Impact of melt ponds on Arctic sea ice simulations from 1990 to 2007

The extent and thickness of the Arctic sea ice cover has decreased dramatically in the past few decades with minima in sea ice extent in September 2007 and 2011 and climate models did not predict this decline. One of the processes poorly represented in sea ice models is the formation and evolution of melt ponds. Melt ponds form on Arctic sea ice during the melting season and their presence affects the heat and mass balances of the ice cover, mainly by decreasing the value of the surface albedo by up to 20%. We have developed a melt pond model suitable for forecasting the presence of melt ponds based on sea ice conditions. This model has been incorporated into the Los Alamos CICE sea ice model, the sea ice component of several IPCC climate models. Simulations for the period 1990 to 2007 are in good agreement with observed ice concentration. In comparison to simulations without ponds, the September ice volume is nearly 40% lower. Sensitivity studies within the range of uncertainty reveal that, of the parameters pertinent to the present melt pond parameterization and for our prescribed atmospheric and oceanic forcing, variations of optical properties and the amount of snowfall have the strongest impact on sea ice extent and volume. We conclude that melt ponds will play an increasingly important role in the melting of the Arctic ice cover and their incorporation in the sea ice component of Global Circulation Models is essential for accurate future sea ice forecasts

The Community Climate System Model version 4

Author Posting. © American Meteorological Society, 2011. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 24 (2011): 4973–4991, doi:10.1175/2011JCLI4083.1.The fourth version of the Community Climate System Model (CCSM4) was recently completed and released to the climate community. This paper describes developments to all CCSM components, and documents fully coupled preindustrial control runs compared to the previous version, CCSM3. Using the standard atmosphere and land resolution of 1° results in the sea surface temperature biases in the major upwelling regions being comparable to the 1.4°-resolution CCSM3. Two changes to the deep convection scheme in the atmosphere component result in CCSM4 producing El Niño–Southern Oscillation variability with a much more realistic frequency distribution than in CCSM3, although the amplitude is too large compared to observations. These changes also improve the Madden–Julian oscillation and the frequency distribution of tropical precipitation. A new overflow parameterization in the ocean component leads to an improved simulation of the Gulf Stream path and the North Atlantic Ocean meridional overturning circulation. Changes to the CCSM4 land component lead to a much improved annual cycle of water storage, especially in the tropics. The CCSM4 sea ice component uses much more realistic albedos than CCSM3, and for several reasons the Arctic sea ice concentration is improved in CCSM4. An ensemble of twentieth-century simulations produces a good match to the observed September Arctic sea ice extent from 1979 to 2005. The CCSM4 ensemble mean increase in globally averaged surface temperature between 1850 and 2005 is larger than the observed increase by about 0.4°C. This is consistent with the fact that CCSM4 does not include a representation of the indirect effects of aerosols, although other factors may come into play. The CCSM4 still has significant biases, such as the mean precipitation distribution in the tropical Pacific Ocean, too much low cloud in the Arctic, and the latitudinal distributions of shortwave and longwave cloud forcings.National Science Foundation, which sponsors NCAR and the CCSM Project. The project is also sponsored by the U.S. Department of Energy (DOE). Thanks are also due to the many other software engineers and scientists who worked on developing CCSM4, and to the Computational and Information Systems Laboratory at NCAR, which provided the computing resources through the Climate Simulation Laboratory. Hunke was supported within theClimate, Ocean and Sea Ice Modeling project at Los Alamos National Laboratory, which is funded by the Biological and Environmental Research division of the DOE Office of Science. The Los Alamos National Laboratory is operated by theDOENationalNuclear Security Administration under Contract DE-AC52-06NA25396. Raschwas supported by theDOEOffice of Science, Earth System Modeling Program, which is part of the DOE Climate Change Research Program. The Pacific Northwest National Laboratory is operated forDOEbyBattelle Memorial Institute under Contract DE-AC06-76RLO 1830. Worley was supported by the Climate Change Research Division of the Office of Biological and Environmental Research and by the Office ofAdvanced Scientific Computing Research, both in the DOE Office of Science, under Contract DE-AC05-00OR22725 with UT-Batelle, LLC

Toward quantifying the increasing role oceanic heat in sea ice loss in the new Arctic

Author Posting. © American Meteorological Society, 2015. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Bulletin of the American Meteorological Society 96 (2015): 2079–2105, doi:10.1175/BAMS-D-13-00177.1.The loss of Arctic sea ice has emerged as a leading signal of global warming. This, together with acknowledged impacts on other components of the Earth system, has led to the term “the new Arctic.” Global coupled climate models predict that ice loss will continue through the twenty-first century, with implications for governance, economics, security, and global weather. A wide range in model projections reflects the complex, highly coupled interactions between the polar atmosphere, ocean, and cryosphere, including teleconnections to lower latitudes. This paper summarizes our present understanding of how heat reaches the ice base from the original sources—inflows of Atlantic and Pacific Water, river discharge, and summer sensible heat and shortwave radiative fluxes at the ocean/ice surface—and speculates on how such processes may change in the new Arctic. The complexity of the coupled Arctic system, and the logistic and technological challenges of working in the Arctic Ocean, require a coordinated interdisciplinary and international program that will not only improve understanding of this critical component of global climate but will also provide opportunities to develop human resources with the skills required to tackle related problems in complex climate systems. We propose a research strategy with components that include 1) improved mapping of the upper- and middepth Arctic Ocean, 2) enhanced quantification of important process, 3) expanded long-term monitoring at key heat-flux locations, and 4) development of numerical capabilities that focus on parameterization of heat-flux mechanisms and their interactions.2016-06-0

