Vegetation indices, especially the normalized difference vegetation index (NDVI), are widely used in urban vegetation assessments. However, estimating the vegetation abundance in urban scenes using the NDVI has constraints due to the complex spectral signature related to the urban structure, materials and other factors compared to natural ground surfaces. This paper employs the 3D discrete anisotropic radiative transfer (DART) model to simulate the spectro-directional reflectance of synthetic urban scenes with various urban geometries and building materials using a flux-tracking method under shaded and sunlit conditions. The NDVI is calculated using the spectral radiance in the red (0.6545 μm) and near-infrared bands (0.865 μm). The effects of the urban material heterogeneity and 3D structure on the NDVI, and the performance of three NDVI-based fractional vegetation cover (FVC) inversion algorithms, are evaluated. The results show that the effects of the building material heterogeneity on the NDVI are negligible under sunlit conditions but not negligible under shaded conditions. The NDVI value of building components within synthetic scenes is approximately zero. The shaded road exhibits a higher NDVI value in comparison to the illuminated road because of scattering from adjacent pixels. In order to correct the effects of scattering caused by building geometry, the reflectance of the Landsat 8/OLI image is corrected using the sky view factor (SVF) and then used to calculate the FVC. Jilin-1 satellite images with high spatial resolution (0.5 m) are used to extract the vegetation cover and then aggregated to 30 m spatial resolution to calculate the FVC for validation. The results show that the RMSE is up to 0.050 after correction, while the RMSE is 0.169 before correction. This study makes a contribution to the understanding of the effects of the urban 3D structure and material reflectance on the NDVI and provides insights into the retrieval of the FVC in different urban scenes.@en

The West Sahel is facing significant threats to its vegetation and wildlife due to the land degradation and habitat fragmentation. It is crucial to assess the regional vegetation greenness dynamics in order to comprehensively evaluate the effectiveness of protection in the natural reserves. This study analyzes the vegetation greenness trends and the driving factors in the Dosso Partial Faunal Reserve in Niger and nearby unprotected regions—one of the most important habitats for endemic African fauna—using satellite time series data from 2001 to 2020. An overall vegetation browning trend was observed throughout the entire region with significant spatial variability. Vegetation browning dominated in the Dosso Reserve with 17.7% of the area showing a significant trend, while the area with significant greening was 6.8%. In a comparison, the nearby unprotected regions to the north and the east were found to be dominated by vegetation browning and greening, respectively. These results suggest that the vegetation protection practice was not fully effective throughout the Dosso Reserve. The dominant drivers were also diagnosed using the Random Forest model-based method and the Partial Dependence Plot tool, showing that water availability (expressed as soil moisture) and land use/land cover change were the most critical factors affecting vegetation greenness in the study region. Specifically, soil moisture stress and specific land management practices associated with logging, grazing, and land clearing appeared to dominate vegetation browning in the Dosso Reserve. In contrast, the vegetation greening in the central Dosso Reserve and the nearby unprotected region to the east was probably caused by the increase in shrubland/forest, which was related to the effective implementation of protection. These findings improve our understanding of how regional vegetation greenness dynamics respond to environmental changes in the Dosso Reserve and also highlight the need for more effective conservation planning and implementation to ensure sustainable socio-ecological development in the West Sahel.

@en

