Feasibility study for an anisotropic full waveform inversion of cross‐well seismic data

Anisotropy is often observed due to the thin layering or aligned micro‐structures, like small fractures. At the scale of cross‐well tomography, the anisotropic effects cannot be neglected. In this paper, we propose a method of full‐wave inversion for transversely isotropic media and we test its robustness against structured noisy data. Optimization inversion techniques based on a least‐square formalism are used. In this framework, analytical expressions of the misfit function gradient, based on the adjoint technique in the time domain, allow one to solve the inverse problem with a high number of parameters and for a completely heterogeneous medium. The wave propagation equation for transversely isotropic media with vertical symmetry axis is solved using the finite difference method on the cylindrical system of coordinates. This system allows one to model the 3D propagation in a 2D medium with a revolution symmetry. In case of approximately horizontal layering, this approximation is sufficient. The full‐wave inversion method is applied to a crosswell synthetic 2‐component (radial and vertical) dataset generated using a 2D model with three different anisotropic regions. Complex noise has been added to these synthetic observed data. This noise is Gaussian and has the same amplitude f?k spectrum as the data. Part of the noise is localized as a coda of arrivals, the other part is not localized. Five parameter fields are estimated, (vertical) P‐wave velocity, (vertical) S‐wave velocity, volumetric mass and the Thomsen anisotropic parameters epsilon and delta. Horizontal exponential correlations have been used. The results show that the full‐wave inversion of cross‐well data is relatively robust for high‐level noise even for second‐order parameters such as Thomsen epsilon and delta anisotropic parameters.     

In situ oxygen uptake rates by coastal sediments under the influence of the Rhône River (NW Mediterranean Sea)

The influence of riverine inputs on biogeochemical cycling and organic matter recycling in sediments on the continental shelf off the Rhône River mouth (NW Mediterranean Sea) was investigated by measuring sediment oxygen uptake rates using a combination of in situ and laboratory techniques. Four stations were investigated during two cruises in June 2001 and June 2002, with depths ranging from 9 to 192 m and over a distance to the Rhône River mouth ranging from 4 to 36 km. Diffusive oxygen uptake (DOU) rates were determined using an in situ sediment microprofiler and total oxygen uptake (TOU) rates were measured using sediment core incubations. There was good agreement between these two techniques which indicates that the non-diffusive fraction of the oxygen flux was minimal at the investigated stations. DOU rates ranged from 3.7±0.4 mmol O₂ m⁻² d⁻¹ at the continental shelf break to 19.3±0.5 mmol O₂ m⁻² d⁻¹ in front of the Rhône River mouth. Sediment oxygen uptake rates mostly decreased with increasing depth and with distance from the Rhône mouth. The highest oxygen uptake rate was observed at 63 m on the Rhône prodelta, corresponding to intense remineralization of organic matter. This oxygen uptake rate was much larger than expected for the increasing bathymetry, which indicates that biogeochemical cycles and benthic deposition are largely influenced by the Rhône River inputs. This functioning was also supported by the detailed spatial distribution of total organic carbon (TOC), total nitrogen (TN) and C/N atomic ratio in surficial sediments. Sediments of the Rhône prodelta are enriched in organic carbon (2–2.2%) relative to the continental shelf sediments (<1%) and showed C/N ratios exceeding Redfield stoichiometry for fresh marine organic matter. A positive exponential correlation was found between DOU and TOC contents (r²=0.98, n=4). South-westward of the Rhône River mouth, sediments contained highly degraded organic matter of both terrestrial and marine origin, due to direct inputs from the Rhône River, sedimentation of marine organic matter and organic material redeposition after resuspension events.     

Climate impacts of recent multidecadal changes in Atlantic Ocean Sea Surface Temperature: a multimodel comparison

