Combination of GRACE monthly gravity field solutions from different processing strategies

We combine the publicly available GRACE monthly gravity field time series to produce gravity fields with reduced systematic errors. We first compare the monthly gravity fields in the spatial domain in terms of signal and noise. Then, we combine the individual gravity fields with comparable signal content, but diverse noise characteristics. We test five different weighting schemes: equal weights, non-iterative coefficient-wise, order-wise, or field-wise weights, and iterative field-wise weights applying variance component estimation (VCE). The combined solutions are evaluated in terms of signal and noise in the spectral and spatial domains. Compared to the individual contributions, they in general show lower noise. In case the noise characteristics of the individual solutions differ significantly, the weighted means are less noisy, compared to the arithmetic mean: The non-seasonal variability over the oceans is reduced by up to 7.7% and the root mean square (RMS) of the residuals of mass change estimates within Antarctic drainage basins is reduced by 18.1% on average. The field-wise weighting schemes in general show better performance, compared to the order- or coefficient-wise weighting schemes. The combination of the full set of considered time series results in lower noise levels, compared to the combination of a subset consisting of the official GRACE Science Data System gravity fields only: The RMS of coefficient-wise anomalies is smaller by up to 22.4% and the non-seasonal variability over the oceans by 25.4%. This study was performed in the frame of the European Gravity Service for Improved Emergency Management (EGSIEM; http://www.egsiem.eu) project. The gravity fields provided by the EGSIEM scientific combination service (ftp://ftp.aiub.unibe.ch/EGSIEM/) are combined, based on the weights derived by VCE as described in this article.     

Rockfall risk management using a pre-failure deformation database

Terrestrial laser scanning (TLS) monitoring has been used to estimate the location, volume, and kinematics of a variety of small magnitude rockfalls before failure (1–1000 m³ range), and in some cases, potential failure time has been assessed through the application of inverse velocity methods. However, our current understanding of rock slope pre-failure behavior for this magnitude range and prediction ability is based on observations of a small number of failure case histories. In this study, a pre-failure deformation database was constructed for rockfall volumes exceeding 0.1 m³, observed over a 1252-day study interval at the Goldpan rock slope, British Columbia, Canada, in order to better understand the pre-failure behavior of rock slopes and provide an empirical means of estimating temporal failure ranges. Repeated TLS datasets were acquired at an average scanning interval of 2–3 months. A total of 90 rockfall events were recorded at this site, during this time period, of which 64 (71%) exhibited measurable deformation prior to failure. Classification of rockfalls by volume suggests that a scale dependency may exist, as deformation was detected for a greater proportion of rockfalls >?5 m³ (92%) than for smaller rockfalls in the range of 0.1–0.5 m³ (61%). A lower rate of pre-failure deformation detection was also reported for planar sliding failures as compared with wedge or toppling failures, suggesting that deformation was less easily detected for these failure types. This study proposes and implements a framework for rockfall assessment and forecasting that does not require continuous monitoring of deformation.     

Revisiting the thermodynamics of hardening plasticity for unsaturated soils

A thermodynamically consistent extension of the constitutive equations of saturated soils to unsaturated conditions is often worked out through the use of a unique ‘effective’ interstitial pressure, accounting equivalently for the pressures of the saturating fluids acting separately on the internal solid walls of the pore network. The natural candidate for this effective interstitial pressure is the space averaged interstitial pressure. In contrast experimental observations have revealed that, at least, a pair of stress state variables was needed for a suitable framework to describe stress–strain–strength behaviour of unsaturated soils. The thermodynamics analysis presented here shows that the most general approach to the behaviour of unsaturated soils actually requires three stress state variables: the suction, which is required to describe the invasion of the soil by the liquid water phase through the retention curve; two effective stresses, which are required to describe the soil deformation at water saturation held constant. However a simple assumption related to the plastic flow rule leads to the final need of only a Bishop-like effective stress to formulate the stress–strain constitutive equation describing the soil deformation, while the retention properties still involve the suction and possibly the deformation. Commonly accepted models for unsaturated soils, that is the Barcelona Basic Model and any approach based on the use of an effective averaged interstitial pressure, appear as special extreme cases of the thermodynamic formulation proposed here.     

Computation of the space and time evolution of equilibrium-line altitudes on Andean glaciers (10°N–55°S)

A previous study of Fox [Fox, A.N. 1993. Snowline altitude and climate at present and during the Last Pleistocene Glacial Maximum in the Central Andes (5°–28°S). Ph.D. Thesis. Cornell University.] showed that for a fixed 0 °C isotherm altitude, the equilibrium-line altitude (ELA) of the Peruvian and Bolivian glaciers from 5 to 20°S can be expressed based on a log–normal expression of local mid-annual rainfall amount (P). In order to extrapolate the function to the whole Andes (10°N to 55°S) a local 0 °C isotherm altitude is introduced. Two applications of this generalised function are presented. One concerns the space evolution of mean inter-annual ELA for three decades (1961–1990) over the whole South American continent. A high-resolution data set (grid data: 10′ for latitude/longitude) of mean monthly air surface temperature and precipitation is used. Mean annual values over the 1961–1990 period were calculated. On each grid element, the mean annual 0 °C isotherm altitude is determined from an altitudinal temperature gradient and mean annual temperature (T) at ground level. The 0 °C isotherm altitude is then associated with the annual precipitation amount to compute the ELA. Using computed ELA and the digital terrain elevation model GTOPO30, we determine the extent of the glacierised area in Andean regions under modern climatic conditions. The other application concerns the ELA time evolution on Zongo Glacier (Bolivia), where inter-annual ELA variations are computed from 1995 to 1999. For both applications, the computed values of ELA are in good agreement with those derived from glacier mass balance measurements.     

Dating cleavage formation in slates and phyllites with the 40Ar/39Ar laser microprobe: an example from the western New England Appalachians, USA

Determinations of the absolute age of cleavage formation can provide fundamental information about the evolution of orogenic belts. However, when applied to cleavages in slates and phyllites, conventional dating methods are complicated by problems related to mineral separation and the presence of multiple cleavage generations. In situ high-spatial-resolution ⁴⁰Ar/³⁹Ar laser microprobe geochronology and microstructural observations indicate that the age of cleavage formation in slates and phyllites can be constrained by analysing zones of tightly packed cleavage domains. Three regionally developed cleavages (S₂, S₃, and S₄) are present in the northern Taconic Allochthon of Vermont and New York. Representative samples were studied from a variety of localities where these cleavages, which are defined by white micas, are well developed. In the suite of samples, only S₃ and S₄ are expressed as domains that are sufficiently wide and spatially isolated in thin section to permit quantitative ⁴⁰Ar/³⁹Ar geochronology. Mean ⁴⁰Ar/³⁹Ar laser microprobe ages for these domains are 370.7 ± 1.0 Myr for S₃ and 345.5 ± 1.7 Myr for S₄. Because estimates of the Ar closure temperature for white micas are substantially higher than the inferred growth temperatures of the micas defining S₃ and S₄, these values are interpreted as periods since cleavage formation. This interpretation is consistent with independent geochronological constraints on the age of the Acadian orogeny in the region.     

Arsenic Geochemistry and Hydrostratigraphy in Midwestern U.S. Glacial Deposits

Arsenic concentrations exceeding the U.S. EPA's 10 μg/L standard are common in glacial aquifers in the midwestern United States. Previous studies have indicated that arsenic occurs naturally in these aquifers in association with metal-(hydr)oxides and is released to groundwater under reducing conditions generated by microbial oxidation of organic matter. Despite this delineation of the arsenic source and mechanism of arsenic mobilization, identification of arsenic-impacted aquifers is hindered by the heterogeneous and discontinuous nature of glacial sediments. In much of the Midwest, the hydrostratigraphy of glacial deposits is not sufficiently characterized to predict where elevated arsenic concentrations are likely to occur. This case study from southeast Wisconsin presents a detailed characterization of local stratigraphy, hydrostratigraphy, and geochemistry of the Pleistocene glacial deposits and underlying Silurian dolomite. Analyses of a single core, water chemistry data, and well construction reports enabled identification of two aquifers separated by an organic-rich aquitard. The upper, unconfined aquifer provides potable water, whereas arsenic generally exceeds 10 μg/L in the deeper aquifer. Although coring and detailed hydrostratigraphic characterization are often considered impractical, our results demonstrate that a single core improved interpretation of the complex lithology and hydrostratigraphy. This detailed characterization of hydrostratigraphy facilitated development of well construction guidelines and lays the ground work for further studies of the complex interactions among aquifer sediments, hydrogeology, water chemistry, and microbiology that lead to elevated arsenic in groundwater.     

Osmium,tungsten, and chromium isotopes in sediments and in Ni‐rich spinel at the K‐T boundary: Signature of a chondritic impactor

Abstract— It is now established that a large extraterrestrial object hit the Earth at the end of the Cretaceous period, about 65 Ma ago. We have investigated Re‐Os, Hf‐W, and Mn‐Cr isotope systems in sediments from the Cretaceous and the Paleogene in order to characterize the type of impactor. Within the Cretaceous‐Tertiary (K‐T) boundary layer, extraterrestrial material is mixed with terrestrial material, causing a dilution of the extraterrestrial isotope signature that is difficult to quantify. A phase essentially composed of Ni‐rich spinel, formed in the atmosphere mainly from melted projectile material, is likely to contain the extraterrestrial isotopic signature of the impactor. We show that the analysis of spinel is indeed the best approach to determine the initial isotope composition of the impactor, and that W and Cr isotopes confirm that the projectile was a carbonaceous chondrite.     

Analyse régionale des simulations climatiques du modèle canadien de circulation générale de l'atmosphère pour le territoire québécois

Abstract

Current understanding of the possible nature of climatic change at the regional scale is limited by the spatial resolution of General Circulation Models (GCM). The use of GCM outputs without correction linked to the spatial variability of the variables can bring significant errors in their utilization at the regional scale. The potential of the Canadian GCM for regional applications in Quebec has been analysed by comparison to the climatic normals of temperature and precipitation, measured over the Quebec climatological network, on an annual and seasonal basis. This analysis has been undertaken with the support of a geographical information system (GIS) (PAMAP). In summary, a difference between the climatic normal and the GCM output has been estimated at 20% for temperature and 30% for precipitation. We present an analysis of a corrected regionalized scenario for the province of Quebec of the possible climatic change simulated by the Canadian GCM under the hypothesis of a doubling of atmospheric CO₂. Results show an increase of the annual average temperature of 4° C for summer and 6°C for winter, associated with an average increase of 80 mm (10%) in annual precipitation, reaching 25% in some regions.