Integration of Arrival-Time Datasets for Consistent Quality Control: A Case Study of Amphibious Experiments along the Middle America Trench

We have integrated waveform and arrival-onset data collected in Costa Rica as part of the National Science Foundation (NSF)-sponsored Costa Rica Seismogenic Zone Experiment (CRSEIZE) and along central Costa Rica and Nicaragua as part of the German SFB 574 program. The five arrays, composed of different sensor types (one-and three-component land and ocean bottom seismometers and hydrophones), were archived using different software packages (Antelope and SEISAN) and were automatically and manually picked using various quality criteria resulting in a disparate set of pick weights. We evaluate pick quality using automated arrival detection and picking algorithm based on the wavelet transform and Akaike information criterion picker. The consistency of the arrival information over various scales provides a basis for assigning a quality to the analyst pick. Approximately 31% of P arrival times and 26% of S times have been classified as high-quality picks (quality 0-1). An additional 21% of P times and 27% of S arrivals are good quality (quality 2-3). The revised quality picks are mapped directly into new pick weights for inversion studies. We explore the effect of new weighting and removal of poor data by relocating hypocenters through a minimum 1D velocity model and conducting double-difference local earthquake tomography (LET). Analysis of the hypocenter relocation and seismic velocity tomography results suggest that using the improved quality determinations have a greater effect on improving sharpness in the velocity images than on the magnitude of hypocentral movement

Source parameters of the great Sumatran megathrust earthquakes of 1797 and 1833 inferred from coral microatolls

Large uplifts and tilts occurred on the Sumatran outer arc islands between 0.5° and 3.3°S during great historical earthquakes in 1797 and 1833, as judged from relative sea level changes recorded by annually banded coral heads. Coral data for these two earthquakes are most complete along a 160-km length of the Mentawai islands between 3.2° and 2°S. Uplift there was as great as 0.8 m in 1797 and 2.8 m in 1833. Uplift in 1797 extended 370 km, between 3.2° and 0.5°S. The pattern and magnitude of uplift imply megathrust ruptures corresponding to moment magnitudes (M_w) in the range 8.5 to 8.7. The region of uplift in 1833 ranges from 2° to at least 3.2°S and, judging from historical reports of shaking and tsunamis, perhaps as far as 5°S. The patterns and magnitude of uplift and tilt in 1833 are similar to those experienced farther north, between 0.5° and 3°N, during the giant Nias-Simeulue megathrust earthquake of 2005; the outer arc islands rose as much as 3 m and tilted toward the mainland. Elastic dislocation forward modeling of the coral data yields megathrust ruptures with moment magnitudes ranging from 8.6 to 8.9. Sparse accounts at Padang, along the mainland west coast at latitude 1°S, imply tsunami runups of at least 5 m in 1797 and 3–4 m in 1833. Tsunamis simulated from the pattern of coral uplift are roughly consistent with these reports. The tsunami modeling further indicates that the Indian Ocean tsunamis of both 1797 and 1833, unlike that of 2004, were directed mainly south of the Indian subcontinent. Between about 0.7° and 2.1°S, the lack of vintage 1797 and 1833 coral heads in the intertidal zone demonstrates that interseismic submergence has now nearly equals coseismic emergence that accompanied those earthquakes. The interseismic strains accumulated along this reach of the megathrust have thus approached or exceeded the levels relieved in 1797 and 1833

Interplate seismicity at the CRISP drilling site: The 2002 Mw=6.4 Osa Earthquake at the southeastern end of the Middle America Trench

We investigate potential relations between variations in seafloor relief and age of the incoming plate and interplate seismicity. Westward from Osa Peninsula in Costa Rica, a major change in the character of the incoming Cocos Plate is displayed by abrupt lateral variations in seafloor depth and thermal structure. Here a Mw 6.4 thrust earthquake was followed by three aftershock clusters in June 2002. Initial relocations indicate that the main shock occurred fairly trenchward of most large earthquakes along the Middle America Trench off central Costa Rica. The earthquake sequence occurred while a temporary network of OBH and land stations ∼80 km to the northwest were deployed. By adding readings from permanent local stations, we obtain uncommon P wave coverage of a large subduction zone earthquake. We relocate this catalog using a nonlinear probabilistic approach within both, a 1-D and a 3-D P wave velocity models. The main shock occurred ∼25 km from the trench and probably along the plate interface at 5–10 km depth. We analyze teleseismic data to further constrain the rupture process of the main shock. The best depth estimates indicate that most of the seismic energy was radiated at shallow depth below the continental slope, supporting the nucleation of the Osa earthquake at ∼6 km depth. The location and depth coincide with the plate boundary imaged in prestack depth-migrated reflection lines shot near the nucleation area. Aftershocks propagated downdip to the area of a 1999 Mw 6.9 sequence and partially overlapped it. The results indicate that underthrusting of the young and buoyant Cocos Ridge has created conditions for interplate seismogenesis shallower and closer to the trench axis than elsewhere along the central Costa Rica margin

Audit and feedback and clinical practice guideline adherence: Making feedback actionable

