scholarly journals Mapping technique of climate fields between GCM's and ice models

2010 ◽  
Vol 3 (1) ◽  
pp. 13-41 ◽  
Author(s):  
T. J. Reerink ◽  
M. A. Kliphuis ◽  
R. S. W. van de Wal

Abstract. Here, we present a mapping method OBLIMAP, which projects and interpolates fields like surface temperature, surface mass balance, and surface height between a geographical based coordinate system of a General Circulation Model (GCM) and a rectangular based Ice Model (IM). We derive an oblique stereographic projection and its inverse, which holds for any area at the Earth's surface, and which can be combined with two different interpolation methods. The first one is suited to interpolate the projected fields of a coarse GCM grid on a fine meshed IM grid. The second one is appropriate for the opposite case. Both grids are allowed to be arbitrary and irregularly spaced. Therefore the OBLIMAP technique is suitable for any GCM-IM combination. After a first scan of the GCM grid coordinates and the specification of the IM grid, fast mapping of various fields is possible. To and fro (GCM-IM-GCM) mapping tests with the Climate Community System Model (CCSM) at T42 resolution (~313 km) and the Regional Atmospheric Climate Model (RACMO) at ~11 km and ~55 km, show average temperature differences of less than 0.1 K with small standard deviations. OBLIMAP, available at GMD, is an accurate, robust and well-documented mapping method for coupling an IM with a GCM or to map state of the art initial and forcing fields available at geographical coordinates to any local IM grid with an optimal centered oblique projection. Currently, the oblique stereographic and the oblique Lambert azimuthal equal-area projections for both the sphere and the ellipsoid are implemented in OBLIMAP.

2009 ◽  
Vol 2 (2) ◽  
pp. 935-999
Author(s):  
T. J. Reerink ◽  
M. A. Kliphuis ◽  
R. S. W. van de Wal

Abstract. Here, we present a mapping method OBLIMAP, which projects and interpolates fields like surface temperature, surface mass balance, and surface height between a geographical based coordinate system of a Global Circulation Model (GCM) and a rectangular based Ice Model (IM). We derive an oblique stereographic projection and its inverse, which holds for any area at the Earth surface, and can be combined with two different interpolation methods. The first one is suited to interpolate the projected fields of a coarse GCM grid on a fine meshed IM grid. The second one is appropriate for the opposite case. Both grids are allowed to be arbitrary and irregularly spaced. Therefore the OBLIMAP technique is suitable for any GCM-IM combination. After a first scan of the GCM grid coordinates and the specification of the IM grid, fast mapping of various fields is possible. To and fro (GCM-IM-GCM) mapping tests with the Climate Community System Model (CCSM) at T42 resolution (~313 km) and the Regional Atmospheric Climate Model (RACMO) at ~11 km, show average temperature differences of less than 0.1 K with small standard deviations. OBLIMAP, available at GMD, is an accurate, robust and well-documented mapping method for coupling an IM with a GCM or to map state of the art initial and forcing fields available at geographical coordinates to any local IM grid with an optimal centered projection.


2018 ◽  
Vol 12 (4) ◽  
pp. 1499-1510 ◽  
Author(s):  
Marcus Lofverstrom ◽  
Johan Liakka

Abstract. Coupled climate–ice sheet simulations have been growing in popularity in recent years. Experiments of this type are however challenging as ice sheets evolve over multi-millennial timescales, which is beyond the practical integration limit of most Earth system models. A common method to increase model throughput is to trade resolution for computational efficiency (compromise accuracy for speed). Here we analyze how the resolution of an atmospheric general circulation model (AGCM) influences the simulation quality in a stand-alone ice sheet model. Four identical AGCM simulations of the Last Glacial Maximum (LGM) were run at different horizontal resolutions: T85 (1.4∘), T42 (2.8∘), T31 (3.8∘), and T21 (5.6∘). These simulations were subsequently used as forcing of an ice sheet model. While the T85 climate forcing reproduces the LGM ice sheets to a high accuracy, the intermediate resolution cases (T42 and T31) fail to build the Eurasian ice sheet. The T21 case fails in both Eurasia and North America. Sensitivity experiments using different surface mass balance parameterizations improve the simulations of the Eurasian ice sheet in the T42 case, but the compromise is a substantial ice buildup in Siberia. The T31 and T21 cases do not improve in the same way in Eurasia, though the latter simulates the continent-wide Laurentide ice sheet in North America. The difficulty to reproduce the LGM ice sheets in the T21 case is in broad agreement with previous studies using low-resolution atmospheric models, and is caused by a substantial deterioration of the model climate between the T31 and T21 resolutions. It is speculated that this deficiency may demonstrate a fundamental problem with using low-resolution atmospheric models in these types of experiments.


