Modifying the Jackson index to quantify the relationship between geology,landscape structure,and water transit time in steep wet headwaters

The relationship between stream water mean transit time (MTT), catchment geology, and landscape structure is still poorly characterized. Here, we present a new simple index that builds on the Jackson, Bitew, and Du (2014) index that focuses specifically on permeability contrasts at the soil–bedrock interface and digital elevation model-based physical flow path measurements to identify broad landscape trends of moisture redistribution in the subsurface of steep wet headwater catchments. We use this index to explore the relationship between geology, landscape structure, and water transit time through the lens of landscape anisotropy. We hypothesize that catchments with a greater tendency to shed water laterally will correlate with younger stream water MTT and catchments with a greater tendency to infiltrate water vertically will correlate with older stream water MTT. We tested the new index at eight geologically diverse Pacific Rim catchments in Oregon, Japan, and New Zealand. The new index explained 77% of the variability in measured stream water MTT across these varied sites. These findings suggest that critical zone anisotropy and catchment form are first-order controls on the time scales over which catchments store and release their water and that a simple index may usefully capture this relationship.     

Spatial variation of crop yield response to climate change in East Africa

There is general consensus that the impacts of climate change on agriculture will add significantly to the development challenges of ensuring food security and reducing poverty, particularly in Africa. While these changes will influence agriculture at a broad scale, regional or country-level assessments can miss critical detail. We use high-resolution methods to generate characteristic daily weather data for a combination of different future emission scenarios and climate models to drive detailed simulation models of the maize and bean crops. For the East African region, there is considerable spatial and temporal variation in this crop response. We evaluate the response of maize and beans to a changing climate, as a prelude to detailed targeting of options that can help smallholder households adapt. The results argue strongly against the idea of large, spatially contiguous development domains for identifying and implementing adaptation options, particularly in regions with large variations in topography and current average temperatures. Rather, they underline the importance of localised, community-based efforts to increase local adaptive capacity, take advantage of changes that may lead to increased crop and livestock productivity where this is possible, and to buffer the situations where increased stresses are likely.     

Increasing impacts of climate change upon ecosystems with increasing global mean temperature rise

In a meta-analysis we integrate peer-reviewed studies that provide quantified estimates of future projected ecosystem changes related to quantified projected local or global climate changes. In an advance on previous analyses, we reference all studies to a common pre-industrial base-line for temperature, employing up-scaling techniques where necessary, detailing how impacts have been projected on every continent, in the oceans, and for the globe, for a wide range of ecosystem types and taxa. Dramatic and substantive projected increases of climate change impacts upon ecosystems are revealed with increasing annual global mean temperature rise above the pre-industrial mean (ΔT_g). Substantial negative impacts are commonly projected as ΔT_g reaches and exceeds 2°C, especially in biodiversity hotspots. Compliance with the ultimate objective of the United Nations Framework Convention on Climate Change (Article 2) requires that greenhouse gas concentrations be stabilized within a time frame “sufficient to allow ecosystems to adapt naturally to climate change”. Unless ΔT_g is constrained to below 2°C at most, results here imply that it will be difficult to achieve compliance. This underscores the need to limit greenhouse gas emissions by accelerating mitigation efforts and by protecting existing ecosystems from greenhouse-gas producing land use change processes such as deforestation.     

Teleconnections in a warmer climate: the pliocene perspective

Migrations toward altered sea surface temperature (SST) patterns in the Indo-Pacific region are present in the recent observational record and in future global warming projections. These SSTs are in the form of ??permanent?? El Ni?o-like (herein termed ??El Padre??) and Indian Ocean Dipole (IOD)-like patterns. The Early Pliocene Warm Period, which bears similarity to future warming projections, may have also exhibited these Indo-Pacific SST patterns, as suggested by regional terrestrial paleo-climatic data and general circulation model studies. The ability to corroborate this assessment with paleo-data reconstructions is an advantage of the warm Pliocene period that is not afforded by future warming scenarios. Thus, the Pliocene period provides us with a warm-climate perspective and test bed for understanding potential changes to future atmospheric interactions given these altered SST states. This study specifically assesses how atmospheric teleconnections from El Padre/IOD SST patterns are generated and propagate to create the regional climate signals of the Pliocene period, as these signals may be representative of future regional climatic changes as well. To do this, we construct a holistic diagnostic rubric that allows us to examine atmospheric teleconnections, both energetically and dynamically, as produced by a general circulation model. We incorporate KE??, a diagnostic adapted from the eddy kinetic energy generation field, to assess the available energy transferred to these teleconnections. Using this methodology, we found that relative to our Modern Control experiments, weaker atmospheric teleconnections prevail under warm Pliocene conditions, although pathways of propagation still appear directed toward the southwestern United States from our tropical Pacific sector forcing. Propagation directly emanating from the Indian Ocean forcing sector appears to be largely blocked, although indirect teleconnective pathways appear traversing the Asian continent toward the North Pacific. The changes in the atmospheric circulation of Indian Ocean region in response to the underlying specified SST forcing (and indicated by Pliocene paleo-data) may have a host of implications for energy transfer out of and into the region, including interactions with the Asian jet stream and changes to the seasonal monsoon cycle. These interactions warrant further study in both past and future warm climate scenarios.     

