Developing the next-generation climate system models: challenges and achievements

Although climate models have been improving in accuracy and efficiency over the past few decades, it now seems that these incremental improvements may be slowing. As tera/petascale computing becomes massively parallel, our legacy codes are less suitable, and even with the increased resolution that we are now beginning to use, these models cannot represent the multiscale nature of the climate system. This paper argues that it may be time to reconsider the use of adaptive mesh refinement for weather and climate forecasting in order to achieve good scaling and representation of the wide range of spatial scales in the atmosphere and ocean. Furthermore, the challenge of introducing living organisms and human responses into climate system models is only just beginning to be tackled. We do not yet have a clear framework in which to approach the problem, but it is likely to cover such a huge number of different scales and processes that radically different methods may have to be considered. The challenges of multiscale modelling and petascale computing provide an opportunity to consider a fresh approach to numerical modelling of the climate (or Earth) system, which takes advantage of the computational fluid dynamics developments in other fields and brings new perspectives on how to incorporate Earth system processes. This paper reviews some of the current issues in climate (and, by implication, Earth) system modelling, and asks the question whether a new generation of models is needed to tackle these problems.

Surface wave processes in air-sea interaction

We review briefly recent progress on understanding the role of surface waves on the marine atmospheric boundary layer and the ocean mixed layer and give a global perspective on these processes by analysing ERA-40 data. Ocean surface waves interact with the marine atmospheric boundary layer in two broad regimes: (i) the conventional wind-driven wave regime, when fast winds blow over slower moving waves, and (ii) a wave-driven wind regime when long wavelength swell propagates under low winds, and generates a wave-driven jet in the lower part of the marine boundary layer. Analysis of ERA-40 data indicates that the wave-driven wind regime is as prevalent as the conventional wind-driven regime. Ocean surface waves also change profoundly mixing in the ocean mixed layer through generation of Langmuir circulation. Results from large-eddy simulation are used here to develop a scaling for the resulting Langmuir turbulence, which is a necessary step in developing a parametrization of the process. ERA-40 data is then used to show that the Langmuir regime is the predominant regime over much of the global ocean, providing a compelling motivation for parameterising this process in ocean general circulation models.

Shifts in rhizosphere microbial communities and enzyme activity of Poa alpina across an alpine chronosequence

This study quantifies the influence of Poa alpina on the soil microbial community in primary succession of alpine ecosystems, and whether these effects are controlled by the successional stage. Four successional sites representative of four stages of grassland development (initial, 4 years (non-vegetated); pioneer, 20 years; transition, 75 years; mature, 9500 years old) on the Rotmoos glacier foreland, Austria, were sampled. The size, composition and activity of the microbial community in the rhizosphere and bulk soil were characterized using the chloroform-fumigation extraction procedure, phospholipid fatty acid (PLFA) analysis and measurements of the enzymes beta-glucosidase, beta-xylosidase, N-acetyl-beta-glucosaminidase, leucine aminopeptidase, acid phosphatase and sulfatase. The interplay between the host plant and the successional stage was quantified using principal component (PCA) and multidimensional scaling analyses. Correlation analyses were applied to evaluate the relationship between soil factors (C-org, N-t, C/N ratio, pH, ammonium, phosphorus, potassium) and microbial properties in the bulk soil. In the pioneer stage microbial colonization of the rhizosphere of P. alpina was dependent on the reservoir of microbial species in the bulk soil. As a consequence, the rhizosphere and bulk soil were similar in microbial biomass (ninhydrin-reactive nitrogen (NHR-N)), community composition (PLFA), and enzyme activity. In the transition and mature grassland stage, more benign soil conditions stimulated microbial growth (NHR-N, total amount of PLFA, bacterial PLFA, Gram-positive bacteria, Gram-negative bacteria), and microbial diversity (Shannon index H) in the rhizosphere either directly or indirectly through enhanced carbon allocation. In the same period, the rhizosphere microflora shifted from a G(-) to a more G(+), and from a fungal to a more bacteria-dominated community. Rhizosphere beta-xylosidase, N-acetyl-beta-glucosaminidase, and sulfatase activity peaked in the mature grassland soil, whereas rhizosphere leucine aminopeptidase, beta-glucosidase, and phosphatase activity were highest in the transition stage, probably because of enhanced carbon and nutrient allocation into the rhizosphere due to better growth conditions. Soil organic matter appeared to be the most important driver of microbial colonization in the bulk soil. The decrease in soil pH and soil C/N ratio mediated the shifts in the soil microbial community composition (bacPLFA, bacPLFA/fungPLFA, G(-), G(+)/G(-)). The activities of beta-glucosidase, beta-xylosidase and phosphatase were related to soil ammonium and phosphorus, indicating that higher decomposition rates enhanced the nutrient availability in the bulk soil. We conclude that the major determinants of the microllora vary along the successional gradient: in the pioneer stage the rhizosphere microflora was primarily determined by the harsh soil environment; under more favourable environmental conditions, however, the host plant selected for a specific microbial community that was related to the dynamic interplay between soil properties and carbon supply. (C) 2004 Elsevier Ltd. All rights reserved.

