Explore the vast range of titles available.

The Global Land Surface Satellite (GLASS) product suite currently contains 12 products, including leaf area index, fraction of absorbed photosynthetically active radiation, fraction of green vegetation coverage, gross primary production, broadband albedo, broadband longwave emissivity, downward shortwave radiation and photosynthetically active radiation, land surface temperature, downward and upwelling thermal radiation, all-wave net radiation, and evapotranspiration. These products are generated from the Advanced Very High Resolution Radiometer and Moderate Resolution Imaging Spectroradiometer satellite data. Their unique features include long-term temporal coverage (many from 1981 to the present), high spatial resolutions of the surface radiation products (1 km and 0.05°), spatial continuities without missing pixels, and high quality and accuracy based on extensive validation using in situ measurements and intercomparisons with other existing satellite products. Moreover, the GLASS products are based on robust algorithms that have been published in peer-reviewed literature. Herein, we provide an overview of the algorithm development, product characteristics, and some preliminary applications of these products. We also describe the next steps, such as improving the existing GLASS products, generating more climate data records (CDRs), broadening product dissemination, and fostering their wider utilization. The GLASS products are freely available to the public.

The second version of the coupled Norwegian Earth System Model (NorESM2) is presented and evaluated. NorESM2 is based on the second version of the Community Earth System Model (CESM2) and shares with CESM2 the computer code infrastructure and many Earth system model components. However, NorESM2 employs entirely different ocean and ocean biogeochemistry models. The atmosphere component of NorESM2 (CAM-Nor) includes a different module for aerosol physics and chemistry, including interactions with cloud and radiation; additionally, CAM-Nor includes improvements in the formulation of local dry and moist energy conservation, in local and global angular momentum conservation, and in the computations for deep convection and air–sea fluxes. The surface components of NorESM2 have minor changes in the albedo calculations and to land and sea-ice models.We present results from simulations with NorESM2 that were carried out for the sixth phase of the Coupled Model Intercomparison Project (CMIP6). Two versions of the model are used: one with lower (∼ 2∘) atmosphere–land resolution and one with medium (∼ 1∘) atmosphere–land resolution. The stability of the pre-industrial climate and the sensitivity of the model to abrupt and gradual quadrupling of CO2 are assessed, along with the ability of the model to simulate the historical climate under the CMIP6 forcings. Compared to observations and reanalyses, NorESM2 represents an improvement over previous versions of NorESM in most aspects. NorESM2 appears less sensitive to greenhouse gas forcing than its predecessors, with an estimated equilibrium climate sensitivity of 2.5 K in both resolutions on a 150-year time frame; however, this estimate increases with the time window and the climate sensitivity at equilibration is much higher. We also consider the model response to future scenarios as defined by selected Shared Socioeconomic Pathways (SSPs) from the Scenario Model Intercomparison Project defined under CMIP6. Under the four scenarios (SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5), the warming in the period 2090–2099 compared to 1850–1879 reaches 1.3, 2.2, 3.0, and 3.9 K in NorESM2-LM, and 1.3, 2.1, 3.1, and 3.9 K in NorESM-MM, robustly similar in both resolutions. NorESM2-LM shows a rather satisfactory evolution of recent sea-ice area. In NorESM2-LM, an ice-free Arctic Ocean is only avoided in the SSP1-2.6 scenario.

The Snow, Ice, and Aerosol Radiative (SNICAR) model has been used in various capacities over the last 15 years to model the spectral albedo of snow with light-absorbing constituents (LACs). Recent studies have extended the model to include an adding-doubling two-stream solver and representations of non-spherical ice particles; carbon dioxide snow; snow algae; and new types of mineral dust, volcanic ash, and brown carbon. New options also exist for ice refractive indices and solar-zenith-angle-dependent surface spectral irradiances used to derive broadband albedo. The model spectral range was also extended deeper into the ultraviolet for studies of extraterrestrial and high-altitude cryospheric surfaces. Until now, however, these improvements and capabilities have not been merged into a unified code base. Here, we document the formulation and evaluation of the publicly available SNICAR-ADv3 source code, web-based model, and accompanying library of constituent optical properties. The use of non-spherical ice grains, which scatter less strongly into the forward direction, reduces the simulated albedo perturbations from LACs by ∼9%–31%, depending on which of the three available non-spherical shapes are applied. The model compares very well against measurements of snow albedo from seven studies, though key properties affecting snow albedo are not fully constrained with measurements, including ice effective grain size of the top sub-millimeter of the snowpack, mixing state of LACs with respect to ice grains, and site-specific LAC optical properties. The new default ice refractive indices produce extremely high pure snow albedo (>0.99) in the blue and ultraviolet part of the spectrum, with such values only measured in Antarctica so far. More work is needed particularly in the representation of snow algae, including experimental verification of how different pigment expressions and algal cell concentrations affect snow albedo. Representations and measurements of the influence of liquid water on spectral snow albedo are also needed.

