Recovery of groundwater in the Sanriku region contaminated by the tsunami inundation from the 2011 Tohoku earthquake

Many wells in the Sanriku region used as sources for water supply systems were heavily contaminated by the tsunami of the 2011 great Tohoku earthquake on March 11 in 2011. To better understand the nature of the groundwater contamination by the tsunami inundation and to clarify the recovery process of contaminated groundwater at the study wells, groundwater monitoring has been conducted once or twice yearly since early summer in 2011. High and abnormal values of electric conductivity (EC), chloride ion concentration (CIC), Na⁺, Ca⁺, heavy metal ions, and heavier isotopes of the contaminated groundwater were also obtained in April and June 2011. The chemical elements have rapidly and exponentially decreased as a result of effective pumping of the contaminated groundwater from the study wells and because of abundant rainfall in 2011. In April 2015 (about 4 years after the tsunami inundation), the CIC and EC of the contaminated groundwater of two study wells in Minamisanriku town had reached pre-inundation values. The estimated residence times of groundwater of the two study wells were 105–118 days in the full-day pumping stage and 910–1000 days in the daytime-only pumping stage.     

Evaluation of site effects on ground motions based on equivalent linear site response analysis and liquefaction potential in Chennai,south India

We study local site effects with detailed geotechnical and geophysical site characterization to evaluate the site-specific seismic hazard for the seismic microzonation of the Chennai city in South India. A Maximum Credible Earthquake (MCE) of magnitude 6.0 is considered based on the available seismotectonic and geological information of the study area. We synthesized strong ground motion records for this target event using stochastic finite-fault technique, based on a dynamic corner frequency approach, at different sites in the city, with the model parameters for the source, site, and path (attenuation) most appropriately selected for this region. We tested the influence of several model parameters on the characteristics of ground motion through simulations and found that stress drop largely influences both the amplitude and frequency of ground motion. To minimize its influence, we estimated stress drop after finite bandwidth correction, as expected from an M6 earthquake in Indian peninsula shield for accurately predicting the level of ground motion. Estimates of shear wave velocity averaged over the top 30 m of soil (V_S30) are obtained from multichannel analysis of surface wave (MASW) at 210 sites at depths of 30 to 60 m below the ground surface. Using these V_S30 values, along with the available geotechnical information and synthetic ground motion database obtained, equivalent linear one-dimensional site response analysis that approximates the nonlinear soil behavior within the linear analysis framework was performed using the computer program SHAKE2000. Fundamental natural frequency, Peak Ground Acceleration (PGA) at surface and rock levels, response spectrum at surface level for different damping coefficients, and amplification factors are presented at different sites of the city. Liquefaction study was done based on the V_S30 and PGA values obtained. The major findings suggest show that the northeast part of the city is characterized by (i) low V_S30 values (<?200 m/s) associated with alluvial deposits, (ii) relatively high PGA value, at the surface, of about 0.24 g, and (iii) factor of safety and liquefaction below unity at three sites (no. 12, no. 37, and no. 70). Thus, this part of the city is expected to experience damage for the expected M6 target event.     

Experimental study of headcut erosion in cohesive soils under different consolidation types and hydraulic parameters

Headcut is a change in stream channel elevation, where there is concentrated flow. Most of the past studies focus on non-cohesive soils, although many problems on the streams occur because of cohesive beds and banks. In this study, eight samples of cohesive soils, with a different composition of silt and clay, for different waterfall heights and flow velocity under long- and short-term natural consolidation conditions were tested. In one of the tests, a sand layer was settled on the headcut bed to investigate its impact on headcut erosion. By increasing clay content, the headcut will remain vertical as it moves backward. Result showed that the effect of clay content reduction was more noticeable under the long-term consolidation condition. In general, the effect of clay percentage variation on the measured parameters is much higher than the effects of waterfall height or flow velocity, and the effect of consolidation type and adding a sand layer on the measured values is much higher than the effect of clay percentage variation on the waterfall height and flow velocity.     

New classes of spectral densities for lattice processes and random fields built from simple univariate margins