Flash droughts tend to cause severe damage to agriculture due to their characteristics of sudden onset and rapid intensification. Early detection of the response of vegetation to flash droughts is of utmost importance in mitigating the effects of flash droughts, as it can provide a scientific basis for establishing an early warning system. The commonly used method of determining the response time of vegetation to flash drought, based on the response time index or the correlation between the precipitation anomaly and vegetation growth anomaly, leads to the late detection of irreversible drought effects on vegetation, which may not be sufficient for use in analyzing the response of vegetation to flash drought for early earning. The evapotranspiration-based (ET-based) drought indices are an effective indicator for identifying and monitoring flash drought. This study proposes a novel approach that applies cross-spectral analysis to an ET-based drought index, i.e., Evaporative Stress Anomaly Index (ESAI), as the forcing and a vegetation-based drought index, i.e., Normalized Vegetation Anomaly Index (NVAI), as the response, both from medium-resolution remote sensing data, to estimate the time lag of the response of vegetation vitality status to flash drought. An experiment on the novel method was carried out in North China during March–September for the period of 2001–2020 using remote sensing products at 1 km spatial resolution. The results show that the average time lag of the response of vegetation to water availability during flash droughts estimated by the cross-spectral analysis over North China in 2001–2020 was 5.9 days, which is shorter than the results measured by the widely used response time index (26.5 days). The main difference between the phase lag from the cross-spectral analysis method and the response time from the response time index method lies in the fundamental processes behind the definitions of the vegetation response in the two methods, i.e., a subtle and dynamic fluctuation signature in the response signal (vegetation-based drought index) that correlates with the fluctuation in the forcing signal (ET-based drought index) versus an irreversible impact indicated by a negative NDVI anomaly. The time lag of the response of vegetation to flash droughts varied with vegetation types and irrigation conditions. The average time lag for rainfed cropland, irrigated cropland, grassland, and forest in North China was 5.4, 5.8, 6.1, and 6.9 days, respectively. Forests have a longer response time to flash droughts than grasses and crops due to their deeper root systems, and irrigation can mitigate the impacts of flash droughts. Our method, based on cross-spectral analysis and the ET-based drought index, is innovative and can provide an earlier warning of impending drought impacts, rather than waiting for the irreversible impacts to occur. The information detected at an earlier stage of flash droughts can help decision makers in developing more effective and timely strategies to mitigate the impact of flash droughts on ecosystems.

@en

Improving irrigation water management is a key concern for the agricultural sector, and it requires extensive and comprehensive tools that provide a complete knowledge of crop water use and requirements. This study presents a novel methodology to explicitly estimate daily gross and net crop water requirements, actual crop water use, and irrigation efficiency of center pivot irrigation systems, by mainly utilizing the Sentinel-2 MultiSpectral Instrument (MSI) imagery at the farm scale. ETMonitor model is adapted to estimate actual water use (as the sum of canopy transpiration and evaporation of water intercepted by canopy and evaporation from soil) at daily/10-m resolution, benefiting from the high-resolution Sentinel-2 data and thus to assess the irrigation efficiency at the farm scale. The gross irrigation water requirement is estimated from the net crop water requirement and the water loss, including the water droplet evaporation directly into the air during application before droplets fall on the canopy and canopy interception loss. The method was applied to a pilot farmland with two major crops (wheat and potato) in the Inner Mongolia Autonomous Region of China, where modern equipment and appropriate irrigation methods are deployed for efficient water use. The estimated actual crop water use showed good agreement with the ground observations, e.g. the determination coefficients range from 0.67 to 0.81 and root mean square errors range from 0.56 mm/day to 1.24 mm/day for wheat and potato when comparing the estimated evapotranspiration with the measurement by the eddy covariance system. It also showed that the losses of total irrigated volume were 25.4% for wheat and 23.7% for potato, respectively, and found that the water allocation was insufficient to meet the water requirement in this irrigated area. This suggests that the amount of water applied was insufficient to meet the crop water requirement and the inherent water losses in the center pivot irrigation system, which imply the necessity to improve the irrigation practice to use the water more efficiently.

@en

Accurate and continuous estimation of surface albedo is vital for assessing and understanding land–surface–atmosphere interactions. We developed a method for estimating instantaneous all-sky at-surface shortwave upwelling radiance and albedo over the Tibetan Plateau. The method accounts for the complex interplay of topography and atmospheric interactions and aims to mitigate the occurrence of data gaps. Employing an RTLSR-kernel-driven model, we retrieved surface shortwave albedo with a 1 km resolution, incorporating direct, isotropic diffuse; circumsolar diffuse; and surrounding terrain irradiance into the all-sky solar surface irradiance. The at-surface upwelling radiance and surface shortwave albedo estimates were in satisfactory agreement with ground observations at four stations in the Tibetan Plateau, with RMSE values of 56.5 W/m² and 0.0422, 67.6 W/m² and 0.0545, 98.6 W/m² and 0.0992, and 78.0 98.6 W/m² and 0.639. This comparison indicated an improved accuracy of at-surface upwelling radiance and surface albedo and significantly reduced data gaps. Valid observations increased substantially in comparison to the MCD43A2 data product, with the new method achieving an increase ranging from 40% to 200% at the four stations. Our study demonstrates that by integrating terrain, cloud properties, and radiative transfer modeling, the accuracy and completeness of retrieved surface albedo and radiance in complex terrains can be effectively improved.

