Lateglacial and Holocene palaeoenvironments at Sierra La Madera,northwestern Mexico: Sharp signals of Younger Dryas and North American monsoon

During the last glacial termination, the climate system experienced intense global variations whose causes and impacts are not fully defined, particularly for low latitudes. The northwestern Mexico Sky Islands present a climate-sensitive ideal setting to record palaeoecological and palaeoclimatic changes due to their physiographic complexity and location in the ecotone between temperate and tropical ecosystems. High-resolution pollen analysis and a detailed sedimentological study were conducted at the Ciénega Tonibabi tropical thorn scrub site. The 15 540–0 cal a bp nearly continuous record shows that the North Atlantic Ocean did have a cold and humid climatic influence during the glacial stages of the end of the Pleistocene, including a sharp pulse during the Younger Dryas. However, a shift to the Pacific Ocean influence occurred during the Holocene, which led to the development of the El Niño conditions prevailing today. Colder and warmer phases follow one another with higher or lower winter precipitation, including a sharp Bølling–Allerød and development and intensification of the North American monsoon. They are reflected in hydrological changes as well as in the advances, retreats and intermingling of coniferous forests and tropical thorn scrub.     

Dynamical constraints on the Local Group from the CMB and 2MRS dipoles

We place constraints on the dynamics of the Local Group (LG) by comparing the dipole of the cosmic microwave background (CMB) with the peculiar velocity induced by the Two Micron All-Sky Redshift Survey galaxy sample. The analysis is limited by the lack of surveyed galaxies behind the zone of avoidance (ZoA). We therefore allow for a component of the LG velocity due to unknown mass concentrations behind the ZoA, as well as for an unknown transverse velocity of the Milky Way relative to the Andromeda galaxy. We infer extra motion along the direction of the Galactic Centre (where Galactic confusion and dust obscuration peaks) at the 95 per cent significance level. With a future survey of the ZoA it might be possible to constrain the transverse velocity of the Milky Way relative to Andromeda.     

Towards an accurate real-time locator of infrasonic sources

Infrasonic signals propagate from an atmospheric source via media with stochastic and fast space-varying conditions. Hence, their travel time, the amplitude at sensor recordings and even manifestation in the so-called “shadow zones” are random. Therefore, the traditional least-squares technique for locating infrasonic sources is often not effective, and the problem for the best solution must be formulated in probabilistic terms. Recently, a series of papers has been published about Bayesian Infrasonic Source Localization (BISL) method based on the computation of the posterior probability density function (PPDF) of the source location, as a convolution of a priori probability distribution function (APDF) of the propagation model parameters with likelihood function (LF) of observations. The present study is devoted to the further development of BISL for higher accuracy and stability of the source location results and decreasing of computational load. We critically analyse previous algorithms and propose several new ones. First of all, we describe the general PPDF formulation and demonstrate that this relatively slow algorithm might be among the most accurate algorithms, provided the adequate APDF and LF are used. Then, we suggest using summation instead of integration in a general PPDF calculation for increased robustness, but this leads us to the 3D space-time optimization problem. Two different forms of APDF approximation are considered and applied for the PPDF calculation in our study. One of them is previously suggested, but not yet properly used is the so-called “celerity-range histograms” (CRHs). Another is the outcome from previous findings of linear mean travel time for the four first infrasonic phases in the overlapping consecutive distance ranges. This stochastic model is extended here to the regional distance of 1000 km, and the APDF introduced is the probabilistic form of the junction between this travel time model and range-dependent probability distributions of the phase arrival time picks. To illustrate the improvements in both computation time and location accuracy achieved, we compare location results for the new algorithms, previously published BISL-type algorithms and the least-squares location technique. This comparison is provided via a case study of different typical spatial data distributions and statistical experiment using the database of 36 ground-truth explosions from the Utah Test and Training Range (UTTR) recorded during the US summer season at USArray transportable seismic stations when they were near the site between 2006 and 2008.     

Numerical study on basin-edge effects in the seismic response of the Gubbio valley,Central Italy

The Gubbio basin in Central Italy is a intermountain basin of extensional tectonic origin, typical of Central and Southern Apennines, characterized by moderate seismicity. The strongest recorded event within the area is a magnitude 5.7 earthquake which occurred on 29 April 1984 along the Gubbio fault, bordering the eastern side of the basin. The main objective of this study is to analyze the features of earthquake ground motion as related to basin-edge effects, by performing physics-based numerical simulations of the 1984 earthquake through a high-performance spectral element code. The simulated ground motions are found in reasonable agreement with the recorded motions when using the kinematic source model developed by Ameri et al. (Bull Seismol Soc Am 99:647–663, 2009), with a rise-time equal to 1 s and a nucleation point located in the middle of the fault. Pronounced differences were noted between records from the basin and adjacent sites at outcropping bedrock, owing to both the strong impedance contrast between soft alluvial sites and bedrock formations (lithostratigraphic amplification), as well as lateral discontinuities related to the 2D/3D geometry of the basin (generation of surface waves). Since the fault was located beneath the basin, 1D amplification effects were found to be more relevant than those associated with the generation of surface waves from the basin edge. Finally, an envelope delay spectrum was computed for the simulated ground motions, showing that surface waves are excited in the frequency band of 0.2–0.8 Hz with a significant increase of ground motion duration within the basin.     