The use of earth observation techniques to improve catchment-scale pollution predictions

Remote sensing can potentially provide information useful in improving pollution transport modelling in agricultural catchments. Realisation of this potential will depend on the availability of the raw data, development of information extraction techniques, and the impact of the assimilation of the derived information into models. High spatial resolution hyperspectral imagery of a farm near Hereford, UK is analysed. A technique is described to automatically identify the soil and vegetation endmembers within a field, enabling vegetation fractional cover estimation. The aerially-acquired laser altimetry is used to produce digital elevation models of the site. At the subfield scale the hypothesis that higher resolution topography will make a substantial difference to contaminant transport is tested using the AGricultural Non-Point Source (AGNPS) model. Slope aspect and direction information are extracted from the topography at different resolutions to study the effects on soil erosion, deposition, runoff and nutrient losses. Field-scale models are often used to model drainage water, nitrate and runoff/sediment loss, but the demanding input data requirements make scaling up to catchment level difficult. By determining the input range of spatial variables gathered from EO data, and comparing the response of models to the range of variation measured, the critical model inputs can be identified. Response surfaces to variation in these inputs constrain uncertainty in model predictions and are presented. Although optical earth observation analysis can provide fractional vegetation cover, cloud cover and semi-random weather patterns can hinder data acquisition in Northern Europe. A Spring and Autumn cloud cover analysis is carried out over seven UK sites close to agricultural districts, using historic satellite image metadata, climate modelling and historic ground weather observations. Results are assessed in terms of probability of acquisition probability and implications for future earth observation missions. (C) 2003 Elsevier Ltd. All rights reserved.

Shifts in rhizosphere microbial communities and enzyme activity of Poa alpina across an alpine chronosequence

