This report is for a capstone project for Cornell MPS, concentration Geospatial Applications, to illustrate mastery of a complex geospatial analysis. It used a wide variety of techniques and tools. The project investigated an important topic in the context of the urban environment under climate variability: the spatial-temporal distribution of vegetation phenology in Beijing, focusing on the impacts of local climate zones (LCZ) and urbanization levels. Using MODIS MCD12Q2 data, we analyzed the start of season (SOS), end of season (EOS), and length of season (LOS) across different urbanization gradients and LCZs. Results reveal that urbanized regions exhibit earlier SOS and longer LOS, while natural areas have a later EOS and shorter LOS. LCZs with lower density buildings tend to have a more extended growing season. These findings underscore the significance of urbanization and LCZ in shaping vegetation phenology, providing valuable insights into climate change adaptation and urban ecological management.