During the twentieth century sea surface temperatures in the Atlantic Ocean exhibited prominent multidecadal variations. The source of such variations has yet to be rigorously established—but the question of their impact on climate can be investigated. Here we report on a set of multimodel experiments to examine the impact of patterns of warming in the North Atlantic, and cooling in the South Atlantic, derived from observations, that is characteristic of the positive phase of the Atlantic Multidecadal Oscillation (AMO). The experiments were carried out with six atmospheric General Circulation Models (including two versions of one model), and a major goal was to assess the extent to which key climate impacts are consistent between the different models. The major climate impacts are found over North and South America, with the strongest impacts over land found over the United States and northern parts of South America. These responses appear to be driven by a combination of an off-equatorial Gill response to diabatic heating over the Caribbean due to increased rainfall within the region and a Northward shift in the Inter Tropical Convergence Zone (ITCZ) due to the anomalous cross-equatorial SST gradient. The majority of the models show warmer US land temperatures and reduced Mean Sea Level Pressure during summer (JJA) in response to a warmer North Atlantic and a cooler South Atlantic, in line with observations. However the majority of models show no significant impact on US rainfall during summer. Over northern South America, all models show reduced rainfall in southern hemisphere winter (JJA), whilst in Summer (DJF) there is a generally an increase in rainfall. However, there is a large spread amongst the models in the magnitude of the rainfall anomalies over land. Away from the Americas, there are no consistent significant modelled responses. In particular there are no significant changes in the North Atlantic Oscillation (NAO) over the North Atlantic and Europe in Winter (DJF). Additionally, the observed Sahel drying signal in African rainfall is not seen in the modelled responses. Suggesting that, in contrast to some studies, the Atlantic Multidecadal Oscillation was not the primary driver of recent reductions in Sahel rainfall.     

A Note on k - ε Modelling of Vegetation Canopy Air-Flows

The k - turbulence model is a standard of computational software packages for engineering, yet its application to canopy turbulence has not received comparable attention. This is probably due to the additional source (and/or sink) terms, whose parameterization remained uncertain. This model must include source terms for both turbulent kinetic energy (k) and the viscous dissipation rate (), to account for vegetation wake turbulence budget. In this note, we show how Kolmogorov's relation allows for an analytical solution to be calculated within the portion of a dense and homogeneous canopy where the mixing length does not vary. By substitution within model equations, this solution allows for a set of constraints on source term model coefficients to be derived.Those constraints should meet both Reynolds averaged Navier–Stokes equationsand large-eddy simulation sub-grid scale turbulence modelling requirements.Although originating from within a limited portion of the canopy, the predictedcoefficients values must be valid elsewhere in order to make the model capable of predicting the whole canopy-layer flow with a single set of constants.     

Weathering of sandstone lotus petals at the Angkor site: a 1,000-year stone durability trial

This article is based on field investigations carried out at Ta Keo and Banteay Srei temples (Angkor, Cambodia), which have similar environmental histories and display the same recurrent decorative pattern (lotus petal) in three different sandstone lithotypes. Based on two original scales of mechanical and chemical weathering, the degree of deterioration of c.2000 sandstone petals was visually assessed, resulting in the proposal of a durability scale. An explanatory scheme was provided, based on laboratory analyses of stone materials (ICP-AES, XRD, SEM-EDS, polarizing microscopy and mercury injection). The key drivers of sandstone durability are both mineralogical and petrophysical. The remarkable preservation of the pink sandstone is due to the quartzose nature of its clasts, the non-swelling clay content of its abundant matrix (kaolinite) and its porosity structure (large connected pores and probable non-connected porosity). The overall good preservation of the green sandstone is due to its extremely low porosity and specific surface area, associated with its strong cementation. The severe deterioration of the grey sandstone is due to its laminated structure (oriented biotites), the susceptibility to hygric swelling and solution of its binding agents (chlorite/smectite and calcite) and its porosity structure, which evolves during a two-stage weathering sequence. The first stage is characterized by a bimodal porosity due to the partial clogging of macropores by swelling clays, and the second stage by a substantial increase in porosity and the penetration of salts (barite and anhydrite). Conversely, oxidation phenomena associated with weathering can lead to a drastic reduction of porosity.     