BACKGROUND: As a strategy for improving clinical practice guideline (CPG) adherence, audit and feedback (A&F) has been found to be variably effective, yet A&F research has not investigated the impact of feedback characteristics on its effectiveness. This paper explores how high performing facilities (HPF) and low performing facilities (LPF) differ in the way they use clinical audit data for feedback purposes. METHOD: Descriptive, qualitative, cross-sectional study of a purposeful sample of six Veterans Affairs Medical Centers (VAMCs) with high and low adherence to six CPGs, as measured by external chart review audits. One-hundred and two employees involved with outpatient CPG implementation across the six facilities participated in one-hour semi-structured interviews where they discussed strategies, facilitators and barriers to implementing CPGs. Interviews were analyzed using techniques from the grounded theory method. RESULTS: High performers provided timely, individualized, non-punitive feedback to providers, whereas low performers were more variable in their timeliness and non-punitiveness and relied on more standardized, facility-level reports. The concept of actionable feedback emerged as the core category from the data, around which timeliness, individualization, non-punitiveness, and customizability can be hierarchically ordered. CONCLUSION: Facilities with a successful record of guideline adherence tend to deliver more timely, individualized and non-punitive feedback to providers about their adherence than facilities with a poor record of guideline adherence. Consistent with findings from organizational research, feedback intervention characteristics may influence the feedback's effectiveness at changing desired behaviors

Seismogenic zone structure of the southern Middle America Trench, Costa Rica

The shallow seismogenic portion of subduction zones generates damaging large and great earthquakes. This study provides structural constraints on the seismogenic zone of the Middle America Trench offshore central Costa Rica and insights into the physical and mechanical characteristics controlling seismogenesis. We have located ~300 events that occurred following the MW 6.9, 20 August 1999, Quepos, Costa Rica, underthrusting earthquake using a three-dimensional velocity model and arrival time data recorded by a temporary local network of land and ocean bottom seismometers. We use aftershock locations to define the geometry and characteristics of the seismogenic zone in this region. These events define a plane dipping at 19° that marks the interface between the Cocos Plate and the Panama Block. The majority of aftershocks occur below 10 km and above 30 km depth below sea level, corresponding to 30–35 km and 95 km from the trench axis, respectively. Relative event relocation produces a seismicity pattern similar to that obtained using absolute locations, increasing confidence in the geometry of the seismogenic zone. The aftershock locations spatially correlate with the downdip extension of the oceanic Quepos Plateau and reflect the structure of the main shock rupture asperity. This strengthens an earlier argument that the 1999 Quepos earthquake ruptured specific bathymetric highs on the downgoing plate. We believe that subduction of this highly disrupted seafloor has established a set of conditions which presently limit the seismogenic zone to be between 10 and 35 km below sea level

Geodetic and seismic constraints on some seismogenic zone processes in Costa Rica

New seismic and geodetic data from Costa Rica provide insight into seismogenic zone processes in Central America, where the Cocos and Caribbean plates converge. Seismic data are from combined land and ocean bottom deployments in the Nicoya peninsula in northern Costa Rica and near the Osa peninsula in southern Costa Rica. In Nicoya, inversion of GPS data suggests two locked patches centered at 14 ± 2 and 39 ± 6 km depth. Interplate microseismicity is concentrated in the more freely slipping intermediate zone, suggesting that small interseismic earthquakes may not accurately outline the updip limit of the seismogenic zone, the rupture zone for future large earthquakes, at least over the short (∼1 year) observation period. We also estimate northwest motion of a coastal “sliver block” at 8 ± 3 mm/yr, probably related to oblique convergence. In the Osa region to the south, convergence is orthogonal to the trench. Cocos-Caribbean relative motion is partitioned here, with ∼8 cm/yr on the Cocos-Panama block boundary (including a component of permanent shortening across the Fila Costeña fold and thrust belt) and ∼1 cm/yr on the Panama block–Caribbean boundary. The GPS data suggest that the Cocos plate–Panama block boundary is completely locked from ∼10–50 km depth. This large locked zone, as well as associated forearc and back-arc deformation, may be related to subduction of the shallow Cocos Ridge and/or younger lithosphere compared to Nicoya, with consequent higher coupling and compressive stress in the direction of plate convergence

The Great Sumatra-Andaman Earthquake of 26 December 2004

The two largest earthquakes of the past 40 years ruptured a 1600-kilometer-long portion of the fault boundary between the Indo-Australian and southeastern Eurasian plates on 26 December 2004 [seismic moment magnitude (M_w) = 9.1 to 9.3] and 28 March 2005 (M_w = 8.6). The first event generated a tsunami that caused more than 283,000 deaths. Fault slip of up to 15 meters occurred near Banda Aceh, Sumatra, but to the north, along the Nicobar and Andaman Islands, rapid slip was much smaller. Tsunami and geodetic observations indicate that additional slow slip occurred in the north over a time scale of 50 minutes or longer

Imaging crustal structure in South-Central Costa Rica with Receiver Functions

