skip to main content


Title: Assessing uncertainties in sea ice extent climate indicators
The uncertainties in sea ice extent (total area covered by sea ice with concentration>15%) derived from passive microwave sensors are assessed in two ways. Absolute uncertainty (accuracy) is evaluated based on the comparison of the extent between several products. There are clear biases between the extent from the different products that are of the order of 500 000 to 1×106 km2 depending on the season and hemisphere. These biases are due to differences in the algorithm sensitivity to ice edge conditions and the spatial resolution of different sensors. Relative uncertainty is assessed by examining extents from the National Snow and Ice Data Center Sea Ice Index product. The largest source of uncertainty,∼100 000 km2, is between near-real-time and final products due to different input source data and different processing and quality control. For consistent processing, the uncertainty is assessed using different input source data and by varying concentration algorithm parameters. This yields a relative uncertainty of 30 000–70 000 km2. The Arctic minimum extent uncertainty is∼40 000 km2. Uncertainties in comparing with earlier parts of the record may be higher due to sensor transitions. For the first time, this study provides a quantitative estimate of sea ice extent uncertainty.  more » « less
Award ID(s):
1748953
NSF-PAR ID:
10128750
Author(s) / Creator(s):
Date Published:
Journal Name:
Environmental letters
Volume:
14
ISSN:
0013-9300
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Ocean colour is recognised as an Essential Climate Variable (ECV) by the Global Climate Observing System (GCOS); and spectrally-resolved water-leaving radiances (or remote-sensing reflectances) in the visible domain, and chlorophyll-a concentration are identified as required ECV products. Time series of the products at the global scale and at high spatial resolution, derived from ocean-colour data, are key to studying the dynamics of phytoplankton at seasonal and inter-annual scales; their role in marine biogeochemistry; the global carbon cycle; the modulation of how phytoplankton distribute solar-induced heat in the upper layers of the ocean; and the response of the marine ecosystem to climate variability and change. However, generating a long time series of these products from ocean-colour data is not a trivial task: algorithms that are best suited for climate studies have to be selected from a number that are available for atmospheric correction of the satellite signal and for retrieval of chlorophyll-a concentration; since satellites have a finite life span, data from multiple sensors have to be merged to create a single time series, and any uncorrected inter-sensor biases could introduce artefacts in the series, e.g., different sensors monitor radiances at different wavebands such that producing a consistent time series of reflectances is not straightforward. Another requirement is that the products have to be validated against in situ observations. Furthermore, the uncertainties in the products have to be quantified, ideally on a pixel-by-pixel basis, to facilitate applications and interpretations that are consistent with the quality of the data. This paper outlines an approach that was adopted for generating an ocean-colour time series for climate studies, using data from the MERIS (MEdium spectral Resolution Imaging Spectrometer) sensor of the European Space Agency; the SeaWiFS (Sea-viewing Wide-Field-of-view Sensor) and MODIS-Aqua (Moderate-resolution Imaging Spectroradiometer-Aqua) sensors from the National Aeronautics and Space Administration (USA); and VIIRS (Visible and Infrared Imaging Radiometer Suite) from the National Oceanic and Atmospheric Administration (USA). The time series now covers the period from late 1997 to end of 2018. To ensure that the products meet, as well as possible, the requirements of the user community, marine-ecosystem modellers, and remote-sensing scientists were consulted at the outset on their immediate and longer-term requirements as well as on their expectations of ocean-colour data for use in climate research. Taking the user requirements into account, a series of objective criteria were established, against which available algorithms for processing ocean-colour data were evaluated and ranked. The algorithms that performed best with respect to the climate user requirements were selected to process data from the satellite sensors. Remote-sensing reflectance data from MODIS-Aqua, MERIS, and VIIRS were band-shifted to match the wavebands of SeaWiFS. Overlapping data were used to correct for mean biases between sensors at every pixel. The remote-sensing reflectance data derived from the sensors were merged, and the selected in-water algorithm was applied to the merged data to generate maps of chlorophyll concentration, inherent optical properties at SeaWiFS wavelengths, and the diffuse attenuation coefficient at 490 nm. The merged products were validated against in situ observations. The uncertainties established on the basis of comparisons with in situ data were combined with an optical classification of the remote-sensing reflectance data using a fuzzy-logic approach, and were used to generate uncertainties (root mean square difference and bias) for each product at each pixel. 
    more » « less
  2. Accurate multidecadal radiative flux records are vital to understand Arctic amplification and constrain climate model uncertainties. Uncertainty in the NASA Clouds and the Earth’s Radiant Energy System (CERES)-derived irradiances is larger over sea ice than any other surface type and comes from several sources. The year-long Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition in the central Arctic provides a rare opportunity to explore uncertainty in CERES-derived radiative fluxes. First, a systematic and statistically robust assessment of surface shortwave and longwave fluxes was conducted using in situ measurements from MOSAiC flux stations. The CERES Synoptic 1degree (SYN1deg) product overestimates the downwelling shortwave flux by +11.40 Wm–2 and underestimates the upwelling shortwave flux by –15.70 Wm–2 and downwelling longwave fluxes by –12.58 Wm–2 at the surface during summer. In addition, large differences are found in the upwelling longwave flux when the surface approaches the melting point (approximately 0°C). The biases in downwelling shortwave and longwave fluxes suggest that the atmosphere represented in CERES is too optically thin. The large negative bias in upwelling shortwave flux can be attributed in large part to lower surface albedo (–0.15) in satellite footprint relative to surface sensors. Additionally, the results show that the spectral surface albedo used in SYN1deg overestimates albedo in visible and mid-infrared bands. A series of radiative transfer model perturbation experiments are performed to quantify the factors contributing to the differences. The CERES-MOSAiC broadband albedo differences (approximately 20 Wm–2) explain a larger portion of the upwelling shortwave flux difference than the spectral albedo shape differences (approximately 3 Wm–2). In addition, the differences between perturbation experiments using hourly and monthly MOSAiC surface albedo suggest that approximately 25% of the sea ice surface albedo variability is explained by factors not correlated with daily sea ice concentration variability. Biases in net shortwave and longwave flux can be reduced to less than half by adjusting both albedo and cloud inputs toward observed values. The results indicate that improvements in the surface albedo and cloud data would substantially reduce the uncertainty in the Arctic surface radiation budget derived from CERES data products. 
    more » « less
  3. Abstract. The role of clouds in the Arctic radiation budget is not well understood. Ground-based and airborne measurements provide valuable data to test and improve our understanding. However, the ground-based measurements are intrinsically sparse, and the airborne observations are snapshots in time and space. Passive remote sensing measurements from satellite sensors offer high spatial coverage and an evolving time series, having lengths potentially of decades. However, detecting clouds by passive satellite remote sensing sensors is challenging over the Arctic because of the brightness of snow and ice in the ultraviolet and visible spectral regions and because of the small brightness temperature contrast to the surface. Consequently, the quality of the resulting cloud data products needs to be assessed quantitatively. In this study, we validate the cloud data products retrieved from the Advanced Very High Resolution Radiometer (AVHRR) post meridiem (PM) data from the polar-orbiting NOAA-19 satellite and compare them with those derived from the ground-based instruments during the sunlit months. The AVHRR cloud data products by the European Space Agency (ESA) Cloud Climate Change Initiative (Cloud_CCI) project uses the observations in the visible and IR bands to determine cloud properties. The ground-based measurements from four high-latitude sites have been selected for this investigation: Hyytiälä (61.84∘ N, 24.29∘ E), North Slope of Alaska (NSA; 71.32∘ N, 156.61∘ W), Ny-Ålesund (Ny-Å; 78.92∘ N, 11.93∘ E), and Summit (72.59∘ N, 38.42∘ W). The liquid water path (LWP) ground-based data are retrieved from microwave radiometers, while the cloud top height (CTH) has been determined from the integrated lidar–radar measurements. The quality of the satellite products, cloud mask and cloud optical depth (COD), has been assessed using data from NSA, whereas LWP and CTH have been investigated over Hyytiälä, NSA, Ny-Å, and Summit. The Cloud_CCI COD results for liquid water clouds are in better agreement with the NSA radiometer data than those for ice clouds. For liquid water clouds, the Cloud_CCI COD is underestimated roughly by 3 optical depth (OD) units. When ice clouds are included, the underestimation increases to about 5 OD units. The Cloud_CCI LWP is overestimated over Hyytiälä by ≈7 g m−2, over NSA by ≈16 g m−2, and over Ny-Å by ≈24 g m−2. Over Summit, CCI LWP is overestimated for values ≤20 g m−2 and underestimated for values >20 g m−2. Overall the results of the CCI LWP retrievals are within the ground-based instrument uncertainties. To understand the effects of multi-layer clouds on the CTH retrievals, the statistics are compared between the single-layer clouds and all types (single-layer + multi-layer). For CTH retrievals, the Cloud_CCI product overestimates the CTH for single-layer clouds. When the multi-layer clouds are included (i.e., all types), the observed CTH overestimation becomes an underestimation of about 360–420 m. The CTH results over Summit station showed the highest biases compared to the other three sites. To understand the scale-dependent differences between the satellite and ground-based data, the Bland–Altman method is applied. This method does not identify any scale-dependent differences for all the selected cloud parameters except for the retrievals over the Summit station. In summary, the Cloud_CCI cloud data products investigated agree reasonably well with those retrieved from ground-based measurements made at the four high-latitude sites.

     
    more » « less
  4. Abstract

    We outline a new and improved uncertainty analysis for the Goddard Institute for Space Studies Surface Temperature product version 4 (GISTEMP v4). Historical spatial variations in surface temperature anomalies are derived from historical weather station data and ocean data from ships, buoys, and other sensors. Uncertainties arise from measurement uncertainty, changes in spatial coverage of the station record, and systematic biases due to technology shifts and land cover changes. Previously published uncertainty estimates for GISTEMP included only the effect of incomplete station coverage. Here, we update this term using currently available spatial distributions of source data, state‐of‐the‐art reanalyses, and incorporate independently derived estimates for ocean data processing, station homogenization, and other structural biases. The resulting 95% uncertainties are near 0.05 °C in the global annual mean for the last 50 years and increase going back further in time reaching 0.15 °C in 1880. In addition, we quantify the benefits and inherent uncertainty due to the GISTEMP interpolation and averaging method. We use the total uncertainties to estimate the probability for each record year in the GISTEMP to actually be the true record year (to that date) and conclude with 87% likelihood that 2016 was indeed the hottest year of the instrumental period (so far).

     
    more » « less
  5. Sea level rise (SLR) may impose substantial economic costs to coastal communities worldwide, but characterizing its global impact remains challenging because SLR costs depend heavily on natural characteristics and human investments at each location – including topography, the spatial distribution of assets, and local adaptation decisions. To date, several impact models have been developed to estimate the global costs of SLR. Yet, the limited availability of open-source and modular platforms that easily ingest up-to-date socioeconomic and physical data sources restricts the ability of existing systems to incorporate new insights transparently. In this paper, we present a modular, open-source platform designed to address this need, providing end-to-end transparency from global input data to a scalable least-cost optimization framework that estimates adaptation and net SLR costs for nearly 10 000 global coastline segments and administrative regions. Our approach accounts both for uncertainty in the magnitude of global mean sea level (g.m.s.l.) rise and spatial variability in local relative sea level rise. Using this platform, we evaluate costs across 230 possible socioeconomic and SLR trajectories in the 21st century. According to the latest Intergovernmental Panel on Climate Change Assessment Report (AR6), g.m.s.l. is likely to rise during the 21st century by 0.40–0.69 m if late-century warming reaches 2 ∘C and by 0.58–0.91 m with 4 ∘C of warming (Fox-Kemper et al., 2021). With no forward-looking adaptation, we estimate that annual costs of sea level rise associated with a 2 ∘C scenario will likely fall between USD 1.2 and 4.0 trillion (0.1 % and 1.2 % of GDP, respectively) by 2100, depending on socioeconomic and sea level rise trajectories. Cost-effective, proactive adaptation would provide substantial benefits, lowering these values to between USD 110 and USD 530 billion (0.02 and 0.06 %) under an optimal adaptation scenario. For the likely SLR trajectories associated with 4 ∘C warming, these costs range from USD 3.1 to 6.9 trillion (0.3 % and 2.0 %) with no forward-looking adaptation and USD 200 billion to USD 750 billion (0.04 % to 0.09 %) under optimal adaptation. The Intergovernmental Panel on Climate Change (IPCC) notes that deeply uncertain physical processes like marine ice cliff instability could drive substantially higher global sea level rise, potentially approaching 2.0 m by 2100 in very high emission scenarios. Accordingly, we also model the impacts of 1.5 and 2.0 m g.m.s.l. rises by 2100; the associated annual cost estimates range from USD 11.2 to 30.6 trillion (1.2 % and 7.6 %) under no forward-looking adaptation and USD 420 billion to 1.5 trillion (0.08 % to 0.20 %) under optimal adaptation. Our modeling platform used to generate these estimates is publicly available in an effort to spur research collaboration and support decision-making, with segment-level physical and socioeconomic input characteristics provided at https://doi.org/10.5281/zenodo.7693868 (Bolliger et al., 2023a) and model results at https://doi.org/10.5281/zenodo.7693869 (Bolliger et al., 2023b).

     
    more » « less