易破碎粒状材料本构研究

粒状材料被广泛地应用到土木工程的各个领域。大部分粒状材料在外力作用下较容易产生颗粒破碎,颗粒破碎对材料的力学性能有很大的影响。为了更好地描述易破碎粒状材料的力学特性,基于弹塑性力学和临界状态土力学开发了一种能够考虑颗粒破碎效应的本构模型。此本构模型能够考虑在剪应力和压应力作用下引起的颗粒破碎对粒状材料力学性能的影响。为了考虑颗粒破碎的影响,基于对Cambria 砂的高压试验,得出了临界状态线位置与消耗塑性功之间的关系,并称之为“破碎方程”。此本构模型拥有双屈服面,分别考虑剪切和等向压缩产生的塑性变形。通过试验结果与数值计算结果的对比得出,新的本构模型能较好地描述易破碎粒状材料的力学特性     

Diagenesis and the reconstruction of paleoenvironments: A method to restore original δO values of carbonate and phosphate from fossil tooth enamel

Intra-tooth δ¹⁸O variations within the carbonate (δ¹⁸Oc) and phosphate (δ¹⁸Op) components of tooth apatite were measured for Miocene and Pliocene hypsodont mammals from Afghanistan, Greece and Chad in order to evaluate the resistance of enamel to diagenetic alteration. Application of water-apatite interaction models suggest that the different kinetic behaviours of the phosphate-water and carbonate-water systems can be used to detect subtle oxygen isotope disequilibria in fossil enamel when intra-individual variations are considered. Selective alteration of the oxygen isotope composition from the carbonate component of Afghan and Greek enamels suggests inorganic isotopic exchange processes. Microbially-induced isotopic exchange for phosphate is demonstrated for the first time in enamel samples from Chad, in association with extensive recrystallization. In Chad, δ¹⁸Op values were derived from partial isotopic exchange with fossil groundwater during early diagenesis. Mass balance calculations using average carbonate content in enamel as a proxy for recrystallization, and the lowest δ¹⁸Op value of dentine as a proxy for the isotopic composition of the diagenetic fluid, indicate that diagenesis can alter δ¹⁸Op by as much as 3‰ in some enamel samples. This diagenetic alteration is also responsible for a decrease in intra-individual variations of up to 1‰ in affected specimens. The effects of diagenesis on δ¹⁸Op values of fossil enamel are not systematic, however, and can only be estimated if sequential δ¹⁸Op and δ¹⁸Oc analyses are performed on fossil enamel and dentine. Reconstruction of large temporal- or spatial-scale paleoclimates based on δ¹⁸Op analyses from mammalian teeth cannot be considered valid if enamel has been affected by bacterial activity or if the data cannot be corrected for diagenetic effects.     

Nanoscale zero-valent iron flakes for groundwater treatment

Even today the remediation of organic contaminant source zones poses significant technical and economic challenges. Nanoscale zero-valent iron (NZVI) injections have proved to be a promising approach especially for source zone treatment. We present the development and the characterization of a new kind of NZVI with several advantages on the basis of laboratory experiments, model simulations and a field test. The developed NZVI particles are manufactured by milling, consist of 85 % Fe(0) and exhibit a flake-like shape with a thickness of <100 nm. The mass normalized perchloroethylene (PCE) dechlorination rate constant was 4.1 × 10^?3 L/g h compared to 4.0 × 10^?4 L/g h for a commercially available reference product. A transport distance of at least 190 cm in quartz sand with a grain size of 0.2–0.8 mm and Fe(0) concentrations between 6 and 160 g/kg (sand) were achieved without significant indications of clogging. The particles showed only a low acute toxicity and had no longterm inhibitory effects on dechlorinating microorganisms. During a field test 280 kg of the iron flakes was injected to a depth of 10–12 m into quaternary sand layers with hydraulic conductivities ranging between 10^?4 and 10^?5 m/s. Fe(0) concentrations of 1 g/kg (sand) or more [up to 100 g/kg (sand)] were achieved in 80 % of the targeted area. The iron flakes have so far remained reactive for more than 1 year and caused a PCE concentration decrease from 20.000–30.000 to 100–200 µg/L. Integration of particle transport processes into the OpenGeoSys model code proved suitable for site-specific 3D prediction and optimization of iron flake injections.     