Measurement of surface water runoff from plots of two different sizes

Intensities and amounts of water infiltration and runoff on sloping land are governed by the rainfall pattern and soil hydraulic conductivity, as well as by the microtopography and soil surface conditions. These components are closely interrelated and occur simultaneously, and their particular contribution may change during a rainfall event, or their effects may vary at different field scales. The scale effect on the process of infiltration/runoff was studied under natural field and rainfall conditions for two plot sizes: small plots of 0·25 m² and large plots of 50 m². The measurements were carried out in the central region of Chile in a piedmont most recently used as natural pastureland. Three blocks, each having one large plot and five small plots, were established. Cumulative rainfall and runoff quantities were sampled every 5 min. Significant variations in runoff responses to rainfall rates were found for the two plot sizes. On average, large plots yielded only 40% of runoff quantities produced on small plots per unit area. This difference between plot sizes was observed even during periods of continuous runoff. Copyright © 2002 John Wiley & Sons, Ltd.     

Structure and dynamics of undercurrents in the western boundary current of the Bay of Bengal

Ocean Dynamics - The structure and variability of undercurrents in the East India Coastal Current (EICC), which is the western boundary current system in the Bay of Bengal (BoB), and the mechanisms...     

Century-scale nitrogen and phosphorus controls of the carbon cycle

In recent decades, humans have become a very important force in the Earth system, demonstrating that emissions (gaseous, liquid, and solid) are the cause of many of our environmental issues. These emissions are responsible for major global reorganizations of the biogeochemical cycles. The oceans are now a net sink of atmospheric CO₂, whereas in their preindustrial state they were a source; the trophic state of the coastal oceans is progressively moving toward increased heterotrophy; and the terrestrial realm is now vacillating between trophic states, whereas in preindustrial times it was autotrophic. In this paper, we present model calculations that underscore the role of human-induced perturbations in changing Earth's climate, specifically the role of anthropogenic nitrogen and phosphorus in controlling processes in the global carbon cycle since the year 1850 with projections to the year 2035. Our studies show that since the late 1940's emissions of nitrogen and phosphorus have been sequestered in the terrestrial living phytomass and groundwater. This nutrient-enhanced fertilization of terrestrial biota, coupled with rising atmospheric CO₂ and global temperature, has induced a sink of anthropogenic CO₂ that roughly balances the emission of CO₂ owing to land use change. In the year 2000, for example, the model-calculated terrestrial biotic sink was 1730 Mtons C/year, while the emission of CO₂ from changes in land use was 1820 Mtons C/year, a net flux of 90 Mtons C/year emitted to the atmosphere. In the global aquatic environment, enhanced terrestrial inputs of biotically reactive phosphorus (about 8.5 Mtons P/year) and inorganic nitrogen (about 54 Mtons N/year), have induced increased new production and burial of organic carbon in marine sediments, which is a small sink of anthropogenic CO₂. It is predicted that the response of the global land reservoirs of C, N, and P to sustained anthropogenic perturbations will be maintained in the same direction of change over the range of projected scenarios of global population increase and temperature change for the next 35 years. The magnitude of change is significantly larger when the global temperature increase is maximum, especially with respect to the processes of remobilization of the biotically important nutrients nitrogen and phosphorus.     

Coupled silicon–oxygen isotope fractionation traces Archaean silicification

Silica alteration zones and cherts are a conspicuous feature of Archaean greenstone belts worldwide and provide evidence of extensive mobilisation of silica in the marine environment of the early Earth. In order to understand the process(es) of silicification we measured the silicon and oxygen isotope composition of sections of variably silicified basalts and overlying bedded cherts from the Theespruit, Hooggenoeg and Kromberg Formations of the Barberton Greenstone Belt, South Africa.The δ³⁰Si and δ¹⁸O values of bulk rock increase with increasing amount of silicification from unsilicified basalts (?0.64‰ < δ³⁰Si < ?0.01‰ and + 8.6‰ < δ¹⁸O < + 11.9‰) to silicified basalts (δ³⁰Si and δ¹⁸O values as high as + 0.81‰ and + 15.6‰, respectively). Cherts generally have positive isotope ratios (+ 0.21‰ < δ³⁰Si < + 1.05‰ and + 10.9 < δ¹⁸O < + 17.1), except two cherts, which have negative δ³⁰Si values, but high δ¹⁸O (up to + 19.5‰).The pronounced positive correlations between δ³⁰Si, δ¹⁸O and SiO₂ imply that the isotope variation is driven by the silicification process which coevally introduced both ¹⁸O and ³⁰Si into the basalts. The oxygen isotope variation in the basalts from about 8.6‰ to 15.6‰ is likely to represent temperature-dependent isotope fractionation during alteration. Our proposed model for the observed silicon isotope variation relies on a temperature-controlled basalt dissolution vs. silica deposition process.     