This study quantifies the influence of Poa alpina on the soil microbial community in primary succession of alpine ecosystems, and whether these effects are controlled by the successional stage. Four successional sites representative of four stages of grassland development (initial, 4 years (non-vegetated); pioneer, 20 years; transition, 75 years; mature, 9500 years old) on the Rotmoos glacier foreland, Austria, were sampled. The size, composition and activity of the microbial community in the rhizosphere and bulk soil were characterized using the chloroform-fumigation extraction procedure, phospholipid fatty acid (PLFA) analysis and measurements of the enzymes beta-glucosidase, beta-xylosidase, N-acetyl-beta-glucosaminidase, leucine aminopeptidase, acid phosphatase and sulfatase. The interplay between the host plant and the successional stage was quantified using principal component (PCA) and multidimensional scaling analyses. Correlation analyses were applied to evaluate the relationship between soil factors (C-org, N-t, C/N ratio, pH, ammonium, phosphorus, potassium) and microbial properties in the bulk soil. In the pioneer stage microbial colonization of the rhizosphere of P. alpina was dependent on the reservoir of microbial species in the bulk soil. As a consequence, the rhizosphere and bulk soil were similar in microbial biomass (ninhydrin-reactive nitrogen (NHR-N)), community composition (PLFA), and enzyme activity. In the transition and mature grassland stage, more benign soil conditions stimulated microbial growth (NHR-N, total amount of PLFA, bacterial PLFA, Gram-positive bacteria, Gram-negative bacteria), and microbial diversity (Shannon index H) in the rhizosphere either directly or indirectly through enhanced carbon allocation. In the same period, the rhizosphere microflora shifted from a G(-) to a more G(+), and from a fungal to a more bacteria-dominated community. Rhizosphere beta-xylosidase, N-acetyl-beta-glucosaminidase, and sulfatase activity peaked in the mature grassland soil, whereas rhizosphere leucine aminopeptidase, beta-glucosidase, and phosphatase activity were highest in the transition stage, probably because of enhanced carbon and nutrient allocation into the rhizosphere due to better growth conditions. Soil organic matter appeared to be the most important driver of microbial colonization in the bulk soil. The decrease in soil pH and soil C/N ratio mediated the shifts in the soil microbial community composition (bacPLFA, bacPLFA/fungPLFA, G(-), G(+)/G(-)). The activities of beta-glucosidase, beta-xylosidase and phosphatase were related to soil ammonium and phosphorus, indicating that higher decomposition rates enhanced the nutrient availability in the bulk soil. We conclude that the major determinants of the microllora vary along the successional gradient: in the pioneer stage the rhizosphere microflora was primarily determined by the harsh soil environment; under more favourable environmental conditions, however, the host plant selected for a specific microbial community that was related to the dynamic interplay between soil properties and carbon supply. (C) 2004 Elsevier Ltd. All rights reserved.

Stratospheric temperature and radiative forcing response to 11-year solar cycle changes in irradiance and ozone

The 11-yr solar cycle temperature response to spectrally resolved solar irradiance changes and associated ozone changes is calculated using a fixed dynamical heating (FDH) model. Imposed ozone changes are from satellite observations, in contrast to some earlier studies. A maximum of 1.6 K is found in the equatorial upper stratosphere and a secondary maximum of 0.4 K in the equatorial lower stratosphere, forming a double peak in the vertical. The upper maximum is primarily due to the irradiance changes while the lower maximum is due to the imposed ozone changes. The results compare well with analyses using the 40-yr ECMWF Re-Analysis (ERA-40) and NCEP/NCAR datasets. The equatorial lower stratospheric structure is reproduced even though, by definition, the FDH calculations exclude dynamically driven temperature changes, suggesting an important role for an indirect dynamical effect through ozone redistribution. The results also suggest that differences between the Stratospheric Sounding Unit (SSU)/Microwave Sounding Unit (MSU) and ERA-40 estimates of the solar cycle signal can be explained by the poor vertical resolution of the SSU/MSU measurements. The adjusted radiative forcing of climate change is also investigated. The forcing due to irradiance changes was 0.14 W m−2, which is only 78% of the value obtained by employing the standard method of simple scaling of the total solar irradiance (TSI) change. The difference arises because much of the change in TSI is at wavelengths where ozone absorbs strongly. The forcing due to the ozone change was only 0.004 W m−2 owing to strong compensation between negative shortwave and positive longwave forcings.

The Effects of Stratification on Flow Separation

Separation of stratified flow over a two-dimensional hill is inhibited or facilitated by acceleration or deceleration of the flow just outside the attached boundary layer. In this note, an expression is derived for this acceleration or deceleration in terms of streamline curvature and stratification. The expression is valid for linear as well as nonlinear deformation of the flow. For hills of vanishing aspect ratio a linear theory can be derived and a full regime diagram for separation can be constructed. For hills of finite aspect ratio scaling relationships can be derived that indicate the presence of a critical aspect ratio, proportional to the stratification, above which separation will occur as well as a second critical aspect ratio above which separation will always occur irrespective of stratification.

