Soil moisture mediates association between the winter North Atlantic Oscillation and summer growth in the Park Grass Experiment

Several aspects of terrestrial ecosystems are known to be associated with the North Atlantic Oscillation (NAO) through effects of the NAO on winter climate, but recently the winter NAO has also been shown to be correlated with the following summer climate, including drought. Since drought is a major factor determining grassland primary productivity, the hypothesis was tested that the winter NAO is associated with summer herbage growth through soil moisture availability, using data from the Park Grass Experiment at Rothamsted, UK between 1960 and 1999. The herbage growth rate, mean daily rainfall, mean daily potential evapotranspiration (PE) and the mean and maximum potential soil moisture deficit (PSMD) were calculated between the two annual cuts in early summer and autumn for the unlimed, unfertilized plots. Mean and maximum PSMD were more highly correlated than rainfall or PE with herbage growth rate. Regression analysis showed that the natural logarithm of the herbage growth rate approximately halved for a 250 mm increase in maximum PSMD over the range 50-485 mm. The maximum PSMD was moderately correlated with the preceding winter NAO, with a positive winter NAO index associated with greater maximum PSMD. A positive winter NAO index was also associated with low herbage growth rate, accounting for 22% of the interannual variation in the growth rate. It was concluded that the association between the winter NAO and summer herbage growth rate is mediated by the PSMD in summer.

The meteorology of the Western Indian Ocean, and the influence of the east African highlands

This paper reviews the meteorology of the Western Indian Ocean and uses a state–of–the–art atmospheric general circulation model to investigate the influence of the East African Highlands on the climate of the Indian Ocean and its surrounding regions. The new 44–year re–analysis produced by the European Centre for Medium range Weather Forecasts (ECMWF) has been used to construct a new climatology of the Western Indian Ocean. A brief overview of the seasonal cycle of the Western Indian Ocean is presented which emphasizes the importance of the geography of the Indian Ocean basin for controlling the meteorology of the Western Indian Ocean. The principal modes of inter–annual variability are described, associated with El Niño and the Indian Ocean Dipole or Zonal Mode, and the basic characteristics of the subseasonal weather over the Western Indian Ocean are presented, including new statistics on cyclone tracks derived from the ECMWF re–analyses. Sensitivity experiments, in which the orographic effects of East Africa are removed, have shown that the East African Highlands, although not very high, play a significant role in the climate of Africa, India and Southeast Asia, and in the heat, salinity and momentum forcing of the Western Indian Ocean. The hydrological cycle over Africa is systematically enhanced in all seasons by the presence of the East African Highlands, and during the Asian summer monsoon there is a major redistribution of the rainfall across India and Southeast Asia. The implied impact of the East African Highlands on the ocean is substantial. The East African Highlands systematically freshen the tropical Indian Ocean, and act to focus the monsoon winds along the coast, leading to greater upwelling and cooler sea–surface temperatures.

Evidence of trends in daily climate extremes over southern and west Africa

[ 1] There has been a paucity of information on trends in daily climate and climate extremes, especially from developing countries. We report the results of the analysis of daily temperature ( maximum and minimum) and precipitation data from 14 south and west African countries over the period 1961 - 2000. Data were subject to quality control and processing into indices of climate extremes for release to the global community. Temperature extremes show patterns consistent with warming over most of the regions analyzed, with a large proportion of stations showing statistically significant trends for all temperature indices. Over 1961 to 2000, the regionally averaged occurrence of extreme cold ( fifth percentile) days and nights has decreased by - 3.7 and - 6.0 days/decade, respectively. Over the same period, the occurrence of extreme hot (95th percentile) days and nights has increased by 8.2 and 8.6 days/decade, respectively. The average duration of warm ( cold) has increased ( decreased) by 2.4 (0.5) days/decade and warm spells. Overall, it appears that the hot tails of the distributions of daily maximum temperature have changed more than the cold tails; for minimum temperatures, hot tails show greater changes in the NW of the region, while cold tails have changed more in the SE and east. The diurnal temperature range (DTR) does not exhibit a consistent trend across the region, with many neighboring stations showing opposite trends. However, the DTR shows consistent increases in a zone across Namibia, Botswana, Zambia, and Mozambique, coinciding with more rapid increases in maximum temperature than minimum temperature extremes. Most precipitation indices do not exhibit consistent or statistically significant trends across the region. Regionally averaged total precipitation has decreased but is not statistically significant. At the same time, there has been a statistically significant increase in regionally averaged daily rainfall intensity and dry spell duration. While the majority of stations also show increasing trends for these two indices, only a few of these are statistically significant. There are increasing trends in regionally averaged rainfall on extreme precipitation days and in maximum annual 5-day and 1-day rainfall, but only trends for the latter are statistically significant.

