Cloud computing and virtualization within the regional climate model and evaluation system

The Regional Climate Model Evaluation System (RCMES) facilitates the rapid, flexible inclusion of NASA observations into climate model evaluations. RCMES provides two fundamental components. A database (RCMED) is a scalable point-oriented cloud database used to elastically store remote sensing observations and to make them available using a space time query interface. The analysis toolkit (RCMET) is a Python-based toolkit that can be delivered as a cloud virtual machine, or as an installer package deployed using Python Buildout to users in order to allow for temporal and spatial regridding, metrics calculation (RMSE, bias, PDFs, etc.) and end-user visualization. RCMET is available to users in an “offline”, lone scientist mode based on a virtual machine dynamically constructed with model outputs and observations to evaluate; or on an institution’s computational cluster seated close to the observations and model outputs. We have leveraged RCMES within the content of the Coordinated Regional Downscaling Experiment (CORDEX) project, working with the University of Cape Town and other institutions to compare the model output to NASA remote sensing data; in addition we are also working with the North American Regional Climate Change Assessment Program (NARCCAP). In this paper we explain the contribution of cloud computing to RCMES’s specifically describing studies of various cloud databases we evaluated for RCMED, and virtualization toolkits for RCMET, and their potential strengths in delivering user-created dynamic regional climate model evaluation virtual machines for our users.     

Post-peak,fluid-mediated modification of granulite facies zircon and monazite in the Trivandrum Block,southern India

The quarry at Kottavattom in the Trivandrum Block of southern India contains spectacular examples of fluid-assisted alteration of high-grade metamorphic rocks. Garnet-biotite gneiss has undergone a change in mineral assemblage to form submetre scale orthopyroxene-bearing patches, later retrogressed to form an amphibole-bearing lithology. These patches, often referred to as arrested or incipient charnockite, crosscut the original metamorphic foliation and are typically attributed to passage of a low aH₂O fluid through the rock. Whilst this conversion is recognised as a late stage process, little detailed chronological work exists to link it temporally to metamorphism in the region. Zircon and monazite analysed from Kottavattom not only record metamorphism in the Trivandrum Block but also show internal, lobate textures crosscutting the original zoning, consistent with fluid-aided coupled dissolution-reprecipitation during formation of the orthopyroxene-bearing patches. High-grade metamorphism at the quarry occurred between the formation of metamorphic monazite at ~585 Ma and the growth of metamorphic zircon at ~523 Ma. The fluid-assisted alteration of the garnet-biotite gneiss is poorly recorded by altered zircon with only minimal resetting of the U–Pb system, whereas monazite has in some cases undergone complete U–Pb resetting and records an age for fluid infiltration at ~495 Ma. The fluid event therefore places the formation of the altered patches at least 25 Myr after the zircon crystallisation in the garnet-biotite gneiss. The most likely fluid composition causing the modification and U–Pb resetting of zircon and monazite is locally derived hypersaline brine.     

Significance of ancient sulfide PGE and Re–Os signatures in the mantle beneath Calatrava,Central Spain

Spinel lherzolite and wehrlite xenoliths from the Cenozoic Calatrava volcanic field carry the geochemical imprint of metasomatic agents that have affected the subcontinental lithospheric mantle beneath Central Iberia. Some xenoliths (mainly wehrlites) were enriched in REE, Sr, P, and CO₂ by silicic-carbonate-rich metasomatic melts/fluids, while others record the effects of subduction-related hydrous silicate fluids that have precipitated amphibole and induced high Ti/Eu in primary clinopyroxene. The petrographic observations and geochemical data suggest that interstitial glass in the xenoliths represent the quenched products of Si-rich melts that infiltrated the mantle peridotite shortly before the entrainment of the xenoliths in the host magmas that erupted ca 2 million years ago. During their infiltration, the metasomatic melts reacted with peridotite, resulting in silica enrichment, while remobilizing grains of iron-rich monosulfide solid solution (Fe-rich Mss) initially enclosed in, or intergranular to, primary olivine and pyroxenes. In situ laser ablation inductively coupled plasma-mass spectrometry analysis of single sulfide grains reveals that the Fe-rich Mss in glass shows platinum-group element (PGE) patterns and ¹⁸⁷Os/¹⁸⁸Os compositions identical to the Fe-rich Mss occurring as inclusions in, or at grain boundaries of primary silicates. Moreover, independent of its microstructural position, Fe-rich Mss exhibits PGE and ¹⁸⁷Os/¹⁸⁸Os signatures typical of Mss either residual after partial melting or crystallized directly from sulfide melts. Our findings reveal that young metasomatic melt(s)/fluid(s) may carry remobilized sulfides with PGE and Os-isotopic signatures identical to those of texturally older sulfides in the peridotite xenolith. These sulfides thus still provide useful information about the timing and nature of older magmatic events in the subcontinental mantle.     