The three-dimensional vortical nature of atmospheric and oceanic turbulent flows

Using a novel numerical method at unprecedented resolution, we demonstrate that structures of small to intermediate scale in rotating, stratified flows are intrinsically three-dimensional. Such flows are characterized by vortices (spinning volumes of fluid), regions of large vorticity gradients, and filamentary structures at all scales. It is found that such structures have predominantly three-dimensional dynamics below a horizontal scale LLR, where LR is the so-called Rossby radius of deformation, equal to the characteristic vertical scale of the fluid H divided by the ratio of the rotational and buoyancy frequencies f/N. The breakdown of two-dimensional dynamics at these scales is attributed to the so-called "tall-column instability" [D. G. Dritschel and M. de la Torre Juárez, J. Fluid. Mech. 328, 129 (1996)], which is active on columnar vortices that are tall after scaling by f/N, or, equivalently, that are narrow compared with LR. Moreover, this instability eventually leads to a simple relationship between typical vertical and horizontal scales: for each vertical wave number (apart from the vertically averaged, barotropic component of the flow) the average horizontal wave number is equal to f/N times the vertical wave number. The practical implication is that three-dimensional modeling is essential to capture the behavior of rotating, stratified fluids. Two-dimensional models are not valid for scales below LR. ©1999 American Institute of Physics.

The cumulus-capped boundary layer. II: Interface fluxes

This paper considers the relationship between the mean temperature and humidity profiles and the fluxes of heat and moisture at cloud base and the base of the inversion in the cumulus-capped boundary layer. The relationships derived are based on an approximate form of the scalar-flux budget and the scaling properties of the turbulent kinetic energy (TKE) budget. The scalar-flux budget gives a relationship between the change in the virtual potential temperature across either the cloud base transition zone or the inversion and the flux at the base of the layer. The scaling properties of the TKE budget lead to a relationship between the heat and moisture fluxes and the mean subsaturation through the liquid-water flux. The 'jump relation' for the virtual potential temperature at cloud base shows the close connection between the cumulus mass flux in the cumulus-capped boundary layer and the entrainment velocity in the dry-convective boundary layer. Gravity waves are shown to be an important feature of the inversion.

Current changes in tropical precipitation

Current changes in tropical precipitation from satellite data and climate models are assessed. Wet and dry regions of the tropics are defined as the highest 30% and lowest 70% of monthly precipitation values. Observed tropical ocean trends in the wet regime (1.8%/decade) and the dry regions (−2.6%/decade) according to the Global Precipitation Climatology Project (GPCP) over the period including Special Sensor Microwave Imager (SSM/I) data (1988–2008), where GPCP is believed to be more reliable, are of smaller magnitude than when including the entire time series (1979–2008) and closer to model simulations than previous comparisons. Analysing changes in extreme precipitation using daily data within the wet regions, an increase in the frequency of the heaviest 6% of events with warming for the SSM/I observations and model ensemble mean is identified. The SSM/I data indicate an increased frequency of the heaviest events with warming, several times larger than the expected Clausius–Clapeyron scaling and at the upper limit of the substantial range in responses in the model simulations.

Effectiveness of landmark analysis for establishing locality in p2p networks

Locality to other nodes on a peer-to-peer overlay network can be established by means of a set of landmarks shared among the participating nodes. Each node independently collects a set of latency measures to landmark nodes, which are used as a multi-dimensional feature vector. Each peer node uses the feature vector to generate a unique scalar index which is correlated to its topological locality. A popular dimensionality reduction technique is the space filling Hilbert’s curve, as it possesses good locality preserving properties. However, there exists little comparison between Hilbert’s curve and other techniques for dimensionality reduction. This work carries out a quantitative analysis of their properties. Linear and non-linear techniques for scaling the landmark vectors to a single dimension are investigated. Hilbert’s curve, Sammon’s mapping and Principal Component Analysis have been used to generate a 1d space with locality preserving properties. This work provides empirical evidence to support the use of Hilbert’s curve in the context of locality preservation when generating peer identifiers by means of landmark vector analysis. A comparative analysis is carried out with an artificial 2d network model and with a realistic network topology model with a typical power-law distribution of node connectivity in the Internet. Nearest neighbour analysis confirms Hilbert’s curve to be very effective in both artificial and realistic network topologies. Nevertheless, the results in the realistic network model show that there is scope for improvements and better techniques to preserve locality information are required.