Indian Ocean-monsoon coupled interactions and impending monsoon droughts

Monsoon droughts over the Indian subcontinent emanate from failures in the seasonal (June-September) monsoon rains. While prolonged dry-spells ("monsoon-breaks'') pervade on sub-seasonal/intra-seasonal time-scales, the underlying causes for these long-lasting anomalies remain elusive. Based on analyses of a suite of observed data sets, we report an ocean-atmosphere dynamical coupling on intra-seasonal time-scales, in the tropical Indian Ocean, which is pivotal in forcing extended monsoon-breaks and causing droughts over the subcontinent. This coupling involves a feedback between the monsoonal flow and thermocline depth in the Equatorial Eastern Indian Ocean (EEIO), in which an anomaly of the summer monsoon circulation induces downwelling and maintains a higher-than-normal heat-content. The near-equatorial anomalies induce strong and sustained suppression of monsoon rainfall over the subcontinent. It is concluded that the intra-seasonal evolution of the ocean-monsoon coupled system is a vital key to unlocking the dynamics of monsoon droughts.

Recent trends in sea level pressure in the Indian Ocean region

During the second half of the twentieth century the Indian Ocean exhibited a rapid rise in sea surface temperatures (SST). It has been argued - largely on the basis of experiments with atmospheric GCMs - that this rapid warming was an important cause of remote changes in climate, in particular an increasing trend in the North Atlantic Oscillation Index and decreases in African rainfall. Here however we present evidence that the Indian Ocean warming was associated with local increases in sea level pressure (SLP). These increases are inconsistent with results from experiments in which an atmospheric GCM is forced by historical SST, which show robust decreases in SLP. The clear discrepancy between the observed and simulated trends in SLP suggests that the response of some atmospheric GCMs to the Indian Ocean warming may not provide a reliable guide to the behaviour of the real world.

The effect of soil phosphorus on particulate phosphorus in land runoff

Accumulation of surplus phosphorus (P) in the soil and the resulting increased transport of P in land runoff contribute to freshwater eutrophication. The effects of increasing soil P (19–194 mg Olsen-P (OP) kg−1) on the concentrations of particulate P (PP), and sorption properties (Qmax, k and EPCo) of suspended solids (SS) in overland flow from 15 unreplicated field plots established on a dispersive arable soil were measured over three monitoring periods under natural rainfall. Concentrations of PP in plot runoff increased linearly at a rate of 2.6 μg litre−1 per mg OP kg−1 of soil, but this rate was approximately 50% of the rate of increase in dissolved P (< 0.45 μm). Concentrations of SS in runoff were similar across all plots and contained a greater P sorption capacity (mean + 57%) than the soil because of enrichment with fine silt and clay (0.45–20 μm). As soil P increased, the P enrichment ratio of the SS declined exponentially, and the values of P saturation (Psat; 15–42%) and equilibrium P concentration (EPCo; 0.7–5.5 mg litre−1) in the SS fell within narrower ranges compared with the soils (6–74% and 0.1–10 mg litre−1, respectively). When OP was < 100 mg kg−1, Psat and EPCo values in the SS were smaller than those in the soil and vice-versa, suggesting that eroding particles from soils with both average and high P fertility would release P on entering the local (Rosemaund) stream. Increasing soil OP from average to high P fertility increased the P content of the SS by approximately 10%, but had no significant (P > 0.05) effect on the Psat, or EPCo, of the SS. Management options to reduce soil P status as a means of reducing P losses in land runoff and minimizing eutrophication risk may therefore have more limited effect than is currently assumed in catchment management.

Global climate and ocean circulation on an aquaplanet ocean-atmosphere general circulation model

A low resolution coupled ocean-atmosphere general circulation model OAGCM is used to study the characteristics of the large scale ocean circulation and its climatic impacts in a series of global coupled aquaplanet experiments. Three configurations, designed to produce fundamentally different ocean circulation regimes, are considered. The first has no obstruction to zonal flow, the second contains a low barrier that blocks zonal flow in the ocean at all latitudes, creating a single enclosed basin, whilst the third contains a gap in the barrier to allow circumglobal flow at high southern latitudes. Warm greenhouse climates with a global average air surface temperature of around 27C result in all cases. Equator to pole temperature gradients are shallower than that of a current climate simulation. Whilst changes in the land configuration cause regional changes in temperature, winds and rainfall, heat transports within the system are little affected. Inhibition of all ocean transport on the aquaplanet leads to a reduction in global mean surface temperature of 8C, along with a sharpening of the meridional temperature gradient. This results from a reduction in global atmospheric water vapour content and an increase in tropical albedo, both of which act to reduce global surface temperatures. Fitting a simple radiative model to the atmospheric characteristics of the OAGCM solutions suggests that a simpler atmosphere model, with radiative parameters chosen a priori based on the changing surface configuration, would have produced qualitatively different results. This implies that studies with reduced complexity atmospheres need to be guided by more complex OAGCM results on a case by case basis.

Simulating current global river runoff with a global hydrological model: model revisions, validation, and sensitivity analysis

Global hydrological models (GHMs) model the land surface hydrologic dynamics of continental-scale river basins. Here we describe one such GHM, the Macro-scale - Probability-Distributed Moisture model.09 (Mac-PDM.09). The model has undergone a number of revisions since it was last applied in the hydrological literature. This paper serves to provide a detailed description of the latest version of the model. The main revisions include the following: (1) the ability for the model to be run for n repetitions, which provides more robust estimates of extreme hydrological behaviour, (2) the ability of the model to use a gridded field of coefficient of variation (CV) of daily rainfall for the stochastic disaggregation of monthly precipitation to daily precipitation, and (3) the model can now be forced with daily input climate data as well as monthly input climate data. We demonstrate the effects that each of these three revisions has on simulated runoff relative to before the revisions were applied. Importantly, we show that when Mac-PDM.09 is forced with monthly input data, it results in a negative runoff bias relative to when daily forcings are applied, for regions of the globe where the day-to-day variability in relative humidity is high. The runoff bias can be up to - 80% for a small selection of catchments but the absolute magnitude of the bias may be small. As such, we recommend future applications of Mac-PDM.09 that use monthly climate forcings acknowledge the bias as a limitation of the model. The performance of Mac-PDM.09 is evaluated by validating simulated runoff against observed runoff for 50 catchments. We also present a sensitivity analysis that demonstrates that simulated runoff is considerably more sensitive to method of PE calculation than to perturbations in soil moisture and field capacity parameters.

Biogenic nitric oxide emissions from soils: impact on NOx and ozone over West Africa during AMMA (African Monsoon Multidisciplinary Analysis) - observational study

Chemical and meteorological parameters measured on board the Facility for Airborne Atmospheric Measurements (FAAM) BAe 146 Atmospheric Research Aircraft during the African Monsoon Multidisciplinary Analysis (AMMA) campaign are presented to show the impact of NOx emissions from recently wetted soils in West Africa. NO emissions from soils have been previously observed in many geographical areas with different types of soil/vegetation cover during small scale studies and have been inferred at large scales from satellite measurements of NOx. This study is the first dedicated to showing the emissions of NOx at an intermediate scale between local surface sites and continental satellite measurements. The measurements reveal pronounced mesoscale variations in NOx concentrations closely linked to spatial patterns of antecedent rainfall. Fluxes required to maintain the NOx concentrations observed by the BAe-146 in a number of cases studies and for a range of assumed OH concentrations (1×106 to 1×107 molecules cm−3) are calculated to be in the range 8.4 to 36.1 ng N m−2 s−1. These values are comparable to the range of fluxes from 0.5 to 28 ng N m−2 s−1 reported from small scale field studies in a variety of non-nutrient rich tropical and sub-tropical locations reported in the review of Davidson and Kingerlee (1997). The fluxes calculated in the present study have been scaled up to cover the area of the Sahel bounded by 10 to 20 N and 10 E to 20 W giving an estimated emission of 0.03 to 0.30 Tg N from this area for July and August 2006. The observed chemical data also suggest that the NOx emitted from soils is taking part in ozone formation as ozone concentrations exhibit similar fine scale structure to the NOx, with enhancements over the wet soils. Such variability can not be explained on the basis of transport from other areas. Delon et al. (2008) is a companion paper to this one which models the impact of soil NOx emissions on the NOx and ozone concentration over West Africa during AMMA. It employs an artificial neural network to define the emissions of NOx from soils, integrated into a coupled chemistry-dynamics model. The results are compared to the observed data presented in this paper. Here we compare fluxes deduced from the observed data with the model-derived values from Delon et al. (2008).

Factors governing cellular convection in orographic precipitation

The development of shallow cellular convection in warm orographic clouds is investigated through idealized numerical simulations of moist flow over topography using a cloud-resolving numerical model. Buoyant instability, a necessary element for moist convection, is found to be diagnosed most accurately through analysis of the moist Brunt–Väisälä frequency (N_m) rather than the vertical profile of θ_e. In statically unstable orographic clouds (N_m^2) < 0), additional environmental and terrain-related factors are shown to have major effects on the amount of cellularity that occurs in 2D simulations. One of these factors, the basic-state wind shear, may suppress convection in 2D yet allow for longitudinal convective roll circulations in 3D. The presence of convective structures within an orographic cloud substantially enhanced the maximum rainfall rates, precipitation efficiencies, and precipitation accumulations in all simulations.