The TROPOspheric Monitoring Instrument (TROPOMI) on board the Sentinel 5 Precursor (S5-P) satellite provides methane (CH4) measurements with high accuracy and exceptional temporal and spatial resolution and sampling. TROPOMI CH4 measurements are highly valuable to constrain emissions inventories and for trend analysis, with strict requirements on the data quality. This study describes the improvements that we have implemented to retrieve CH4 from TROPOMI using the RemoTeC full-physics algorithm. The updated retrieval algorithm features a constant regularization scheme of the inversion that stabilizes the retrieval and yields less scatter in the data and includes a higher resolution surface altitude database. We have tested the impact of three state-of-the-art molecular spectroscopic databases (HITRAN 2008, HITRAN 2016 and Scientific Exploitation of Operational Missions – Improved Atmospheric Spectroscopy Databases SEOM-IAS) and found that SEOM-IAS provides the best fitting results. The most relevant update in the TROPOMI XCH4 data product is the implementation of an a posteriori correction fully independent of any reference data that is more accurate and corrects for the underestimation at low surface albedo scenes and the overestimation at high surface albedo scenes. After applying the correction, the albedo dependence is removed to a large extent in the TROPOMI versus satellite (Greenhouse gases Observing SATellite – GOSAT) and TROPOMI versus ground-based observations (Total Carbon Column Observing Network – TCCON) comparison, which is an independent verification of the correction scheme. We validate 2 years of TROPOMI CH4 data that show the good agreement of the updated TROPOMI CH4 with TCCON (−3.4 ± 5.6 ppb) and GOSAT (−10.3 ± 16.8 ppb) (mean bias and standard deviation). Low- and high-albedo scenes as well as snow-covered scenes are the most challenging for the CH4 retrieval algorithm, and although the a posteriori correction accounts for most of the bias, there is a need to further investigate the underlying cause.

The Aerosol Robotic Network (AERONET) Version 3 (V3) aerosol retrieval algorithm is described, which is based on the Version 2 (V2) algorithm with numerous updates. Comparisons of V3 aerosol retrievals to those of V2 are presented, along with a new approach to estimate uncertainties in many of the retrieved aerosol parameters. Changes in the V3 aerosol retrieval algorithm include (1) a new polarized radiative transfer code (RTC), which replaced the scalar RTC of V2, (2) detailed characterization of gas absorption by adding NO2 and H2O to specify total gas absorption in the atmospheric column, specification of vertical profiles of all the atmospheric species, (3) new bidirectional reflectance distribution function (BRDF) parameters for land sites adopted from the MODIS BRDF/Albedo product, (4) a new version of the extraterrestrial solar flux spectrum, and (5) a new temperature correction procedure of both direct Sun and sky radiance measurements. The potential effect of each change in V3 on single scattering albedo (SSA) retrievals was analyzed. The operational almucantar retrievals of V2 versus V3 were compared for four AERONET sites: GSFC, Mezaira, Mongu, and Kanpur. Analysis showed very good agreement in retrieved parameters of the size distributions. Comparisons of SSA retrievals for dust aerosols (Mezaira) showed a good agreement in 440 nm SSA, while for longer wavelengths V3 SSAs are systematically higher than those of V2, with the largest mean difference at 675 nm due to cumulative effects of both extraterrestrial solar flux and BRDF changes. For non-dust aerosols, the largest SSA deviation is at 675 nm due to differences in extraterrestrial solar flux spectrums used in each version. Further, the SSA 675 nm mean differences are very different for weakly (GSFC) and strongly (Mongu) absorbing aerosols, which is explained by the lower sensitivity to a bias in aerosol scattering optical depth by less absorbing aerosols. A new hybrid (HYB) sky radiance measurement scan is introduced and discussed. The HYB combines features of scans in two different planes to maximize the range of scattering angles and achieve scan symmetry, thereby allowing for cloud screening and spatial averaging, which is an advantage over the principal plane scan that lacks robust symmetry. We show that due to an extended range of scattering angles, HYB SSA retrievals for dust aerosols exhibit smaller variability with solar zenith angles (SZAs) than those of almucantar (ALM), which allows extension of HYB SSA retrievals to SZAs less than 50∘ to as small as 25∘. The comparison of SSA retrievals from closely time-matched HYB and ALM scans in the 50 to 75∘ SZA range showed good agreement with the differences below ∼0.005. We also present an approach to estimate retrieval uncertainties which utilizes the variability in retrieved parameters generated by perturbing both measurements and auxiliary input parameters as a proxy for retrieval uncertainty. The perturbations in measurements and auxiliary inputs are assumed as estimated biases in aerosol optical depth (AOD), radiometric calibration of sky radiances combined with solar spectral irradiance, and surface reflectance. For each set of Level 2 Sun/sky radiometer observations, 27 inputs corresponding to 27 combinations of biases were produced and separately inverted to generate the following statistics of the inversion results: average, standard deviation, minimum and maximum values. From these statistics, standard deviation (labeled U27) is used as a proxy for estimated uncertainty, and a lookup table (LUT) approach was implemented to reduce the computational time. The U27 climatological LUT was generated from the entire AERONET almucantar (1993–2018) and hybrid (2014–2018) scan databases by binning U27s in AOD (440 nm), Angström exponent (AE, 440–870 nm), and SSA (440, 675, 870, 1020 nm). Using this LUT approach, the uncertainty estimates U27 for each individual V3 Level 2 retrieval can be obtained by interpolation using the corresponding measured and inverted combination of AOD, AE, and SSA.

