@cmcc.it
Ocean modeling and Data Assimilation Section
Euro-Mediterranean Center on Climate Chnage
Geophysics, Modeling and Simulation
Scopus Publications
Scholar Citations
Scholar h-index
Scholar i10-index
Lorenzo Zampieri, David Clemens‐Sewall, Anne Sledd, Nils Hutter, and Marika Holland
American Geophysical Union (AGU)
AbstractModels struggle to accurately simulate observed sea ice thickness changes, which could be partially due to inadequate representation of thermodynamic processes. We analyzed co‐located winter observations of the Arctic sea ice from the Multidisciplinary Drifting Observatory for the Study of the Arctic Climate for evaluating and improving thermodynamic processes in sea ice models, aiming to enable more accurate predictions of the warming climate system. We model the sea ice and snow heat conduction for observed transects forced by realistic boundary conditions to understand the impact of the non‐resolved meter‐scale snow and sea ice thickness heterogeneity on horizontal heat conduction. Neglecting horizontal processes causes underestimating the conductive heat flux of 10% or more. Furthermore, comparing model results to independent temperature observations reveals a ∼5 K surface temperature overestimation over ice thinner than 1 m, attributed to shortcomings in parameterizing surface turbulent and radiative fluxes rather than the conduction. Assessing the model deficiencies and parameterizing these unresolved processes is required for improved sea ice representation.
Clare Eayrs, Won Sang Lee, Emilia Jin, Jean-François Lemieux, François Massonnet, Martin Vancoppenolle, Lorenzo Zampieri, Luke G. Bennetts, Ed Blockley, Eui-Seok Chung,et al.
American Meteorological Society
Lorenzo Zampieri, Gabriele Arduini, Marika Holland, Sarah P. E. Keeley, Kristian Mogensen, Matthew D. Shupe, and Steffen Tietsche
American Meteorological Society
Abstract 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. Significance Statement 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.
Longjiang Mu, Lars Nerger, Jan Streffing, Qi Tang, Bimochan Niraula, Lorenzo Zampieri, Svetlana N. Loza, and Helge F. Goessling
American Geophysical Union (AGU)
A new version of the AWI Coupled Prediction System is developed based on the Alfred Wegener Institute Climate Model v3.0. Both the ocean and the atmosphere models are upgraded or replaced, reducing the computation time by a factor of 5 at a given resolution. This allowed us to increase the ensemble size from 12 to 30, maintaining a similar resolution in both model components. The online coupled data assimilation scheme now additionally utilizes sea‐surface salinity and sea‐level anomaly as well as temperature and salinity profile observations. Results from the data assimilation demonstrate that the sea‐ice and ocean states are reasonably constrained. In particular, the temperature and salinity profile assimilation has mitigated systematic errors in the deeper ocean, although issues remain over polar regions where strong atmosphere‐ocean‐ice interaction occurs. One‐year‐long sea‐ice forecasts initialized on 1 January, 1 April, 1 July and 1 October from 2003 to 2019 are described. To correct systematic forecast errors, sea‐ice concentration from 2011 to 2019 is calibrated by trend‐adjusted quantile mapping using the preceding forecasts from 2003 to 2010. The sea‐ice edge raw forecast skill is within the range of operational global subseasonal‐to‐seasonal forecast systems, outperforming a climatological benchmark for about 2 weeks in the Arctic and about 3 weeks in the Antarctic. The calibration is much more effective in the Arctic: Calibrated sea‐ice edge forecasts outperform climatology for about 45 days in the Arctic but only 27 days in the Antarctic. Both the raw and the calibrated forecast skill exhibit strong seasonal variations.
Jan Streffing, Dmitry Sidorenko, Tido Semmler, Lorenzo Zampieri, Patrick Scholz, Miguel Andrés-Martínez, Nikolay Koldunov, Thomas Rackow, Joakim Kjellsson, Helge Goessling,et al.
Copernicus GmbH
Abstract. We developed a new version of the Alfred Wegener Institute Climate Model (AWI-CM3), which has higher skills in representing the observed climatology and better computational efficiency than its predecessors. Its ocean component FESOM2 (Finite-volumE Sea ice–Ocean Model) has the multi-resolution functionality typical of unstructured-mesh models while still featuring a scalability and efficiency similar to regular-grid models. The atmospheric component OpenIFS (CY43R3) enables the use of the latest developments in the numerical-weather-prediction community in climate sciences. In this paper we describe the coupling of the model components and evaluate the model performance on a variable-resolution (25–125 km) ocean mesh and a 61 km atmosphere grid, which serves as a reference and starting point for other ongoing research activities with AWI-CM3. This includes the exploration of high and variable resolution and the development of a full Earth system model as well as the creation of a new sea ice prediction system. At this early development stage and with the given coarse to medium resolutions, the model already features above-CMIP6-average skills (where CMIP6 denotes Coupled Model Intercomparison Project phase 6) in representing the climatology and competitive model throughput. Finally we identify remaining biases and suggest further improvements to be made to the model.
Gabriele Arduini, Sarah Keeley, Jonathan J. Day, Irina Sandu, Lorenzo Zampieri, and Gianpaolo Balsamo
American Geophysical Union (AGU)
Correctly representing the snow on sea‐ice has great potential to improve cryosphere‐atmosphere coupling in forecasting and monitoring (e.g., reanalysis) applications, via improved modeling of surface temperature, albedo and emissivity. This can also enhance the all‐weather all‐surface coupled data assimilation for atmospheric satellite radiances. Using wintertime observations from two Arctic field campaigns, SHEBA and N‐ICE2015, and satellite data, we explore the merits of different approaches to represent the snow over sea‐ice in a set of 5‐day coupled forecasts. Results show that representing the snow insulation effects is essential for capturing the wintertime surface temperature variability over sea‐ice and its response to changes in the atmospheric forcing. Modeling the snow over sea‐ice improves the representation of strong cooling events, reduces surface temperature biases in clear‐sky conditions and improves the simulation of surface‐based temperature inversions. In clear‐sky conditions, when using a multi‐layer snow scheme the root‐mean‐squared error in the surface temperature is reduced by about 60% for both N‐ICE2015 and SHEBA. This study also highlights the role of compensating errors in different components of the surface energy budget in the Arctic boundary layer. During warm air intrusions, errors in the surface temperature increase when cloud phase and cloud radiative processes are misrepresented in the model, inducing large errors in the net radiative energy at the surface. This work indicates that numerical weather prediction systems can fully benefit from a better representation of snow over sea‐ice, for example, with multi‐layer snow schemes, combined with improvements to other boundary layer processes including mixed phase clouds.
Jan Streffing, Tido Semmler, Lorenzo Zampieri, and Thomas Jung
American Meteorological Society
AbstractThe impact of Arctic sea ice decline on the weather and climate in mid-latitudes is still much debated, with observation suggesting a strong and models a much weaker link. In this study, we use the atmospheric model OpenIFS, in a set of model experiments following the protocol outlined in the Polar Amplification Model Intercomparison Project (PAMIP), to investigate whether the simulated atmospheric response to future changes in Arctic sea ice fundamentally depends on model resolution. More specifically, we increase the horizontal resolution of the model from 125km to 39km with 91 vertical levels; in a second step resolution is further increased to 16km with 137 levels in the vertical. The model does produce a response to sea ice decline with a weaker mid latitude Atlantic jet and increased blocking in the high latitude Atlantic, but no sensitivity to resolution can be detected with 100 members. Furthermore we find that the ensemble convergence toward the mean is not impacted by the model resolutions considered here.
Lorenzo Zampieri, Frank Kauker, Jörg Fröhle, Hiroshi Sumata, Elizabeth C. Hunke, and Helge F. Goessling
American Geophysical Union (AGU)
We have equipped the unstructured‐mesh global sea‐ice and ocean model FESOM2 with a set of physical parameterizations derived from the single‐column sea‐ice model Icepack. The update has substantially broadened the range of physical processes that can be represented by the model. The new features are directly implemented on the unstructured FESOM2 mesh, and thereby benefit from the flexibility that comes with it in terms of spatial resolution. A subset of the parameter space of three model configurations, with increasing complexity, has been calibrated with an iterative Green's function optimization method to test the impact of the model update on the sea‐ice representation. Furthermore, to explore the sensitivity of the results to different atmospheric forcings, each model configuration was calibrated separately for the NCEP‐CFSR/CFSv2 and ERA5 forcings. The results suggest that a complex model formulation leads to a better agreement between modeled and the observed sea‐ice concentration and snow thickness, while differences are smaller for sea‐ice thickness and drift speed. However, the choice of the atmospheric forcing also impacts the agreement of the FESOM2 simulations and observations, with NCEP‐CFSR/CFSv2 being particularly beneficial for the simulated sea‐ice concentration and ERA5 for sea‐ice drift speed. In this respect, our results indicate that parameter calibration can better compensate for differences among atmospheric forcings in a simpler model (i.e., sea‐ice has no heat capacity) than in more realistic formulations with a prognostic sea‐ice thickness distribution and sea ice enthalpy.
Longjiang Mu, Lars Nerger, Qi Tang, Svetlana N. Loza, Dmitry Sidorenko, Qiang Wang, Tido Semmler, Lorenzo Zampieri, Martin Losch, and Helge F. Goessling
American Geophysical Union (AGU)
This paper describes and evaluates the assimilation component of a seamless sea ice prediction system, which is developed based on the fully coupled Alfred Wegener Institute, Helmholtz Center for Polar and Marine Research Climate Model (AWI‐CM, v1.1). Its ocean/ice component with unstructured‐mesh discretization and smoothly varying spatial resolution enables seamless sea ice prediction across a wide range of space and time scales. The model is complemented with the Parallel Data Assimilation Framework to assimilate observations in the ocean/ice component with an Ensemble Kalman Filter. The focus here is on the data assimilation of the prediction system. First, the performance of the system is tested in a perfect‐model setting with synthetic observations. The system exhibits no drift for multivariate assimilation, which is a prerequisite for the robustness of the system. Second, real observational data for sea ice concentration, thickness, drift, and sea surface temperature are assimilated. The analysis results are evaluated against independent in situ observations and reanalysis data. Further experiments that assimilate different combinations of variables are conducted to understand their individual impacts on the model state. In particular, assimilating sea ice drift improves the sea ice thickness estimate, and assimilating sea surface temperature is able to avert a circulation bias of the free‐running model in the Arctic Ocean at middepth. Finally, we present preliminary results obtained with an extended system where the atmosphere is constrained by nudging toward reanalysis data, revealing challenges that still need to be overcome to adapt the ocean/ice assimilation. We consider this system a prototype on the way toward strongly coupled data assimilation across all model components.
Lorenzo Zampieri and Helge F. Goessling
American Geophysical Union (AGU)
To counteract global warming, a geoengineering approach that aims at intervening in the Arctic ice‐albedo feedback has been proposed. A large number of wind‐driven pumps shall spread seawater on the surface in winter to enhance ice growth, allowing more ice to survive the summer melt. We test this idea with a coupled climate model by modifying the surface exchange processes such that the physical effect of the pumps is simulated. Based on experiments with RCP 8.5 scenario forcing, we find that it is possible to keep the late‐summer sea ice cover at the current extent for the next ∼60 years. The increased ice extent is accompanied by significant Arctic late‐summer cooling by ∼1.3 K on average north of the polar circle (2021–2060). However, this cooling is not conveyed to lower latitudes. Moreover, the Arctic experiences substantial winter warming in regions with active pumps. The global annual‐mean near‐surface air temperature is reduced by only 0.02 K (2021–2060). Our results cast doubt on the potential of sea ice targeted geoengineering to mitigate climate change.
Lorenzo Zampieri, Helge F. Goessling, and Thomas Jung
American Geophysical Union (AGU)
Coupled subseasonal forecast systems with dynamical sea ice have the potential of providing important predictive information in polar regions. Here, we evaluate the ability of operational ensemble prediction systems to predict the location of the sea ice edge in Antarctica. Compared to the Arctic, Antarctica shows on average a 30% lower skill, with only one system remaining more skillful than a climatological benchmark up to ∼30 days ahead. Skill tends to be highest in the west Antarctic sector during the early freezing season. Most of the systems tend to overestimate the sea ice edge extent and fail to capture the onset of the melting season. All the forecast systems exhibit large initial errors. We conclude that subseasonal sea ice predictions could provide marginal support for decision‐making only in selected seasons and regions of the Southern Ocean. However, major progress is possible through investments in model development, forecast initialization and calibration.
Nils Hutter, Lorenzo Zampieri, and Martin Losch
Copernicus GmbH
Abstract. Leads and pressure ridges are dominant features of the Arctic sea ice cover. Not only do they affect heat loss and surface drag, but they also provide insight into the underlying physics of sea ice deformation. Due to their elongated shape they are referred to as linear kinematic features (LKFs). This paper introduces two methods that detect and track LKFs in sea ice deformation data and establish an LKF data set for the entire observing period of the RADARSAT Geophysical Processor System (RGPS). Both algorithms are available as open-source code and applicable to any gridded sea ice drift and deformation data. The LKF detection algorithm classifies pixels with higher deformation rates compared to the immediate environment as LKF pixels, divides the binary LKF map into small segments, and reconnects multiple segments into individual LKFs based on their distance and orientation relative to each other. The tracking algorithm uses sea ice drift information to estimate a first guess of LKF distribution and identifies tracked features by the degree of overlap between detected features and the first guess. An optimization of the parameters of both algorithms, as well as an extensive evaluation of both algorithms against handpicked features in a reference data set, is presented. A LKF data set is derived from RGPS deformation data for the years from 1996 to 2008 that enables a comprehensive description of LKFs. LKF densities and LKF intersection angles derived from this data set agree well with previous estimates. Further, a stretched exponential distribution of LKF length, an exponential tail in the distribution of LKF lifetimes, and a strong link to atmospheric drivers, here Arctic cyclones, are derived from the data set. Both algorithms are applied to output of a numerical sea ice model to compare the LKF intersection angles in a high-resolution Arctic sea ice simulation with the LKF data set.
Lorenzo Zampieri, Helge F. Goessling, and Thomas Jung
American Geophysical Union (AGU)
With retreating sea ice and increasing human activities in the Arctic come a growing need for reliable sea ice forecasts up to months ahead. We exploit the subseasonal‐to‐seasonal prediction database and provide the first thorough assessment of the skill of operational forecast systems in predicting the location of the Arctic sea ice edge on these time scales. We find large differences in skill between the systems, with some showing a lack of predictive skill even at short weather time scales and the best producing skillful forecasts more than 1.5 months ahead. This highlights that the area of subseasonal prediction in the Arctic is in an early stage but also that the prospects are bright, especially for late summer forecasts. To fully exploit this potential, it is argued that it will be imperative to reduce systematic model errors and develop advanced data assimilation capacity.