Quasi arithmetic and Archimedean functionals are used to build new classes of spectral densities for processes defined on any d-dimensional lattice \mathbbZ^d{\mathbb{Z}^d} and random fields defined on the d-dimensional Euclidean space \mathbbR^d{\mathbb{R}^d}, given simple margins. We discuss the mathematical features of the proposed constructions, and show rigorously as well as through examples, that these new classes of spectra generalize celebrated classes introduced in the literature. Additionally, we obtain permissible spectral densities as linear combinations of quasi arithmetic or Archimedean functionals, whose associated correlation functions may attain negative values or oscillate between positive and negative ones. We finally show that these new classes of spectral densities can be used for nonseparable processes that are not necessarily diagonally symmetric.     

From global to local: Testing the potential of cross-scaling in global data sets

The present study investigates the potential of readily available and easily accessible global data sets to understand regional/local level interactions in wetland systems. The biogeographical zones of India were used a base-frame to select three sites. The study well fits the interests of National Wetland Committee of India to investigate and document fundamental information on wetland extent/distribution. The national partnership with SACON represents this interest. SACON commenced the inland wetland inventory module at national scale using geospatial data, although the provincial scale analysis is underway. In addition, the global irrigated area mapping (GIAM-IWMI) project generated multi-scalar spatial outputs for irrigated/rain-fed areas. With the existing information base, a multi-level geospatial analysis using Arc GIS algorithmic modelling was used to derive comprehensive appraisal of wetland systems complementing the data from GIAM and SACON. It was observed that the overlap between the two layers was 58 percent for Gujarat and 10 percent in Tamil Nadu. In Krishna basin the wetland’s cover 1.04 million hectare excluding the rice agro-ecosystem. The difference in the biogeography of the case sites governs the gradient of information derived from both data layers. Additionally, the global lakes and wetlands database (GLWD) database added thematic information on coastal wetlands. In summary we describe the cross-scaling the global data layers to compliment the regional/national level monitoring assignments.     

Simultaneous estimation of cross-validation errors in least squares collocation applied for statistical testing and evaluation of the noise variance components

The cross-validation technique is a popular method to assess and improve the quality of prediction by least squares collocation (LSC). We present a formula for direct estimation of the vector of cross-validation errors (CVEs) in LSC which is much faster than element-wise CVE computation. We show that a quadratic form of CVEs follows Chi-squared distribution. Furthermore, a posteriori noise variance factor is derived by the quadratic form of CVEs. In order to detect blunders in the observations, estimated standardized CVE is proposed as the test statistic which can be applied when noise variances are known or unknown. We use LSC together with the methods proposed in this research for interpolation of crustal subsidence in the northern coast of the Gulf of Mexico. The results show that after detection and removing outliers, the root mean square (RMS) of CVEs and estimated noise standard deviation are reduced about 51 and 59%, respectively. In addition, RMS of LSC prediction error at data points and RMS of estimated noise of observations are decreased by 39 and 67%, respectively. However, RMS of LSC prediction error on a regular grid of interpolation points covering the area is only reduced about 4% which is a consequence of sparse distribution of data points for this case study. The influence of gross errors on LSC prediction results is also investigated by lower cutoff CVEs. It is indicated that after elimination of outliers, RMS of this type of errors is also reduced by 19.5% for a 5 km radius of vicinity. We propose a method using standardized CVEs for classification of dataset into three groups with presumed different noise variances. The noise variance components for each of the groups are estimated using restricted maximum-likelihood method via Fisher scoring technique. Finally, LSC assessment measures were computed for the estimated heterogeneous noise variance model and compared with those of the homogeneous model. The advantage of the proposed method is the reduction in estimated noise levels for those groups with the fewer number of noisy data points.     

Temporal variability in summertime bottom hypoxia in shallow areas of Mobile Bay,Alabama

