Nothing Special   »   [go: up one dir, main page]

You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (280)

Search Parameters:
Keywords = land surface albedo

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
15 pages, 1077 KiB  
Technical Note
Quantifying Annual Glacier Mass Change and Its Influence on the Runoff of the Tuotuo River
by Lin Liu, Xueyu Zhang and Zhimin Zhang
Remote Sens. 2024, 16(20), 3898; https://doi.org/10.3390/rs16203898 - 20 Oct 2024
Viewed by 530
Abstract
Glacier meltwater is an indispensable water supply for billions of people living in the catchments of major Asian rivers. However, the role of glaciers on river runoff regulation is seldom investigated due to the lack of annual glacier mass balance observation. In this [...] Read more.
Glacier meltwater is an indispensable water supply for billions of people living in the catchments of major Asian rivers. However, the role of glaciers on river runoff regulation is seldom investigated due to the lack of annual glacier mass balance observation. In this study, we employed an albedo-based model with a daily land surface albedo dataset to derive the annual glacier mass balance over the Tuotuo River Basin (TRB). During 2000–2022, an annual glacier mass balance range of −0.89 ± 0.08 to 0.11 ± 0.11 m w.e. was estimated. By comparing with river runoff records from the hydrometric station, the contribution of glacier mass change to river runoff was calculated to be 0.00–31.14% for the studied period, with a mean value of 9.97%. Moreover, we found that the mean contribution in drought years is 20.07%, which is approximately five times that in wet years (4.30%) and twice that in average years (9.49%). Therefore, our results verify that mountain glaciers act as a significant buffer against drought in the TRB, at least during the 2000–2022 period. Full article
Show Figures

Figure 1

Figure 1
<p>The geographic location of the Tuotuo River Basin (Olivine Yellow) and the three study sites (A, B, and C). Glacier boundaries were obtained from the second Chinese glacier inventory. The location of the Tuotuohe hydrometric station is indicated as a red five-pointed star.</p>
Full article ">Figure 2
<p>The extracted annual minimum regional-average surface albedo time series for the three study sites in the TRB between 2000 and 2022. The red vertical dashed line indicates the year 2012.</p>
Full article ">Figure 3
<p>The observed annual runoff of the Tuotuo River between 2000 and 2022.</p>
Full article ">Figure 4
<p>Annual basin-wide average precipitation (<b>a</b>) and evaporation (<b>b</b>) for the Tuotuo River Basin between 2000 and 2022.</p>
Full article ">Figure 5
<p>Comparison between annual glacier mass change and annual river runoff for the Tuotuo River Basin during the period of 2000–2022.</p>
Full article ">Figure 6
<p>The amount of annual precipitation minus annual evaporation in 2000–2022.</p>
Full article ">Figure 7
<p>Relationship between river runoff and the difference in precipitation and evaporation over the Tuotuo River Basin during the period of 2000–2022.</p>
Full article ">
21 pages, 7794 KiB  
Article
Spatial–Temporal Variations and Driving Factors of the Albedo of the Qilian Mountains from 2001 to 2022
by Huazhu Xue, Haojie Zhang, Zhanliang Yuan, Qianqian Ma, Hao Wang and Zhi Li
Atmosphere 2024, 15(9), 1081; https://doi.org/10.3390/atmos15091081 - 6 Sep 2024
Viewed by 620
Abstract
Surface albedo plays a pivotal role in the Earth’s energy balance and climate. This study conducted an analysis of the spatial distribution patterns and temporal evolution of albedo, normalized difference vegetation index (NDVI), normalized difference snow index snow cover (NSC), and land surface [...] Read more.
Surface albedo plays a pivotal role in the Earth’s energy balance and climate. This study conducted an analysis of the spatial distribution patterns and temporal evolution of albedo, normalized difference vegetation index (NDVI), normalized difference snow index snow cover (NSC), and land surface temperature (LST) within the Qilian Mountains (QLMs) from 2001 to 2022. This study evaluated the spatiotemporal correlations of albedo with NSC, NDVI, and LST at various temporal scales. Additionally, the study quantified the driving forces and relative contributions of topographic and natural factors to the albedo variation of the QLMs using geographic detectors. The findings revealed the following insights: (1) Approximately 22.8% of the QLMs exhibited significant changes in albedo. The annual average albedo and NSC exhibited a minor decline with rates of −0.00037 and −0.05083 (Sen’s slope), respectively. Conversely, LST displayed a marginal increase at a rate of 0.00564, while NDVI experienced a notable increase at a rate of 0.00178. (2) The seasonal fluctuations of NSC, LST, and vegetation collectively influenced the overall albedo changes in the Qilian Mountains. Notably, the highly similar trends and significant correlations between albedo and NSC, whether in intra-annual monthly variations, multi-year monthly anomalies, or regional multi-year mean trends, indicate that the changes in snow albedo reflected by NSC played a major role. Additionally, the area proportion and corresponding average elevation of PSI (permanent snow and ice regions) slightly increased, potentially suggesting a slow upward shift of the high mountain snowline in the QLMs. (3) NDVI, land cover type (LCT), and the Digital Elevation Model (DEM, which means elevation) played key roles in shaping the spatial pattern of albedo. Additionally, the spatial distribution of albedo was most significantly influenced by the interaction between slope and NDVI. Full article
(This article belongs to the Special Issue Vegetation and Climate Relationships (3rd Edition))
Show Figures

Figure 1

Figure 1
<p>Geographic location and elevation of the QLMs.</p>
Full article ">Figure 2
<p>(<b>a</b>) Slope, (<b>b</b>) aspect, and (<b>c</b>) land cover types of the QLMs (the explanation of the abbreviation is included in <a href="#atmosphere-15-01081-t001" class="html-table">Table 1</a>).</p>
Full article ">Figure 3
<p>The multi-year average (<b>a</b>) albedo, (<b>b</b>) NSC, (<b>c</b>) NDVI, and (<b>d</b>) LST.</p>
Full article ">Figure 4
<p>Trends in the annual average albedo of the QLMs from 2001 to 2022 in relation to (<b>a</b>) annual average NSC, (<b>b</b>) annual average NDVI, (<b>c</b>) annual average LST; (<b>d</b>) trends in the average elevation and area percentage of PSI regions in the QLMs.</p>
Full article ">Figure 5
<p>The spatial distribution of significant or non-significant changes in (<b>a</b>) albedo, (<b>b</b>) NSC, (<b>c</b>) NDVI, and (<b>d</b>) LST in the QLMs from 2001 to 2022.</p>
Full article ">Figure 6
<p>Comparison of the multi-year monthly average values of albedo with (<b>a</b>) NSC, (<b>b</b>) NDVI, and (<b>c</b>) LST in the QLM region from 2001 to 2022.</p>
Full article ">Figure 7
<p>Monthly albedo anomalies (<b>a</b>) and monthly NSC anomalies, (<b>b</b>) and monthly NDVI anomalies, (<b>c</b>) and monthly LST anomalies in the QLMs from 2001 to 2022.</p>
Full article ">Figure 8
<p>(<b>a</b>–<b>e</b>) Changes in explanatory power (q values) in 2001, 2006, 2011, 2016, and 2021.</p>
Full article ">Figure 9
<p>Striking differences in the driving factors (At a confidence level of 95%, “Y” indicates a significant difference in the spatial distribution of albedo due to the two factors, while “N” indicates the opposite).</p>
Full article ">Figure 10
<p>(<b>a</b>–<b>e</b>) Changes in interactive explanatory power (q values) in 2001, 2006, 2011, 2016, and 2021; (<b>f</b>) average interactive explanatory power (q values) of 5 years (Bi: Enhance, bivariate, ENL: Enhance, nonlinear. The annotations inside parentheses indicate a higher frequency of occurrence of interaction types within five years. Without annotations, it indicates that the interaction types remained consistent over the 5 years).</p>
Full article ">
29 pages, 19031 KiB  
Article
Directional Applicability Analysis of Albedo Retrieval Using Prior BRDF Knowledge
by Hu Zhang, Qianrui Xi, Junqin Xie, Xiaoning Zhang, Lei Chen, Yi Lian, Hongtao Cao, Yan Liu, Lei Cui and Yadong Dong
Remote Sens. 2024, 16(15), 2744; https://doi.org/10.3390/rs16152744 - 26 Jul 2024
Viewed by 629
Abstract
Surface albedo measures the proportion of incoming solar radiation reflected by the Earth’s surface. Accurate albedo retrieval from remote sensing data usually requires sufficient multi-angular observations to account for the surface reflectance anisotropy. However, most middle and high-resolution remote sensing satellites lack the [...] Read more.
Surface albedo measures the proportion of incoming solar radiation reflected by the Earth’s surface. Accurate albedo retrieval from remote sensing data usually requires sufficient multi-angular observations to account for the surface reflectance anisotropy. However, most middle and high-resolution remote sensing satellites lack the capability to acquire sufficient multi-angular observations. Existing algorithms for retrieving surface albedo from single-direction reflectance typically rely on land cover types and vegetation indices to extract the corresponding prior knowledge of surface anisotropic reflectance from coarse-resolution Bidirectional Reflectance Distribution Function (BRDF) products. This study introduces an algorithm for retrieving albedo from directional reflectance based on a 3 × 3 BRDF archetype database established using the 2015 global time-series Moderate Resolution Imaging Spectro-radiometer (MODIS) BRDF product. For different directions, BRDF archetypes are applied to the simulated MODIS directional reflectance to retrieve albedo. By comparing the retrieved albedos with the MODIS albedo, the BRDF archetype that yields the smallest Root Mean Squared Error (RMSE) is selected as the prior BRDF for the direction. A lookup table (LUT) that contains the optimal BRDF archetypes for albedo retrieval under various observational geometries is established. The impact of the number of BRDF archetypes on the accuracy of albedo is analyzed according to the 2020 MODIS BRDF. The LUT is applied to the MODIS BRDF within specific BRDF archetype classes to validate its applicability under different anisotropic reflectance characteristics. The applicability of the LUT across different data types is further evaluated using simulated reflectance or real multi-angular measurements. The results indicate that (1) for any direction, a specific BRDF archetype can retrieve a high-accuracy albedo from directional reflectance. The optimal BRDF archetype varies with the observation direction. (2) Compared to the prior BRDF knowledge obtained through averaging method, the BRDF archetype LUT based on the 3 × 3 BRDF archetype database can more accurately retrieve the surface albedo. (3) The BRDF archetype LUT effectively eliminates the influence of surface anisotropic reflectance characteristics in albedo retrieval across different scales and types of data. Full article
Show Figures

