Orthotropic anisotropy: on contributions of elasticity parameters to a difference in quasi-P-wave-squared velocities resulted from propagation in two orthogonal symmetry planes

We investigate the dependence of quasi P-wave phase velocity propagating in orthotropic media on particular elasticity parameters. Specifically, due to mathematical facilitation, we consider the squared-velocity difference, , resulted from propagation in two mutually perpendicular symmetry planes. In the context of the effective medium theory, may be viewed as a parameter evaluating the influence of cracks – embedded in the background medium – parallel to one or both aforementioned planes. Our investigation is both theoretical and numerical. Based on Christoffel's equations, we propose two accurate approximations of . Due to them, we interpret the aforementioned squared-velocity difference as being twice more dependent on , than on . To describe the magnitude of the dependence, we consider the proportions between the partial derivatives of . Further, it occurs that is influenced by the ratio of vertically propagating quasi P-wave to vertically propagating quasi S-wave. Anomalously high might be caused by the low P/S ratio, which in turn can be an indicator of the presence of gas in natural fractures or aligned porosity. Also, we carry out numerical sensitivity study, according to which is approximately twice more dependent on than on , twice more sensitive to than to , and equally dependent on as on . The dependence on and can be neglected, especially for small phase angles. We verify the approximations and perform the sensitivity study, using eight examples of the elasticity tensors.     

运用尾波干涉技术监测瑞昌—阳新MS4.6地震后地壳介质变化

基于2011年瑞昌—阳新M_S4.6地震序列事件波形及其观测报告,采用双差重定位法、波形互相关法和尾波干涉法,筛选出了同时被黄梅台和九江台记录且波形互相关系数≥0.920的两组重复地震对:D₁（2011-09-11 05:44 M_L2.3和2011-10-24 08:06 M_L2.2）和D₂（2011-09-18 07:06 M_L2.8和2011-10-2408:06 M_L2.2）,发现瑞昌—阳新M_S4.6地震后1.5个月内（2011-09-11—10-24）震源区S波早期尾波波速增加。      

Correction to: Accelerated reproduction of 2-D periodic waves

Ocean Dynamics - A Correction to this paper has been published: https://doi.org/ https://doi.org/10.1007/s10236-021-01450-3     

Zur Bestimmung des gesamten anorganischen Kohlenstoffes in natürlichen Gewässern durch Titration mit Salzsäure

The knowledge of total inorganic carbon concentration (c) is important for characterizing natural waters. It is usually measured by the titration alkalinity (“m-value”) and pH which depend on temperature and ionic strength. This paper demonstrates that Ca (and Mg) can influence the calculation of from titration alkalinity, too. Errors result from neglecting this influence. In such cases the share of ion pair CaCO amounts to more than 50 % of . General relationships among the influencing factors are given by tables calculated with the help of a BASIC computer programme for calculation from titration alkalinity, pH and Ca concentration.     

Effects of variations in fluid properties and fracture geometry on dispersion,anisotropy, and reflection in media with planar fractures

Our objective was to discover the effect of variations in fluid properties and fracture geometry on the velocity of seismic wave propagation in fluid‐saturated media with parallel planar fractures. We used numerical models calculated by analytical solutions to examine the behaviour of P‐wave phase velocity dispersion in the normal direction to layering, in non‐porous and porous media with planar fractures. We also examined the anisotropy of low frequency phase and group velocities of fast and slow P‐waves and angular‐dependent reflection coefficients in media with planar fractures, under conditions of saturation by fluids with varying bulk moduli, densities, and fracture apertures. We defined several parametre , , and characterising dispersion, characterising anisotropy, characterising the difference between fast and slow modes, and R₀ and characterising reflection. Our results show that the behaviour of dispersion shows wider stopbands in the case of gas saturation. Concavity indicator of dispersion for gas saturation was greater than that for liquid saturation and is usually greater than one. Anisotropy is more sensitive to bulk modulus contrast than to density contrast between the solid and the fluid, and is more sensitive to density contrast than to bulk modulus contrast. The case of gas saturation usually had a greater negative R₀ and a greater value of compared with those of brine and heavy and light oil saturations. Our results are helpful in distinguishing fluid types saturating geophysical fractures and estimating the aperture and spacing of planar fractures. In seismic exploration, bulk modulus and fluid density can provide useful information in distinguishing among brine, oil, and gas; fracture geometry is important to estimate the permeability of reservoirs.     

