Diabatic balance model for the equatorial atmosphere

Using an asymptotic expansion, a balance model is derived for the shallow-water equations (SWE) on the equatorial beta-plane that is valid for planetary-scale equatorial dynamics and includes Kelvin waves. In contrast to many theories of tropical dynamics, neither a strict balance between diabatic heating and vertical motion nor a small Froude number is required. Instead, the expansion is based on the smallness of the ratio of meridional to zonal length scales, which can also be interpreted as a separation in time scale. The leading-order model is characterized by a semigeostrophic balance between the zonal wind and meridional pressure gradient, while the meridional wind v vanishes; the model is thus asymptotically nondivergent, and the nonzero correction to v can be found at the next order. Importantly for applications, the diagnostic balance relations are linear for winds when inferring the wind field from mass observations and the winds can be diagnosed without direct observations of diabatic heating. The accuracy of the model is investigated through a set of numerical examples. These examples show that the diagnostic balance relations can remain valid even when the dynamics do not, and the balance dynamics can capture the slow behavior of a rapidly varying solution.

Bollène-2002 experiment: radar quantitative precipitation estimation in the Cévennes–Vivarais region, France

The Bollène-2002 Experiment was aimed at developing the use of a radar volume-scanning strategy for conducting radar rainfall estimations in the mountainous regions of France. A developmental radar processing system, called Traitements Régionalisés et Adaptatifs de Données Radar pour l’Hydrologie (Regionalized and Adaptive Radar Data Processing for Hydrological Applications), has been built and several algorithms were specifically produced as part of this project. These algorithms include 1) a clutter identification technique based on the pulse-to-pulse variability of reflectivity Z for noncoherent radar, 2) a coupled procedure for determining a rain partition between convective and widespread rainfall R and the associated normalized vertical profiles of reflectivity, and 3) a method for calculating reflectivity at ground level from reflectivities measured aloft. Several radar processing strategies, including nonadaptive, time-adaptive, and space–time-adaptive variants, have been implemented to assess the performance of these new algorithms. Reference rainfall data were derived from a careful analysis of rain gauge datasets furnished by the Cévennes–Vivarais Mediterranean Hydrometeorological Observatory. The assessment criteria for five intense and long-lasting Mediterranean rain events have proven that good quantitative precipitation estimates can be obtained from radar data alone within 100-km range by using well-sited, well-maintained radar systems and sophisticated, physically based data-processing systems. The basic requirements entail performing accurate electronic calibration and stability verification, determining the radar detection domain, achieving efficient clutter elimination, and capturing the vertical structure(s) of reflectivity for the target event. Radar performance was shown to depend on type of rainfall, with better results obtained with deep convective rain systems (Nash coefficients of roughly 0.90 for point radar–rain gauge comparisons at the event time step), as opposed to shallow convective and frontal rain systems (Nash coefficients in the 0.6–0.8 range). In comparison with time-adaptive strategies, the space–time-adaptive strategy yields a very significant reduction in the radar–rain gauge bias while the level of scatter remains basically unchanged. Because the Z–R relationships have not been optimized in this study, results are attributed to an improved processing of spatial variations in the vertical profile of reflectivity. The two main recommendations for future work consist of adapting the rain separation method for radar network operations and documenting Z–R relationships conditional on rainfall type.

Timing and depositional environments of a Middle Pleistocene glaciation of northeast England: New evidence from Warren House Gill, County Durham

At various times during the Quaternary, north-eastern England was a zone of confluence between dynamic ice lobes sourced from the Pennines, northern Scotland, the Cheviots, and Scandinavia. The region thus has some of the most complex exposures of Middle to Late Pleistocene sediments in Britain, with both interglacial and glacial sediments deposited in terrestrial and marine settings. We investigated sedimentary sequences exposed on the coastline of County Durham at Warren House Gill, and present a new model of British and Fennoscandian Ice Sheet interaction in the North Sea Basin during the Middle Pleistocene. The stratigraphy at Warren House Gill consists of a lower diamicton and upper estuarine sediments, both part of the Warren House Formation. They are separated from the overlying Weichselian Blackhall and Horden tills by a substantial unconformity. The lower diamicton of the Warren House Formation is re-interpreted here as an MIS 8 to 12 glaciomarine deposit containing ice-rafted lithics from north-eastern Scotland and the northeast North Sea, and is renamed the ‘Ash Gill Member’. It is dated by lithological comparison to the Easington Raised Beach, Middle Pleistocene Amino Acid Racemisation values, and indirectly by optically stimulated luminescence. The overlying shallow subaqueous sediments were deposited in an estuarine environment by suspension settling and bottom current activity. They are named the ‘Whitesides Member’, and form the uppermost member of the Warren House Formation. During glaciation, ice-rafted material was deposited in a marine embayment. There is no evidence of a grounded, onshore Scandinavian ice sheet in County Durham during MIS 6, which has long been held as the accepted stratigraphy. This has major implications for the currently accepted British Quaternary Stratigraphy. Combined with recent work on the Middle Pleistocene North Sea Drift from Norfolk, which is now suggested to have been deposited by a Scottish ice sheet, the presence of a Scandinavian ice sheet in eastern England at any time during the Quaternary is becoming increasingly doubtful.

