The triggering of orographic rainbands by small-scale topography

The triggering of convective orographic rainbands by small-scale topographic features is investigated through observations of a banded precipitation event over the Oregon Coastal Range and simulations using a cloud-resolving numerical model. A quasi-idealized simulation of the observed event reproduces the bands in the radar observations, indicating the model’s ability to capture the physics of the band-formation process. Additional idealized simulations reinforce that the bands are triggered by lee waves past small-scale topographic obstacles just upstream of the nominal leading edge of the orographic cloud. Whether a topographic obstacle in this region is able to trigger a strong rainband depends on the phase of its lee wave at cloud entry. Convective growth only occurs downstream of obstacles that give rise to lee-wave-induced displacements that create positive vertical velocity anomalies w_c and nearly zero buoyancy anomalies b_c as air parcels undergo saturation. This relationship is quantified through a simple analytic condition involving w_c, b_c, and the static stability N_m^2 of the cloud mass. Once convection is triggered, horizontal buoyancy gradients in the cross-flow direction generate circulations that align the bands parallel to the flow direction.

Ground deformation at Soufriere Hills Volcano, Montserrat during 1998-2000 measured by radar interferometry and GPS

We examine the motion of the ground surface on the Soufriere Hills Volcano, Montserrat between 1998 and 2000 using radar interferometry (InSAR). To minimise the effects of variable atmospheric water vapour on the InSAR measurements we use independently-derived measurements of the radar path delay from six continuous GPS receivers. The surfaces providing a measurable inter-ferometric signal are those on pyroclastic flow deposits, mainly emplaced in 1997. Three types of surface motion can be discriminated. Firstly, the surfaces of thick, valley-filling deposits subsided at rates of 150-120 mm/year in the year after emplacement to 50-30 mm/year two years later. This must be due to contraction and settling effects during cooling. The second type is the near-field motion localised within about one kilometre of the dome. Both subsidence and uplift events are seen and though the former could be due to surface gravitational effects, the latter may reflect shallow (< 1 km) pressurisation effects within the conduit/dome. Far-field motions of the surface away from the deeply buried valleys are interpreted as crustal strains. Because the flux of magma to the surface stopped from March 1998 to November 1999 and then resumed from November 1999 through 2000, we use InSAR data from these two periods to test the crustal strain behaviour of three models of magma supply: open, depleting and unbalanced. The InSAR observations of strain gradients of 75-80 mm/year/krn uplift during the period of quiescence on the western side of the volcano are consistent with an unbalanced model in which magma supply into a crustal magma chamber continues during quiescence, raising chamber pressure that is then released upon resumption of effusion. GPS motion vectors agree qualitatively with the InSAR displacements but are of smaller magnitude. The discrepancy may be due to inaccurate compensation for atmospheric delays in the InSAR data. (c) 2005 Elsevier B.V. All rights reserved.

The zonal asymmetry of the Southern Hemisphere storm-track

Atmospheric general circulation model experiments have been performed to investigate how the significant zonal asymmetry in the Southern Hemisphere (SH) winter storm track is forced by sea surface temperature (SST) and orography. An experiment with zonally symmetric tropical SSTs expands the SH upper-tropospheric storm track poleward and eastward and destroys its spiral structure. Diagnosis suggests that these aspects of the observed storm track result from Rossby wave propagation from a wave source in the Indian Ocean region associated with the monsoon there. The lower-tropospheric storm track is not sensitive to this forcing. However, an experiment with zonally symmetric midlatitude SSTs exhibits a marked reduction in the magnitude of the maximum intensity of the lower-tropospheric storm track associated with reduced SST gradients in the western Indian Ocean. Experiments without the elevation of the South African Plateau or the Andes show reductions in the intensity of the major storm track downstream of them due to reduced cyclogenesis associated with the topography. These results suggest that the zonal asymmetry of the SH winter storm track is mainly established by stationary waves excited by zonal asymmetry in tropical SST in the upper troposphere and by local SST gradients in the lower troposphere, and that it is modified through cyclogenesis associated with the topography of South Africa and South America.

The cumulus-capped boundary layer. I: Modelling transports in the cloud layer

