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 (287)

Search Parameters:
Keywords = land surface albedo

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
22 pages, 9741 KiB  
Article
Assessing Green Strategies for Urban Cooling in the Development of Nusantara Capital City, Indonesia
by Radyan Putra Pradana, Vinayak Bhanage, Faiz Rohman Fajary, Wahidullah Hussainzada, Mochamad Riam Badriana, Han Soo Lee, Tetsu Kubota, Hideyo Nimiya and I Dewa Gede Arya Putra
Climate 2025, 13(2), 30; https://doi.org/10.3390/cli13020030 - 31 Jan 2025
Viewed by 886
Abstract
The relocation of Indonesia’s capital to Nusantara in East Kalimantan has raised concerns about microclimatic impacts resulting from proposed land use and land cover (LULC) changes. This study explored strategies to mitigate these impacts by using dynamical downscaling with the Weather Research and [...] Read more.
The relocation of Indonesia’s capital to Nusantara in East Kalimantan has raised concerns about microclimatic impacts resulting from proposed land use and land cover (LULC) changes. This study explored strategies to mitigate these impacts by using dynamical downscaling with the Weather Research and Forecasting model integrated with the urban canopy model (WRF-UCM). Numerical experiments at a 1 km spatial resolution were used to evaluate the impacts of green and mitigation strategies on the proposed master plan. In this process, five scenarios were analyzed, incorporating varying proportions of blue–green spaces and modifications to building walls and roof albedos. Among them, scenario 5, with 65% blue–green spaces, exhibited the highest cooling potential, reducing average urban surface temperatures by approximately 2 °C. In contrast, scenario 4, which allocated equal shares of built-up areas and mixed forests (50% each), achieved a more modest reduction of approximately 1 °C. The adoption of nature-based solutions and sustainable urban planning in Nusantara underscores the feasibility of climate-resilient urban development. This framework could inspire other cities worldwide, showcasing how urban growth can align with environmental sustainability. Full article
(This article belongs to the Special Issue Applications of Smart Technologies in Climate Risk and Adaptation)
Show Figures

Figure 1

Figure 1
<p>Study area map representing (<b>a</b>) Indonesia and the domains for the Weather Research and Forecast (WRF) model and the location of the Nusantara capital city and (<b>b</b>) the location of the government center core area (KIPP), the main area of Nusantara (KIKN), and the entire area, including the Nusantara capital city future development plan (IKN) and the actual land use and land cover (LULC) classes across the Nusantara region.</p>
Full article ">Figure 2
<p>Variations in the LULC data for the KIKN area used for the WRF numerical simulations: (<b>a</b>) scenario 1 (before development); (<b>b</b>) scenario 2, representing the baseline (50% greenery and 50% urban); (<b>c</b>) scenario 3 (50% grasslands and 50% urban); (<b>d</b>) scenario 4 (50% mixed forest and 50% urban); (<b>e</b>) scenario 5 (65% greenery and 35% urban); (<b>f</b>) scenario 6 (35% greenery and 65% urban).</p>
Full article ">Figure 3
<p>Spatial distributions of LULC (left column) and category-wise surface air temperature probability distributions (right column) at 01:00 and 16:00 local time for the Nusantara area domain for different scenarios: scenario 1 (<b>a</b>,<b>b</b>), scenario 2 (<b>c</b>,<b>d</b>), scenario 3 (<b>e</b>,<b>f</b>), scenario 4 (<b>g</b>,<b>h</b>), scenario 5 (<b>i</b>,<b>j</b>), and scenario 6 (<b>k</b>,<b>l</b>). The dashed lines indicate the north–south (NS) and east–west (EW) cross-sections used in <a href="#climate-13-00030-f004" class="html-fig">Figure 4</a>.</p>
Full article ">Figure 4
<p>Differences in surface air temperature (Temp) and wind speed (WS) at 01:00 and 16:00 local time at the east–west (EW) and north–south (NS) cross-sections before (scenario 1) and after (scenarios 2–6) Nusantara city development; Temp and WS differences (<b>a</b>) between scenario 1 and scenario 2, (<b>b</b>) between scenario 1 and scenario 3, (<b>c</b>) between scenario 1 and scenario 4, (<b>d</b>) between scenario 1 and scenario 5, and (<b>e</b>) between scenario 1 and scenario 6. <a href="#climate-13-00030-f003" class="html-fig">Figure 3</a> shows the locations of the EW and NS cross-sections. The mitigation measures were incorporated into scenarios 2–6 by adjusting the albedo values to 0.8 for roofs and 0.7 for walls. (For references to color blocks in this figure legend, please see <a href="#climate-13-00030-f003" class="html-fig">Figure 3</a>).</p>
Full article ">Figure 5
<p>Comparisons of the simulated surface air temperature (Temp) and wind speed (WS) from scenario 2 at 01:00 and 16:00 local time along (<b>A</b>) EW1, (<b>B</b>) EW2, (<b>C</b>) NS1, and (<b>D</b>) NS2 cross-sections and (<b>E</b>) the LULC pattern of scenario 2 with the locations of the cross-sections.</p>
Full article ">Figure 6
<p>Comparison of the average hourly surface air temperature over computational domain 3 for various scenarios (<b>a</b>) with water bodies, (<b>b</b>) without water bodies, and (<b>c</b>) the difference in air temperature between (<b>b</b>) and (<b>a</b>).</p>
Full article ">Figure A1
<p>(<b>a</b>) Observed monthly variations in rainfall and surface air temperature in the study area from Jan 2016–Dec 2020, and (<b>b</b>) the highest average surface air temperature occurred on 21 October 2020 during the five-year period, as indicated by the red circle.</p>
Full article ">Figure A2
<p>Paired comparison of simulated surface air temperature with wind speed before and after Nusantara city development: scenario 2 (<b>a</b>,<b>b</b>), scenario 3 (<b>c</b>,<b>d</b>), scenario 4 (<b>e</b>,<b>f</b>), scenario 5 (<b>g</b>,<b>h</b>), and scenario 6 (<b>i</b>,<b>j</b>). (For references to color blocks in this figure legend, please see <a href="#climate-13-00030-f003" class="html-fig">Figure 3</a>).</p>
Full article ">Figure A2 Cont.
<p>Paired comparison of simulated surface air temperature with wind speed before and after Nusantara city development: scenario 2 (<b>a</b>,<b>b</b>), scenario 3 (<b>c</b>,<b>d</b>), scenario 4 (<b>e</b>,<b>f</b>), scenario 5 (<b>g</b>,<b>h</b>), and scenario 6 (<b>i</b>,<b>j</b>). (For references to color blocks in this figure legend, please see <a href="#climate-13-00030-f003" class="html-fig">Figure 3</a>).</p>
Full article ">Figure A3
<p>Spatial wind patterns at 01:00 and 16:00 local time on 21 October 2020. (For references to color blocks in this figure legend, please see <a href="#climate-13-00030-f003" class="html-fig">Figure 3</a>).</p>
Full article ">
49 pages, 68388 KiB  
Article
Improved Stereophotogrammetric and Multi-View Shape-from-Shading DTMs of Occator Crater and Its Interior Cryovolcanism-Related Bright Spots
by Alicia Neesemann, Stephan van Gasselt, Ralf Jaumann, Julie C. Castillo-Rogez, Carol A. Raymond, Sebastian H. G. Walter and Frank Postberg
Remote Sens. 2025, 17(3), 437; https://doi.org/10.3390/rs17030437 - 27 Jan 2025
Viewed by 464
Abstract
Over the course of NASA’s Dawn Discovery mission, the onboard framing camera mapped Ceres across a wide wavelength spectrum at varying polar science orbits and altitudes. With increasing resolution, the uniqueness of the 92 km wide, young Occator crater became evident. Its central [...] Read more.
Over the course of NASA’s Dawn Discovery mission, the onboard framing camera mapped Ceres across a wide wavelength spectrum at varying polar science orbits and altitudes. With increasing resolution, the uniqueness of the 92 km wide, young Occator crater became evident. Its central cryovolcanic dome, Cerealia Tholus, and especially the associated bright carbonate and ammonium chloride deposits—named Cerealia Facula and the thinner, more dispersed Vinalia Faculae—are the surface expressions of a deep brine reservoir beneath Occator. Understandably, this made this crater the target for future sample return mission studies. The planning and preparation for this kind of mission require the characterization of potential landing sites based on the most accurate topography and orthorectified image data. In this work, we demonstrate the capabilities of the freely available and open-source USGS Integrated Software for Imagers and Spectrometers (ISIS 3) and Ames Stereo Pipeline (ASP 2.7) in creating high-quality image data products as well as stereophotogrammetric (SPG) and multi-view shape-from-shading (SfS) digital terrain models (DTMs) of the aforementioned spectroscopically challenging features. The main data products of our work are four new DTMs, including one SPG and one SfS DTM based on High-Altitude Mapping Orbit (HAMO) (CSH/CXJ) and one SPG and one SfS DTM based on Low-Altitude Mapping Orbit (LAMO) (CSL/CXL), along with selected Extended Mission Orbit 7 (XMO7) framing camera (FC) data. The SPG and SfS DTMs were calculated to a GSD of 1 and 0.5 px, corresponding to 136 m (HAMO SPG), 68 m (HAMO SfS), 34 m (LAMO SPG), and 17 m (LAMO SfS). Finally, we show that the SPG and SfS approaches we used yield consistent results even in the presence of high albedo differences and highlight how our new DTMs differ from those previously created and published by the German Aerospace Center (DLR) and the Jet Propulsion Laboratory (JPL). Full article
Show Figures

