Tectonic significance of an earthquake sequence in the Zacoalco half-graben, Jalisco, Mexico

We studied a sequence of small earthquakes that occurred during the months of April and May of 1997, in Jalisco, southwestern Mexico. The earthquakes were located along a set of active faults that form the Zacoalco half-graben (La Lima fault system), west of Lake Chapala, within the rift–rift–rift triple junction. A total of 33 events were located, with magnitudes ranging from 1.5 to 3.5, recorded by a portable array of broadband seismographs. We identified two groups of events: one corresponding to a shallow normal fault, synthetic to La Lima fault system, and another group associated with a deeper fault. The events that occurred on the synthetic fault show normal faulting oriented on a NW–SE plane, dipping shallowly towards the SW. The other group of mechanisms showed either a normal fault oriented NW–SE and dipping steeply to the NE, or a very shallow-dipping normal fault, dipping to the SW. Earthquake distribution and fault plane solutions suggest that the Zacoalco half-graben developed from blocks that rotate as slip occurs on listric faults. These mechanisms could represent the type of motion expected for larger earthquakes in the area, like the one that occurred in 1568.     

Contributing Areas to Domestic Wells in Dipping Sedimentary Rocks under Extreme Recharge Events

We use particle tracking to determine contributing areas (CAs) to wells for transient flow models that simulate cyclic domestic pumping and extreme recharge events in a small synthetic watershed underlain by dipping sedimentary rocks. The CAs consist of strike-oriented bands at locations where the water table intersects high-hydraulic conductivity beds, and from which groundwater flows to the pumping well. Factors that affect the size and location of the CAs include topographic flow directions, rock dip direction, cross-bed fracture density, and position of the well relative to streams. For an effective fracture porosity (n_e) of 10⁻⁴, the fastest advective travel times from CAs to wells are only a few hours. These results indicate that wells in this type of geologic setting can be highly vulnerable to contaminants or pathogens flushed into the subsurface during extreme recharge events. Increasing n_e to 10⁻³ results in modestly smaller CAs and delayed well vulnerability due to slower travel times. CAs determined for steady-state models of the same setting, but with long-term average recharge and pumping rates, are smaller than CAs in the models with extreme recharge. Also, the earliest-arriving particles arrive at the wells later in the steady-state models than in the extreme-recharge models. The results highlight the importance of characterizing geologic structure, simulating plausible effective porosities, and simulating pumping and recharge transience when determining CAs in fractured rock aquifers to assess well vulnerability under extreme precipitation events.     

The effect of lithology on valley width,terrace distribution,and bedload provenance in a tectonically stable catchment with flat‐lying stratigraphy

How rock resistance or erodibility affects fluvial landforms and processes is an outstanding question in geomorphology that has recently garnered attention owing to the recognition that the erosion rates of bedrock channels largely set the pace of landscape evolution. In this work, we evaluate valley width, terrace distribution, and bedload provenance in terms of reach scale variation in lithology in the study reach and discuss the implications for landscape evolution in a catchment with relatively flat‐lying stratigraphy and very little uplift. A reach of the Buffalo National River in Arkansas was partitioned into lithologic reaches and the mechanical and chemical resistance of the main lithologies making up the catchment was measured. Valley width and the spatial distribution of terraces were compared among the different lithologic reaches. The surface grain size and provenance of coarse (2–90 mm) sediment of both modern gravel bars and older terrace deposits that make up the former bedload were measured and defined. The results demonstrate a strong impact of lithology upon valley width, terrace distribution, and bedload provenance and therefore, upon landscape evolution processes. Channel down‐cutting through different lithologies creates variable patterns of resistance across catchments and continents. Particularly in post‐tectonic and non‐tectonic landscapes, the variation in resistance that arises from the exhumation of different rocks in channel longitudinal profiles can impact local base levels, initiating responses that can be propagated through channel networks. The rate at which that response is transmitted through channels is potentially amplified and/or mitigated by differences between the resistance of channel beds and bedload sediment loads. In the study reach, variation in lithologic resistance influences the prevalence of lateral and vertical processes, thus producing a spatial pattern of terraces that reflects rock type rather than climate, regional base level change, or hydrologic variability. Copyright © 2017 John Wiley & Sons, Ltd.     

The Magma Feeding System of the Klyuchevskaya Group of Volcanoes (Kamchatka)

Doklady Earth Sciences - Using observations from a large-scale seismic station network, we built a seismic tomography model of velocity anomalies for longitudinal (P-wave) and transverse (S-wave)...     

Estimates of Lithospheric Failure Cycle Parameters from Regional Earthquake Catalogues

Izvestiya, Physics of the Solid Earth - Abstract—Based on the generalized frequency–magnitude relation for earthquakes, the failure cycle parameters are estimated from regional seismic...     

Back front of seismicity induced by non‐linear pore pressure diffusion

Borehole fluid injections are accompanied by microseismic activity not only during but also after termination of the fluid injection. Previously, this phenomenon has been analysed, assuming that the main triggering mechanism is governed by a linear pressure diffusion in a hydraulically isotropic medium. In this context the so‐called back front of seismicity has been introduced, which allows to characterize the hydraulic transport from the spatiotemporal distribution of post‐injection induced events. However, rocks are generally anisotropic, and in addition, fluid injections can strongly enhance permeability. In this case, permeability becomes a function of pressure. For such situations, we carry out a comprehensive study about the behaviour and parametrization of the back front. Based on a model of a factorized anisotropic pressure dependence of permeability, we present an approach to reconstruct the principal components of the diffusivity tensor. We apply this approach to real microseismic data and show that the back front characterizes the least hydraulic transport. To investigate the back front of non‐linear pore‐fluid pressure diffusion, we numerically consider a power‐law and an exponential‐dependent diffusivity. To account for a post‐injection enhanced hydraulic state of the rock, we introduce a model of a frozen (i.e., nearly unchanged after the stimulation) medium diffusivity and generate synthetic seismicity. We find that, for a weak non‐linearity and 3D exponential diffusion, the linear diffusion back front is still applicable. This finding is in agreement with microseismic data from Ogachi and Fenton Hill. However, for a strong non‐linear fluid–rock interaction such as hydraulic fracturing, the back front can significantly deviate from a time dependence of a linear diffusion back front. This is demonstrated for a data set from the Horn River Basin. Hence, the behaviour of the back front is a strong indicator of a non‐linear fluid–rock interaction.     

Algal δ15N values detect a wastewater effluent plume in nearshore and offshore surface waters and three-dimensionally model the plume across a coral reef on Maui, Hawai'i, USA

The coral reef at Kahekili, Maui is located ~300 m south of the Lahaina Wastewater Reclamation Facility which uses four Class V injection wells to dispose of 3-5 million gallons of wastewater effluent daily. Prior research documented that the wastewater effluent percolates into the nearshore region of Kahekili. To determine if the wastewater effluent was detectable in the surface waters offshore, we used algal bioassays from the nearshore region to 100 m offshore and throughout the water column from the surface to the benthos. These algal bioassays documented that significantly more wastewater effluent was detected in the surface rather than the benthic waters and allowed us to generate a three-dimensional model of the wastewater plume in the Kahekili coastal region. Samples located over freshwater seeps had the highest δ(15)N values (~30-35‰) and the effluent was detected in surface samples 500 m south and 100 m offshore of the freshwater seeps (~8-11‰).