Scalar-flux budgets have been obtained from large-eddy simulations (LESs) of the cumulus-capped boundary layer. Parametrizations of the terms in the budgets are discussed, and two parametrizations for the transport term in the cloud layer are proposed. It is shown that these lead to two models for scalar transports by shallow cumulus convection. One is equivalent to the subsidence detrainment form of convective tendencies obtained from mass-flux parametrizations of cumulus convection. The second is a flux-gradient relationship that is similar in form to the non-local parametrizations of turbulent transports in the dry-convective boundary layer. Using the fluxes of liquid-water potential temperature and total water content from the LES, it is shown that both models are reasonable diagnostic relations between fluxes and the vertical gradients of the mean fields. The LESs used in this study are for steady-state convection and it is possible to treat the fluxes of conserved thermodynamic variables as independent, and ignore the effects of condensation. It is argued that a parametrization of cumulus transports in a model of the cumulus-capped boundary layer should also include an explicit representation of condensation. A simple parametrization of the liquid-water flux in terms of conserved variables is also derived.

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.

Ocean circulation and climate in an idealised Pangean OAGCM

An idealised Pangean configuration is integrated in a coupled ocean atmosphere general circulation model to investigate the form of the ocean circulation and its impacts on the large scale climate system. A vigorous, hemispherically symmetric overturning is found, driven by deep water formation at high latitudes. Whilst the peak mass transport is around 100Sv, a low vertical temperature gradient in the ocean means that the maximum heat transport is only 1.2PW. The geographical change in the coupled model is found to produce a global average warming of 2°C, despite an increase in global surface albedo. This occurs through changes in the atmospheric water vapour and cloud distributions. There is also reduction in the equator-pole temperature gradient, largely attributable to the same causes, avoiding the paradox of low meridional temperature gradients without increased polar heat transport.

The impact of global warming on the tropical Pacific Ocean and El Niño

The El Niño–Southern Oscillation (ENSO) is a naturally occurring fluctuation that originates in the tropical Pacific region and affects ecosystems, agriculture, freshwater supplies, hurricanes and other severe weather events worldwide. Under the influence of global warming, the mean climate of the Pacific region will probably undergo significant changes. The tropical easterly trade winds are expected to weaken; surface ocean temperatures are expected to warm fastest near the equator and more slowly farther away; the equatorial thermocline that marks the transition between the wind-mixed upper ocean and deeper layers is expected to shoal; and the temperature gradients across the thermocline are expected to become steeper. Year-to-year ENSO variability is controlled by a delicate balance of amplifying and damping feedbacks, and one or more of the physical processes that are responsible for determining the characteristics of ENSO will probably be modified by climate change. Therefore, despite considerable progress in our understanding of the impact of climate change on many of the processes that contribute to El Niño variability, it is not yet possible to say whether ENSO activity will be enhanced or damped, or if the frequency of events will change.

The role of lateral ocean physics in the upper ocean thermal balance of a coupled ocean-atmosphere GCM

The sensitivity of the upper ocean thermal balance of an ocean-atmosphere coupled GCM to lateral ocean physics is assessed. Three 40-year simulations are performed using horizontal mixing, isopycnal mixing, and isopycnal mixing plus eddy induced advection. The thermal adjustment of the coupled system is quite different between the simulations, confirming the major role of ocean mixing on the heat balance of climate. The initial adjustment phase of the upper ocean (SST) is used to diagnose the physical mechanisms involved in each parametrisation. When the lateral ocean physics is modified, significant changes of SST are seen, mainly in the southern ocean. A heat budget of the annual mixed layer (defined as the “bowl”) shows that these changes are due to a modified heat transfer between the bowl and the ocean interior. This modified heat intake of the ocean interior is directly due to the modified lateral ocean physics. In isopycnal diffusion, this heat exchange, especially marked at mid-latitudes, is both due to an increased effective surface of diffusion and to the sign of the isopycnal gradients of temperature at the base of the bowl. As this gradient is proportional to the isopycnal gradient of salinity, this confirms the strong role of salinity in the thermal balance of the coupled system. The eddy induced advection also leads to increased exchanges between the bowl and the ocean interior. This is both due to the shape of the bowl and again to the existence of a salinity structure. The lateral ocean physics is shown to be a significant contributor to the exchanges between the diabatic and the adiabatic parts of the ocean.

Impacts of a pesticide on pollinator species richness at different spatial scales

