[1] Large, sustained well water level changes (>10 cm) in response to distant (more than hundreds of kilometers) earthquakes have proven enigmatic for over 30 years. Here we use high sampling rates at a well near Grants Pass, Oregon, to perform the first simultaneous analysis of both the dynamic response of water level and sustained changes, or steps. We observe a factor of 40 increase in the ratio of water level amplitude to seismic wave ground velocity during a sudden coseismic step. On the basis of this observation we propose a new model for coseismic pore pressure steps in which a temporary barrier deposited by groundwater flow is entrained and removed by the more rapid flow induced by the seismic waves. In hydrothermal areas, this mechanism could lead to 4 Â 10 À2 MPa pressure changes and triggered seismicity.

Earthquakes have been observed to affect hydrological systems in a variety of ways--water well levels can change dramatically, streams can become fuller and spring discharges can increase at the time of earthquakes. Distant earthquakes may even increase the permeability in faults. Most of these hydrological observations can be explained by some form of permeability increase. Here we use the response of water well levels to solid Earth tides to measure permeability over a 20-year period. At the time of each of seven earthquakes in Southern California, we observe transient changes of up to 24 degrees in the phase of the water level response to the dilatational volumetric strain of the semidiurnal tidal components of wells at the Piñon Flat Observatory in Southern California. After the earthquakes, the phase gradually returns to the background value at a rate of less than 0.1 degrees per day. We use a model of axisymmetric flow driven by an imposed head oscillation through a single, laterally extensive, confined, homogeneous and isotropic aquifer to relate the phase response to aquifer properties. We interpret the changes in phase response as due to changes in permeability. At the time of the earthquakes, the permeability at the site increases by a factor as high as three. The permeability increase depends roughly linearly on the amplitude of seismic-wave peak ground velocity in the range of 0.21-2.1 cm s(-1). Such permeability increases are of interest to hydrologists and oil reservoir engineers as they affect fluid flow and might determine long-term evolution of hydrological and oil-bearing systems. They may also be interesting to seismologists, as the resulting pore pressure changes can affect earthquakes by changing normal stresses on faults.

Approximately 0.4% of explosive volcanic eruptions occur within a few days of large, distant earthquakes. This many "triggered" eruptions is much greater than expected by chance. Several mechanisms have been proposed to explain triggering through changes in magma overpressure, including the growth of bubbles, the advection of large pressures by rising bubbles, and overturn of magma chambers. Alternatively, triggered eruptions may occur through failure of rocks surrounding stored magma. All these mechanisms require a process that enhances small static stress changes caused by earthquakes or that can convert (the larger) transient, dynamic strains into permanent changes in pressure. All proposed processes, in addition to viscoelastic relaxation of stresses, can result in delayed triggering of eruptions, though quantifying the connection between earthquakes and delayed, triggered eruptions is much more challenging. Mud volcanoes and geysers also respond to distant earthquakes. Mud volcanoes, that discharge mud from depths greater than many hundreds of meters, may be triggered by liquefaction caused by shaking, and may thus be similar to small mud volcanoes that originate within a few meters of the surface. Changes in permeability of the matrix surrounding main geyser conduits, by opening or creating new fractures, may explain the observed changes in their eruption frequency.

Seismologists have never directly observed rupture in Earth’s interior. Instead, they glean information from seismic waves, geodetic measurements, and numerical experiments.

[1] Oscillations in stress, such as those created by earthquakes, can increase permeability and fluid mobility in geologic media. In natural systems, strain amplitudes as small as 10À6 can increase discharge in streams and springs, change the water level in wells, and enhance production from petroleum reservoirs. Enhanced permeability typically recovers to prestimulated values over a period of months to years. Mechanisms that can change permeability at such small stresses include unblocking pores, either by breaking up permeability-limiting colloidal deposits or by mobilizing droplets and bubbles trapped in pores by capillary forces. The recovery time over which permeability returns to the prestimulated value is governed by the time to reblock pores, or for geochemical processes to seal pores. Monitoring permeability in geothermal systems where there is abundant seismicity, and the response of flow to local and regional earthquakes, would help test some of the proposed mechanisms and identify controls on permeability and its evolution.

Earthquakes occur as a result of global plate motion. However, this simple picture is far from complete. Some plate boundaries glide past each other smoothly, while others are punctuated by catastrophic failures. Some earthquakes stop after only a few hundred metres while others continue rupturing for a thousand kilometres. Earthquakes are sometimes triggered by other large earthquakes thousands of kilometres away. We address these questions by dissecting the observable phenomena and separating out the quantifiable features for comparison across events. We begin with a discussion of stress in the crust followed by an overview of earthquake phenomenology, focusing on the parameters that are readily measured by current seismic techniques. We briefly discuss how these parameters are related to the amplitude and frequencies of the elastic waves measured by seismometers as well as direct geodetic measurements of the Earth's deformation. We then review the major processes thought to be active during the rupture and discuss their relation to the observable parameters. We then take a longer range view by discussing how earthquakes interact as a complex system. Finally, we combine subjects to approach the key issue of earthquake initiation. This concluding discussion will require using the processes introduced in the study of rupture as well as some novel mechanisms. As our observational database improves, our computational ability accelerates and our laboratories become more refined, the next few decades promise to bring more insights on earthquakes and perhaps some answers.

[1] We report on the topographic roughness measurements of five exhumed faults and thirteen surface earthquake ruptures over a large range of scales: from 50 mm to 50 km. We used three scanner devices (LiDAR, laser profilometer, white light interferometer), spanning complementary scale ranges from 50 mm to 10 m, to measure the 3-D topography of the same objects, i.e., five exhumed slip surfaces (Vuache-Sillingy, Bolu, Corona Heights, Dixie Valley, Magnola). A consistent geometrical property, i.e., self-affinity, emerges as the morphology of the slip surfaces shows at first order, a linear behavior on a log-log plot where axes are fault roughness and spatial length scale, covering five decades of length-scales. The observed fault roughness is scale dependent, with an anisotropic self-affine behavior described by four parameters: two power law exponents H, constant among all the faults studied but slightly anisotropic (H k = 0.58 AE 0.07 in the slip direction and H ? = 0.81 AE 0.04 perpendicular to it), and two pre-factors showing variability over the faults studied. For larger scales between 200 m and 50 km, we have analyzed the 2-D roughness of the surface rupture of thirteen major continental earthquakes. These ruptures show geometrical properties consistent with the slip-perpendicular behavior of the smaller-scale measurements. Our analysis suggests that the inherent non-alignment between the exposed traces and the along or normal slip direction results in sampling the slip-perpendicular geometry. Although a data gap exists between the scanned fault scarps and rupture traces, the measurements are consistent within the error bars with a single geometrical description, i.e., consistent dimensionality, over nine decades of length scales.

The majority of earthquakes are aftershocks, yet aftershock physics is not well understood. Many studies suggest that static stress changes trigger aftershocks, but recent work suggests that shaking (dynamic stresses) may also play a role. Here we measure the decay of aftershocks as a function of distance from magnitude 2-6 mainshocks in order to clarify the aftershock triggering process. We find that for short times after the mainshock, when low background seismicity rates allow for good aftershock detection, the decay is well fitted by a single inverse power law over distances of 0.2-50 km. The consistency of the trend indicates that the same triggering mechanism is working over the entire range. As static stress changes at the more distant aftershocks are negligible, this suggests that dynamic stresses may be triggering all of these aftershocks. We infer that the observed aftershock density is consistent with the probability of triggering aftershocks being nearly proportional to seismic wave amplitude. The data are not fitted well by models that combine static stress change with the evolution of frictionally locked faults.