Coseismic responses and the mechanism behind MW 5.1 earthquake of March 14, 2005 in the Koyna–Warna region, India

An earthquake of M_w 5.1 occurred on March 14, 2005, in the seismically active Koyna–Warna region in western India, the site known for the largest reservoir triggered seismicity (RTS) in the world. For more than four decades, earthquakes with M  4.0 have occurred in this region at regular intervals. Impoundment of reservoirs and changes in lake levels can trigger earthquakes by two processes of stress modifications, namely direct loading effect of the reservoir and diffusion through various faults and fractures. In this paper we analysed the reservoir water level data at Koyna and Warna reservoirs prior to the occurrence of the March 14, 2005 earthquake, to explain the dominant mechanism behind its occurrence and its correlation with the observed coseismic changes. We conclude that the diffusion process, not the reservoir load effect, is the dominating mechanism triggering earthquakes in the region. The coseismic changes in deep well water levels sensitive to earth tides are found to be to the order of 1–12 cm.     

Uncertainty as the overlap of alternate conditional distributions

An important task in modern geostatistics is the assessment and quantification of resource and reserve uncertainty. This uncertainty is valuable support information for many management decisions. Uncertainty at specific locations and uncertainty in the global resource is of interest. There are many different methods to build models of uncertainty, including Kriging, Cokriging, and Inverse Distance. Each method leads to different results. A method is proposed to combine local uncertainties predicted by different models to obtain a combined measure of uncertainty that combines good features of each alternative. The new estimator is the overlap of alternate conditional distributions.     

Accurate local upscaling with variable compact multipoint transmissibility calculations

We propose a new single-phase local upscaling method that uses spatially varying multipoint transmissibility calculations. The method is demonstrated on two-dimensional Cartesian and adaptive Cartesian grids. For each cell face in the coarse upscaled grid, we create a local fine grid region surrounding the face on which we solve two generic local flow problems. The multipoint stencils used to calculate the fluxes across coarse grid cell faces involve the six neighboring pressure values. They are required to honor the two generic flow problems. The remaining degrees of freedom are used to maximize compactness and to ensure that the flux approximation is as close as possible to being two-point. The resulting multipoint flux approximations are spatially varying (a subset of the six neighbors is adaptively chosen) and reduce to two-point expressions in cases without full-tensor anisotropy. Numerical tests show that the method significantly improves upscaling accuracy as compared to commonly used local methods and also compares favorably with a local–global upscaling method.     

Geochemistry of amphibolitized eclogites and cross-cutting tonalitic–trondhjemitic dykes in the Metamorphic Kimi Complex in East Rhodope (N.E. Greece): implications for partial melting at the base of a thickened crust

In the ultra-high pressure Metamorphic Kimi Complex widespread tonalitic–trondhjemitic dykes, with an intrusion age ca. 65–63 Ma, cross-cut boudins and layers of amphibolitized eclogites. Geochemical investigation proclaims the tied genetic relationship of the amphibolitized eclogites and the associated tonalitic–trondhjemitic dykes. The major and trace element contents and rare earth element patterns of the amphibolitized eclogites indicate formation of their protoliths by fractional crystallization of tholeiitic magmas in a back-arc environment. The tonalites and trondhjemites are characterized by moderate to high Sr contents (>130 ppm), and low Y (<8.2 ppm) and heavy rare earth element contents (Yb content of 0.19–0.88 ppm). The chemical composition of the tonalitic and trondhjemitic dykes are best explained by partial melting of a tholeiitic source like the amphibolitized eclogites with residual garnet and amphibole, at the base of a thickened crust during Early Tertiary subduction/accretion at the southern margins of the European continent.     

Surface water–groundwater interaction and chemistry in a mineral-armored hydrothermal outflow channel, Yellowstone National Park, USA

Small quantities of groundwater interact with hydrothermal surface water to drive in-stream geochemical processes in a silica-armored hot-spring outflow channel in Yellowstone National Park, USA. The objective of this study was to characterize the hydrology and geochemistry of this unique system in order to (1) learn more about the Yellowstone Plateau’s subsurface water mixing between meteoric and hydrothermal waters and (2) learn more about the chemical and physical processes that lead to accumulation of streambed cements, i.e., streambed armor. A combination of hydrological, geochemical, mineralogical, microscopic, and petrographic techniques were used to identify groundwater and surface-water exchange. Interaction could be identified in winter because of differences in surface water and groundwater composition but interaction at other times of the year cannot be ruled out. Dissolved constituents originating from groundwater (e.g., Fe(II) and Mg) were traced downstream until oxidation and/or subsequent precipitation with silica removed them, particularly where high affinity substrates like cyanobacterial surfaces were present. Because the stream lies in a relatively flat drainage basin and is fed mainly by a seasonally relatively stable hot spring, this system allowed study of the chemical processes along a stream without the obscuring effects of sedimentation.     

