Coupling on the MHT - Supplementary material
|
|
- Sophia Paul
- 5 years ago
- Views:
Transcription
1 Coupling on the MHT - Supplementary material February, S. Map of the GPS network Figure S: Map indicating the names and position of the continuous GPS stations used in this study to determine the pattern of coupling on the MHT. S. Derivation of secular velocities from the GPS time series S.. Noise model for the inversion Several studies of GPS time series have established that the daily estimates of GPS positions are temporally correlated (Langbein and Johnson, 997; Zhang et al., 997; Mao et al., 999;
2 Williams, 3a; Williams et al., 4). Assuming a purely white noise model is therefore incorrect and although it doesn t affect much the value of the final parameters inverted for, it results in a dramatic underestimation of their uncertainties. We thus add to the white noise in our GPS time series a component of colored noise, i.e. a noise that has a power spectrum of the form: P κ (f) f κ, (S) where f is the temporal frequency, and κ is called the spectral index (Mandelbrot and Van Ness, 968). The spectral index is estimated for each time series by first fitting equation () (see main paper) to the time series assuming a white noise, and computing a periodogram of the fit s residuals. The spectral index is then estimated by fitting to the power spectrum a combination of white and colored noise (figure S): P (f) = P + P c f κ, (S) where P and P c are the respective amplitudes of the white and the colored noise. Once the spectral index is estimated, we build the covariance matrix of the data as the sum of white and colored noise covariance matrices. The relative amplitudes of both noises are estimated by a Maximum Likelihood Estimation (MLE) method (Williams et al., 4). The covariance matrix for the white noise is the usual diagonal covariance matrix C w = diag(σ, σ,..., σ n), where σ i is the standard deviation of data point number i. The colored noise covariance matrix C κ is built following an adaptation of the method described in Williams (3a): C κ = t s TT T, where t s is the sampling interval (so t s = day for GPS time series), and the matrix T is defined as: ψ... ψ ψ... T = ψ ψ ψ..., (S4) ψ n ψ n ψ n... ψ where the quantities ψ n are defined by the recurrence formula: { ψ = ψ n+ = n κ/ n+ ψ n. The rows and columns corresponding to times with no data are then removed from the covariance matrix. The final data covariance matrix is given by (S3) C D = a C w + b C κ,. (S)
3 SIMR north! Residuals power spectrum DAMA east! Residuals power spectrum Spectral index:! =!.47 Spectral index:! =!.87! (a) SIMR, north component JMSM north! Residuals power spectrum Spectral index:! =!.93 Spectral index:! =!.6 (c) ODRE, north component (d) JMSM, north component IISC north! Residuals power spectrum CHLM north! Residuals power spectrum f (yr!) f (yr )! Spectral index:! =!.6 Spectral index:! =!.78 (b) DAMA, east component ODRE north! Residuals power spectrum f (yr!) f (yr ) f (yr!) f (yr!) (e) IISC, north component (f) CHLM, north component Figure S: Power spectra of the noise (blue curve) for a sample of time series and fit assuming a combination of white noise and colored noise (red curve). The spectral index κ of the colored noise is indicated on each plot. where a and b are the parameters to be estimated by MLE, measuring respectively the amplitude of white and colored noise. Assuming a Gaussian distribution of the uncertainties on GPS positions, the likelihood that has to be maximized with respect to a and b is then likelihood(cd ) = rt CD r e, (π)n/ (det CD )/ (S6) where r is vector of residuals of the fit and N is the number of daily GPS positions available. The fit and residuals on some time series are shown on figure S3. 3
4 SIMR! North displacement and fit DAMA! East displacement and fit!r = !r =. 998 SIMR! North displacement residuals ODRE! North displacement and fit 8 JMSM! North displacement and fit 6 (b) DAMA, east component!r =.!r = ODRE! North displacement residuals JMSM! North displacement residuals 4 (a) SIMR, north component (c) ODRE, north component (d) JMSM, north component IISC! North displacement and fit CHLM! North displacement and fit 4 3!r =.44!r =.68 6 IISC! North displacement residuals 8 9 CHLM! North displacement residuals DAMA! East displacement residuals (e) IISC, north component (f) CHLM, north component Figure S3: Fits and residuals of the continuous GPS time series at some stations. For each time series, the upper plot shows the raw data (blue curve) and the fit (green curve) with equation (). The value of the reduced chi square χr of each fit is indicated on the upper plot s lower right corner of the corresponding figure. S.. Uncertainties due to unmodeled steps in the time series Steps in the time series can be of many different origins, being actually tectonic, environmental or coming from equipment malfunction, human error, etc. (Williams, 3b). The 4
5 ones large enough to be detected are included in the model (equation ()), but smaller ones remain unnoticed and affect the estimates of model parameters and their uncertainties. Therefore, those uncertainties have to be adjusted accordingly. For convenience, we will assume that those unmodeled steps account for all the errors on the model. Those steps are assumed to happen at a frequency ν, and to have a random Gaussian amplitude N (, σx). The standard deviation on the secular velocity due to those steps is then (Williams, 3b) σ v = σ x ν, (S7) T where T is the length of the time series. In the case of our GPS time series, the amplitude of the steps that were actually detected was always greater than. time the median value of the uncertainties on the daily positions in the time series. We hence take σ x = σ D, where. denotes the median value and σ D is the uncertainty on daily positions of the time series. We estimate ν through the following considerations. First, the steps that were large enough to be detected in the time series happened on average once every years. Assuming that the smaller the steps are, the more frequent they would be, the value for ν should be greater than /. On the other hand, a value of ν overestimated (ν > in this case) results in larger uncertainties on the secular velocity, and eventually leads to values of a reduced chi square smaller than when one fits the Euler pole of the Indian plate in the ITRF reference frame (see section 3.3), indicating that the uncertainties on the GPS velocities are probably overestimated. As a result, we chose a value of ν = /3, which gives the final formula for the uncertainties on the secular velocity due to unmodeled steps: σ v = σ D 3T. (S8) The velocities and corresponding uncertainties that we obtain at the GPS stations used in this paper as well as at the DORIS stations COLA and EVEB are given in table S. S.3 Slip resolution and Laplacian The result of our inversion should be assessed in view of its resolution. This information is contained in the resolution matrix: R = ( G T C d G T + Λ T Λ ) G T C d G T, (S9) where G is the Green s matrix defined in equation () from the main paper, C d is the data covariance matrix and Λ is the Laplacian matrix. The diagonal of R tells how well the slip value on each patch can be retrieved by the inversion. However, it doesn t express how each patch correlates with its neighbors. This information is contained in each of the individual columns of R: column number i is the vector of parameters (i.e. the slip on each patch) returned by the inversion from an input dataset corresponding to a unit slip on patch i and no slip on other patches. Usually, what the inversion returns is slip on a more or less
6 spread area centered on patch i. The characteristic size of this area is estimated by fitting a bell curve to the slip on the patches as a function of distance to patch i (Lohman, 4), and taking the standard deviation of that bell curve. Namely, for each patch i, we find the distance w i that minimizes the quantity: N p ( χ i = j= R ji R ii ) d ij e w i, (S) where N p is the number of patches on the fault, R ji is the value of the coefficient (j, i) of the resolution matrix R (row j and column i), and d ij is the distance between patches i and j. This idea of an estimate of the resolution scale on each patch is also used in order to more efficiently smooth our model by weighting the Laplacian according to the resolution on each patch. Since the Laplacian matrix is not yet available (this is what we try to determine), we compute a first resolution matrix using the Moore-Penrose pseudoinverse matrix (Aster et al., ), keeping only the singular values larger than % of the maximum one. We then compute how far each patch correlates with its neighbors with the method previously described applied to this resolution matrix. Finally, each line of the Laplacian matrix is weighted by the decimal logarithm of the resolution size on the corresponding patch. S.4 Supplementary figures on the pattern of coupling on the MHT S.4. Laplacian smoothing On figure S4 we test how different values of the Laplacian smoothing affect the estimate of the moment deficit accumulated every year. Weights assigned to the Laplacian too small Station ID Site name Latitude ( N) Longitude ( E) Elevation (m) BAN Bangalore DGAR Diego Garcia GUAO Guao GUAM Guam Observatory HYDE Hyderabad IISC Indian Inst. Science KUNM Kunming LHAS Lhasa LHAZ Lhasa POL Poligan IVTAN SELE Selezaschita TAIW Taipei URUM Urumqi WUHN Wuhan Table S: List of IGS sites included in the daily regional processing. 6
7 Velocities in ITRF (mm/yr) Time of operation Station lon ( E) lat ( N) V e V n V u Init. End DAMA ± ± ±.37 Nov 997 current GUMB ± ±.34.7 ±.9 Nov 997 current SIMR ± ±.8.68 ±.7 Nov 997 Apr. BRNG ± ± ± 3. Mar 4 May 9 BRN ± ± ±.48 May 9 current CHLM ± ± ±. Mar 4 current JMSM ± ± ±.36 Ma. 4 current KKN ± ±.4.3 ±.3 Jan 4 current KLDN ± ± ±.7 Apr 4 current MSTG ± ±..79 ± 3.7 Apr 4 Sept 4 a MST ± ±.3.8 ± 4.9 Oct 9 current ODRE ± ±.37. ±. Mar 4 current SIM ± ±..4 ±.8 Mar 4 current SRGK ± ±.9 4. ± 3. Mar Feb 7 TPLJ ±. 3.4 ±.3. ±.3 Mar 4 current BMCL ± ± ±.9 Mar 7 current DLPA ±.63 ±..3 ±.68 May 7 current GRHI ± ± ±.64 May 7 current JMLA ± ±.4. ±.6 May 7 current NPGJ ± ±.74.6 ±.63 May 7 current BYNA ±.38 6 ± ±.39 May 8 current DNGD ± ±.8.48 ±.8 May 8 current DRCL ± ±..64 ±.34 Mar 8 current GNTW ± ±.63.7 ± 3. Apr 8 current RMJT ± ±.48.4 ± 4.36 Oct 8 current RMTE ± ± ±. Sep 8 current SMKT ±.7.84 ± ±.83 May 8 current SYBC ±.8.93 ± ±.88 Oct 8 current CUOM ±.9.38 ±..3 ±. Oct 6 current JRGR ±.9.39 ±. 3.3 ±.7 Mar 7 current XGBA ±. 8. ± ± 4. Mar 7 Sep 7 b YARE ± ± ±.86 Oct 6 current ZHXZ ± ±.4.87 ±.88 Oct 6 current MALD ± ±.49. ±.93 Jul 999 May 6 HYDE ± ± ±.8 Sept current IISC ± ±.3. ±.6 Oct 997 current COLA ± ± ± 3.78 Jan 993 Sep 4 EVEB ± ± ± 3.64 May 993 current Table S: Estimates of the secular velocity at the continuous GPS stations in ITRF and dates of operation of each station. The uncertainties on the velocities indicated are the -σ uncertainties. See text for details on the derivation of those quantities. Gaps in the time series are not unfrequent, and one should keep in mind that they are not indicated in this table. a A -day campaign measurement has also been done with a different antenna on the station s monument in October 9. b 4 additional points in May 9 made the positions at this station exploitable. 7
8 (λ <.8) lead to models featuring locked patches only underneath data points, right next to creeping patches. Besides being unphysical and resulting in very high reduced chi squares, such models are highly dependent on the data spatial distribution and must then be rejected. A smoothing too large (λ > ) tends to lead to a fault locked further at depth, and with a very smooth locked-creeping transition, which doesn t fit the data anymore (reduced χ > 3 on figure S4). Within the range of Laplacian weight.8 < λ <, the moment deficit accumulated each year remains within the uncertainties determined by the inversion. Moment deficit rate (Nm/yr) 7. x " continuous GPS " campaign GPS " leveling " total Moment deficit rate " of different data sets Laplacian weight (! parameter) Figure S4: Variation of the χ of the fit and of the moment deficit rate for different values of the weight attributed to the Laplacian in the inversion. The black curve shows the moment deficit accumulated every year as a function of the weight attributed to the Laplacian. The dashed black line and grey shaded area represent the rate of moment deficit with uncertainties derived in this study, i.e M = 6.6 ±.4 9 Nm/yr. The green, red and blue curves respectively represent the value of the χ of the fit to the continuous GPS, campaign and leveling data. S.4. Direction of extension of the Tibetan plateau Figure S shows the sensitivity of the long term velocity and the moment accumulation rate estimated in this study to the direction chosen for the extension of the Tibetan Plateau. The direction N98.E has been chosen because it is the one that affects the least the estimates of the long term velocities (it is the most perpendicular direction to those velocities, i.e. it is the direction onto which the sum of the projections of the East and West long term velocities 8
9 reaches a minimum). But there is no real reason to prevent this direction from varying by a few degrees from the N98.E azimuth. Figure S shows that even by changing this direction by, the final values of the parameters remain within their estimated uncertainties. Far velocity (mm/yr) v west v east 9 9 x 9 Moment deficit (Nm/yr) Azimuth of extension (N o E) Figure S: Impact of the azimuth selected for the extension of the Tibetan plateau on the long term East and West velocities (upper plot) and the moment deficit rate (lower plot). The solid thick lines represent the values of the parameters with respect the azimuth, the horizontal dashed lines and filled area of corresponding colors are the values with - σ uncertainties that we retained in this study (corresponding to an azimuth of N98.E): V e = 7.8 ±. mm/yr, V w =. ± mm/yr and M = 6.6 ±.4 9 Nm/yr. S.4.3 Recurrence time of large eartquakes Assuming that the moment deficit of M = 6.6 ±.4 9 Nm/yr computed in the main paper was released through earthquakes following a Gutenberg-Richter distribution up to a maximum magnitude, the recurrence time of those largest earthquakes (corresponding to those largest magnitudes) is plotted on figure S6(a). The black lines (solid, dashed and dotted) correspond to M = Nm/yr, with different percentages of this moment deficit being released seismically, while the grey surrounding lines show the extent corresponding to the uncertainties on M. This plot shows that earthquakes as large as the 9 Assam earthquake, whose moment magnitude is estimated at M w 8. (Ambraseys and Douglas, 4; Chen and Molnar, 977), could happen as often as once every 7 years within the borders of Nepal. As far as frequency is concerned, this would be the worst case 9
10 scenario where all the moment deficit accumulated was released seismically in earthquakes whose magnitude wouldn t exceed 8.. However, too many parameters remain unknown to make any accurate estimation on the return period of major earthquakes. Should the actual b-value of the seismicity distribution in Nepal slightly differ from, equation (8) shows that those estimates would be significantly affected. Another unknown parameter is the largest possible earthquake magnitude in Nepal which has a paramount effect, as shown on figure S6(b). Indeed if the seismicity on the MHT doesn t go beyond those M w 8. earthquakes, they would indeed have a period of return of about 7 years. But if we assume that the MHT can produce earthquakes up to M w 9., then the return period of M w 8. earthquakes would become of the order of 6 years. Recurrence period (yr)! =.! =.8! = Maximum M w of earthquakes (a) Recurrence time of the largest possible earthquakes in Nepal assuming a release of a proportion α of the accumulated moment by a seismicity following a Gutenberg-Richter distribution with b =. 9 T REC (M w = 8.) (yr) Maximum M w of earthquakes (b) Reccurence time of M w 8. earthquakes as a function of the largest possible earthquakes happening in Nepal, for a moment accumulation of M = Nm/yr released entirely seismically. Figure S6: Estimations on the recurrence time of earthquakes.
11 References Ambraseys, J. J., and J. Douglas, Magnitude calibration of north indian earthquakes, Geophysical Journal International, 9, 6 6, 4. Aster, R. C., C. H. Thurber, and B. Borchers, Parameter Estimation and Inverse Problems, Academic Press,. Chen, W.-P., and P. Molnar, Seismic moments of major earthquakes and the average rate of slip in central asia, Journal of Geophysical Research, 8(), 94969, 977. Langbein, J., and H. Johnson, Correlated errors in geodetic time series: Implications for time-dependent deformation, Journal of Geophysical Research,, 9 63, 997. Lohman, R. B., The inversion of geodetic data for earthquake parameters, Ph.D. thesis, California Institute of Technology, 4. Mandelbrot, B., and J. Van Ness, Fractional brownian motions, fractional noises, and applications, SIAM Rev.,, 4 439, 968. Mao, A., C. G. A. Harrison, and T. H. Dixon, Noise in gps coordinate time series, Journal of Geophysical Research, 4, , 999. Williams, S. D., Y. Bock, P. Fang, P. Jamason, R. M. Nikolaidis, L. Prawirodirdjo, M. Miller, and D. J. Johnson, Error analysis of continuous gps position time series, Journal of Geophysical Research, 9, 4. Williams, S. D. P., The effect of coloured noise on the uncertainties of rates estimated from geodetic time series, Journal of Geodesy, 76, , 3a. Williams, S. D. P., Offsets in global positioning system time series, Journal of Geophysical Research, 8, doi:.9/jb,6, 3b. Zhang, J., Y. Bock, H. Johnson, P. Fang, S. Williams, J. Genrich, S. Wdowinski, and J. Behr, Southern california permanent gps geodetic array: Error analysis of daily position estimates and site velocities, J. Geophys. Res.,, 8,3 8,, 997.
Lab 9: Satellite Geodesy (35 points)
Lab 9: Satellite Geodesy (35 points) Here you will work with GPS Time Series data to explore plate motion and deformation in California. This lab modifies an exercise found here: http://www.unavco.org:8080/cws/pbonucleus/draftresources/sanandreas/
More informationNoise Characteristics in High Precision GPS Positioning
Noise Characteristics in High Precision GPS Positioning A.R. Amiri-Simkooei, C.C.J.M. Tiberius, P.J.G. Teunissen, Delft Institute of Earth Observation and Space systems (DEOS), Delft University of Technology,
More informationThe problem (1/2) GPS velocity fields in plate boundary zones are very smooth. What does this smoothness hide?
Block models The problem (1/2) GPS velocity fields in plate boundary zones are very smooth Figure from Tom Herring, MIT What does this smoothness hide? Continuous deformation? Rigid block motions, with
More informationSeasonal and long-term vertical deformation in the Nepal Himalaya constrained by GPS and GRACE measurements
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 117,, doi:10.1029/2011jb008925, 2012 Seasonal and long-term vertical deformation in the Nepal Himalaya constrained by GPS and GRACE measurements Yuning Fu 1 and Jeffrey
More informationto: Interseismic strain accumulation and the earthquake potential on the southern San
Supplementary material to: Interseismic strain accumulation and the earthquake potential on the southern San Andreas fault system by Yuri Fialko Methods The San Bernardino-Coachella Valley segment of the
More informationTransient signal detection using GPS measurements: Transient inflation at Akutan volcano, Alaska, during early 2008
GEOPHYSICAL RESEARCH LETTERS, VOL. 38,, doi:10.1029/2011gl046904, 2011 Transient signal detection using GPS measurements: Transient inflation at Akutan volcano, Alaska, during early 2008 Kang Hyeun Ji
More informationJournal of Geophysical Research Letters Supporting Information for
Journal of Geophysical Research Letters Supporting Information for InSAR observations of strain accumulation and fault creep along the Chaman Fault system, Pakistan and Afghanistan H. Fattahi 1, F. Amelung
More informationSupplementary Material
1 Supplementary Material 2 3 4 Interseismic, megathrust earthquakes and seismic swarms along the Chilean subduction zone (38-18 S) 5 6 7 8 9 11 12 13 14 1 GPS data set We combined in a single data set
More informationGPS Strain & Earthquakes Unit 5: 2014 South Napa earthquake GPS strain analysis student exercise
GPS Strain & Earthquakes Unit 5: 2014 South Napa earthquake GPS strain analysis student exercise Strain Analysis Introduction Name: The earthquake cycle can be viewed as a process of slow strain accumulation
More informationGNSS measurements in Ulaanbaatar and its vicinity
GNSS measurements in Ulaanbaatar and its vicinity Erdenezul Danzansan Institute of Astronomy and Geophysics MAS, Ulaanbaatar, Mongolia erdenezul@iag.ac.mn Andrei Miroshnichenko Institute of the Earth s
More informationSUPPLEMENTARY INFORMATION
doi: 1.138/nature962 Data processing A network of 5 continuously recording GPS stations (LAGU, CHIN, ENAP, GUAD and JUAN) was installed after the earthquake (Figure 1, main text). The data considered in
More informationRegional Geodesy. Shimon Wdowinski. MARGINS-RCL Workshop Lithospheric Rupture in the Gulf of California Salton Trough Region. University of Miami
MARGINS-RCL Workshop Lithospheric Rupture in the Gulf of California Salton Trough Region Regional Geodesy Shimon Wdowinski University of Miami Rowena Lohman, Kim Outerbridge, Tom Rockwell, and Gina Schmalze
More informationJournal of Geophysical Research (Solid Earth) Supporting Information for
Journal of Geophysical Research (Solid Earth) Supporting Information for Postseismic Relocking of the Subduction Megathrust Following the 2007 Pisco, Peru earthquake D.Remy (a), H.Perfettini (b), N.Cotte
More informationElectronic supplement for Forearc motion and deformation between El Salvador and Nicaragua: GPS, seismic, structural, and paleomagnetic observations
DR2011053 Electronic supplement for Forearc motion and deformation between El Salvador and Nicaragua: GPS, seismic, structural, and paleomagnetic observations by D. Alvarado et al., Lithosphere, April,
More informationStrain accumulation studies between Antarctica and India through geodetic tying of two continents from continuous GPS measurements
Indian Journal of Marine Sciences Vol. 37(4), December 2008, pp. 404-411 Strain accumulation studies between Antarctica and India through geodetic tying of two continents from continuous GPS measurements
More informationMegathrust earthquakes: How large? How destructive? How often? Jean-Philippe Avouac California Institute of Technology
Megathrust earthquakes: How large? How destructive? How often? Jean-Philippe Avouac California Institute of Technology World seismicity (data source: USGS) and velocities relative to ITRF1997 (Sella et
More informationInquiry: Sumatran earthquakes with GPS Earth Science Education
Inquiry: Sumatran earthquakes with GPS Earth Science Education www.earthobservatory.sg Preparation: Before doing this investigation, complete two introductory investigations using GPS data from UNAVCO
More informationPrinciples of the Global Positioning System Lecture 23
12.540 Principles of the Global Positioning System Lecture 23 Prof. Thomas Herring 05/12/03 12.540 Lec 23 1 OVERVIEW Run through HW 03. Final lecture: Application area and GPS results Analysis of contemporary
More informationSupporting Information for Interseismic Coupling on the Main Frontal Thrust
JOURNAL OF GEOPHYSICAL RESEARCH Supporting Information for Interseismic Coupling on the Main Frontal Thrust DOI:./ V. L. Stevens and JP Avouac, 4 5 6 7 8 9 Contents of this file. Tables S to S.. Figures
More informationGAMINGRE 8/1/ of 7
FYE 09/30/92 JULY 92 0.00 254,550.00 0.00 0 0 0 0 0 0 0 0 0 254,550.00 0.00 0.00 0.00 0.00 254,550.00 AUG 10,616,710.31 5,299.95 845,656.83 84,565.68 61,084.86 23,480.82 339,734.73 135,893.89 67,946.95
More informationActual Continuous Kinematic Model (ACKIM) of the Earth s Crust based on ITRF2014
Actual Continuous Kinematic Model (ACKIM) of the Earth s Crust based on ITRF2014 Hermann Drewes Technische Universität München Joint Scientific Assembly of the International Association of Geodesy and
More informationNational Report of Greece to EUREF 2016
National Report of Greece to EUREF 2016 M.Gianniou, E. Mitropoulou, D. Mastoris National Cadastre and Mapping Agency S.A. Mesogion Ave. 288, 15562 Athens, Greece 1. Introduction During the last months
More informationSeismogeodesy for rapid earthquake and tsunami characterization
Seismogeodesy for rapid earthquake and tsunami characterization Yehuda Bock Scripps Orbit and Permanent Array Center Scripps Institution of Oceanography READI & NOAA-NASA Tsunami Early Warning Projects
More informationPrinciples of the Global Positioning System Lecture 11
12.540 Principles of the Global Positioning System Lecture 11 Prof. Thomas Herring http://geoweb.mit.edu/~tah/12.540 Statistical approach to estimation Summary Look at estimation from statistical point
More informationGPS Strain & Earthquakes Unit 3: Finding location and velocity data for PBO GPS stations
GPS Strain & Earthquakes Unit 3: Finding location and velocity data for PBO GPS stations Original activity by Vince Cronin (Baylor University). Revisions by Beth Pratt-Sitaula (UNAVCO). Analyzing the velocities
More informationVelocities of GPS sites throughout Japan are derived from combination of SINEX
DR2009039 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 Appendix : GPS data analysis and data interpretation GPS data analysis and transient deformation Velocities of GPS sites throughout
More informationFrom Earthquakes to Mountains: the Earth s Crust in Motion
Dean s Freshman Honor Seminar (SCI 110) Eric Calais Associate Professor of Geophysics Purdue University Department of Earth and Atmospheric Sciences ecalais@purdue.edu http:/www.eas.purdue.edu/~calais
More informationESTIMATES OF HORIZONTAL DISPLACEMENTS ASSOCIATED WITH THE 1999 TAIWAN EARTHQUAKE
ESTIMATES OF HORIZONTAL DISPLACEMENTS ASSOCIATED WITH THE 1999 TAIWAN EARTHQUAKE C. C. Chang Department of Surveying and Mapping Engineering Chung Cheng Institute of Technology, Taiwan, ROC ABSTRACT A
More informationA TESTABLE FIVE-YEAR FORECAST OF MODERATE AND LARGE EARTHQUAKES. Yan Y. Kagan 1,David D. Jackson 1, and Yufang Rong 2
Printed: September 1, 2005 A TESTABLE FIVE-YEAR FORECAST OF MODERATE AND LARGE EARTHQUAKES IN SOUTHERN CALIFORNIA BASED ON SMOOTHED SEISMICITY Yan Y. Kagan 1,David D. Jackson 1, and Yufang Rong 2 1 Department
More informationEarthquake distribution is not random: very narrow deforming zones (= plate boundaries) versus large areas with no earthquakes (= rigid plate
Earthquake distribution is not random: very narrow deforming zones (= plate boundaries) versus large areas with no earthquakes (= rigid plate interiors) Tectonic plates and their boundaries today -- continents
More informationStrain accumulation at Yucca Mountain, Nevada,
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 104, NO. B8, PAGES 17,627-17,631, AUGUST 10, 1999 Strain accumulation at Yucca Mountain, Nevada, 1983-1998 J. C. Savage, J.L. Svarc, and W. H. Prescott U.S. Geological
More informationFORMAT FOR TECHNICAL PAPERS
TOWARDS A MODERNIZED GEODETIC DATUM FOR NEPAL: OPTIONS FOR DEVELOPING AN ACCURATE TERRESTRIAL REFERENCE FRAME FOLLOWING THE APRIL 25, 2015 MW7.8 GORKHA EARTHQUAKE Chris PEARSON 1, Niraj MANANDHA 2 ABSTRACT
More informationBasics of the modelling of the ground deformations produced by an earthquake. EO Summer School 2014 Frascati August 13 Pierre Briole
Basics of the modelling of the ground deformations produced by an earthquake EO Summer School 2014 Frascati August 13 Pierre Briole Content Earthquakes and faults Examples of SAR interferograms of earthquakes
More informationJackson County 2013 Weather Data
Jackson County 2013 Weather Data 61 Years of Weather Data Recorded at the UF/IFAS Marianna North Florida Research and Education Center Doug Mayo Jackson County Extension Director 1952-2008 Rainfall Data
More informationPlate Boundary Observatory Working Group for the Central and Northern San Andreas Fault System PBO-WG-CNSA
Plate Boundary Observatory Working Group for the Central and Northern San Andreas Fault System PBO-WG-CNSA Introduction Our proposal focuses on the San Andreas fault system in central and northern California.
More informationThe International Terrestrial Reference System and ETRS89: Part I : General concepts
The International Terrestrial Reference System and ETRS89: Part I : General concepts Zuheir ALTAMIMI Laboratoire de Recherche en Géodésie Institut national de l information géographique et forestière (IGN),
More informationSupporting Information for Break of slope in earthquake-size distribution reveals creep rate along the San Andreas fault system
GEOPHYSICAL RESEARCH LETTERS Supporting Information for Break of slope in earthquake-size distribution reveals creep rate along the San Andreas fault system Inessa Vorobieva, 1,2 Peter Shebalin, 1,2 Clément
More informationGeodesy (InSAR, GPS, Gravity) and Big Earthquakes
Geodesy (InSAR, GPS, Gravity) and Big Earthquakes Mathew Pritchard Teh-Ru A. Song Yuri Fialko Luis Rivera Mark Simons UJNR Earthquake Research Panel, Morioka, Japan - Nov 6, 2002 Goals Accurate and high
More informationSolid Earth Deformations Induced by the Sumatra Earthquakes of : GPS Detection of Co-Seismic Displacements
Solid Earth Deformations Induced by the Sumatra Earthquakes of 2-25: GPS Detection of Co-Seismic Displacements and Tsunami-Induced Loading H.-P. Plag, G. Blewitt, C. Kreemer, W.C. Haond Nevada Bureau of
More informationGPS measurements of present-day uplift in the Southern Tibet
LETTER Earth Planets Space, 52, 735 739, 2000 GPS measurements of present-day uplift in the Southern Tibet Caijun Xu, Jingnan Liu, Chenghua Song, Weiping Jiang, and Chuang Shi School of Geoscience and
More informationKinematic Waveform Inversion Study of Regional Earthquakes in Southwest Iberia
Kinematic Waveform Inversion Study of Regional Earthquakes in Southwest Iberia Ana Domingues Under supervision of Prof. João Fonseca and Dr. Susana Custódio Dep. Physics, IST, Lisbon, Portugal November
More informationEstimating crustal deformations using GNSS. Riccardo Barzaghi DICA-Politecnico di Milano
Estimating crustal deformations using GNSS Riccardo Barzaghi DICA-Politecnico di Milano GNSS non-permanent campaigns The local campaigns require a dense GNSS network in a quite small area a careful distribution
More informationAssessing the solution quality of the earthquake location problem
Semesterarbeit HS 2008 Department of Earth Sciences ETH Zurich Assessing the solution quality of the earthquake location problem Sabrina Schönholzer Supervisor: Technical Supervisor: Prof. Dr. E. Kissling
More informationEarthquake Stress Drops in Southern California
Earthquake Stress Drops in Southern California Peter Shearer IGPP/SIO/U.C. San Diego September 11, 2009 Earthquake Research Institute Lots of data for big earthquakes (rupture dimensions, slip history,
More information3.3 Analysis Coordinator
1. Introduction 2. Analysis of ITRF2014 post-seismic parameterization In this report we outline the activities of the Analysis Coordinator during 2015. The main activities were analysis of the ITRF2014
More informationDetermination of Current Velocity Field (Rate) of North Anatolian Fault in Izmit-Sapanca Segment
Determination of Current Velocity Field (Rate) of North Anatolian Fault in Izmit-Sapanca Segment Cetin MEKIK, Bulent TURGUT, Haluk OZENER, Turkish Republic INTRODUCTION Turkey is geologically the part
More informationExercise 6. Solar Panel Orientation EXERCISE OBJECTIVE DISCUSSION OUTLINE. Introduction to the importance of solar panel orientation DISCUSSION
Exercise 6 Solar Panel Orientation EXERCISE OBJECTIVE When you have completed this exercise, you will understand how the solar illumination at any location on Earth varies over the course of a year. You
More informationSurface Anomalies Prior to Earthquakes
Surface Anomalies Prior to Earthquakes Habibeh Valizadeh, Shattri B. Mansor Husaini Omar and Farid Azad Department of Civil Engineering Universiti Putra Malaysia Serdang, Selangor Malaysia shattri@eng.upm.edu.my
More informationVertical to Horizontal (V/H) Ratios for Large Megathrust Subduction Zone Earthquakes
Vertical to Horizontal (V/H) Ratios for Large Megathrust Subduction Zone Earthquakes N.J. Gregor Consultant, Oakland, California, USA N.A. Abrahamson University of California, Berkeley, USA K.O. Addo BC
More informationpeak half-hourly New South Wales
Forecasting long-term peak half-hourly electricity demand for New South Wales Dr Shu Fan B.S., M.S., Ph.D. Professor Rob J Hyndman B.Sc. (Hons), Ph.D., A.Stat. Business & Economic Forecasting Unit Report
More informationThe Size and Duration of the Sumatra-Andaman Earthquake from Far-Field Static Offsets
The Size and Duration of the Sumatra-Andaman Earthquake from Far-Field Static Offsets P. Banerjee, 1 F. F. Pollitz, 2 R. Bürgmann 3 * 1 Wadia Institute of Himalayan Geology, Dehra Dun, 248001, India. 2
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION doi: 10.1038/ngeo739 Supplementary Information to variability and distributed deformation in the Marmara Sea fault system Tobias Hergert 1 and Oliver Heidbach 1,* 1 Geophysical
More informationChiang Rai Province CC Threat overview AAS1109 Mekong ARCC
Chiang Rai Province CC Threat overview AAS1109 Mekong ARCC This threat overview relies on projections of future climate change in the Mekong Basin for the period 2045-2069 compared to a baseline of 1980-2005.
More informationOngoing Oroclinal bending in the Cascadia forearc and its relation to concave-outboard plate margin geometry
GSA DATA REPOSITORY 2019058 Ongoing Oroclinal bending in the Cascadia forearc and its relation to concave-outboard plate margin geometry Theron Finley 1,2, Kristin Morell 1,3, Lucinda Leonard 1, Christine
More informationCase Study of Australia
Case Study of Australia John Dawson, Geoscience Australia 1 Australia s Geodetic eco-system Geocentric Coordinates Observing Infrastructure Height Datum Geoid Models Tools Services Standards 2 Australia
More informationSeparating Tectonic, Magmatic, Hydrological, and Landslide Signals in GPS Measurements near Lake Tahoe, Nevada-California
Separating Tectonic, Magmatic, Hydrological, and Landslide Signals in GPS Measurements near Lake Tahoe, Nevada-California Geoffrey Blewitt, Corné Kreemer, William C. Hammond, & Hans-Peter Plag NV Geodetic
More informationRegional Bias in Estimates of Earthquake Ms Due to Surface-Wave Path Effects
Bulletin of the Seismological Society of America, Vol. 84, No. 2, pp. 377-382, April 1994 Regional Bias in Estimates of Earthquake Ms Due to Surface-Wave Path Effects by Rachel E. Abercrombie Abstract
More informationComparison Figures from the New 22-Year Daily Eddy Dataset (January April 2015)
Comparison Figures from the New 22-Year Daily Eddy Dataset (January 1993 - April 2015) The figures on the following pages were constructed from the new version of the eddy dataset that is available online
More information= (G T G) 1 G T d. m L2
The importance of the Vp/Vs ratio in determining the error propagation and the resolution in linear AVA inversion M. Aleardi, A. Mazzotti Earth Sciences Department, University of Pisa, Italy Introduction.
More informationWidespread Ground Motion Distribution Caused by Rupture Directivity during the 2015 Gorkha, Nepal Earthquake
Widespread Ground Motion Distribution Caused by Rupture Directivity during the 2015 Gorkha, Nepal Earthquake Kazuki Koketsu 1, Hiroe Miyake 2, Srinagesh Davuluri 3 and Soma Nath Sapkota 4 1. Corresponding
More informationRealizing a geodetic reference frame using GNSS in the presence of crustal deformations: The case of Greece
Realizing a geodetic reference frame using GNSS in the presence of crustal deformations: The case of Greece M. Gianniou National Cadastre and Mapping Agency S.A. Outline 1. Introduction 2. Deformation
More informationEstimating fault slip rates, locking distribution, elastic/viscous properites of lithosphere/asthenosphere. Kaj M. Johnson Indiana University
3D Viscoelastic Earthquake Cycle Models Estimating fault slip rates, locking distribution, elastic/viscous properites of lithosphere/asthenosphere Kaj M. Johnson Indiana University In collaboration with:
More informationobservational data analysis (Supplementary Material)
1 Assessing the equatorial long-wave approximation: asymptotics and observational data analysis (Supplementary Material) 3 by H. Reed Ogrosky and Samuel N. Stechmann 4 1. Kelvin and MRG waves 5 6 The analysis
More informationSupplementary appendix
Supplementary appendix This appendix formed part of the original submission and has been peer reviewed. We post it as supplied by the authors. Supplement to: Lowe R, Stewart-Ibarra AM, Petrova D, et al.
More informationKinematics of the Southern California Fault System Constrained by GPS Measurements
Title Page Kinematics of the Southern California Fault System Constrained by GPS Measurements Brendan Meade and Bradford Hager Three basic questions Large historical earthquakes One basic question How
More informationRegional deformation and kinematics from GPS data
Regional deformation and kinematics from GPS data Jessica Murray, Jerry Svarc, Elizabeth Hearn, and Wayne Thatcher U. S. Geological Survey Acknowledgements: Rob McCaffrey, Portland State University UCERF3
More informationDEFORMATION KINEMATICS OF TIBETAN PLATEAU DETERMINED FROM GPS OBSERVATIONS
DEFORMATION KINEMATICS OF TIBETAN PLATEAU DETERMINED FROM GPS OBSERVATIONS Jinwei Ren Institute of Geology, China Seismological Bureau, Beijing 100029 China Tel: (10)62009095; Fax: (10)62009003; email:
More informationSegmentation in episodic tremor and slip all along Cascadia
Segmentation in episodic tremor and slip all along Cascadia Michael R. Brudzinski and Richard M. Allen Geology 35 (10) 907-910, 2007, doi: 10.1130/G23740A.1 Data Repository: Methods for Automated Data
More informationThe effect and correction of non-tectonic crustal deformation for continuous GPS position time series
The effect and correction of non-tectonic crustal deformation for continuous GPS position time series WANG Min 1,2) SHEN Zheng-Kang 1,3) DONG Da-Nan 4) 1) State Key Laboratory of Earthquake Dynamics, Institute
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lecture 16, 9 Nov. 2017 www.geosc.psu.edu/courses/geosc508 Energy Balance of dynamic rupture Crack tip stress field Frictional Rupture Fronts Meet in the lab (522
More informationSolid Earth Deformations Induced by the Sumatra Earthquakes of : GPS Detection of Co-Seismic Displacements
Solid Earth Deformations Induced by the Sumatra Earthquakes of 2-25: GPS Detection of Co-Seismic Displacements and Tsunami-Induced Loading H.-P. Plag, G. Blewitt, C. Kreemer, W.C. Haond Nevada Bureau of
More informationTesting for fault activity at Yucca Mountain, Nevada, using independent GPS results from the BARGEN network
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 33, L14302, doi:10.1029/2006gl026140, 2006 Testing for fault activity at Yucca Mountain, Nevada, using independent GPS results from the BARGEN
More informationPreliminary slip model of M9 Tohoku earthquake from strongmotion stations in Japan - an extreme application of ISOLA code.
Preliminary slip model of M9 Tohoku earthquake from strongmotion stations in Japan - an extreme application of ISOLA code. J. Zahradnik 1), F. Gallovic 1), E. Sokos 2) G-A. Tselentis 2) 1) Charles University
More informationRapid Determination of Earthquake Magnitude using GPS for Tsunami Warning Systems: An Opportunity for IGS to Make a Difference
Rapid Determination of Earthquake Magnitude using GPS for Tsunami Warning Systems: An Opportunity for IGS to Make a Difference Geoffrey Blewitt, 1 Corné Kreemer, 1 William C. Hammond, 1 Hans-Peter Plag,
More informationWhat happened before the last five strong earthquakes in Greece: Facts and open questions
86 Proc. Jpn. Acad., Ser. B 82 (2006) [Vol. 82, What happened before the last five strong earthquakes in Greece: Facts and open questions By Panayiotis A. VAROTSOS ) Solid State Section and Solid Earth
More informationPrinciples of the Global Positioning System Lecture 18" Mathematical models in GPS" Mathematical models used in GPS"
12.540 Principles of the Global Positioning System Lecture 18" Prof. Thomas Herring" Room 54-820A; 253-5941" tah@mit.edu" http://geoweb.mit.edu/~tah/12.540 " Mathematical models in GPS" Review assignment
More informationNew satellite mission for improving the Terrestrial Reference Frame: means and impacts
Fourth Swarm science meeting and geodetic missions workshop ESA, 20-24 March 2017, Banff, Alberta, Canada New satellite mission for improving the Terrestrial Reference Frame: means and impacts Richard
More informationGPS time series and sea level
GPS time series and sea level M. Poutanen 1, H. Koivula 1, M. Tervo 1,2, K. Kahma 3, M. Ollikainen 1, H. Virtanen 1 1 Finnish Geodetic Institute, 2 University of Helsinki, 3 Finnish Institute of Marine
More informationTime Series Analysis
Time Series Analysis A time series is a sequence of observations made: 1) over a continuous time interval, 2) of successive measurements across that interval, 3) using equal spacing between consecutive
More informationPrediction of Great Japan Earthquake 11 March of 2011 by analysis of microseismic noise. Does Japan approach the next Mega-EQ?
Prediction of Great Japan Earthquake 11 March of 2011 by analysis of microseismic noise. Does Japan approach the next Mega-EQ? Alexey Lyubushin lyubushin@yandex.ru, http://alexeylyubushin.narod.ru/ Institute
More informationGEORED Project: GNSS Geodesy Network for Geodynamics Research in Colombia, South America. Héctor Mora-Páez
GEORED Project: GNSS Geodesy Network for Geodynamics Research in Colombia, South America. Héctor Mora-Páez Colombian Geological Survey UNAVCO SCIENCE MEETING Feb 27 March 1, 2012 Boulder, CO, TECTONIC
More informationLecture 14: Strain Examples. GEOS 655 Tectonic Geodesy Jeff Freymueller
Lecture 14: Strain Examples GEOS 655 Tectonic Geodesy Jeff Freymueller A Worked Example Consider this case of pure shear deformation, and two vectors dx 1 and dx 2. How do they rotate? We ll look at vector
More informationLateral extrusion and tectonic escape in Ilan Plain of northeastern Taiwan
Lateral extrusion and tectonic escape in Ilan Plain of northeastern Taiwan Angelier, J., Chang, T.Y., Hu, J.C., Chang, C.P., Siame, L., Lee, J.C., Deffontaines, B., Chu, H.T, Lu, C.Y., Does extrusion occur
More informationIntroduction: the Abruzzo earthquake The network and the processing strategies. displacements estimation at earthquake epoch. horizontal displacements
The Abruzzo earthquake: temporal and spatial analysis of the first geodetic results L. Biagi, S. Caldera, D. Dominici, F. Sansò Politecnico di Milano Università degli studi de L Aquila Outline Introduction:
More information25th Seismic Research Review - Nuclear Explosion Monitoring: Building the Knowledge Base
KARATAU, MAKANCHI, AND LOCATION IN WESTERN CHINA Lee K. Steck, George E. Randall, Michael L. Begnaud, Charlotte A. Rowe, and Aaron A. Velasco Los Alamos National Laboratory Sponsored by National Nuclear
More informationEarthquake. What is it? Can we predict it?
Earthquake What is it? Can we predict it? What is an earthquake? Earthquake is the vibration (shaking) and/or displacement of the ground produced by the sudden release of energy. Rocks under stress accumulate
More informationCodal provisions of seismic hazard in Northeast India
Codal provisions of seismic hazard in Northeast India Sandip Das 1, Vinay K. Gupta 1, * and Ishwer D. Gupta 2 1 Department of Civil Engineering, Indian Institute of Technology, Kanpur 208 016, India 2
More informationON THE RELATION BETWEEN GPS STRAIN FIELD AND ACTIVE FAULTS IN THE EASTERN BALTIC REGION
ON THE RELATION BETWEEN GPS STRAIN FIELD AND ACTIVE FAULTS IN THE EASTERN BALTIC REGION B.A. Assinovskaya 1, V.L. Gorshkov 2, N.V. Shcherbakova 2, N.M. Panas 1 1 Geophysical Service RAS, seismic station
More informationFULL MOMENT TENSOR ANALYSIS USING FIRST MOTION DATA AT THE GEYSERS GEOTHERMAL FIELD
PROCEEDINGS, Thirty-Eighth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 2013 SGP-TR-198 FULL MOMENT TENSOR ANALYSIS USING FIRST MOTION DATA AT
More informationpeak half-hourly South Australia
Forecasting long-term peak half-hourly electricity demand for South Australia Dr Shu Fan B.S., M.S., Ph.D. Professor Rob J Hyndman B.Sc. (Hons), Ph.D., A.Stat. Business & Economic Forecasting Unit Report
More informationErrors in the Data Measurements
Voyager PLS Error Analysis Errors in the Data Measurements For all purposes in this Error Analysis memo, it will be assumed that currents are measured and displayed in femtoamps. This is the natural range
More informationLarge surface wave of the 2004 Sumatra-Andaman earthquake captured by the very long baseline kinematic analysis of 1-Hz GPS data
Earth Planets Space, 58, 153 157, 2006 Large surface wave of the 2004 Sumatra-Andaman earthquake captured by the very long baseline kinematic analysis of 1-Hz GPS data Yusaku Ohta, Irwan Meilano, Takeshi
More informationAutomatic Moment Tensor Analyses, In-Situ Stress Estimation and Temporal Stress Changes at The Geysers EGS Demonstration Project
PROCEEDINGS, 42nd Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 13-15, 2017 SGP-TR-212 Automatic Moment Tensor Analyses, In-Situ Stress Estimation and
More informationGPS study of N-S trending Karaburun Belt (Turkey) and its E-W trending eastern part
Presented at the FIG Working Week 2017, GPS study of N-S trending Karaburun Belt (Turkey) and its E-W trending eastern part May 29 - June 2, 2017 in Helsinki, Finland Muzaffer KAHVECİ Selcuk University,
More informationPractical considerations for determining Euler Pole Parameters for the terrestrial reference frames in the United States
Practical considerations for determining Euler Pole Parameters for the terrestrial reference frames in the United States Jacob M. HECK and Dru A. SMITH, United States Key words: GNSS/GPS, Positioning,
More informationTectonic deformations in Greece and the operation of HEPOS network
Tectonic deformations in Greece and the operation of HEPOS network M. Gianniou KTIMATOLOGIO S.A. (Hellenic Cadastre) Abstract Today, networks of permanent reference stations are broadly used for the realization
More informationEarthquake stress drop estimates: What are they telling us?
Earthquake stress drop estimates: What are they telling us? Peter Shearer IGPP/SIO/U.C. San Diego October 27, 2014 SCEC Community Stress Model Workshop Lots of data for big earthquakes (rupture dimensions,
More informationSimulated and Observed Scaling in Earthquakes Kasey Schultz Physics 219B Final Project December 6, 2013
Simulated and Observed Scaling in Earthquakes Kasey Schultz Physics 219B Final Project December 6, 2013 Abstract Earthquakes do not fit into the class of models we discussed in Physics 219B. Earthquakes
More informationSupporting Online Material for
www.sciencemag.org/cgi/content/full/1131692/dc1 Supporting Online Material for Localized Temporal Change of the Earth s Inner Core Boundary This PDF file includes: Materials and Methods Figs. S1 to S3
More informationInversion of Earthquake Rupture Process:Theory and Applications
Inversion of Earthquake Rupture Process:Theory and Applications Yun-tai CHEN 12 * Yong ZHANG 12 Li-sheng XU 2 1School of the Earth and Space Sciences, Peking University, Beijing 100871 2Institute of Geophysics,
More information