K–Ar ages of meteorites: Clues to parent-body thermal histories

Whereas most radiometric chronometers give formation ages of individual meteorites >4.5 Ga ago, the K–Ar chronometer rarely gives times of meteorite formation. Instead, K–Ar ages obtained by the ³⁹Ar–⁴⁰Ar technique span the entire age of the solar system and typically measure the diverse thermal histories of meteorites or their parent objects, as produced by internal parent body metamorphism or impact heating. This paper briefly explains the Ar–Ar dating technique. It then reviews Ar–Ar ages of several different types of meteorites, representing at least 16 different parent bodies, and discusses the likely thermal histories these ages represent. Ar–Ar ages of ordinary (H, L, and LL) chondrites, R chondrites, and enstatite meteorites yield cooling times following internal parent body metamorphism extending over ∼200 Ma after parent body formation, consistent with parent bodies of ∼100 km diameter. For a suite of H-chondrites, Ar–Ar and U–Pb ages anti-correlate with the degree of metamorphism, consistent with increasing metamorphic temperatures and longer cooling times at greater depths within the parent body. In contrast, acapulcoites–lodranites, although metamorphosed to higher temperatures than chondrites, give Ar–Ar ages which cluster tightly at ∼4.51 Ga. Ar–Ar ages of silicate from IAB iron meteorites give a continual distribution across ∼4.53–4.32 Ga, whereas silicate from IIE iron meteorites give Ar–Ar ages of either ∼4.5 Ga or ∼3.7 Ga. Both of these parent bodies suffered early, intense collisional heating and mixing. Comparison of Ar–Ar and I–Xe ages for silicate from three other iron meteorites also suggests very early collisional heating and mixing. Most mesosiderites show Ar–Ar ages of ∼3.9 Ga, and their significantly sloped age spectra and Ar diffusion properties, as well as Ni diffusion profiles in metal, indicate very deep burial after collisional mixing and cooling at a very slow rate of ∼0.2 °C/Ma. Ar–Ar ages of a large number of brecciated eucrites range over ∼3.4–4.1 Ga, similar to ages of many lunar highland rocks. These ages on both bodies were reset by large impact heating events, possibly initiated by movements of the giant planets. Many impact-heated chondrites show impact-reset Ar–Ar ages of either >3.5 Ga or <1.0 Ga, and generally only chondrites show these younger ages. The younger ages may represent orbital evolution times in the asteroid belt prior to ejection into Earth-crossing orbits. Among martian meteorites, Ar–Ar ages of nakhlites are similar to ages obtained from other radiometric chronometers, but apparent Ar–Ar ages of younger shergottites are almost always older than igneous crystallization ages, because of the presence of excess (parentless) ⁴⁰Ar. This excess ⁴⁰Ar derives from shock-implanted martian atmosphere or from radiogenic ⁴⁰Ar inherited from the melt. Differences between meteorite ages obtained from other chronometers (e.g., I–Xe and U–Pb) and the oldest measured Ar–Ar ages are consistent with previous suggestions that the ⁴⁰K decay parameters in common use are incorrect and that the K–Ar age of a 4500 Ma meteorite should be possibly increased, but by no more than ∼20 Ma.     

The transition to microbial photosynthesis in hot spring ecosystems

Even casual observations of continental hot springs reveal that photosynthesis has its limits. In an effort to explore the transition to photosynthesis, field measurements of temperature and pH were made at 996 hot spring locations at Yellowstone National Park ranging from 14° to 94 °C and pH from 0.8 to 9.7. In addition, sulfide measurements were made in 426 of these locations showing concentrations up to 8820 μg L^? 1 total sulfide. These data indicate that the previously established upper temperature (73–75 °C) for the transition to photosynthesis is reached in many basic hot springs, but that the transition occurs at lower temperature with decreasing pH below ~ 6.5. As an example, no strong evidence for photosynthesis was found above 45 °C at pH ~ 2. In several locations, photosynthesis appears to be suppressed despite temperatures and pH values that permit photosynthesis elsewhere. Sulfide concentrations may be responsible for the suppression of photosynthesis at these sites. Total sulfide concentrations were observed to decrease downstream in hot spring outflow channels. Abiotic processes (degassing, oxidation, mineral precipitation, etc.) are too slow to account for these decreases, suggesting an explanation from microbial sulfide oxidation that is supported by field experiments. Microbial sulfide oxidation may determine the ultimate suitability of some hot springs for microbial photosynthesis.     

Modeling locally varying anisotropy of CO2 emissions in the United States