Forests of the tropical eastern Andean flank during the middle pleistocene

Inter-bedded volcanic and organic sediments from Erazo (Ecuador) indicate the presence of four different forest assemblages on the eastern Andean flank during the middle Pleistocene. Radiometric dates (40Ar–39Ar) obtained fromthe volcanic ash indicate that deposition occurred between 620,000 and 192,000 years ago. Examination of the organic sediment composition and the fossil pollen, wood and charcoal it contains provides insight into depositional environment, vegetation assemblage and fire history. The high organic content and abundance of macro fossils found throughout the sediment suggest that during the period of deposition the local environment was either a swamp or a shallow water body. The correlation of fire activity (peaks in charcoal abundance) with volcanic ash deposits through most of the record suggests that volcanoes were the main source of ignition. The low abundance of grass (typically b10%) throughout the sedimentary sequence along with the low abundance of other taxa indicative of open vegetation suggests the persistence of forest at Erazo. Four types of forest assemblage were identified (with the first taxa as the most dominant): i) Alnus-Arecaceae, ii) Miconia- Melastomataceae/Combretaceae-Moraceae/Urticaceae, iii) Arecaceae-Alnus, and iv) Podocarpus with Oreopanax sp. and Melastomataceae/Combretaceae. Changes in the forest floristic composition indicate high vegetation turnover and reassortment of taxa between upper and lower montane forests during the middle Pleistocene as well as the persistence of forest cover.

Deglacial changes of the southern margin of the southern westerly winds revealed by terrestrial records from SW Patagonia (52 degrees S)

Much of the ongoing discussion regarding synchrony or bipolar asynchrony of paleoclimate events has centered on the timing and structure of the last glacial termination in the southern mid- latitudes, in particular the southwestern Patagonian region (50�e55�S). Its location adjacent to the Drake Passage andnear the southern margin of the southern westerly winds (SWW) allows examining the postulated links between the Southern Oceane SWW coupled system and tmospheric CO2 variations through the last glacial termination. Results from two sites located in the Última Esperanza area (52�S) allow us to infer SWW-driven changes in hydrologic balance during this critical time interval. These findings indicate peatland development under temperate/wet conditions between 14,600 and 14,900 cal yr BP, followed by cooling and a lake transgressive phase that led to a shallow lake during the early part of the Antarctic Cold Reversal (ACR, 13,600-14,600 cal yr BP), followed in turn by a deeper lake and modest warming during Younger Dryas time (YD, 11,800-13,000 cal yr BP), superseded by terrestrialization and forest expansion at the beginning of the Holocene. We propose that the SWW (i) strengthened and shifted northward during ACR time causing a precipitation rise in northwestern and southwestern Patagonia coeval with mid- and high-latitude cooling and a halt in the deglacial atmospheric CO2 rise; (ii) shifted southward during YD time causing a precipitation decline/increase in NW/SW Patagonia, respectively, high-latitude warming, and invigorated CO2 release from the Southern Ocean; (iii) became weaker between 10,000 and 11,500 cal yr BP causing a precipitation decline throughout Patagonia, concurrent with peak mid- and high-latitude temperatures and atmospheric CO2 concentrations.

Foehn jets over the Larsen C ice shelf, Antarctica