Figure 1

Figure 1
<p>CSH CYCLE 1 photometrically corrected RGB orthomosaic of the Occator crater (F5IR, F2GREEN, F8BLUE). The map is a stereographic projection with the projection and image center at 22.879°N<sub><math display="inline"><semantics> <mi>φ</mi> </semantics></math></sub>/239.429°E (19.865°N<sub><math display="inline"><semantics> <mi>ψ</mi> </semantics></math></sub>). RGB values were limited to R: 0.0257–0.0339; G: 0.0555–0.0869; and B: 0.0552–0.0815. (<b>a</b>) Due to the high albedo difference, bright deposits within Occator appear overexposed in the applied histogram stretch. To get an idea about their shape and distribution, we compiled a CSL/CXL RGB orthomosaic with a histogram stretch optimized for the bright deposits, shown in <a href="#remotesensing-17-00437-f002" class="html-fig">Figure 2</a>. (<b>b</b>) Extent of our CSL/CXL ASP SPG DTM. (<b>c</b>,<b>d</b>) Areas for which we calculated the highest resolution DTM (CSL/CXL ASP SfS DTMs).</p>
Full article ">Figure 2
<p>Photometrically corrected CSL/CXL RGB orthomosaic of Occator’s interior. The mosaic was compiled from 13 individual images per channel (see <a href="#remotesensing-17-00437-t0A2" class="html-table">Table A2</a>) using F5IR, F2GREEN, and F8BLUE as the three RGB bands. In contrast to <a href="#remotesensing-17-00437-f001" class="html-fig">Figure 1</a>, images were photometrically corrected based on the 482 × 446 km ellipsoid, not on the DTM, to preserve the topography-related brightness variations and morphology, respectively. RGB values were limited to R: 0.0016–0.6634; G: 0.0003–1.5579; and B: 0.0000–1.1953. The figure is a stereographic projection centered at 23.05°N<sub><math display="inline"><semantics> <mi>φ</mi> </semantics></math></sub>/241.05°E (20.02°N<sub><math display="inline"><semantics> <mi>ψ</mi> </semantics></math></sub>). To the left, we see Cerealia Facula with its fractured dome, Cerealia Tholus, and the orange-colored exposed hydrated sodium chloride [<a href="#B24-remotesensing-17-00437" class="html-bibr">24</a>]. To the right, we see various smaller and thinner bright deposits collectively named Vinalia Faculae (<b>a</b>–<b>e</b>).</p>
Full article ">Figure 3
<p>Number of acquired FC2 images and geometries of polar mapping orbits during Dawn’s mission at Ceres. Left: Number of images acquired by the FC2 during Dawn’s mission to Ceres. Left: During the course of Dawn’s mission at Ceres, 67,969 images in total were acquired by the FC2. Except for images acquired for the purpose of the camera calibration and orientation of the spacecraft, as well as for the search of moons orbiting Ceres, by far the most images mapping Ceres directly were taken during the HAMO and LAMO. Ancillary image acquisition was carried out by the FC1 in order to increase spatial coverage during Dawn’s time-limited final mission phase, XM2, but is not included in this figure (see <a href="#remotesensing-17-00437-t001" class="html-table">Table 1</a>). Right: Indicated orbits correspond to the median distance to Ceres’ center during the different orbit phases. 360° corresponds to the period between 2015/01/01 and 2018/12/31. Geometries of Dawn’s highly elliptical 2nd extended mission orbits (XMO5–XMO7) flown during the final mission, extended mission XM2, are not shown in this figure for scale reasons.</p>
Full article ">Figure 4
<p>Detection of erroneous pixels. In the upper two rows, we present average images created for the eight individual FC filters and two average flat-field images for the F1 and F7 filters taken while the front door was closed and the calibration lamp (callamp) on. Note that images used in this context were not photometrically but only radiometrically corrected. The east–west shading therefore stems from the illumination conditions during image acquisition and not from camera shading. The five static erroneous pixel clusters recognized in all the average images are marked by an ‘x’ in the upper left subfigure. In the bottom row, we present magnifications of the five erroneous pixel clusters to illustrate their extent, marked by the dashed outlines. A 3 sigma histogram stretch was applied to each average image.</p>
Full article ">Figure 5
<p>Basic ISIS 3 pre-processing workflow. Gray: functions; red: temporary textfiles (here, GCP.net files); green: temporary raster files; blue: final photometrically corrected, orthorectified F1CLEAR image; yellow: F1CLEAR orthomosaic, HAMO-based SPG DTM, and reconstructed SPICE kernels. The asterisk in campt* stands for a script we wrote that reads out the lat/lon values at a specific sample/line position, converts them into cartesian coordinates (the ApprioryX, Y, and Z values) and creates a <span class="html-italic">qtie</span>-readable GCP netfile. A detailed description of the flowchart is given in <a href="#sec4dot4-remotesensing-17-00437" class="html-sec">Section 4.4</a>.</p>
Full article ">Figure 6
<p>ASP SPG processing workflow. This figure is a version of Figure 14.1 from the ASP 2.7 documentation (<a href="https://stereopipeline.readthedocs.io/en/latest/correlation.html" target="_blank">https://stereopipeline.readthedocs.io/en/latest/correlation.html</a>, accessed on 23 Decemeber 2024) that we have modified. The parameters specified in the stereo.txt file and passed to the <span class="html-italic">stereo</span> command are listed in Appendix <a href="#remotesensing-17-00437-t0A8" class="html-table">Table A8</a>. Other parameters passed to the <span class="html-italic">point2dem</span> command to remove additional erroneous points from the point cloud during the DTM generation are described in <a href="#sec4dot5-remotesensing-17-00437" class="html-sec">Section 4.5</a>.</p>
Full article ">Figure 7
<p>Latitude- and longitude-dependent deviation between the CSH/CXJ and CSL/CXL ASP SPG and SfS DTMs created in our study. For the upper two plots, the mean latitude and longitude values were calculated from the extent of area 1, while they were calculated for area 4 in the four lower plots.</p>
Full article ">Figure 8
<p>Latitude- and longitude-dependent deviation between the HAMO and LAMO SPG DTMs. Mean latitude and longitude values were calculated for the extent of area 3.</p>
Full article ">Figure 9
<p>Latitude- and longitude-dependent deviation of our 4 new CSH/CXJ and CSL/CXL SPG and SfS DTMs and the JPL HAMO/LAMO SPC DTM [<a href="#B58-remotesensing-17-00437" class="html-bibr">58</a>,<a href="#B59-remotesensing-17-00437" class="html-bibr">59</a>] from the DLR HAMO SPG DTM [<a href="#B55-remotesensing-17-00437" class="html-bibr">55</a>,<a href="#B56-remotesensing-17-00437" class="html-bibr">56</a>].</p>
Full article ">Figure 10
<p>Visual comparison of the seven DTMs for the example of the fresh crater located at 14.281°N<sub><math display="inline"><semantics> <mi>φ</mi> </semantics></math></sub>/233.489°E (12.295°N<sub><math display="inline"><semantics> <mi>ψ</mi> </semantics></math></sub>). Images used to create the orthomosaics are listed in <a href="#remotesensing-17-00437-t0A5" class="html-table">Table A5</a>.</p>
Full article ">Figure 11
<p>Topographic profiles of the fresh crater. Subfigures (<b>a</b>,<b>b</b>) (left and middle panel) are stereographic projections centered at 14.281°N<sub><math display="inline"><semantics> <mi>φ</mi> </semantics></math></sub>/233.489°E (12.295°N<sub><math display="inline"><semantics> <mi>ψ</mi> </semantics></math></sub>). (<b>a</b>) Photometrically corrected FC2 F1CLEAR close-up view of the fresh crater. (<b>b</b>) Photometrically corrected RGB (F5IR, F2GREEN, F8BLUE) orthomosaic of the fresh crater. In total, we derived 52 profiles at 3 degree intervals between 24–90°, 144–177°, and 240–288° for each of the seven DTMs extending from the crater center. (<b>c</b>) (right panel) Average topographic profiles for each of the seven DTMs. <sup>a</sup> [<a href="#B59-remotesensing-17-00437" class="html-bibr">59</a>], <sup>b</sup> [<a href="#B58-remotesensing-17-00437" class="html-bibr">58</a>], <sup>c</sup> [<a href="#B106-remotesensing-17-00437" class="html-bibr">106</a>], <sup>d</sup> [<a href="#B104-remotesensing-17-00437" class="html-bibr">104</a>], <sup>e</sup> [<a href="#B56-remotesensing-17-00437" class="html-bibr">56</a>], <sup>f</sup> [<a href="#B55-remotesensing-17-00437" class="html-bibr">55</a>].</p>
Full article ">Figure 12
<p>Detailed comparison of the topographic profiles of the fresh crater. As a reference (black solid line), we used our new CSL/CXL ASP SfS DTM, as it has the highest effective resolution and the highest d/D ratio of 0.255 and plotted it together with the profiles extracted from the other six DTMs (<b>a</b>–<b>f</b>). Black and gray triangles mark the inflection points (the highest elevation of the rim crest) of our reference profile and the other profiles. Note that the highest congruence exists between profiles taken from our CSL/CXL ASP SfS, our CSL/CXL ASP SPG, and the DLR CSL/CXL SPG DTMs.</p>
Full article ">Figure 13
<p>Visual comparison of the seven DTMs for the example of the Cerealia Tholus at 22.626°N<sub><math display="inline"><semantics> <mi>φ</mi> </semantics></math></sub>/239.581°E (19.648°N<sub><math display="inline"><semantics> <mi>ψ</mi> </semantics></math></sub>). Images used to create the orthomosaics are listed in <a href="#remotesensing-17-00437-t0A6" class="html-table">Table A6</a>.</p>
Full article ">Figure 14
<p>Overview of the location of the topographic profile across the Cerealia Tholus. We created three different maps, each with their advantages and disadvantages, in order to show which specific surface features are covered by our topographic profile. Our 25,498 m long profile goes from west to east while crossing the highest elevations (the Lohri, Cerealia, and Kekri Tholus) within Occator. (<b>a</b>) Semi-transparent, color-coded CSL/CXL/XMO7 ASP SfS DTM overlaid on the corresponding hillshade model. Topography contour lines are plotted in 100 m intervals. (<b>b</b>) CSL/CXL RGB color composite of Cerealia Facula. (<b>c</b>) Generated slope map overlaid on a curvature map. This combined map highlights aspects of the surface’s shape or features, such as the circular and other rather subparallel fault systems as well as numerous little mounds, at a detailed level.</p>
Full article ">Figure 15
<p>Topographic profiles across the Cerealia Tholus.</p>
Full article ">Figure 16
<p>Topographic profiles across Vinalia Faculae. (<b>a</b>) Topographic profiles across Vinalia Faculae, extracted from our CSL/CXL SPG (black line) and SfS (dark grey line) DTMs, as well as the HAMO/LAMO-based SPC DTM (blue line) from the JPL. Additionally, the albedo along the profile line was extracted based on the photometrically corrected CSL/CXL F1CLEAR orthomosaic included in this study. (<b>b</b>) Deviations between the lower resolution yet more robust CSL/CXL ASP SPG DTM, the CSL/CXL ASP SfS DTM, and the HAMO/LAMO-based SPC DTM from the JPL. (<b>c</b>) The CSL/CXL ASP SfS DTM of Vinalia Faculae, represented as elevations above the 482 × 446 km ellipsoid. (<b>d</b>) Photometrically corrected CSL/CXL F1CLEAR orthomosaic of Vinalia Faculae. Both maps are stereographic projections, centered at 23.292°N<sub><math display="inline"><semantics> <mi>φ</mi> </semantics></math></sub>/242.487°E (20.234°N<sub><math display="inline"><semantics> <mi>ψ</mi> </semantics></math></sub>). The course of the topographic profiles shown in panel (<b>a</b>) is indicated by a black-and-green dashed line.</p>
Full article ">
18 pages, 11390 KiB  
Article
Quantifying the Carbon Reduction Potential of Urban Parks Under Extreme Heat Events Using Interpretable Machine Learning: A Case Study of Jinan, China
by Lemin Yu, Wenru Li, Changhui Zheng and Xiaowen Lin
Atmosphere 2025, 16(1), 79; https://doi.org/10.3390/atmos16010079 - 14 Jan 2025
Viewed by 543
Abstract
Greenhouse gas emissions are primary drivers of climate change, and the intensification of extreme heat and urban heat island effects poses serious threats to urban ecosystems, public health, and energy consumption. This study systematically evaluated the carbon reduction potential of 369 urban parks [...] Read more.
Greenhouse gas emissions are primary drivers of climate change, and the intensification of extreme heat and urban heat island effects poses serious threats to urban ecosystems, public health, and energy consumption. This study systematically evaluated the carbon reduction potential of 369 urban parks in Jinan during extreme heat events using land surface temperature (LST) retrieval, combined with CatBoost + SHAP machine learning methods. Results indicate that the LST in Jinan ranged from 1.77 °C to 59.44 °C, and 278 parks exhibited significant cooling effects, collectively saving 2943 tons of CO2 per day—offsetting 11.28% of the city’s fossil fuel emissions. Small parks, such as community parks, demonstrated higher carbon-saving efficiency (CSE), while large ecological parks showed greater carbon-saving intensity (CSI). CSE was strongly correlated with vegetation coverage and surrounding population density, with efficiency increasing when the vegetation index was within 0.3–0.7 and population density ranged 0–5000 or 15,000–22,500 people. CSI was influenced by evapotranspiration and park geometric form, increasing significantly when the park area exceeded 250 hectares or evapotranspiration ranged 2.5–6.0. However, elevation and albedo negatively impacted both metrics, with the lowest CSI observed when elevation exceeded 150 m or albedo surpassed 18%. Full article
(This article belongs to the Special Issue Urban Impact on the Low Atmosphere Processes)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>–<b>c</b>) Study area and distribution of urban park types.</p>
Full article ">Figure 2
<p>Research flowchart.</p>
Full article ">Figure 3
<p>(<b>a</b>) Cooling curves for the park. (<b>b</b>–<b>d</b>) Carbon-saving modeling for the outdoors.</p>
Full article ">Figure 4
<p>Cross-verification of RTE-based LST.</p>
Full article ">Figure 5
<p>(<b>a</b>) Land surface temperature spatial distribution; (<b>b</b>) CSI distribution; (<b>c</b>) CSE distribution.</p>
Full article ">Figure 6
<p>(<b>a</b>) Contribution to carbon neutrality; (<b>b</b>) carbon offset by park type.</p>
Full article ">Figure 7
<p>Correlation coefficients between carbon reduction indicators of urban parks and internal and external factors (*** indicates a significance level of 0.001, ** indicates a significance level of 0.01, * indicates a significance level of 0.05).</p>
Full article ">Figure 8
<p>Summary map of factors affecting carbon reduction in park green spaces.</p>
Full article ">Figure 9
<p>Single-factor feature dependency plot. (<b>a</b>) Park_Area; (<b>b</b>) Park_DEM; (<b>c</b>) Park_ET; (<b>d</b>) Park_Albedo; (<b>e</b>) Buffer_Tree; (<b>f</b>) Park_Area(CSE); (<b>g</b>) Park_DEM(CSE); (<b>h</b>) Park_NDVI(CSE); (<b>i</b>) Buffer_Tree(CSE); (<b>j</b>) Buffer_POP(CSE).</p>
Full article ">Figure 10
<p>Partial factor bivariate dependency plot. (<b>a</b>) Carbon reduction intensity. (<b>b</b>) Carbon reduction efficiency.</p>
Full article ">
30 pages, 60239 KiB  
Article
Retrieval and Evaluation of Global Surface Albedo Based on AVHRR GAC Data of the Last 40 Years
by Shaopeng Li, Xiongxin Xiao, Christoph Neuhaus and Stefan Wunderle
Remote Sens. 2025, 17(1), 117; https://doi.org/10.3390/rs17010117 - 1 Jan 2025
Viewed by 829
Abstract
In this study, the global land surface albedo namely GAC43 was retrieved for the years 1979 to 2020 using Advanced Very High Resolution Radiometer (AVHRR) global area coverage (GAC) data onboard National Oceanic and Atmospheric Administration (NOAA) and Meteorological Operational (MetOp) satellites. We [...] Read more.
In this study, the global land surface albedo namely GAC43 was retrieved for the years 1979 to 2020 using Advanced Very High Resolution Radiometer (AVHRR) global area coverage (GAC) data onboard National Oceanic and Atmospheric Administration (NOAA) and Meteorological Operational (MetOp) satellites. We provide a comprehensive retrieval process of the GAC43 albedo, followed by a comprehensive assessment against in situ measurements and three widely used satellite-based albedo products, the third edition of the CM SAF cLoud, Albedo and surface RAdiation (CLARA-A3), the Copernicus Climate Change Service (C3S) albedo product, and MODIS BRDF/albedo product (MCD43). Our quantitative evaluations indicate that GAC43 demonstrates the best stability, with a linear trend of ±0.002 per decade at nearly all pseudo invariant calibration sites (PICS) from 1982 to 2020. In contrast, CLARA-A3 exhibits significant noise before the 2000s due to the limited availability of observations, while C3S shows substantial biases during the same period due to imperfect sensors intercalibrations. Extensive validation at globally distributed homogeneous sites shows that GAC43 has comparable accuracy to C3S, with an overall RMSE of approximately 0.03, but a smaller positive bias of 0.012. Comparatively, MCD43C3 shows the lowest RMSE (~0.023) and minimal bias, while CLARA-A3 displays the highest RMSE (~0.042) and bias (0.02). Furthermore, GAC43, CLARA-A3, and C3S exhibit overestimation in forests, with positive biases exceeding 0.023 and RMSEs of at least 0.028. In contrast, MCD43C3 shows negligible bias and a smaller RMSE of 0.015. For grasslands and shrublands, GAC43 and MCD43C3 demonstrate comparable estimation uncertainties of approximately 0.023, with close positive biases near 0.09, whereas C3S and CLARA-A3 exhibit higher RMSEs and biases exceeding 0.032 and 0.022, respectively. All four albedo products show significant RMSEs around 0.035 over croplands but achieve the highest estimation accuracy better than 0.020 over deserts. It is worth noting that significant biases are typically attributed to insufficient spatial representativeness of the measurement sites. Globally, GAC43 and C3S exhibit similar spatial distribution patterns across most land surface conditions, including an overestimation compared to MCD43C3 and an underestimation compared to CLARA-A3 in forested areas. In addition, GAC43, C3S, and CLARA-A3 estimate higher albedo values than MCD43C3 in low-vegetation regions, such as croplands, grasslands, savannas, and woody savannas. Besides the fact that the new GAC43 product shows the best stability covering the last 40 years, one has to consider the higher proportion of backup inversions before 2000. Overall, GAC43 offers a promising long-term and consistent albedo with good accuracy for future studies such as global climate change, energy balance, and land management policy. Full article
Show Figures

