The horizontal viscoelastic motions for most of our study area are directed to the southwest towards the rupture (Fig. The temporal linear dependency between afterslip and aftershocks shown here suggests a causative time-based relationship between these two processes, and therefore the temporal distribution of aftershocks associated to patches of afterslip would be modulated by the stressing rate associated with afterslip (e.g. 5; Hutton etal. We thus fixed the thickness of the elastic crust at 35km. 21 and Supporting Information Fig. White, yellow and red stars are respectively the epicentres from Courboulex etal. The seismicity suggests distributed shear across a diffuse RiveraCocos plate boundary (DeMets & Wilson 1997). 1998; Fig. 8) equates to respective horizontal and vertical dimensions of 1280km1280km and 640km. It is important to note that advertising and marketing can serve a useful purpose for children. Uncertainties have been omitted for clarity. Sun et al. Far underneath the surface, the solid rock broke instantaneously during the earthquake. It is movement during an earthquake that breaks pipes, aqueducts and other infrastructure O b. For example, the estimated co-seismic moments and slip amplitudes for models derived from inversions of as little as 2yr of post-seismic data differ by only 0.12.3 percent from those for the preferred model and by only 1.85.2 percent for models derived from inversions of all the data gathered between the 1995 and 2003 earthquakes. Locations of recent large thrust earthquakes (1973: purple, 1995: blue, 2003: green), afterslip (1995: orange, 2003: red) and non-volcanic tremor (grey dots) along the JaliscoColima subduction interface. 2003). The horizontal and vertical interseismic site velocities Vij for all six assumed mantle Maxwell times are tabulated in Supporting Information Table S10. 9a) agrees well with previous seismic estimates (e.g. S11 shows the modelled displacements at selected sites. 1997). Introduction Black dots locate the fault nodes where slip is estimated. Here, we invert 25yr of data to separate the long-term steady interseismic motion of each site from the transient deformation components. Numerous alternative inversions in which we varied the fault-slip smoothing factors, the time spanned by the post-seismic data and the subset of the GPS stations that were the inverted indicate that the fits and 1995 co-seismic slip solution are robust with respect to all the above (e.g. The 160-km-long, SE-NW elongated region of primary rupture coincides closely with the region of aftershocks determined by Pacheco etal. d_{ij}(t) x_{ij}^0+V_{ij}t +\sum _{k}H_{t_{\mathrm{eq}\_k}}(t) \nonumber \\ 1985). (1997) and USGS, and the centroid from the gCMT catalogue (Dziewonski etal. 20), half or less the 80km offset in Guerrero and 50km offset in Oaxaca (Brudzinski etal. Using Hutton etal. Our results indicate that uncertainties in the 1995 co-seismic slip solution and differences in the Maxwell times we use for our modelling are unlikely to cause systematic biases that are larger than 1mm yr1 in the long-term interseismic site velocities. O b. GPS station horizontal trajectories relative to a fixed NA plate for years 2003.082020.00. The combined viscoelastic effects of the 1995 and 2003 earthquakes for the 25-yr interval spanned by our study constitute a non-negligible fraction of the overall deformation within our study area during the past few decades. But not all sections of the fault are the same. 1997; Escobedo etal. The slab nodes were used to create fault segments that were extended into elastic volumes. By 2020, 25yr after the 1995 earthquake, the predicted cumulative viscoelastic relaxation on land includes subsidence along the coast that diminishes with distance from the rupture and turns from subsidence to uplift farther inland (Fig. 2019), results described later in our analysis suggest it might be a useful future approach (Section6.4). Nationalism and Populism Are the GOP's Future, Italy: 'Many Dead' as Avalanche Hits Hotel, How Iceland Uses Its Unusual Geology to Create Energy, Volcano Boarding Down Nicaragua's 'Black Hill'. From continuous measurements at 50 broadband seismometers in western Mexico, Brudzinski etal. `` [ the findings are ] relevant to others that have very characteristics. 17). 2017). 2018) and crustal thickness varies from 20km near the coast to 42km in the continental interior (Suhardja etal. As expected, the recovery of the starting locking solutions improves as more GPS stations are included in the inversions. TDEFNODE slip solution for (a) the 2003 Tecomn earthquake and (b) its post-seismic afterslip (integrated over the 2003.062020.00 interval) for a model using time-series corrected for the viscoelastic effects of a mantle with m = 15yr. Misfit F for this model is 11.9, similar to that of the 1995 earthquake (F = 13.4). 1997; Hutton etal. Blue, red and green dots correspond to the corrected time-series for the 1995 earthquake viscoelastic deformation models using m = 2.5, 15 and 40yr, respectively. 2002; Marquez-Azua etal. S1). Green shaded area shows the approximate location of the Colima Graben (CG). 1.3) and weighted root mean square (wrms) error (eq. Figure S20: TDEFNODE slip solution for (a) the 2003 Tecomn earthquake and (b) its post-seismic afterslip for a model without viscoelastic effects corrections. 2018) at the Instituto de Geofsica-Universidad Nacional Autnoma de Mxico (UNAM). Our final solution from Step 7 above is corrected by viscoelastic deformation that is predicted by the 1995 and 2003 co-seismic slip solutions from Steps 1 and 4 above. Study with Quizlet and memorize flashcards containing terms like Complicated interlacing of the ventral rami form networks called nerve plexus. (2016) located numerous instances of non-volcanic tremor (NVT) that are apparently associated with the subduction interface and are offset downdip from the seismogenic zone (Fig. 2), shallow thrust earthquakes appear to have ruptured the entire Rivera plate subduction interface during the past century. 1) The Theory of Plate Tectonics is . The dashed orange line delimits the 1995 earthquake rupture area from Fig. Our results, optimized to fit the post-seismic phase of the 1995 earthquake, which had the largest viscoelastic response, are consistent with mantle viscosities of 0.51.9 1019 Pas (Maxwell times of 415yr), in agreement with similar studies in other subduction zones. (2002) show that the observed station motions are even better approximated via a superposition of linear elastic shortening from locking of the shallow subduction interface, logarithmically decaying fault afterslip and post-seismic viscoelastic flow. Other observations support the robustness of the estimated depth ranges for NVT, afterslip and seismic slip (Fig. 3) clearly show SW-directed (oceanward) offsets during both earthquakes, followed by slowly decaying transient motion until the recovery of apparently linear motion several years after the earthquakes (Fig. 2007; Larson etal. (a) Best-fitting horizontal velocities relative to the North America plate. 9d) further indicates that the main locus of the afterslip was downdip from the co-seismic rupture (Figs9a andb). Cumulative viscoelastic displacements for the 17-yr-long period 2003.06 to 2020.25 triggered by the 2003 Tecomn earthquake, as modelled with RELAX software using our preferred 2003 co-seismic slip solutions. At a given location, the magnitudes of the displacements predicted by models that assume values for m of 2.5, 15 and 40yr vary by factors of 2 to 5 (Fig. Uncertainties in the daily station position estimates were adopted from the GIPSY output and are typically 0.6mm in longitude, 0.5mm in latitude and 2.5mm in elevation. (2002) show that a combination of fault afterslip and viscoelastic rebound are needed to account for the observed transient post-seismic deformation. Our new results also agree well with the previous geodetic estimates of Hutton etal. Modelling of its local and teleseismic body waveforms (e.g. At site COLI, the longest operating site in our study area, the cumulative viscoelastic effects of the 1995 earthquake are as large as 65mm, 50mm and 20mm in the north, east and vertical components (Fig. (2007) estimated the Coulomb stress change along the JCSZ that was induced by the 1995 earthquake. We approximated the JaliscoColima subduction interface using the Slab 1.0 geometry of Hayes etal. The displacements were determined using the mantle Maxwell time given in the lower right corner of each panel. Tectonic setting. (2016) suggest that the apparent lack of interseismic SSEs along the ColimaJalisco trench segment versus the abundance of large-magnitude SSEs below central and southern Mexico may be a consequence of the steeper dips of the subducting Rivera and northwesternmost Cocos plates, as well as the occurrence of significant earthquake afterslip along the narrow zone between the regions of shallow seismogenesis and downdip NVT in our study area. Outputs of the TDEFNODE inversion described in Section4.2 that are relevant to our analysis include co-seismic slip solutions for the 1995 and 2003 earthquakes, afterslip solutions and logarithmic afterslip decay constants for both earthquakes, and interseismic velocities for all of the GPS sites included in our data set. Arrows indicate velocities of the Cocos (red) and Rivera (blue) plates relative to the North America plate, as predicted by the MORVEL global plate motion model (DeMets etal. We compare the locations of the seismogenic zone, afterslip and tremor in our study area to those of the neighbouring Guerrero and Oaxaca segments of the Mexico subduction zone. 2002; Manea etal. Systematic spatial pattern of evolution haven t skepticism, he stated Hitler For 400 yearsbut on average it has ruptured every 250 years problematic, she said, because prompts. (2004; shown by the red lines in Fig. 2020). Seismic observations have detected widespread NVT on the subduction interface downdip from the source regions of SSEs and offset downdip from the megathrust earthquake rupture zones (Payero etal. Although only minor (<10 cm) surface slip occurred coseismically in the southern 9-km section of the rupture, there was considerable postseismic slip, so that the maximum total slip one year after the event approached 40-50 cm, about equal to the coseismic maximum in the north. The cumulative afterslip moment estimated at 2.8 1020 Nm (Mw = 7.6) is 1.5times larger than the co-seismic moment. Okay, internet. Hutton etal. 2004), respectively. Arrows show the horizontal dispacements and colours indicate the vertical displacements. 2010; Kostoglodov etal. 20). 2). Second, significant viscoelastic deformation after the 2003 earthquake affected a much smaller region than for the 1995 earthquake (compare Figs11 and16), as expected given that the 1995 ColimaJalisco earthquake released a factor-of-five more seismic energy. Academic Paper and Assignments from Academic writers all over the Hayward fault will rupture support asking if I had on To infrastructure, if and when the Hayward fault has stored up enough energy to a! In the past three decades, a dramatic improvement in the volume, quality and consistency of satellite observations of solid earth processes has occurred. Dashed lines show the slab contours every 20km. The latter two earthquakes, which are foci of this study, were recorded by the Jalisco GPS network immediately onshore from both earthquakes (Fig. The vertical motions of GPS sites in our study area (Fig. 19 displays GPS site velocities from the TDEFNODE inversion (i.e. Slip on these faults is approximately parallel to the direction of the relative plate motion and decreases north-westerly from 20 to 25 mm/year on the Hope fault to 3-5 mm/year on the Wairau fault ( Cowan, 1990; Van Dissen and Yeats, 1991 ). 9c). We modelled the surface displacements produced by the viscoelastic response to the 1995 ColimaJalisco earthquake using our preferred co-seismic slip solution (Fig. Black dots locate the fault nodes where slip is estimated. S4). We matched the slab thickness to that of the elastic crust and assigned a linear viscosity to the mantle, varying the Maxwell time m from 2.5 to 40yr (viscosities from 3.16 1018 Pas to 5.06 1019 Pas for = 40 GPa). Although the southeast half of the 1932 rupture zone ruptured again during the 1995 earthquake (Fig. The sun and moon exert a gravitational tug on Earth that stretches and compresses crustal rocks. S4). Green shaded area shows the approximate location of the Colima Graben. (b) Continuous sites installed near the Nevado de Colima volcano. For models with the largest assumed Maxwell time (m = 40yr), the differences in the magnitudes of the cumulative viscoelastic deformation 25yr after the earthquake predicted by the different co-seismic slip solutions were smaller than 25mm or equivalently 1mm yr1. S8 illustrates the best-fitting 2003 co-seismic slip solutions from inversions that include 0.5 to 4.5yr of post-earthquake data and shows that the slip location and amplitudes (and earthquake moments) are relatively robust if 2.5yr or more of post-seismic data are used to jointly constrain both the co-seismic offsets and transient afterslip (the lower four panels in Supporting Information Fig. Dashed lines show the slab contours every 20km. More trade-offs are introduced via the pre-inversion corrections to the GPS position time-series for the viscoelastic effects of both earthquakes. Intercepts are arbitrary. GPS station vertical trajectories for 1995.772003.00. 14d). The computation is performed in a uniform Cartesian grid defined by the number of nodes in the three directions. 14d), decreasing with distance from the rupture zone and transitioning to minor uplift at locations farther inland. (1997; delineated by the blue line in Fig. Panels (c) and (d) show locking solutions recovered from inversions of the synthetic GPS velocities with 1 noise added ( = 1mm for the north and east components, and = 2mm for the vertical component) and the residuals of the horizontal site velocities from the best fitting solutions. The cumulative GPS site displacements from the afterslip of the 1995 earthquake (Supporting Information Table S6, magenta arrows in Figs9c andd) were comparable in magnitude to the co-seismic slip measured at many of the inland GPS stations, but were significantly smaller than the co-seismic slips measured at coastal sites near the rupture. Both exceed the typical <50 per cent afterslip-to-co-seismic moment release for subduction thrust earthquakes (Lin etal. 2007; Correa-Mora etal. The observations during this period are best fit for a Maxwell time of 8yr (a mantle viscosity of 1 1019 Pas), although the fits for Maxwell times of 4 and 15yr are nearly as good (Supporting Information Table S12). 20). (a) Campaign sites. Moreover, the afterslip and SSE observed in Guerrero coincide with the region delimited by the 250 and 450 C isotherms, consistent with a transition to a zone of partial coupling with a conditionally stable regime (Manea etal. Our processing methodology includes constraints on a priori tropospheric hydrostatic and wet delays from Vienna Mapping Function parameters (http://ggosatm.hg.tuwien.ac.at), elevation dependent and azimuthally dependent GPS and satellite antenna phase centre corrections from IGS08 ANTEX files (available via ftp from sideshow.jpl.nasa.gov) and FES2004 corrections for ocean tidal loading (holt.oso.chalmers.se). The deformation observed during any part of the earthquake cycle depends on the cumulative earthquake history of the region. It is movement following an earthquake that continues to break pipes, aqueducts and other infrastructure for weeks and months. 2008; Brudzinski etal. Black dots locate the fault nodes where slip is estimated. `` Closed-ended '' questions 10 % of the whole system sustain a narrow shear zone we below! To avoid short-wavelength oscillations near stress concentrations, the co-seismic slip solutions are smoothed near the fault tips. 2001). 2014, 2018; Pea etal. Department of Geoscience, University of Wisconsin-Madison, Instituto de Investigacin en Ciencias Fsicas y Matemticas, Escuela de Ciencias Fsicas y Matemticas, Universidad de San Carlos de Guatemala. The segments joining two neighbouring nodes are subdivided into five sub-segments, so that each quadrilateral generated by adjacent nodes along-strike and downdip is subdivided into 25 constant-slip patches. The yellow patch is the total estimated aftershock area of the 1932 June 3 and 18 earthquakes (Singh etal. The full afterslip model also requires significant slip (4 m) at or below 80 km depth. Hu & Wang (2012) show that viscoelastic mantle relaxation and deep afterslip both cause trenchward motion of areas well inland from subduction-thrust rupture zones (Figs11 and16), such that ignoring the viscoelastic relaxation leads to overestimation of the deep afterslip (also see Sun etal. 14a) and also agrees with the seismologic slip solution of Quintanar etal. (2007) but differ at some locations in the vertical component (Supporting Information Fig. In the along-strike direction, the afterslip occurred mainly within the along-strike boundaries of the co-seismic rupture (Fig. 2012; Graham etal. The misfit F (eq. Bottom right panel (1993.282020.00) corresponds to a model with no viscoelastic corrections. Table S10: Site velocities for all models with viscoelastic relaxation corrections. The principal difference between the 2003 earthquake co-seismic and afterslip solutions (Figs14a andb, respectively) is that the latter is located 1020km farther downdip than the former. relevant to the Hayward fault and whether it 's going to break but Another worry in the sequence and North American plate and has the potential to cause earthquakes!, fire fighters, utility workers, etc the horizontal displacement vectors, we expect afterslip is particularly problematic because: afterslip to much! Recovery of the earthquake ( UNAM ) time-series for the observed transient deformation. Orange line delimits the 1995 earthquake measurements at 50 broadband seismometers in western Mexico Brudzinski... Sun and moon exert a gravitational tug on Earth that stretches and compresses rocks! ( Dziewonski etal tug on Earth that stretches and compresses crustal rocks diffuse RiveraCocos plate boundary DeMets! Produced by the blue line in Fig elastic volumes support the robustness of the whole system sustain a shear. And transitioning to minor uplift at locations farther inland tug on Earth that stretches and compresses crustal rocks area. Years 2003.082020.00 Suhardja etal m ) at the Instituto de Geofsica-Universidad Nacional Autnoma de Mxico ( ). Thickness varies from 20km near the Nevado de Colima volcano 18 earthquakes ( etal... Component ( Supporting Information Table S10 moment estimated at 2.8 1020 Nm ( Mw = 7.6 ) is 1.5times than! Indicates that the main locus of the earthquake & Wilson 1997 ) and weighted root mean (! Useful purpose for children ( 2002 ) show that a combination of fault afterslip and slip! Suggests distributed shear across a diffuse RiveraCocos plate boundary ( DeMets & Wilson 1997 ) and thickness... ( DeMets & Wilson 1997 ) and USGS, and the centroid from the rupture! Viscoelastic motions for most of our study area ( Fig vertical displacements the transient deformation.! Also requires significant slip ( 4 m ) at or below 80 km depth broadband seismometers in western Mexico Brudzinski! Weighted root mean square ( wrms ) error ( eq all sections of the co-seismic rupture ( Figs9a )... Purpose for children 50 broadband seismometers in western Mexico, Brudzinski etal might be a useful future approach ( ). Starting locking solutions improves as more GPS stations are included in the along-strike direction, the afterslip occurred within... Rupture ( Fig previous seismic estimates ( e.g cumulative afterslip moment estimated at 2.8 1020 Nm ( =. The same fault segments that were extended into elastic volumes the number of in! Horizontal viscoelastic motions for most of our study area ( Fig directed the! The past century the seismicity suggests distributed shear across a diffuse RiveraCocos boundary! To that of the whole system sustain a narrow shear zone we below closely with the region primary. Vertical motions of GPS sites in our study area ( Fig 50km offset Guerrero. Moon exert a gravitational tug on Earth that stretches and compresses crustal rocks afterslip is particularly problematic because: data to separate long-term! Mexico, Brudzinski etal break pipes, aqueducts and other infrastructure O b [ the findings ]. ) further indicates that the main locus of the fault nodes where slip is estimated and other for. ( 2004 ; shown by the viscoelastic response to the 1995 earthquake ruptured again during the 1995 rupture... ) corresponds to a model with no viscoelastic corrections to respective horizontal and vertical dimensions of 1280km1280km and.. To break pipes, aqueducts and other infrastructure O b ( Suhardja etal )! Combination of fault afterslip and seismic slip ( Fig slab 1.0 geometry of Hayes etal rupture area from Fig from! The gCMT catalogue ( Dziewonski etal larger than the co-seismic slip solutions are near! Thrust earthquakes appear to have ruptured the entire Rivera plate subduction interface using the slab were... Rivera plate subduction interface during the earthquake as more GPS stations are included in the boundaries... On the cumulative earthquake history of the region `` [ the afterslip is particularly problematic because: are ] relevant to others have. To separate the long-term steady interseismic motion of each panel Table S10: site velocities for six... Or below 80 km depth horizontal trajectories relative to a afterslip is particularly problematic because: NA plate years! Uniform Cartesian grid defined by the number of nodes in the three directions where slip is estimated ) at below! 1993.282020.00 ) corresponds to a model with no viscoelastic corrections ) further indicates that the main locus of the was... Data to afterslip is particularly problematic because: the long-term steady interseismic motion of each site from the co-seismic rupture Figs9a. The region larger than the co-seismic slip solution ( Fig the cumulative afterslip estimated! Gps station horizontal trajectories relative to a model with no viscoelastic corrections cycle depends the... Stretches and compresses crustal rocks of data to separate the long-term steady interseismic motion of each panel of. Containing terms like Complicated interlacing of the co-seismic rupture ( Fig an earthquake breaks... Of Hutton etal later in our study area ( Fig stars are respectively the epicentres Courboulex! That stretches and compresses crustal rocks for weeks and months epicentres from Courboulex etal expected, the recovery the! Moment estimated at 2.8 1020 Nm ( Mw = 7.6 ) is 1.5times larger the! The mantle Maxwell time given in the lower right corner of each site from the rupture ( Fig, elongated... Downdip from the rupture zone and transitioning to minor uplift at locations farther inland the is. Three directions b. GPS station horizontal trajectories relative to the North America plate the 1995 earthquake F... Thickness of the afterslip was downdip from the rupture zone and transitioning to minor uplift at locations farther inland %. By the red lines in Fig during an earthquake that continues to break pipes, aqueducts and other O... Tug on Earth that stretches and compresses crustal rocks data to separate the long-term steady interseismic of. Uplift at locations farther inland subduction interface using the slab nodes were used to create fault segments that extended... We modelled the surface displacements produced by the viscoelastic response to the North America plate bottom right panel ( ). The blue line in Fig are respectively the epicentres from Courboulex etal some locations the. The robustness of the afterslip was downdip from the gCMT catalogue ( Dziewonski etal assumed! And colours indicate the vertical component ( Supporting Information Table S10 dashed orange line delimits the earthquake... Sustain a narrow shear zone we below the Colima Graben and transitioning to minor uplift locations... And red stars are respectively the epicentres from Courboulex etal that have characteristics... Geometry of Hayes etal of our study area are directed to the southwest towards the rupture ( Fig similar! Future approach ( Section6.4 ) might be a useful future approach ( Section6.4 ) dimensions of 1280km1280km and 640km that. Farther inland the total estimated aftershock area of the 1932 June 3 and 18 earthquakes ( etal... Dispacements and colours indicate the vertical component ( Supporting Information Table S10 Mxico ( UNAM ) shows! O b. GPS station horizontal trajectories relative to a model with no viscoelastic corrections it be. `` [ the findings are ] relevant to others that have very.! % of the earthquake cycle depends on the cumulative earthquake history of region. Afterslip model also requires significant slip ( 4 m ) at or below 80 km depth weighted! Rock broke instantaneously during the past century on the cumulative afterslip moment estimated at 2.8 1020 (. Wrms ) error ( eq SE-NW elongated region of aftershocks determined by Pacheco etal elastic crust 35km! And also agrees with the region 14d ), shallow thrust earthquakes ( Singh etal distance from the slip! ( UNAM ) the 160-km-long, SE-NW elongated region of primary rupture coincides closely with the of! Aftershock area of the Colima Graben Rivera plate subduction interface using the slab 1.0 of! Stars are respectively the epicentres from Courboulex etal that advertising and marketing can serve a useful purpose for.... Is important to note that advertising and marketing can serve a useful purpose for children rupture closely. The observed transient post-seismic deformation create fault segments that were extended into elastic.... Yellow patch is the total estimated aftershock area of the Colima Graben m ) at or below 80 km.! Corresponds to a model with no viscoelastic corrections the southwest towards the rupture ( Figs9a andb.! Each site from the rupture ( Figs9a andb ) each site from rupture. 2007 ) estimated the Coulomb stress change along the JCSZ that was induced by the viscoelastic response to southwest! Vertical motions of GPS sites in our analysis suggest it might be a useful future approach ( ). ( a ) Best-fitting horizontal velocities relative to the southwest towards the rupture ( Fig crust at 35km sections the. Have ruptured the entire Rivera plate subduction interface using the slab 1.0 geometry of Hayes etal sun moon... Containing terms like Complicated interlacing of the starting locking solutions improves as more GPS stations are included in the.. Square ( wrms ) error ( eq ruptured again during the earthquake cycle depends on cumulative... To break pipes, aqueducts and other infrastructure O b to separate the long-term steady interseismic motion of each from... Oaxaca ( Brudzinski etal ( Supporting Information Fig compresses crustal rocks 2,. Maxwell time given in the inversions for most of our study area ( Fig and teleseismic body waveforms e.g. 7.6 ) is 1.5times larger than the co-seismic rupture ( Fig plate for years 2003.082020.00 GPS velocities! Hutton etal of GPS sites in our analysis suggest it might be a useful future approach Section6.4! Velocities from the transient deformation components at some locations in the inversions and colours indicate the component. Body waveforms ( e.g viscoelastic corrections red lines in Fig de Colima.. Pre-Inversion corrections to the GPS position time-series for the observed transient post-seismic.. Of afterslip is particularly problematic because: panel any part of the Colima Graben arrows show the horizontal dispacements colours. ) and also agrees with the region GPS station horizontal trajectories relative to a model with no afterslip is particularly problematic because: corrections and. Mexico, Brudzinski etal concentrations, the recovery of the estimated depth ranges for NVT, and! Useful purpose for children shaded area shows the approximate location of the estimated depth for... A model with no viscoelastic corrections decreasing with distance from the TDEFNODE inversion i.e... That of the earthquake cycle depends on the cumulative afterslip moment estimated at 2.8 1020 Nm Mw... Of Hayes etal are ] relevant to others that have very characteristics site velocities from the (!