GT0 Explosion Sources for IMS Infrasound Calibration: Charge Design and Yield Estimation from Near-source Observations

Three large-scale on-surface explosions were conducted by the Geophysical Institute of Israel (GII) at the Sayarim Military Range, Negev desert, Israel: about 82 tons of strong high explosives in August 2009, and two explosions of about 10 and 100 tons of ANFO explosives in January 2011. It was a collaborative effort between Israel, CTBTO, USA and several European countries, with the main goal to provide fully controlled ground truth (GT0) infrasound sources, monitored by extensive observations, for calibration of International Monitoring System (IMS) infrasound stations in Europe, Middle East and Asia. In all shots, the explosives were assembled like a pyramid/hemisphere on dry desert alluvium, with a complicated explosion design, different from the ideal homogenous hemisphere used in similar experiments in the past. Strong boosters and an upward charge detonation scheme were applied to provide more energy radiated to the atmosphere. Under these conditions the evaluation of the actual explosion yield, an important source parameter, is crucial for the GT0 calibration experiment. Audio-visual, air-shock and acoustic records were utilized for interpretation of observed unique blast effects, and for determination of blast wave parameters suited for yield estimation and the associated relationships. High-pressure gauges were deployed at 100–600 m to record air-blast properties, evaluate the efficiency of the charge design and energy generation, and provide a reliable estimation of the charge yield. The yield estimators, based on empirical scaled relations for well-known basic air-blast parameters—the peak pressure, impulse and positive phase duration, as well as on the crater dimensions and seismic magnitudes, were analyzed. A novel empirical scaled relationship for the little-known secondary shock delay was developed, consistent for broad ranges of ANFO charges and distances, which facilitates using this stable and reliable air-blast parameter as a new potential yield estimator. The delay data of the 2009 shot with IMI explosives, characterized by much higher detonation velocity, are clearly separated from ANFO data, thus indicating a dependence on explosive type. This unique dual Sayarim explosion experiment (August 2009/January 2011), with the strongest GT0 sources since the establishment of the IMS network, clearly demonstrated the most favorable westward/eastward infrasound propagation up to 3,400/6,250 km according to appropriate summer/winter weather pattern and stratospheric wind directions, respectively, and thus verified empirically common models of infrasound propagation in the atmosphere.     

Vegetation structure controls on snow and soil moisture in restored ponderosa pine forests

Thinning of semi-arid forests to reduce wildfire risk is believed to improve forest health by increasing soil moisture. Increased snowpack, reduced transpiration and reduced rainfall interception are frequently cited mechanisms by which reduced canopy density may increase soil moisture. However, the relative importance of these factors has not been rigorously evaluated in field studies. We measured snow depth, snow water equivalent (SWE) and the spatial and temporal variation in soil moisture at four experimental paired treatment-control thinning sites in high elevation ponderosa pine forest northern Arizona, USA. We compared snow and soil moisture measurements with forest structure metrics derived from aerial imagery and 3-dimensional lidar data to determine the relationship between vegetation structure, snow and soil moisture throughout the annual hydrologic cycle. Soil moisture was consistently and significantly higher in thinned forest plots, even though the treatments were performed 8–11 years before this study. However, we did not find evidence that SWE was higher in thinned forests across a range of snow conditions. Regression tree analysis of soil moisture and vegetation structure data provided some evidence that localized differences in transpiration and interception of precipitation influence the spatial pattern of soil moisture at points in the annual hydrologic cycle when the system is becoming increasingly water limited. However, vegetation structure explained a relatively low amount of the spatial variance (R² < 0.23) in soil moisture. Continuous measurements of soil moisture in depth profiles showed stronger attenuation of soil moisture peaks in thinned sites, suggesting differences in infiltration dynamics may explain the difference in soil moisture between treatments as opposed to overlying vegetation alone. Our results show limited support for commonly cited relationships between vegetation structure, snow and soil moisture and indicate that future research is needed to understand how reduction in tree density alters soil hydraulic properties.