@en

Introduction: The vegetation dynamics of the Sahel-Sudan-Guinea region in Africa, one of the largest transition zones between arid and humid zones, is of great significance for understanding regional ecosystem changes. However, a time-unvarying trend based on linear assumption challenges the overall understanding of vegetation greenness evolution and of tracking a complex ecosystem response to climate in the Sahel-Sudan-Guinea region. Methods: This study first applied the ensemble empirical mode decomposition (EEMD) method to detect the time-varying trends in vegetation greenness based on normalized difference vegetation index (NDVI) data in the region during 2001–2020, and then identified the dominant climatic drivers of NDVI trends by employing explainable machine learning framework. Results: The study revealed an overall vegetation greening but a significant nonlinear spatio-temporal evolution characteristic over the region. Trend reversals, i.e., browning-to-greening and greening-to-browning, were dominant in approximately 60% of the study area. The browning-to-greening reversal was primarily observed in the southern Sahel, Congo Basin north of the Equator, and East Africa, with a breakpoint around 2008, while the greening-to-browning reversal was mainly observed in West Africa, with a breakpoint around 2011. The sustained greening primarily took place in northern Sahel, Central African Republic and South Sudan; while sustained browning clustered in central West Africa and Uganda, mainly in agricultural lands. Furthermore, the combination of Random Forest (RF) algorithm and the SHapley Additive exPlanations (SHAP) method could robustly model and reveal the relationships between the observed trends in NDVI and in climatic variables, also detected by applying EEMD. The results suggested that air temperature and precipitation were the most important climatic drivers controlling the NDVI trends across the Sahel-Sudan-Guinea region. The NDVI trends were more likely to have negative correlations with solar radiation and vapor pressure deficit in arid areas, while they could have positive correlations in humid areas. The study also found that large-scale climate changes induced by sea surface temperature (SST) anomalies had strong relationships with trend reversals in vegetation greenness at a sub-continental scale. These findings advanced the understanding of the impacts of climatic drivers on vegetation greenness evolution in the Sahel-Sudan-Guinea region.

@en

Trees are the most important natural factor to alleviate the urban heat island effect and the latent heat flux (LE) they release contributes significantly to urban cooling. In this study, the model ENVI-met was used to study the influence of building geometry on the LE exchanged by trees at the block scale in compact urban areas. The building density (BD), building height (BH) and sky view factor (SVF) were used to characterize building geometry. The sensitivity of LE to building geometry was estimated by multi-linear regression analysis. The following conclusions were drawn:(1) the LE of trees is sensitive to building geometry, higher during daytime than night-time in winter and the higher during night-time than daytime in summer; (2) LE changes as expected across the seasons, with LE larger in summer; (3) The shade affects the LE of a tree by influencing the solar irradiance; (4) In winter the average LE of a single tree is larger with 0.06 than 0.12 fractional vegetation cover (FVC). This study can provide useful leads towards further research to explore latent heat exchanges by trees at the street scale.

@en

The global validation for SMAP (Soil Moisture Active Passive) Level-4 surface soil moisture using well-established core validation sites does not comprehensively account for all landscapes on earth’s surface due to the diverse nature of their intrinsic characteristics. Due to the inhibitive cost of implementing a standard validation site, this study presents an alternative sampling approach suitable for localized validation of SMAP data in non-complex terrains. It involves clustering a large heterogeneous study area to smaller units of non-complex terrains where landscape-defining characteristics are largely homogeneous, thus permitting the computation of areal soil moisture as a simple arithmetic mean of near-linear point measurements. This allows optimization of limited resources (a hand-held moisture sensor with site-specific calibrations) to balance the need for spatial representativeness of samples in the sampling unit and the need for temporal proximity of sampled measurements to the aggregation time of the satellite product being validated. For ease of movement, transect sampling is implemented along access roads that run across the sampling units to allow sufficient measurement replications within a reasonably short time. Validations with four different landscapes in Kenya show a good agreement between in situ measurements and SMAP with R2 of 0.76, 0.72, 0.80, and 0.82, and biases of −0.0246, +0.0113, 0.0004, and + 0.0035 m3 m−3, respectively for Mawego, Kuresoi, Sotik and Kapsuser sites. These results only marginally differ from those obtained with a spatially distributed sampling method, indicating the potential of the proposed sampling design for time and cost effectiveness in validations at non-complex terrains. An analysis of the temporal variability of SMAP soil moisture in the watershed is also presented, with an assessment of its significance in the selection of sampling sites for validation. Particularly, the concept of temporal stability of soil moisture as a basis for clustering validation sites is evaluated.@en

