Challenges in simulating sea ice in Earth System Models

2012 ◽  
Vol 3 (6) ◽  
pp. 509-526 ◽  
Author(s):  
Dirk Notz
Keyword(s):  
Sea Ice ◽  
Earth System ◽  
System Models ◽  
Earth System Models

scholarly journals Geometric remapping of particle distributions in the Discrete Element Model for Sea Ice (DEMSI v0.0)

2021 ◽  
Author(s):  
Adrian K. Turner ◽  
Kara J. Peterson ◽  
Dan Bolintineanu
Keyword(s):  
Sea Ice ◽  
Element Model ◽  
Discrete Element ◽  
Discrete Element Model ◽  
Earth System ◽  
System Models ◽  
Ice Edge ◽  
Ice Floes ◽  
Lagrangian Motion ◽  
Earth System Models

Abstract. A new sea ice dynamical core, the Discrete Element Model for Sea Ice (DEMSI), is under development for use in coupled Earth system models. DEMSI is based on the discrete element method, which models collections of ice floes as interacting Lagrangian particles. In basin-scale sea ice simulations the Lagrangian motion results in significant convergence and ridging, which requires periodic remapping of sea ice variables from a deformed particle configuration back to an undeformed initial distribution. At the resolution required for Earth system models we cannot resolve individual sea ice floes, so we adopt the sub-gridscale thickness distribution used in continuum sea ice models. This choice leads to a series of hierarchical tracers depending on ice fractional area or concentration that must be remapped consistently. The circular discrete elements employed in DEMSI help improve the computational efficiency at the cost of increased complexity in the effective element area definitions for sea ice cover that are required for the accurate enforcement of conservation. An additional challenge is the accurate remapping of element values along the ice edge, the location of which varies due to the Lagrangian motion of the particles. In this paper we describe a particle-to-particle remapping approach based on well-established geometric remapping ideas that enforces conservation, bounds-preservation, and compatibility between associated tracer quantities, while also robustly managing remapping at the ice edge. One element of the remapping algorithm is a novel optimization-based flux correction that enforces concentration bounds in the case of non-uniform motion. We demonstrate the accuracy and utility of the algorithm in a series of numerical test cases.

scholarly journals Intercomparison and improvement of two-stream shortwave radiative transfer schemes in Earth system models for a unified treatment of cryospheric surfaces

2019 ◽  
Vol 13 (9) ◽  
pp. 2325-2343 ◽  
Author(s):  
Cheng Dang ◽  
Charles S. Zender ◽  
Mark G. Flanner
Keyword(s):  
Radiative Transfer ◽  
Sea Ice ◽  
Hybrid Model ◽  
Incident Beam ◽  
Radiative Properties ◽  
Transfer Model ◽  
Earth System ◽  
System Models ◽  
Near Ir ◽  
Earth System Models

Abstract. Snow is an important climate regulator because it greatly increases the surface albedo of middle and high latitudes of the Earth. Earth system models (ESMs) often adopt two-stream approximations with different radiative transfer techniques, the same snow therefore has different solar radiative properties depending whether it is on land or on sea ice. Here we intercompare three two-stream algorithms widely used in snow models, improve their predictions at large zenith angles, and introduce a hybrid model suitable for all cryospheric surfaces in ESMs. The algorithms are those employed by the SNow ICe and Aerosol Radiative (SNICAR) module used in land models, dEdd–AD used in Icepack, the column physics used in the Los Alamos sea ice model CICE and MPAS-Seaice, and a two-stream discrete-ordinate (2SD) model. Compared with a 16-stream benchmark model, the errors in snow visible albedo for a direct-incident beam from all three two-stream models are small (<±0.005) and increase as snow shallows, especially for aged snow. The errors in direct near-infrared (near-IR) albedo are small (<±0.005) for solar zenith angles θ<75∘, and increase as θ increases. For diffuse incidence under cloudy skies, dEdd–AD produces the most accurate snow albedo for both visible and near-IR (<±0.0002) with the lowest underestimate (−0.01) for melting thin snow. SNICAR performs similarly to dEdd–AD for visible albedos, with a slightly larger underestimate (−0.02), while it overestimates the near-IR albedo by an order of magnitude more (up to 0.04). 2SD overestimates both visible and near-IR albedo by up to 0.03. We develop a new parameterization that adjusts the underestimated direct near-IR albedo and overestimated direct near-IR heating persistent across all two-stream models for θ>75∘. These results are incorporated in a hybrid model SNICAR-AD, which can now serve as a unified solar radiative transfer model for snow in ESM land, land ice, and sea ice components.