This paper addresses temporal variability in bottom hypoxia in broad shallow areas of Mobile Bay, Alabama. Time-series data collected in the summer of 2004 from one station (mean depth of 4 m) exhibit bottom dissolved oxygen (DO) variations associated with various time scales of hours to days. Despite a large velocity shear, stratification was strong enough to suppress vertical mixing most of the time. Bottom DO was closely related to the vertical salinity gradient (ΔS). Hypoxia seldom occurred when ΔS (over 2.5 m) was <2 psu and occurred almost all the time when ΔS was >8 psu in the absence of extreme events like hurricanes. Oxygen balance between vertical mixing and total oxygen demand was considered for bottom water from which oxygen demand and diffusive oxygen flux were estimated. The estimated decay rates at 20°C ranging between 0.175–0.322 d⁻¹ and the corresponding oxygen consumption as large as 7.4 g O₂ m⁻² d⁻¹ fall at the upper limit of previously reported ranges. The diffusive oxygen flux and the corresponding vertical diffusivity estimated for well mixed conditions range between 8.6–9.5 g O₂ m⁻² d⁻¹ and 2.6–2.9 m² d⁻¹, respectively. Mobile Bay hypoxia is likely to be associated with a large oxygen demand, supported by both water column and sediment oxygen demands, so that oxygen supply from surface water during destratification events would be quickly exhausted to return to hypoxic conditions within a few hours to days after destratification events are terminated.     

Geochemical characteristics of island-arc volcanic rocks in the Nan-Nam Pat-Phetchabun zone,northern Thailand

Late Permian-Early Triassic (P₂-T₁) volcanic rocks distributed on the eastern side of ocean-ridge and oceanic-island basalts in the Nan-Uttaradit zone were analyzed from aspects of petrographic characteristics, rock assemblage, REE, trace elements, geotectonic setting, etc., indicating that those volcanic rocks possess the characteristic features of island-arc volcanic rocks. The volcanic rock assemblage is basalt-basaltic andesite-andesite. The volcanic rocks are sub-alkaline, dominated by calc-alkaline series, with tholeiite series coming next. The chemical composition of the volcanic rocks is characterized by low TiO₂ and K₂O and high Al₂O₃ and Na₂O. Their REE patterns are of the flat, weak LREE-enrichment right-inclined type. The trace elements are characterized by the enrichment of large cation elements such as K, Rb and Ba, common enrichment of U and Th, and depletion of Nb, Ta, Zr and Hf. The petrochemical plot falls within the field of volcanic rocks, in consistency with the plot of island-arc volcanic rocks in the Jinsha River zone of China. This island-arc volcanic zone, together with the ocean-ridge/oceanic island type volcanic rocks in the Nan-Uttaradit zone, constitutes the ocean-ridge volcanic rock-island-arc magmatic rock zones which are distributed in pairs, indicating that the oceanic crust of the Nan-Uttaradit zone once was of eastward subduction. This work is of great significance in exploring the evolution of paleo-Tethys in the Nan-Uttaradit zone.     

Assessment of three temperature reconstruction methods in the virtual reality of a climate simulation

The performance of statistical climate reconstruction methods in the pre-instrumental period is uncertain, as they are calibrated in a short instrumental period but applied to much longer reconstructions time spans. Here, the virtual reality created by a climate simulation of the past millennium with the model ECHO-G is used as a test bed of three methods to reconstruct the annual Northern Hemisphere temperature. The methods are Composite plus Scaling, the inverse regression method of Mann et al. (Nature 392:779–787, 1998) and a direct principal-components regression method. The testing methodology is based on the construction of pseudo-proxies derived from the climate model output, the application of each of these methods to pseudo-proxy timeseries, and the comparison of their result with the simulated mean temperature. Different structures of the noise have been used to construct pseudo-proxies, ranging from the simulated grid-point precipitation. Also, one sparse and one denser pseudo-proxy network, co-located with two real networks, have been considered. All three methods underestimate the simulated variations of the Northern Hemisphere temperature, but the Composite plus Scaling method clearly displays a better performance and is robust against the different noise models and network size. The most relevant factor determining the skill of the reconstruction appears to be the network size, whereas the different noise models tend to yield similar results.     

Segmentation algorithm for long time series analysis

Time series analysis is an important issue in the earth science-related engineering applications such as hydrology, meteorology and environmetrics. Inconsistency and nonhomogeneity that might arise in a time series yield segments with different statistical characteristics. In this study, an algorithm based on the first order statistical moment (average) of a time series is developed and applied on five time series with length ranging from 84 items to nearly 1,300. Comparison to the existing segmentation algorithms proves the applicability and usefulness of the proposed algorithm in long hydrometeorological and geophysical time series analysis.