Previously unknown foehn jets have been identified to the east of the Antarctic Peninsula (AP) above the Larsen C Ice Shelf. These jets have major implications for the east coast of the AP, a region of rapid climatic warming and where two large sections of ice shelf have collapsed in recent years. During three foehn events across the AP, leeside warming and drying is seen in new aircraft observations and simulated well by the Met Office Unified Model (MetUM) at ∼1.5 km grid spacing. In case A, weak southwesterly flow and an elevated upwind inversion characterise a highly nonlinear flow regime with upwind flow blocking. In case C strong northwesterly winds characterise a relatively linear case with little upwind flow blocking. Case B resides somewhere between the two in flow regime linearity. The foehn jets – apparent in aircraft observations where available and MetUM simulations of all three cases – are mesoscale features (up to 60 km in width) originating from the mouths of leeside inlets. Through back trajectory analysis they are identified as a type of gap flow. In cases A and B the jets are distinct, being strongly accelerated relative to the background flow, and confined to low levels above the Larsen C Ice Shelf. They resemble the ‘shallow foehn’ of the Alps. Case C resembles a case of ‘deep foehn’, with the jets less distinct. The foehn jets are considerably cooler and moister relative to adjacent regions of calmer foehn air. This is due to a dampened foehn effect in the jet regions: in case A the jets have lower upwind source regions, and in the more linear case C there is less diabatic warming and precipitation along jet trajectories due to the reduced orographic uplift across the mountain passes.

Hydrology and climatology at Laguna La Gaiba, lowland Bolivia: complex responses to climatic forcings over the last 25,000 years

Diatom, geochemical and isotopic data provide a record of environmental change in Laguna La Gaiba, lowland Bolivia (17°450S, 57°350W), over the last ca. 25 000 years. High-resolution diatom analysis around the Last Glacial–Interglacial Transition provides new insights into this period of change. The full and late glacial lake was generally quite shallow, but with evidence of periodic flooding. At about 13 100 cal a BP, just before the start of the Younger Dryas chronozone, the diatoms indicate shallower water conditions, but there is a marked change at about 12 200 cal a BP indicating the onset of a period of high variability, with rising water levels punctuated by periodic drying. From ca. 11 800 to 10 000 cal a BP, stable, deeper water conditions persisted. There is evidence for drying in the early to middle Holocene, but not as pronounced as that reported from elsewhere in the southern hemisphere tropics of South America. This was followed by the onset of wetter conditions in the late Holocene consistent with insolation forcing. Conditions very similar to present were established about 2100 cal a BP. A complex response to both insolation forcing and millennial-scale events originating in the North Atlantic is noted.

The development of a data-driven application benchmarking approach to performance modelling

Performance modelling is a useful tool in the lifeycle of high performance scientific software, such as weather and climate models, especially as a means of ensuring efficient use of available computing resources. In particular, sufficiently accurate performance prediction could reduce the effort and experimental computer time required when porting and optimising a climate model to a new machine. In this paper, traditional techniques are used to predict the computation time of a simple shallow water model which is illustrative of the computation (and communication) involved in climate models. These models are compared with real execution data gathered on AMD Opteron-based systems, including several phases of the U.K. academic community HPC resource, HECToR. Some success is had in relating source code to achieved performance for the K10 series of Opterons, but the method is found to be inadequate for the next-generation Interlagos processor. The experience leads to the investigation of a data-driven application benchmarking approach to performance modelling. Results for an early version of the approach are presented using the shallow model as an example.

Geophysical characterization of the complex dynamics of groundwater and seawater exchange in a highly stressed aquifer system linked to a coastal lagoon (SE Spain)

Anthropogenic pressure influences the two-way interactions between shallow aquifers and coastal lagoons. Aquifer overexploitation may lead to seawater intrusion, and aquifer recharge from rainfall plus irrigation may, in turn, increase the groundwater discharge into the lagoon. We analyse the evolution, since the 1950s up to the present, of the interactions between the Campo de Cartagena Quaternary aquifer and the Mar Menor coastal lagoon (SE Spain). This is a very heterogeneous and anisotropic detrital aquifer, where aquifer–lagoon interface has a very irregular geometry. Using electrical resistivity tomography, we clearly identified the freshwater–saltwater transition zone and detected areas affected by seawater intrusion. Severity of the intrusion was spatially variable and significantly related to the density of irrigation wells in 1950s–1960s, suggesting the role of groundwater overexploitation. We distinguish two different mechanisms by which water from the sea invades the land: (a) horizontal advance of the interface due to a wide exploitation area and (b) vertical rise (upconing) caused by local intensive pumping. In general, shallow parts of the geophysical profiles show higher electrical resistivity associated with freshwater mainly coming from irrigation return flows, with water resources mostly from deep confined aquifers and imported from Tagus river, 400 km north. This indicates a likely reversal of the former seawater intrusion process.

Impact of atmospheric transport on the evolution of microphysical and optical properties of Saharan dust

