A series solution for surface motion ampli-fication due to underground group cavities: Incident P waves

Introduction In recent years some big cities in China commonly build subways in central urban areas to improve crowded surface transportation, and the subways usually are large-scale group cavities composed by two parallel cavities. Construction of large-scale group cavities will certainly affect the design ground motion near the cavities, and further affect seismic safety of the existing build-ings nearby. So it is of significance to estimate this effect on seismic design of the buildings. H…     

The Multiscale TROPIcal CatchmentS critical zone observatory M-TROPICS dataset III: Hydro-geochemical monitoring of the Mule Hole catchment,south India

Despite the importance of tropical ecosystems for climate regulation, biodiversity, water and nutrient cycles, only a few Critical Zone Observatories (CZOs) are located in the tropics. Among these, most are in humid climates, while very few data exist for semi-arid and sub-humid climates, due to the difficulty of estimating hydro-geochemical balances in catchments with ephemeral streams. We contribute to fill this gap by presenting a meteorological and hydro-geochemical dataset acquired at the Mule Hole catchment (4.1 km²), a pristine dry deciduous forest located in a biosphere reserve in south India. The dataset consists of time series of variables related to (i) meteorology, including rainfall, air temperature, relative humidity, wind speed and direction, and global radiation, (ii) hydrology, including water level and discharge at the catchment outlet, (iii) hydrogeology, including manual (monthly) and/or automated (from 15 min to hourly) groundwater levels in nine piezometers and (iv) geochemistry, including suspended sediment content in the stream and chemical composition of rainfall (event based), groundwater (monthly sampling) and stream water (storm events, 15 min to hourly frequency with an automatic sampler). The time series extend from 2003 to 2019. Measurement errors are minimized by frequent calibration of sensors and quality checks, both in the field and in the laboratory. Despite these precautions, several data gaps exist, due to occasional access restriction to the site and instrument destruction by wildlife. Results show that large seasonal and interannual variations of climatic conditions were reflected in the large variations of stream flow and groundwater recharge, as well as in water chemical composition. Notably, they reveal a long-term evolution of groundwater storage, suggesting hydrogeological cycles on a decadal scale. This dataset, alone or in combination with other data, has already allowed to better understand water and element cycling in tropical dry forests, and the role of forest diversity on biogeochemical cycles. As tropical ecosystems are underrepresented by Critical Zone Observatories, we expect this data note to be valuable for the global scientific community.     

Scattering and reflection of SH waves around a slope on an elastic wedged space

The scattering and reflection of SH waves by a slope on an elastic wedged space is investigated. A series solution is obtained by using the wave function expansion method. The slope on a wedged space is divided into two subregions by an artificial, auxiliary circular arc. The wave fields with unknown complex coefficients within each sub-region are derived. Applying Graf addition theorem, the scattered waves in the sub-regions are expressed in a global coordinate system. Fourier transform is adopted to derive a consistent form of standing waves in the inner region using the orthogonality of the cosine functions. The boundary-valued problem is solved by stress and displacement continuity along the artificial, auxiliary arc to obtain the unknown complex coefficients. Parametric studies are next performed to investigate how the topography from the slope on the wedged space will affect the scattering and diffraction, and hence the amplification and de-amplification of the SH waves. Numerical results show that the surface motions on the slope of the wedged space is influenced greatly by the topography. Amplification of the surface motions near the slope vertex is significant. The corresponding phases along the wedged space surfaces are consistent with the direction that the SH waves are propagating.     

A grain-size correction for metal pollution indexes in river sediments

According to the List of Hazardous Substances of the Agency of the Toxic Substances and Disease Registry(ATSDR 2017),some metallic elements such as Pb or Cd are still considered as the most polluting elements in the world.These elements accumulate in sediments,and there are various methods available that differentiate lithogenic sources from anthropogenic sources.For that,the natural geochemical background is required,however,its definition is far from unanimous:it can be a global or local value,single or statistically estimated value,depending on the studies.Our study is focused on the Suba
e River in the state of Bahia,Northeastern Brazil,which was historically contaminated by lead metallurgy.The river sediments were sampled at 21 locations:one at the discharge point of the plant's rainwater basin,ten upstream,and ten downstream.The total contents of Fe,Ti,Mn,Pb,Zn,Cu,Cr,particle-size distribution,and organic matter were analyzed.The conventional geochemical indexes,i.e.,the pollution degree(mCd),pollution load index,and enrichment factor indicated contamination by Zn,Pb,Cr,and Cu.Then,a new enrichment factor is introduced,assuming that the sediments come from erosion of soils in the watershed.After obtaining the natural concentrations of metals in the clay,silt,and sand fractions of uncontaminated reference soils,we corrected the pollution indexes by calculating a new individual background in each sediment sample,according to its grain size distribution.This new approach provided more precise indexes in the assessment of sediment pollution,by highlighting higher contamination of Zn and Pb(around 50%)and at the same time the absence of Cr and Cu contamination,two metals not involved with the metallurgical activity.     