A model-based assessment of the effects of projected climate change on the water resources of Jordan

This paper is concerned with the quantification of the likely effect of anthropogenic climate change on the water resources of Jordan by the end of the twenty-first century. Specifically, a suite of hydrological models are used in conjunction with modelled outcomes from a regional climate model, HadRM3, and a weather generator to determine how future flows in the upper River Jordan and in the Wadi Faynan may change. The results indicate that groundwater will play an important role in the water security of the country as irrigation demands increase. Given future projections of reduced winter rainfall and increased near-surface air temperatures, the already low groundwater recharge will decrease further. Interestingly, the modelled discharge at the Wadi Faynan indicates that extreme flood flows will increase in magnitude, despite a decrease in the mean annual rainfall. Simulations projected no increase in flood magnitude in the upper River Jordan. Discussion focuses on the utility of the modelling framework, the problems of making quantitative forecasts and the implications of reduced water availability in Jordan.

Statistical downscaling of river flows

An extensive statistical ‘downscaling’ study is done to relate large-scale climate information from a general circulation model (GCM) to local-scale river flows in SW France for 51 gauging stations ranging from nival (snow-dominated) to pluvial (rainfall-dominated) river-systems. This study helps to select the appropriate statistical method at a given spatial and temporal scale to downscale hydrology for future climate change impact assessment of hydrological resources. The four proposed statistical downscaling models use large-scale predictors (derived from climate model outputs or reanalysis data) that characterize precipitation and evaporation processes in the hydrological cycle to estimate summary flow statistics. The four statistical models used are generalized linear (GLM) and additive (GAM) models, aggregated boosted trees (ABT) and multi-layer perceptron neural networks (ANN). These four models were each applied at two different spatial scales, namely at that of a single flow-gauging station (local downscaling) and that of a group of flow-gauging stations having the same hydrological behaviour (regional downscaling). For each statistical model and each spatial resolution, three temporal resolutions were considered, namely the daily mean flows, the summary statistics of fortnightly flows and a daily ‘integrated approach’. The results show that flow sensitivity to atmospheric factors is significantly different between nival and pluvial hydrological systems which are mainly influenced, respectively, by shortwave solar radiations and atmospheric temperature. The non-linear models (i.e. GAM, ABT and ANN) performed better than the linear GLM when simulating fortnightly flow percentiles. The aggregated boosted trees method showed higher and less variable R2 values to downscale the hydrological variability in both nival and pluvial regimes. Based on GCM cnrm-cm3 and scenarios A2 and A1B, future relative changes of fortnightly median flows were projected based on the regional downscaling approach. The results suggest a global decrease of flow in both pluvial and nival regimes, especially in spring, summer and autumn, whatever the considered scenario. The discussion considers the performance of each statistical method for downscaling flow at different spatial and temporal scales as well as the relationship between atmospheric processes and flow variability.