Saharan dust affects the climate by altering the radiation balance and by depositing minerals to the Atlantic Ocean. Both are dependent on particle size. We present aircraft measurements comprising 42 profiles of size distribution (0.1–300 µm), representing freshly uplifted dust, regional aged dust, and dust in the Saharan Air Layer (SAL) over the Canary Islands. The mean effective diameter of dust in SAL profiles is 4.5 µm smaller than that in freshly uplifted dust, while the vertical structure changes from a low shallow layer (0–1.5 km) to a well-mixed deep Saharan dust layer (0–5 km). Size distributions show a loss of 60 to 90% of particles larger than 30 µm 12 h after uplift. The single scattering albedo (SSA) increases from 0.92 to 0.94 to 0.95 between fresh, aged, and SAL profiles: this is enough to alter heating rates by 26%. Some fresh dust close to the surface shows SSA as low as 0.85

Investigation of discrepancies in satellite rainfall estimates over Ethiopia

Tropical Applications of Meteorology Using Satellite and Ground-Based Observations (TAMSAT) rainfall estimates are used extensively across Africa for operational rainfall monitoring and food security applications; thus, regional evaluations of TAMSAT are essential to ensure its reliability. This study assesses the performance of TAMSAT rainfall estimates, along with the African Rainfall Climatology (ARC), version 2; the Tropical Rainfall Measuring Mission (TRMM) 3B42 product; and the Climate Prediction Center morphing technique (CMORPH), against a dense rain gauge network over a mountainous region of Ethiopia. Overall, TAMSAT exhibits good skill in detecting rainy events but underestimates rainfall amount, while ARC underestimates both rainfall amount and rainy event frequency. Meanwhile, TRMM consistently performs best in detecting rainy events and capturing the mean rainfall and seasonal variability, while CMORPH tends to overdetect rainy events. Moreover, the mean difference in daily rainfall between the products and rain gauges shows increasing underestimation with increasing elevation. However, the distribution in satellite–gauge differences demon- strates that although 75% of retrievals underestimate rainfall, up to 25% overestimate rainfall over all eleva- tions. Case studies using high-resolution simulations suggest underestimation in the satellite algorithms is likely due to shallow convection with warm cloud-top temperatures in addition to beam-filling effects in microwave- based retrievals from localized convective cells. The overestimation by IR-based algorithms is attributed to nonraining cirrus with cold cloud-top temperatures. These results stress the importance of understanding re- gional precipitation systems causing uncertainties in satellite rainfall estimates with a view toward using this knowledge to improve rainfall algorithms.

On the cause of a magnetospheric flux transfer event

We present a detailed investigation of a magnetospheric flux transfer event (FTE) seen by the Active Magnetospheric Tracer Explorer (AMPTE) UKS and IRM satellites around 1046 UT on October 28, 1984. This event has been discussed many times previously in the literature and has been cited as support for a variety of theories of FTE formation. We make use of a model developed to reproduce ion precipitations seen in the cusp ionosphere. The analysis confirms that the FTE is well explained as a brief excursion into an open low-latitude boundary layer (LLBL), as predicted by two theories of magnetospheric FTEs: namely, that they are bulges in the open LLBL due to reconnection rate enhancements or that they are indentations of the magnetopause by magnetosheath pressure increases (but in the presence of ongoing steady reconnection). The indentation of the inner edge of the open LLBL that these two models seek to explain is found to be shallow for this event. The ion model reproduces the continuous evolution of the ion distribution function between the sheath-like population at the event center and the surrounding magnetospheric populations; it also provides an explanation of the high-pressure core of the event as comprising field lines that were reconnected considerably earlier than those that are draped over it to give the event boundary layer. The magnetopause transition parameter is used to isolate a field rotation on the boundaries of the core, which is subjected to the tangential stress balance test. The test identifies this to be a convecting structure, which is neither a rotational discontinuity (RD) nor a contact discontinuity, but could possibly be a slow shock. In addition, evidence for ion reflection off a weak RD on the magnetospheric side of this structure is found. The event structure is consistent in many ways with features predicted for the open LLBL by analytic MHD theories and by MHD and hybrid simulations. The de Hoffman-Teller velocity of the structure is significantly different from that of the magnetosheath flow, indicating that it is not an indentation caused by a high-pressure pulse in the sheath but is consistent with the motion of newly opened field lines (different from the sheath flow because of the magnetic tension force) deduced from the best fit to the ion data. However, we cannot here rule out the possibility that the sheath flow pattern has changed in the long interval between the two satellites observing the FTE and subsequently emerging into the magnetosheath; thus this test is not conclusive in this particular case. Analysis of the fitted elapsed time since reconnection shows that the core of the event was reconnected in one pulse and the event boundary layer was reconnected in a subsequent pulse. Between these two pulses is a period of very low (but nonzero) reconnection rate, which lasts about 14 mins. Thus the analysis supports, but does not definitively verify, the concept that the FTE is a partial passage into an open LLBL caused by a traveling bulge in that layer produced by a pulse in reconnection rate.