Around 90% of the oceanic and inland waters’ reflectance registered in satellite detectors comes from the atmospheric contribution. Hence the water-leaving radiances in the Near-InfraRed (NIR) region are above the zero value over inland waters because of sediments and dissolved organic particles, this radiance cannot be ignored. To accurately retrieve water quality parameters from water-leaving reflectance, atmospheric correction is the most important step. This study evaluated five reliable atmospheric correction algorithms (AC) known as: (ACOLITE, C2RCC, iCOR, 6SV, and Sen2Cor) against optical in-situ measurements collected above the water in Qiandao Lake, China using Sentinel-2 Multi-Spectral Imager. 60 in-situ water samples and optical measurements (range 400–900 nm) above the water were collected at different points in Qiandao Lake. The spectra measurements were used to validate the atmospheric correction processors. All ACs that were evaluated showed high levels of uncertainty. ACOLITE and ICOR performed the best statistics with root mean square differences (RMSD) (0.006 sr⁻¹) while Sen2Cor achieved the lowest RMSD (0.023 sr⁻¹) across the different modules. ACOLITE, had a better performance when applied to meso- and hypereutrophic waters, compared with oligotrophic, while C2RCC performs better at the wavelength of 833 nm (0.007 sr⁻¹). Finally, 6S performs better at the wavelength of 665 nm (0.015 sr⁻¹). This study introduces insights and addresses a significant research gap in the field of atmospheric correction for satellite imagery over inland waters. Prior studies have primarily focused on atmospheric correction algorithms for coastal and open ocean environments while few studies focused on the unique characteristics and challenges associated with inland water bodies. The findings of this study are crucial for researchers, remote sensing experts, and environmental scientists working with Sentinel-2A imagery, as it enables them to make more accurate and reliable interpretations of water quality and other environmental parameters derived from satellite data.

@en

Accurate knowledge of the at-surface solar irradiance (SSI) is essential for retrieving surface and atmospheric properties using satellite measurements of backscattered and reflected radiance. The latter is affected by surface-atmosphere interactions, including the effects of terrain. The SSI is affected by the same processes. This study proposes a method to estimate the components of instantaneous SSI: direct, isotropic and circumsolar diffuse, and terrain irradiance, which is expected to improve the simultaneous retrieval of aerosol optical depth (AOD) and surface reflectance. The method takes into account the coupled effects of topography and atmosphere by combining parameterization and the lookup table (LUT) approaches. The method was applied to rugged terrain over the Tibetan Plateau using Moderate Resolution Imaging Spectrometer (MODIS) atmosphere and surface data, the fifth generation European Centre for Medium-Range Weather Forecasts reanalysis (ERA5) data, Cloud-Aerosol Lidar With Orthogonal Polarization (CALIOP) aerosol data, and a digital elevation model (DEM). The results showed that the SSI estimates were in satisfactory agreement with ground observations at four stations over the Tibetan Plateau (TP) in 2018 with R2 values of 0.61, 0.44, 0.41, and 0.49, respectively, and root mean square error (RMSE) of 205.7, 176.9, 186.0, and 201.2 W/m2, respectively. Estimations of the diffuse irradiance were evaluated separately against the only available in situ observations at the Dali Station, and the results were better than our SSI estimates with R2, RMSE, and relative bias (BIAS) being 0.71, 94.98 W/m2, and 31%, respectively. The isotropic and circumsolar diffuse irradiances accounted for 37.57% and 7.68% of the total annual SSI, respectively, while diffuse irradiance accounted for 46.48% of the total annual SSI. Under clear skies, every 0.1 increase in AOD caused about a 35-W/m2 increase in diffuse irradiance and a decrease of about 25 W/m2 of SSI.

@en