In a spatial property modeling context, the variables of interest to be modeled often display complex nonlinear features. Techniques to incorporate these nonlinear features, such as multiple point statistics or cummulants, are often complex with input parameters that are difficult to infer. The methodology proposed in this paper uses a classical vector-based definition of locally varying anisotropy to characterize nonlinear features and incorporate locally varying anisotropy into numerical property models. The required input is an exhaustive field of anisotropy orientation and magnitude. The methodology consists of (1) using the shortest path distance between locations to define the covariance between points in space (2) multidimensional scaling of the domain to ensure positive definite kriging equations and (3) estimation or simulation with kriging or sequential Gaussian simulation. The only additional parameter required when kriging or simulating with locally varying anisotropy is the number of dimensions to retain in multidimensional scaling. The methodology is demonstrated on a CO₂ emissions data set for the United States in 2002 and shows an improvement in cross validation results as well as a visual reproduction of nonlinear features.     

The Hf-Nd isotopic composition of marine sediments

This paper presents new major and trace-element data and Lu-Hf and Sm-Nd isotopic compositions for representative suites of marine sediment samples from 14 drill sites outboard of the world’s major subduction zones. These suites and samples were chosen to represent the global range in lithology, Lu/Hf ratios, and sediment flux in subducting sediments worldwide. The data reported here represent the most comprehensive data set on subducting sediments and define the Hf-Nd isotopic variations that occur in oceanic sediments and constrain the processes that caused them.Using new marine sediment data presented here, in conjunction with published data, we derive a new Terrestrial Array given by the equation, ε_Hf = 1.55 × ε_Nd + 1.21. This array was calculated using >3400 present-day Hf and Nd isotope values. The steeper slope and smaller y-intercept of this array, compared to the original expression (ε_Hf = 1.36 × ε_Nd + 2.89; Vervoort et al., 1999) reflects the use of present day values and the unradiogenic Hf of old continental samples included in the array.In order to examine the Hf-Nd isotopic variations in marine sediments, we have classified our samples into 5 groups based on lithology and major and trace-element geochemical compositions: turbidites, terrigenous clays, and volcaniclastic, hydrothermal and hydrogenetic sediments. Compositions along the Terrestrial Array are largely controlled by terrigenous material derived from the continents and delivered to the ocean basins via turbidites, volcaniclastic sediments, and volcanic inputs from magmatic arcs. Compositions below the Terrestrial Array derive from unradiogenic Hf in zircon-rich turbidites. The anomalous compositions above the Terrestrial Array largely reflect the decoupled behavior of Hf and Nd during continental weathering and delivery to the ocean. Both terrigenous and hydrogenetic clays possess anomalously radiogenic Hf, reflecting terrestrial sedimentary and weathering processes on the one hand and marine inheritance on the other. This probably occurs during complementary processes involving preferential retention of unradiogenic Hf on the continents in the form of zircon and release of radiogenic Hf from the breakdown of easily weathered, high Lu-Hf phases such as apatite.     

Optical dating of tsunami-laid sand from an Oregon coastal lake

Optical ages for five samples of tsunami-laid sand from an Oregon coastal lake were determined using an infrared optical-dating method on K-feldspar separates and, as a test of accuracy, compared to ages determined by AMS ¹⁴C dating of detrital plant fragments found in the same beds. Two optical ages were about 20% younger than calibrated ¹⁴C ages of about 3.1 and 4.3 ka. Correction of the optical ages using measured anomalous fading rates brings them into agreement with the ¹⁴C ages. The approach used holds significant promise for improving the accuracy of infrared optical-dating methods.Luminescence data for the other three samples result in optical age limits much greater than the ¹⁴C ages. These data provide a textbook demonstration of the correlation between scatter in the luminescence intensity of individual sample aliquots and their normalization values that is expected when the samples contain sand grains not adequately exposed to daylight just prior to or during deposition and burial. Thus, the data for these three samples suggest that the tsunamis eroded young and old sand deposits before dropping the sand in the lake.     

Be isotope and B-Be investigations of the historic eruptions of Mt. Vesuvius

Cosmogenic and radioactive ¹⁰Be, stable ⁹Be and B concentrations have been determined for four alkaline lavas from the recent cycle of historical activity at Mt. Vesuvius, in the Central Campania Province (Italy). The goal of this study was to use the Be isotope and Be-B systematics of these lavas in a manner analogous to that used in regions unequivocally related to active subduction, in hopes of being able to document a subduction origin for the Vesuvian lavas. Four lavas measured to date have low ¹⁰Be concentrations as well as low B/Be ratios. While the low ¹⁰Be concentrations could reflect subduction and incorporation of old sediments and/or only contributions from the basaltic part of the subducting plate, the combination of low ¹⁰Be and low B/Be ratios must be interpreted as the absence of a subduction signature in these lavas. Unfortunately, the absence of a subduction signature in the lavas cannot be used unequivocally to argue against recent subduction in the Central Campania region. Subduction of an unusually young, hot slab which has lost its B through prograde metamorphic reactions at shallow levels could explain the absence of a subduction signature, as could extensive crustal contamination. In addition, recent studies in Java show that alkaline lavas in this region of active subduction never show a subduction signature, even when erupting in close proximity to calc-alkaline and tholeiitic lavas which do; by analogy, the Vesuvian lavas studied may have been generated in a part of the sub-arc mantle which either does not experience or does not preserve the chemical signature of subduction modification. The present data set does not allow us to make definitive statements about the role of subduction in the origin of the Central Campania lavas.