Validating the Flash Code: Vortex-Dominated Flows

As a component of the Flash Center’s validation program, we compare FLASH simulation results with experimental results from Los Alamos National Laboratory. The flow of interest involves the lateral interaction between a planar M _a = 1.2 shock wave with a cylinder of gaseous sulfur hexafluoride (SF₆) in air, and in particular the development of primary and secondary instabilities after the passage of the shock. While the overall evolution of the flow is comparable in the simulations and experiments, small-scale features are difficult to match. We focus on the sensitivity of numerical results to simulation parameters.     

Petrogenesis of extension-related alkaline volcanism in Karaburhan (Sivrihisar–Eskisehir), NW Anatolia, Turkey

Alkaline lavas were erupted as phonolites and trachytes around Karaburhan (Sivrihisar–Eskisehir, NW Anatolia) within the Izmir–Ankara–Erzincan suture zone. These volcanic rocks were emplaced as domes, close and parallel to the ophiolite thrust line. According to ⁴⁰Ar/³⁹Ar geochronological analyses of sanidine crystals from the phonolites, the age of the alkaline volcanics is 25 Ma (Late Oligocene–Early Miocene).The flow-textured phonolites are porphyritic and consist mainly of sanidine, clinopyroxene, and feldspathoid crystals. The clinopyroxenes show compositional zoning, with aegirine (Na_0.82–0.96Fe⁺³_0.68–0.83) rims and aegirine–augite cores (containing calcium, magnesium, and Fe⁺²). Some aegirine–augites are replaced with sodium-, calcium-, and magnesium-rich amphibole (hastingsite). Feldspathoid (hauyne) crystals enriched with elemental Na and Ca have been almost completely altered to zeolite and carbonate minerals. The fine-grained trachytes with a trachytic texture consist of feldspar (oligoclase and sanidine) phenocrystals and clinopyroxene microphenocrystals within a groundmass made up largely of alkali feldspar microlites.Although there are some differences in their element patterns, the phonolites and trachytes exhibit enrichment in LILEs (Sr, K, Rb, Ba, Th) and LREEs (La, Ce, Pr, Nd) and negative anomalies in Nb and Ta. These geochemical characteristics indicate a lithospheric mantle enriched by fluids extracted from the subduction component. In addition, the high ⁸⁷Sr/⁸⁶Sr (0.706358–0.708052) and low ¹⁴³Nd/¹⁴⁴Nd (0.512546–0.512646) isotope concentrations of the alkaline lavas reflect a mantle source that has undergone metasomatism by subduction-derived fluids. Petrogenetic modeling indicates that the alkaline lavas generated from the subduction-modified lithospheric mantle have undergone assimilation, fractional crystallization, and crustal contamination, acquiring high Pb, Ba, Rb, and Sr contents and Pb isotopic compositions during their ascent through the thickened crust in an extensional setting.     

Ca isotope cycling in a forested ecosystem

Reports of large Ca isotope fractionations between trees and soils prompted this study of a Boreal forest ecosystem near La Ronge, Saskatchewan, to improve understanding of this phenomenon. The results on five tree species (black spruce, trembling aspen, white spruce, jack pine, balsam poplar) confirm that nutrient Ca uptake by plants favors the light isotopes, thus driving residual Ca in plant available soil pools towards enrichment in the heavy isotopes. Substantial within-tree fraction occurs in tissues formed along the transpiration stream, with low δ⁴⁴Ca values in fine roots (2 mm), intermediate values in stemwood, and high values in foliage. Separation factors between different plant tissues are similar between species, but the initial fractionation step in the tips of the fine roots is species specific, and/or sensitive to the local soil environment. Soil water δ⁴⁴Ca values appear to increase with depth to at least 35 cm below the top of the forest floor, which is close to the deepest level of fine roots. The heavy plant fractionated signature of Ca in the finely rooted upper soils filters downward where it is retained on ion exchange sites, leached into groundwater, and discharged into surface waters.The relationship between Ca uptake by tree fine roots and the pattern of δ⁴⁴Ca enrichment with soil depth was modeled for two Ca pools: the forest floor (litter) and the underlying (upper B) mineral soil. Six study plots were investigated along two hillside toposequences trending upwards from a first order stream. We used allometric equations describing the Ca distribution in boreal tree species to calculate weighted average δ⁴⁴Ca values for the stands in each plot and estimate Ca uptake rates. The δ⁴⁴Ca value of precipitation was measured, and soil weathering signatures deduced, by acid leaching of lower B mineral soils. Steady state equations were used to derive a set of model Ca fluxes and fractionation factors for each plot. The model reproduces the increase in δ⁴⁴Ca with depth found in forest floor and upper B soil waters. Transient model runs show that the forest Ca cycle is sensitive to changes in plant Ca uptake rate, such as would occur during ontogeny or disturbance. Accordingly, secular records of δ⁴⁴Ca in tree ring cellulose have the potential to monitor changes in the forest Ca cycle through time, thus providing a new tool for evaluating natural and anthropogenic impacts on forest health. Another model run shows that by changing the size of the isotope fractionation factor and adjusting for differences in forest productivity, that the range in Ca isotope fractionation in forested ecosystems reported in the literature, thus far, is reproduced. As a quantitative tool, the Ca cycling model produces a reasonable set of relative Ca fluxes for the La Ronge site, consistent with Environment Canada’s measurements for wet deposition in the region and simulated Ca release from soil mineral weathering using the PROFILE model. But the sensitivity of the model is limited by the small range of fractionation observed in this boreal shield setting of ∼1‰, which limits accuracy. If the model were applied to a site with a greater range in δ⁴⁴Ca values among the principal Ca fluxes, it is capable of producing robust and reliable estimations of Ca fluxes that are otherwise difficult to measure in forested ecosystems.     

