Target-level waveform inversion: a prospective application of the convolution-type representation for the acoustic wavefield

Nowadays, full-waveform inversion, based on fitting the measured surface data with modelled data, has become the preferred approach to recover detailed physical parameters from the subsurface. However, its application is computationally expensive for large inversion domains. Furthermore, when the subsurface has a complex geological setting, the inversion process requires an appropriate pre-conditioning scheme to retrieve the medium parameters for the desired target area in a reliable manner. One way of dealing with both aspects is by waveform inversion schemes in a target-oriented fashion. Therefore, we propose a prospective application of the convolution-type representation for the acoustic wavefield in the frequency–space domain formulated as a target-oriented waveform inversion method. Our approach aims at matching the observed and modelled upgoing wavefields at a target depth level in the subsurface, where the seismic wavefields, generated by sources distributed above this level, are available. The forward modelling is performed by combining the convolution-type representation for the acoustic wavefield with solving the two-way acoustic wave-equation in the frequency–space domain for the target area. We evaluate the effectiveness of our inversion method by comparing it with the full-domain full-waveform inversion process through some numerical examples using synthetic data from a horizontal well acquisition geometry, where the sources are located at the surface and the receivers are located along a horizontal well at the target level. Our proposed inversion method requires less computational effort and, for this particular acquisition, it has proven to provide more accurate estimates of the target zone below a complex overburden compared to both full-domain full-waveform inversion process and local full-waveform inversion after applying interferometry by multidimensional deconvolution to get local-impulse responses.     

C-O bonds in carbonate minerals: A new kind of paleothermometer

The abundance of the doubly substituted CO₂ isotopologue, ¹³C¹⁸O¹⁶O, in CO₂ produced by phosphoric acid digestion of synthetic, inorganic calcite and natural, biogenic aragonite is proportional to the concentration of ¹³C-¹⁸O bonds in reactant carbonate, and the concentration of these bonds is a function of the temperature of carbonate growth. This proportionality can be described between 1 and 50 °C by the function: Δ₄₇ = 0.0592 · 10⁶ · T⁻² − 0.02, where Δ₄₇ is the enrichment, in per mil, of ¹³C¹⁸O¹⁶O in CO₂ relative to the amount expected for a stochastic (random) distribution of isotopes among all CO₂ isotopologues, and T is the temperature in Kelvin. This relationship can be used for a new kind of carbonate paleothermometry, where the temperature-dependent property of interest is the state of ordering of ¹³C and ¹⁸O in the carbonate lattice (i.e., bound together vs. separated into different CO₃²⁻ units), and not the bulk δ¹⁸O or δ¹³C values. Current analytical methods limit precision of this thermometer to ca. ± 2 °C, 1σ. A key feature of this thermometer is that it is thermodynamically based, like the traditional carbonate-water paleothermometer, and so is suitable for interpolation and even modest extrapolation, yet is rigorously independent of the δ¹⁸O of water and δ¹³C of DIC from which carbonate grew. Thus, this technique can be applied to parts of the geological record where the stable isotope compositions of waters are unknown. Moreover, simultaneous determinations of Δ₄₇ and δ¹⁸O for carbonates will constrain the δ¹⁸O of water from which they grew.     

Carbonate clumped isotope thermometry of deep-sea corals and implications for vital effects

Here we calibrate the carbonate clumped isotope thermometer in modern deep-sea corals. We examined 11 specimens of three species of deep-sea corals and one species of a surface coral spanning a total range in growth temperature of 2-25 °C. External standard errors for individual measurements ranged from 0.005‰ to 0.011‰ (average: 0.0074‰) which corresponds to ∼1-2 °C. External standard errors for replicate measurements of Δ₄₇ in corals ranged from 0.002‰ to 0.014‰ (average: 0.0072‰) which corresponds to 0.4-2.8 °C. We find that skeletal carbonate from deep-sea corals shows the same relationship of Δ₄₇ (the measure of ¹³C-¹⁸O ordering) to temperature as does inorganic calcite. In contrast, the δ¹³ C and δ¹⁸O values of these carbonates (measured simultaneously with Δ₄₇ for every sample) differ markedly from equilibrium with seawater; i.e., these samples exhibit pronounced ‘vital effects’ in their bulk isotopic compositions. We explore several reasons why the clumped isotope compositions of deep-sea coral skeletons exhibit no evidence of a vital effect despite having large conventional isotopic vital effects.     