Figure 1

Figure 1
<p>The flowchart for albedo retrieval from directional reflectance based on BRDF archetypes. The red, blue and green lines in the validation section represent the processing workflows for different validation data.</p>
Full article ">Figure 2
<p>The shapes of BRDF archetypes (red line) on the PP at an SZA of 45° for the NIR band. (<b>a</b>–<b>i</b>) refer to the nine BRDF archetype classes. The gray lines refer to 100 normalized MODIS BRDF selected from each BRDF archetype class randomly.</p>
Full article ">Figure 3
<p>The viewing zenith angles (<b>a</b>) and azimuth angles (<b>b</b>) of the LUT. The radius represents the zenith angle, and the polar angle represents the azimuth angle. Each point represents a direction, and different colors represent the magnitudes of the angles.</p>
Full article ">Figure 4
<p>Angular sampling of multi-angular observations. (<b>a</b>,<b>b</b>) refer to the MODIS observations within 2021.101–2021.116 and 305–320, (<b>c</b>) shows the angular distribution of POLDER data named ‘brdf_ndvi03_0634_2286.txt’, and (<b>d</b>) represents the angular distribution pattern of ground measurements named ‘Parabola.1987.ifc3-site36.inp’. Solid dots represent the locations of the view, and the red open circles refer to the locations of the sun.</p>
Full article ">Figure 5
<p>The comparison of albedos retrieved from different BRDF archetypes and directional reflectance with MODIS albedo in the NIR band. (<b>a</b>–<b>i</b>) represent the inversion results for the nine BRDF archetypes, respectively. The observation is positioned with an SZA of 45° and a VZA of 55° in the backward direction of the PP. The color represents the density of overlapping points.</p>
Full article ">Figure 6
<p>The comparison between directional reflectance (<b>a</b>–<b>d</b>) or albedo (<b>e</b>–<b>h</b>) retrieved from the BRDF archetype with the least RMSE and MODIS albedo in the NIR band. (<b>a</b>,<b>e</b>) represent the direction with an VZA of 55° in the backward direction of PP; (<b>b</b>,<b>f</b>) represent the forward direction of 45° in PP; (<b>c</b>,<b>g</b>) represent the nadir direction; (<b>d</b>,<b>h</b>) represent the direction of 60° in CPP. The color represents the density of overlapping points.</p>
Full article ">Figure 7
<p>The distribution of RMSE<sub>r</sub> (<b>a</b>–<b>d</b>) and RMSE<sub>a</sub> (<b>e</b>–<b>h</b>) in the red band over the viewing hemisphere under SZA of 5° (<b>a</b>,<b>e</b>), 30° (<b>b</b>,<b>f</b>), 45° (<b>c</b>,<b>g</b>), and 60° (<b>d</b>,<b>h</b>). The radius represents the zenith angle, and the polar angle represents the azimuth angle.</p>
Full article ">Figure 8
<p>The distribution of RMSE<sub>r</sub> (<b>a</b>–<b>d</b>) and RMSE<sub>a</sub> (<b>e</b>–<b>h</b>) in the NIR band over the viewing hemisphere under SZA of 5° (<b>a</b>,<b>e</b>), 30° (<b>b</b>,<b>f</b>), 45° (<b>c</b>,<b>g</b>), and 60° (<b>d</b>,<b>h</b>).</p>
Full article ">Figure 9
<p>The BRDF archetype LUTs for the red (<b>a</b>–<b>d</b>) and NIR (<b>e</b>–<b>h</b>) bands for retrieving BSA based on directional reflectance under SZA of 5° (<b>a</b>,<b>e</b>), 30° (<b>b</b>,<b>f</b>), 45° (<b>c</b>,<b>g</b>), and 60° (<b>d</b>,<b>h</b>). Different colors represent different BRDF archetypes.</p>
Full article ">Figure 10
<p>The BRDF archetype LUTs for the red (<b>a</b>–<b>d</b>) and NIR (<b>e</b>–<b>h</b>) bands for retrieving WSA based on directional reflectance under SZA of 5° (<b>a</b>,<b>e</b>), 30° (<b>b</b>,<b>f</b>), 45° (<b>c</b>,<b>g</b>), and 60° (<b>d</b>,<b>h</b>). Different colors represent different BRDF archetypes.</p>
Full article ">Figure 11
<p>The 3D pattern of mean BRDF at an SZA of 30°. (<b>a</b>) is the red band, and (<b>b</b>) is the NIR band. Colors represent the magnitude of reflectance.</p>
Full article ">Figure 12
<p>The distribution of RMSE<sub>a</sub> based on the mean BRDF in the red (<b>a</b>–<b>d</b>) and NIR (<b>e</b>–<b>h</b>) bands over the viewing hemisphere under SZA of 5° (<b>a</b>,<b>e</b>), 30° (<b>b</b>,<b>f</b>), 45° (<b>c</b>,<b>g</b>), and 60° (<b>d</b>,<b>h</b>).</p>
Full article ">Figure 13
<p>The proportions of directions with RMSE less than 0.025 and 0.045 in the red and NIR bands. (<b>a</b>,<b>b</b>) refer to the BSA and WSA, respectively. The RMSEs are calculated based on directional reflectance, mean BRDF, and LUTs established using 6 × 1, 2 × 2, 3 × 3, and 5 × 5 BRDF archetypes.</p>
Full article ">Figure 14
<p>Validation based on 2020 MODIS BRDF product. The distribution of RMSE<sub>a</sub> in the red (<b>a</b>–<b>d</b>) and NIR (<b>e</b>–<b>h</b>) bands over the viewing hemisphere under SZA of 5° (<b>a</b>,<b>e</b>), 30° (<b>b</b>,<b>f</b>), 45° (<b>c</b>,<b>g</b>), and 60° (<b>d</b>,<b>h</b>).</p>
Full article ">Figure 15
<p>The distribution of RMSE<sub>a</sub> over the viewing hemisphere within each BRDF archetype class at an SZA of 45°. (<b>a</b>–<b>i</b>) represent the nine BRDF archetype classes, respectively.</p>
Full article ">Figure 15 Cont.
<p>The distribution of RMSE<sub>a</sub> over the viewing hemisphere within each BRDF archetype class at an SZA of 45°. (<b>a</b>–<b>i</b>) represent the nine BRDF archetype classes, respectively.</p>
Full article ">Figure 16
<p>The distribution of RMSE<sub>r</sub> over the viewing hemisphere within each BRDF archetype class at an SZA of 45°. (<b>a</b>–<b>i</b>) represent the nine BRDF archetype classes, respectively.</p>
Full article ">Figure 17
<p>Accuracy evaluation of albedo retrieval using LUT based on multi-angular data simulated by PROSAIL. (<b>a</b>–<b>c</b>) refer to the RMSE<sub>r</sub> and (<b>d</b>–<b>f</b>) refer to the RMSE<sub>a</sub> over the viewing hemisphere under SZA of 15°, 45°, and 60°.</p>
Full article ">Figure 18
<p>The comparison between simulated MODIS directional reflectance (<b>a</b>,<b>c</b>,<b>e</b>,<b>f</b>) or albedo (<b>b</b>,<b>d</b>,<b>f</b>,<b>h</b>) retrieved from the BRDF archetype LUTs and MODIS albedo. (<b>a</b>,<b>b</b>,<b>e</b>,<b>f</b>) refer to the red band, and (<b>c</b>,<b>d</b>,<b>g</b>,<b>h</b>) refer to the NIR band. The color represents the density of overlapping points.</p>
Full article ">Figure 19
<p>The comparison between POLDER observations (<b>a</b>,<b>c</b>) or albedo (<b>b</b>,<b>d</b>) retrieved from the BRDF archetype LUTs and POLDER albedo based on multi-angular observations. (<b>a</b>,<b>b</b>) refer to the red band, and (<b>c</b>,<b>d</b>) refer to the NIR band.</p>
Full article ">Figure 20
<p>The comparison between ground observations (<b>a</b>,<b>c</b>) or albedo (<b>b</b>,<b>d</b>) retrieved from the BRDF archetype LUTs and albedo based on multi-angular observations. (<b>a</b>,<b>b</b>) refer to the red band, and (<b>c</b>,<b>d</b>) refer to the NIR band.</p>
Full article ">
21 pages, 3645 KiB  
Article
Evaluating the Performance of the Enhanced Ross-Li Models in Characterizing BRDF/Albedo/NBAR Characteristics for Various Land Cover Types in the POLDER Database
by Anxin Ding, Ziti Jiao, Alexander Kokhanovsky, Xiaoning Zhang, Jing Guo, Ping Zhao, Mingming Zhang, Hailan Jiang and Kaijian Xu
Remote Sens. 2024, 16(12), 2119; https://doi.org/10.3390/rs16122119 - 11 Jun 2024
Viewed by 824
Abstract
The latest versions of the Ross-Li model include kernels that represent isotropic reflection of the surface, describe backward reflection of soil and vegetation systems, characterize strong forward reflection of snow, and adequately consider the hotspot effect (i.e., RossThick-LiSparseReciprocalChen-Snow, RTLSRCS), theoretically able to effectively [...] Read more.
The latest versions of the Ross-Li model include kernels that represent isotropic reflection of the surface, describe backward reflection of soil and vegetation systems, characterize strong forward reflection of snow, and adequately consider the hotspot effect (i.e., RossThick-LiSparseReciprocalChen-Snow, RTLSRCS), theoretically able to effectively characterize BRDF/Albedo/NBAR features for various land surface types. However, a systematic evaluation of the RTLSRCS model is still lacking for various land cover types. In this paper, we conducted a thorough assessment of the RTLSRCS and RossThick-LiSparseReciprocalChen (RTLSRC) models in characterizing BRDF/Albedo/NBAR characteristics by using the global POLDER BRDF database. The primary highlights of this paper include the following: (1) Both models demonstrate high accuracy in characterizing the BRDF characteristics across 16 IGBP types. However, the accuracy of the RTLSRC model is notably reduced for land cover types with high reflectance and strong forward reflection characteristics, such as Snow and Ice (SI), Deciduous Needleleaf Forests (DNF), and Barren or Sparsely Vegetated (BSV). In contrast, the RTLSRCS model shows a significant improvement in accuracy for these land cover types. (2) These two models exhibit highly consistent albedo inversion across various land cover types (R2 > 0.9), particularly in black-sky and blue-sky albedo, except for SI. However, significant differences in white-sky albedo inversion persist between these two models for Evergreen Needleleaf Forests (ENF), Evergreen Broadleaf Forests (EBF), Urban Areas (UA), and SI (p < 0.05). (3) The NBAR values inverted by these two models are nearly identical across the other 15 land cover types. However, the consistency of NBAR results is relatively poor for SI. The RTLSRC model tends to overestimate compared to the RTLSRCS model, with a noticeable bias of approximately 0.024. This study holds significant importance for understanding different versions of Ross-Li models and improving the accuracy of satellite BRDF/Albedo/NBAR products. Full article
(This article belongs to the Special Issue Remote Sensing of Surface BRDF and Albedo)
Show Figures