Pesticides are an important potential cause of biodiversity and pollinator decline. Little is known about the impacts of pesticides on wild pollinators in the field. Insect pollinators were sampled in an agricultural system in Italy with the aim of detecting the impacts of pesticide use. The insecticide fenitrothion was over 150 times greater in toxicity than other pesticides used in the area, so sampling was set up around its application. Species richness of wild bees, bumblebees and butterflies were sampled at three spatial scales to assess responses to pesticide application: (i) the ‘field’ scale along pesticide drift gradients; (ii) the ‘landscape’ scale sampling in different crops within the area and (iii) the ‘regional’ scale comparing two river basins with contrasting agricultural intensity. At the field scale, the interaction between the application regime of the insecticide and the point in the season was important for species richness. Wild bee species richness appeared to be unaffected by one insecticide application, but declined after two and three applications. At the landscape scale, the species richness of wild bees declined in vine fields where the insecticide was applied, but did not decline in maize or uncultivated fields. At the regional scale, lower bumblebee and butterfly species richness was found in the more intensively farmed basin with higher pesticide loads. Our results suggest that wild bees are an insect pollinator group at particular risk from pesticide use. Further investigation is needed on how the type, quantity and timing of pesticide application impacts pollinators.

Evaluation of the S(T) assimilation method with the Argo dataset

Recent observations from the Argo dataset of temperature and salinity profiles are used to evaluate a series of 3-year data assimilation experiments in a global ice–ocean general circulation model. The experiments are designed to evaluate a new data assimilation system whereby salinity is assimilated along isotherms, S(T ). In addition, the role of a balancing salinity increment to maintain water mass properties is investigated. This balancing increment is found to effectively prevent spurious mixing in tropical regions induced by univariate temperature assimilation, allowing the correction of isotherm geometries without adversely influencing temperature–salinity relationships. In addition, the balancing increment is able to correct a fresh bias associated with a weak subtropical gyre in the North Atlantic using only temperature observations. The S(T ) assimilation method is found to provide an important improvement over conventional depth level assimilation, with lower root-mean-squared forecast errors over the upper 500 m in the tropical Atlantic and Pacific Oceans. An additional set of experiments is performed whereby Argo data are withheld and used for independent evaluation. The most significant improvements from Argo assimilation are found in less well-observed regions (Indian, South Atlantic and South Pacific Oceans). When Argo salinity data are assimilated in addition to temperature, improvements to modelled temperature fields are obtained due to corrections to model density gradients and the resulting circulation. It is found that observations from the Argo array provide an invaluable tool for both correcting modelled water mass properties through data assimilation and for evaluating the assimilation methods themselves.

Impact of hydrographic data assimilation on the modelled Atlantic meridional overturning circulation

Here we make an initial step toward the development of an ocean assimilation system that can constrain the modelled Atlantic Meridional Overturning Circulation (AMOC) to support climate predictions. A detailed comparison is presented of 1° and 1/4° resolution global model simulations with and without sequential data assimilation, to the observations and transport estimates from the RAPID mooring array across 26.5° N in the Atlantic. Comparisons of modelled water properties with the observations from the merged RAPID boundary arrays demonstrate the ability of in situ data assimilation to accurately constrain the east-west density gradient between these mooring arrays. However, the presence of an unconstrained "western boundary wedge" between Abaco Island and the RAPID mooring site WB2 (16 km offshore) leads to the intensification of an erroneous southwards flow in this region when in situ data are assimilated. The result is an overly intense southward upper mid-ocean transport (0–1100 m) as compared to the estimates derived from the RAPID array. Correction of upper layer zonal density gradients is found to compensate mostly for a weak subtropical gyre circulation in the free model run (i.e. with no assimilation). Despite the important changes to the density structure and transports in the upper layer imposed by the assimilation, very little change is found in the amplitude and sub-seasonal variability of the AMOC. This shows that assimilation of upper layer density information projects mainly on the gyre circulation with little effect on the AMOC at 26° N due to the absence of corrections to density gradients below 2000 m (the maximum depth of Argo). The sensitivity to initial conditions was explored through two additional experiments using a climatological initial condition. These experiments showed that the weak bias in gyre intensity in the control simulation (without data assimilation) develops over a period of about 6 months, but does so independently from the overturning, with no change to the AMOC. However, differences in the properties and volume transport of North Atlantic Deep Water (NADW) persisted throughout the 3 year simulations resulting in a difference of 3 Sv in AMOC intensity. The persistence of these dense water anomalies and their influence on the AMOC is promising for the development of decadal forecasting capabilities. The results suggest that the deeper waters must be accurately reproduced in order to constrain the AMOC.

Land degradation and vegetation distribution in Chott El Beida wetland, Algeria

