Landsat 8 OLI Broadband Albedo Validation in Antarctica and Greenland

The albedo is a fundamental component of the processes that govern the energy budget, and particularly important in the context of climate change. However, a satellite-based high-resolution (30 m) albedo product which can be used in the polar regions up to 82.5° latitude during the summer seasons is lacking. To cover this gap, in this study we calculate satellite-based broadband albedo from Landsat 8 OLI and validate it against broadband albedo measurements from in situ stations located on the Antarctic and Greenland icesheets. The model to derive the albedo from raw satellite data includes an atmospheric and topographic correction and conversion from narrow-band to broadband albedo, and at each step different options were taken into account, in order to provide the best combination of corrections. Results, after being cleaned from anomalous data, show a good agreement with in situ albedo measurements, with a mean absolute error between in situ and satellite albedo of 0.021, a root mean square error of 0.026, a standard deviation of 0.015, a correlation coefficient of 0.995 (p < 0.01) and a bias estimate of −0.005. Considering the structure of the model, it could be applied to data from previous sensors of the Landsat family and help construct a record to analyze albedo variations in the polar regions.

ENSO and Light-Absorbing Impurities and Their Impact on Snow Albedo in the Sierra Nevada de Santa Marta, Colombia

Snow albedo is an important variable in the coupled atmosphere-earth system at the global level. Moreover, studying its behavior allows us to know the state of the cryosphere. The Sierra Nevada de Santa Marta (SNSM) is a glacier area and the northernmost tropical (10.82° N, 73.75° W) region in South America. It has a height of up to 5775 m.a.sl., which is the second highest mountain in the world near the marine coast. We analyzed variations in snow albedo related to snow cover, snowfall, temperature, light-absorbing impurities such as blank carbon (BC), organic carbon (OC) and dust, and El Niño—Southern Oscillation (ENSO) phenomenon through 20 years (2000–2020). We mainly use daily data from the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the Terra and Aqua NASA satellites. Results showed through correlations that snow albedo has decreased due to Land Surface Temperature (55%, p < 0.001), a positive phase of ENSO (42%, p < 0.001) and dust (37%, p < 0.01) in the SNSM. Additionally, a dust negative effect was more evident on the southern side (up to 49%, p < 0.001) of the SNSM. Backward trajectories by the NOAA HYSPLIT model suggest that dust sources would be soil erosion in the surrounding region. Results can help recognize the influence of ENSO and dust in the glacier decrease of the SNSM.

Evaluation of the Snow Albedo Retrieved from the Snow Kernel Improved the Ross-Roujean BRDF Model

The original kernel-driven bidirectional reflectance distribution function (BRDF) models were developed based on soil-vegetation systems. To further improve the ability of the models to characterize the snow surface scattering properties, a snow kernel was derived from the asymptotic radiative transfer (ART) model and used in the kernel-driven BRDF model framework. However, there is a need to further evaluate the influence of using this snow kernel to improve the original kernel-driven models in snow albedo retrieval applications. The aim of this study is to perform such an evaluation using a variety of snow BRDF data. The RossThick-Roujean (RTR) model is used as a framework for taking in the new snow kernel (hereafter named the RTS model) since the Roujean geometric-optical (GO) kernel captures a neglectable hotspot effect and represents a more prominent dome-shaped BRDF, especially at a small solar zenith angle (SZA). We obtained the following results: (1) The RTR model has difficulties in reconstructing the snow BRDF shape, especially at large SZAs, which tends to underestimate the reflectance in the forward direction and overestimate reflectance in the backward direction for various data sources. In comparison, the RTS model performs very well in fitting snow BRDF data and shows high accuracy for all data. (2) The RTR model retrieved snow albedos at SZAs = 30°–70° are underestimated by 0.71% and 0.69% in the red and near-infrared (NIR) bands, respectively, compared with the simulation results of the bicontinuous photon tracking (bic-PT) model, which serve as “real” values. However, the albedo retrieved by the RTS model is significantly improved and generally agrees well with the simulation results of the bic-PT model, although the improved model still somewhat underestimates the albedo by 0.01% in the red band and overestimates the albedo by 0.05% in the NIR band, respectively, at SZAs = 30°–70°, which may be negligible. (3) The albedo derived by these two models shows a high correlation (R2 > 0.9) between the field-measured and Polarization and Directionality of the Earth's Reflectances (POLDER) data, especially for the black-sky albedo. However, the albedo derived using the RTR model is significantly underestimated compared with the RTS model. The RTR model underestimates the black-sky albedo (white-sky albedo) retrievals by 0.62% (1.51%) and 0.93% (2.08%) in the red and NIR bands, respectively, for the field-measured data. The shortwave black-sky and white-sky albedos derived using the RTR model for the POLDER data are underestimated by 1.43% and 1.54%, respectively, compared with the RTS model. These results indicate that the snow kernel in the kernel-driven BRDF model frame is more accurate in snow albedo retrievals and has the potential for application in the field of the regional and global energy budget.

A machine learning approach to estimation of downward solar radiation from satellite-derived data products: An application over a semi-arid ecosystem in the U.S

Shortwave solar radiation is an important component of the surface energy balance and provides the principal source of energy for terrestrial ecosystems. This paper presents a machine learning approach in the form of a random forest (RF) model for estimating daily downward solar radiation flux at the land surface over complex terrain using MODIS (MODerate Resolution Imaging Spectroradiometer) remote sensing data. The model-building technique makes use of a unique network of 16 solar flux measurements in the semi-arid Reynolds Creek Experimental Watershed and Critical Zone Observatory, in southwest Idaho, USA. Based on a composite RF model built on daily observations from all 16 sites in the watershed, the model simulation of downward solar radiation matches well with the observation data (r2 = 0.96). To evaluate model performance, RF models were built from 12 of 16 sites selected at random and validated against the observations at the remaining four sites. Overall root mean square errors (RMSE), bias, and mean absolute error (MAE) are small (range: 37.17 W/m2-81.27 W/m2, -48.31 W/m2-15.67 W/m2, and 26.56 W/m2-63.77 W/m2, respectively). When extrapolated to the entire watershed, spatiotemporal patterns of solar flux are largely consistent with expected trends in this watershed. We also explored significant predictors of downward solar flux in order to reveal important properties and processes controlling downward solar radiation. Based on the composite RF model built on all 16 sites, the three most important predictors to estimate downward solar radiation include the black sky albedo (BSA) near infrared band (0.858 μm), BSA visible band (0.3-0.7 μm), and clear day coverage. This study has important implications for improving the ability to derive downward solar radiation through a fusion of multiple remote sensing datasets and can potentially capture spatiotemporally varying trends in solar radiation that is useful for land surface hydrologic and terrestrial ecosystem modeling.