Characterizing the degree of amplitude‐variation‐with‐offset nonlinearity in seismic physical modelling reflection data

The nonlinearity of the seismic amplitude‐variation‐with‐offset response is investigated with physical modelling data. Nonlinearity in amplitude‐variation‐with‐offset becomes important in the presence of large relative changes in acoustic and elastic medium properties. A procedure for pre‐processing physical modelling reflection data is enacted on the reflection from a water‐plexiglas boundary. The resulting picked and processed amplitudes are compared with the exact solutions of the plane‐wave Zoeppritz equations, as well as approximations that are first, second, and third order in , , and . In the low angle range of 0°–20°, the third‐order plane‐wave approximation is sufficient to capture the nonlinearity of the amplitude‐variation‐with‐offset response of a liquid‐solid boundary with , , and ρ contrasts of 1485–2745 m/s, 0–1380 m/s, and 1.00–1.19 gm/cc respectively, to an accuracy value of roughly 1%. This is in contrast to the linear Aki–Richards approximation, which is in error by as much as 25% in the same angle range. Even‐order nonlinear corrective terms are observed to be primarily involved in correcting the angle dependence of , whereas the odd‐order nonlinear terms are involved in determining the absolute amplitude‐variation‐with‐offset magnitudes.     

Correction to: Sea-state contributions to sea-level variability in the European Seas

Ocean Dynamics - A Correction to this paper has been published: https://doi.org/10.1007/s10236-020-01434-9     

Correction to: Benchmarking the software packages to model and assess the seismic response of unreinforced masonry existing buildings through nonlinear static analyses

Bulletin of Earthquake Engineering - A correction to this paper has been published: https://doi.org/10.1007/s10518-021-01105-0     

Regulating Environmental Factors of Nutrients Release from Wheat Straw Biochar for Sustainable Agriculture

One mechanism by which biochar application enhances soil nutrient availability is through direct nutrients release from biochar. However, factors controlling the release processes are poorly understood. In this study, the effects of pH, biochar to water ratio, temperature, ionic strength, and equilibration time on the release of PO, NO, NH, K⁺, Na⁺, Ca²⁺, and Mg²⁺ from biochar were evaluated in simulated experiments. The release of PO, K⁺, Ca²⁺, and Mg²⁺ was significantly affected by extraction pH, suggesting that their release from biochar was pH dependent or an H⁺‐consuming process. Correlation analysis indicated that PO and Ca²⁺, PO and Mg²⁺, and Ca²⁺ and Mg²⁺ were co‐solubilized with increasing soil acidity. To a lesser extent, the recovery of the nutrients was also affected by the ratio of biochar to water: more nutrients were soluble with more water supply. In contrast, the release of Na was not affected by pH while the concentration increased with decreasing biochar to water ratio. Meanwhile, other factors (temperature, ionic strength, and equilibration time) had less effect on nutrient release from biochar. Under the influence of pH, the patterns of NO and PO release from biochar were different: extractable NO concentration was not affected by the pH but more PO was released in strongly acidic conditions. Our data suggested that P was mainly retained in inorganic forms while N was in organic forms in biochar. We conclude that environmental factors have marked influences on nutrients release from biochar.     

Correction to: Initial time dependence of wind- and density-driven Lagrangian residual velocity in a tide-dominated bay

Ocean Dynamics - A Correction to this paper has been published: https://doi.org/10.1007/s10236-021-01459-8     

Diel patterns in coastal‐stream nitrate concentrations linked to evapotranspiration in the riparian zone of a low‐relief,agricultural catchment

Evapotranspiration (ET) can cause diel fluctuations in the elevation of the water table and the stage in adjacent streams. The diel fluctuations of water levels change head gradients throughout the day, causing specific discharge through near‐stream sediment to fluctuate at the same time scale. In a previous study, we showed that specific discharge controls the residence time of groundwater in streambed sediment that, in turn, exerted the primary control on removal from groundwater passing through the streambed. In this study, we examine the magnitude of diel specific discharge patterns through the streambed driven by ET in the riparian zone with a transient numerical saturated–unsaturated groundwater flow model. On the basis of a first‐order kinetic model for removal, we predicted diel fluctuations in stream concentrations. Model results indicated that ET drove a diel pattern in specific discharge through the streambed and riparian zone (the removal zones). Because specific discharge is inversely proportional to groundwater travel time through the removal zones and travel time determines the extent of removal, diel changes in ET can result in a diel pattern in concentration in the stream. The model predictions generally matched observations made during summertime base‐flow conditions in a small coastal plain stream in Virginia. A more complicated pattern was observed following a seasonal drawdown period, where source components to the stream changed during the receding limb of the hydrograph and resulted in diel fluctuations being superimposed over a multi‐day trend in concentrations. Copyright © 2013 John Wiley & Sons, Ltd.     