Toward a combined seasonal weather and crop productivity forecasting system: Determination of the working spatial scale

A methodology is presented for the development of a combined seasonal weather and crop productivity forecasting system. The first stage of the methodology is the determination of the spatial scale(s) on which the system could operate; this determination has been made for the case of groundnut production in India. Rainfall is a dominant climatic determinant of groundnut yield in India. The relationship between yield and rainfall has been explored using data from 1966 to 1995. On the all-India scale, seasonal rainfall explains 52% of the variance in yield. On the subdivisional scale, correlations vary between variance r(2) = 0.62 (significance level p < 10(-4)) and a negative correlation with r(2) = 0.1 (p = 0.13). The spatial structure of the relationship between rainfall and groundnut yield has been explored using empirical orthogonal function (EOF) analysis. A coherent, large-scale pattern emerges for both rainfall and yield. On the subdivisional scale (similar to 300 km), the first principal component (PC) of rainfall is correlated well with the first PC of yield (r(2) = 0.53, p < 10(-4)), demonstrating that the large-scale patterns picked out by the EOFs are related. The physical significance of this result is demonstrated. Use of larger averaging areas for the EOF analysis resulted in lower and (over time) less robust correlations. Because of this loss of detail when using larger spatial scales, the subdivisional scale is suggested as an upper limit on the spatial scale for the proposed forecasting system. Further, district-level EOFs of the yield data demonstrate the validity of upscaling these data to the subdivisional scale. Similar patterns have been produced using data on both of these scales, and the first PCs are very highly correlated (r(2) = 0.96). Hence, a working spatial scale has been identified, typical of that used in seasonal weather forecasting, that can form the basis of crop modeling work for the case of groundnut production in India. Last, the change in correlation between yield and seasonal rainfall during the study period has been examined using seasonal totals and monthly EOFs. A further link between yield and subseasonal variability is demonstrated via analysis of dynamical data.