2011 ◽  
Vol 24 (3) ◽  
pp. 609-623 ◽  
Author(s):  
Sebastian H. Mernild ◽  
Glen E. Liston ◽  
Christopher A. Hiemstra ◽  
Jens H. Christensen ◽  
Martin Stendel ◽  
...  

Abstract A regional atmospheric model, the HIRHAM4 regional climate model (RCM) using boundary conditions from the ECHAM5 atmosphere–ocean general circulation model (AOGCM), was downscaled to a 500-m gridcell increment using SnowModel to simulate 131 yr (1950–2080) of hydrologic cycle evolution in west Greenland’s Kangerlussuaq drainage. Projected changes in the Greenland Ice Sheet (GrIS) surface mass balance (SMB) and runoff are relevant for potential hydropower production and prediction of ecosystem changes in sensitive Kangerlussuaq Fjord systems. Mean annual surface air temperatures and precipitation in the Kangerlussuaq area were simulated to increase by 3.4°C and 95 mm water equivalent (w.eq.), respectively, between 1950 and 2080. The local Kangerlussuaq warming was less than the average warming of 4.8°C simulated for the entire GrIS. The Kangerlussuaq SMB loss increased by an average of 0.3 km3 because of a 0.4 km3 rise in precipitation, 0.1 km3 rise in evaporation and sublimation, and 0.6 km3 gain in runoff (1950–2080). By 2080, the spring runoff season begins approximately three weeks earlier. The average modeled SMB and runoff is approximately −0.1 and 1.2 km3 yr−1, respectively, indicating that ∼10% of the Kangerlussuaq runoff is explained by the GrIS SMB net loss. The cumulative net volume loss (1950–2080) from SMB was 15.9 km3, and runoff was 151.2 km3 w.eq. This runoff volume is expected to have important hydrodynamic and ecological impacts on the stratified salinity in the Kangerlussuaq Fjord and on the transport of freshwater to the ocean.


2006 ◽  
Vol 19 (16) ◽  
pp. 3903-3931 ◽  
Author(s):  
H. Schmidt ◽  
G. P. Brasseur ◽  
M. Charron ◽  
E. Manzini ◽  
M. A. Giorgetta ◽  
...  

Abstract This paper introduces the three-dimensional Hamburg Model of the Neutral and Ionized Atmosphere (HAMMONIA), which treats atmospheric dynamics, radiation, and chemistry interactively for the height range from the earth’s surface to the thermosphere (approximately 250 km). It is based on the latest version of the ECHAM atmospheric general circulation model of the Max Planck Institute for Meteorology in Hamburg, Germany, which is extended to include important radiative and dynamical processes of the upper atmosphere and is coupled to a chemistry module containing 48 compounds. The model is applied to study the effects of natural and anthropogenic climate forcing on the atmosphere, represented, on the one hand, by the 11-yr solar cycle and, on the other hand, by a doubling of the present-day concentration of carbon dioxide. The numerical experiments are analyzed with the focus on the effects on temperature and chemical composition in the mesopause region. Results include a temperature response to the solar cycle by 2 to 10 K in the mesopause region with the largest values occurring slightly above the summer mesopause. Ozone in the secondary maximum increases by up to 20% for solar maximum conditions. Changes in winds are in general small. In the case of a doubling of carbon dioxide the simulation indicates a cooling of the atmosphere everywhere above the tropopause but by the smallest values around the mesopause. It is shown that the temperature response up to the mesopause is strongly influenced by changes in dynamics. During Northern Hemisphere summer, dynamical processes alone would lead to an almost global warming of up to 3 K in the uppermost mesosphere.