The DOE E3SM Coupled Model Version 1: Overview and Evaluation at Standard Resolution

This work documents the first version of the U.S. Department of Energy (DOE) new Energy Exascale Earth System Model (E3SMv1). We focus on the standard resolution of the fully coupled physical model designed to address DOE mission-relevant water cycle questions. Its components include atmosphere and land (110-km grid spacing), ocean and sea ice (60 km in the midlatitudes and 30 km at the equator and poles), and river transport (55 km) models. This base configuration will also serve as a foundation for additional configurations exploring higher horizontal resolution as well as augmented capabilities in the form of biogeochemistry and cryosphere configurations. The performance of E3SMv1 is evaluated by means of a standard set of Coupled Model Intercomparison Project Phase 6 (CMIP6) Diagnosis, Evaluation, and Characterization of Klima simulations consisting of a long preindustrial control, historical simulations (ensembles of fully coupled and prescribed SSTs) as well as idealized CO2 forcing simulations. The model performs well overall with biases typical of other CMIP-class models, although the simulated Atlantic Meridional Overturning Circulation is weaker than many CMIP-class models. While the E3SMv1 historical ensemble captures the bulk of the observed warming between preindustrial (1850) and present day, the trajectory of the warming diverges from observations in the second half of the twentieth century with a period of delayed warming followed by an excessive warming trend. Using a two-layer energy balance model, we attribute this divergence to the model’s strong aerosol-related effective radiative forcing (ERFari+aci = -1.65 W/m2) and high equilibrium climate sensitivity (ECS = 5.3 K).Plain Language SummaryThe U.S. Department of Energy funded the development of a new state-of-the-art Earth system model for research and applications relevant to its mission. The Energy Exascale Earth System Model version 1 (E3SMv1) consists of five interacting components for the global atmosphere, land surface, ocean, sea ice, and rivers. Three of these components (ocean, sea ice, and river) are new and have not been coupled into an Earth system model previously. The atmosphere and land surface components were created by extending existing components part of the Community Earth System Model, Version 1. E3SMv1’s capabilities are demonstrated by performing a set of standardized simulation experiments described by the Coupled Model Intercomparison Project Phase 6 (CMIP6) Diagnosis, Evaluation, and Characterization of Klima protocol at standard horizontal spatial resolution of approximately 1° latitude and longitude. The model reproduces global and regional climate features well compared to observations. Simulated warming between 1850 and 2015 matches observations, but the model is too cold by about 0.5 °C between 1960 and 1990 and later warms at a rate greater than observed. A thermodynamic analysis of the model’s response to greenhouse gas and aerosol radiative affects may explain the reasons for the discrepancy.Key PointsThis work documents E3SMv1, the first version of the U.S. DOE Energy Exascale Earth System ModelThe performance of E3SMv1 is documented with a set of standard CMIP6 DECK and historical simulations comprising nearly 3,000 yearsE3SMv1 has a high equilibrium climate sensitivity (5.3 K) and strong aerosol-related effective radiative forcing (-1.65 W/m2)Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/151288/1/jame20860_am.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/151288/2/jame20860.pd

Comparison of elastic-viscous-plastic and viscous-plastic dynamics models using a high resolution Arctic sea ice model

A nonlinear viscous-plastic (VP) rheology proposed by Hibler (1979) has been demonstrated to be the most suitable of the rheologies commonly used for modeling sea ice dynamics (Kreyscher et al., 1997). However, the presence of a huge range of effective viscosities hinders numerical implementations of this model, particularly on high resolution grids or when the ice model is coupled to an ocean or atmosphere model. Hunke and Dukowicz (1997) have modified the VP model by including elastic waves as a numerical regularization in the case of zero strain rate. This modification (EVP) allows an efficient, fully explicit discretization that adapts well to parallel architectures. We present a comparison of EVP and VP dynamics model results from two 5-year simulations of Arctic sea ice, obtained with a high resolution sea ice model. The purpose of the comparison is to determine how differently the two dynamics models behave, and to decide whether the elastic-viscous-plastic model is preferable for high resolution climate simulations, considering its high efficiency in parallel computation. Results from the first year of this experiment (1990) are discussed in detail in Hunke and Zhang (1997). The model equations were solved on a 300x360 mesh with a resolution of about 18 km and a 4 hr timestep, and driven by 1990-1994 ECMWF atmospheric data above and 1992 ocean model output below. The two runs are identical except for the ice rheology; the EVP and VP dynamics models used here are documented in Hunke and Dukowicz (1997) and Zhang et al. (1997), respectively. For further details about the model specifications and experiment design, see Hunke and Zhang (1997).DOE CHAMMP Progra