The aim of this study is to explore the environmental factors that determine plant Community distribution in northeast Algeria. This paper provides a quantitative analysis of the vegetation-environment relationships for a study site in the Cholt El Beida wetland, a RAMSAR site in Setif, Algeria. Sixty vegetation plots were sampled and analysed using TWINSPAN and Detrended Correspondence Analysis (DCA) in order to identify the principal vegetation communities and determine the environmental gradients associated with these. 127 species belonging to 41 families and 114 genera were recorded. Six of the recorded species were endemic representing 4.7% of the total species. The richest families were Compositae, Gramineae, Cruciferae and Chenopodiaceae. Therophytes and hemicryptophytes were the most frequent life forms. the Mediterranean floristic element is dominant and is represented by 39 species. The samples were classified into four main community types. The principal DCA axes represent gradients of soil salinity, moisture and anthropogenic pressure. The use of classification in combination with ordination techniques resulted in a good discrimination between plant communities and a greater understanding of controlling environmental factors. The methodology adopted can be employed for improving baseline information on plant community ecology and distribution in often critically endangered Mediterranean wetland areas. (C) 2008 Elsevier Ltd. All rights reserved.

Taxonomical vs. functional responses of bee communities to fire in two contrasting climatic regions

Valuable insights into mechanisms of community responses to environmental change can be gained by analysing in tandem the variation in functional and taxonomic composition along environmental gradients. We assess the changes in species and functional trait composition (i.e. dominant traits and functional diversity) of diverse bee communities in contrasting fire-driven systems in two climatic regions: Mediterranean (scrub habitats in Israel) and temperate (chestnut forests in southern Switzerland). In both climatic regions, there were shifts in species diversity and composition related to post-fire age. In the temperate region, functional composition responded markedly to fire; however, in the Mediterranean, the taxonomic response to fire was not matched by functional replacement. These results suggest that greater functional stability to fire in the Mediterranean is achieved by replacement of functionally similar species (i.e. functional redundancy) which dominate under different environmental conditions in the heterogeneous landscapes of the region. In contrast, the greater functional response in the temperate region was attributed to a more rapid post-fire vegetation recovery and shorter time-window when favourable habitat was available relative to the Mediterranean. Bee traits can be used to predict the functional responses of bee communities to environmental changes in habitats of conservation importance in different regions with distinct disturbance regimes. However, predictions cannot be generalized from one climatic region to another where distinct habitat configurations occur.

Phylogeny and the hierarchical organization of plant diversity

R. H. Whittaker's idea that plant diversity can be divided into a hierarchy of spatial components from alpha at the within-habitat scale through beta for the turnover of species between habitats to gamma along regional gradients implies the underlying existence of alpha, beta, and gamma niches. We explore the hypothesis that the evolution of a, (3, and gamma niches is also hierarchical, with traits that define the a niche being labile, while those defining a and 7 niches are conservative. At the a level we find support for the hypothesis in the lack of close significant phylogenetic relationship between meadow species that have similar a niches. In a second test, a niche overlap based on a variety of traits is compared between congeners and noncongeners in several communities; here, too, there is no evidence of a correlation between a niche and phylogeny. To test whether beta and gamma niches evolve conservatively, we reconstructed the evolution of relevant traits on evolutionary trees for 14 different clades. Tests against null models revealed a number of instances, including some in island radiations, in which habitat (beta niche) and elevational maximum (an aspect of the gamma niche) showed evolutionary conservatism.

Macrophyte community structure and species occurrence in relation to environmental determinants in the ephemeral aquatic habitats of Gavdos, Greece

The aims of this study were to explore the environmental factors that determine the distribution of plant communities in temporary rock pools and provide a quantitative analysis of vegetation-environment relationships for five study sites on the island of Gavdos, southwest of Crete, Greece. Data from 99 rock pools were collected and analysed using Two-Way Indicator Species Analysis (TWINSPAN), Detrended Correspondence Analysis (DCA) and Canonical Correspondence Analysis (CCA) to identify the principal communities and environmental gradients that are linked to community distribution. A total of 46 species belonging to 21 families were recorded within the study area. The dominant families were Labiatae, Gramineae and Compositae while therophytes and chamaephytes were the most frequent life forms. The samples were classified into six community types using TWINSPAN, which were also corroborated by CCA analysis. The principal gradients for vegetation distribution, identified by CCA, were associated with water storage and water retention ability, as expressed by pool perimeter and water depth. Generalised Additive Models (GAMs) were employed to identify responses of four dominant rock pool species to water depth. The resulting species response curves showed niche differentiation in the cases of Callitriche pulchra and Tillaea vaillantii and revealed competition between Zannichellia pedunculata and Chara vulgaris. The use of classification in combination with ordination techniques resulted in a good discrimination between plant communities. Generalised Additive Models are a powerful tool in investigating species response curves to environmental gradients. The methodology adopted can be employed for improving baseline information on plant community ecology and distribution in Mediterranean ephemeral pools.