Coupled atmosphere-biophysics-hydrology models for environmental modeling

The formulation and implementation of LEAF-2, the Land Ecosystem–Atmosphere Feedback model, which comprises the representation of land–surface processes in the Regional Atmospheric Modeling System (RAMS), is described. LEAF-2 is a prognostic model for the temperature and water content of soil, snow cover, vegetation, and canopy air, and includes turbulent and radiative exchanges between these components and with the atmosphere. Subdivision of a RAMS surface grid cell into multiple areas of distinct land-use types is allowed, with each subgrid area, or patch, containing its own LEAF-2 model, and each patch interacts with the overlying atmospheric column with a weight proportional to its fractional area in the grid cell. A description is also given of TOPMODEL, a land hydrology model that represents surface and subsurface downslope lateral transport of groundwater. Details of the incorporation of a modified form of TOPMODEL into LEAF-2 are presented. Sensitivity tests of the coupled system are presented that demonstrate the potential importance of the patch representation and of lateral water transport in idealized model simulations. Independent studies that have applied LEAF-2 and verified its performance against observational data are cited. Linkage of RAMS and TOPMODEL through LEAF-2 creates a modeling system that can be used to explore the coupled atmosphere–biophysical–hydrologic response to altered climate forcing at local watershed and regional basin scales.

Cloud-resolving model simulations with one- and two-way couplings via the weak temperature gradient approximation

A cloud-resolving model is modified to implement the weak temperature gradient approximation in order to simulate the interactions between tropical convection and the large-scale tropical circulation. The instantaneous domain-mean potential temperature is relaxed toward a reference profile obtained from a radiative–convective equilibrium simulation of the cloud-resolving model. For homogeneous surface conditions, the model state at equilibrium is a large-scale circulation with its descending branch in the simulated column. This is similar to the equilibrium state found in some other studies, but not all. For this model, the development of such a circulation is insensitive to the relaxation profile and the initial conditions. Two columns of the cloud-resolving model are fully coupled by relaxing the instantaneous domain-mean potential temperature in both columns toward each other. This configuration is energetically closed in contrast to the reference-column configuration. No mean large-scale circulation develops over homogeneous surface conditions, regardless of the relative area of the two columns. The sensitivity to nonuniform surface conditions is similar to that obtained in the reference-column configuration if the two simulated columns have very different areas, but it is markedly weaker for columns of comparable area. The weaker sensitivity can be understood as being a consequence of a formulation for which the energy budget is closed. The reference-column configuration has been used to study the convection in a local region under the influence of a large-scale circulation. The extension to a two-column configuration is proposed as a methodology for studying the influence on local convection of changes in remote convection.

The ENSO-Asian monsoon interaction in a coupled ocean-atmosphere GCM

In this study, the authors evaluate the (El Niño–Southern Oscillation) ENSO–Asian monsoon interaction in a version of the Hadley Centre coupled ocean–atmosphere general circulation model (CGCM) known as HadCM3. The main focus is on two evolving anomalous anticyclones: one located over the south Indian Ocean (SIO) and the other over the western North Pacific (WNP). These two anomalous anticyclones are closely related to the developing and decaying phases of the ENSO and play a crucial role in linking the Asian monsoon to ENSO. It is found that the HadCM3 can well simulate the main features of the evolution of both anomalous anticyclones and the related SST dipoles, in association with the different phases of the ENSO cycle. By using the simulated results, the authors examine the relationship between the WNP/SIO anomalous anticyclones and the ENSO cycle, in particular the biennial component of the relationship. It is found that a strong El Niño event tends to be followed by a more rapid decay and is much more likely to become a La Niña event in the subsequent winter. The twin anomalous anticyclones in the western Pacific in the summer of a decaying El Niño are crucial for the transition from an El Niño into a La Niña. The El Niño (La Niña) events, especially the strong ones, strengthen significantly the correspondence between the SIO anticyclonic (cyclonic) anomaly in the preceding autumn and WNP anticyclonic (cyclonic) anomaly in the subsequent spring, and favor the persistence of the WNP anomaly from spring to summer. The present results suggest that both El Niño (La Niña) and the SIO/WNP anticyclonic (cyclonic) anomalies are closely tied with the tropospheric biennial oscillation (TBO). In addition, variability in the East Asian summer monsoon, which is dominated by the internal atmospheric variability, seems to be responsible for the appearance of the WNP anticyclonic anomaly through an upper-tropospheric meridional teleconnection pattern over the western and central Pacific.

