Modern Kinematics of the Cordillera Blanca Fault, Northern-Central Peru, from GPS Geodesy

Abstract:

The Cordillera Blanca batholith is an intrusive emplacement of Miocene age (5 Ma) that lies over the northern-central Peruvian Andes. Its western flank is bounded by a normal fault that extends for about 220 km in a NW-SE direction and dips west at a low angle. There is no historical evidence of large earthquakes associated with this fault, but geological studies indicate Quaternary or Holocene activity with as much as 2.5 m of vertical slip. Other related events in the neighborhood of the Cordillera Blanca Fault are the 1942 Quiches (M 6) earthquake and the 1970 Chimbote earthquake (Mw 8). The Chimbote event caused the highest death toll of all times in Peru (70,000 casualties) as it triggered a massive mudslide that swept all along its path and buried the town of Yungay. The ice-and-snow covered, elevated peaks (up to 18,000 ft.) in the Cordillera Banca mountains, as well as the existance of the active bounding normal fault system, represent a big hazard to this region. This study presents the first geodetic estimates of the modern kinematics of the Cordillera Blanca Fault and uses GPS measurements carried out between 2014 and 2017, as well as baseline changes observed between two cGPS stations across its southern end.

Slidecast:

https://vimeo.com/278051014

Transportable Array Alaska – Power and Telemetry System Design and Performance

Abstract:

The Alaska Transportable Array (ATA) is a broadband seismic network made up of 280 new and existing stations that uniformly cover Alaska and north-west Canada. Over five field season from 2013-2017, 194 new stations were installed. In order to meet the performance requirements for the network, new power systems and telemetry systems were designed and produced on a large scale. The majority of ATA stations are installed in completely off grid locations with no access to an established electrical grid or telecommunications network. These off grid stations are fully autonomous and use solar panels and batteries for power and satellite modems to transmit data in real time. The stations also integrate additional instrument packages that include pressure sensors, infrasound sensors and weather stations. ATA stations use a dual battery chemistry power system made up of rechargeable lithium batteries and sealed lead acid batteries to minimize the system’s weight and maximize its energy storage. Savings in weight and volume significantly reduce the cost and complexity of the logistics needed to move the equipment from where it’s fabricated to where it operates. The autonomous ATA stations use high bandwidth satellite modems on the Inmarsat satellite network to transmit station data in near real time. The modems are configured to transmit station data in real time during the summer when incoming solar energy can be harvested with solar panels to keep the station batteries charged. In the winter months, the modems are operated in a lower power mode to maximize battery life yet still deliver data with less than 1 hour of latency. Additionally, autonomous ATA stations use very low power and low bandwidth Iridium modems to provide power system data as well as backup station state of health data in the event of a failure of the primary satellite modem. This state of health data helps operators troubleshoot the cause of outages and improves station service planning and success.

Slidecast:

https://vimeo.com/278026574

Uncertainty Estimation of Moment Tensor Source Types

Abstract:

A moment tensor is a symmetric matrix that expresses the source for a seismic event. Uncertainty characterization of moment tensors is vital for any interpretation about the moment tensor, such as whether an event is likely to be an earthquake or not. We provide a method for characterizing and visualizing the uncertainty for a full moment tensor. Our uncertainty summary has four components: (1) variation in waveform misfit for the best-fitting moment tensor at each lune point; (2) probability density p(v,w) for moment tensor source type; (3) confidence curve Pcon(V); (4) confidence parameter Pav, which is the area under the confidence curve, with large values representing high concentration of probability near the best-fitting moment tensor. These characterizations are facilitated by a uniform parameterization of moment tensors with fixed magnitude. The parameters describe the source type, with v and w, and the orientation with strike angle, cosine dip angle, and rake angle. The parameterization is uniform in the sense that a uniform distribution of 5-tuples in the coordinate domain corresponds to a uniform distribution in the moment tensor space. Stated otherwise, volumes (i.e. 5-volumes) in the coordinate domain are proportional to the corresponding volumes of moment tensors. Uncertainty in source type and be derived from the probability density p(v,w) in source type, which is depicted easily on the vw rectangular domain. We discuss how this approach could be applied to event screening.

Slidecast:

https://vimeo.com/278017773

Rayleigh Wave Group-Velocity across the Dominican Republic and Puerto Rico from Ambient Noise Tomography