Figure 1

Figure 1
<p>The global distribution of the selected POLDER pixels, with red and blue points representing pixels from the 2006 and 2008 POLDER BRDF databases, respectively.</p>
Full article ">Figure 2
<p>The fitting accuracy of the global POLDER BRDF database for 16 land cover types using the RTLSRC and RTLSRCS models. The metrics include the R<sup>2</sup> (<b>a</b>), RMSE (<b>b</b>), NRMSE (<b>c</b>) and MRE (<b>d</b>).</p>
Full article ">Figure 3
<p>The average values and standard deviation of BRDF parameters for the RTLSRC and RTLSRCS models for 16 land cover types.</p>
Full article ">Figure 3 Cont.
<p>The average values and standard deviation of BRDF parameters for the RTLSRC and RTLSRCS models for 16 land cover types.</p>
Full article ">Figure 4
<p>The comparison results of the RTLSRC and RTLSRCS models for the inversion of BSA and WSA across 16 land cover types, with the results from the RTLSRCS model serving as reference data. The metrics include the R<sup>2</sup> (<b>a</b>), RMSE (<b>b</b>), bias (<b>c</b>), NRMSE (<b>d</b>) and MRE (<b>e</b>).</p>
Full article ">Figure 5
<p>The T-test (<b>a</b>) and <span class="html-italic">p</span> value (<b>b</b>) results comparing the RTLSRC and RTLSRCS models for the inversion of BSA and WSA across 16 land cover types, with the results from the RTLSRCS model serving as the reference data.</p>
Full article ">Figure 6
<p>The average values and standard deviation of NBAR values for the RTLSRC and RTLSRCS models for 16 land cover types.</p>
Full article ">Figure 7
<p>The variations in the average values and standard deviation of the NBAR values inverted by the RTLSRC and RTLSRCS models vary with the SZA for SI.</p>
Full article ">
26 pages, 10154 KiB  
Article
Retrieval of At-Surface Upwelling Radiance and Albedo by Parameterizing Cloud Scattering and Transmittance over Rugged Terrain
by Junru Jia, Massimo Menenti, Li Jia, Qiting Chen and Anlun Xu
Remote Sens. 2024, 16(10), 1723; https://doi.org/10.3390/rs16101723 - 13 May 2024
Viewed by 1118
Abstract
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 [...] Read more.
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/m2 and 0.0422, 67.6 W/m2 and 0.0545, 98.6 W/m2 and 0.0992, and 78.0 98.6 W/m2 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. Full article
Show Figures

Figure 1

