scholarly journals LIVVkit 2.1: automated and extensible ice sheet model validation

2019 ◽  
Vol 12 (3) ◽  
pp. 1067-1086 ◽  
Author(s):  
Katherine J. Evans ◽  
Joseph H. Kennedy ◽  
Dan Lu ◽  
Mary M. Forrester ◽  
Stephen Price ◽  
...  

Abstract. A collection of scientific analyses, metrics, and visualizations for robust validation of ice sheet models is presented using the Land Ice Verification and Validation toolkit (LIVVkit), version 2.1, and the LIVVkit Extensions repository (LEX), version 0.1. This software collection targets stand-alone ice sheet or coupled Earth system models, and handles datasets and analyses that require high-performance computing and storage. LIVVkit aims to enable efficient and fully reproducible workflows for postprocessing, analysis, and visualization of observational and model-derived datasets in a shareable format, whereby all data, methodologies, and output are distributed to users for evaluation. Extending from the initial LIVVkit software framework, we demonstrate Greenland ice sheet simulation validation metrics using the coupled Community Earth System Model (CESM) as well as an idealized stand-alone high-resolution Community Ice Sheet Model, version 2 (CISM2), coupled to the Albany/FELIX velocity solver (CISM-Albany or CISM-A). As one example of the capability, LIVVkit analyzes the degree to which models capture the surface mass balance (SMB) and identifies potential sources of bias, using recently available in situ and remotely sensed data as comparison. Related fields within atmosphere and land surface models, e.g., surface temperature, radiation, and cloud cover, are also diagnosed. Applied to the CESM1.0, LIVVkit identifies a positive SMB bias that is focused largely around Greenland's southwest region that is due to insufficient ablation.

2018 ◽  
Author(s):  
Katherine J. Evans ◽  
Joseph H. Kennedy ◽  
Dan Lu ◽  
Mary M. Forrester ◽  
Stephen Price ◽  
...  

Abstract. A collection of scientific analyses, metrics, and visualizations for robust validation of ice sheet models is presented using the LIVVkit package, version 2.1. This software collection targets stand-alone ice sheet or coupled Earth system models, and handles datasets and operations that require high-performance computing and storage. LIVVkit aims to enable efficient and fully reproducible workflows for post-processing, analysis, and visualization of observational and model-derived datasets in a shareable format, whereby all data, methodologies, and output are distributed to users for evaluation. We demonstrate LIVVkit validation for a Greenland ice sheet simulation using the coupled Community Earth System Model, CESM, as well as an idealized stand-alone high-resolution ice sheet model, CISM-Albany. As one example of the capability, LIVVkit analyzes the degree to which models capture the surface mass balance (SMB) and identifies potential sources of bias, using recently available in-situ and remotely sensed data as comparison. Related fields within atmosphere and land surface models, e.g. surface temperature, radiation, and cloud cover, are also diagnosed. Applied to the CESM1.0, LIVVkit identifies a positive SMB bias that is focused largely around Greenland's southwest region that is due to insufficient ablation.


2013 ◽  
Vol 26 (20) ◽  
pp. 7793-7812 ◽  
Author(s):  
Miren Vizcaíno ◽  
William H. Lipscomb ◽  
William J. Sacks ◽  
Jan H. van Angelen ◽  
Bert Wouters ◽  
...  

Abstract The modeling of the surface mass balance (SMB) of the Greenland Ice Sheet (GIS) requires high-resolution models in order to capture the observed large gradients in the steep marginal areas. Until now, global climate models have not been considered suitable to model ice sheet SMB owing to model biases and insufficient resolution. This study analyzes the GIS SMB simulated for the period 1850–2005 by the Community Earth System Model (CESM), which includes a new ice sheet component with multiple elevation classes for SMB calculations. The model is evaluated against observational data and output from the regional model Regional Atmospheric Climate Model version 2 (RACMO2). Because of a lack of major climate biases, a sophisticated calculation of snow processes (including surface albedo evolution) and an adequate downscaling technique, CESM is able to realistically simulate GIS surface climate and SMB. CESM SMB agrees reasonably well with in situ data from 475 locations (r = 0.80) and output from RACMO2 (r = 0.79). The simulated mean SMB for 1960–2005 is 359 ± 120 Gt yr−1 in the range of estimates from regional climate models. The simulated seasonal mass variability is comparable with mass observations from the Gravity Recovery and Climate Experiment (GRACE), with synchronous annual maximum (May) and minimum (August–September) and similar amplitudes of the seasonal cycle. CESM is able to simulate the bands of precipitation maxima along the southeast and northwest margins, but absolute precipitation rates are underestimated along the southeastern margin and overestimated in the high interior. The model correctly simulates the major ablation areas. Total refreezing represents 35% of the available liquid water (the sum of rain and melt).