2015 ◽  
Vol 15 (10) ◽  
pp. 5537-5555 ◽  
Author(s):  
R. Eichinger ◽  
P. Jöckel ◽  
S. Brinkop ◽  
M. Werner ◽  
S. Lossow

Abstract. This modelling study aims at an improved understanding of the processes that determine the water vapour budget in the stratosphere by means of the investigation of water isotope ratios. An additional (and separate from the actual) hydrological cycle has been introduced into the chemistry–climate model EMAC, including the water isotopologues HDO and H218O and their physical fractionation processes. Additionally an explicit computation of the contribution of methane oxidation to H2O and HDO has been incorporated. The model expansions allow detailed analyses of water vapour and its isotope ratio with respect to deuterium throughout the stratosphere and in the transition region to the troposphere. In order to assure the correct representation of the water isotopologues in the model's hydrological cycle, the expanded system has been evaluated in several steps. The physical fractionation effects have been evaluated by comparison of the simulated isotopic composition of precipitation with measurements from a ground-based network (GNIP) and with the results from the isotopologue-enabled general circulation model ECHAM5-wiso. The model's representation of the chemical HDO precursor CH3D in the stratosphere has been confirmed by a comparison with chemical transport models (1-D, CHEM2D) and measurements from radiosonde flights. Finally, the simulated stratospheric HDO and the isotopic composition of water vapour have been evaluated, with respect to retrievals from three different satellite instruments (MIPAS, ACE-FTS, SMR). Discrepancies in stratospheric water vapour isotope ratios between two of the three satellite retrievals can now partly be explained.


2014 ◽  
Vol 7 (3) ◽  
pp. 883-907 ◽  
Author(s):  
R. Fischer ◽  
S. Nowicki ◽  
M. Kelley ◽  
G. A. Schmidt

Abstract. The method of elevation classes, in which the ice surface model is run at multiple elevations within each grid cell, has proven to be a useful way for a low-resolution atmosphere inside a general circulation model (GCM) to produce high-resolution downscaled surface mass balance fields for use in one-way studies coupling atmospheres and ice flow models. Past uses of elevation classes have failed to conserve mass and energy because the transformation used to regrid to the atmosphere was inconsistent with the transformation used to downscale to the ice model. This would cause problems for two-way coupling. A strategy that resolves this conservation issue has been designed and is presented here. The approach identifies three grids between which data must be regridded and five transformations between those grids required by a typical coupled atmosphere–ice flow model. This paper develops a theoretical framework for the problem and shows how each of these transformations may be achieved in a consistent, conservative manner. These transformations are implemented in Glint2, a library used to couple atmosphere models with ice models. Source code and documentation are available for download. Confounding real-world issues are discussed, including the use of projections for ice modeling, how to handle dynamically changing ice geometry, and modifications required for finite element ice models.


2009 ◽  
Vol 137 (9) ◽  
pp. 2851-2868 ◽  
Author(s):  
Masaru Inatsu ◽  
Masahide Kimoto

Abstract This study newly developed the interactively nested climate model (INCL) using a general circulation model (GCM) interactively nested with a regional atmospheric model (RAM). One interactive experiment with finer RAM topography and another with coarser topography, as well as offline versions of each experiment, were performed to investigate the effects of subsynoptic-scale eddies and subsynoptic-scale mountains in northeast Asia on the larger-scale climate, using the GCM with T42 atmosphere and the RAM with 40-km mesh size in the INCL system. The subsynoptic-scale eddy effect restrictively increased synoptic-scale eddy activity within the RAM domain. In contrast, subsynoptic-scale mountains had the effect of robust anticyclonic circulation around the Sea of Japan and effectively forced larger-scale circulation. The effect was positively fed back to the mean field and amplified the anticyclonic circulation accompanied by suppressed storm activity in northeast Asia. The results suggest that subsynoptic-scale mountains affect not only subsynoptic-scale eddies but also the global climate.