Mountains are important suppliers of freshwater to downstream areas, affecting large populations in particular in High Mountain Asia (HMA). Yet, the propagation of water from HMA headwaters to downstream areas is not fully understood, as interactions in the mountain water cycle between the cryo-, hydro- and biosphere remain elusive. We review the definition of blue and green water fluxes as liquid water that contributes to runoff at the outlet of the selected domain (blue) and water lost to the atmosphere through vapor fluxes, that is evaporation from water, ground, and interception plus transpiration (green) and propose to add the term white water to account for the (often neglected) evaporation and sublimation from snow and ice. We provide an assessment of models that can simulate the cryo-hydro-biosphere continuum and the interactions between spheres in high mountain catchments, going beyond disciplinary separations. Land surface models are uniquely able to account for such complexity, since they solve the coupled fluxes of water, energy, and carbon between the land surface and atmosphere. Due to the mechanistic nature of such models, specific variables can be compared systematically to independent remote sensing observations–providing vital insights into model accuracy and enabling the understanding of the complex watersheds of HMA. We discuss recent developments in spaceborne earth observation products that have the potential to support catchment modeling in high mountain regions. We then present a pilot study application of the mechanistic land surface model Tethys & Chloris to a glacierized watershed in the Nepalese Himalayas and discuss the use of high-resolution earth observation data to constrain the meteorological forcing uncertainty and validate model results. We use these insights to highlight the remaining challenges and future opportunities that remote sensing data presents for land surface modeling in HMA.

@en

In the Tibetan Plateau (TP) region, the foreseeable increase in air temperature may have profound and complex effects on the local hydrological cycle, and is likely to increase water loss from the land surface to the atmosphere through evapotranspiration (ET). Quantifying ET and its regulatory mechanisms are major challenges for understanding the water cycle and land-atmosphere interactions in the TP region. We evaluated the performance of several Earth observation-based ET datasets in the TP region, and explored the spatiotemporal variation of ET in the same region. The accuracy of different global ET datasets was evaluated, and ETMonitor and PML-V2 provide the best accuracy with overall high correlation, low bias, and low root mean square error. ETMonitor ET is also the only product with both high spatial (~1 km) and temporal (daily) resolution. ETMonitor ET may reflect the effect of mountain topography on ET better than other global products, i.e., ET values are higher in the humid valleys with denser vegetation cover and higher soil moisture, and ET values are lower on the mountain slopes at higher elevations with less vegetation cover and colder climate. Other ET products failed to capture the spatial patterns of ET in the mountainous regions, and this suggests that the spatial resolution is not the only dominant factor leading to the poorer performance of these ET products in the mountain regions of the TP. The results show that multi-year average ET is 339 mm/yr in the TP region during 2000-2021, which accounts for about 51% of the total precipitation in the TP region. From 2000 to 2021, ET over the Tibetan Plateau shows an overall increasing trend with large spatial variability.

@en

Water depth, a fundamental characteristic of a lake, is important for understanding climatic, ecological, and hydrological processes. However, lake water depth data are still scarce due to the high cost of in-situ measurements and the limitations of remote sensing observations. In this study, a novel method was developed to estimate time series of pixel-wise water depths of lakes that have ever exposed their bottom by remote sensing observations. Lake water depths were calculated as the difference between the elevations of the dynamic water surface and the historical lakebed elevations using optical images and DEM data. The method was applied in the Sahel-Sudano-Guinean region of Africa where complex climatic conditions and rare in-situ measurements. Experiments showed that the proposed method could get consistent water depths compared with the HydroLAKES data, i.e. with a MAE of 0.86 m and a RMSE of 1.69 m, and water surface elevations similar to the estimates derived from ICESat/ICESat-2 measurements with a MAE of 3.79 m and a RMSE of 5.92 m. The method can provide pixel-wise information on lake water depth at high temporal frequency, and is expected to provide an efficient solution to gather essential information on lakes.

@en

Glaciers are crucial water resources in the Third Pole (the Tibetan Plateau and its surroundings) and are shrinking in response to climate change. Glacier albedo is an expression of glacier interactions with climate and dust/black carbon, and albedo reduction enhances glacier mass loss, but its changes and potential drivers remain poorly quantified. We leverage satellite observations to explore the variability of glacier albedo and understand its sensitivity to potential drivers and its future evolution. We find that glacier albedo has declined during 2001–2020, but high interannual variability is also an important signal. These variations are highly sensitive to air temperature and snow conditions and to nearby dust/black carbon emission sources. Future changes to these drivers will lead to further decreases of 2.9%–12.5% in glacier albedo by 2100 under different warming scenarios. These findings highlight the importance of albedo in glacier future evolution and the urgency of action to mitigate climate warming.