A fuzzy approach for the network congestion problem

In the recent years, the unpredictable growth of the Internet has moreover pointed out the congestion problem, one of the problems that historicallyha ve affected the network. This paper deals with the design and the evaluation of a congestion control algorithm which adopts a FuzzyCon troller. The analogyb etween Proportional Integral (PI) regulators and Fuzzycon trollers is discussed and a method to determine the scaling factors of the Fuzzycon troller is presented. It is shown that the Fuzzycon troller outperforms the PI under traffic conditions which are different from those related to the operating point considered in the design.

Marine ecosystem models for earth systems applications: the MarQUEST experience

The MarQUEST (Marine Biogeochemistry and Ecosystem Modelling Initiative in QUEST) project was established to develop improved descriptions of marine biogeochemistry, suited for the next generation of Earth system models. We review progress in these areas providing insight on the advances that have been made as well as identifying remaining key outstanding gaps for the development of the marine component of next generation Earth system models. The following issues are discussed and where appropriate results are presented; the choice of model structure, scaling processes from physiology to functional types, the ecosystem model sensitivity to changes in the physical environment, the role of the coastal ocean and new methods for the evaluation and comparison of ecosystem and biogeochemistry models. We make recommendations as to where future investment in marine ecosystem modelling should be focused, highlighting a generic software framework for model development, improved hydrodynamic models, and better parameterisation of new and existing models, reanalysis tools and ensemble simulations. The final challenge is to ensure that experimental/observational scientists are stakeholders in the models and vice versa.

On the steadiness of separating meandering currents

The existence of inertial steady currents that separate from a coast and meander afterward is investigated. By integrating the zonal momentum equation over a suitable area, it is shown that retroflecting currents cannot be steady in a reduced gravity or in a barotropic model of the ocean. Even friction cannot negate this conclusion. Previous literature on this subject, notably the discrepancy between several articles by Nof and Pichevin on the unsteadiness of retroflecting currents and steady solutions presented in other papers, is critically discussed. For more general separating current systems, a local analysis of the zonal momentum balance shows that given a coastal current with a specific zonal momentum structure, an inertial, steady, separating current is unlikely, and the only analytical solution provided in the literature is shown to be inconsistent. In a basin-wide view of these separating current systems, a scaling analysis reveals that steady separation is impossible when the interior flow is nondissipative (e.g., linear Sverdrup-like). These findings point to the possibility that a large part of the variability in the world’s oceans is due to the separation process rather than to instability of a free jet.

Instability of surface-temperature filaments in strain and shear

The effects of uniform straining and shearing on the stability of a surface quasi-geostrophic temperature filament are investigated. Straining is shown to stabilize perturbations for wide filaments but only for a finite time until the filament thins to a critical width, after which some perturbations can grow. No filament can be stabilized in practice, since there are perturbations that can grow large for any strain rate. The optimally growing perturbations, defined as solutions that reach a certain threshold amplitude first, are found numerically for a wide range of parameter values. The radii of the vortices formed through nonlinear roll-up are found to be proportional to θ/s, where θ is the temperature anomaly of the filament and s the strain rate, and are not dependent on the initial size of the filament. Shearing is shown to reduce the normal-mode growth rates, but it cannot stabilize them completely when there are temperature discontinuities in the basic state; smooth filaments can be stabilized completely by shearing and a simple scaling argument provides the shear rate required. Copyright © 2010 Royal Meteorological Society