2011 ◽  
Vol 4 (4) ◽  
pp. 3047-3065
Author(s):  
R. S. Smith

Abstract. FAMOUS is an ocean-atmosphere general circulation model of low resolution, based on version 4.5 of the UK MetOffice Unified Model. Here we update the model description to account for changes in the model as it is used in the CMIP5 EMIC model intercomparison project (EMICmip) and a number of other studies. Most of these changes correct errors found in the code. The EMICmip version of the model (XFXWB) has a better-conserved water budget and additional cooling in some high latitude areas, but otherwise has a similar climatology to previous versions of FAMOUS. A variant of XFXWB is also described, with changes to the dynamics at the top of the model which improve the model climatology (XFHCC).


2017 ◽  
Vol 24 (4) ◽  
pp. 681-694 ◽  
Author(s):  
Yuxin Zhao ◽  
Xiong Deng ◽  
Shaoqing Zhang ◽  
Zhengyu Liu ◽  
Chang Liu ◽  
...  

Abstract. Climate signals are the results of interactions of multiple timescale media such as the atmosphere and ocean in the coupled earth system. Coupled data assimilation (CDA) pursues balanced and coherent climate analysis and prediction initialization by incorporating observations from multiple media into a coupled model. In practice, an observational time window (OTW) is usually used to collect measured data for an assimilation cycle to increase observational samples that are sequentially assimilated with their original error scales. Given different timescales of characteristic variability in different media, what are the optimal OTWs for the coupled media so that climate signals can be most accurately recovered by CDA? With a simple coupled model that simulates typical scale interactions in the climate system and twin CDA experiments, we address this issue here. Results show that in each coupled medium, an optimal OTW can provide maximal observational information that best fits the characteristic variability of the medium during the data blending process. Maintaining correct scale interactions, the resulting CDA improves the analysis of climate signals greatly. These simple model results provide a guideline for when the real observations are assimilated into a coupled general circulation model for improving climate analysis and prediction initialization by accurately recovering important characteristic variability such as sub-diurnal in the atmosphere and diurnal in the ocean.


2007 ◽  
Vol 4 (5) ◽  
pp. 3413-3440 ◽  
Author(s):  
E. P. Maurer ◽  
H. G. Hidalgo

Abstract. Downscaling of climate model data is essential to most impact analysis. We compare two methods of statistical downscaling to produce continuous, gridded time series of precipitation and surface air temperature at a 1/8-degree (approximately 140 km² per grid cell) resolution over the western U.S. We use NCEP/NCAR Reanalysis data from 1950–1999 as a surrogate General Circulation Model (GCM). The two methods included are constructed analogues (CA) and a bias correction and spatial downscaling (BCSD), both of which have been shown to be skillful in different settings, and BCSD has been used extensively in hydrologic impact analysis. Both methods use the coarse scale Reanalysis fields of precipitation and temperature as predictors of the corresponding fine scale fields. CA downscales daily large-scale data directly and BCSD downscales monthly data, with a random resampling technique to generate daily values. The methods produce comparable skill in producing downscaled, gridded fields of precipitation and temperatures at a monthly and seasonal level. For daily precipitation, both methods exhibit some skill in reproducing both observed wet and dry extremes and the difference between the methods is not significant, reflecting the general low skill in daily precipitation variability in the reanalysis data. For low temperature extremes, the CA method produces greater downscaling skill than BCSD for fall and winter seasons. For high temperature extremes, CA demonstrates higher skill than BCSD in summer. We find that the choice of most appropriate downscaling technique depends on the variables, seasons, and regions of interest, on the availability of daily data, and whether the day to day correspondence of weather from the GCM needs to be reproduced for some applications. The ability to produce skillful downscaled daily data depends primarily on the ability of the climate model to show daily skill.


Sign in / Sign up

Export Citation Format

Share Document