@en

Spatiotemporal residual noise in terrestrial earth observation products, often caused by unfavorable atmospheric conditions, impedes their broad applications. Most users prefer to use gap-filled remote sensing products with time series reconstruction (TSR) algorithms. Applying currently available implementations of TSR to large-volume datasets is time-consuming and challenging for non-professional users with limited computation or storage resources. This study introduces a new open-source software package entitled ‘HANTS-GEE’ that implements a well-known and robust TSR algorithm, i.e. Harmonic ANalysis of Time Series (HANTS), on the Google Earth Engine (GEE) platform for scalable reconstruction of terrestrial earth observation data. Reconstruction tasks can be conducted on user-defined spatiotemporal extents when raw datasets are available on GEE. According to site-based and regional-based case evaluation, the new tool can effectively eliminate cloud contamination in the time series of earth observation data. Compared with traditional PC-based HANTS implementation, the HANTS-GEE provides quite consistent reconstruction results for most terrestrial vegetated sites. The HANTS-GEE can provide scalable reconstruction services with accelerated processing speed and reduced internet data transmission volume, promoting algorithm usage by much broader user communities. To our knowledge, the software package is the first tool to support full-stack TSR processing for popular open-access satellite sensors on cloud platforms.

@en

The potential drivers of vegetation changes in the Sudano-Sahelian region of Africa remain poorly understood due to complex interactions between climatic and anthropogenic processes. In this study, we analyzed the vegetation greenness trends in relation to rainfall variability that we considered the essence of climatic effects on vegetation in a well-known water-limited environment by using time series of satellite data in the Sudano-Sahelian region during 2001–2020. We quantified in more detail the relative contributions of rainfall variability (climatic factor), land use/land cover (LULC) change, and fire occurrence change (non-climatic factors) to vegetation greenness trends in selected sub-regions. The results showed that vegetation greening was widespread (26.9% of the total study area), while vegetation browning was more clustered in central West Africa (5% of the total study area). About half of the vegetation greening area can be explained by long-term rainfall variability during 2001–2020, but most of the area characterized by a browning trend was unrelated to rainfall variability. An analysis of the relative importance showed that LULC changes had significant local effects on vegetation greenness and that these changes were characterized by a strong spatial heterogeneity in specific sub-regions. Gains in cropland and natural vegetation related to positive land management were probably the dominant drivers of greening in Senegal and Ethiopia. Also, the combined impacts of rainfall variability and LULC changes contributed to greening trends in the arid zone, particularly in Mali and Sudan. In contrast, vegetation browning in central West Africa appeared to be driven by cropland gain and natural vegetation loss associated with extensive agricultural production activities. Furthermore, we found that repeated fires for agricultural expansion in central West Africa intensified vegetation browning. These results advanced our understanding of vegetation dynamics in response to climatic and non-climatic factors in Sudano-Sahelian drylands characterized by increasing pressures on land resources.

@en

The global climate has undergone unequivocal warming [...]@en

Urban morphology affects the sensible heat flux and net radiation exchange which can alter urban heat mitigation plans. This study first parameterized the geometric effects on the net radiation, and then calculated the net radiation and sensible heat flux in the urban landscape of Hong Kong. Considering that the sensible heat flux is the main heat sink in compact urban areas, this study proposes a Normalized Urban Sensible Heat Mitigation Index (NUSHMI) based on the ratio of the net radiation and sensible heat flux. Overall, there is major difference in the dependence of net radiation and sensible heat flux on geometric parameters. Net radiation Rn, reaches an optimal value, either maximum or minimum depending on the parameters of SVF and a standard deviation of building height σ_h, at intermediate parameter values, which suggests a guideline relevant to urban design targeting the mitigation of urban climate. Contrariwise, sensible heat flux decreases or increases, again depending on SVF and σ_h, is being considered, with increasing values of the same parameters. For example, R_n, reaches a minimum value for a Sky View Factor (SVF) between 0.5 and 0.6, while it reaches a maximum value for a standard deviation of building height σ_h between 20 and 30 m. These two results suggest that radiative forcing, i.e. R_n, can be minimized by urban space with SVF around 0.55 and σ_h around 25 m. The relationships between sensible heat flux and SVF or σ_h do not show multiple minima or maxima (as with R_n), with the exception of building density, which could also be applied as a guideline in urban design. The results based on the proposed NUSHMI indicated the NUSHMI reaches the highest values when building density is about 0.7 and building height is about 80 m and when the building height standard deviation within an area is about 10 m to 20 m. These findings revealed how the urban morphology affects the surface heat flux exchange between urban canopy and atmosphere boundary layer, and can help to design an efficient urban landscape towards urban heat mitigation for highly compacted cities, e.g. controlling the building density, height, and the height deviation. This combination of urban geometric parameters identifies an urban configuration maximizing the dissipation of absorbed radiant energy as sensible heat. It should be noted, however, that heat load upon buildings would be reduced at the price of maximizing heat dissipation within the built-up space.