Figure 1

Figure 1
<p>Local solar times and solar zenith angles of equator observations for all AVHRR-carrying NOAA and MetOp satellites used to generate GAC43 albedo products as shown in (<b>a</b>,<b>b</b>), respectively. SZA &gt; 90° indicates night conditions.</p>
Full article ">Figure 2
<p>Globally distributed sites with homogeneous characteristics and corresponding land cover types defined by the IGBP from the MCD12C1 product. Purple squares located in the desert are used to evaluate temporal stability, while other sites are utilized for direct validations.</p>
Full article ">Figure 3
<p>Flowchart for this study.</p>
Full article ">Figure 4
<p>The performance of full inversion and full and backup inversion at various IGBP land cover types.</p>
Full article ">Figure 5
<p>The performance of the GAC43 albedo with full inversions at various land cover types, where panels (<b>a</b>–<b>h</b>) represent the land cover types of BSV, CRO, DBF, EBF, ENF, GRA, OSH and WSA, respectively. In the plots, the red solid line represents the 1:1 line, and the green dotted line and purple solid lines represent the limits of deviation ±0.02 and ±0.04, respectively.</p>
Full article ">Figure 6
<p>Google Earth <sup>TM</sup> images were used to visually illustrate the heterogeneity surrounding selected homogeneous sites representing various land cover types: (<b>a</b>) EBF, (<b>b</b>) BSV, (<b>c</b>) CRO and (<b>d</b>) GRA, as defined by the MCD12C1 IGBP classification. The red circle in each image denotes a radius of 2.5 km.</p>
Full article ">Figure 7
<p>Inter-comparison performance among four satellite-based albedo products. The top four subfigures (<b>a</b>–<b>d</b>) show the accuracy of all available matching samples between in situ measurements and estimated albedo values derived from satellite products, while the bottom four subfigures (<b>e</b>–<b>h</b>) give the performance of that using same samples.</p>
Full article ">Figure 8
<p>The performance of four satellite-based albedo products using same samples across various land surface types, evaluated in terms of (<b>a</b>) RMSE and (<b>b</b>) bias, respectively. The <span class="html-italic">x</span>-axis represents the land cover type classified as forest, grassland or shrublands, cropland, and desert, and corresponding available samples.</p>
Full article ">Figure 9
<p>The temporal performance of four satellite-based albedo products related to in situ measurements, and each subplot represents one case of different land cover surface, including (<b>a</b>) EBF, (<b>b</b>) ENF, (<b>c</b>) DBF, (<b>d</b>) GRA, and (<b>e</b>) CRO, respectively. The grey shaded areas depict situations with snow cover.</p>
Full article ">Figure 10
<p>Spatial distributions of GAC43 BSA in July 2013 are shown in subgraph (<b>a</b>), with corresponding differences from (<b>b</b>) CLARA-A3, (<b>c</b>) C3S, and (<b>d</b>) MCD43C3 in the same month, respectively.</p>
Full article ">Figure 11
<p>Percentage difference in BSA values between (<b>a</b>) GAC43 and CLARA-A3, (<b>b</b>) GAC43 and C3S, and (<b>c</b>) GAC43 and MCD43C3 in July 2013.</p>
Full article ">Figure 12
<p>The scattering plots between GAC43 BSA and (<b>a</b>) CLARA-A3 BSA, (<b>b</b>) C3S BSA, and (<b>c</b>) MCD43C3 BSA using all snow-free monthly pixels in July 2013, where the red lines indicate 1:1.</p>
Full article ">Figure 13
<p>The monthly BSA for the four satellite-based products across various land cover types in July 2013, where panels (<b>a</b>–<b>i</b>) represent the land cover types of CRO, DBF, DNF, EBF, ENF, GRA, MF, SAV and WSA, respectively. In the plots, the bottom values of each albedo product are the median of all corresponding land cover estimates. The top values match available samples.</p>
Full article ">Figure 14
<p>Monthly BSA from GAC43, MCD43C3, C3S, and CALRA-A3 at three randomly selected PICS sites: (<b>a</b>) Arabia 2, 20.19°N, 51.63°E; (<b>b</b>) Libya 3, 23.22°N, 23.23°E; and (<b>c</b>) Sudan 1, 22.11°N, 28.11°E, all characterized by BSV land surfaces as defined by IGBP.</p>
Full article ">Figure 15
<p>Box plots of the slope per decade for GAC43, CLARA-A3, C3S, and MCD43C3 at all PICS sites, where (<b>a</b>–<b>d</b>) represent the corresponding statistics during 1982–1990, 1991–2000, 2001–2010 and 2011–2020, respectively, and three dashed grey lines represent the 75%, 50%, and 25% quantiles. Red dotted lines indicate the horizontal line where slope is 0.</p>
Full article ">Figure 16
<p>Percentage of full inversions for the years 2004, 2008, 2012, and 2016 based on GAC43 (<b>top</b>) and MCD43A3 (<b>bottom</b>).</p>
Full article ">Figure 17
<p>Percentage of full inversions of GAC43 at various continents from 1979 to 2020.</p>
Full article ">Figure A1
<p>Spatial distributions of GAC43 BSA in July 2004 are shown in subgraph (<b>a</b>), with corresponding differences from (<b>b</b>) CLARA-A3, (<b>c</b>) C3S, and (<b>d</b>) MCD43C3 in the same month, respectively.</p>
Full article ">Figure A2
<p>Spatial distributions of GAC43 BSA in July 2008 are shown in subgraph (<b>a</b>), with corresponding differences from (<b>b</b>) CLARA-A3, (<b>c</b>) C3S, and (<b>d</b>) MCD43C3 in the same month, respectively.</p>
Full article ">Figure A3
<p>Spatial distributions of GAC43 BSA in July 2012 are shown in subgraph (<b>a</b>), with corresponding differences from (<b>b</b>) CLARA-A3, (<b>c</b>) C3S, and (<b>d</b>) MCD43C3 in the same month, respectively.</p>
Full article ">Figure A4
<p>Spatial distributions of GAC43 BSA in July 2016 are shown in subgraph (<b>a</b>), with corresponding differences from (<b>b</b>) CLARA-A3, (<b>c</b>) C3S, and (<b>d</b>) MCD43C3 in the same month, respectively.</p>
Full article ">Figure A5
<p>Percentages of full inversions for the years between 1979 and 2020 based on GAC43 data record.</p>
Full article ">
21 pages, 5602 KiB  
Article
Quantitative Inversion of Martian Hydrous Minerals Based on LSTM-1DCNN Model
by Xinbao Liu, Ming Jin, Xiangnan Liu, Zhiming Yang, Zengqian Hou and Xiaozhong Ding
Remote Sens. 2025, 17(1), 94; https://doi.org/10.3390/rs17010094 - 30 Dec 2024
Viewed by 686
Abstract
Hydrous minerals are significant indicators of the ancient aqueous environment on Mars, and orbital hyperspectral data are one of the most effective tools for obtaining information about the distribution of hydrous minerals on the Martian surface. However, prolonged weathering, erosion, and other external [...] Read more.
Hydrous minerals are significant indicators of the ancient aqueous environment on Mars, and orbital hyperspectral data are one of the most effective tools for obtaining information about the distribution of hydrous minerals on the Martian surface. However, prolonged weathering, erosion, and other external forces result in complex mixing effects, often weakening the spectral absorption features of individual minerals. This study proposes a quantitative inversion method for Martian hydrous minerals by integrating a radiative transfer model with a deep learning network. Based on the physics of the Hapke radiative transfer model, the single-scattering albedo spectra of mineral end members were obtained. Additionally, the Linear Spectral Mixture Model was employed to generate a large number of fully constrained mineral mixture samples, providing theoretical support for experimental data. An LSTM-1DCNN model was trained to establish a data-driven quantitative inversion framework. CRISM data were applied to the Eberswalde Crater region to retrieve the abundances of 21 hydrous minerals, including tremolite, opal, and serpentine. The average abundance of hydrous minerals was calculated to be 0.018, with a total area proportion of approximately 8%. Additionally, by analyzing the distribution areas of hydrous silicates, hydrous sulfates, and hydrous hydroxides, the water activity history of the region was inferred. The results align with findings from related studies and mineral spectral index results. By incorporating deep learning into traditional mixing models, this study identifies the distribution of various low-abundance hydrous minerals, enhancing the accuracy of Martian hydrous mineral inversion. It is expected to provide valuable references for the selection of landing sites for Tianwen-3 and support the smooth implementation of China’s Mars exploration mission. Full article
(This article belongs to the Special Issue Planetary Geologic Mapping and Remote Sensing (Second Edition))
Show Figures