2014 ◽  
Vol 27 (1) ◽  
pp. 215-226 ◽  
Author(s):  
Miren Vizcaíno ◽  
William H. Lipscomb ◽  
William J. Sacks ◽  
Michiel van den Broeke

Abstract This study presents the first twenty-first-century projections of surface mass balance (SMB) changes for the Greenland Ice Sheet (GIS) with the Community Earth System Model (CESM), which includes a new ice sheet component. For glaciated surfaces, CESM includes a sophisticated calculation of energy fluxes, surface albedo, and snowpack hydrology (melt, percolation, refreezing, etc.). To efficiently resolve the high SMB gradients at the ice sheet margins and provide surface forcing at the scale needed by ice sheet models, the SMB is calculated at multiple elevations and interpolated to a finer 5-km ice sheet grid. During a twenty-first-century simulation driven by representative concentration pathway 8.5 (RCP8.5) forcing, the SMB decreases from 372 ± 100 Gt yr−1 in 1980–99 to −78 ± 143 Gt yr−1 in 2080–99. The 2080–99 near-surface temperatures over the GIS increase by 4.7 K (annual mean) with respect to 1980–99, only 1.3 times the global increase (+3.7 K). Snowfall increases by 18%, while surface melt doubles. The ablation area increases from 9% of the GIS in 1980–99 to 28% in 2080–99. Over the ablation areas, summer downward longwave radiation and turbulent fluxes increase, while incoming shortwave radiation decreases owing to increased cloud cover. The reduction in GIS-averaged July albedo from 0.78 in 1980–99 to 0.75 in 2080–99 increases the absorbed solar radiation in this month by 12%. Summer warming is strongest in the north and east of Greenland owing to reduced sea ice cover. In the ablation area, summer temperature increases are smaller due to frequent periods of surface melt.


2018 ◽  
Author(s):  
Leonardus van Kampenhout ◽  
Alan M. Rhoades ◽  
Adam R. Herrington ◽  
Colin M. Zarzycki ◽  
Jan T. M. Lenaerts ◽  
...  

Abstract. In this study, the resolution dependence of the simulated Greenland Ice Sheet surface mass balance in the variable-resolution Community Earth System Model (VR-CESM) is investigated. Coupled atmosphere-land simulations are performed on three regionally refined grids over Greenland at 1° (~111 km), 0.5°(~55 km), and 0.25° (~28 km), maintaining a quasi-uniform resolution of 1° (~111 km) over the rest of the globe. The SMB in the accumulation zone is significantly improved compared to airborne radar and in-situ observations, with a general wetting at the margins and a drying in the interior GrIS. Total precipitation decreases with resolution, which is in line with best-available regional climate model results. In the ablation zone, VR-CESM starts developing a positive SMB bias in some locations. Potential driving mechanisms are proposed, amongst which are diversions in large scale circulation, changes in cloud cover, and changes in summer snowfall. Overall, our results demonstrate that VR-CESM is a viable new tool in the cryospheric sciences and can be used to dynamically downscale future scenarios and/or be interactively coupled to dynamical ice sheet models.


2021 ◽  
Author(s):  
Adam Schneider ◽  
Stephen Price ◽  
Jonathan Wolfe ◽  
Charles Zender