An Isotope Dilution ICP‐MS Method for the Determination of Mg/Ca and Sr/Ca Ratios in Calcium Carbonate

Mg/Ca and Sr/Ca ratios in calcium carbonate are important components of many palaeoclimate studies. We present an isotope dilution method relying on a single mixed spike containing ²⁵Mg, ⁴³Ca and ⁸⁷Sr. Dozens of samples per day, as small as 10 μg of carbonate, could be dissolved, spiked and run in an ICP‐MS with a precision of 0.8% (2 RSD). Two instruments types, a sector field and a quadrupole ICP‐MS, were compared. The best long term precision found was 0.4% (2 RSD), although this increased by up to a factor of two when samples of very different Mg or Sr content were run together in the same sequence. Long term averages for the two instruments concurred. No matrix effects were detected for a range of Ca concentrations between 0.2 and 2 mmol l^‐1. Accuracy, tested by measuring synthetic standard solutions, was 0.8% with some systematic trends. We demonstrate the strength of this isotope dilution method for (a) obtaining accurate results for sample sets that present a broad Mg and Sr range and (b) testing solid carbonates as candidate reference materials for interlaboratory consistency. Mg/Ca and Sr/Ca results for reference materials were in good agreement with values from the literature.     

ROSA: A High-cadence,Synchronized Multi-camera Solar Imaging System

The Rapid Oscillations in the Solar Atmosphere (ROSA) instrument is a synchronized, six-camera high-cadence solar imaging instrument developed by Queen’s University Belfast. The system is available on the Dunn Solar Telescope at the National Solar Observatory in Sunspot, New Mexico, USA, as a common-user instrument. Consisting of six 1k × 1k Peltier-cooled frame-transfer CCD cameras with very low noise (0.02 – 15 e s⁻¹ pixel⁻¹), each ROSA camera is capable of full-chip readout speeds in excess of 30 Hz, or 200 Hz when the CCD is windowed. Combining multiple cameras and fast readout rates, ROSA will accumulate approximately 12 TB of data per 8 hours observing. Following successful commissioning during August 2008, ROSA will allow for multi-wavelength studies of the solar atmosphere at a high temporal resolution.     

Propagation of drought through groundwater—a new approach using linear reservoir theory

The effect of drought on groundwater heads and discharge is often complex and poorly understood. Therefore the propagation of a drought from groundwater recharge to discharge and the influence of aquifer characteristics on the propagation was analysed by tracking a drought in recharge through a linear reservoir. The recharge was defined as a sinusoid function with a period of 1 year. The decrease in recharge owing to drought was simulated by multiplying the recharge during 1 year with a drought fraction between 0 and 1, which represents a decrease in the recharge of 100 to 0%, respectively. The droughts were identified using the threshold level approach, with a threshold that is constant in time. For this case analytical formulations were derived, which express the drought duration and deficit in the groundwater discharge in terms of the decrease in recharge, the reservoir coefficient that characterizes aquifer properties and the height of the threshold level. The results showed that the delay in the groundwater system caused a shift of the main part of the decrease in recharge from the high‐flow to the low‐flow period. This resulted in an increase in drought deficit for discharge compared with the drought deficit for recharge. Also the development of multiyear droughts caused an increase in drought deficit. The attenuation in the groundwater system caused a decrease in drought deficit. In most cases the net effect of these processes was an increase of drought deficit as a result of the propagation through groundwater. Only for small droughts the deficit decreased from recharge to discharge. The amount of increase or decrease depends on the reservoir coefficient and the severity of the drought. Under most conditions a maximum in the drought deficit occurred for a reservoir coefficient of around 200 days. Copyright © 2003 John Wiley & Sons, Ltd.     