Deriving global flood hazard maps of fluvial floods through a physical model cascade

Global flood hazard maps can be used in the assessment of flood risk in a number of different applications, including (re)insurance and large scale flood preparedness. Such global hazard maps can be generated using large scale physically based models of rainfall-runoff and river routing, when used in conjunction with a number of post-processing methods. In this study, the European Centre for Medium Range Weather Forecasts (ECMWF) land surface model is coupled to ERA-Interim reanalysis meteorological forcing data, and resultant runoff is passed to a river routing algorithm which simulates floodplains and flood flow across the global land area. The global hazard map is based on a 30 yr (1979–2010) simulation period. A Gumbel distribution is fitted to the annual maxima flows to derive a number of flood return periods. The return periods are calculated initially for a 25×25 km grid, which is then reprojected onto a 1×1 km grid to derive maps of higher resolution and estimate flooded fractional area for the individual 25×25 km cells. Several global and regional maps of flood return periods ranging from 2 to 500 yr are presented. The results compare reasonably to a benchmark data set of global flood hazard. The developed methodology can be applied to other datasets on a global or regional scale.

A mechanism for Atlantic multidecadal variability in the Kiel Climate Model

Atlantic Multidecadal Variability (AMV) is investigated in a millennial control simulation with the Kiel Climate Model (KCM), a coupled atmosphere–ocean–sea ice model. An oscillatory mode with approximately 60 years period and characteristics similar to observations is identified with the aid of three-dimensional temperature and salinity joint empirical orthogonal function analysis. The mode explains 30 % of variability on centennial and shorter timescales in the upper 2,000 m of the North Atlantic. It is associated with changes in the Atlantic Meridional Overturning Circulation (AMOC) of ±1–2 Sv and Atlantic Sea Surface Temperature (SST) of ±0.2 °C. AMV in KCM results from an out-of-phase interaction between horizontal and vertical ocean circulation, coupled through Irminger Sea convection. Wintertime convection in this region is mainly controlled by salinity anomalies transported by the Subpolar Gyre (SPG). Increased (decreased) dense water formation in this region leads to a stronger (weaker) AMOC after 15 years, and this in turn leads to a weaker (stronger) SPG after another 15 years. The key role of salinity variations in the subpolar North Atlantic for AMV is confirmed in a 1,000 year long simulation with salinity restored to model climatology: No low frequency variations in convection are simulated, and the 60 year mode of variability is absent.

Updated estimate of aerosol direct radiative forcing from satellite observations and comparison against the Hadley Centre climate model

The fourth assessment report of the Intergovernmental Panel on Climate Change (IPCC) includes a comparison of observation-based and modeling-based estimates of the aerosol direct radiative forcing. In this comparison, satellite-based studies suggest a more negative aerosol direct radiative forcing than modeling studies. A previous satellite-based study, part of the IPCC comparison, uses aerosol optical depths and accumulation-mode fractions retrieved by the Moderate Resolution Imaging Spectroradiometer (MODIS) at collection 4. The latest version of MODIS products, named collection 5, improves aerosol retrievals. Using these products, the direct forcing in the shortwave spectrum defined with respect to present-day natural aerosols is now estimated at −1.30 and −0.65 Wm−2 on a global clear-sky and all-sky average, respectively, for 2002. These values are still significantly more negative than the numbers reported by modeling studies. By accounting for differences between present-day natural and preindustrial aerosol concentrations, sampling biases, and investigating the impact of differences in the zonal distribution of anthropogenic aerosols, good agreement is reached between the direct forcing derived from MODIS and the Hadley Centre climate model HadGEM2-A over clear-sky oceans. Results also suggest that satellite estimates of anthropogenic aerosol optical depth over land should be coupled with a robust validation strategy in order to refine the observation-based estimate of aerosol direct radiative forcing. In addition, the complex problem of deriving the aerosol direct radiative forcing when aerosols are located above cloud still needs to be addressed.

An evaluation of the spatiotemporal structure of large-scale European drought within the HiGEM climate model