Abstract:

The eastern North America-Caribbean (NA-CAR) plate boundary near the islands of Hispaniola and Puerto Rico is a complex transition zone in which strain is accommodated by two transform fault systems and oblique subduction. In 2013, scientists from Baylor University, the Autonomous University of Santo Domingo, and the Puerto Rico Seismic Network deployed 16 broadband stations on the Dominican Republic to expand the local permanent network. The goal of the Greater Antilles Seismic Program (GrASP) is to combine its data with that from permanent networks in Puerto Rico, Haiti, Cuba, the Cayman Islands, and Jamaica to develop a better understanding of the crust and upper mantle structure in the Northeastern Caribbean (Greater Antilles). One important goal of GrASP is to develop robust velocity models that can be used to improve earthquake location and seismic hazard efforts. In this study, we focus on obtaining Rayleigh wave group velocity maps from ambient noise tomography. Here we use vertical-component broadband data recorded at 53 stations between 2010 to present, to obtain Green’s functions between 1165 pairs of stations. From these, we obtain dispersion curves by the application of FTAN methods with phase-matched filtering. Group velocity maps are generated between 4 to 40 s. One-dimensional shear wave velocity models are generated for selected sites by convolving group velocity results with depth kernels produced by forward modeling. The Markov Chain Monte Carlo method is applied to estimate the maximum likelihood S wave velocity model and, specifically, Moho depth. Results show strong correlations with large-scale geological and tectonic features for periods between 4 – 40 s, such as the Cordillera Central in both the Dominican Republic and Puerto Rico, the Mona Passage, and the NA-CAR subduction zone. Preliminary 1D shear wave velocity models suggest that the Moho underneath the Dominican Republic and Puerto Rico is in the range of 28 to 32 km.

Slidecast:

https://vimeo.com/278019459

Seismo-Geodetic Monitoring of the Marmara Seismic Gap

Abstract:

The North Anatolian Fault Zone in the Sea of Marmara did not generate a M>7.0 earthquake since 1766. This fault section stores ~20 mm annual slip deficit and therefore is expected to accommodate at least one at most three M>7.0 in near future. In this study, we continuously monitor this critically strained fault section using seismo-geodetic stations that are equipped with 100Hz sampling seismographs and 1Hz sampling GPS recorders. This configuration allows covering a broad spectral band and is sensitive to both fast/slow tectonic motions at large/small temporal and spatial scale, from milliseconds to years, from centimeters to tens of kilometers. Therefore, recorded seismo-geodetic data will be used to identify (1) along-fault variation of the slip deficit, (2) fault segmentation, (3) interaction between slip-deficit and background seismicity, (4) pre-seismic seismo-geodetic symptoms and (5) co-seismic slip in case of M 7.0 earthquake(s).

Slidecast:

https://vimeo.com/278051069

Geodetic Characterization of Kinematic Segmentation Along the Presqu’ile Du Sud / Enriquillo Fault System in Southern Haiti

Abstract:

The ~3 M inhabitant capital region of Haiti, severely affected by the devastating January 12, 2010, M7.0 earthquake, continues to expand at a fast rate. Accurate characterization of regional earthquake sources is key to inform urban development and construction practices through improved regional seismic hazard estimates. We recently showed, using Global Positioning System (GPS) measurements and simple elastic models, that seismogenic strain accumulation in southern Haiti involves an overlooked component of shortening on a south-dipping reverse fault along the southern edge of the Cul-de-Sac basin at 9+-3 mm/yr, in addition to the well-known component of left-lateral strike-slip motion at 6+-2 mm/yr (Symithe and Calais, Tectonophysics, 2016; Calais et al., C.R. Geosciences, 2016). Here we use a significantly improved regional GPS velocity field to revisit the issue of kinematic segmentation in southern Haiti, with a transition from purely strike-slip motion on the quasi-vertical Presqu’Ile du Sud fault, into highly oblique motion on a south-dipping fault along the southern edge of the Cul-de-Sac basin. Reverse faulting on such a south-dipping plate-boundary fault implies ground shaking in Port-au-Prince twice stronger than expected if the major fault was purely strike-slip, as currently assumed in the seismic hazard map integrated into the Haiti National Building Code.

Slidecast:

https://vimeo.com/278050961