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.     

Assessment on variability of extreme climate events for the Upper Thames River basin in Canada

Climate change may affect magnitude and frequency of regional extreme events with possibility of serious impacts on the existing infrastructure systems. This study investigates how the current spatial and temporal variations of extreme events are affected by climate change in the Upper Thames River basin, Ontario, Canada. A weather generator model is implemented to obtain daily time series of three climate variables for two future climate scenarios. The daily time series are disaggregated into hourly to capture characteristics of intense and rapidly changing storms. The maximum annual precipitation events for five short durations, 6‐, 12‐, 24‐, 48‐, and 72‐h durations, at each station are extracted from the generated hourly data. The frequency and seasonality analyses are conducted to investigate the temporal and spatial variability of extreme precipitation events corresponding to each duration. In addition, this study investigates the impacts of increase in temperature using reliability, resilience, and vulnerability. The results indicate that the extreme precipitation events under climate change will occur earlier than in the past. In addition, episodes of extremely high temperature may last longer up to 19·7% than under the no‐change climate scenario. This study points out that the revision of the design storms (e.g. 100‐ or 250‐year return period) is warranted for the west and the south east region of the basin. Copyright © 2011 John Wiley & Sons, Ltd.     

Estimation of methane fluxes from bottom sediments of Lake Baikal

Methane bubble fluxes in gas flares from bottom sediments in Lake Baikal were estimated for the first time using hydroacoustic methods. Earlier work has demonstrated the occurrence of gas seeps both inside and outside of areas where bottom simulating reflectors were identified in seismic profiles. Fluxes ranged from 14 to 216 tons per year, with the flux for the entire area of the central and southern basins ranging from 1,400 to 2,800 tons per year. Comparison with other water bodies showed that fluxes from the most intensive Baikal flares were similar to those in the Norwegian and Okhotsk seas. Gas hydrates decompose at the lower boundary of the gas hydrate stability zone due to sedimentation. Calculation of the amount of methane produced due to sedimentation gave a total of between 2,600 and 14,000 tons per year for the central and southern basins of the lake. Based on rough estimation, the total flux from shallow- and deep-water gas seeps is similar to the amount of methane produced due to sedimentation. This suggests that gas hydrates possibly occupy much more than 10?% of the pore volume near the base of the gas hydrate stability zone, or that there are other reasons for gas hydrate dissociation and bubble flux from these bottom sediments.     

Geochemical characteristics and evaluation of minor and trace elements pollution in groundwater of Tuticorin city,Tamil Nadu,India using geospatial techniques

In this work the major and trace elements pollution levels in groundwater of Tuticorin city of Tamil Nadu was investigated. This coastal city with numerous large and small scale industries has increasingly affected the water quality by dumping effluents. The mean major and trace elements levels in groundwater in descending order: Si > Sr > Zn > Fe > Rb > Se > Al > Mn > Cu > As > Cr > Ba > V > Ni > Pb > Mo > Co > Sb > Cd > Ag. The concentrations of As, Se, Pb, Ba, Al, Fe and V (0.084, 0.150, 0.057, 0.837, 0.357 and 0.052) in some of the groundwater samples exceed the standard value of world health organization (WHO, 2004). Interrelationship between the elements was done by correlation matrix and factor analysis which indicates the contribution from industrial and irrigation fields. Spatial interpolation of inverse distance method (IDW) maps of all the elements suggested that discharge of effluents from chemical factory (STERLITE, TAC, SPIC, HEAVY WATER PLANT, NILA SEA FOODS), municipal wastewater, fertilizers were added to the natural sources. Results of this research suggests that proper measures should be taken to protect the vital groundwater resources in the Tuticorin city.     

Mid-Miocene (post 12 Ma) displacement along the central Karakoram fault zone in the Nubra Valley,Ladakh, India from spot LA-ICPMS U/Pb zircon ages of granites

The Karakoram fault zone is a prominent right lateral fault that connects the frontal thrust of the North Pamir with the Indus suture zone near Mount Kailas. Its nature and age of initiation is controversial. In the Nubra valley, Ladakh, India, a Karakoram range granite is thrust over Cretaceous magmatic arc rocks and this thrust is cut by a western strand of the Karakoram fault zone. Three different lithologies from this granite gave weighted mean zircon U/Pb ages of 12.92±0.77 Ma, 12.41±0.43 Ma, and 11.72±0.31 Ma. The ages indicate a relatively short intrusive history of about 1 Ma for the phases: the geochemistry is practically identical to the Pangong leucogranites in the same tectonic block. The Karakoram fault zone in this area is thus less than ~12 Ma old which supports a post middle Miocene (Serravallian) age of Karakoram fault initiation in this area.     

Subduction-related 200 Ma Talun metagranite,SE Taiwan: an age constraint for palaeo-Pacific plate subduction beneath South China Block during the Mesozoic

Various tectonic models have been proposed to account for the widely distributed igneous activities in the southeastern part of the South China Block (SCB) during the Triassic–Jurassic period. One of the major contending debates is on the timing of initiation of the palaeo-Pacific plate subduction under the SCB, due to lack of unequivocal evidence for arc magmatism during the period in this region.