<p>Since 1993, nearly 10 percent of the observed rise in global mean sea level can be attributed to the coincident increase in surface mass loss from the Greenland Ice Sheet (GrIS) (Meredith et al., 2019; WCRP, 2018). To determine the GrIS surface mass balance (SMB), defined as the ice sheet’s annual net (surface) mass increase due to snow accumulation minus ablation, a climate model can be coupled to a snowpack model, which enables simulating relevant hydrologic processes including precipitation, phase changes, and runoff. Recent developments within the Energy Exascale Earth System Model (E3SM) include an active ice sheet component. To explore GrIS snowpack conditions relevant to present-day climate, we conduct simulations demonstrating the evolution of SMB and accumulation of snowpack depth, first in E3SM’s land component (ELM). After forcing ELM’s surface condition using 20th century atmospheric reanalysis, we couple ELM to E3SM’s atmosphere component (EAM) and simulate both atmospheric and snowpack conditions over a fixed GrIS geometry. Finally, we activate the MPAS-Albany Land Ice model (MALI), which enables prognostic SMB calculations including elevation-change feedbacks. We find broad agreement in the spatial patterns of GrIS SMB compared to regional climate model (RACMO) and Community Earth System Model (CESM) simulations. We provide insights regarding the use of a statistical downscaling method, which involves using multiple elevation classes with time-varying areal coverages within ELM grid-cells. Within this dynamic system, we can begin investigating elevation feedbacks, where the atmospheric temperature lapse rate allows the SMB to accelerate both positively and negatively in a rapidly changing climate.</p><p>References</p><ul><li>Meredith, M., M. Sommerkorn, S. Cassotta, C. Derksen, A. Ekaykin, A. Hollowed, G. Kofinas, A. Mackintosh, J. Melbourne-Thomas, M.M.C. Muelbert, G. Ottersen, H. Pritchard, and E.A.G. Schuur, 2019: Polar Regions. In: IPCC Special Report on the Ocean and Cryosphere in a Changing Climate [H.-O. Pörtner, D.C. Roberts, V. Masson-Delmotte, P. Zhai, M. Tignor, E. Poloczanska, K. Mintenbeck, A. Alegría, M. Nicolai, A. Okem, J. Petzold, B. Rama, N.M. Weyer (eds.)]. In press.</li> <li>WCRP Global Sea Level Budget Group: Global sea-level budget 1993–present, Earth Syst. Sci. Data, 10, 1551–1590, https://doi.org/10.5194/essd-10-1551-2018, 2018.</li> </ul>


2020 ◽  
Vol 14 (4) ◽  
pp. 1425-1435 ◽  
Author(s):  
Brice Noël ◽  
Leonardus van Kampenhout ◽  
Willem Jan van de Berg ◽  
Jan T. M. Lenaerts ◽  
Bert Wouters ◽  
...  

Abstract. We present a reconstruction of historical (1950–2014) surface mass balance (SMB) of the Greenland ice sheet (GrIS) using a high-resolution regional climate model (RACMO2; ∼11 km) to dynamically downscale the climate of the Community Earth System Model version 2 (CESM2; ∼111 km). After further statistical downscaling to 1 km spatial resolution, evaluation using in situ SMB measurements and remotely sensed GrIS mass change shows good agreement. Comparison with an ensemble of previously conducted RACMO2 simulations forced by climate reanalysis demonstrates that the current product realistically represents the long-term average and variability of individual SMB components and captures the recent increase in meltwater runoff that accelerated GrIS mass loss. This means that, for the first time, climate forcing from an Earth system model (CESM2), which assimilates no observations, can be used without additional corrections to reconstruct the historical GrIS SMB and its recent decline that initiated mass loss in the 1990s. This paves the way for attribution studies of future GrIS mass loss projections and contribution to sea level rise.


2019 ◽  
Vol 13 (12) ◽  
pp. 3193-3208 ◽  
Author(s):  
Raymond Sellevold ◽  
Leonardus van Kampenhout ◽  
Jan T. M. Lenaerts ◽  
Brice Noël ◽  
William H. Lipscomb ◽  
...  

Abstract. The modeling of ice sheets in Earth system models (ESMs) is an active area of research with applications to future sea level rise projections and paleoclimate studies. A major challenge for surface mass balance (SMB) modeling with ESMs arises from their coarse resolution. This paper evaluates the elevation class (EC) method as an SMB downscaling alternative to the dynamical downscaling of regional climate models. To this end, we compare EC-simulated elevation-dependent surface energy and mass balance gradients from the Community Earth System Model 1.0 (CESM1.0) with those from the regional climate model RACMO2.3. The EC implementation in CESM1.0 combines prognostic snow albedo, a multilayer snow model, and elevation corrections for two atmospheric forcing variables: temperature and humidity. Despite making no corrections for incoming radiation and precipitation, we find that the EC method in CESM1.0 yields similar SMB gradients to RACMO2.3, in part due to compensating biases in snowfall, surface melt, and refreezing gradients. We discuss the sensitivity of the results to the lapse rate used for the temperature correction. We also evaluate the impact of the EC method on the climate simulated by the ESM and find minor cooling over the Greenland ice sheet and Barents and Greenland seas, which compensates for a warm bias in the ESM due to topographic smoothing. Based on our diagnostic procedure to evaluate the EC method, we make several recommendations for future implementations.


2020 ◽  
Vol 14 (7) ◽  
pp. 2253-2265
Author(s):  
Jan T. M. Lenaerts ◽  
M. Drew Camron ◽  
Christopher R. Wyburn-Powell ◽  
Jennifer E. Kay