We compare the characteristics of synthetic European droughts generated by the HiGEM1 coupled climate model run with present day atmospheric composition with observed drought events extracted from the CRU TS3 data set. The results demonstrate consistency in both the rate of drought occurrence and the spatiotemporal structure of the events. Estimates of the probability density functions for event area, duration and severity are shown to be similar with confidence > 90%. Encouragingly, HiGEM is shown to replicate the extreme tails of the observed distributions and thus the most damaging European drought events. The soil moisture state is shown to play an important role in drought development. Once a large-scale drought has been initiated it is found to be 50% more likely to continue if the local soil moisture is below the 40th percentile. In response to increased concentrations of atmospheric CO2, the modelled droughts are found to increase in duration, area and severity. The drought response can be largely attributed to temperature driven changes in relative humidity. 1 HiGEM is based on the latest climate configuration of the Met Office Hadley Centre Unified Model (HadGEM1) with the horizontal resolution increased to 1.25 x 0.83 degrees in longitude and latitude in the atmosphere and 1/3 x 1/3 degrees in the ocean.

Chemistry‐climate model simulations of spring Antarctic ozone

Coupled chemistry‐climate model simulations covering the recent past and continuing throughout the 21st century have been completed with a range of different models. Common forcings are used for the halogen amounts and greenhouse gas concentrations, as expected under the Montreal Protocol (with amendments) and Intergovernmental Panel on Climate Change A1b Scenario. The simulations of the Antarctic ozone hole are compared using commonly used diagnostics: the minimum ozone, the maximum area of ozone below 220 DU, and the ozone mass deficit below 220 DU. Despite the fact that the processes responsible for ozone depletion are reasonably well understood, a wide range of results is obtained. Comparisons with observations indicate that one of the reasons for the model underprediction in ozone hole area is the tendency for models to underpredict, by up to 35%, the area of low temperatures responsible for polar stratospheric cloud formation. Models also typically have species gradients that are too weak at the edge of the polar vortex, suggesting that there is too much mixing of air across the vortex edge. Other models show a high bias in total column ozone which restricts the size of the ozone hole (defined by a 220 DU threshold). The results of those models which agree best with observations are examined in more detail. For several models the ozone hole does not disappear this century but a small ozone hole of up to three million square kilometers continues to occur in most springs even after 2070.

Secular trends in daily precipitation characteristics: greenhouse gas simulation with a coupled AOGCM

Secular trends of daily precipitation characteristics are considered in the transient climate change experiment with a coupled atmosphere-ocean general circulation model ECHAM4/OPYC3 for 1900-2099. The climate forcing is due to increasing concentrations of the greenhouse gases in the atmosphere. Mean daily precipitation, precipitation intensity, probability of wet days and parameters of the gamma distribution are analyzed. Particular attention is paid to the changes of heavy precipitation, Analysis of the annual mean precipitation trends for 1900-1999 revealed general agreement with observations with significant positive trends in mean precipitation over continental areas. In the 2000-2099 period precipitation trend patterns followed the tendency obtained for 1900-1999 but with significantly increased magnitudes. Unlike the annual mean precipitation trends for which negative values were found for some continental areas, the mean precipitation intensity and scale parameter of the fitted gamma distribution increased over all land territories . Negative trends in the number of wet days were found over most of the land areas except high latitudes in the Northern Hemisphere. The shape parameter of the gamma distribution in general revealed a slight negative trend in the areas of the precipitation increase. Investigation of daily precipitation revealed an unproportional increase of heavy precipitation events for the land areas including local maxima in Europe and the eastern United States.

Impact of global warming on the interannual and interdecadal climate modes in a coupled GCM

