Subduction zone megathrusts produce the majority of the world's largest earthquakes. To understand the processes that control seismicity here, it is important to improve our knowledge on the subduction interface characteristics and its spatial variations. Nicoya Peninsula, Costa Rica, extends the continental landmass ~50 km towards the trench, making it a very suitable place to study interface activity from right on the top of the seismogenic zone of the Middle America Subduction Zone (MASZ). We contribute to and utilize an earthquake catalog of 8765 analyst-picked events to determine the spatial variability in the earthquake frequency-magnitude distribution (FMD) in this region. After initial detection, magnitude determination and location, the events are precisely relocated using a locally derived 3-D seismic compressional and shear wave velocity model (DeShon et al., 2006). After restricting the dataset to events nearest the interface and with low formal error (horizontal location error < 5 km), we retain a subset of 3226 events that best resolves interface activity. Beneath Nicoya, we determine the spatial variability and mean FMD of the interface, and focus on the relative relationship of small-to-large earthquakes, termed b-value. Across the region, the overall b-value (1.18 ± 0.04) is higher than the global average (b~1), and much larger than the global subduction zone average (b~0.6). Significant variation in b-value is observed along the active plate interface. A well resolved zone of lower b is observed at and offshore central Nicoya coast, in a previously determined locked patch using deformation observed from Global Positioning System (GPS). Conversely, high b-values prevail over the subducted portion of the Fisher ridge, which likely ruptured in the 1990 Gulf of Nicoya Mw 7.0 earthquake. Observed regions of low b-value approximately corresponds to more strongly-locked segments of the subduction interface resulting in higher differential stress, which may be released in the next large interface earthquake in this part of the MASZ. Across the region the b-value is found to vary inversely with the degree of interface locking. Thus, it is proposed that if sufficient data exist, spatial b-value mapping can be used as a proxy to determine interface locking. This method is especially useful along the subduction megathrust, which is generally offshore making geodetic measurements difficult.

<p>The potential of variations of <i>b-</i>values in the G-R relation, log<i>N=a-bM</i> as earthquake precursors for small events (rockbursts) in Zinkgruvan mine, Sweden and for tectonic (large) earthquakes in the Andaman-Sumatra region were investigated. </p><p>The temporal frequency-magnitude distribution, <i>b(t)</i>, of rockbursts in Zinkgruvan mine was examined using high quality data recorded during the period November 1996 to April 2004 with magnitude ranges from <i>M</i><i>w</i>= -2.4 to 2.6. A sliding time-window was applied to compute <i>b</i>-values. The windows contain 50 events and were shifted with steps of 5 events. The results indicated that <i>b</i>-values significantly drop preceding rockbursts of magnitude <i>M</i><i>w</i>≥1.6.</p><p>Temporal and spatial variations of <i>b</i>-values were also examined for tectonic earthquakes, magnitude <i>M</i><i>w</i>≥4.1, in the Andaman-Sumatra region. Earthquake data from the ISC, IDC, NEIC and HVRD earthquake catalogs for a period from 01/01/1995 to 12/26/2004 were used for analysis. Spatial variations of <i>b</i> were calculated from circular areas containing 50 events, with nodes on a 0.5° x 0.5° grid. The analysis shows that <i>b(t)</i> estimates using data from different catalogs<i> </i>are comparable and that large earthquakes are preceded by a drop in <i>b(t)</i> of about 0.3~1. The distribution of stress deduced from <i>b</i>-value mapping shows that large earthquakes occurred in the high stress, i.e. low <i>b</i>-value, areas.</p><p>Aftershock sequences of the <i>M</i><i>w</i>=9, December 26, 2004 and the <i>M</i><i>w</i>=8.7, March 28, 2005 shocks were investigated by using the same methods. Results from aftershock sequences show similar behaviour as for the large and presumed independent main events.</p><p>The observed variations of <i>b</i>-values with time and in space support the hypothesis that <i>b</i>-values have a precursory potential. The method can be used for a wide range of earthquake magnitude, from microearthquakes (<i>M</i><i>w</i><3) to giant tectonic shocks (<i>M</i><i>w</i>~9) and for both of independent shocks and aftershocks.</p>

The potential of variations of b-values in the G-R relation, logN=a-bM as earthquake precursors for small events (rockbursts) in Zinkgruvan mine, Sweden and for tectonic (large) earthquakes in the Andaman-Sumatra region were investigated. The temporal frequency-magnitude distribution, b(t), of rockbursts in Zinkgruvan mine was examined using high quality data recorded during the period November 1996 to April 2004 with magnitude ranges from Mw= -2.4 to 2.6. A sliding time-window was applied to compute b-values. The windows contain 50 events and were shifted with steps of 5 events. The results indicated that b-values significantly drop preceding rockbursts of magnitude Mw≥1.6. Temporal and spatial variations of b-values were also examined for tectonic earthquakes, magnitude Mw≥4.1, in the Andaman-Sumatra region. Earthquake data from the ISC, IDC, NEIC and HVRD earthquake catalogs for a period from 01/01/1995 to 12/26/2004 were used for analysis. Spatial variations of b were calculated from circular areas containing 50 events, with nodes on a 0.5° x 0.5° grid. The analysis shows that b(t) estimates using data from different catalogs are comparable and that large earthquakes are preceded by a drop in b(t) of about 0.3~1. The distribution of stress deduced from b-value mapping shows that large earthquakes occurred in the high stress, i.e. low b-value, areas. Aftershock sequences of the Mw=9, December 26, 2004 and the Mw=8.7, March 28, 2005 shocks were investigated by using the same methods. Results from aftershock sequences show similar behaviour as for the large and presumed independent main events. The observed variations of b-values with time and in space support the hypothesis that b-values have a precursory potential. The method can be used for a wide range of earthquake magnitude, from microearthquakes (Mw&lt;3) to giant tectonic shocks (Mw~9) and for both of independent shocks and aftershocks.