Abstract. The dominant mass input component of the Greenland Ice Sheet (GrIS) is precipitation, whose amounts and phase are poorly constrained by observations. Here we use spaceborne radar observations from CloudSat to map the precipitation frequency and phase on the GrIS, and we use those observations, in combination with a satellite simulator to enable direct comparison between observations and model, to evaluate present-day precipitation frequency in the Community Earth System Model (CESM). The observations show that substantial variability of snowfall frequency over the GrIS exists, that snowfall occurs throughout the year, and that snowfall frequency peaks in spring and fall. Rainfall is rare over the GrIS and only occurs in regions under 2000 m elevation and in the peak summer season. Although CESM overestimates the rainfall frequency, it reproduces the spatial and seasonal variability of precipitation frequency reasonably well. Driven by the high-emission, worst-case Representative Concentration Pathway (RCP) 8.5 scenario, CESM indicates that rainfall frequency will increase considerably across the GrIS, and will occur at higher elevations, potentially exposing a much larger GrIS area to rain and associated meltwater refreezing, firn warming, and reduced storage capacity. This technique can be applied to evaluate precipitation frequency in other climate models and can aid in planning future satellite campaigns.


2019 ◽  
Vol 13 (6) ◽  
pp. 1547-1564 ◽  
Author(s):  
Leonardus van Kampenhout ◽  
Alan M. Rhoades ◽  
Adam R. Herrington ◽  
Colin M. Zarzycki ◽  
Jan T. M. Lenaerts ◽  
...  

Abstract. In this study, the resolution dependence of the simulated Greenland ice sheet surface mass balance (GrIS SMB) in the variable-resolution Community Earth System Model (VR-CESM) is investigated. Coupled atmosphere–land simulations are performed on two regionally refined grids over Greenland at 0.5∘ (∼55 km) and 0.25∘ (∼28 km), maintaining a quasi-uniform resolution of 1∘ (∼111 km) over the rest of the globe. On the refined grids, the SMB in the accumulation zone is significantly improved compared to airborne radar and in situ observations, with a general wetting (more snowfall) at the margins and a drying (less snowfall) in the interior GrIS. Total GrIS precipitation decreases with resolution, which is in line with best-available regional climate model results. In the ablation zone, CESM starts developing a positive SMB bias with increased resolution in some basins, notably in the east and the north. The mismatch in ablation is linked to changes in cloud cover in VR-CESM, and a reduced effectiveness of the elevation classes subgrid parametrization in CESM. Overall, our pilot study introduces VR-CESM as a new tool in the cryospheric sciences, which could be used to dynamically downscale SMB in scenario simulations and to force dynamical ice sheet models through the CESM coupling framework.


2013 ◽  
Vol 26 (19) ◽  
pp. 7352-7371 ◽  
Author(s):  
William H. Lipscomb ◽  
Jeremy G. Fyke ◽  
Miren Vizcaíno ◽  
William J. Sacks ◽  
Jon Wolfe ◽  
...  

Abstract The Glimmer Community Ice Sheet Model (Glimmer-CISM) has been implemented in the Community Earth System Model (CESM). Glimmer-CISM is forced by a surface mass balance (SMB) computed in multiple elevation classes in the CESM land model and downscaled to the ice sheet grid. Ice sheet evolution is governed by the shallow-ice approximation with thermomechanical coupling and basal sliding. This paper describes and evaluates the initial model implementation for the Greenland Ice Sheet (GIS). The ice sheet model was spun up using the SMB from a coupled CESM simulation with preindustrial forcing. The model's sensitivity to three key ice sheet parameters was explored by running an ensemble of 100 GIS simulations to quasi equilibrium and ranking each simulation based on multiple diagnostics. With reasonable parameter choices, the steady-state GIS geometry is broadly consistent with observations. The simulated ice sheet is too thick and extensive, however, in some marginal regions where the SMB is anomalously positive. The top-ranking simulations were continued using surface forcing from CESM simulations of the twentieth century (1850–2005) and twenty-first century (2005–2100, with RCP8.5 forcing). In these simulations the GIS loses mass, with a resulting global-mean sea level rise of 16 mm during 1850–2005 and 60 mm during 2005–2100. This mass loss is caused mainly by increased ablation near the ice sheet margin, offset by reduced ice discharge to the ocean. Projected sea level rise is sensitive to the initial geometry, showing the importance of realistic geometry in the spun-up ice sheet.


Sign in / Sign up

Export Citation Format

Share Document