Frequency‐domain acoustic‐elastic coupled waveform inversion using the Gauss‐Newton conjugate gradient method

We developed a frequency‐domain acoustic‐elastic coupled waveform inversion based on the Gauss‐Newton conjugate gradient method. Despite the use of a high‐performance computer system and a state‐of‐the‐art parallel computation algorithm, it remained computationally prohibitive to calculate the approximate Hessian explicitly for a large‐scale inverse problem. Therefore, we adopted the conjugate gradient least‐squares algorithm, which is frequently used for geophysical inverse problems, to implement the Gauss‐Newton method so that the approximate Hessian is calculated implicitly. Thus, there was no need to store the Hessian matrix. By simultaneously back‐propagating multi‐components consisting of the pressure and displacements, we could efficiently extract information on the subsurface structures. To verify our algorithm, we applied it to synthetic data sets generated from the Marmousi‐2 model and the modified SEG/EAGE salt model. We also extended our algorithm to the ocean‐bottom cable environment and verified it using ocean‐bottom cable data generated from the Marmousi‐2 model. With the assumption of a hard seafloor, we recovered both the P‐wave velocity of complicated subsurface structures as well as the S‐wave velocity. Although the inversion of the S‐wave velocity is not feasible for the high Poisson's ratios used to simulate a soft seafloor, several strategies exist to treat this problem. Our example using multi‐component data showed some promise in mitigating the soft seafloor effect. However, this issue still remains open.     

Thermal anomalies associated with shallow gas hydrates in the K-2 mud volcano, Lake Baikal

Thermal measurements and hydrate mapping in the vicinity of the K-2 mud volcano in Lake Baikal have revealed a particular type of association of thermal anomalies (29–121?mW?m^–2) near hydrate-forming layers. Detailed coring within K-2 showed that hydrates are restricted to two distinct zones at sub-bottom depths exceeding 70–300?cm. Temperature data from stations with hydrate recovery and degassing features all display low thermal gradients. Otherwise, the thermal gradients within the mud volcano are generally increased. These findings imply a more complicated thermal regime than often assumed for mud volcanoes, with important roles for both fluids and hydrates. The coexistence of neighbouring low and high thermal anomalies is interpreted to result from discharging and recharging fluid activity, rather than hydrate thermodynamics. It is suggested that hydrates play a key role in controlling the fluid circulation pattern at an early stage. At a later stage, the inflow of undersaturated lake water would favour the dissolution of structure I hydrates and the formation of structure II hydrates, the latter having been observed on top of structure I hydrates in the K-2 mud volcano.     

Reduction of future disaster damages by learning from disaster experiences

This paper examines the effect of a country’s own past disaster experiences and nearby countries’ past experiences on subsequent disaster damage. We use global disaster data from 1990 to 2010, which include disaster-related death tolls for both natural and technological disasters that are further divided into sub-categories. Overall, we find evidence of a reduction effect of past disaster damage on future disaster damage. More detailed analyses show that an adaptation effect seems to be present for certain combinations of disaster types and levels of economic development. The results show that a country’s own experiences reduce future damage for natural disasters but that the marginal effect is larger for lower-income countries. On the other hand, for technological disasters, a robust impact of experiences was found only in higher-income countries. In terms of the disaster experiences of nearby countries, the adaptation effect was found only for natural disasters, and the marginal impact was relatively higher for higher-income countries.     

Acoustic full waveform inversion of synthetic land and marine data in the Laplace domain

Elastic waves, such as Rayleigh and mode‐converted waves, together with amplitude versus offset variations, serve as noise in full waveform inversion using the acoustic approximation. Heavy preprocessing must be applied to remove elastic effects to invert land or marine data using the acoustic inversion method in the time or frequency domains. Full waveform inversion using the elastic wave equation should be one alternative; however, multi‐parameter inversion is expensive and sensitive to the starting velocity model. We implement full acoustic waveform inversion of synthetic land and marine data in the Laplace domain with minimum preprocessing (i.e., muting) to remove elastic effects. The damping in the Laplace transform can be thought of as an automatic time windowing. Numerical examples show that Laplace‐domain acoustic inversion can yield correct smooth velocity models even with the noise originating from elastic waves. This offers the opportunity to develop an accurate smooth starting model for subsequent inversion in the frequency domain.     

Organohalogen compounds in deep-sea fishes from the western North Pacific, off-Tohoku, Japan: Contamination status and bioaccumulation profiles

Twelve species of deep-sea fishes collected in 2005 from the western North Pacific, off-Tohoku, Japan were analyzed for organohalogen compounds. Among the compounds analyzed, concentrations of DDTs and PCBs (up to 23,000 and 12,400 ng/g lipid wt, respectively) were the highest. The present study is the foremost to report the occurrence of brominated flame retardants such as PBDEs and HBCDs in deep-sea organisms from the North Pacific region. Significant positive correlations found between δ¹⁵N (‰) and PCBs, DDTs and PBDEs suggest the high biomagnification potential of these contaminants in food web. The large variation in δ¹³C (‰) values observed between the species indicate multiple sources of carbon in the food web and specific accumulation of hydrophobic organohalogen compounds in benthic dwelling carnivore species like snubnosed eel. The results obtained in this study highlight the usefulness of deep-sea fishes as sentinel species to monitor the deep-sea environment.     