scholarly journals A Variational Method for Sea Ice Ridging in Earth System Models

2019 ◽  
Vol 11 (3) ◽  
pp. 771-805 ◽  
Author(s):  
Andrew F. Roberts ◽  
Elizabeth C. Hunke ◽  
Samy M. Kamal ◽  
William H. Lipscomb ◽  
Christopher Horvat ◽  
...  
Keyword(s):  
Variational Method ◽  
Sea Ice ◽  
Earth System ◽  
System Models ◽  
Earth System Models

scholarly journals Mid-Holocene Antarctic sea-ice increase driven by marine ice sheet retreat

2021 ◽  
Vol 17 (1) ◽  
pp. 1-19
Author(s):  
Kate E. Ashley ◽  
Robert McKay ◽  
Johan Etourneau ◽  
Francisco J. Jimenez-Espejo ◽  
Alan Condron ◽  
...  
Keyword(s):  
Sea Ice ◽  
Climate Models ◽  
Global Climate Models ◽  
Earth System ◽  
Ice Shelf ◽  
Sea Ice Extent ◽  
Glacial Meltwater ◽  
System Models ◽  
Antarctic Sea Ice ◽  
Earth System Models

Abstract. Over recent decades Antarctic sea-ice extent has increased, alongside widespread ice shelf thinning and freshening of waters along the Antarctic margin. In contrast, Earth system models generally simulate a decrease in sea ice. Circulation of water masses beneath large-cavity ice shelves is not included in current Earth System models and may be a driver of this phenomena. We examine a Holocene sediment core off East Antarctica that records the Neoglacial transition, the last major baseline shift of Antarctic sea ice, and part of a late-Holocene global cooling trend. We provide a multi-proxy record of Holocene glacial meltwater input, sediment transport, and sea-ice variability. Our record, supported by high-resolution ocean modelling, shows that a rapid Antarctic sea-ice increase during the mid-Holocene (∼ 4.5 ka) occurred against a backdrop of increasing glacial meltwater input and gradual climate warming. We suggest that mid-Holocene ice shelf cavity expansion led to cooling of surface waters and sea-ice growth that slowed basal ice shelf melting. Incorporating this feedback mechanism into global climate models will be important for future projections of Antarctic changes.

Sea ice in Earth system models

Sea Ice ◽  
2016 ◽  
pp. 304-325 ◽  
Author(s):  
Dirk Notz ◽  
Cecilia M. Bitz
Keyword(s):  
Sea Ice ◽  
Earth System ◽  
System Models ◽  
Earth System Models

scholarly journals The evaluation of Earth System Models: discussion summary

2011 ◽  
Vol 6 ◽  
pp. 216-221
Author(s):  
Sönke Zaehle ◽  
Colin Prentice ◽  
Sarah Cornell
Keyword(s):  
Earth System ◽  
System Models ◽  
Earth System Models

scholarly journals Using field observations to inform thermal hydrology models of permafrost dynamics with ATS (v0.83)

2015 ◽  
Vol 8 (4) ◽  
pp. 3235-3292 ◽  
Author(s):  
A. L. Atchley ◽  
S. L. Painter ◽  
D. R. Harp ◽  
E. T. Coon ◽  
C. J. Wilson ◽  
...  
Keyword(s):  
Carbon Sink ◽  
Thermal Response ◽  
Field Measurements ◽  
Early Summer ◽  
Active Layer Thickness ◽  
Earth System ◽  
Model Refinement ◽  
System Models ◽  
Developed Surface ◽  
Earth System Models

Abstract. Climate change is profoundly transforming the carbon-rich Arctic tundra landscape, potentially moving it from a carbon sink to a carbon source by increasing the thickness of soil that thaws on a seasonal basis. However, the modeling capability and precise parameterizations of the physical characteristics needed to estimate projected active layer thickness (ALT) are limited in Earth System Models (ESMs). In particular, discrepancies in spatial scale between field measurements and Earth System Models challenge validation and parameterization of hydrothermal models. A recently developed surface/subsurface model for permafrost thermal hydrology, the Advanced Terrestrial Simulator (ATS), is used in combination with field measurements to calibrate and identify fine scale controls of ALT in ice wedge polygon tundra in Barrow, Alaska. An iterative model refinement procedure that cycles between borehole temperature and snow cover measurements and simulations functions to evaluate and parameterize different model processes necessary to simulate freeze/thaw processes and ALT formation. After model refinement and calibration, reasonable matches between simulated and measured soil temperatures are obtained, with the largest errors occurring during early summer above ice wedges (e.g. troughs). The results suggest that properly constructed and calibrated one-dimensional thermal hydrology models have the potential to provide reasonable representation of the subsurface thermal response and can be used to infer model input parameters and process representations. The models for soil thermal conductivity and snow distribution were found to be the most sensitive process representations. However, information on lateral flow and snowpack evolution might be needed to constrain model representations of surface hydrology and snow depth.

