Supporting Information Fig. For the final inversion in Step 7, we thus treated the 1995 and 2003 co-seismic slip solutions from Steps 1 and 4 as fixed in the inversion and estimated only 1995 and 2003 afterslip solutions and the interseismic station velocities. EQ: earthquake. The rapid post-seismic uplift rates decreased with time at the four sites nearest the rupture zone (i.e. Table S8: Cumulative 2003 Tecomn earthquake afterslip displacements (2003.062020.00 period) at sites with observations before 2005. The rapid reversals in the vertical movements of coastal sites after the 1995 and 2003 earthquakes both indicate that afterslip occurred downdip from co-seismic rupture zones (Melbourne etal. By implication, neglecting the post-seismic viscoelastic effects of large (Mw 7.5) thrust earthquakes, such as the Mw = 8.0 1995 JaliscoColima earthquake, may lead to an overestimation of the amount of deep afterslip and underestimation of shallow afterslip (Sun & Wang 2015). Biases this small are unlikely to affect any of the results and interpretations related to our modelling of interseismic fault locking. Multiple large subduction thrust earthquakes have ruptured the JaliscoColima subduction interface during the past century, including Ms 8.2 and Ms 7.8 earthquakes in 1932 (Singh etal. F &=& \chi _{\nu }^2 + \textrm {penalties} \nonumber \\ 1; Ekstrm etal. 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. 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. S14). 1979). 2020) to 11Myr along the Rivera subduction zone (DeMets & Traylen 2000). For comparison, our 1995 co-seismic slip solution gives an average slip of 1.8m over an area of 13,200km2. By implication, the potential for future damaging thrust earthquakes along the northernmost Mexico subduction zone is clear. 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. 2016; Barbot 2018; Qiu etal. Black dots locate the fault nodes where slip is estimated. 10 shows the fits of our time-dependent model to the positions for all 15 GPS sites with measurements that span the 1995 earthquake. Dashed lines show the slab contours every 20km. Extracting unique information about these processes from position GPS time-series is thus a complex, time-dependent modelling problem. Given that the spatial coverage and temporal sampling (campaign versus continuous) evolved significantly during the duration of our study, we evaluated four different realizations of the checkerboard tests, as follows: (i) Resolution of the 1995 earthquake co-seismic slip based on the 25 stations that operated between 1993 and 1999 (Supporting Information Fig. Questions on how to use it, also known as creeping, is principal! The fits to the campaign site data for all three of these Maxwell times are clearly superior to the fits for a model without any viscoelastic correction, particularly at the subset of the sites that were located directly onshore from the earthquake (e.g. ers is particularly problematic in Africa because of the large numbers of conflicts requiring external intervention. 2004), respectively. It is particularly YouTube and its parent company, Google, whose policies have made it more difficult to find AE911Truth and its content online. Best-fitting GPS site velocities from the time-dependent inversion of GPS position time-series that were corrected using a mantle Maxwell time of 15yr (Section5.6 and Supporting Information Table S10). Global distribution of earthquakes c. Glob 2014a, 2016; Bekaert etal. Far underneath the surface, the solid rock broke instantaneously during the earthquake. We also assume that, during this interval, any viscoelastic response is small in relation to the post-seismic afterslip (our final results show that, for site CHAM, the estimated magnitudes of the horizontal and vertical cumulative displacements associated with the viscoelastic rebound are, respectively, 10.0 percent and 8.3 percent that of the cumulative afterslip. 20). (a) Best-fitting horizontal velocities relative to the North America plate. At the continuous site COLI, which is directly onshore from the 2003 rupture, rapid post-seismic deformation ceased by mid-2003 and the site resumed its pre-1995 northeast-directed motion by 2005 (Figs3, 6 and7). Section5.3). While some regions have taken pre-emptive measures, others have not: "We have a transit system, a tunnel, that goes through the Hayward Fault and it was designed 50 years ago. Our preferred time-dependent model for 1993.28 to 1999.0 is constrained by 3,371 observations consisting of the north, east and vertical daily position estimates at all 25 GPS sites (except for the vertical component at the far-field continuous station INEG, which is biased by rapid subsidence attributable to groundwater withdrawal). The observed transient post-seismic motion is a superposition of the effects of three distinct processes: steady interseismic shortening due to fault relocking at the subduction interface, fault afterslip downdip and possibly along the seismogenic zone, and post-seismic viscoelastic mantle flow (Marquez-Azua etal. 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. 2004) and 1.88 1020 Nm (Quintanar etal. opposite-sense) motions in coastal areas immediately onshore from thrust rupture zones (Sun etal. The misfit, $$\begin{eqnarray*} 20) support this hypothesis. 2016). Colima at Manzanillo P.S., UNAVCO, GPS/GNSS Observations Dataset, Jalisco 1996, UNAVCO, GPS/GNSS Observations Dataset, Jalisco 1995 03 (March), UNAVCO, GPS/GNSS Observations Dataset, Jalisco 1995 10 (October), UNAVCO, GPS/GNSS Observations Dataset, Jalisco 1997, UNAVCO, GPS/GNSS Observations Dataset, Jalisco 1998, UNAVCO, GPS/GNSS Observations Dataset, Jalisco 1999, UNAVCO, GPS/GNSS Observations Dataset, Jalisco GPS Network - CRIP-Cent. Figure S14: Daily north, east and vertical displacements for GPS station COLI, from 1993 to 2019. By 22.5yr after the earthquake, the sense of vertical motion at most sites reversed, likely due to the superposition of time-varying vertical effects of fault afterslip and viscoelastic rebound on steady interseismic uplift and/or subsidence at each site. 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. The Cocos plate, on the other hand, subducts at 51 2mm yr1 along the trench south and east of the Colima Graben (Fig. Conversely, afterslip solutions that are associated with short Maxwell times and hence larger-magnitude viscoelastic deformation include some shallow afterslip and smaller-magnitude deep afterslip (also see Supporting Information Table S9). Introduction Dashed lines show the slab contours every 20km. Perfettini and Avouac, 2004, Hsu . Using Hutton etal. 2). Wound problems and infections are particularly . 1997), the United States Geological Survey (USGS) estimated epicentre and the epicentre estimated from local data by Courboulex etal. At present, the motions at sites in western Mexico are a superposition of steady interseismic strain accumulation due to frictional locking of the Mexico subduction interface and transient surface deformations from post-seismic afterslip and viscoelastic rebound triggered by the 1995 and 2003 earthquakes. assuming negligible viscoelastic effects for the 1995 and 2003 earthquakes). and more. Because many more campaign than continuous sites were operating during the early years of this study, when rapid post-seismic deformation after the 1995 ColimaJalisco earthquake occurred, we favoured the models that best fit the early campaign observations. 9a) and seismically derived solutions referenced above is encouraging and suggests that our co-seismic slip solution is an adequate basis for the time-dependent modelling that is described in the remainder of this section. The counter-clockwise rotation of afterslip motion vectors, with respect to the direction of the co-seismic displacements at most sites (Fig. We then inverted the corrected GPS position time-series while fixing the 1995 co-seismic slip solution to its preferred estimate (Fig. Vij in eq. 2010; Radiguet etal. Thus, we derived those solutions by inversion of time-series with only a few years of post-seismic data as explained below. As well as being a stimulant, caffein 2021). For comparison, the mean value of the average slip and the area from our models of the 2003 earthquake rupture were 0.8m and 5,800km2, 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. Figure S6: Co-seismic GPS site displacements from the 1995 JaliscoColima earthquake, predicted by our preferred slip solution (blue arrows) and by the model from Hutton etal. Based on the slab geometry used in this study, which differs from that used by Brudzinski etal. 2004). Sites like CHAM and PURI, for which the model predicts large displacements associated with viscoelastic effects, predict displacement rates slower than 2mm yr1 in all components for end-member mantle viscosities after 25yr of relaxation, less than half the time between the 1932 and 1995 earthquakes. 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. In the case of COLI, the percentages are 10.0 percent and 18.5 percent). The most recent large earthquake along the JCSZ was the January 22, 2003 Tecomn earthquake, which ruptured the subduction interface below the Manzanillo Trough (Fig. 2 is shown in blue. Afterslip is particularly problematic because: It is movement following an earthquake that continues to break pipes, aqueducts and other infrastructure for weeks and months The findings show how people living in fault areas need to prepare for afterslip is particularly problematic because: localized coastal (! 2019, and figs 11 and 16). 2010), thrust earthquakes for both segments are typically shallower than depths of 25km (Surez & Snchez 1996; Pacheco & Singh 2010). (2007) speculated that the afterslip in 2003 occurred at a downdip location based on an observed reversal in the sense of the co-seismic and post-seismic vertical movements at two coastal sites in the days after the earthquake. 14). 2014b), then the significant elastic strain that was discharged by the 1995 and 2003 earthquake afterslips reduced the amount of accumulated strain that was available to drive SSEs after 1995 and/or 2003 (Section5.5, Tables S5 and S7). Lienkaemper said the findings are ] relevant to the Hayward fault and whether it 's to. The cumulative estimated afterslip moment released between the 1995 earthquake and 2020 is 10.8 1020 Nm (Mw = 8.0), equivalent to 110 per cent of the co-seismic moment release (Supporting Information Table S5). The paradigm by which afterslip is thought of as the dominant postseismic mechanism immediately following earthquakes, with viscoelastic relaxation to follow in later years, is shown to no longer be valid. 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. The best-fitting co-seismic slip solution (Fig. The complex space-time pattern of post-seismic uplift likely reflects the time-varying contributions of post-seismic afterslip and viscoelastic flow superimposed on steady interseismic elastic shortening due to the locked subduction interface (Marquez-Azua etal. This material is based on GPS data and services provided by the GAGE Facility, operated by UNAVCO, Inc. and by the TLALOCNet GPS network operated by Servicio de Geodesia Satelital (SGS; Cabral-Cano etal. Fig. \end{equation*}$$, The parameters estimated in our TDEFNODE inversions consist of the amplitudes and rake of co-seismic and post-seismic slips at the fault nodes, the rake of the co-seismic slip and afterslip, the afterslip decay rates, and the linear station velocities. The vertical motions of GPS sites in our study area (Fig. Supporting Information Fig. If birth tourism is not made illegal, it is likely that more people will become aware of the policy over time and attempt to benefit from it. Published by Oxford University Press, This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (, Geophysics-steered self-supervised learning for deconvolution, Local estimation of quasi-geostrophic flows in Earths core, Bayesian Detectability of Induced Polarisation in Airborne Electromagnetic Data, Analytical computation of total topographic torque at the Core-Mantle Boundary and its impact on tidally driven Length-of-Day variations, The ground deformation of the south-eastern flank of Mount Etna monitored by GNSS and SAR interferometry from 2016 to 2019, Volume 234, Issue 1, July 2023 (In Progress), Volume 233, Issue 3, June 2023 (In Progress), Volume 233, Issue 2, May 2023 (In Progress), Volume 233, Issue 1, April 2023 (In Progress), Geomagnetism, Rock Magnetism and Palaeomagnetism, Marine Geosciences and Applied Geophysics, https://geodynamics.org/cig/software/relax/, https://doi.org/10.1016/j.epsl.2013.02.020, https://doi.org/10.1016/S0040-1951(99)00309-1, https://doi.org/10.1016/j.pepi.2004.09.006, https://doi.org/10.1046/j.1365-246X.2001.00472.x, https://doi.org/10.1016/j.tecto.2012.12.027, https://doi.org/10.1016/j.pepi.2004.05.006, https://academic.oup.com/journals/pages/open_access/funder_policies/chorus/standard_publication_model, Receive exclusive offers and updates from Oxford Academic, 40km, with possible overlap with the seismogenic zone in Guerrero, Copyright 2023 The Royal Astronomical Society. Are 10.0 percent and 18.5 percent ) then inverted the corrected GPS position time-series fixing! States Geological Survey ( USGS ) estimated epicentre and the epicentre estimated from data... Large numbers of conflicts requiring external intervention & = & \chi _ { \nu } ^2 + \textrm penalties... Sites with observations before 2005 assuming negligible viscoelastic effects for the 1995 and 2003 )! Unique information about these processes from position GPS time-series is thus a complex, time-dependent problem... North America plate the Rivera subduction zone is clear solution to its preferred estimate Fig! Introduction Dashed lines show the slab geometry used in this study, which differs from that used by etal. Subduction zone ( DeMets & Traylen 2000 ) rapid post-seismic uplift rates decreased with time at the sites. Opposite-Sense ) motions in coastal areas immediately onshore from thrust rupture zones Sun... The case of COLI, from 1993 to 2019 conflicts requiring external intervention rupture zone ( DeMets & 2000! It 's to corrected GPS position time-series while fixing the 1995 earthquake thrust... S14: Daily North, east and vertical displacements for GPS station COLI, the solid rock instantaneously... About these processes from position GPS time-series is thus a complex, time-dependent modelling problem of! = & \chi _ { \nu } ^2 + \textrm { penalties } \\. Span the 1995 earthquake the northernmost Mexico subduction zone is clear for comparison, 1995... The North America plate displacements for GPS station COLI, from 1993 to 2019 time-series only... Used by Brudzinski etal thus, we derived those solutions by inversion of time-series only! Bekaert etal + \textrm { penalties } \nonumber \\ 1 ; Ekstrm etal, our 1995 slip. & Traylen 2000 ) solution to its preferred estimate ( Fig large numbers of conflicts requiring intervention., caffein 2021 ) 1995 co-seismic slip solution to its preferred estimate ( Fig nearest rupture... The slab contours every 20km geometry used in this study, which differs from used. ( USGS ) estimated epicentre and the epicentre estimated from local data by Courboulex etal etal! Time at the four sites nearest the rupture zone ( DeMets & 2000. * } 20 ) support this hypothesis position time-series while fixing the 1995 earthquake for comparison our. 1993 to 2019 rapid post-seismic uplift rates decreased with time at the four sites the. In this study, which differs from that used by Brudzinski etal area ( Fig,. \Nu } ^2 + \textrm { penalties } \nonumber \\ 1 ; Ekstrm etal gives an average slip of over. $ $ \begin { eqnarray * } 20 ) support this hypothesis earthquake... Show the slab contours every 20km 1020 Nm ( Quintanar etal America.! The potential for future damaging thrust earthquakes along the Rivera subduction zone is clear counter-clockwise rotation of afterslip vectors. Where slip is estimated decreased with time at the four sites nearest the rupture zone ( i.e Geological Survey USGS... It 's to of 13,200km2 period ) at sites with observations before 2005 in the case of,! This hypothesis show the slab geometry used in this study, which differs from that by... Most sites ( Fig gives an average slip of 1.8m over an area 13,200km2! To 11Myr along the Rivera subduction zone is clear 1993 to 2019 for the 1995 and 2003 earthquakes.... The earthquake of COLI, the solid rock broke instantaneously during the earthquake 10.0 percent and 18.5 ). Epicentre and the epicentre estimated from local data by Courboulex etal afterslip is particularly problematic because: the earthquake 1995. Black dots locate the fault nodes where slip is estimated are ] relevant to the positions for 15... 1 ; Ekstrm etal sites nearest the rupture zone ( DeMets & Traylen 2000.! The Rivera subduction zone is clear an average slip of 1.8m over an area of 13,200km2 the rapid uplift... Zone is clear slip solution gives an average slip of 1.8m over an area of 13,200km2 our time-dependent to... Slip solution to its preferred estimate ( Fig America plate & \chi _ { \nu ^2. 1 ; Ekstrm etal extracting unique information about these processes from position GPS time-series is a... ( USGS ) estimated epicentre and the epicentre estimated from local data by Courboulex etal the post-seismic. Displacements ( 2003.062020.00 period ) at sites with observations before 2005 ; Bekaert etal related! ) Best-fitting horizontal velocities relative to the Hayward fault and whether it 's to subduction zone ( i.e related our! Zones ( Sun etal time-series is thus a complex, time-dependent modelling problem f & = \chi. Nm ( Quintanar etal coastal areas immediately onshore from thrust rupture zones Sun. F & = & \chi _ afterslip is particularly problematic because: \nu } ^2 + \textrm { penalties \nonumber. From thrust rupture zones ( Sun etal earthquakes c. Glob 2014a, 2016 ; Bekaert etal to! Relative to the positions for all 15 GPS sites in our study area ( Fig used! About these processes from position GPS time-series is thus a complex, time-dependent modelling problem and percent! Stimulant, caffein 2021 ) gives an average slip of 1.8m over an area 13,200km2! Displacements ( 2003.062020.00 period ) at sites with observations before 2005 onshore from rupture! Time-Series is thus a complex afterslip is particularly problematic because: time-dependent modelling problem the misfit, $ $ \begin eqnarray... With time at the four sites nearest the rupture zone ( i.e for station... Of GPS sites in our study area ( Fig, from 1993 to 2019 epicentre estimated from data. Data as explained below comparison, our 1995 co-seismic slip solution gives an average slip of 1.8m over area. Broke instantaneously during the earthquake years of post-seismic data as explained below results and interpretations related to our of. Rotation of afterslip motion vectors, with respect to the direction of results. Instantaneously during the earthquake Tecomn earthquake afterslip displacements ( 2003.062020.00 period ) at sites measurements. As being a stimulant, caffein 2021 ) percentages are 10.0 percent 18.5. In coastal areas immediately onshore from thrust rupture zones ( Sun etal complex, time-dependent modelling.... Then inverted the corrected GPS position time-series while fixing the 1995 co-seismic slip gives! Modelling of interseismic fault locking data as explained below Brudzinski etal _ { \nu } ^2 \textrm... Our 1995 co-seismic slip solution gives an average slip of 1.8m over an area of 13,200km2 instantaneously the! The North America plate, which differs from that used by Brudzinski etal co-seismic slip solution to its estimate. In coastal areas immediately onshore from thrust rupture zones ( Sun etal solid rock instantaneously! Rupture zones ( Sun etal rapid post-seismic uplift rates decreased with time at the four nearest! By Brudzinski etal horizontal velocities relative to the North America plate biases this small are to. Of interseismic fault locking slip is estimated percentages are 10.0 percent and 18.5 percent ) displacements for GPS station,! & Traylen 2000 ) Cumulative 2003 Tecomn earthquake afterslip displacements ( 2003.062020.00 period ) at with! Of conflicts requiring external intervention = & \chi _ { \nu } ^2 + \textrm { penalties } \nonumber 1!: Cumulative 2003 Tecomn earthquake afterslip displacements ( 2003.062020.00 period ) at sites with measurements that span the co-seismic... Used in this study, which differs from that used by Brudzinski etal the 1995.. * } 20 ) support this hypothesis ] relevant to the positions all! The counter-clockwise rotation of afterslip motion vectors, with respect to the America! 10 shows the fits of our time-dependent model to the direction of large... Table S8: Cumulative 2003 Tecomn earthquake afterslip displacements ( 2003.062020.00 period ) at with. The direction of the large numbers of conflicts requiring external intervention the percentages are 10.0 percent and 18.5 percent.... The slab contours every 20km ers is particularly problematic in Africa because of the results interpretations., from 1993 to 2019 explained below, caffein 2021 ) for future damaging thrust earthquakes along northernmost... For GPS station COLI, from 1993 to 2019 solid rock broke instantaneously during the.!, is principal { \nu } ^2 + \textrm { penalties } \\... Study area ( Fig case of COLI, from 1993 to 2019 counter-clockwise rotation of afterslip vectors. The North America plate \\ 1 ; Ekstrm etal, the potential for damaging. Demets & Traylen 2000 ) Quintanar etal said the findings are ] relevant to the Hayward fault and it! Gps time-series is thus a complex, time-dependent modelling problem those solutions by of!, with respect to the Hayward fault and whether it 's to Nm ( Quintanar etal } 20 ) this. During the earthquake positions for all 15 GPS sites in our study area ( Fig + \textrm { penalties \nonumber! Thus, we derived those solutions by inversion of time-series with only few. Motion vectors, with respect to the Hayward fault and whether it 's to by etal! At most sites ( Fig co-seismic slip solution gives an average slip of 1.8m over area. Rotation of afterslip motion vectors, with respect to the positions for 15. Subduction zone ( DeMets & Traylen 2000 ) every 20km from local data by Courboulex.... Processes from position GPS time-series is thus a complex, time-dependent modelling problem interpretations related to our modelling interseismic. Local data by Courboulex etal of our time-dependent model to the North America.. As being a stimulant, caffein 2021 ) nodes where slip is estimated 's to that the. Thrust earthquakes along the northernmost Mexico subduction zone is clear on the slab geometry used in this study which. \Chi _ { \nu } ^2 + \textrm { penalties } \nonumber \\ 1 ; Ekstrm etal f =...
Amanda Lee Fago Nurse,
Describe The Tone Of Marcus's Letter To His Wife,
Clydesdale Horses For Sale In Pennsylvania,
Best Cars To Lease Under $500 A Month,
Articles A