Catalytic Reduction of Hexaminecobalt(III) by Pitch‐based Spherical Activated Carbon (PBSAC)

The wet ammonia (NH₃) desulfurization process can be retrofitted to remove nitric oxide (NO) and sulfur dioxide (SO₂) simultaneously by adding soluble cobalt(II) salt into the aqueous ammonia solution. Activated carbon is used as a catalyst to regenerate hexaminecobalt(II), Co(NH₃), so that NO removal efficiency can be maintained at a high level for a long time. In this study, the catalytic performance of pitch‐based spherical activated carbon (PBSAC) in the simultaneous removal of NO and SO₂ with this wet ammonia scrubbing process has been studied systematically. Experiments have been performed in a batch stirred cell to test the catalytic characteristics of PBSAC in the catalytic reduction of hexaminecobalt(III), Co(NH₃). The experimental results show that PBSAC is a much better catalyst in the catalytic reduction of Co(NH₃) than palm shell activated carbon (PSAC). The Co(NH₃) reduction reaction rate increases with PBSAC when the PBSAC dose is below 7.5 g/L. The Co(NH₃) reduction rate increases with its initial concentration. Best Co(NH₃) conversion is gained at a pH range of 2.0–6.0. A high temperature is favorable to such reaction. The intrinsic activation energy of 51.00 kJ/mol for the Co(NH₃) reduction catalyzed by PBSAC has been obtained. The experiments manifest that the simultaneous elimination of NO and SO₂ by the hexaminecobalt solution coupled with catalytic regeneration of hexaminecobalt(II) can maintain a NO removal efficiency of 90% for a long time.     

Time‐lapse pre‐stack seismic data registration and inversion for CO2 sequestration study at Cranfield

Pre‐stack seismic data are indicative of subsurface elastic properties within the amplitude versus offset characteristic and can be used to detect elastic rock property changes caused by injection. We perform time‐lapse pre‐stack 3‐D seismic data analysis for monitoring sequestration at Cranfield. The time‐lapse amplitude differences of Cranfield datasets are found entangled with time‐shifts. To disentangle these two characters, we apply a local‐correlation‐based warping method to register the time‐lapse pre‐stack datasets, which can effectively separate the time‐shift from the time‐lapse seismic amplitude difference without changing the original amplitudes. We demonstrate the effectiveness of our registration method by evaluating the inverted elastic properties. These inverted time‐lapse elastic properties can be reliably used for monitoring plumes.     

Spatial average aquifer recharge through atmospheric chloride mass balance and its uncertainty in continental Spain

The atmospheric chloride mass balance (CMB) method allows spatial evaluations of the average diffuse aquifer recharge by rainfall () in large and varied territories when long‐term steady conditions can be assumed. Often, the distributed average CMB variables necessary to calculate have to be estimated from the available variable‐length data series, which may be of suboptimal quality and spatial coverage. This paper explains the use of these data and the reliability of the results in continental Spain, chosen as a large and varied territory. The CMB variables have been regionalized by ordinary kriging at the same 4976 nodes of a 10 km × 10 km grid. Nodal values vary from 14 to 810 mm year^–1, 90% ranging from 30 to 300 mm year^–1. The recharge‐to‐precipitation ratios vary from 0.03 in low‐permeability formations and semiarid areas to 0.65 in some carbonate massifs. Integrated average results for the whole of continental Spain yield a potential aquifer recharge of 64 km³ year^?1, the net recharge over permeable formations (40% of the territory) being 32 km³ year^?1. Two main sources of uncertainty affecting (given by the coefficient of variation, CV), induced by the inherent natural variability of the variables (CV_R) and from mapping (), have been segregated. The average CV_R is 0.13 and could be improved with longer data series. The average is 0.07 and may be decreased with better data coverage. The estimates were compared with other regional and local recharge estimates, being 4% and 1% higher, respectively. Copyright © 2012 John Wiley & Sons, Ltd.     

Chemical dynamics of N-containing ionic species in a boreal forest snowcover during the spring melt period