In this study, we investigated the impact of global warming on the variabilities of large-scale interannual and interdecadal climate modes and teleconnection patterns with two long-term integrations of the coupled general circulation model of ECHAM4/OPYC3 at the Max-Planck-Institute for Meteorology, Hamburg. One is the control (CTRL) run with fixed present-day concentrations of greenhouse gases. The other experiment is a simulation of transient greenhouse warming, named GHG run. In the GHG run the averaged geopotential height at 500 hPa is increased significantly, and a negative phase of the Pacific/North American (PNA) teleconnection-like distribution pattern is intensified. The standard deviation over the tropics (high latitudes) is enhanced (reduced) on the interdecadal time scales and reduced (enhanced) on the interannual time scales in the GHG run. Except for an interdecadal mode related to the Southern Oscillation (SO) in the GHG run, the spatial variation patterns are similar for different (interannual + interdecadal, interannual, and interdecadal) time scales in the GHG and CTRL runs. Spatial distributions of the teleconnection patterns on the interannual and interdecadal time scales in the GHG run are also similar to those in the CTRL run. But some teleconnection patterns show linear trends and changes of variances and frequencies in the GHG run. Apart from the positive linear trend of the SO, the interdecadal modulation to the El Niño/SO cycle is enhanced during the GHG 2040 ∼ 2099. This is the result of an enhancement of the Walker circulation during that period. La Niña events intensify and El Niño events relatively weaken during the GHG 2070 ∼ 2090. It is interesting to note that with increasing greenhouse gas concentrations the relation between the SO and the PNA pattern is reversed significantly from a negative to a positive correlation on the interdecadal time scales and weakened on the interannual time scales. This suggests that the increase of the greenhouse gas concentrations will trigger the nonstationary correlation between the SO and the PNA pattern both on the interdecadal and interannual time scales.

A low-order model investigation of the analysis of gravity waves in the ensemble Kalman filter.

The behavior of the ensemble Kalman filter (EnKF) is examined in the context of a model that exhibits a nonlinear chaotic (slow) vortical mode coupled to a linear (fast) gravity wave of a given amplitude and frequency. It is shown that accurate recovery of both modes is enhanced when covariances between fast and slow normal-mode variables (which reflect the slaving relations inherent in balanced dynamics) are modeled correctly. More ensemble members are needed to recover the fast, linear gravity wave than the slow, vortical motion. Although the EnKF tends to diverge in the analysis of the gravity wave, the filter divergence is stable and does not lead to a great loss of accuracy. Consequently, provided the ensemble is large enough and observations are made that reflect both time scales, the EnKF is able to recover both time scales more accurately than optimal interpolation (OI), which uses a static error covariance matrix. For OI it is also found to be problematic to observe the state at a frequency that is a subharmonic of the gravity wave frequency, a problem that is in part overcome by the EnKF.However, error in themodeled gravity wave parameters can be detrimental to the performance of the EnKF and remove its implied advantages, suggesting that a modified algorithm or a method for accounting for model error is needed.

Assessment of temperature, trace species, and ozone in chemistry-climate model simulations of the recent past

Simulations of the stratosphere from thirteen coupled chemistry-climate models (CCMs) are evaluated to provide guidance for the interpretation of ozone predictions made by the same CCMs. The focus of the evaluation is on how well the fields and processes that are important for determining the ozone distribution are represented in the simulations of the recent past. The core period of the evaluation is from 1980 to 1999 but long-term trends are compared for an extended period (1960–2004). Comparisons of polar high-latitude temperatures show that most CCMs have only small biases in the Northern Hemisphere in winter and spring, but still have cold biases in the Southern Hemisphere spring below 10 hPa. Most CCMs display the correct stratospheric response of polar temperatures to wave forcing in the Northern, but not in the Southern Hemisphere. Global long-term stratospheric temperature trends are in reasonable agreement with satellite and radiosonde observations. Comparisons of simulations of methane, mean age of air, and propagation of the annual cycle in water vapor show a wide spread in the results, indicating differences in transport. However, for around half the models there is reasonable agreement with observations. In these models the mean age of air and the water vapor tape recorder signal are generally better than reported in previous model intercomparisons. Comparisons of the water vapor and inorganic chlorine (Cly) fields also show a large intermodel spread. Differences in tropical water vapor mixing ratios in the lower stratosphere are primarily related to biases in the simulated tropical tropopause temperatures and not transport. The spread in Cly, which is largest in the polar lower stratosphere, appears to be primarily related to transport differences. In general the amplitude and phase of the annual cycle in total ozone is well simulated apart from the southern high latitudes. Most CCMs show reasonable agreement with observed total ozone trends and variability on a global scale, but a greater spread in the ozone trends in polar regions in spring, especially in the Arctic. In conclusion, despite the wide range of skills in representing different processes assessed here, there is sufficient agreement between the majority of the CCMs and the observations that some confidence can be placed in their predictions.

