Geology Reports⌕ Search

Geology topics

M. Cocco

Publications and source records attributed to M. Cocco.

8 recordsLinked to original sources
USGS2007

A global search inversion for earthquake kinematic rupture history: Application to the 2000 western Tottori, Japan earthquake

[1] We present a two-stage nonlinear technique to invert strong motions records and geodetic data to retrieve the rupture history of an earthquake on a finite fault. To account for the actual rupture complexity, the fault parameters are spatially variable peak slip velocity, slip direction, rupture time and risetime. The unknown parameters are given at the nodes of the subfaults, whereas the parameters within a subfault are allowed to vary through a bilinear interpolation of the nodal values. The forward modeling is performed with a discrete wave number technique, whose Green's functions include the complete response of the vertically varying Earth structure. During the first stage, an algorithm based on the heat-bath simulated annealing generates an ensemble of models that efficiently sample the good data-fitting regions of parameter space. In the second stage (appraisal), the algorithm performs a statistical analysis of the model ensemble and computes a weighted mean model and its standard deviation. This technique, rather than simply looking at the best model, extracts the most stable features of the earthquake rupture that are consistent with the data and gives an estimate of the variability of each model parameter. We present some synthetic tests to show the effectiveness of the method and its robustness to uncertainty of the adopted crustal model. Finally, we apply this inverse technique to the well recorded 2000 western Tottori, Japan, earthquake ( Mw 6.6); we confirm that the rupture process is characterized by large slip (3-4 m) at very shallow depths but, differently from previous studies, we imaged a new slip patch (2-2.5 m) located deeper, between 14 and 18 km depth.

Tottori↗
USGS2005

A frictional population model of seismicity rate change

We study models of seismicity rate changes caused by the application of a static stress perturbation to a population of faults and discuss our results with respect to the model proposed by Dieterich (1994). These models assume distribution of nucleation sites (e.g., faults) obeying rate-state frictional relations that fail at constant rate under tectonic loading alone, and predicts a positive static stress step at time to will cause an immediate increased seismicity rate that decays according to Omori's law. We show one way in which the Dieterich model may be constructed from simple general idead, illustratted using numerically computed synthetic seismicity and mathematical formulation. We show that seismicity rate change predicted by these models (1) depend on the particular relationship between the clock-advanced failure and fault maturity, (2) are largest for the faults closest to failure at to, (3) depend strongly on which state evolution law faults obey, and (4) are insensitive to some types of population hetrogeneity. We also find that if individual faults fail repeatedly and populations are finite, at timescales much longer than typical aftershock durations, quiescence follows at seismicity rate increase regardless of the specific frictional relations. For the examined models the quiescence duration is comparable to the ratio of stress change to stressing rate ????/??,which occurs after a time comparable to the average recurrence interval of the individual faults in the population and repeats in the absence of any new load may pertubations; this simple model may partly explain observations of repeated clustering of earthquakes. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Imaging the complexity of an active normal fault system: The 1997 Colfiorito (central Italy) case study

Six moderate magnitude earthquakes (5 < Mw < 6) ruptured normal fault segments of the southern sector of the North Apennine belt (central Italy) in the 1997 Colfiorito earthquake sequence. We study the progressive activation of adjacent and nearby parallel faults of this complex normal fault system using ???1650 earthquake locations obtained by applying a double-difference location method, using travel time picks and waveform cross-correlation measurements. The lateral extent of the fault segments range from 5 to 10 km and make up a broad, ???45 km long, NW trending fault system. The geometry of each segment is quite simple and consists of planar faults gently dipping toward SW with an average dip of 40??-45??. The fault planes are not listric but maintain a constant dip through the entire seismogenic volume, down to 8 km depth. We observe the activation of faults on the hanging wall and the absence of seismicity in the footwall of the structure. The observed fault segmentation appears to be due to the lateral heterogeneity of the upper crust: preexisting thrusts inherited from Neogene's compressional tectonic intersect the active normal faults and control their maximum length. The stress tensor obtained by inverting the six main shock focal mechanisms of the sequence is in agreement with the tectonic stress active in the inner chain of the Apennine, revealing a clear NE trending extension direction. Aftershock focal mechanisms show a consistent extensional kinematics, 70% of which are mechanically consistent with the main shock stress field.