Figure 1
<p>A DEM map of the third pole showing the study area and locations of the ground stations (black dots).</p>
Full article ">Figure 2
<p>Schematic diagram of the workflow for the retrieval of the surface albedo in rugged terrain.</p>
Full article ">Figure 3
<p>Comparison of the monthly distribution of the sum of the number of valid daily surface reflectance observations (ValObs) between our method and MODIS MCD43A2 that can be used for successful retrieval of daily surface albedo within each 16-day moving window at the four stations on the Tibetan Plateau in 2018: (<b>a</b>) Dali Station, (<b>b</b>) MAWORS Station, (<b>c</b>) NAMORS Station, and (<b>d</b>) QOMS Station.</p>
Full article ">Figure 4
<p>Validation of the instantaneous at-surface upwelling estimated using the method developed in this study when CF &lt; 50% by comparison with the ground observations at Dali, NAMORS, MAWORS, and QOMS stations on the Tibetan Plateau in 2018 during the MODIS Terra and Aqua overpass time.</p>
Full article ">Figure 5
<p>The MAPE for at-surface upwelling radiance estimates at different cloud fractions at Dali, MAWORS, NAMORS, and QOMS stations on the Tibetan Plateau in 2018.</p>
Full article ">Figure 6
<p>(<b>a</b>) The estimates of instantaneous at-surface upwelling radiance of the Tibetan Plateau (MYD Swath: UTC 2018-01-02-0715) and (<b>b</b>) MODIS true color composite images of bands 4, 3, and 1 at the same time.</p>
Full article ">Figure 7
<p>Comparison of clear day results of the estimated broadband surface albedo using the new method proposed in this paper, the MODIS shortwave WSA and BSA and the ground measured albedo at Dali, MAWORS, NAMORS, and QOMS stations on the Tibetan Plateau in 2018.</p>
Full article ">Figure 7 Cont.
<p>Comparison of clear day results of the estimated broadband surface albedo using the new method proposed in this paper, the MODIS shortwave WSA and BSA and the ground measured albedo at Dali, MAWORS, NAMORS, and QOMS stations on the Tibetan Plateau in 2018.</p>
Full article ">Figure 8
<p>The MAPE between the estimated surface albedo by our method and the ground measurements of albedo at different cloud fractions at Dali, MAWORS, NAMORS, and QOMS stations in 2018.</p>
Full article ">Figure 9
<p>Spatial distribution over an area in the surrounding of the MAWORS station on the Tibetan Plateau: (<b>a</b>) elevation (m), (<b>b</b>) slope (°), (<b>c</b>) aspect (°), (<b>d</b>) MCD43A3 black-sky albedo, (<b>e</b>) MCD43A3 white-sky albedo, and (<b>f</b>) our retrievals of surface broadband albedo on 2018.01.01.</p>
Full article ">Figure 9 Cont.
<p>Spatial distribution over an area in the surrounding of the MAWORS station on the Tibetan Plateau: (<b>a</b>) elevation (m), (<b>b</b>) slope (°), (<b>c</b>) aspect (°), (<b>d</b>) MCD43A3 black-sky albedo, (<b>e</b>) MCD43A3 white-sky albedo, and (<b>f</b>) our retrievals of surface broadband albedo on 2018.01.01.</p>
Full article ">Figure 10
<p>The time series of the estimated broadband surface albedo by our method, MODIS shortwave WSA and BSA, and ground-measured albedo at the MAWORS site on the Tibetan Plateau in 2018.</p>
Full article ">
16 pages, 8104 KiB  
Technical Note
Diurnal Asymmetry Effects of Photovoltaic Power Plants on Land Surface Temperature in Gobi Deserts
by Xubang Wang, Qianru Zhou, Yong Zhang, Xiang Liu, Jianquan Liu, Shengyun Chen, Xinxin Wang and Jihua Wu
Remote Sens. 2024, 16(10), 1711; https://doi.org/10.3390/rs16101711 - 11 May 2024
Viewed by 1356
Abstract
The global expansion of photovoltaic (PV) power plants, especially in ecologically fragile regions like the Gobi Desert, highlights the suitability of such areas for large-scale PV development. The most direct impact of PV development in the Gobi Desert is temperature change that results [...] Read more.
The global expansion of photovoltaic (PV) power plants, especially in ecologically fragile regions like the Gobi Desert, highlights the suitability of such areas for large-scale PV development. The most direct impact of PV development in the Gobi Desert is temperature change that results from the land-use-induced albedo changes; however, the detailed and systemic understanding of the effects of PV expansion on land surface temperature remains limited. This study focuses on the 16 largest PV plants in the Chinese Gobi Desert, utilizing remote sensing data to assess their effects on land surface temperature. Our result showed a cooling effect during the daytime (−0.69 ± 0.10 °C), but a warming effect during the nighttime (0.23 ± 0.05 °C); the overall effect on the daily mean was a cooling effect (−0.22 ± 0.05 °C). Seasonal variations were observed, with the most significant cooling effect in autumn and the weakest in summer. The PV area was the most significant factor which influenced the temperature variation across PV plants. Our findings enrich our understanding of the environmental effects arising from the construction of PV plants and provide vital information for the design and management of increasingly renewable electricity systems globally. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Study area. (<b>a</b>) Location of the Gobi region and photovoltaic power plants in China; (<b>b</b>) Dunhuang PV power plant; (<b>c</b>) Google Earth satellite imagery of Dunhuang PV power plant. The green line and blue line in (<b>b</b>,<b>c</b>) indicate in (1 km buffer) and out (15 km buffer) of the Dunhuang plant.</p>
Full article ">Figure 2
<p>Workflow of this study. In the top left image, the red rectangle represents the PV panel, the green and blue lines represent the buffer, and the orange color represents the Gobi surface. In the top right image, the colored image represents LST data. In the lower image, ***: statistically significant at <span class="html-italic">p</span> &lt; 0.001 levels.</p>
Full article ">Figure 3
<p>Effects of photovoltaic power plant on LST of (<b>a</b>) daytime period and (<b>b</b>) nighttime period.</p>
Full article ">Figure 4
<p>The PV power plant effects on the annual means of LST of (<b>a</b>) daytime period and (<b>b</b>) nighttime period in all photovoltaic power plants, the black line represents the extent of the PV plant.</p>
Full article ">Figure 5
<p>Effect of photovoltaic power plants on LST across all plants. (<b>a</b>) Diurnal variations in effects (ΔLST); the background violin plot characterizes the distribution of plants in each diurnal period effect, while white dots represent the mean value. Statistical difference was tested by one-sample <span class="html-italic">t</span>-test between each period effect and zero (μ = 0) and independent two-sample <span class="html-italic">t</span>-test between daytime and nighttime period. ***: statistically significant at <span class="html-italic">p</span> &lt; 0.001 levels, respectively. Seasonal variation in effects (ΔLST) separated into (<b>b</b>) daily mean, (<b>c</b>) daytime period, and (<b>d</b>) nighttime period. Statistical difference was tested by Kruskal–Wallis analysis and Dunn’s test as a post hoc analysis to investigate pairwise differences between seasons. The boxes represent the interquartile range, the lines inside the boxes represent the medians, and the whiskers denote the lowest and highest values within 1.5 times the interquartile range. Lowercase letters denote significant differences between seasons. Colored dots represent each plant data.</p>
Full article ">Figure 6
<p>Effects of PV power plant on monthly LSTs across all plants, separated into (<b>a</b>) daily mean and (<b>b</b>) daytime and nighttime period.</p>
Full article ">Figure 7
<p>Factors that influenced the effects of PV power plant on LST, include area, mean annual temperature (MAT), mean annual precipitation (MAP), solar radiation (Rs), wind speed (Ws), and water vapor pressure (Vp). Estimate effect sizes with 95% confidence intervals are derived from the weighted average standardized coefficients of models with ΔAICc &lt; 4. The relative importance of factors on (<b>a</b>) daily mean, (<b>b</b>) daytime period, and (<b>c</b>) nighttime period, as estimated by linear models. Blue lines indicate negative effects, and red lines indicate positive effects. *: statistically significant at <span class="html-italic">p</span> &lt; 0.05 level. Model-averaged importance of the predictors and the <span class="html-italic">p</span>-value of each factor are shown in <a href="#remotesensing-16-01711-f008" class="html-fig">Figure 8</a> and <a href="#app1-remotesensing-16-01711" class="html-app">Table S3</a>.</p>
Full article ">Figure 8
<p>Importance of each predictor of the PV power plant effects, (<b>a</b>) daily mean, (<b>b</b>) daytime period, and (<b>c</b>) nighttime period. The importance value is based on the sum of Akaike weights derived from model selection using corrected Akaike’s information criteria. Cutoff is set at 0.8 (dash line) to differentiate between essential and nonessential predictors.</p>
Full article ">
23 pages, 6750 KiB  
Article
Assessing Satellite Data’s Role in Substituting Ground Measurements for Urban Surfaces Characterization: A Step towards UHI Mitigation
by Davide Parmeggiani, Francesca Despini, Sofia Costanzini, Malvina Silvestri, Federico Rabuffi, Sergio Teggi and Grazia Ghermandi
Atmosphere 2024, 15(5), 551; https://doi.org/10.3390/atmos15050551 - 29 Apr 2024
Viewed by 946
Abstract
Urban surfaces play a crucial role in shaping the Urban Heat Island (UHI) effect by absorbing and retaining significant solar radiation. This paper explores the potential of high-resolution satellite imagery as an alternative method for characterizing urban surfaces to support UHI mitigation strategies [...] Read more.
Urban surfaces play a crucial role in shaping the Urban Heat Island (UHI) effect by absorbing and retaining significant solar radiation. This paper explores the potential of high-resolution satellite imagery as an alternative method for characterizing urban surfaces to support UHI mitigation strategies in urban redevelopment plans. We utilized Landsat images spanning the past 40 years to analyze trends in Land Surface Temperature (LST). Additionally, WorldView-3 (WV3) imagery was acquired for surface characterization, and the results were compared with ground truth measurements using the ASD FieldSpec 4 spectroradiometer. Our findings revealed a strong correlation between satellite-derived surface reflectance and ground truth measurements across various urban surfaces, with Root Mean Square Error (RMSE) values ranging from 0.01 to 0.14. Optimal characterization was observed for surfaces such as bituminous membranes and parking with cobblestones (RMSE < 0.03), although higher RMSE values were noted for tiled roofs, likely due to aging effects. Regarding surface albedo, the differences between satellite-derived data and ground measurements consistently remained below 12% for all surfaces, with the lowest values observed in high heat-absorbing surfaces like bituminous membranes. Despite challenges on certain surfaces, our study highlights the reliability of satellite-derived data for urban surface characterization, thus providing valuable support for UHI mitigation efforts. Full article
(This article belongs to the Special Issue UHI Analysis and Evaluation with Remote Sensing Data)
Show Figures

Figure 1

Figure 1
<p>Location of the study area.</p>
Full article ">Figure 2
<p>WorldView3 image of the study area and overlook of the chosen ROIs highlighted in yellow. 1: Polyolefin roof, 2: new tiles roof, 3: aged tiles roof, 4: asphalt parking, 5: parking with cobblestones, 6: bituminous membrane.</p>
Full article ">Figure 3
<p>Methodology chart.</p>
Full article ">Figure 4
<p>ASD FieldSpec 4 acquisition. (<b>a</b>) Parking with cobblestones, (<b>b</b>) polyolefin roof, (<b>c</b>) aged tiles roof, (<b>d</b>) asphalt parking, (<b>e</b>) new tiles roof, (<b>f</b>) bituminous membrane.</p>
Full article ">Figure 5
<p>Median, maximum, and minimum LST values for summer period from 1985 to 2023, and number of available summer images.</p>
Full article ">Figure 6
<p>WV3 image (<b>left</b>), LST map [°C] (<b>center</b>), and albedo map (<b>right</b>) of the study area. Boxes highlighted and depicted in the figure with zoom are used to visualize critical areas.</p>
Full article ">Figure 7
<p>Comparison between the spectral signatures of bituminous membrane achieved by the ASD Fieldspec 4, the Jasco Spectrophotometer, and the resampled WV3 imagery.</p>
Full article ">Figure 8
<p>Comparative spectral signature of the six acquired ROI with WV3 Surface Reflectance and ASD Fieldspec 4 measurements resampled with WV3 spectrum.</p>
Full article ">
20 pages, 9422 KiB  
Article
Impact of Wildfires on Land Surface Cold Season Climate in the Northern High-Latitudes: A Study on Changes in Vegetation, Snow Dynamics, Albedo, and Radiative Forcing
by Melissa Linares and Wenge Ni-Meister
Remote Sens. 2024, 16(8), 1461; https://doi.org/10.3390/rs16081461 - 20 Apr 2024
Viewed by 1336
Abstract
Anthropogenic climate change is increasing the occurrence of wildfires, especially in northern high latitudes, leading to a shift in land surface climate. This study aims to determine the predominant climatic effects of fires in boreal forests to assess their impact on vegetation composition, [...] Read more.
Anthropogenic climate change is increasing the occurrence of wildfires, especially in northern high latitudes, leading to a shift in land surface climate. This study aims to determine the predominant climatic effects of fires in boreal forests to assess their impact on vegetation composition, surface albedo, and snow dynamics. The influence of fire-induced changes on Earth’s radiative forcing is investigated, while considering variations in burn severity and postfire vegetation structure. Six burn sites are explored in central Alaska’s boreal region, alongside six control sites, by utilizing Moderate Resolution Imaging Spectroradiometer (MODIS)-derived albedo, Leaf Area Index (LAI), snowmelt timing data, AmeriFlux radiation, National Land Cover Database (NLCD) land cover, and Monitoring Trends in Burn Severity (MTBS) data. Key findings reveal significant postfire shifts in land cover at each site, mainly from high- to low-stature vegetation. A continuous increase in postfire surface albedo and negative surface shortwave forcing was noted even after 12 years postfire, particularly during the spring and at high-severity burn areas. Results indicate that the cooling effect from increased albedo during the snow season may surpass the warming effects of earlier snowmelt. The overall climate impact of fires depends on burn severity and vegetation composition. Full article
(This article belongs to the Special Issue Remote Sensing of Solar Radiation Absorbed by Land Surfaces)
Show Figures

Figure 1