The simulation of the opposing fluxes of latent heat and CO2 over various land-use types: coupling a gas exchange model to a mesoscale atmospheric model

A mesoscale meteorological model (FOOT3DK) is coupled with a gas exchange model to simulate surface fluxes of CO2 and H2O under field conditions. The gas exchange model consists of a C3 single leaf photosynthesis sub-model and an extended big leaf (sun/shade) sub-model that divides the canopy into sunlit and shaded fractions. Simulated CO2 fluxes of the stand-alone version of the gas exchange model correspond well to eddy-covariance measurements at a test site in a rural area in the west of Germany. The coupled FOOT3DK/gas exchange model is validated for the diurnal cycle at singular grid points, and delivers realistic fluxes with respect to their order of magnitude and to the general daily course. Compared to the Jarvis-based big leaf scheme, simulations of latent heat fluxes with a photosynthesis-based scheme for stomatal conductance are more realistic. As expected, flux averages are strongly influenced by the underlying land cover. While the simulated net ecosystem exchange is highly correlated with leaf area index, this correlation is much weaker for the latent heat flux. Photosynthetic CO2 uptake is associated with transpirational water loss via the stomata, and the resulting opposing surface fluxes of CO2 and H2O are reproduced with the model approach. Over vegetated surfaces it is shown that the coupling of a photosynthesis-based gas exchange model with the land-surface scheme of a mesoscale model results in more realistic simulated latent heat fluxes.

European storminess and associated circulation weather types: future changes deduced from a multi-model ensemble of GCM simulations

A range of possible changes in the frequency and characteristics of European wind storms under future climate conditions was investigated on the basis of a multi-model ensemble of 9 coupled global climate model (GCM) simulations for the 20th and 21st centuries following the IPCC SRES A1B scenario. A multi-model approach allowed an estimation of the (un)certainties of the climate change signals. General changes in large-scale atmospheric flow were analysed, the occurrence of wind storms was quantified, and atmospheric features associated with wind storm events were considered. Identified storm days were investigated according to atmospheric circulation, associated pressure patterns, cyclone tracks and wind speed patterns. Validation against reanalysis data revealed that the GCMs are in general capable of realistically reproducing characteristics of European circulation weather types (CWTs) and wind storms. Results are given with respect to frequency of occurrence, storm-associated flow conditions, cyclone tracks and specific wind speed patterns. Under anthropogenic climate change conditions (SRES A1B scenario), increased frequency of westerly flow during winter is detected over the central European investigation area. In the ensemble mean, the number of detected wind storm days increases between 19 and 33% for 2 different measures of storminess, only 1 GCM revealed less storm days. The increased number of storm days detected in most models is disproportionately high compared to the related CWT changes. The mean intensity of cyclones associated with storm days in the ensemble mean increases by about 10 (±10)% in the Eastern Atlantic, near the British Isles and in the North Sea. Accordingly, wind speeds associated with storm events increase significantly by about 5 (±5)% over large parts of central Europe, mainly on days with westerly flow. The basic conclusions of this work remain valid if different ensemble contructions are considered, leaving out an outlier model or including multiple runs of one particular model.

Averaging, slaving and balance dynamics in a simple atmospheric model

We report numerical results from a study of balance dynamics using a simple model of atmospheric motion that is designed to help address the question of why balance dynamics is so stable. The non-autonomous Hamiltonian model has a chaotic slow degree of freedom (representing vortical modes) coupled to one or two linear fast oscillators (representing inertia-gravity waves). The system is said to be balanced when the fast and slow degrees of freedom are separated. We find adiabatic invariants that drift slowly in time. This drift is consistent with a random-walk behaviour at a speed which qualitatively scales, even for modest time scale separations, as the upper bound given by Neishtadt’s and Nekhoroshev’s theorems. Moreover, a similar type of scaling is observed for solutions obtained using a singular perturbation (‘slaving’) technique in resonant cases where Nekhoroshev’s theorem does not apply. We present evidence that the smaller Lyapunov exponents of the system scale exponentially as well. The results suggest that the observed stability of nearly-slow motion is a consequence of the approximate adiabatic invariance of the fast motion.