Figure 1

Figure 1
<p>Overview map of the Eberswalde crater on Mars and the study area FRT000060dd. The Eberswalde crater underwent significant post-impact modifications, leaving only the northeastern crater rim intact. The FRT000060dd image is located northwest of the Eberswalde crater. The western part of the crater features a prominent valley where a river connects to the alluvial fan within the image coverage, forming a relatively complete aqueous landform. The base map data are CTX images.</p>
Full article ">Figure 2
<p>Technology roadmap.</p>
Full article ">Figure 3
<p>LSTM model structure diagram and Bi-LSTM layer structure diagram, where the LSTM layer neuron structure is referenced from Li et al. [<a href="#B32-remotesensing-17-00094" class="html-bibr">32</a>].</p>
Full article ">Figure 4
<p>1DCNN structure and data chart.</p>
Full article ">Figure 5
<p>Model accuracy, precision, recall, and F1 score trends as <span class="html-italic">T</span><sub>1</sub> changes.</p>
Full article ">Figure 6
<p>Qualitative network mineral identification accuracy chart.</p>
Full article ">Figure 7
<p>Scatter plot of the test set for multi-batch training of the quantitative network. Panels (<b>a</b>–<b>f</b>) represent the results for 1DCNN training batches of 100, 150, 200, 250, 300, and 350, respectively.</p>
Full article ">Figure 8
<p>Box plot of predicted mineral abundance within the FRT000060dd map area.</p>
Full article ">Figure 9
<p>Statistical chart of the pixels occupied by major minerals.</p>
Full article ">Figure 10
<p>Comparison of model inversion results and CRISM data. (<b>a</b>) A false-color RGB image (red = 2.5295 μm, green = 1.5066 μm, blue = 1.0800 μm). The red arrow indicates the pixel location. In (<b>b</b>), the solid black line represents the orbital spectral curve of the pixel, whereas the dashed black line shows the simulated spectral curve based on the inversion results. The other dashed, colored lines correspond to the end-member spectra. The variations observed in the 1.4–1.55 μm range, marked by the arrow, result from instrumental errors [<a href="#B33-remotesensing-17-00094" class="html-bibr">33</a>], whereas the sharp peaks in the 1.9–2.1 μm range, also marked by an arrow, are due to atmospheric correction related to CO<sub>2</sub> [<a href="#B10-remotesensing-17-00094" class="html-bibr">10</a>].</p>
Full article ">Figure 11
<p>The distribution map of hydrous mineral types: (<b>a</b>) hydrous sulfate minerals, which have been highlighted with a red box, (<b>b</b>) hydrous hydroxide minerals, and (<b>c</b>) hydrous silicate minerals. The bottom figure is the gray-scale data for the 1.1652 μm.</p>
Full article ">Figure 12
<p>A comparison between mineral abundance inversion results and mineral spectral indices. Panels (<b>a</b>,<b>b</b>) represent the Olivine abundance and Olivine Index. The red rectangular region in Panel (<b>a</b>) represents a potential olivine distribution area. Panel (<b>c</b>) illustrates the scatter plot corresponding to the data presented in Panels (<b>a</b>,<b>b</b>); the Pearson Correlation Coefficient (ρ) between the two variables in Panel (<b>c</b>) is 0.666, indicating a moderately strong positive correlation. Panels (<b>d</b>,<b>e</b>) represent the Pyroxene abundance and Pyroxene Index. Panel (<b>f</b>) illustrates the scatter plot corresponding to the data presented in Panels (<b>d</b>,<b>e</b>); the ρ between the two variables in Panel (<b>f</b>) is 0.335, indicating a certain positive correlation.</p>
Full article ">Figure 13
<p>Horizontal and vertical elevation distribution maps of the three types of hydrous minerals. (<b>a</b>) The horizontal distribution of three types of hydrous minerals: hydrous sulfate minerals, hydrous hydroxide minerals, and hydrous silicate minerals. (<b>b</b>,<b>c</b>) The AA’ and BB’ cross-sectional profiles corresponding to (<b>a</b>), respectively, showing the elevation distribution of different minerals along and near the section lines.</p>
Full article ">
28 pages, 13027 KiB  
Article
From Fields to Microclimate: Assessing the Influence of Agricultural Landscape Structure on Vegetation Cover and Local Climate in Central Europe
by Jan Kuntzman and Jakub Brom
Remote Sens. 2025, 17(1), 6; https://doi.org/10.3390/rs17010006 - 24 Dec 2024
Viewed by 534
Abstract
Agricultural intensification through simplification and specialization has homogenized diverse landscapes, reducing their heterogeneity and complexity. While the negative impact of large, simplified fields on biodiversity has been well-documented, the role of landscape structure in mitigating climatic extremes and stabilizing climate is becoming increasingly [...] Read more.
Agricultural intensification through simplification and specialization has homogenized diverse landscapes, reducing their heterogeneity and complexity. While the negative impact of large, simplified fields on biodiversity has been well-documented, the role of landscape structure in mitigating climatic extremes and stabilizing climate is becoming increasingly important. Despite considerable knowledge of landscape cover types, understanding of how landscape structure influences climatic characteristics remains limited. To explore this further, we studied an area along the Czech–Austrian border, where socio-political factors have created stark contrasts in landscape structure, despite a similar topography. Using Land Parcel Information System (LPIS) data, we analyzed the landscape structure on both sides and processed eight Landsat 8 and 9 OLI/TIRS scenes from the 2022 vegetation season to calculate spectral indices (NDVI, NDMI) and microclimatic features (surface temperature, albedo, and energy fluxes). Our findings revealed significant differences between the two regions. Czech fields, with their larger, simpler structure and lower edge density, can amplify local climatic extremes. In contrast, the distribution of values on the Austrian side was more even, likely due to the greater diversity of cultivated crops, a more spatially diverse landscape, and a balanced spread of agricultural activities over time. In light of climate change and biodiversity conservation, these results emphasize the need to protect and restore landscape complexity to enhance resilience and environmental stability. Full article
Show Figures