Journal of Geophysical Research B: Solid Earth↗
USGS2001

Solving the dynamic rupture problem with different numerical approaches and constitutive laws

We study the dynamic initiation, propagation and arrest of a 2-D in-plane shear rupture by solving the elastodynamic equation by using both a boundary integral equation method and a finite difference approach. For both methods we adopt different constitutive laws: a slip-weakening (SW) law, with constant weakening rate, and rate- and state-dependent friction laws (Dieterich-Ruina). Our numerical procedures allow the use of heterogeneous distributions of constitutive parameters along the fault for both formulations. We first compare the two solution methods with an SW law, emphasizing the required stability conditions to achieve a good resolution of the cohesive zone and to avoid artificial complexity in the solutions. Our modelling results show that the two methods provide very similar time histories of dynamic source parameters. We point out that, if a careful control of resolution and stability is performed, the two methods yield identical solutions. We have also compared the rupture evolution resulting from an SW and a rate- and state-dependent friction law. This comparison shows that despite the different constitutive formulations, a similar behaviour is simulated during the rupture propagation and arrest. We also observe a crack tip bifurcation and a jump in rupture velocity (approaching the P-wave speed) with the Dieterich-Ruina (DR) law. The rupture arrest at a barrier (high strength zone) and the barrier-healing mechanism are also reproduced by this law. However, this constitutive formulation allows the simulation of a more general and complex variety of rupture behaviours. By assuming different heterogeneous distributions of the initial constitutive parameters, we are able to model a barrier-healing as well as a self-healing process. This result suggests that if the heterogeneity of the constitutive parameters is taken into account, the different healing mechanisms can be simulated. We also study the nucleation phase duration Tn, defined as the time necessary for the crack to reach the half-length Ic. We compare the Tn values resulting from distinct simulations calculated using different constitutive laws and different sets of constitutive parameters. Our results confirm that the DR law provides a different description of the nucleation process than the SW law adopted in this study. We emphasize that the DR law yields a complete description of the rupture process, which includes the most prominent features of SW.

Geophysical Journal International↗
USGS1999

Stress transfer by the 1988-1989 M=5.3 and 5.4 Lake Elsman foreshocks to the Loma Prieta fault: Unclamping at the site of peak mainshock slip

We study the stress transferred by the June 27, 1988, M = 5.3 and August 8, 1989, M = 5.4 Lake Elsman earthquakes, the largest events to strike within 15 km of the future Loma Prieta rupture zone during 74 years before the 1989 M = 6.9 Loma Prieta earthquake. We find that the first Lake Elsman event brought the rupture plane of the second event 0.3–1.6 bars (0.03–0.16 MPa) closer to Coulomb failure but that the Lake Elsman events did not bring the future Loma Prieta hypocentral zone closer to failure. Instead, the Lake Elsman earthquakes are calculated to have reduced the normal stress on (or “undamped”) the Loma Prieta rupture surface by 0.5–1.0 bar (0.05–0.10 MPa) at the site where the greatest slip subsequently occurred in the Loma Prieta earthquake. This association between the sites of peak unclamping and slip suggests that the Lake Elsman events did indeed influence the Loma Prieta rupture process. Unclamping the fault would have locally lowered the resistance to sliding. Such an effect could have been enhanced if the lowered normal stress permitted fluid infusion into the undamped part of the fault. Although less well recorded, the M L = 5.0 1964 and M L = 5.3 1967 Corralitos events struck within 10 km of the southwest end of the future Loma Prieta rupture. No similar relationship between the normal stress change and subsequent Loma Prieta slip is observed, although the high-slip patch southwest of the Loma Prieta epicenter corresponds roughly to the site of calculated Coulomb stress increase for a low coefficient of friction. The Lake Elsman-Loma Prieta result is similar to that for the 1987 M = 6.2 Elmore Ranch and M = 6.7 Superstition Hills earthquakes, suggesting that foreshocks might influence the distribution of mainshock slip rather than the site of mainshock nucleation.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Two-way coupling between Vesuvius eruptions and southern Apennine earthquakes, Italy, by elastic stress transfer