Emplacement ages and sources of kimberlites and related rocks in southern Africa: U–Pb ages and Sr–Nd isotopes of groundmass perovskite

Groundmass perovskite has been dated by LA-ICPMS in 135 kimberlites and related rocks from 110 localities across southern Africa. Sr and/or Nd isotopes have been analysed by LA-MC-ICPMS in a subset of these and integrated with published data. The age distribution shows peaks at 1,600–1,800, 1,000–1,200, 500–800 and 50–130 Ma. The major “bloom” of Group I kimberlites at ca 90 ± 10 Ma was preceded by a slow build-up in magmatic activity from ca 180 Ma. The main pulse of Group II kimberlites at 120–130 Ma was a distinct episode within this build-up. Comparison of the isotopic data with seismic tomography images suggests that metasomatized subcontinental lithospheric mantle (SCLM) with very low ε _Nd and high ⁸⁷Sr/⁸⁶Sr, (the isotopic signature of Group II kimberlites) was focused in low-Vs zones along translithospheric structures. Such metasomatized zones existed as early as 1,800 Ma, but were only sporadically tapped until the magmatic build-up began at ca 180 Ma, and contributed little to the kimberlitic magmas after ca 110 Ma. We suggest that these metasomatized volumes resided in the deep SCLM and that their low-melting point components were “burned off” by rising temperatures, presumably during an asthenospheric upwelling that led to SCLM thinning and a rise in the ambient geotherm between 120 and 90 Ma. The younger Group I kimberlites therefore rarely interacted with such SCLM, but had improved access to shallower volumes of differently metasomatized, ancient SCLM with low ⁸⁷Sr/⁸⁶Sr and intermediate ε _Nd (0–5). The kimberlite compositions therefore reflect the evolution of the SCLM of southern Africa, with metasomatic-enrichment events from as early as 1.8 Ga, through a major thermal and compositional change at ca 110 Ma, and the major kimberlite “bloom” around 90 Ma.     

Bayesian modelling the retreat of the Irish Sea Ice Stream

We present an 8000‐year history spanning 650 km of ice margin retreat for the largest marine‐terminating ice stream draining the former British–Irish Ice Sheet. Bayesian modelling of the geochronological data shows the ISIS expanded 34.0–25.3 ka, accelerating into the Celtic Sea to reach maximum limits 25.3–24.5 ka before a collapse with rapid marginal retreat to the northern Irish Sea Basin (ISB). This retreat was rapid and driven by climatic warming, sea‐level rise, mega‐tidal amplitudes and reactivation of meridional circulation in the North Atlantic. The retreat, though rapid, is uneven, with the stepped retreat pattern possibly a function of the passage of the ice stream between normal and adverse ice bed gradients and changing ice stream geometry. Initially, wide calving margins and adverse slopes encouraged rapid retreat (～550 m a^?1) that slowed (～100 m a^?1) at the topographic constriction and bathymetric high between southern Ireland and Wales before rates increased (～200 m a^?1) across adverse bed slopes and wider and deeper basin configuration in the northern ISB. These data point to the importance of the ice bed slope and lateral extent in predicting the longer‐term (>1000 a) patterns and rates of ice‐marginal retreat during phases of rapid collapse, which has implications for the modelling of projected rapid retreat of present‐day ice streams. Copyright © 2013 John Wiley & Sons, Ltd.     

Implementing quality control on a random number stream to improve a stochastic weather generator