Design and optimisation of a large-area process-based model for annual crops

The formulation of a new process-based crop model, the general large-area model (GLAM) for annual crops is presented. The model has been designed to operate on spatial scales commensurate with those of global and regional climate models. It aims to simulate the impact of climate on crop yield. Procedures for model parameter determination and optimisation are described, and demonstrated for the prediction of groundnut (i.e. peanut; Arachis hypogaea L.) yields across India for the period 1966-1989. Optimal parameters (e.g. extinction coefficient, transpiration efficiency, rate of change of harvest index) were stable over space and time, provided the estimate of the yield technology trend was based on the full 24-year period. The model has two location-specific parameters, the planting date, and the yield gap parameter. The latter varies spatially and is determined by calibration. The optimal value varies slightly when different input data are used. The model was tested using a historical data set on a 2.5degrees x 2.5degrees grid to simulate yields. Three sites are examined in detail-grid cells from Gujarat in the west, Andhra Pradesh towards the south, and Uttar Pradesh in the north. Agreement between observed and modelled yield was variable, with correlation coefficients of 0.74, 0.42 and 0, respectively. Skill was highest where the climate signal was greatest, and correlations were comparable to or greater than correlations with seasonal mean rainfall. Yields from all 35 cells were aggregated to simulate all-India yield. The correlation coefficient between observed and simulated yields was 0.76, and the root mean square error was 8.4% of the mean yield. The model can be easily extended to any annual crop for the investigation of the impacts of climate variability (or change) on crop yield over large areas. (C) 2004 Elsevier B.V. All rights reserved.

Simulation of crop yields using ERA-40: Limits to skill and nonstationarity in weather-yield relationships

Reanalysis data provide an excellent test bed for impacts prediction systems. because they represent an upper limit on the skill of climate models. Indian groundnut (Arachis hypogaea L.) yields have been simulated using the General Large-Area Model (GLAM) for annual crops and the European Centre for Medium-Range Weather Forecasts (ECMWF) 40-yr reanalysis (ERA-40). The ability of ERA-40 to represent the Indian summer monsoon has been examined. The ability of GLAM. when driven with daily ERA-40 data, to model both observed yields and observed relationships between subseasonal weather and yield has been assessed. Mean yields "were simulated well across much of India. Correlations between observed and modeled yields, where these are significant. are comparable to correlations between observed yields and ERA-40 rainfall. Uncertainties due to the input planting window, crop duration, and weather data have been examined. A reduction in the root-mean-square error of simulated yields was achieved by applying bias correction techniques to the precipitation. The stability of the relationship between weather and yield over time has been examined. Weather-yield correlations vary on decadal time scales. and this has direct implications for the accuracy of yield simulations. Analysis of the skewness of both detrended yields and precipitation suggest that nonclimatic factors are partly responsible for this nonstationarity. Evidence from other studies, including data on cereal and pulse yields, indicates that this result is not particular to groundnut yield. The detection and modeling of nonstationary weather-yield relationships emerges from this study as an important part of the process of understanding and predicting the impacts of climate variability and change on crop yields.