Atmospheric model systems, such as those used for weather forecast and reanalysis production, often have significant and systematic errors in their representation of the Arctic surface energy budget and its components. The newly available observation data of the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition (2019/2020) enable a range of model analyses and validation in order to advance our understanding of potential model deficiencies. In the present study, we analyze deficiencies in the surface radiative energy budget over Arctic sea ice in the ERA5 global atmospheric reanalysis by comparing against the winter MOSAiC campaign data, as well as, a pan-Arctic level-2 MODIS ice surface temperature remote sensing product. We find that ERA5 can simulate the timing of radiatively clear periods, though it is not able to distinguish the two observed radiative Arctic winter states, radiatively clear and opaquely cloudy, in the distribution of the net surface radiative budget. The ERA5 surface temperature over Arctic sea ice has a conditional error with a positive bias in radiatively clear conditions and a negative bias in opaquely cloudy conditions. The mean surface temperature error is 4°C for radiatively clear situations at MOSAiC and up to 15°C in some parts of the Arctic. The spatial variability of the surface temperature, given by 4 observation sites at MOSAiC, is not captured by ERA5 due to its spatial resolution but represented in the level-2 satellite product. The sensitivity analysis of possible error sources, using satellite products of snow depth and sea ice thickness, shows that the positive surface temperature errors during radiatively clear events are, to a large extent, caused by insufficient sea ice thickness and snow depth representation in the reanalysis system. A positive bias characterizes regions with ice thickness greater than 1.5 m, while the negative bias for thinner ice is partly compensated by the effect of snow.
Atmospheric reanalyses are widely used to estimate the past atmospheric near-surface state over sea ice. They provide boundary conditions for sea ice and ocean numerical simulations and relevant information for studying polar variability and anthropogenic climate change. Previous research revealed the existence of large near-surface temperature biases (mostly warm) over the Arctic sea ice in the current generation of atmospheric reanalyses, which is linked to a poor representation of the snow over the sea ice and the stably stratified boundary layer in the forecast models used to produce the reanalyses. These errors can compromise the employment of reanalysis products in support of polar research. Here, we train a fully connected neural network that learns from remote sensing infrared temperature observations to correct the existing generation of uncoupled atmospheric reanalyses (ERA5, JRA-55) based on a set of sea ice and atmospheric predictors, which are themselves reanalysis products. The advantages of the proposed correction scheme over previous calibration attempts are the consideration of the synoptic weather and cloud state, compatibility of the predictors with the mechanism responsible for the bias, and a self-emerging seasonality and multidecadal trend consistent with the declining sea ice state in the Arctic. The correction leads on average to a 27% temperature bias reduction for ERA5 and 7% for JRA-55 if compared to independent in situ observations from the MOSAiC campaign (respectively, 32% and 10% under clear-sky conditions). These improvements can be beneficial for forced sea ice and ocean simulations, which rely on reanalyses surface fields as boundary conditions.
This study illustrates a novel method based on machine learning for reducing the systematic surface temperature errors that characterize multiple atmospheric reanalyses in sea ice–covered regions of the Arctic under clear-sky conditions. The correction applied to the temperature field is consistent with the local weather and the sea ice and snow conditions, meaning that it responds to seasonal changes in sea ice cover as well as to its long-term decline due to global warming. The corrected reanalysis temperature can be employed to support polar research activities, and in particular to better simulate the evolution of the interacting sea ice and ocean system within numerical models.
- NSF-PAR ID:
- 10415612
- Publisher / Repository:
- American Meteorological Society
- Date Published:
- Journal Name:
- Monthly Weather Review
- Volume:
- 151
- Issue:
- 6
- ISSN:
- 0027-0644
- Page Range / eLocation ID:
- p. 1443-1458
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
more » « less
Abstract The Iceland and Greenland Seas are a crucial region for the climate system, being the headwaters of the lower limb of the Atlantic Meridional Overturning Circulation. Investigating the atmosphere–ocean–ice processes in this region often necessitates the use of meteorological reanalyses—a representation of the atmospheric state based on the assimilation of observations into a numerical weather prediction system. Knowing the quality of reanalysis products is vital for their proper use. Here we evaluate the surface‐layer meteorology and surface turbulent fluxes in winter and spring for the latest reanalysis from the European Centre for Medium‐Range Weather Forecasts, i.e., ERA5. In situ observations from a meteorological buoy, a research vessel, and a research aircraft during the Iceland–Greenland Seas Project provide unparalleled coverage of this climatically important region. The observations are independent of ERA5. They allow a comprehensive evaluation of the surface meteorology and fluxes of these subpolar seas and, for the first time, a specific focus on the marginal ice zone. Over the ice‐free ocean, ERA5 generally compares well to the observations of surface‐layer meteorology and turbulent fluxes. However, over the marginal ice zone, the correspondence is noticeably less accurate: for example, the root‐mean‐square errors are significantly higher for surface temperature, wind speed, and surface sensible heat flux. The primary reason for the difference in reanalysis quality is an overly smooth sea‐ice distribution in the surface boundary conditions used in ERA5. Particularly over the marginal ice zone, unrepresented variability and uncertainties in how to parameterize surface exchange compromise the quality of the reanalyses. A parallel evaluation of higher‐resolution forecast fields from the Met Office's Unified Model corroborates these findings.
-
A Lagrangian Snow‐Evolution System for Sea‐Ice Applications (SnowModel‐LG): Part I—Model Descriptionmore » « less
Abstract A Lagrangian snow‐evolution model (SnowModel‐LG) was used to produce daily, pan‐Arctic, snow‐on‐sea‐ice, snow property distributions on a 25 × 25‐km grid, from 1 August 1980 through 31 July 2018 (38 years). The model was forced with NASA's Modern Era Retrospective‐Analysis for Research and Applications‐Version 2 (MERRA‐2) and European Centre for Medium‐Range Weather Forecasts (ECMWF) ReAnalysis‐5th Generation (ERA5) atmospheric reanalyses, and National Snow and Ice Data Center (NSIDC) sea ice parcel concentration and trajectory data sets (approximately 61,000, 14 × 14‐km parcels). The simulations performed full surface and internal energy and mass balances within a multilayer snowpack evolution system. Processes and features accounted for included rainfall, snowfall, sublimation from static‐surfaces and blowing‐snow, snow melt, snow density evolution, snow temperature profiles, energy and mass transfers within the snowpack, superimposed ice, and ice dynamics. The simulations produced horizontal snow spatial structures that likely exist in the natural system but have not been revealed in previous studies spanning these spatial and temporal domains. Blowing‐snow sublimation made a significant contribution to the snowpack mass budget. The superimposed ice layer was minimal and decreased over the last four decades. Snow carryover to the next accumulation season was minimal and sensitive to the melt‐season atmospheric forcing (e.g., the average summer melt period was 3 weeks or 50% longer with ERA5 forcing than MERRA‐2 forcing). Observed ice dynamics controlled the ice parcel age (in days), and ice age exerted a first‐order control on snow property evolution.
-
more » « less
Abstract Convectively coupled equatorial waves are a significant source of atmospheric variability in the tropics. Current numerical models continue to struggle in simulating the coupled diabatic heating fields that are responsible for the development and maintenance of these waves. This study investigates how the diabatic fields associated with Mixed Rossby–Gravity waves (MRGs) are represented in four reanalysis products by using a unique observational dataset from the TRMM‐KWAJEX (Tropical Rainfall Measuring Mission‐Kwajalein Experiment) field campaign. These reanalyses include ERA5, Japanese 55‐year Reanalysis (JRA‐55), Climate Forecast System Reanalysis (CFSR), and Modern‐Era Retrospective Analysis for Research and Applications (MERRA). We found that all four reanalyses captured the MRG structures in winds and temperature, and to a lesser degree in the humidity field except in the boundary layer. However, only the ERA5 and MERRA reanalyses captured the gradual rise and succession of the diabatic heating from boundary layer turbulence, shallow convection, cumulus congestus, and deep convection within the waves. ERA5 is the only product that also captured the gradual rise of the subgrid‐scale vertical transport of moist static energy. All reanalysis products underestimated the diabatic heating from cumulus congestus. Results provide observational basis on what aspects of MRG can be trusted and what cannot in the reanalysis products.
-
more » « less
Abstract The surface downwelling longwave radiation component (LW↓) is crucial for the determination of the surface energy budget and has significant implications for the resilience of ice surfaces in the polar regions. Accurate model evaluation of this radiation component requires knowledge about the phase, vertical distribution, and associated temperature of water in the atmosphere, all of which control the LW↓ signal measured at the surface. In this study, we examine the LW↓ model errors found in the Antarctic Mesoscale Prediction System (AMPS) operational forecast model and the ERA5 model relative to observations from the ARM West Antarctic Radiation Experiment (AWARE) campaign at McMurdo Station and the West Antarctic Ice Sheet (WAIS) Divide. The errors are calculated separately for observed clear-sky conditions, ice-cloud occurrences, and liquid-bearing cloud-layer (LBCL) occurrences. The analysis results show a tendency in both models at each site to underestimate the LW↓ during clear-sky conditions, high error variability (standard deviations > 20 W m−2) during any type of cloud occurrence, and negative LW↓ biases when LBCLs are observed (bias magnitudes >15 W m−2 in tenuous LBCL cases and >43 W m−2 in optically thick/opaque LBCLs instances). We suggest that a generally dry and liquid-deficient atmosphere responsible for the identified LW↓ biases in both models is the result of excessive ice formation and growth, which could stem from the model initial and lateral boundary conditions, microphysics scheme, aerosol representation, and/or limited vertical resolution.