The influence of bed forms on reference concentration and suspension under waves and currents

Hydrodynamic, suspension and bed-form measurements were made 2 km off the Dutch coast near Noordwijk aan Zee in ∼14 m water depth for a period of 32 days in 2003. Tidal currents were just able to suspend sand at the bed at peak spring tide but most suspension and transport occurred as a result of the combination of waves and currents. Burst-average (17 min) sand concentration profiles (C¯-profiles) from an acoustic backscatter instrument were used to define the (varying) location of the sea-bed, following the method used by Green et al. [Green, M.O., Dolphin, T.J., Swales, A., Vincent, C.E., 1999. Transport of mixed-size sediments in a tidal channel. Coastal Sediments ‘99, edited by N.C. Kraus, and W.G. McDougal, ASCE, Long Island, New York, pp. 644–658]. Reference concentrations at the sea-bed (C₀) and at 1 cm (C₁) were examined in relation to both the hydrodynamic conditions and the type of bed forms present. The C₀ predictive equations of Green and Black [Green, M.O., Black, K.P., Suspended sediment reference concentration under waves: field measurements and critical analysis of two predictive models, Coastal Engineering, 38, 115–141, 1999](short-wave ripples) and Nielsen [Nielsen, P., Suspended sediment concentrations under waves, Coastal Engineering, 10, 23–31, 1986](all bed forms; includes ripple steepness), both of which require knowledge of the bed-form type, were not as successful in explaining the variance in our C₀ data as a regression of C₀ against the skin-friction Shields parameter θ′_cw that ignored bed-form type (73% of variance explained). The values of the reference concentration C₁ were compared with the Lee et al. [Lee, G.-H., Dade, W.B., Friedrichs, C.T., Vincent, C.E., Examination of Reference Concentration Under Waves and Currents on the Inner Shelf., Journal of Geophysical Research, 109, 1–10, 2004] equation which predicts C₁ from the product of the Shields parameter and the inverse Rouse parameter; 51% of the variance in C₁ was explained.     

Closed‐Form Approximations for Two‐Dimensional Groundwater Age Patterns in a Fresh Water Lens

Simple closed‐form approximations are presented for calculating the steady‐state groundwater age distribution in two‐dimensional vertical cross sections of idealized fresh water lenses overlying salt water, for aquifers that are vertically semi‐infinite and of finite thickness. The approximations are developed on the basis of existing one‐dimensional analytical solutions for travel‐time calculation in fresh water lenses and approximate streamline formulations. The two‐dimensional age distributions based on the closed‐form solutions match convincingly with numerical simulations. As expected, notable deviations from the numerical solution are encountered at the groundwater flow divide and when submarine groundwater discharge occurs. Ratios of recharge over hydraulic conductivities are varied to explore how the magnitude of the deviations changes, and it is found that the approximate closed‐form solutions perform well over a range of conditions found in natural systems.     

On the Use of Hydrological Models and Satellite Data to Study the Water Budget of River Basins Affected by Human Activities: Examples from the Garonne Basin of France

Natural and anthropogenic forcing factors and their changes significantly impact water resources in many river basins around the world. Information on such changes can be derived from fine scale in situ and satellite observations, used in combination with hydrological models. The latter need to account for hydrological changes caused by human activities to correctly estimate the actual water resource. In this study, we consider the catchment area of the Garonne river (in France) to investigate the capabilities of space-based observations and up-to-date hydrological modeling in estimating water resources of a river basin modified by human activities and a changing climate. Using the ISBA–MODCOU and SWAT hydrological models, we find that the water resources of the Garonne basin display a negative climate trend since 1960. The snow component of the two models is validated using the moderate-resolution imaging spectroradiometer snow cover extent climatology. Crop sowing dates based on remote sensing studies are also considered in the validation procedure. Use of this dataset improves the simulated evapotranspiration and river discharge amounts when compared to conventional data. Finally, we investigate the benefit of using the MAELIA multi-agent model that accounts for a realistic agricultural and management scenario. Among other results, we find that changes in crop systems have significant impacts on water uptake for agriculture. This work constitutes a basis for the construction of a future modeling framework of the sociological and hydrological system of the Garonne river region.     