A study of the changes in the ionic loads of NO, NH, SO and H⁺ in a boreal forest snowpack at Lake Laflamme, Québec was carried out using hydrological and chemical data from field lysimeters. The results showed that depletion of the N-containing species occurs periodically in the snowpack during meltwater discharge. Rain-on-snow events led to in-pack losses of NO and NH at a rate of 130 μeq m^?2 day^?1 and 101·3 μeq m^?2day^?1 respectively. On dry days, however, dry deposition and deposition of organic debris from the canopy resulted in increases of 183·3 μeq m^?2day^?1 for NO and 4·5 μeq m^?2day^?1 for NH in the pack. In contrast, SO₄^2? showed continual in-pack increases due to deposition of 5·0 μeq m^?2day^?1 for wet days and 92·6 μeq m^?2day^?1 for dry days. The depletion of NO and NH is due to microbiological uptake of these nutrients during periods when the free water content of the pack is high. Controlled melts in a laboratory snowmelt simulator containing snow and organic matter from the forest canopy at Lake Laflamme showed losses of NO and NH similar to those observed in the field. As the microbiological uptake proceeds at a rate comparable to that of ionic load increases in the pack by dry deposition, models of the chemical dynamics of snowmelt should take the former into account in any system where organic content of the snowpack is appreciable.     

Adsorption of Complex Pollutants from Aqueous Solutions by Nanocomposite Materials

In view of water pollutants becoming more complex, both anionic and cationic pollutants need to be removed. The multi‐pollutants simultaneous removal is paid more and more attention. Hence, development composite materials for treatment complex wastewater are the aim of this study. In this research, iron–nickel nanoparticles deposited onto aluminum oxide (α‐Al₂O₃) and carbon nanotubes (CNTs) to form nanocomposite materials Fe–Ni/Al₂O₃ and Fe–Ni/CNTs, respectively, were used as adsorbents. The adsorption capacities of Fe–Ni/Al₂O₃ and Fe–Ni/CNTs for AO7, HSeO, and Pb²⁺ were observed to be 5.46, 8.28, 27.02, and 25.6 mg/g, 15.29 and 17.12 mg/g, separately. The composite materials with negative charges were superior in adsorption of anionic pollutants. Using orthogonal experimental design and analysis of variance to co‐treat dye AO7, HSeO and Pb²⁺ in aqueous solutions, seven testing factors were included: (1) adsorbent types, (2) amount of iron, (3) solution pHs, (4) AO7 concentrations, (5) Pb²⁺ concentrations, (6) HSeO concentrations and (7) reaction time. The experimental results showed that the removal of complex pollutants AO7, HSeO, and Pb²⁺ on Fe–Ni/CNTs could reach up to 90% in the optimal treatment conditions. When using Fe–Ni/CNTs as the adsorbent, the sorption isothermals were well fitted in the Freundlich isotherm, and R² could reach up to 0.98.     

Multi-parameter reflection waveform inversion for acoustic transversely isotropic media with a vertical symmetry axis

Full waveform inversion in transversely isotropic media with a vertical symmetry axis provides an opportunity to better match the data at the near and far offsets. However, multi-parameter full waveform inversion, in general, suffers from serious cycle-skipping and trade-off problems. Reflection waveform inversion can help us recover a background model by projecting the residuals of the reflected wavefield along the reflection wavepath. Thus, we extend reflection waveform inversion to acoustic transversely isotropic media with a vertical symmetry axis utilizing the proper parameterization for reduced parameter trade-off. From a radiation patterns analysis, an acoustic transversely isotropic media with a vertical symmetry axis is better described by a combination of the normal-moveout velocity and the anisotropic parameters η and δ for reflection waveform inversion applications. We design a three-stage inversion strategy to construct the optimal resulting model. In the first stage, we only invert for the background by matching the simulated reflected wavefield from the perturbations of and δ with the observed reflected wavefield. In the second stage, the background and η are optimized simultaneously and the far-offset reflected wavefield mainly contribute to their updates. We perform Born modelling to compute the reflected wavefield for the two stages of reflection waveform inversion. In the third stage, we perform full waveform inversion for the acoustic transversely isotropic media with a vertical symmetry axis to delineate the high-wavenumber structures. For this stage, the medium is described by a combination of the horizontal velocity , η and ε instead of , η and δ. The acoustic multi-parameter full waveform inversion utilizes the diving waves to improve the background as well as utilizes reflection for high-resolution information. Finally, we test our inversion algorithm on the modified Sigsbee 2A model (a salt free part) and a two-dimensional line from a three-dimensional ocean bottom cable dataset. The results demonstrate that the proposed reflection waveform inversion approach can recover the background model for acoustic transversely isotropic media with a vertical symmetry axis starting from an isotropic model. This recovered background model can mitigate the cycle skipping of full waveform inversion and help the inversion recover higher resolution structures.     