scholarly journals Evolution of the Arabian Sea Upwelling from the Last Millennium to the Future as Simulated by Earth System Models

Climate ◽  
2021 ◽  
Vol 9 (5) ◽  
pp. 72
Author(s):  
Xing Yi ◽  
Birgit Hünicke ◽  
Eduardo Zorita
Keyword(s):  
Arabian Sea ◽  
Positive Trend ◽  
Last Millennium ◽  
Earth System ◽  
System Models ◽  
Sediment Records ◽  
The Future ◽  
Water Stratification ◽  
Arabian Sea Upwelling ◽  
Earth System Models

Arabian Sea upwelling in the past has been generally studied based on the sediment records. We apply two earth system models and analyze the simulated water vertical velocity to investigate coastal upwelling in the western Arabian Sea over the last millennium. In addition, two models with slightly different configurations are also employed to study the upwelling in the 21st century under the strongest and the weakest greenhouse gas emission scenarios. With a negative long-term trend caused by the orbital forcing of the models, the upwelling over the last millennium is found to be closely correlated with the sea surface temperature, the Indian summer Monsoon and the sediment records. The future upwelling under the Representative Concentration Pathway (RCP) 8.5 scenario reveals a negative trend, in contrast with the positive trend displayed by the upwelling favorable along-shore winds. Therefore, it is likely that other factors, like water stratification in the upper ocean layers caused by the stronger surface warming, overrides the effect from the upwelling favorable wind. No significant trend is found for the upwelling under the RCP2.6 scenario, which is likely due to a compensation between the opposing effects of the increase in upwelling favorable winds and the water stratification.

GFDL’s ESM2 Global Coupled Climate–Carbon Earth System Models. Part I: Physical Formulation and Baseline Simulation Characteristics

2012 ◽  
Vol 25 (19) ◽  
pp. 6646-6665 ◽  
Author(s):  
John P. Dunne ◽  
Jasmin G. John ◽  
Alistair J. Adcroft ◽  
Stephen M. Griffies ◽  
Robert W. Hallberg ◽  
...  
Keyword(s):  
Climate Changes ◽  
Southern Oscillation ◽  
Heat Content ◽  
Ocean Dynamics ◽  
Thermocline Depth ◽  
Earth System ◽  
System Models ◽  
Model Version ◽  
Earth System Models

Abstract The physical climate formulation and simulation characteristics of two new global coupled carbon–climate Earth System Models, ESM2M and ESM2G, are described. These models demonstrate similar climate fidelity as the Geophysical Fluid Dynamics Laboratory’s previous Climate Model version 2.1 (CM2.1) while incorporating explicit and consistent carbon dynamics. The two models differ exclusively in the physical ocean component; ESM2M uses Modular Ocean Model version 4p1 with vertical pressure layers while ESM2G uses Generalized Ocean Layer Dynamics with a bulk mixed layer and interior isopycnal layers. Differences in the ocean mean state include the thermocline depth being relatively deep in ESM2M and relatively shallow in ESM2G compared to observations. The crucial role of ocean dynamics on climate variability is highlighted in El Niño–Southern Oscillation being overly strong in ESM2M and overly weak in ESM2G relative to observations. Thus, while ESM2G might better represent climate changes relating to total heat content variability given its lack of long-term drift, gyre circulation, and ventilation in the North Pacific, tropical Atlantic, and Indian Oceans, and depth structure in the overturning and abyssal flows, ESM2M might better represent climate changes relating to surface circulation given its superior surface temperature, salinity, and height patterns, tropical Pacific circulation and variability, and Southern Ocean dynamics. The overall assessment is that neither model is fundamentally superior to the other, and that both models achieve sufficient fidelity to allow meaningful climate and earth system modeling applications. This affords the ability to assess the role of ocean configuration on earth system interactions in the context of two state-of-the-art coupled carbon–climate models.

Sign in / Sign up

Export Citation Format

Share Document