Nitrogen isotopes in ophiolitic metagabbros: A re-evaluation of modern nitrogen fluxes in subduction zones and implication for the early Earth atmosphere

Nitrogen contents and isotope compositions together with major and trace element concentrations were determined in a sequence of metagabbros from the western Alps (Europe) in order to constrain the evolution and behavior of N during hydrothermal alteration on the seafloor and progressive dehydration during subduction in a cold slab environment (8 °C/km). The rocks investigated include: (i) low-strain metagabbros that equilibrated under greenschist to amphibolite facies (Chenaillet Massif), blueschist facies (Queyras region) and eclogite facies (Monviso massif) conditions and (ii) highly-strained mylonites and associated eclogitic veins from the Monviso Massif. In all samples, nitrogen (2.6-55 ppm) occurs as bound ammonium () substituting for K or Na-Ca in minerals. Cu concentrations show a large variation, from 73.2 to 6.4 ppm, and are used as an index of hydrothermal alteration on the seafloor because of Cu fluid-mobility at relatively high temperature (>300 °C). In low-strain metagabbros, δ¹⁵N values of +0.8‰ to +8.1‰ are negatively correlated with Cu concentrations. Eclogitic mylonites and veins display Cu concentrations lower than 11 ppm and show a δ¹⁵N-Cu relationship that does not match the δ¹⁵N-Cu correlation found in low-strain rocks. This δ¹⁵N-Cu correlation preserved in low-strain rocks is best interpreted by leaching of Cu-N compounds, possibly of the form Cu(NH₃)₂²⁺, during hydrothermal alteration. Recognition that the different types of low-strain metagabbros show the same δ¹⁵N-Cu correlation indicates that fluid release during subduction zone metamorphism did not modify the original N and Cu contents of the parent hydrothermally-altered metagabbros. In contrast, the low Cu content present in eclogitic veins and mylonites implies that ductile deformation and veining were accompanied either by a loss of copper or that externally-derived nitrogen was added to the system.We estimate the global annual flux of N subducted by metagabbros as 4.2 (±2.0) × 10¹¹ g/yr. This value is about half that of sedimentary rocks, which suggests that gabbros carry a significant portion of the subducted nitrogen. The net budget between subducted N and that outgassed at volcanic arcs indicates that ∼80% of the subducted N is not recycled to the surface. On a global scale, the total amount of N buried to the mantle via subduction zones is estimated to be three times higher than that released from the mantle via mid-ocean ridges, arc and intraplate volcanoes and back-arc basins. This implies that N contained in Earth surface reservoirs, mainly in the atmosphere, is progressively transferred and sequestered into the mantle, with a net flux of ∼9.6 × 10¹¹ g/yr. Assuming a constant flux of subducted N over the Earth’s history indicates that an amount equivalent to the present atmospheric N may have been sequestered into the silicate Earth over a period of 4 billion years.     

Velocity model estimation for depth imaging: Comparison of three tomography methods on a 2D real data set

The estimation of velocity models is still crucial in seismic reflection imaging as it controls the quality of the depth‐migrated image, which is the basis of geological interpretation. Among the numerous existing methods for velocity determination, tomographic methods are very attractive for their efficiency and ability to retrieve heterogeneities of the medium. We present three tomographic methods in order to estimate heterogeneous velocity models from 2D prestack PP reflection data: a traveltime tomography in the time‐migrated domain, a traveltime and slope tomography in the non‐migrated time domain, and a slope tomography in the depth‐migrated domain. The first method (traveltime tomography in the time domain) is based on continuous picked events, whereas the two slope tomographic methods, one in the time domain and the other in the depth domain, are based on locally coherent events, with no assumptions about reflector geometry or the unknown velocity field. The purpose of this paper is not to describe in detail the theoretical basis and implementation of the methods, but to apply and compare their output using the same marine real data set. Based on the estimated velocity models, the migrated images and the common‐image gathers from the three processing routes, the relative strengths and weaknesses of the methods are discussed. Finally, similarities are indicated and potential alternative approaches are proposed.