Figure 1
<p>The study sites are in the center of Alaska, covering six burn sites (red) and their corresponding control sites (blue).</p>
Full article ">Figure 2
<p>Flowchart of Data Processing and Analysis.</p>
Full article ">Figure 3
<p>Column one displays burn severity percent distributions for each site, indicating the extent of fire impact. Column two features corresponding burn severity maps visually representing the distribution of burn severity classes. Columns three and four contrast the NLCD 2001 prefire and NLCD 2016 postfire land cover maps, illustrating changes in vegetation over time. The fifth column provides the NLCD 2001 land cover maps for control sites.</p>
Full article ">Figure 4
<p>2001 and 2016 NLCD Land Cover Percent for Burn Sites. The stacked bars show the percentage of different land cover types at each burn site location in 2001 and 2016, derived from the National Land Cover Database (NLCD). Land cover categories include evergreen forest, deciduous forest, shrub/scrub, emergent herbaceous wetlands, woody wetlands, grassland/herbaceous, open water, and developed/low-intensity areas. The total absolute change percentage and vegetation cover density classifications (dense vs. sparse) are provided for each burn site between the two time periods. This allows for visualizing the impact of fires on shifting land cover compositions and forest density at these locations over 15 years.</p>
Full article ">Figure 5
<p>The summer months (June, July, and August) and Winter Months (January, February, and March) LAI values for all study sites and control sites over two decades. The red dashed line represents the wildfires.</p>
Full article ">Figure 6
<p>The mean difference in snowmelt dates (black solid) between each site and its control site from 2001 to 2018. Trendlines, represented in blue, depicts the prefire trajectory of changes in snowmelt timing, while the red trendline indicates the general postfire trajectory of these changes. The slope of each trendline is provided in the legend for reference.</p>
Full article ">Figure 7
<p>Time series of albedo for different burn severity classes across all burn sites from 2007 to 2021. The black dotted line marks the fire.</p>
Full article ">Figure 8
<p>Monthly mean differences (black solid) in albedo between fire and control sites. Pre- (dotted blue line) and postfire (dotted red line) trendlines indicate the general trajectory of changes in albedo difference over the time series, with the slope of each trendline provided in the legend.</p>
Full article ">Figure 9
<p>Daily mean Surface Shortwave Forcing (SSF) across various burn severity classes spanning from 2010 to 2017.</p>
Full article ">
17 pages, 1035 KiB  
Article
Modification and Validation of the Soil–Snow Module in the INM RAS Climate Model
by Alexey Chernenkov, Evgeny Volodin, Sergey Kostrykin, Maria Tarasevich and Vasilisa Vorobyeva
Atmosphere 2024, 15(4), 422; https://doi.org/10.3390/atmos15040422 - 29 Mar 2024
Cited by 1 | Viewed by 1240
Abstract
This paper describes the modification of a simple land snow cover module of the INM RAS climate model. The possible liquid water and refreezing of meltwater in the snow layer are taken into account by the proposed parameterization. This is particularly important for [...] Read more.
This paper describes the modification of a simple land snow cover module of the INM RAS climate model. The possible liquid water and refreezing of meltwater in the snow layer are taken into account by the proposed parameterization. This is particularly important for modelling the transition season, as this phenomenon is mainly observed during the formation and melting of the snow cover when the surface temperature fluctuates around 0 °C. The snow density evolution simulation is also added. This parameterization is implemented in the INM-CM snow module and verified on observation data using the ESM-SnowMIP-like protocol. As a result, the INM-CM mean climate snow melt periods are refined, particularly in middle and high latitudes. The snow-covered area according to the model is also improved. In the future, a modified version of the land snow module can be used, coupled with a snow albedo model that takes into account snow metamorphism. This module can also be applied to sea ice snow. Full article
(This article belongs to the Section Atmospheric Techniques, Instruments, and Modeling)
Show Figures

Figure 1

