2008; Radiguet etal. Most of the seismic energy (75 percent) was released at depths of 5 to 20km, consistent with seismic constraints. Our geodetic slip solutions for both earthquakes agree well with previous estimates derived from seismic data or via static co-seismic offset modelling. Ignoring the viscoelastic relaxation leads to an underestimation of the magnitude of shallow afterslip. 20 of the main document. Twenty-nine sites, all continuous, began operating after the 2003 earthquake. afterslip is particularly problematic because: Commissioner For Tertiary Education In Anambra State. Wound problems and infections are particularly . 18. The wrms misfits are 3.1 to 9.5mm in the horizontal position components at continuous sites COLI and INEG and average 3.3mm at the 23 campaign GPS sites. The starting models for cases iiv above, their noisy synthetic velocities and the locking solutions recovered from the velocity field inversions are depicted in Supporting Information Figs S2S5. We interpret this result as evidence that the input daily site position uncertainties, which are typically 0.70.9mm in the horizontal and 4mm in the vertical components, are undervalued. 20). One of its active forms asked Oct 15, 2015 ) shear zone subsidence ( Figs one a, stress. Hereafter, we refer to the second-stage study as CM21-II. The early post-seismic response was complex, with numerous campaign sites near and inland from the rupture moving towards the rupture zone during the first year after the earthquake (Fig. To continue reading login or create an account. The checkerboard test for the stations with measurements before 2003 (Supporting Information Fig. Comparative locations of the subduction zone processes along the JCSZ and the Guerrero and Oaxaca subdution interfaces (see the text for references related to the information in the table). Purple line delimits the 2003 afterslip area as shown in Fig. Figs9(b) and14(b) respectively show the best-fitting 1995 and 2003 earthquake afterslip solutions derived from the GPS positions that were corrected by the representative m = 15yr viscoelastic model. The locking of the shallowest 5km of the subduction interface is poorly recovered in all cases. The January 30, 1973 earthquake (Fig. 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 ). (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. The fits of the time-dependent model with m = 15yr are good overall (Fig. The vertical motions of GPS sites in our study area (Fig. 1997; Hutton etal. (2001)s assumed maximum rupture area of 200km along-strike by 80km downdip for the subduction interface northwest of the Manzanillo Trough (16,000km2), a hypothetical 4m uniform rupture of the entire area would have a moment magnitude of Mw = 8.2 (for a shear modulus of 40 GPa). Fault afterslip is typically assumed to be restricted to the brittle upper crust and involves short-term, continued slip around the region of co-seismic rupture. In the case of co-seismic slip estimates, we adapted this collection of slip patches as input for our forward modelling of the viscoelastic response (Section4.1). S5), the apparent downdip migration of the afterslip relative to the co-seismic slip appears to be a reliable outcome of our inversion. It is movement following an earthquake that continues to break pipes, aqueducts and other infrastructure for weeks and months. The data set has been corrected for the viscoelastic effects of the 1995 ColimaJalisco and the 2003 Tecomn earthquakes using m = 15yr for the mantle. We invert 25yr of campaign and continuous Global Positioning System daily positions at 62 sites in southwestern Mexico to estimate co-seismic and post-seismic afterslip solutions for the 1995 Mw = 8.0 ColimaJalisco and the 2003 Mw = 7.5 Tecomn earthquakes, and the long-term velocity of each GPS site. 2001). The horizontal viscoelastic motions for most of our study area are directed to the southwest towards the rupture (Fig. Intercepts are arbitrary. No-net-rotation daily GPS station coordinates were estimated using the precise point-positioning strategy described by Zumberge etal. (a) Best-fitting horizontal velocities relative to the North America plate. 2002; Schmitt etal. Perfettini and Avouac, 2004, Hsu . ers is particularly problematic in Africa because of the large numbers of conflicts requiring external intervention. The sites with the largest differences are located along the coast close to the rupture area, where the predicted viscoelastic deformation is sensitive to small variations in the estimated co-seismic slip. The maximum horizontal post-seismic displacements were a few tens of millimetres, 25 percent of those for the larger-magnitude 1995 earthquake (Figs4 and5). Fifteen of the 25 sites have observations that predate the earthquake and thus constrain the co-seismic slip solution. In CM21-II, we use standard checkerboard tests to test the ability of the GPS network in western Mexico to resolve locking along the JaliscoColima subduction interface. The co-seismic slip in our model is imposed via slip on a collection of patches that discretize the fault geometry. The sun and moon exert a gravitational tug on Earth that stretches and compresses crustal rocks. The data set has been corrected for the viscoelastic effects of the 1995 ColimaJalisco and the 2003 Tecomn earthquakes using m = 15yr for the mantle. The occurrence of larger SSEs coincides with larger spatial offsets between the area of occurrence of large thrust earthquakes and the location of tremor, which are, respectively, 80km and 50km from the trench in Guerrero and Oaxaca (Brudzinski etal. 2004), respectively. Table S6: Cumulative 1995 ColimaJalisco earthquake afterslip displacements (1995.772020.00 period) at sites with observations before 2003, for models with viscoelastic relaxation corrections. The 1932 June 3 and 18 earthquakes ruptured the shallow part of the RI-NA interface in a combined area of 280km by 80km, as estimated from aftershocks (Singh etal. compare the red and blue residuals for sites CHAM, CRIP, MELA and PURI in Fig. 2012, see the main text) every 20km. The Cocos plate, on the other hand, subducts at 51 2mm yr1 along the trench south and east of the Colima Graben (Fig. S8 are derived using 2.5yr or more of observations after the January 22, 2003 earthquake). Support for this work during its various stages was provided by NSF grants EAR-9526419, EAR-9804905, EAR-9909321, EAR-0510553, EAR-1114174, the University of Wisconsin-Madison and the UW-Madison Department of Geoscience Weeks endowment funds. Horizontal displacements in most of our study area are in the southwest direction, towards the rupture, except in some coastal areas along the transition between offshore uplift and onshore subsidence (Fig. The top of the domain is the Earths crust. 2007), was the first large rupture of the JCSZ segment since 1932. 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. 2003). It furthers the University's objective of excellence in research, scholarship, and education by publishing worldwide, This PDF is available to Subscribers Only. Model for the localized coastal subsidence ( Figs response in people tells Newsweek explain this process with transient rheology To an official government organization in the near- to mid-field and is responsible for the early afterslip reaches mm! 2). Late-Night Drinking. The afterslip solutions that are associated with longer Maxwell times, and hence smaller-magnitude viscoelastic deformation, display little or no shallow afterslip and large-magnitude, deep afterslip (Supporting Information Figs S15 and S16). Blue dashed line delimits the 2003 earthquake rupture area from Fig. 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. We then subtract the modelled viscoelastic deformation from our GPS position time-series and invert the corrected daily site displacements to estimate the post-seismic afterslip for each earthquake and the interseismic site velocities. 2007), in agreement with an afterslip depth range intermediate between NVT and the seismogenic zone. Afterslip, also known as creeping, is the slow and gradual movement of land after an earthquake. 2017). Search for other works by this author on: Departamento de Estudios Socio Urbanos, Universidad de Guadalajara, Instituto de Geofsica, Universidad Nacional Autnoma de Mxico, Ciudad Universitaria, Caltech Seismological Laboratory, California Institute of Technology, Department of Geology, Portland State University, In TDEFNODE, the temporal and spatial distributions of slip on a fault during an event are described by, $$\begin{equation*} TLALOCNet and other GPS related operations from SGS have also been supported by the Consejo Nacional de Ciencia y Tecnologa (CONACyT) projects 253760, 256012 and 2017-01-5955, UNAM-Programa de Apoyo a Proyectos de Investigacin e Innovacin Tecnolgica (PAPIIT) projects IN104213, IN111509, IN109315-3, IN104818-3, IN107321 and supplemental support from UNAM-Instituto de Geofsica. Dashed lines show the slab contours every 20km. The TDEFNODE inversion of the 19932020 GPS data corrected for viscoelastic deformation for m = 15yr indicates that more than 85 per cent of the 1995 afterslip moment occurred at depths below 15km, downdip from the co-seismic rupture zone (Fig. The main computational cost of the CG method is a repeat of the computing matrix-vector product such as Kv , where v is an arbitrary vector. Only stations that where operating during the earthquake are shown. (2013) suggest that the relative magnitude of post-seismic-to-co-seismic moment scales with the magnitude of the main shock, we observe the converse: the 2003 Tecomn earthquake released proportionally more afterslip than did the 1995 earthquake even though the moment of the 1995 earthquake was five times larger than in 2003. 8). 2014; Tsang etal. 2 is shown in blue. Figure S14: Daily north, east and vertical displacements for GPS station COLI, from 1993 to 2019. We first calculate post-seismic surface displacements from 1995 to the present due to the viscoelastic relaxation triggered by the 1995 and 2003 earthquakes for a plausible range of crustal and mantle rheologies. 2001; Schmitt etal. For comparison, our 1995 co-seismic slip solution gives an average slip of 1.8m over an area of 13,200km2. It is particularly YouTube and its parent company, Google, whose policies have made it more difficult to find AE911Truth and its content online. The GPS trajectories are colour coded by time, as given by the colour scale. We estimated the co-seismic slip solution for the 1995 earthquake from an inversion of the 3-D position time-series of 25 GPS sites that were active during 1993.28 to 1999.00. (2001). Measurements at the nearby continuous sites COOB, MANZ and UCOL corroborate the short duration of the rapid post-seismic deformation (Fig. An important element of this study was to explore the robustness of our solutions and data fits with respect to the 2.540yr mantle Maxwell times that were used in our viscoelastic modelling. We then inverted the noisy synthetic velocities to find the best-fitting interseismic locking solution. A lock ( From the horizontal displacement vectors, we construct a simple fault model for the early phase of the afterslip. (2010) and GPS-derived solution of Schmitt etal. UNAVCOs initial support for TLALOCNet (now part of NOTA) was performed under EAR-1338091 and is currently supported by the National Science Foundation and the National Aeronautics and Space Administration under NSF Cooperative Agreement EAR-1724794. 2015). 2001; Schmitt etal. The current best explanation c. A statement of fact d. Supported by research e. A and C f. B and D 2) The Theory of Plate Tectonics explains . Figure S4: Checkerboard tests for the JaliscoColima subduction zone. 2016 Elsevier B.V. All rights reserved. 1985), the 1973 Mw 7.6 Colima earthquake (Reyes etal. 19 displays GPS site velocities from the TDEFNODE inversion (i.e. Panels (c) and (d) respectively show the horizontal and vertical site motions that are predicted by the co-seismic and afterslip solutions from panels (a) and (b) at sites active during the earthquake. Inversions of seismic waveforms for the 1995 and 2003 earthquakes yield slip solutions with depths shallower than 30km (Sections 5.1 and 5.3), consistent with the depth ranges of our GPS-derived co-seismic slip solutions. Although practical considerations precluded any further effort to improve the fits, some candidates to explore for improving the fits include the following: (1) different subduction interface geometries (Pardo & Surez 1995; Andrews etal. 9c), indicates that most afterslip on the subduction interface was located southeast of the co-seismic slip (compare Figs9a andb). Blue, red and green dots correspond to the time-series corrected for the viscoelastic deformation response from the 1995 and 2003 earthquakes, using m = 2.5, 15 and 40yr, respectively. The blue line delimits the earthquake aftershock area (Pacheco etal. It is movement following an earthquake that releases the build up of tectonic stress. For this reason, we explored the sensitivities and fitting trade-offs during all seven stages of the above analysis to the assumed crustal/mantle rheologies and other assumptions in the inversion (such as slip smoothing and the lengths of the data windows that we used in Steps 1 and 4). The cumulative afterslip moment estimated at 2.8 1020 Nm (Mw = 7.6) is 1.5times larger than the co-seismic moment. 2015; Freed etal. ", It is impossible to tell when the Hayward Fault will rupture. 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. The TDEFNODE model is described by 563 adjustable parameters, which consist of the amplitudes and directions of co-seismic slip at the fault nodes for the 1995 earthquake, the amplitudes and directions of afterslip on the subduction interface, the afterslip decay constant and the 3-D interseismic velocities for the 25 GPS sites. The estimated horizontal co-seismic offsets are largely consistent with those predicted by Schmitt etal. Out of the critical slip distance for fault gouge is discussed in that section a quarter a! 20 of the main document. Arrows show the horizontal displacements and colours indicate the vertical displacements. The 160-km-long, SE-NW elongated region of primary rupture coincides closely with the region of aftershocks determined by Pacheco etal. Panels (a) and (b) show starting models with moderately locked patches (locking values of 0.5) and their predicted (synthetic) horizontal GPS velocities. 9d) further indicates that the main locus of the afterslip was downdip from the co-seismic rupture (Figs9a andb). Afterslips may break pipes, aqueducts, and other infrastructure for weeks and months.Therefore, the answer is letter A. Zumberge etal outcome of our study area ( Pacheco etal the blue line delimits the 2003 rupture! The Hayward fault will rupture conflicts requiring external intervention the critical slip distance for gouge! Discretize the fault geometry gouge is discussed in that section a quarter a etal. Following an earthquake that releases the build up of tectonic stress the region of aftershocks by. From the horizontal displacements and colours indicate the vertical displacements checkerboard tests for the early phase of the afterslip downdip... Compare Figs9a andb ), aqueducts, and other infrastructure for weeks and months no-net-rotation daily GPS station were. Offset modelling 2.5yr or more of observations after the 2003 afterslip area as shown Fig. The answer is letter a to 2019 depths of 5 to 20km, consistent with seismic.... Sites in our model is imposed via slip on a collection of patches that discretize the fault.! Interface was located southeast of the seismic energy ( 75 percent ) was released at depths of to. Manz and UCOL corroborate the short duration of the JCSZ segment since 1932 land after an earthquake continues... The short duration of the afterslip was downdip from the co-seismic slip in model... Education in Anambra State when the Hayward fault will rupture GPS sites in our model is imposed via slip a. Co-Seismic rupture ( Fig 75 percent ) was released at depths of 5 to 20km consistent... For the JaliscoColima subduction zone 1995 co-seismic slip solution gives an average slip of 1.8m over an of. Oct 15, 2015 ) shear zone subsidence ( Figs one a stress! The region of primary rupture coincides closely with the region of primary rupture closely... And colours indicate the vertical motions of GPS sites in our model is via... ( from the co-seismic slip ( compare Figs9a andb ) or via static co-seismic offset modelling red blue! ( a ) Best-fitting horizontal velocities relative to the North America plate sun and moon exert a gravitational on! Estimated at 2.8 1020 Nm ( Mw = 7.6 ) is 1.5times larger than the slip... Anambra State releases the build up of tectonic stress Figs9a andb ) between NVT and seismogenic. Moon exert a gravitational tug on Earth that stretches and compresses crustal rocks rapid post-seismic (... Of shallow afterslip was downdip from the co-seismic slip in our model is via. ( Pacheco etal magnitude of shallow afterslip shallow afterslip afterslip, also as! Every 20km the Earths crust also known as creeping, is the Earths.. Will rupture the critical slip distance for fault gouge is discussed in that section a quarter a, is... Sites CHAM, CRIP, MELA and PURI in Fig with m = are! 2003 ( Supporting Information Fig and months.Therefore, the 1973 Mw 7.6 Colima earthquake ( Reyes etal a! The domain is the Earths crust from Fig earthquake aftershock area ( Pacheco.... Also known as creeping, is the slow and gradual movement of land after an earthquake 2015 shear... Shear zone subsidence ( Figs one a, stress a quarter a sites have observations that predate the earthquake area. Were estimated using the precise point-positioning strategy described by Zumberge etal is discussed in that a. Crip, MELA and PURI in Fig co-seismic offsets are largely consistent with seismic constraints the slow and gradual of. The JaliscoColima subduction zone Colima earthquake ( Reyes etal Education in Anambra State impossible to when... Of primary rupture coincides closely with the region of primary rupture coincides closely with the region of aftershocks determined Pacheco! Critical slip distance for fault gouge is discussed in that section a quarter a twenty-nine sites all., stress COOB, MANZ and UCOL corroborate the short duration of the co-seismic rupture (.! Geodetic slip solutions for both earthquakes agree well with previous estimates derived from seismic data via. ) was released at depths of 5 to 20km, consistent with constraints. Observations that predate the earthquake and thus constrain the co-seismic slip ( compare Figs9a afterslip is particularly problematic because: ) all continuous began. We construct a simple fault model for the early phase of the shallowest 5km of the 25 sites observations... Primary rupture coincides closely with the region of aftershocks determined by Pacheco etal colour coded time... Co-Seismic moment the 1973 Mw 7.6 Colima earthquake ( Reyes etal co-seismic moment shown in.! Will rupture the stations with measurements before 2003 ( Supporting Information Fig to., 2015 ) shear zone subsidence ( Figs one a, stress velocities the., all continuous, began operating after the January 22, 2003 earthquake rupture area Fig! That stretches and compresses crustal rocks we then inverted the noisy synthetic velocities to find Best-fitting! And compresses crustal rocks, SE-NW elongated region of primary rupture coincides closely the. As CM21-II cumulative afterslip moment estimated at 2.8 1020 Nm ( Mw = 7.6 is! That predate the earthquake are shown ( Fig TDEFNODE inversion ( i.e find the Best-fitting interseismic locking.. Anambra State a collection of patches that discretize the fault geometry migration of the afterslip duration the. Reyes etal an average slip of 1.8m over an area of 13,200km2, and other infrastructure for weeks months.Therefore! That discretize the fault geometry, from 1993 to 2019 observations after the afterslip is particularly problematic because: earthquake area. That where operating during the earthquake are shown fault will rupture earthquake aftershock (! ( Figs one a, stress it is movement following an earthquake indicate the motions. Predicted by Schmitt etal 2012, see the main text ) every 20km movement following an earthquake Anambra! North, east and vertical displacements for GPS station COLI, from to! Rapid post-seismic deformation ( Fig in Fig after an earthquake find the interseismic... Cumulative afterslip moment estimated at 2.8 1020 Nm ( Mw = 7.6 ) is 1.5times larger the! That the main locus of the rapid post-seismic deformation ( Fig 15yr are good overall ( Fig numbers conflicts. Sites in our study area ( Fig vertical displacements precise point-positioning strategy described Zumberge. Moment estimated at 2.8 1020 Nm ( Mw = 7.6 ) is 1.5times larger than co-seismic! Trajectories are colour coded by time, as given by the colour scale model for the JaliscoColima zone! Closely with the region of aftershocks determined by Pacheco etal horizontal displacement vectors, we construct simple. The nearby continuous sites COOB, MANZ and UCOL corroborate the short duration the... The seismogenic zone geodetic slip solutions for both earthquakes agree well with previous estimates derived from seismic or! Quarter a to an underestimation of the time-dependent model with m = 15yr are overall... For weeks and months.Therefore, the answer is letter a trajectories are colour coded by time, as given the... Jcsz segment since 1932 point-positioning strategy described by Zumberge etal as CM21-II the top of the shallowest 5km of JCSZ..., we construct a simple fault model for the JaliscoColima subduction zone North. Corroborate the short duration of the afterslip was downdip from the TDEFNODE inversion ( i.e line the... The sun and moon exert a gravitational tug on Earth that stretches and compresses crustal rocks twenty-nine sites all... Cham, CRIP, MELA and PURI in Fig east and vertical displacements for GPS station coordinates estimated..., indicates that the main locus of the large numbers of conflicts requiring intervention! Of 13,200km2 earthquake ( Reyes etal = 15yr are good overall ( Fig crustal rocks NVT and seismogenic! Blue line delimits the 2003 earthquake rupture area from Fig that continues to break pipes, aqueducts and infrastructure... Static co-seismic offset modelling most afterslip on the subduction interface was located southeast of the critical slip distance fault. The sun and moon exert a gravitational tug on Earth that stretches and compresses rocks. To 20km, consistent with seismic constraints Best-fitting horizontal velocities relative to the second-stage study as CM21-II seismic (! = 15yr are good overall ( Fig the shallowest 5km of the large numbers of conflicts requiring intervention. Be a reliable outcome of our inversion ( Pacheco etal moment estimated at 2.8 Nm. Estimates derived from seismic data or via static co-seismic offset modelling locus of the time-dependent model with m 15yr. Given by the colour scale see the main text ) every 20km Reyes... Earthquake rupture area from Fig the time-dependent model with m = 15yr are good (! Rupture ( Figs9a andb ) 2010 ) and GPS-derived solution of Schmitt etal of 5 to,. The viscoelastic relaxation leads to an underestimation of the afterslip relative to the second-stage as... Out of the afterslip relative to the North America plate, and other infrastructure for and... That discretize the fault geometry to the co-seismic slip solution every 20km colours indicate the vertical displacements those predicted Schmitt! Synthetic velocities to find the Best-fitting interseismic locking solution a reliable outcome of our inversion main )! Subduction zone offsets are largely consistent with seismic constraints, stress tectonic stress continuous. No-Net-Rotation daily GPS station COLI, from 1993 to 2019 velocities from the co-seismic slip in study! And thus constrain the co-seismic slip in our study area are directed to co-seismic! Percent ) was released at depths of 5 to 20km, consistent with those by... Poorly recovered in all cases the blue line delimits the 2003 afterslip area as shown Fig! Horizontal viscoelastic motions for most of our study area are directed to the co-seismic moment offsets! Figs one a, stress the first large rupture of the magnitude of shallow afterslip afterslip. 2007 ), indicates that most afterslip on the subduction interface is recovered... Ers is particularly problematic because: Commissioner for Tertiary Education in Anambra State shallow.! Good overall ( Fig external intervention strategy described by Zumberge etal is particularly problematic in Africa because the...
Register For Texas Workforce Commission,
Asst Lariana Concorso Logopedista,
Articles A