3D angle-domain common-image gathers for migration velocity analysis

Angle‐domain common‐image gathers (ADCIGs) are an essential tool for migration velocity analysis (MVA). We present a method for computing ADCIGs in 3D from the results of wavefield‐continuation migration. The proposed methodology can be applied before or after the imaging step in a migration procedure. When computed before imaging, 3D ADCIGs are functions of the offset ray parameters (p, p) ; we derive the geometric relationship that links the offset ray parameters to the aperture angle γ and the reflection azimuth φ. When computed after imaging, 3D ADCIGs are directly produced as functions of γ and φ. The mapping of the offset ray parameters (p, p) into the angles (γ, φ) depends on both the local dips and the local interval velocity; therefore, the transformation of ADCIGs computed before imaging into ADCIGs that are functions of the actual angles is difficult in complex structure. By contrast, the computation of ADCIGs after imaging is efficient and accurate even in the presence of complex structure and a heterogeneous velocity function. On the other hand, the estimation of the offset ray parameters (p, p) is less sensitive to velocity errors than the estimation of the angles (γ, φ). When ADCIGs that are functions of the offset ray parameters (p, p) are adequate for the application of interest (e.g. ray‐based tomography), the computation of ADCIGs before imaging might be preferable. Errors in the migration velocity cause the image point in the angle domain to shift along the normal to the apparent geological dip. By assuming stationary rays (i.e. small velocity errors), we derive a quantitative relationship between this normal shift and the traveltime perturbation caused by velocity errors. This relationship can be directly used in an MVA procedure to invert depth errors measured from ADCIGs into migration velocity updates. In this paper, we use it to derive an approximate 3D residual moveout (RMO) function for measuring inconsistencies between the migrated images at different γ and φ. We tested the accuracy of our kinematic analysis on a 3D synthetic data set with steeply dipping reflectors and a vertically varying propagation velocity. The tests confirm the accuracy of our analysis and illustrate the limitations of the straight‐ray approximation underlying our derivation of the 3D RMO function.     

Solutes in overland flow following fire in eucalyptus and pine forests,northern Portugal

Small bounded plots were used to investigate dissolved nutrient losses in overland flow in eucalyptus and pine forests in northern Portugal following understorey fires. Losses of NON, POCa²⁺, Mg²⁺ and K⁺ were measured over 19 months during the first 2–3 years after fire. Solute losses in overland flow increased after fire in eucalyptus and pine terrain due to the mineralization of litter and vegetation and increased overland flow. The data suggest that (i) elevated losses persist for at least 2 and 3 years at the pine and eucalyptus sites respectively and (ii) soluble losses of POK⁺ in a post‐fire cycle may exceed those adsorbed to eroded material. Losses of POalthough relatively small, are potentially the most detrimental to soil fertility and forest productivity because of the limited opportunity for P replenishment and correspondingly high P losses adsorbed to eroded sediment. Copyright © 2000 John Wiley & Sons, Ltd.     

Development of duration‐dependent damage‐based inelastic response spectra

A simple relationship is proposed in this paper to construct damage‐based inelastic response spectra including the effect of ground motion duration that it can be used for damage control in seismic design of structures. This relation is established for three groups of ground motions with short‐duration, moderate‐duration, and long‐duration ranges. To develop the model, the duration effect is included in the cyclic ductility of structures by an energy‐based method, and then strength reduction factors are computed based on this modified ductility (named ). The strength reduction factors were calculated for 44 stiffness‐degrading oscillators having vibration periods between 0.05 and 4.0 s, four ultimate ductility capacities, and five damage levels subjected to 296 earthquake records. The results showed that ductility capacity, damage level, and ground motion duration are effective parameters in the energy dissipation of structures, which affect the spectra. The values of short‐period oscillators (e.g., low‐rise structures) under short‐duration records are generally greater than those under moderate‐duration and long‐duration records. Residual analysis has been made in terms of magnitude and distance to examine the validity of the proposed simple expression. Finally, the introduced spectra were compared with three previously published proposals. Copyright © 2016 John Wiley & Sons, Ltd.