Fe xiii emission lines in active region spectra obtained with the Solar Extreme-Ultraviolet Research Telescope and Spectrograph

Recent fully relativistic calculations of radiative rates and electron impact excitation cross-sections for Fe  xiii are used to generate emission-line ratios involving 3s²3p²–3s3p³ and 3s²3p²–3s²3p3d transitions in the 170–225 and 235–450 Å wavelength ranges covered by the Solar Extreme-Ultraviolet Research Telescope and Spectrograph (SERTS). A comparison of these line ratios with SERTS active region observations from rocket flights in 1989 and 1995 reveals generally very good agreement between theory and experiment. Several new Fe  xiii emission features are identified, at wavelengths of 203.79, 259.94, 288.56 and 290.81 Å. However, major discrepancies between theory and observation remain for several Fe  xiii transitions, as previously found by Landi and others, which cannot be explained by blending. Errors in the adopted atomic data appear to be the most likely explanation, in particular for transitions which have 3s²3p3d ¹D₂ as their upper level. The most useful Fe  xiii electron-density diagnostics in the SERTS spectral regions are assessed, in terms of the line pairs involved being (i) apparently free of atomic physics problems and blends, (ii) close in wavelength to reduce the effects of possible errors in the instrumental intensity calibration, and (iii) very sensitive to changes in N _e over the range  10⁸–10¹¹ cm⁻³  . It is concluded that the ratios which best satisfy these conditions are 200.03/202.04 and 203.17/202.04 for the 170–225 Å wavelength region, and 348.18/320.80, 348.18/368.16, 359.64/348.18 and 359.83/368.16 for 235–450 Å.     

Automated detection and tracking of solar magnetic bright points

Magnetic bright points (MBPs) in the internetwork are among the smallest objects in the solar photosphere and appear bright against the ambient environment. An algorithm is presented that can be used for the automated detection of the MBPs in the spatial and temporal domains. The algorithm works by mapping the lanes through intensity thresholding. A compass search, combined with a study of the intensity gradient across the detected objects, allows the disentanglement of MBPs from bright pixels within the granules. Object growing is implemented to account for any pixels that might have been removed when mapping the lanes. The images are stabilized by locating long-lived objects that may have been missed due to variable light levels and seeing quality. Tests of the algorithm, employing data taken with the Swedish Solar Telescope, reveal that ≈90 per cent of MBPs within a  75 × 75  arcsec² field of view are detected.     

Assessment of aquifer recharge and groundwater availability in a semiarid region of Brazil in the context of an interbasin water transfer scheme

Particularly in arid and semiarid areas, more and more populations rely almost entirely on imported water. However, the extent to which intentional discharge into transiting river systems and unintentional leakage may be augmenting water resources for communities along and down gradient of the water transfer scheme has not previously been subject to research. The objective of this study was to assess both the potential of a large-scale water transfer (WT) scheme to increase groundwater availability by channel transmission losses in a large dryland aquifer system (2,166 km²) in Brazil, and the capability of the receiving streams to transport water downstream under a prolonged drought. An integrated surface-water/groundwater model was developed to improve the estimation of the groundwater resources, considering the spatio-temporal variability of infiltrated rainfall for aquifer recharge. Aquifer recharge from the WT scheme was simulated under prolonged drought conditions, applying an uncertainty analysis of the most influential fluxes and parameters. The annual recharge (66 mm/year) was approximately twice the amount of water abstracted (1990–2016); however, the annual recharge dropped to 13.9 mm/year from 2012 to 2016, a drought period. Under similar drought conditions, the additional recharge (6.89 × 106 m³/year) from the WT scheme did not compensate for the decrease in groundwater head in areas that do not surround the receiving streams. Actually, the additional recharge is counteracted by a decrease of 25% of natural groundwater recharge or an increase of 50% in pumping rate; therefore, WT transmission losses alone would not solve the issue of the unsustainable management of groundwater resources.