@en

The urban heat island effect poses a growing threat to human society, especially in densely populated and developed megacities. With the introduction of the Local Climate Zones (LCZ) framework, new perspectives and findings have been brought to urban heat island studies. This study investigated the cooling effect of vegetation and albedo on the surface urban heat island (SUHI) in the classification system of LCZ during different seasons, using three Chinese megacities as case study areas. Single-factor linear regression and Pearson’s correlation coefficient were applied to analyze the seasonal cooling effect of both albedo and the NDVI on the SUHI within different LCZs. The results show that (1) the variability of the SUHI is reflected in its dominance and intensity within certain LCZs in different cities and in the efficiency of cooling factors; (2) the cooling effect of vegetation is dominant in each season, and the cooling effect produced by albedo within specific seasons can be differentiated by LCZs. This study provides valuable information for the mitigation of the SUHI magnitude in specific regions and at specific times of the year.@en

Information on crop yield is important for food security, in particular under the conditions of climate change and growing population worldwide. We developed a new fully distributed, high spatial resolution, model of biomass accumulation and crop yield applicable to a highly heterogeneous desert-oasis agroecosystem. The bulk of required input data is obtained by retrieving pixel-wise biogeophysical variables from a suite of very diverse satellite data. Both temperature and water stress conditions at field-scale are given full consideration, while the model was designed to strike a balance between model applicability and satisfactory characterization of the heterogeneous desert-oasis system to estimate field-scale yield. The development of this model relies on three main innovations. First, the start and end of the growing season were estimated for each pixel by calibrating the high spatial and temporal resolution observations of Normalized Difference Vegetation Index (NDVI) by Sentinal-2 (S2) MSI (Multi-Spectral Instrument) against limited local phenological information. Second, to monitor crop water stress, account taken of irrigation, a process-based water and energy balance model was applied to estimate the actual evapotranspiration (ET). This requires knowledge of soil water availability, which is characterized by downscaling the ASCAT (Advanced SCATterrometer) soil moisture data product. To capture the dominant features of the eco-hydrological conditions in the desert and oasis agroecosystem, ET was further downscaled from the 1 km resolution. Third, likewise the water stress indicator, the air temperature stress indicator was mapped after characterizing the thermal contrast and heterogeneity of the desert-oasis system, by generating time series of air temperature at 1 km spatial resolution using the MODIS (Moderate Resolution Imaging Spectroradiometer) Land Surface Temperature (LST) data product. In the temporal dimension, gaps were mitigated by applying time series analysis techniques to reconstruct cloud-free time series of LST, NDVI, fAPAR and albedo. These innovations add up to a high resolution characterization of crop response to the geospatial variability of weather and climate forcing in the desert-oasis agroecosystem. The model was applied to the dominant crops, i.e., spring wheat, maize, sunflower, and melon, in the oases of the Shiyang River Basin (northwestern China) characterized by a rather fragmented land use. The high resolution of pixel-wise ecohydrological parameters, i.e., crop phenology, temperature stress and water stress factors successfully reflect differences of crops with different phenology and location in the oases. The relative errors for wheat and maize yields compared to the census data are less than 5% at district level. At the county level, the relative errors of wheat yields of Liangzhou, Minqin, Gulang, Jinchuan, and Yongchang equal to 0.87%, 24.2%, 9.7%, 12.5%, and 7.2%. For maize, the dominant crop, the error on estimated yields was less than 5%, except in Gulang. The relative error on estimated yield for sunflower was less than 10% compared to agricultural census data. The relative error on estimated melon yield was 16%. This performance highlights the applicability of the model to estimate field-scale yields in agroecosystems characterized by fragmented land use.

@en