A benchmark-driven modelling approach for evaluating deployment choices on a multi-core architecture

The complexity of current and emerging architectures provides users with options about how best to use the available resources, but makes predicting performance challenging. In this work a benchmark-driven model is developed for a simple shallow water code on a Cray XE6 system, to explore how deployment choices such as domain decomposition and core affinity affect performance. The resource sharing present in modern multi-core architectures adds various levels of heterogeneity to the system. Shared resources often includes cache, memory, network controllers and in some cases floating point units (as in the AMD Bulldozer), which mean that the access time depends on the mapping of application tasks, and the core's location within the system. Heterogeneity further increases with the use of hardware-accelerators such as GPUs and the Intel Xeon Phi, where many specialist cores are attached to general-purpose cores. This trend for shared resources and non-uniform cores is expected to continue into the exascale era. The complexity of these systems means that various runtime scenarios are possible, and it has been found that under-populating nodes, altering the domain decomposition and non-standard task to core mappings can dramatically alter performance. To find this out, however, is often a process of trial and error. To better inform this process, a performance model was developed for a simple regular grid-based kernel code, shallow. The code comprises two distinct types of work, loop-based array updates and nearest-neighbour halo-exchanges. Separate performance models were developed for each part, both based on a similar methodology. Application specific benchmarks were run to measure performance for different problem sizes under different execution scenarios. These results were then fed into a performance model that derives resource usage for a given deployment scenario, with interpolation between results as necessary.

The hidden organic carbon in deep mineral soils

Aims Current estimates of soil organic carbon (SOC) are based largely on surficial measurements to depths of 0.3 to 1 m. Many of the world’s soils greatly exceed 1 m depth and there are numerous reports of biological activity to depths of many metres. Although SOC storage to depths of up to 8 m has been previously reported, the extent to which SOC is stored at deeper depths in soil profiles is currently unknown. This paper aims to provide the first detailed analysis of these previously unreported stores of SOC. Methods Soils from five sites in the deeply weathered regolith in the Yilgarn Craton of south-western Australia were sampled and analysed for total organic carbon by combustion chromatography. These soils ranged between 5 and 38 m (mean 21 m) depth to bedrock and had been either recently reforested with Pinus pinaster or were under agriculture. Sites had a mean annual rainfall of between 399 and 583 mm yr−1. Results The mean SOC concentration across all sites was 2.30 ± 0.26 % (s.e.), 0.41 ± 0.05 % and 0.23 ± 0.04 % in the surface 0.1, 0.1–0.5 and 0.5 to 1.0 m increments, respectively. The mean value between 1 and 5 m was 0.12 ± 0.01 %, whereas between 5 and 35 m the values decreased from 0.04 ± 0.002 % to 0.03 ± 0.003 %. Mean SOC mass densities for each of the five locations varied from 21.8–37.5 kg C m−2, and were in toto two to five times greater than would be reported with sampling to a depth of 0.5 m. Conclusions This finding may have major implications for estimates of global carbon storage and modelling of the potential global impacts of climate change and land-use change on carbon cycles. The paper demonstrates the need for a reassessment of the current arbitrary shallow soil sampling depths for assessing carbon stocks, a revision of global SOC estimates and elucidation of the composition and fate of deep carbon in response to land use and climate change

A sustainable agricultural landscape for Australia: A review of interlacing carbon sequestration, biodiversity and salinity management in agroforestry systems

Transformation of the south-western Australian landscape from deep-rooted woody vegetation systems to shallow-rooted annual cropping systems has resulted in the severe loss of biodiversity and this loss has been exacerbated by rising ground waters that have mobilised stored salts causing extensive dry land salinity. Since the original plant communities were mostly perennial and deep rooted, the model for sustainable agriculture and landscape water management invariably includes deep rooted trees. Commercial forestry is however only economical in higher rainfall (>700 mm yr−1) areas whereas much of the area where biodiversity is threatened has lower rainfall (300–700 mm yr−1). Agroforestry may provide the opportunity to develop new agricultural landscapes that interlace ecosystem services such as carbon mitigation via carbon sequestration and biofuels, biodiversity restoration, watershed management while maintaining food production. Active markets are developing for some of these ecosystem services, however a lack of predictive metrics and the regulatory environment are impeding the adoption of several ecosystem services. Nonetheless, a clear opportunity exists for four major issues – the maintenance of food and fibre production, salinisation, biodiversity decline and climate change mitigation – to be managed at a meaningful scale and a new, sustainable agricultural landscape to be developed.