Sea ice cools Earth by reducing its absorbed solar energy. We combine radiative transfer modeling with satellite‐derived surface albedo, sea ice, and cloud distributions to quantify the top‐of‐atmosphere sea ice radiative effect (SIRE). Averaged over 1980–2023, Arctic and Antarctic SIREs range from −0.64 to −0.86 W m−2 and −0.85 to −0.98 W m−2, respectively, with different cloud data sets and assumptions of climatological versus annually‐varying clouds. SIRE trends, however, are relatively insensitive to these assumptions. Arctic SIRE has weakened quasi‐linearly at a rate of 0.04–0.05 W m−2 decade−1, implying a 21%–27% reduction in the reflective power of Arctic sea ice since 1980. Antarctic sea ice exhibited a regime change in 2016, resulting in 2016–2023 Antarctic and global SIRE being 0.08–0.12 and 0.22–0.27 W m−2 weaker, respectively, relative to 1980–1988. Global sea ice has therefore lost 13%–15% of its planetary cooling effect since the early/mid 1980s, and the implied global sea ice albedo feedback is 0.24–0.38 W m−2 K−1. Plain Language Summary Ice that forms on the surface of the ocean (sea ice) is highly reflective of sunlight, especially in comparison with ocean water, which has very low reflectivity. Sea ice therefore cools Earth by decreasing the amount of sunlight that it absorbs. We present a measure of this planetary cooling effect and evaluate how it has changed since 1980, around the advent of consistent satellite observations of Earth's surface and atmosphere. At individual locations and times, the cooling effect depends on the amount of incoming sunlight (which varies strongly with season in polar regions), cloud coverage (which masks the underlying surface from sunlight), and the reflectivity of the sea ice and overlying snow. Averaged globally, this effect varies with the areal coverage of sea ice and has therefore weakened with shrinking ice coverage. The planetary cooling effects of Arctic and Antarctic sea ice during 2016–2023 were about 20% and 12% less, respectively, than they were during 1980–1988. Disappearing sea ice is therefore amplifying climate change by causing Earth to absorb roughly an additional 0.3 W m−2 of solar power for each degree Celsius of global warming, a feedback that is stronger than that simulated by most climate models. Key Points We quantify the top‐of‐atmosphere radiative effect of global sea ice with historical surface albedo, sea ice, and cloud data sets The Arctic sea ice radiative effect has weakened at 0.04–0.05 W m−2decade−1, or by about 24%, since 1980 The planetary cooling effect of global sea ice was about 0.25 W m−2 (14%) weaker during 2016–2023 than during 1980–1988

We implement a set of new parameterizations into the widely used Snow, Ice, and Aerosol Radiative (SNICAR) model to account for effects of snow grain shape (spherical vs. nonspherical) and black carbon (BC)–snow mixing state (external vs. internal). We find that nonspherical snow grains lead to higher pure albedo but weaker BC-induced albedo reductions relative to spherical snow grains, while BC–snow internal mixing significantly enhances albedo reductions relative to external mixing. The combination of snow nonsphericity and internal mixing suggests an important interactive effect on BC-induced albedo reduction. Comparisons with observations of clean and BC-contaminated snow albedo show that model simulations accounting for both snow nonsphericity and BC–snow internal mixing perform better than those using the common assumption of spherical snow grains and external mixing. We further apply the updated SNICAR model with comprehensive in situ measurements of BC concentrations in the Tibetan Plateau snowpack to quantify the present-day (2000–2015) BC-induced snow albedo effects from a regional and seasonal perspective. The BC concentrations show distinct and substantial sub-regional and seasonal variations, with higher values in the non-monsoon season and low altitudes. As a result, the BC-induced regional mean snow albedo reductions and surface radiative effects vary by up to an order of magnitude across different sub-regions and seasons, with values of 0.7–30.7 and 1.4–58.4 W m−2 for BC externally mixed with fresh and aged snow spheres, respectively. The BC radiative effects are further complicated by uncertainty in snow grain shape and BC–snow mixing state. BC–snow internal mixing enhances the mean albedo effects over the plateau by 30–60 % relative to external mixing, while nonspherical snow grains decrease the mean albedo effects by up to 31 % relative to spherical grains. Based on this study, extensive measurements and improved model characterization of snow grain shape and aerosol–snow mixing state are urgently needed in order to precisely evaluate BC–snow albedo effects.