For decades, stochastic modellers have used computerized random number generators to produce random numeric sequences fitting a specified statistical distribution. Unfortunately, none of the random number generators we tested satisfactorily produced the target distribution. The result is generated distributions whose mean even diverges from the mean used to generate them, regardless of the length of run. Non‐uniform distributions from short sequences of random numbers are a major problem in stochastic climate generation, because truly uniform distributions are required to produce the intended climate parameter distributions. In order to ensure generation of a representative climate with the stochastic weather generator CLIGEN within a 30‐year run, we tested the climate output resulting from various random number generators. The resulting distributions of climate parameters showed significant departures from the target distributions in all cases. We traced this failure back to the uniform random number generators themselves. This paper proposes a quality control approach to select only those numbers that conform to the expected distribution being retained for subsequent use. The approach is based on goodness‐of‐fit analysis applied to the random numbers generated. Normally distributed deviates are further tested with confidence interval tests on their means and standard deviations. The positive effect of the new approach on the climate characteristics generated and the subsequent deterministic process‐based hydrology and soil erosion modelling are illustrated for four climatologically diverse sites. Copyright © 2007 John Wiley & Sons, Ltd.     

Controls on surface water chemistry in two lake‐watersheds in the Adirondack region of New York: differences in nitrogen solute sources and sinks

The southwestern Adirondack region of New York receives among the highest rates of atmospheric nitrogen (N) deposition in the USA. Atmospheric N deposition to sensitive ecosystems, like the Adirondacks, may increase the acidification of soils through losses of exchangeable nutrient cations, and the acidification of surface waters associated with enhanced mobility of nitrate (NO₃^?). However, watershed attributes, including surficial terrestrial characteristics, in‐lake processing, and geological settings, have been found to complicate the relationships between atmospheric N deposition and N drainage losses. We studied two lake‐watersheds in the southwestern Adirondacks, Grass Pond and Constable Pond, which are located in close proximity (～26 km) and receive similarly high N deposition, but have contrasting watershed attributes (e.g. wetland area, geological settings). Since the difference in the influence of N deposition was minimal, we were able to examine both within‐ and between‐watershed influences of land cover, the contribution of glacial till groundwater inputs, and in‐lake processes on surface water chemistry with particular emphasis on N solutes and dissolved organic carbon (DOC). Monthly samples at seven inlets and one outlet of each lake were collected from May to October in 1999 and 2000. The concentrations of NO₃^? were high at the Grass Pond inlets, especially at two inlets, and NO₃^? was the major N solute at the Grass Pond inlets. The concentrations of likely weathering products (i.e. dissolved Si, Ca²⁺, Mg²⁺, Na⁺) as well as acid neutralizing capacity and pH values, were also particularly high at those two Grass Pond inlets, suggesting a large contribution of groundwater inputs. Dissolved organic N (DON) was the major N solute at the Constable Pond inlets. The higher concentrations of DON and DOC at the Constable Pond inlets were attributed to a large wetland area in the watershed. The DOC/DON ratios were also higher at the Constable Pond inlets, possibly due to a larger proportion of coniferous forest area. Although DON and DOC were strongly related, the stronger relationship of the proportion of wetland area with DOC suggests that additional factors regulate DON. The aggregated representation of watershed physical features (i.e. elevation, watershed area, mean topographic index, hypsometric‐analysis index) was not clearly related to the lake N and DOC chemistry. Despite distinctive differences in inlet N chemistry, NO₃^? and DON concentrations at the outlets of the two lakes were similar. The lower DOC/DON ratios at the lake outlets and at the inlets having upstream ponds suggest the importance of N processing and organic N sources within the lakes. Although an inverse relationship between NO₃^? and DOC/DON has been suggested to be indicative of a N deposition gradient, the existence of this relationship for sites that receive similar atmospheric N deposition suggest that the relationship between NO₃^? and the DOC/DON ratio is derived from environmental and physical factors. Our results suggest that, despite similar wet N deposition at the two watershed sites, N solutes entering lakes were strongly affected by hydrology associated with groundwater contribution and the presence of wetlands, whereas N solutes leaving lakes were strongly influenced by in‐lake processing. Copyright © 2006 John Wiley & Sons, Ltd.