Field scale heterogeneity of redox conditions in till-upscaling to a catchment nitrate model

Point scale studies in different settings of glacial geology show a large local variation of redox conditions. There is a need to develop an upscaling methodology for catchment scale models. This paper describes a study of field-scale heterogeneity of redox-interfaces in a till aquitard within an area of 600?×?600 m. The results showed significant variation of the depths to the redox-interface and thickness of the aquitard. Nitrate was present above the redox-interface but reduced to non-detectable levels a few metres below the interface. An upscaling approach for an area of 92 km² is proposed. Two models are proposed to predict the depth to the redox-interface in the aquitard and the resulting nitrate recharge concentrations to an underlying aquifer. The first model assumes that the depth to the redox-interface reflects the hydraulic head in the aquitard, and the second model assumes that the depth of the redox-interface is randomly distributed according to a log-normal probability distribution function. The upscaling approach using the random redox model estimated recharge concentrations comparable to the observed concentration in the underlying aquifer. The presented upscaling approach is applicable in distributed catchment models where sub-grid variability cannot be represented by the large grids.     

Sediment Organic Carbon in Todos Santos Bay, Baja California, Mexico

Sediment grain size and organic carbon (OC) data collected over the past 50 years, together with δ¹³C values of OC in recently collected samples, were analyzed to improve understanding of sediment OC distribution and abundance in Todos Santos Bay. Sediments in the submarine canyon at the mouth of the bay and in quiet-water locations along the shore are fine grained, high in OC, and have generally low δ¹³C values; sediments in high-energy environments are low in OC and have high δ¹³C values. A bivariate isotopic mixing model indicates that none of the sediments contain >50% terrigenous OC (average ~30%), and that the terrigenous OC content of the sediments is a small proportion of the OC content of local soils. Sediment OC composition is apparently controlled by energy-related sorting and deposition, oxidation of much of the original terrigenous OC, and replacement of some terrigenous OC by marine OC.     

A New Algorithm for Wind-Vector Retrieval From Scatterometers

A new efficient algorithm for retrieving wind-vector solutions from scatterometers is developed based on a criterion of minimum normalized standard deviation (NSD) of wind speed derived from backscatter measurements using a geophysical model function (GMF). Its performance has been evaluated through simulations using QSCAT-1 GMF and the QuikSCAT observational geometry. The present algorithm, named the NSD algorithm, is found to be computationally more efficient (two to three times) besides being at par with the maximum-likelihood estimator (MLE) algorithm in terms of retrieval skill, retrieval errors, and distribution of solutions, on the basis of simulations as well as comparison of limited QuikSCAT-data-derived winds with National Centers for Environmental Prediction and European Centre for Medium-Range Weather Forecasts model winds. Simulation results and analysis of sample QuikSCAT data are presented.      

Modeling study of the effects of the coagulation kernel with van der Waals forces and turbulence on the particle size distribution

In this study, variations in the size distributions due to different assumptions for the coagulation kernel are investigated. In order to evaluate how the coagulation kernel influences the form of the particle size distribution and to describe the process of the Brownian coagulation, the Brownian coagulation kernel, including the van der Waals forces, is compared with the most frequently used coagulation coefficients. Retardation should be considered for interparticle interactions for particles larger than 1μm radius. However, for particle sizes larger than 0.1 μm, the Brownian kernel is not dominant, so the retardation effect can be ignored. The inclusion of the van der Waals forces in the Brownian coagulation kernel caused a faster coagulation process in the small particle size range. Taking account of the turbulent coagulation kernel, the turbulent coagulation kernel becomes more important when the turbulent intensity is higher. The turbulent coagulation kernel affects the large particle size range and ignoring the turbulent coagulation kernel will lead to overestimation of particle number concentration in model simulation. The results of this study indicate that the inclusion of van der Waals forces or the turbulent coagulation kernel in the total coagulation kernel impacts on the modeled particle size distributions and total particle number concentration.