Fatty acid profiles of benthic environment associated with artificial reefs in subtropical Hong Kong

Artificial reefs can enhance habitat heterogeneity, especially in seabed degraded by bottom-dredging and trawling. However, the trophodynamics of such reef systems are not well understood. This study provided baseline data on trophic relationships in the benthic environment associated with artificial reefs in late spring and mid summer of subtropical Hong Kong, using fatty acid profiles as an indicator. Data from sediments collected at the reef base, materials from sediment traps deployed on top and bottom of the reefs, total particulate matter from the water column and oyster tissues from reef surface were subjected to principal component analysis. Results showed variations of fatty acid profiles in the total particulate matter, upper sediment trap and oyster tissue samples collected in the two samplings, indicating seasonal, trophodynamic changes within the reef system. The wastes produced by fish aggregating at the reefs can also contribute a source of biodeposits to the nearby benthic environment.     

Reconstruction of historical nutrient levels in Korean and Japanese coastal areas based on dinoflagellate cyst assemblages

Dinoflagellate cysts acquired from sediment cores were analyzed in order to reconstruct historical nutrient levels in Gamak Bay, Korea and Ariake Bay, Japan. Dinoflagellate cyst assemblages in Gamak Bay were characterized by high proportions of heterotrophic cysts such as Brigantedinium spp., Protoperidinium americanum and Polykrikos cysts, which suggested that nutrients levels may have already been high before 1970s, and then increased further to the hypertrophic conditions of the 1990s. In contrast, dinoflagellate cyst assemblages in Ariake Bay were characterized by high relative abundances of Lingulodinium machaerophorum and Spiniferites spp., which suggested that nutrient levels in Ariake Bay had increased gradually since the mid 1960s, and may have been significantly enhanced by the mid 1980s. Dinoflagellate cyst assemblages reflecting environmental changes in the two bays are contrasting, perhaps due to different nutrient enrichment mechanisms. This suggests that the indicators of nutrient levels encoded in dinoflagellate cyst assemblages may exhibit site-specific information.     

Electrokinetic Removal of Petroleum Hydrocarbon from Residual Clayey Soil Following a Washing Process

This study investigates total petroleum hydrocarbon (TPH) removal from residual clayey soil, after a washing procedure, using an electrokinetic process. Eight electrokinetic experiments were carried out to investigate the characteristics of TPH removal. When 0.1 M MgSO₄ or 0.1 M NaOH was used as an electrolyte, the electric current rapidly increased within the first 100 or 200 h, respectively. A negatively charged soil surface resulted in a more negative zeta potential and greater electroosmotic flow toward the cathode. Therefore, the accumulated electroosmotic flow (EOF) when using 0.1 M NaOH as the anolyte‐purging solution was higher than when using 0.1 M MgSO₄. Although the energy consumption for the two purging solutions was similar, the efficiencies of TPH removal when 0.1 M MgSO₄ and 0.1 M NaOH with surfactant were used were 0 and 39%, respectively, because the electroosmotic flow rate increased with TPH removal efficiency. When 5% isopropyl alcohol (IPA) was used as a circulation solution, the electric current increased but the TPH removal was similar to that using water. In terms of energy consumption, the use of a surfactant‐enhanced electrokinetic process with NaOH as electrolyte was effective in removing TPHs from low‐permeability soil.     

2D Laplace-Domain Waveform Inversion of Field Data Using a Power Objective Function

The wavefield in the Laplace domain has a very small amplitude except only near the source point. In order to deal with this characteristic, the logarithmic objective function has been used in many Laplace domain inversion studies. The Laplace-domain waveform inversion using the logarithmic objective function has fewer local minima than the time- or frequency domain inversion. Recently, the power objective function was suggested as an alternative to the logarithmic objective function in the Laplace domain. Since amplitudes of wavefields are very small generally, a power <1 amplifies the wavefields especially at large offset. Therefore, the power objective function can enhance the Laplace-domain inversion results. In previous studies about synthetic datasets, it is confirmed that the inversion using a power objective function shows a similar result when compared with the inversion using a logarithmic objective function. In this paper, we apply an inversion algorithm using a power objective function to field datasets. We perform the waveform inversion using the power objective function and compare the result obtained by the logarithmic objective function. The Gulf of Mexico dataset is used for the comparison. When we use a power objective function in the inversion algorithm, it is important to choose the appropriate exponent. By testing the various exponents, we can select the range of the exponent from 5 × 10^?3 to 5 × 10^?8 in the Gulf of Mexico dataset. The results obtained from the power objective function with appropriate exponent are very similar to the results of the logarithmic objective function. Even though we do not get better results than the conventional method, we can confirm the possibility of applying the power objective function for field data. In addition, the power objective function shows good results in spite of little difference in the amplitude of the wavefield. Based on these results, we can expect that the power objective function will produce good results from the data with a small amplitude difference. Also, it can partially be utilized at the sections where the amplitude difference is very small.