During the past 1000 years, eruptions of Vesuvius have often been accompanied by large earthquakes in the Apennines 50–60 km to the northeast. Statistical investigations had shown that earthquakes often preceded eruptions, typically by less than a decade, but did not provide a physical explanation for the correlation. Here, we explore elastic stress interaction between earthquakes and eruptions under the hypothesis that small stress changes can promote events when the Apennine normal faults and the Vesuvius magma body are close to failure. We show that earthquakes can promote eruptions by compressing the magma body at depth and opening suitably oriented near-surface conduits. Voiding the magma body in turns brings these same normal faults closer to Coulomb failure, promoting earthquakes. Such a coupling is strongest if the magma reservoir is a dike oriented normal to the regional extension axis, parallel to the Apennines, and the near-surface conduits and fissures are oriented normal to the Apennines. This preferred orientation suggests that the eruptions issuing from such fissures should be most closely linked in time to Apennine earthquakes. Large Apennine earthquakes since 1400 are calculated to have transferred more stress to Vesuvius than all but the largest eruptions have transferred to Apennine faults, which may explain why earthquakes more commonly lead than follow eruptions. A two-way coupling may thus link earthquakes and Vesuvius eruptions along a 100-km-long set of faults. We test the statistical significance of the earthquake-eruption correlation in the two-way coupling zone, and find a correlation significant at the 95% confidence level.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Frictional constraints on crustal faulting

We consider how variations in fault frictional properties affect the phenomenology of earthquake faulting. In particular, we propose that lateral variations in fault friction produce the marked heterogeneity of slip observed in large earthquakes. We model these variations using a rate- and state-dependent friction law, where we differentiate velocity-weakening behavior into two fields: the strong seismic field is very velocity weakening and the weak seismic field is slightly velocity weakening. Similarly, we differentiate velocity-strengthening behavior into two fields: the compliant field is slightly velocity strengthening and the viscous field is very velocity strengthening. The strong seismic field comprises the seismic slip concentrations, or asperities. The two “intermediate” fields, weak seismic and compliant, have frictional velocity dependences that are close to velocity neutral: these fields modulate both the tectonic loading and the dynamic rupture process. During the interseismic period, the weak seismic and compliant regions slip aseismically, while the strong seismic regions remain locked, evolving into stress concentrations that fail only in main shocks. The weak seismic areas exhibit most of the interseismic activity and aftershocks but can also creep seismically. This “mixed” frictional behavior can be obtained from a sufficiently heterogeneous distribution of the critical slip distance. The model also provides a mechanism for rupture arrest: dynamic rupture fronts decelerate as they penetrate into unloaded complaint or weak seismic areas, producing broad areas of accelerated afterslip. Aftershocks occur on both the weak seismic and compliant areas around a fault, but most of the stress is diffused through aseismic slip. Rapid afterslip on these peripheral areas can also produce aftershocks within the main shock rupture area by reloading weak fault areas that slipped in the main shock and then healed. We test this frictional model by comparing the seismicity and the coseismic slip for the 1966 Parkfield, 1979 Coyote Lake, and 1984 Morgan Hill earthquakes. The interevent seismicity and aftershocks appear to occur on fault areas outside the regions of significant slip: these regions are interpreted as either weak seismic or compliant, depending on whether or not they manifest interevent seismicity.

Journal of Geophysical Research B: Solid Earth↗