Permeability controls fluid flow in fault zones and is a proxy for rock damage after an earthquake. We used the tidal response of water level in a deep borehole to track permeability for 18 months in the damage zone of the causative fault of the 2008 moment magnitude 7.9 Wenchuan earthquake. The unusually high measured hydraulic diffusivity of 2.4 × 10(-2) square meters per second implies a major role for water circulation in the fault zone. For most of the observation period, the permeability decreased rapidly as the fault healed. The trend was interrupted by abrupt permeability increases attributable to shaking from remote earthquakes. These direct measurements of the fault zone reveal a process of punctuated recovery as healing and damage interact in the aftermath of a major earthquake.

Quantitative interpretation of the tidal response of water levels measured in wells has long been made either with a model for perfectly confined aquifers or with a model for purely unconfined aquifers. However, many aquifers may be neither totally confined nor purely unconfined at the frequencies of tidal loading but behave somewhere between the two end‐members. Here we present a more general model for the tidal response of groundwater in aquifers with both horizontal flow and vertical leakage. The model has three independent parameters: the transmissivity (T) and storativity (S) of the aquifer and the specific leakage (K′/b′) of the leaking aquitard, where K′ and b′ are the hydraulic conductivity and the thickness of the aquitard, respectively. If T and S are known independently, this model may be used to estimate aquitard leakage from the phase shift and amplitude ratio of water level in wells obtained from tidal analysis. We apply the model to interpret the tidal response of water level in a US Geological Survey (USGS) deep monitoring well installed in the Arbuckle aquifer in Oklahoma, into which massive amount of wastewater coproduced from hydrocarbon exploration has been injected. The analysis shows that the Arbuckle aquifer is leaking significantly at this site. We suggest that the present method may be effective and economical for monitoring leakage in groundwater systems, which bears on the safety of water resources, the security of underground waste repositories, and the outflow of wastewater during deep injection and hydrocarbon extraction.

The 11 March 2011 off the Pacific coast of Tohoku (M w 9.0) Earthquake ruptured a 200 km wide megathrust fault, with average displacements of ∼15-20 m. Early estimates of the co-seismic slip distribution using seismic, geodetic and tsunami observations vary significantly in the placement of slip, particularly in the vicinity of the trench. All methods have difficulty resolving the up-dip extent of rupture; onshore geodetic inversions have limited sensitivity to slip far offshore, seismic inversions have instabilities in seismic moment estimation as subfault segments get very shallow, and tsunami inversions average over the total region of ocean bottom uplift. Seismic wave estimates depend strongly on the velocity structure used in the model, which affects both seismic moment estimation and inferred mapping to slip. We explore these ideas using a least-squares inversion of teleseismic P-waves that yields surprisingly large fault displacements (up to ∼60 m) at shallow depth under a protrusion of the upper plate into the trench. This model provides good prediction of GPS static displacements on Honshu. We emphasize the importance of poorly-constrained rigidity variations with depth for estimating fault displacement near the trench. The possibility of large slip at very shallow depth holds implications for up-dip strain accumulation and tsunamigenic earthquake potential of megathrusts elsewhere.

Hydrogeologic properties of fault zones are critical to faulting processes; however, they are not well understood and difficult to measure in situ, particularly in low‐permeability fractured bedrock formations. Analysis of continuous water level response to Earth tides in monitoring wells provides a method to measure the in situ hydrogeologic properties. We utilize four monitoring wells within the San Andreas Fault zone near Logan Quarry to study the fault zone hydrogeologic architecture by measuring the water level tidal response. The specific storage and permeability inferred from the tidal response suggest that there is a difference in properties at different distances from the fault. The sites closer to the fault have higher specific storage and higher permeability than farther from the fault. This difference of properties might be related to the fault zone fracture distribution decreasing away from the fault. Although permeability channels near faults have been documented before, the difference in specific storage near the fault is a new observation. The inferred compliance contrast is consistent with prior estimates of elastic moduli in the near‐fault environment, but the direct measurements are new. The combination of measured permeability and storage yields a diffusivity of about 10−2 m2/s at all the sites both near and far from the fault as a result of the competing effects of permeability and specific storage. This uniform diffusivity structure suggests that the permeability contrast might not efficiently trap fluids during the interseismic period.