Ridging, strength, and stability in high‐resolution sea ice models

The article of record as published may be found at http://dx.doi.org/10.1029/2005JC003355In multicategory sea ice models the compressive strength of the ice pack is often assumed to be a function of the potential energy of pressure ridges. This assumption, combined with other standard features of ridging schemes, allows the ice strength to change dramatically on short timescales. In high-resolution (~10 km) sea ice models with a typical time step (~1 hour), abrupt strength changes can lead to large internal stress gradients that destabilize the flow. The unstable flow is characterized by large oscillations in ice concentration, thickness, strength, velocity, and strain rates. Straightforward, physically motivated changes in the ridging scheme can reduce the likelihood of abrupt strength changes and improve stability. In simple test problems with flow toward and around topography, stability is significantly enhanced by eliminating the threshold fraction G* in the ridging participation function. Use of an exponential participation function increases the maximum stable time step at 10-km resolution from less than 30 min to about 2 hours. Modifying the redistribution function to build thinner ridges modestly improves stability and also gives better agreement between modeled and observed thickness distributions. Allowing the ice strength to increase linearly with the mean ice thickness improves stability but probably underestimates the maximum stresses.Climate Change Prediction Program (CCPP) and the Scientific Discovery through Advanced Computing (SciDAC) program of the U.S. Department of Energy’s Office of ScienceShelf-Basin Interaction (SBI) program of the U.S. National Science FoundationCCPP and through the International Arctic Research Center’s Arctic Ocean Model Intercomparison Project (AOMIP)Arctic Region Supercomputing Center (ARSC) through the U.S. Department of Defense High Performance Computer Modernization Program (HPCMP

Level-ice melt ponds in the Los Alamos sea ice model, CICE

A new meltpond parameterization has been developed for the CICE sea ice model, taking advantage of the level ice tracer available in the model. The ponds evolve according to physically based process descriptions, assuming a depth-area ratio for changes in pond volume. A novel aspect of the new scheme is that the ponds are carried as tracers on the level ice area of each thickness category, thus limiting their spatial extent based on the simulated sea ice topography. This limiting is meant to approximate the horizontal drainage of melt water into depressions in ice floes. Simulated melt pond processes include collection of liquid melt water and rain into ponds, drainage through permeable sea ice or over the edges of floes, infiltration of snow by pond water, and refreezing of ponds. Furthermore, snow that falls on top of ponds whose top surface has refrozen blocks radiation from penetrating into the ponds and sea ice below. Along with a control simulation, we present a range of sensitivity tests to parameters related to each subprocess described by the parameterization. With the exception of one parameter that alters the albedo of snow-covered pond ice, results are not highly sensitive to these parameters unless an entire process is removed. The snow simulation itself is critical, because the volume of snow deposition and rate of snow melt largely determine the timing and extent of the simulated melt ponds. Nevertheless, compensating effects moderate the model's sensitivity to precipitation changes. For instance, infiltration of the snow by melt water postpones the appearance of ponds and the subsequent acceleration of melting through albedo feedback, while snow on top of refrozen pond ice also reduces the ponds' effect on the radiation budget. By construction, the model simulation of level and ridged ice is also important for this parameterization. We find that as sea ice thins, either through time or when comparing sensitivity tests, the area of level ice increases. This leads to an enhanced thinning feedback in the model, because a greater ice area may be exposed to ponding and further thinning due to lowered albedo. © 2012 Elsevier Ltd

Impact of Sea‐Ice Model Complexity on the Performance of an Unstructured‐Mesh Sea‐Ice/Ocean Model under Different Atmospheric Forcings

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.Plain Language Summary: The role of model complexity in determining the performance of sea‐ice numerical simulations is still not completely understood. Some studies suggest that a more sophisticated description of the sea‐ice physics leads to simulations that agree better with sea‐ice observations. Others, however, fail to establish a link between complex model formulations and improved model performance. Here, we investigate this open question by analyzing a set of sea‐ice simulations performed with a revised and improved sea‐ice model that features substantial modularity in terms of model complexity. Ten model parameters in three different model configurations are optimized to improve the agreement between model results and observations, allowing a fair comparison between model configurations with varying complexity. The model optimization is repeated for two different atmospheric forcings to shed light on the relationship between model complexity and other sources of uncertainty in the sea‐ice simulations, such as those associated with the atmospheric conditions. The results suggest that a more complex formulation of our model can lead to a more appropriate representation of sea ice concentration and snow thickness, while it is less relevant for sea‐ice thickness and drift.Key Points: Increased sea‐ice model complexity can improve the simulated sea‐ice concentration and snow thickness Sea‐ice thickness and drift are only weakly affected by model complexity Parameter calibration can better compensate for differences between atmospheric forcings in a simpler modelBundesministerium für Bildung und Forschung (BMBF) http://dx.doi.org/10.13039/501100002347European Commission (EC) http://dx.doi.org/10.13039/501100000780US Department of Energy (DOE