An array of broadband seismometers transecting the Talamanca Range in southern Costa Rica was operated from 2005 until 2007. In combination with data from a short‐period network near Quepos in central Costa Rica, this data is analyzed by the receiver function method to image the crustal structure in south‐central Costa Rica. Two strong positive signals are seen in the migrated images, interpreted as the Moho (at around 35 km depth) and an intra‐crustal discontinuity (15 km depth). A relatively flat crustal and Moho interface underneath the north‐east flank of the Talamanca Range can be followed for a lateral distance of about 50 km parallel to the trench, with only slight changes in the overall geometry. Closer to the coast, the topography of the discontinuities shows several features, most notably a deeper Moho underneath the Talamanca Mountain Range and volcanic arc. Under the highest part of the mountain ranges, the Moho reaches a depth of about 50 km, which indicates that the mountain ranges are approximately isostatically compensated. Local deviations from the crustal thickness expected for isostatic equilibrium occur under the active volcanic arc and in south Costa Rica. In the transition region between the active volcanic arc and the Talamanca Range, both the Moho and intracrustal discontinuity appear distorted, possibly related to the southern edge of the active volcanic zone and deformation within the southern part of the Central Costa Rica Deformed Belt. Near the volcanoes Irazu and Turrialba, a shallow converter occurs, correlating with a low‐velocity, low‐density body seen in tomography and gravimetry. Applying a grid search for the crustal interface depth and vp/vs ratio cannot constrain vp/vs values well, but points to generally low values (<1.7) in the upper crust. This is consistent with quartz‐rich rocks forming the mountain range

Hydrogeological system of erosional convergent margins and its influence on tectonics and interplate seismogenesis

[1] Fluid distribution in convergent margins is by most accounts closely related to tectonics. This association has been widely studied at accretionary prisms, but at half of the Earth's convergent margins, tectonic erosion grinds down overriding plates, and here fluid distribution and its relation to tectonics remain speculative. Here we present a new conceptual model for the hydrological system of erosional convergent margins. The model is based largely on new data and recently published observations from along the Middle America Trench offshore Nicaragua and Costa Rica, and it is consistent with observations from other erosional margins. The observations indicate that erosional margins possess previously unrecognized distinct hydrogeological systems: Most fluid contained in the sediment pores and liberated by early dehydration reactions drains from the plate boundary through a fractured upper plate to seep at the seafloor across the slope, rather than migrating along the décollement toward the deformation front as described for accretionary prisms. The observations indicate that the relative fluid abundance across the plate-boundary fault zone and fluid migration influence long-term tectonics and the transition from aseismic to seismogenic behavior. The segment of the plate boundary where fluid appears to be more abundant corresponds to the locus of long-term tectonic erosion, where tectonic thinning of the overriding plate causes subsidence and the formation of the continental slope. This correspondence between observations indicates that tectonic erosion is possibly linked to the migration of overpressured fluids into the overriding plate. The presence of overpressured fluids at the plate boundary is compatible with the highest flow rates estimated at slope seeps. The change from aseismic to seismogenic behavior along the plate boundary of the erosional margin begins where the amount of fluid at the fault declines with depth, indicating a control on interplate earthquakes. A previously described similar observation along accreting plate boundaries strongly indicates that fluid abundance exerts a first-order control on interplate seismogenesis at all types of subduction zones. We hypothesize that fluid depletion with depth increases grain-to-grain contact, increasing effective stress on the fault, and modifies fault zone architecture from a thick fault zone to a narrower zone of localized slip

Seismogenic zone structure beneath the Nicoya Peninsula, Costa Rica, from three-dimensional local earthquake P- and S-wave tomography

The subduction plate interface along the Nicoya Peninsula, Costa Rica, generates damaging large (Mw > 7.5) earthquakes. We present hypocenters and 3-D seismic velocity models (VP and VP/VS) calculated using simultaneous inversion of P- and S-wave arrival time data recorded from small magnitude, local earthquakes to elucidate seismogenic zone structure. In this region, interseismic cycle microseismicity does not uniquely define the potential rupture extent of large earthquakes. Plate interface microseismicity extends from 12 to 26 and from 17 to 28 km below sea level beneath the southern and northern Nicoya Peninsula, respectively. Microseismicity offset across the plate suture of East Pacific Rise-derived and Cocos-Nazca Spreading Center-derived oceanic lithosphere is ∼5 km, revising earlier estimates suggesting ∼10 km of offset. Interplate seismicity begins downdip of increased locking along the plate interface imaged using GPS and a region of low VP along the plate interface. The downdip edge of plate interface microseismicity occurs updip of the oceanic slab and continental Moho intersection, possibly due to the onset of ductile behaviour. Slow forearc mantle wedge P-wave velocities suggest 20–30 per cent serpentinization across the Nicoya Peninsula region while calculated VP/VS values suggest 0–10 per cent serpentinization. Interpretation of VP/VS resolution at depth is complicated however due to ray path distribution. We posit that the forearc mantle wedge is regionally serpentinized but may still be able to sustain rupture during the largest seismogenic zone earthquakes