Figure 1

Figure 1
<p>Overview map of the area of interest (AOI). The detail shows the landscape structure on both the Czech and Austrian sides of the AOI.</p>
Full article ">Figure 2
<p>Flowchart diagram of the methods used for data analysis and their relationship. Details of each method are given in the text.</p>
Full article ">Figure 3
<p>Summary of the analysis of differences between locations (Czech Republic and Austria) within the observation dates of the area of interest. Contrasts were calculated as <span class="html-italic">post hoc</span> tests of linear mixed models for individual land cover categories and specific biophysical features using the estimated marginal means method. Note: AL—Farmland, AR—Arable Land, VY—Vineyards, PG—Permanent Grassland, FL—Fallow Land, GF—Grass on Field.</p>
Full article ">Figure 4
<p>Violin plots of individual functional variables for the total Farmland in the Czech Republic (orange color) and Austria (blue color).</p>
Full article ">Figure 5
<p>The difference in individual functional variable data distribution between the Czech Republic (orange color) and Austria (blue color) for the total Farmland area.</p>
Full article ">Figure A1
<p>Violin plots of individual functional variables for the Arable Land in the Czech Republic (orange color) and Austria (blue color).</p>
Full article ">Figure A2
<p>Violin plots of individual functional variables for Vineyards in the Czech Republic (orange color) and Austria (blue color).</p>
Full article ">Figure A3
<p>Violin plots of individual functional variables for Permanent Grassland in the Czech Republic (orange color) and Austria (blue color).</p>
Full article ">Figure A4
<p>Violin plots of individual functional variables for Fallow Land in the Czech Republic (orange color) and Austria (blue color).</p>
Full article ">Figure A5
<p>Violin plots of individual functional variables for Grass on Field in the Czech Republic (orange color) and Austria (blue color).</p>
Full article ">Figure A6
<p>The difference in individual functional variable data distribution between the Czech Republic (orange color) and Austria (blue color) for Arable Land.</p>
Full article ">Figure A7
<p>The difference in individual functional variable data distribution between the Czech Republic (orange color) and Austria (blue color) for Vineyards.</p>
Full article ">Figure A8
<p>The difference in individual functional variable data distribution between the Czech Republic (orange color) and Austria (blue color) for Permanent Grasslands.</p>
Full article ">Figure A9
<p>The difference in biophysical features data distribution between the Czech Republic (orange color) and Austria (blue color) for Fallow Land.</p>
Full article ">Figure A10
<p>The difference in individual functional variables data distribution between the Czech Republic (orange color) and Austria (blue color) for Grass on Field.</p>
Full article ">
26 pages, 6384 KiB  
Review
Research Overview on Urban Heat Islands Driven by Computational Intelligence
by Chao Liu, Siyu Lu, Jiawei Tian, Lirong Yin, Lei Wang and Wenfeng Zheng
Land 2024, 13(12), 2176; https://doi.org/10.3390/land13122176 - 13 Dec 2024
Viewed by 1192
Abstract
In recent years, the intensification of the urban heat island (UHI) effect has become a significant concern as urbanization accelerates. This survey comprehensively explores the current status of surface UHI research, emphasizing the role of land use and land cover changes (LULC) in [...] Read more.
In recent years, the intensification of the urban heat island (UHI) effect has become a significant concern as urbanization accelerates. This survey comprehensively explores the current status of surface UHI research, emphasizing the role of land use and land cover changes (LULC) in urban environments. We conducted a systematic review of 8260 journal articles from the Web of Science database, employing bibliometric analysis and keyword co-occurrence analysis using CiteSpace to identify research hotspots and trends. Our investigation reveals that vegetation cover and land use types are the two most critical factors influencing UHI intensity. We analyze various computational intelligence techniques, including machine learning algorithms, cellular automata, and artificial neural networks, used for simulating urban expansion and predicting UHI effects. The study also examines numerical modeling methods, including the Weather Research and Forecasting (WRF) model, while examining the application of Computational Fluid Dynamics (CFD) in urban microclimate research. Furthermore, we evaluate potential mitigation strategies, considering urban planning approaches, green infrastructure solutions, and the use of high-albedo materials. This comprehensive survey not only highlights the critical relationship between land use dynamics and UHIs but also provides a direction for future research in computational intelligence-driven urban climate studies. Full article
(This article belongs to the Special Issue Geospatial Data in Land Suitability Assessment: 2nd Edition)
Show Figures

Figure 1

Figure 1
<p>The flowchart of the research.</p>
Full article ">Figure 2
<p>Number of academic publications on the topic of “Heat Island Effect” from 2014 to 2024.</p>
Full article ">Figure 3
<p>Keyword co-occurrence network in the heat island effects study.</p>
Full article ">Figure 4
<p>UHI effect.</p>
Full article ">Figure 5
<p>LST performance of various land covers in Shenzhen during the UHI phenomenon.</p>
Full article ">Figure 6
<p>Impact of urban heat island effect.</p>
Full article ">Figure 7
<p>Frequency of usage for input parameters.</p>
Full article ">Figure 8
<p>Distribution of research across countries.</p>
Full article ">
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 751
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 ">
19 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 816
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 832
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 1118
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 1371
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 1734
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 1122
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
Cited by 1 | Viewed by 1857
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 ">
Back to TopTop