In situ and high frequency monitoring of suspended sediment properties using a spectrophotometric sensor

It is well known that sediment properties, including sediment‐associated chemical constituents and sediment physical properties, can exhibit significant variations within and between storm runoff events. However, the number of samples included in suspended sediment studies is often limited by time‐consuming and expensive laboratory procedures after stream water sampling. This restricts high frequency sampling campaigns to a limited number of events and reduces accuracy when aiming to estimate fluxes and loads of sediment‐associated chemical constituents. In this study, we address the potential of a portable ultraviolet–visible spectrophotometer (220–730 nm) to estimate suspended sediment properties in a resource efficient way. Several field deployable spectrophotometers are currently available for in‐stream measurements of environmental variables at high temporal resolution. These instruments have primarily been developed and used to quantify solute concentrations (e.g. dissolved organic carbon and NO₃‐N), total concentrations of dissolved and particulate forms (e.g. total organic carbon) and turbidity. Here we argue that light absorbance values can be calibrated to estimate sediment properties. We present light absorbance data collected at 15‐min intervals in the Weierbach catchment (NW Luxembourg, 0.45 km²) from December 2013 to January 2015. In this proof‐of‐concept study, we performed a local calibration using suspended sediment loss‐on‐ignition (LOI) measurements as an example of suspended sediment property. We assessed the performance of several regression models that relate light absorbance measurements with the percentage weight LOI. The MM‐robust regression method presented the lowest standard error of prediction (0.48%) and was selected for calibration (adjusted r² = 0.76 between observed and predicted values). The model was then used to predict LOI during a storm runoff event in December 2014. This study demonstrates that spectrophotometers can be used to estimate suspended sediment properties at high temporal resolution and for long‐time spans in a simple, non‐destructive and affordable manner. Copyright © 2016 John Wiley & Sons, Ltd.     

Fe(II)-Na(I)-Ca(II) Cation Exchange on Montmorillonite in Chloride Medium: Evidence for Preferential Clay Adsorption of Chloride – Metal Ion Pairs in Seawater

Fe(II)–Ca(II), Fe(II)–Na(I), and Fe(II)–Ca(II)–Na(I) exchange experiments on montmorillonite were performed in chloride background. These experiments show the possible sorption of Fe²⁺ and FeCl⁺ ion pairs in exchange site positions, a result confirmed with 77 K ⁵⁷Fe Mössbauer experiments. The sorption data were modeled and the cation exchange selectivity for Fe(II) were found to be nearly equal to that of Ca(II). Vanselow selectivity coefficients, for Na–Fe²⁺ and Na–FeCl⁺ reactions, were found to be equal to 0.4 (0.5 for Ca²⁺) and 2.3 (2.5 for CaCl⁺) respectively. High affinity of montmorillonite for chloride ion pairs seems to be a common mechanism as first stated by Sposito et al., (Soil Sci. Soc. Am. J. 47, 51–56, 1983a), and should have implications e.g., on the chemistry of suspended particles in seawater. Exchange selectivity coefficients derived from this study and others were used to model experimental data on river water and seawater equilibrated particles. The agreement between simulations and experimental data is very good. The simulation shows the predominance of monovalent ion (Na⁺ and chloride ion pairs) sorption on clay particles in seawater. This sorption of monovalent ions leads to the dispersion of particles in seawater and to the extension of a plume of particles spreading away from river deltas, such as that of the River Amazon.