A topical evaluation and discussion of data movement technologies for data-intensive scientific applications

Transferring large volumes of information from one location to potentially many others that are geographically distributed and across varying networks is still prevalent in modern scientific data systems. This is despite the movement to push computation to the data and to reduce data movement needed to compute answers to challenging scientific problems, to disseminate information to the scientific community, and to acquire data for curation and enrichment. Because of this, it is imperative that decisions made regarding data movement systems and architectures be backed by both analytical rigor, and also by empirical evidence and measurement. The purpose of this study is to expand on the work performed by our research team over the last decade and to take a fresh look at the evaluation of multiple topical data transfer technologies in use cases derived from data-intensive scientific systems and applications in the areas of Earth science. We report on the evaluation of a set of data movement technologies against a set of empirically derived comparison dimensions. Based on this evaluation, we make recommendations towards the selection of appropriate data movement technologies in scientific applications and scenarios.     

Platinum-group elements (Ir,Pd, Pt and Rh) in road-deposited sediments in two urban watersheds,Hawaii

This study represents the first detailed investigation of platinum-group elements (PGEs) in road-deposited sediment (RDS) in Hawaii, USA. Thirty-three sample locations, in two urban watersheds in Honolulu, Oahu, Hawaii were sampled. The <63 μm fraction of RDS was digested with aqua regia, followed by matrix separation with Dowex AG50-X8 cation exchange resin. PGEs were analyzed with inductively coupled plasma-mass spectrometry (ICP-MS) equipped with a desolvating nebulizer. Concentrations of Rh, Pd and Pt in residential streets reached 64, 105 and 506 ng/g, respectively. Maximum enrichment ratios, computed as RDS concentrations relative to baseline values, exceeded 400, indicating a significant anthropogenic signal with the sequence Rh > Pt > Pd. Iridium concentrations were uniformly low <1 ng/g, and enrichment ratios support a geogenic source. Significant interelement PGE correlations (Pd–Pt–Rh), combined with the magnitude of PGE pair-wise ratios (Pt/Pd, Pt/Rh and Pd/Rh), and relative percentages comparable to European RDS and roadside soil in Indiana, USA all suggest an automobile source. Attrition of PGE-loaded automobile catalytic converters and subsequent loss to the environment by exhaust emissions explains the significant environmental signal of PGEs in road environments of Hawaii. Further PGE work is required to quantify urban transport paths as PGEs are known to bioaccumulate, cause cellular damage and may have detrimental human health effects.     

An extreme wave event in eastern Yucatán,Mexico: Evidence of a palaeotsunami event during the Mayan times

The Yucatán Peninsula, Mexico, has typically been considered a tectonically stable region with little significant seismic activity. The region though, is one that is regularly affected by hurricanes. A detailed survey of ca 100 km of the eastern Yucatán and Cozumel coast identified the presence of ridges containing individual boulders measuring >1 m in length. The boulder ridges reach 5 m in height and their origin is associated with extreme wave event activity. Previously modelled tsunami waves from known seismically active zones in the region (Muertos Thrust Belt and South Caribbean Deformed Belt) are not of sufficient scale in the area of the Yucatán Peninsula to have produced the boulder ridges recorded in this study. The occurrence of hurricanes in this region is more common, but two of the most destructive (Hurricane Gilbert 1988 and Hurricane Wilma 2005) produced coastal waves too small to have created the ridges recorded here. In this paper, a new tsunami model with a source area located on the Motagua/Swan Island Fault System has been generated that indicates a tsunami event may have caused the extreme wave events that resulted in the deposition of the boulder ridges.     

Recent progress on combining geomorphological and geochronological data with ice sheet modelling,demonstrated using the last British–Irish Ice Sheet

Palaeo-ice sheets are important analogues for understanding contemporary ice sheets, offering a record of ice sheet behaviour that spans millennia. There are two main approaches to reconstructing palaeo-ice sheets. Empirical reconstructions use the available glacial geological and chronological evidence to estimate ice sheet extent and dynamics but lack direct consideration of ice physics. In contrast, numerically modelled simulations implement ice physics, but often lack direct quantitative comparison with empirical evidence. Despite being long identified as a fruitful scientific endeavour, few ice sheet reconstructions attempt to reconcile the empirical and model-based approaches. To achieve this goal, model-data comparison procedures are required. Here, we compare three numerically modelled simulations of the former British–Irish Ice Sheet with the following lines of evidence: (a) position and shape of former margin positions, recorded by moraines; (b) former ice-flow direction and flow-switching, recorded by flowsets of subglacial bedforms; and (c) the timing of ice-free conditions, recorded by geochronological data. These model–data comparisons provide a useful framework for quantifying the degree of fit between numerical model simulations and empirical constraints. Such tools are vital for reconciling numerical modelling and empirical evidence, the combination of which will lead to more robust palaeo-ice sheet reconstructions with greater explicative and ultimately predictive power.