The 191 ± 10 Ma (N = 5, MSWD = 12) calc-alkalic high-K I-type Talun metagranite occurs in the southern Tailuko belt of the Tananao metamorphic complex, Taiwan. In terms of age, this metagranite belongs to the Early Yanshanian igneous activity in the southeastern part of the SCB. However, its geographic position does not accord with the well-known general oceanward younging trend of the Yansnanian igneous rocks. In view of the large age uncertainty reported, this metagranite is redated in this study. Some zircons of this metagranite are high in U content and are metamict. Zircons with low U contents are analysed by SHRIMP yielding a more precise age of 200 ± 2 Ma (N = 10, MSWD = 4). In particular, the εHf(t) of these dated zircons ranges from +4.5 to +12.9. The metagranite mainly consists of quartz, K-feldspar, plagioclase, with minor amounts of garnet, biotite, zircon, apatite, and pyrrhotite. Chlorite and calcite are secondary phases overprinted by the later tectonic event(s). Its initial Sr isotope compositional range is 0.70473–0.70588, and εNd(t), +2.4 to +3.6. The results demonstrate that the genesis of this metagranite could be attributed to the assimilation-fractionation of a depleted mantle-derived basaltic magma, which was most likely related to arc magmatism. The present study therefore offers key evidence that during the Mesozoic, the palaeo-Pacific plate subduction underneath the SCB would have taken place no later than the very early Jurassic.     

Selecting suitable sites for mountain ginseng (<Emphasis Type="Italic">Panax ginseng</Emphasis>) cultivation by using geographically weighted logistic regression

With the well-being trends to pursue a healthy life, mountain ginseng (Panax ginseng) is rising as one of the most profitable forest products in South Korea. This study was aimed at evaluating a new methodology for identifying suitable sites for mountain ginseng cultivation in the country. Forest vegetation data were collected from 46 sites and the spatial distribution of all sites was analyzed using GIS data for topographic position, landform, solar radiation, and topographic wetness. The physical and chemical properties of the soil samples, including moisture content, pH, organic matter, total nitrogen, exchangeable cations, available phosphorous, and soil texture, were analyzed. The cultivation suitability at each site was assessed based on the environmental conditions using logistic regression (LR) and geographically weighted logistic regression (GWLR) and the results of both methods were compared. The results show that the areas with northern aspect and higher levels of solar radiation, moisture content, total nitrogen, and sand ratio are more likely to be identified as suitable sites for ginseng cultivation. In contrast to the LR, the spatial modeling with the GWLR results in an increase in the model fitness and indicates that a significant portion of spatial autocorrelation in the data decreases. A higher value of the area under the receiver operating characteristic (ROC) curve presents a better prediction accuracy of site suitability by the GWLR. The geographically weighted coefficient estimates of the model are non-stationary, and reveal that different site suitability is associated with the geographical location of the forest stands. The GWLR increases the accuracy of selecting suitable sites by considering the geographical variations in the characteristics of the cultivation sites.     

Long-term seismological sea-bottom monitoring using autonomous bottom stations

Merits and demerits of recording of seismic signals at the bottom of water areas are considered. It is shown that long-term seismological monitoring systems should be placed in the regions of industrial development of the shelf and continental slope and in the areas of high seismic and tsunami hazard of oceans and seas. The results obtained during expeditions of the Shirshov Institute of Oceanology of the Russian Academy of Sciences with the use of broadband bottom seismographs are reported. Autonomous bottom seismographs with long-term operation at the bottom and operative communication via satellite and radio channels are proposed for the formation of a marine seismological network.     

Soft‐sediment deformation structures in Cretaceous non‐marine deposits of southeastern Gyeongsang Basin,Korea: Occurrences and origin

A variety of soft‐sediment deformation structures formed during or shortly after deposition occurs in the Cretaceous Seongpori and Dadaepo Formations of the southeastern Gyeongsang Basin exposed along coastal areas of southeastern Korean Peninsula for 0.5–2 km. These are mostly present in a fluvial plain facies, with interbedded lacustrine deposits. In this study, the features of different kinds of soft‐sediment deformation structures have been interpreted on the basis of sedimentology of structure‐bearing deposits, comparison with normal sedimentary structures, timing and mechanism of deformation, and triggering mechanisms. The soft‐sediment deformation structures can be classified into four morphological groups: (i) load structures (load casts, ball‐and‐pillow structures); (ii) soft‐sediment intrusive structures (dish‐and‐pillars, clastic dykes, sills); (iii) ductile disturbed structures (convolute folds, slump structures); and (iv) brittle deformation structures (syndepositional faulting, dislocated breccia). The most probable triggering mechanisms resulting in these structures were seismic shocks. These interpretations are based on the following field observations: (i) location of the study area within tectonically active fault zone reactivated several times during the Cretaceous; (ii) deformation structures confined to single stratigraphic levels; (iii) lateral continuity and occurrences of various soft‐sediment deformation structures in the deformed level over large areas; (iv) absence of depositional slope to indicate gravity sliding or slumping; and (v) similarity to the structures produced experimentally. The soft‐sediment deformation structures in the study areas are thus interpreted to have been generated by seismic shocks with an estimated magnitude of M > 5, representing an intermittent record of the active tectonic and sedimentary processes during the development and evolution of two formations from the late Early Cretaceous to the Late Cretaceous.     

Application of bivariate frequency analysis to the derivation of rainfall–frequency curves

Bivariate distributions have been recently employed in hydrologic frequency analysis to analyze the joint probabilistic characteristics of multivariate storm events. This study aims to derive practical solutions of application for the bivariate distribution to estimate design rainfalls corresponding to the desired return periods. Using the Gumbel mixed model, this study constructed rainfall–frequency curves at sample stations in Korea which provide joint relationships between amount, duration, and frequency of storm events. Based on comparisons and analyses of the rainfall–frequency curves derived from univariate and bivariate storm frequency analyses, this study found that conditional frequency analysis provides more appropriate estimates of design rainfalls as it more accurately represents the natural relationship between storm properties than the conventional univariate storm frequency analysis.