This paper describes the latest version of the algorithm MAIAC (Multi-Angle Implementation of Atmospheric Correction) used for processing the MODIS (Moderate-resolution Imaging Spectroradiometer) Collection6 data record. Since initial publication in 2011-2012, MAIAC has changed considerably to adapt to global processing and improve cloud/snow detection, aerosol retrievals and atmospheric correction of MODIS data. The main changes include (1) transition from a 25 to 1 km scale for retrieval of the spectral regression coefficient (SRC) which helped to remove occasional blockiness at 25 km scale in the aerosol optical depth (AOD) and in the surface reflectance, (2) continuous improvements of cloud detection, (3) introduction of smoke and dust tests to discriminate absorbing fine- and coarse mode aerosols, (4) adding over-water processing, (5) general optimization of the LUT (LookUp-Table)-based radiative transfer for the global processing, and others. MAIAC provides an interdisciplinary suite of atmospheric and land products, including cloud mask (CM), column water vapor (CWV), AOD at 0.47 and 0.55 m, aerosol type (background, smoke or dust) and fine-mode fraction over water; spectral bidirectional reflectance factors (BRF), parameters of Ross-thick Lisparse (RTLS) bidirectional reflectance distribution function (BRDF) model and instantaneous albedo. For snow-covered surfaces, we provide subpixel snow fraction and snow grain size. All products come in standard HDF4 (software library) format at 1 km resolution, except for BRF, which is also provided at 500 m resolution on a sinusoidal grid adopted by the MODIS Land team. All products are provided on per-observation basis in daily files except for the BRDF/Albedo product, which is reported every 8 days. Because MAIAC uses a time series approach, BRDF/Albedo is naturally gap-filled over land where missing values are filled-in with results from the previous retrieval. While the BRDF model is reported for MODIS Land bands 1-7 and ocean band 8, BRF is reported for both land and ocean bands 1-12. This paper focuses on MAIAC cloud detection, aerosol retrievals and atmospheric correction and describes MCD19 data products and quality assurance (QA) flags.

A 15-yr duration record of mooring observations from the eastern (>70°E) Eurasian Basin (EB) of the Arctic Ocean is used to show and quantify the recently increased oceanic heat flux from intermediate-depth (~150–900 m) warm Atlantic Water (AW) to the surface mixed layer and sea ice. The upward release of AW heat is regulated by the stability of the overlying halocline, which we show has weakened substantially in recent years. Shoaling of the AW has also contributed, with observations in winter 2017–18 showing AW at only 80 m depth, just below the wintertime surface mixed layer, the shallowest in our mooring records. The weakening of the halocline for several months at this time implies that AW heat was linked to winter convection associated with brine rejection during sea ice formation. This resulted in a substantial increase of upward oceanic heat flux during the winter season, from an average of 3–4 W m−2 in 2007–08 to >10 W m−2 in 2016–18. This seasonal AW heat loss in the eastern EB is equivalent to a more than a twofold reduction of winter ice growth. These changes imply a positive feedback as reduced sea ice cover permits increased mixing, augmenting the summer-dominated ice-albedo feedback.

The second edition of the satellite-derived climate data record CLARA (The CM SAF Cloud, Albedo And Surface Radiation dataset from AVHRR data – second edition denoted as CLARA-A2) is described. The data record covers the 34-year period from 1982 until 2015 and consists of cloud, surface albedo and surface radiation budget products derived from the AVHRR (Advanced Very High Resolution Radiometer) sensor carried by polar-orbiting, operational meteorological satellites. The data record is produced by the EUMETSAT Climate Monitoring Satellite Application Facility (CM SAF) project as part of the operational ground segment. Its upgraded content and methodology improvements since edition 1 are described in detail, as are some major validation results. Some of the main improvements to the data record come from a major effort in cleaning and homogenizing the basic AVHRR level-1 radiance record and a systematic use of CALIPSO-CALIOP cloud information for development and validation purposes. Examples of applications studying decadal changes in Arctic summer surface albedo and cloud conditions are provided.