Figure 1
<p>Conceptual flowchart of the modified snow model.</p>
Full article ">Figure 2
<p>Monthly average depth (<b>a</b>) and mass (<b>b</b>) of snow cover during the entire simulation period for the Col de Porte site, according to the INM-CM snow module (blue—original version, red—modified) and observations (black).</p>
Full article ">Figure 3
<p>Annual variation of snow depths averaged over the simulation period, according to the INM-CM snow module (blue—original version, red—modified) and site observations (black).</p>
Full article ">Figure 4
<p>Annual variation of snow mass averaged over the simulation period, according to the INM-CM snow module (blue—original version, red—modified) and site observations (black).</p>
Full article ">Figure 5
<p>Monthly average surface temperatures for Col de Porte (<b>a</b>), Swamp Angel (<b>b</b>), and Weissfluhjoch (<b>c</b>) sites during the entire simulation period, according to the INM-CM snow module (blue—original version, red—modified) and observations (black).</p>
Full article ">Figure 6
<p>Quality metrics (normalized root-mean-square errors, bias, and correlation coefficient) for snow depth, snow water equivalent thickness, and surface temperature.</p>
Full article ">Figure 7
<p>Monthly average snow depth during the simulation period for the Col de Porte site, according to the different INM-CM snow module versions (including the version with the correction of the heat transfer coefficient, <math display="inline"><semantics> <msub> <mi>C</mi> <mi>T</mi> </msub> </semantics></math>) and observation data.</p>
Full article ">Figure 8
<p>Quality metrics (normalized root-mean-square errors, bias, and correlation coefficient) for snow depth and SWE for mountain sites (for the original, modified, and modified with the additional tuning of the INM-CM snow module versions).</p>
Full article ">Figure 9
<p>Monthly average soil temperature at a depth of 10 cm during the simulation period for the Reynolds Mountain East site, according to the INM-CM snow module (blue—original version, red—modified) and observations (black).</p>
Full article ">Figure 10
<p>Quality metrics (normalized root-mean-square errors, bias, and correlation coefficient) for the temperature in the upper soil layers (at a depth of ∼10 cm) according to the different INM-CM snow module versions.</p>
Full article ">Figure 11
<p>(<b>a</b>) Mean monthly SWE according to the INM RAS global climate model with the modified version of the soil–snow module, (<b>b</b>) difference with the original version (INM-CM48), and (<b>c</b>) with the NSIDC-Blended5 reanalysis data (in [mm]; model data averaged over the ensemble and all simulation period; month—May).</p>
Full article ">Figure 12
<p>Annual snow-covered area variation in the Northern Hemisphere (averaged over the ensemble and for 2003–2017) according to the original and modified versions of INM-CM and CAMS reanalysis data.</p>
Full article ">Figure 13
<p>Annual snow-covered area variation in the Northern Hemisphere (averaged over the ensemble and for 1998–2010) according to the original and modified versions of INM-CM and the NSIDC-Blended5 reanalysis data.</p>
Full article ">
12 pages, 6008 KiB  
Article
The Influence of Vegetation on Climate Elements in Northwestern China
by Bicheng Huang, Yu Huang, Dan Wu, Xinyue Bao, Yongping Wu, Guolin Feng and Li Li
Atmosphere 2024, 15(3), 325; https://doi.org/10.3390/atmos15030325 - 5 Mar 2024
Cited by 1 | Viewed by 1520
Abstract
Vegetation plays a crucial role in maintaining the balance between nature, water and soil resources. However, understanding its impact mechanisms in arid and semi-arid areas remains limited. This study aims to analyze the spatial–temporal characteristics of the vegetation leaf area index (LAI) and [...] Read more.
Vegetation plays a crucial role in maintaining the balance between nature, water and soil resources. However, understanding its impact mechanisms in arid and semi-arid areas remains limited. This study aims to analyze the spatial–temporal characteristics of the vegetation leaf area index (LAI) and climate elements in typical regions of northwest China and the correlations between LAI and climate elements; it also aims to explore the influence of regional vegetation growth on climate change. The results reveal significant correlations between LAI and various climate elements. Specifically, within the same region, surface temperature, precipitation, vegetation transpiration, and total evaporation show positive correlations with the LAI, whereas surface albedo shows a negative correlation. Vegetation may affect climate through both heat and water exchange between the land and atmosphere. Increased vegetation leads to the enhanced absorption of solar radiation by the land surface, elevating surface temperature. Increased levels of vegetation also increase vegetation transpiration and total evaporation, increasing the water vapor content in the atmosphere and thus leading to increased surface precipitation. Therefore, vegetation distribution plays a role in climate change, and ecological restoration projects in the northwest region hold significant potential for addressing ecological challenges in its arid and semi-arid areas. Full article
(This article belongs to the Section Climatology)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Geographical location and (<b>b</b>) elevation of the study area.</p>
Full article ">Figure 2
<p>The average annual distribution of (<b>a</b>) high vegetation coverage, (<b>b</b>) low vegetation coverage, (<b>c</b>) high vegetation leaf area indices and (<b>d</b>) low vegetation leaf area indices in the northwest region of China during 1979–2020. Region I, II, III represent the typical areas with highest vegetation coverage, relatively high vegetation coverage, and with low vegetation coverage, respectively.</p>
Full article ">Figure 3
<p>Temporal distribution of surface temperature in Northwest China from 1979 to 2020 (<b>a</b>), spatial distribution of surface temperature in Northwest China (<b>b</b>), standardized anomalous surface temperature in regions I, II and III (<b>c</b>), and temporal distribution of surface temperatures in regions I, II, and III (<b>d</b>–<b>f</b>). Black lines indicate trends and red dotted lines indicate low-pass filtering.</p>
Full article ">Figure 4
<p>Temporal distribution of albedo in Northwest China from 1979 to 2020 (<b>a</b>), spatial distribution of albedo in Northwest China (<b>b</b>), standardized anomalous albedo in regions I, II and III (<b>c</b>), and temporal distribution of albedo in regions I, II, and III (<b>d</b>–<b>f</b>). Black lines indicate trends and red dotted lines indicate low-pass filtering.</p>
Full article ">Figure 5
<p>Temporal distribution of precipitation in Northwest China from 1979 to 2020 (<b>a</b>), spatial distribution of precipitationin Northwest China (<b>b</b>), standardized anomalous precipitation in regions I, II and III (<b>c</b>), and temporal distribution of precipitation in regions I, II, and III (<b>d</b>–<b>f</b>). Black lines indicate trends and red dotted lines indicate low-pass filtering.</p>
Full article ">Figure 6
<p>Temporal distribution of total evaporation in Northwest China from 1979 to 2020 (<b>a</b>), spatial distribution of total evaporation in Northwest China (<b>b</b>), standardized anomalous total evaporation in regions I, II and III (<b>c</b>), and temporal distribution of total evaporation in regions I, II, and III (<b>d</b>–<b>f</b>). Black lines indicate trends and red dotted lines indicate low-pass filtering.</p>
Full article ">Figure 7
<p>Relationship between LAI and surface temperature (<b>a</b>,<b>e</b>), albedo (<b>b</b>,<b>f</b>), surface precipitation (<b>c</b>,<b>g</b>), and total evaporation (<b>d</b>,<b>h</b>) in Northwest China from 1979 to 2020. ((<b>a</b>–<b>d</b>) for region I, (<b>e</b>–<b>h</b>) for region II). The blue dots represent the values of corresponding climate factors for the low vegetation leaf area index in the same months each year.</p>
Full article ">
24 pages, 25102 KiB  
Article
Modeling and Estimating the Land Surface Temperature (LST) Using Remote Sensing and Machine Learning (Case Study: Yazd, Iran)
by Mohammad Mansourmoghaddam, Iman Rousta, Hamidreza Ghafarian Malamiri, Mostafa Sadeghnejad, Jaromir Krzyszczak and Carla Sofia Santos Ferreira
Remote Sens. 2024, 16(3), 454; https://doi.org/10.3390/rs16030454 - 24 Jan 2024
Cited by 7 | Viewed by 4047
Abstract
The pressing issue of global warming is particularly evident in urban areas, where urban thermal islands amplify the warming effect. Understanding land surface temperature (LST) changes is crucial in mitigating and adapting to the effect of urban heat islands, and ultimately addressing the [...] Read more.
The pressing issue of global warming is particularly evident in urban areas, where urban thermal islands amplify the warming effect. Understanding land surface temperature (LST) changes is crucial in mitigating and adapting to the effect of urban heat islands, and ultimately addressing the broader challenge of global warming. This study estimates LST in the city of Yazd, Iran, where field and high-resolution thermal image data are scarce. LST is assessed through surface parameters (indices) available from Landsat-8 satellite images for two contrasting seasons—winter and summer of 2019 and 2020, and then it is estimated for 2021. The LST is modeled using six machine learning algorithms implemented in R software (version 4.0.2). The accuracy of the models is measured using root mean square error (RMSE), mean absolute error (MAE), root mean square logarithmic error (RMSLE), and mean and standard deviation of the different performance indicators. The results show that the gradient boosting model (GBM) machine learning algorithm is the most accurate in estimating LST. The albedo and NDVI are the surface features with the greatest impact on LST for both the summer (with 80.3% and 11.27% of importance) and winter (with 72.74% and 17.21% of importance). The estimated LST for 2021 showed acceptable accuracy for both seasons. The GBM models for each of the seasons are useful for modeling and estimating the LST based on surface parameters using machine learning, and to support decision-making related to spatial variations in urban surface temperatures. The method developed can help to better understand the urban heat island effect and ultimately support mitigation strategies to improve human well-being and enhance resilience to climate change. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The digital elevation model (DEM) and elevation contours of the study area (Yazd County) with cities outlined together with intercity (grade 1) and urban (grade 2) road map (<b>a</b>), the location of Yazd County within the Yazd province (<b>b</b>), and the Yazd province location in Iran (<b>c</b>).</p>
Full article ">Figure 2
<p>Method used for determination of the true Euclidean distance (adapted from [<a href="#B79-remotesensing-16-00454" class="html-bibr">79</a>]).</p>
Full article ">Figure 3
<p>The LST maps for the two seasons (winter in the <b>left</b> panels and summer in the <b>right</b> panels) in 2019 (<b>top</b> panels) and 2020 (<b>bottom</b> panels).</p>
Full article ">Figure 4
<p>The LST images of the 2021 winter (<b>left</b>) and summer (<b>right</b>) based on model predictions.</p>
Full article ">Figure 5
<p>Spatial variation in calculated ground surface features: (<b>a</b>) albedo, (<b>b</b>) NDVI, (<b>c</b>) NDBI, (<b>d</b>) NDBaI, (<b>e</b>) distance (m) from water bodies for each pixel, and (<b>f</b>) distance from mountains (top) and grade 1 and 2 roads (bottom) for each pixel map. Each figure includes maps for 2019 (top), 2020 (middle), and 2021 (bottom) for the winter (left panels) and summer (right panels) seasons.</p>
Full article ">Figure 5 Cont.
<p>Spatial variation in calculated ground surface features: (<b>a</b>) albedo, (<b>b</b>) NDVI, (<b>c</b>) NDBI, (<b>d</b>) NDBaI, (<b>e</b>) distance (m) from water bodies for each pixel, and (<b>f</b>) distance from mountains (top) and grade 1 and 2 roads (bottom) for each pixel map. Each figure includes maps for 2019 (top), 2020 (middle), and 2021 (bottom) for the winter (left panels) and summer (right panels) seasons.</p>
Full article ">Figure 6
<p>Shapley Additive exPlanation (SHAP) summary plot.</p>
Full article ">Figure 7
<p>The map of the estimated LST for winter 2021 (<b>top</b>) and summer (<b>down</b>) 2021 for the Yazd cities (<b>left</b>) and its distance from the 45° polynomial line (DFPL) image (<b>right</b>).</p>
Full article ">Figure 7 Cont.
<p>The map of the estimated LST for winter 2021 (<b>top</b>) and summer (<b>down</b>) 2021 for the Yazd cities (<b>left</b>) and its distance from the 45° polynomial line (DFPL) image (<b>right</b>).</p>
Full article ">Figure 8
<p>Scatterplot for estimated and actual LST values (<b>left</b>), and histogram for actual (LST) and estimated (LSTpre) LST values (<b>right</b>) for winter (<b>top</b>) and summer (<b>down</b>) 2021.</p>
Full article ">Figure 8 Cont.
<p>Scatterplot for estimated and actual LST values (<b>left</b>), and histogram for actual (LST) and estimated (LSTpre) LST values (<b>right</b>) for winter (<b>top</b>) and summer (<b>down</b>) 2021.</p>
Full article ">
25 pages, 83192 KiB  
Article
Preliminary Retrieval and Validation of Aerosol Optical Depths from FY-4B Advanced Geostationary Radiation Imager Images
by Dong Zhou, Qingxin Wang, Siwei Li and Jie Yang
Remote Sens. 2024, 16(2), 372; https://doi.org/10.3390/rs16020372 - 17 Jan 2024
Cited by 1 | Viewed by 1777
Abstract
Fengyun-4B (FY-4B) is the latest Chinese next-generation geostationary meteorological satellite. The Advanced Geostationary Radiation Imager (AGRI) aboard FY-4B is equipped with 15 spectral bands, from visible to infrared, suitable for aerosol optical depth (AOD) retrieval. In this study, an overland AOD retrieval algorithm [...] Read more.
Fengyun-4B (FY-4B) is the latest Chinese next-generation geostationary meteorological satellite. The Advanced Geostationary Radiation Imager (AGRI) aboard FY-4B is equipped with 15 spectral bands, from visible to infrared, suitable for aerosol optical depth (AOD) retrieval. In this study, an overland AOD retrieval algorithm was developed for the FY-4B AGRI. Considering the large directional variation in the FY-4B AGRI reflectances, a bidirectional reflectance distribution function (BRDF) database was built, through which to estimate land surface reflectance/albedo. Seasonal aerosol models, based on four geographical regions in China, were developed between 2016 and 2022 using AERONET aerosol products, to improve their applicability to regional distribution differences and seasonal variations in aerosol types. AGRI AODs were retrieved using this new method over China from September 2022 to August 2023 and validated against ground-based measurements. The AGRI, Advanced Himawari Imager (AHI), and Moderate-Resolution Imaging Spectroradiometer (MODIS) official land aerosol products were also evaluated for comparison purposes. The results showed that the AGRI AOD retrievals were highly consistent with the AERONET AOD measurements, with a correlation coefficient (R) of 0.88, root mean square error (RMSE) of 0.14, and proportion that met an expected error (EE) of 65.04%. Intercomparisons between the AGRI AOD and other operational AOD products showed that the AGRI AOD retrievals achieved better performance results than the AGRI, AHI, and MODIS official AOD products. Moreover, the AGRI AOD retrievals showed high spatial integrity and stable performance at different times and regions, as well as under different aerosol loadings and characteristics. These results demonstrate the robustness of the new aerosol retrieval method and the potential of FY-4B AGRI measurements for the monitoring of aerosols with high accuracy and temporal resolutions. Full article
(This article belongs to the Section Atmospheric Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Study area and the distribution of the AERONET sites for aerosol model building.</p>
Full article ">Figure 2
<p>Flowchart of the FY-4B AGRI AOD retrieval algorithm.</p>
Full article ">Figure 3
<p>Four geographical regions of China.</p>
Full article ">Figure 4
<p>The aerosol optical models of average results for four seasons in QTR, NWR, NR, and SR. (<b>a</b>,<b>f</b>,<b>k</b>,<b>p</b>) are extinction coefficients; (<b>b</b>,<b>g</b>,<b>l</b>,<b>q</b>) are phase functions at 440 nm; (<b>c</b>,<b>h</b>,<b>m</b>,<b>r</b>) are size distributions; (<b>d</b>,<b>i</b>,<b>n</b>,<b>s</b>) are single scattering albedo; (<b>e</b>,<b>j</b>,<b>o</b>,<b>t</b>) are AE indices at 440–675 nm.</p>
Full article ">Figure 5
<p>(<b>a</b>) Spectral response function between the MODIS and AGRI blue wavebands. (<b>b</b>) Scatter plots of different surface features of the MODIS and AGRI. The black line is the 1:1 line.</p>
Full article ">Figure 6
<p>Satellite remote sensing and ground verification scatter plots. AERONET AOD versus (<b>a</b>) AGRI AOD, (<b>b</b>) FY4B_LDA AOD, (<b>c</b>) AHI AOD, (<b>d</b>) MODIS DB AOD, (<b>e</b>) MODIS DT AOD at 550 nm. The black dashed lines and red solid lines are the EE lines and regression lines, respectively.</p>
Full article ">Figure 7
<p>Accuracy evaluation of AGRI, FY4B_LDA, AHI, MODIS DB, and MODIS DT AOD retrievals against ground-based measurements for each site. (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>,<b>m</b>): correlation coefficient (R); (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>,<b>n</b>): root mean square error (RMSE); (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>,<b>o</b>): percentage of AOD within the expected error envelopes (%).</p>
Full article ">Figure 7 Cont.
<p>Accuracy evaluation of AGRI, FY4B_LDA, AHI, MODIS DB, and MODIS DT AOD retrievals against ground-based measurements for each site. (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>,<b>m</b>): correlation coefficient (R); (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>,<b>n</b>): root mean square error (RMSE); (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>,<b>o</b>): percentage of AOD within the expected error envelopes (%).</p>
Full article ">Figure 8
<p>Satellite remote sensing and ground verification scatter plots in different seasons. AERONET AOD versus (<b>a</b>–<b>d</b>) AGRI AOD, (<b>e</b>–<b>h</b>) FY4B_LDA AOD, (<b>i</b>–<b>l</b>) AHI AOD, (<b>m</b>–<b>p</b>) MODIS DB AOD, (<b>q</b>–<b>t</b>) MODIS DT AOD at 550 nm.</p>
Full article ">Figure 8 Cont.
<p>Satellite remote sensing and ground verification scatter plots in different seasons. AERONET AOD versus (<b>a</b>–<b>d</b>) AGRI AOD, (<b>e</b>–<b>h</b>) FY4B_LDA AOD, (<b>i</b>–<b>l</b>) AHI AOD, (<b>m</b>–<b>p</b>) MODIS DB AOD, (<b>q</b>–<b>t</b>) MODIS DT AOD at 550 nm.</p>
Full article ">Figure 9
<p>Accuracy evaluation of AGRI, FY4B_LDA, AHI, MODIS DB, and MODIS DT AOD retrievals against ground-based measurements in different months. (<b>a</b>): correlation coefficient (R); (<b>b</b>): root mean square error (RMSE); (<b>c</b>): mean absolute error (MAE); (<b>d</b>): mean relative error (MRE); (<b>e</b>): percentage of AOD within the expected error envelopes (%); (<b>f</b>): relative mean deviation (RMB);.</p>
Full article ">Figure 10
<p>AOD bias (i.e., AGRI AOD—AERONET AOD) in terms of (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>,<b>m</b>) AERONET 550 nm AOD, (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>,<b>n</b>) AERONET 440–675 nm AE, (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>,<b>o</b>) AERONET 440 nm SSA. The black dots and error bars represent the median and standard deviation of the AOD bias.</p>
Full article ">Figure 10 Cont.
<p>AOD bias (i.e., AGRI AOD—AERONET AOD) in terms of (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>,<b>m</b>) AERONET 550 nm AOD, (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>,<b>n</b>) AERONET 440–675 nm AE, (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>,<b>o</b>) AERONET 440 nm SSA. The black dots and error bars represent the median and standard deviation of the AOD bias.</p>
Full article ">Figure 11
<p>Maps of the AGRI AOD on 29 September 2022, from 01:00 UTC to 08:00 UTC.</p>
Full article ">Figure 12
<p>Maps of the FY4B_LDA AOD on 29 September 2022, from 01:00 UTC to 08:00 UTC.</p>
Full article ">Figure 13
<p>Maps of the AHI AOD on 29 September 2022, from 01:00 UTC to 08:00 UTC.</p>
Full article ">Figure 14
<p>Maps of the MODIS DB AOD and MODIS DT AOD on 29 September 2022.</p>
Full article ">Figure 15
<p>Averaged AOD distributions for AGRI, FY4B_LDA, AHI, MODIS DB, and MODIS DT from September 2022 to August 2023.</p>
Full article ">
23 pages, 24014 KiB  
Article
Characterisation of Morphological Patterns for Land Surface Temperature Distribution in Urban Environments: An Approach to Identify Priority Areas
by Karina Angélica García-Pardo, David Moreno-Rangel, Samuel Domínguez-Amarillo and José Roberto García-Chávez
Climate 2024, 12(1), 4; https://doi.org/10.3390/cli12010004 - 28 Dec 2023
Viewed by 2154
Abstract
The validated influence of urban biophysical structure on environmental processes within urban areas has heightened the emphasis on studies examining morphological patterns to determine precise locations and underlying causes of urban climate conditions. The present study aims to characterise morphological patterns describing the [...] Read more.
The validated influence of urban biophysical structure on environmental processes within urban areas has heightened the emphasis on studies examining morphological patterns to determine precise locations and underlying causes of urban climate conditions. The present study aims to characterise morphological patterns describing the distribution of Land Surface Temperature (LST) based on a prior classification of biophysical variables, including urban density (building intensity and average height), surface characteristics, shortwave solar radiation (broadband albedo), and seasonal variations in vegetation cover (high, medium, and low levels), retrieved from multisource datasets. To describe the distribution of LST, the variables were calculated, classified, and subsequently, analysed individually and collectively concerning winter and summer LST values applied in an urban neighbourhood in Madrid, Spain. The results from the analytical approaches (observation, correlations, and multiple regressions) were compared to define the morphological patterns. The selection of areas resulting from the morphological patterns with the most unfavourable LST values showed agreement of up to 89% in summer and up to 70% for winter, demonstrating the feasibility of the methods applied to identify priority areas for intervention by season. Notably, low and high vegetation levels emerged as pivotal biophysical characteristics influencing LST distribution compared to the other characteristics, emphasising the significance of integrating detailed seasonal vegetation variations in urban analyses. Full article
(This article belongs to the Section Climate Change and Urban Ecosystems)
Show Figures

Figure 1

Figure 1
<p>Methods and data applied in the study [<a href="#B22-climate-12-00004" class="html-bibr">22</a>].</p>
Full article ">Figure 2
<p>Integrated classification of biophysical characteristics: (<b>i</b>) winter and (<b>ii</b>) summer.</p>
Full article ">Figure 3
<p>Classification of LST values into groups according to ±1 K variation. (<b>i</b>) LSTw; (<b>ii</b>) LSTs.</p>
Full article ">Figure 4
<p>Observation analysis by LSTw groups and classified polygons.</p>
Full article ">Figure 5
<p>Observation analysis by LSTs groups and classified polygons.</p>
Full article ">Figure 6
<p>Scatter plots from simple linear regression of variables with significant relationships with LSTw_Winter: (<b>i</b>) AvgH; (<b>ii</b>) VCw; (<b>iii</b>) VLCwL; and (<b>iv</b>) VLCwM.</p>
Full article ">Figure 7
<p>Observation of unusual residuals (marked in red) from variables most related to LSTw: (<b>a</b>) map of unusual polygons/areas identified in the case study according to residual polygons (indicated in red); (<b>i</b>) residual plot from VLCwL; (<b>ii</b>) residual plot from AvgH; (<b>iii</b>) residual plot from VLCwM; and (<b>iv</b>) residual plot from VCw.</p>
Full article ">Figure 8
<p>Scatter plots from simple linear regression of variables with significant relationships with LSTs_Summer: (<b>i</b>) FSI; (<b>ii</b>) AvgH; (<b>iii</b>) VCs; (<b>iv</b>) VLCsH; (<b>v</b>)VLCsM; (<b>vi</b>) VLCsL; and (<b>vii</b>) BBAlvar_s.</p>
Full article ">Figure 9
<p>Observation of unusual residuals (marked in red) from variables most related to LSTs: (<b>a</b>) map of unusual polygons/areas identified in the case study according to residual polygons (indicated in red); (<b>i</b>) residual plot from VLCsH; (<b>ii</b>) residual plot from VCs; (<b>iii</b>) residual plot from AvgH; (<b>iv</b>) residual plot from VLCsM; (<b>v</b>) residual plot from BBAlvar_s; and (<b>vi</b>) residual plot from FSI.</p>
Full article ">Figure 10
<p>Pearson product–moment correlations to address multicollinearity concerns between variables.</p>
Full article ">Figure 11
<p>Priority areas identified from the selection by filter of morphological patterns for Winter.</p>
Full article ">Figure 12
<p>Priority areas identified from the selection by filter of morphological patterns for Summer.</p>
Full article ">Figure 13
<p>Comparison between a sample of potentially better adapted areas to winter and summer conditions compared to a sample of priority areas identified from morphological patterns.</p>
Full article ">
20 pages, 1113 KiB  
Article
MPPT Strategy of Waterborne Bifacial Photovoltaic Power Generation System Based on Economic Model Predictive Control
by Minan Tang, Jinping Li, Jiandong Qiu, Xi Guo, Bo An, Yaqi Zhang and Wenjuan Wang
Energies 2024, 17(1), 152; https://doi.org/10.3390/en17010152 - 27 Dec 2023
Cited by 1 | Viewed by 1009
Abstract
At present, the new energy industry represented by photovoltaics has become the main force to realize the optimization of China’s energy structure and the goal of “double carbon”; with the absence of land resources, the waterborne bifacial photovoltaic has ushered in a new [...] Read more.
At present, the new energy industry represented by photovoltaics has become the main force to realize the optimization of China’s energy structure and the goal of “double carbon”; with the absence of land resources, the waterborne bifacial photovoltaic has ushered in a new opportunity. Therefore, in order to address the problem that the maximum power point tracking (MPPT) of photovoltaics (PV) could not take into account, the dynamic economic performance in the control process, an economic model predictive control (EMPC), is proposed in this work to realize the MPPT of the waterborne bifacial PV power generation system. Firstly, the model of the bifacial PV module is constructed by combining the ray-tracing irradiance model and considering the effect of water surface albedo on the irradiance absorbed by the module. Secondly, the EMPC controller is designed based on the state-space model of the system to maximize the power generation as the economic performance index, and to solve the optimal input variables time by time to achieve a rolling optimization with the operational requirements of the system itself as the constraints. Thirdly, the MATLAB/Simulink (R2022a) simulation experimental results verify that the EMPC strategy could be utilized to achieve MPPT of the waterborne bifacial PV power generation system, according to the changes of environment. Finally, it is also demonstrated that the bifacial PV power generation system that employed the EMPC strategy outperformed the traditional MPPT algorithm, with respect to both output power tracking velocity and accuracy, and the power generation could be improved by about 6% to 14.5%, which significantly enhances the system’s dynamic process economics. Full article
(This article belongs to the Section A2: Solar Energy and Photovoltaic Systems)
Show Figures

Figure 1

Figure 1
<p>The schematic diagram of the bifacial photovoltaics (PV)generation system.</p>
Full article ">Figure 2
<p>The equivalent circuit model of the PV cell.</p>
Full article ">Figure 3
<p>The relationship between the the solar altitude angle and the water surface albedo.</p>
Full article ">Figure 4
<p>The water surface albedo curves for rainy and sunny days.</p>
Full article ">Figure 5
<p>Schematic of the irradiance reception of the waterborne bifacial PV module.</p>
Full article ">Figure 6
<p>Diagram of the difference between EMPC and MPC.</p>
Full article ">Figure 7
<p>Schematic illustration of maximum power point tracking (MPPT) for the waterborne bifacial PV generation system.</p>
Full article ">Figure 8
<p>The diagram of the economic model predictive control (EMPC) structure.</p>
Full article ">Figure 9
<p>The flow chart of the EMPC algorithm.</p>
Full article ">Figure 10
<p>Factors impacting the power generation performance of photovoltaic modules: (<b>a</b>) front irradiance, (<b>b</b>) rear irradiance, (<b>c</b>) ambient temperature.</p>
Full article ">Figure 11
<p>Dynamic response waveforms of waterborne bifacial PV under varying solar irradiance: (<b>a</b>) output power, (<b>b</b>) output current, (<b>c</b>) output voltage.</p>
Full article ">Figure 12
<p>Dynamic response waveforms of waterborne bifacial PV under varying ambient temperature: (<b>a</b>) output power, (<b>b</b>) output current, (<b>c</b>) output voltage.</p>
Full article ">
25 pages, 33171 KiB  
Article
Spatial Estimation of Snow Water Equivalent for Glaciers and Seasonal Snow in Iceland Using Remote Sensing Snow Cover and Albedo
by Andri Gunnarsson and Sigurdur M. Gardarsson
Hydrology 2024, 11(1), 3; https://doi.org/10.3390/hydrology11010003 - 26 Dec 2023
Viewed by 2762
Abstract
Efficient water resource management in glacier- and snow-dominated basins requires accurate estimates of the snow water equivalent (SWE) in late winter and spring and melt onset timing and intensity. To understand the high spatio-temporal variability of snow and glacier ablation, a spatially distributed [...] Read more.
Efficient water resource management in glacier- and snow-dominated basins requires accurate estimates of the snow water equivalent (SWE) in late winter and spring and melt onset timing and intensity. To understand the high spatio-temporal variability of snow and glacier ablation, a spatially distributed energy balance model combining satellite-based retrievals of albedo and snow cover was applied. Incoming short-wave energy, contributing to daily estimates of melt energy, was constrained by remotely sensed surface albedo for snow-covered surfaces. Fractional snow cover was used for non-glaciated areas, as it provides estimates of snow cover for each pixel to better constrain snow melt. Thus, available daily estimates of melt energy in a given area were the product of the possible melt energy and the fractional snow cover of the area or pixel for non-glaciated areas. This provided daily estimates of melt water to determine seasonal snow and glacier ablation in Iceland for the period 2000–2019. Observations from snow pits on land and glacier summer mass balance were used for evaluation, and observations from land and glacier-based automatic weather stations were used to evaluate model inputs for the energy balance model. The results show that the interannual SWE variability was generally high both for seasonal snow and glaciers. For seasonal snow, the largest SWE (>1000 mm) was found in mountainous and alpine areas close to the coast, notably in the East- and Westfjords, Tröllaskaga, and in the vicinity of glacier margins. Lower SWE values were observed in the central highlands, flatter inland areas, and at lower elevations. For glaciers, more SWE (glacier ablation) was associated with lower glacier elevations while less melt was observed at higher elevations. For the impurity-rich bare-ice areas that are exposed annually, observed SWE was more than 3000 mm. Full article
(This article belongs to the Topic Hydrology and Water Resources Management)
Show Figures

Figure 1

Figure 1
<p>Overview of study area. Main catchments are outlined with black while glacier boundaries used in this study are red. The letters near each area refer to the catchment identification (ID). Topographic properties for each area are shown in <a href="#hydrology-11-00003-t001" class="html-table">Table 1</a> and <a href="#hydrology-11-00003-t002" class="html-table">Table 2</a>.</p>
Full article ">Figure 2
<p>Evaluation of the downscaled daily RAV2 modeled incoming solar radiation, incoming long-wave radiation, and air temperature with ground observations. (<b>a</b>–<b>c</b>) show results from weather stations operated on glaciers (GAWS) while (<b>d</b>–<b>f</b>) show results from stations outside of glaciers (AWS). Color shows the normalized (0–1) density distribution of data. Dotted black line shows 1:1, and black line shows the calculated linear fit to the data.</p>
Full article ">Figure 3
<p>Comparison of calculated and observed summer mass balance (<math display="inline"><semantics> <msub> <mi>b</mi> <mi>s</mi> </msub> </semantics></math>) for selected glaciers. (<b>a</b>) comparison for Vatnajökull, (<b>b</b>) comparison for Hofsjökull, (<b>c</b>) comparison for Langjökull and (<b>d</b>) comparison for Mýrdalsjökull. The dotted black line shows 1:1 and the gray line represents the calculated linear fit to the data. The color bar refers to the elevation of each observation point in the comparison with respect to the elevation range of each glacier, i.e., the normalized elevation range of each glacier.</p>
Full article ">Figure 4
<p>(<b>a</b>) Comparison of SWE from 1 April 2015 to 2019 for observed SWE using a CS725 and estimation from the reconstructed SWE. (<b>b</b>) Comparison of observed maximum SWE in spring (mid-March to mid-April) against reconstructed SWE for selected years.</p>
Full article ">Figure 5
<p>Spatial pattern mean of reconstructed SWE and glacier ablation for the period 2000–2019. Glacier ablation is shown as negative values (white to red) and seasonal snow is shown as positive values(white to blue) to better distinguish seasonal snow and glacier ablation as their total magnitudes generally had different ranges.</p>
Full article ">Figure 6
<p>Spatial patterns of mean reconstructed SWE for the period 2000–2019 for individual months from April to September. Note that the color scale is different for each figure.</p>
Full article ">Figure 7
<p>Annual spatial patterns for SWE anomalies for 2000–2019. Red colors denote positive values where melt was above the average, i.e., more melt, while blue colors show less melt.</p>
Full article ">Figure 8
<p>Reconstructed SWE for the main catchments and individual years.</p>
Full article ">Figure 9
<p>Variability of ratio between melt contribution from seasonal snow and glacier ablation. Gray circles show the catchment area ratio between land and glacier. On each box, the central mark indicates the median, and the bottom and top edges of the box indicate the 25th and 75th percentiles, respectively. The whiskers extend to the most extreme data points not considered outliers. Outliers are shown as blue circles.</p>
Full article ">Figure 10
<p>Melt in mm (upper panel) and melt volume anomalies in percent (lower figure) for the main Icelandic glaciers and their sub-areas. A melt anomaly of 100% represents the mean for the period of 2000–2019 while values below 100% represent below-average melt and vice versa. Glaciers and their sub-outlets are listed from the lowest average melt (<b>top</b>) to the highest (<b>bottom</b>).</p>
Full article ">Figure 11
<p>Spatial pattern of mean melt season (AMJJAS) trends in terms of the total change of a least-square fit to the SWE from 2000–2019. Statistically significant trends (<span class="html-italic">p</span> &lt; 0.05) are shown with stipples.</p>
Full article ">
Back to TopTop