[1] We present the first regional surface velocity field for Central America, showing crustal response to interaction of the Cocos and Caribbean plates. Elastic half-space models for interseismic strain accumulation on the dipping subduction plate boundary fit the GPS data well and show strain accumulation offshore and beneath the Nicoya and Osa peninsulas in Costa Rica but not in Nicaragua. Since large subduction zone earthquakes occur in Nicaragua, we suggest that interseismic locking in Nicaragua and some other parts of Central America occurs but is mainly shallow, <20 km depth, too far offshore to be detected by our on-land GPS measurements. Our data also show significant trench-parallel motion for most of the region, generally interpreted as due to oblique convergence and strong mechanical coupling between subducting and overriding plates. However, trench-parallel motion is also observed in central Costa Rica, where plate convergence is normal to the trench, and in the Nicaraguan fore arc, where trench-parallel motion is fast, up to 9 mm a À1 , but mechanical coupling is low. A finite element model of collision (as opposed to subduction) involving the aseismic Cocos Ridge also fits the GPS surface velocity field, most significantly reproducing the pattern of trench-parallel motion. We infer that buoyant, thickened CNS-2-Cocos Ridge crust resists normal subduction and instead acts as an indenter to the Caribbean plate, driving crustal shortening in southern Costa Rica and contributing to trench-parallel fore-arc motion in Costa Rica and perhaps Nicaragua as a type of tectonic escape.

No abstract

A vigorously degassing lava lake appeared inside the Santiago pit crater of Masaya volcano (Nicaragua) in December 2015, after years of degassing with no (or minor) incandescence. Here we present an unprecedented‐long (3 years) and continuous volcanic gas record that instrumentally characterizes the (re)activation of the lava lake. Our results show that, before appearance of the lake, the volcanic gas plume composition became unusually CO2 rich, as testified by high CO2/SO2 ratios (mean: 12.2 ± 6.3) and low H2O/CO2 ratios (mean: 2.3 ± 1.3). The volcanic CO2 flux also peaked in November 2015 (mean: 81.3 ± 40.6 kg/s; maximum: 247 kg/s). Using results of magma degassing models and budgets, we interpret this elevated CO2 degassing as sourced by degassing of a volatile‐rich fast‐overturning (3.6–5.2 m3 s−1) magma, supplying CO2‐rich gas bubbles from minimum equivalent depths of 0.36–1.4 km. We propose this elevated gas bubble supply destabilized the shallow (<1 km) Masaya magma reservoir, leading to upward migration of vesicular (buoyant) resident magma, and ultimately to (re)formation of the lava lake. At onset of lava lake activity on 11 December 2015 (constrained by satellite‐based MODIS thermal observations), the gas emissions transitioned to more SO2‐rich composition, and the SO2 flux increased by a factor ∼40% (11.4 ± 5.2 kg/s) relative to background degassing (8.0 kg/s), confirming faster than normal (4.4 versus ∼3 m3 s−1) shallow magma convection. Based on thermal energy records, we estimate that only ∼0.8 of the 4.4 m3 s−1 of magma actually reached the surface to manifest into a convecting lava lake, suggesting inefficient transport of magma in the near‐surface plumbing system.

This work presents a new database of SO2 and CO2 fluxes from the Southern Central American Volcanic Arc (SCAVA) for the period 2015–2016. We report ∼300 SO2 flux measurements from 10 volcanoes and gas ratios from 11 volcanoes in Costa Rica and Nicaragua representing the most extensive available assessment of this ∼500 km arc segment. The SO2 flux from SCAVA is estimated at 6,240 ± 1,150 T/d, about a factor of three higher than previous estimations (1972–2013). We attribute this increase in part to our more complete assessment of the arc. Another consideration in interpreting the difference is the context of increased volcanic activity, as there were more eruptions in 2015–2016 than in any period since ∼1980. A potential explanation for increased degassing and volcanic activity is a change in crustal stress regime (from compression to extension, opening volcanic conduits) following two large (Mw > 7) earthquakes in the region in 2012. The CO2 flux from the arc is estimated at 22,500 ± 4,900 T/d, which is equal to or greater than estimates of C input into the SCAVA subduction zone. Time‐series data sets for arc degassing need to be improved in temporal and spatial coverage to robustly constrain volatile budgets and tectonic controls. Arc volatile budgets are strongly influenced by short‐lived degassing events and arc systems likely display significant short‐term variations in volatile output, calling for expansion of nascent geochemical monitoring networks to achieve spatial and temporal coverage similar to traditional geophysical networks.

[1] Campaign GPS data from a network in the Nicaraguan forearc show a strong component of arc-parallel motion indicating northwest translation of a nearly rigid forearc sliver. Our measured mean velocity for forearc sites of 15.1 mm yr À1 agrees well with the arc-parallel sliver motion predicted previously by DeMets (2001) derived from closure constraints on oblique convergence between the Cocos and Caribbean plates. The lack of a northeasterly oriented arc-normal component of motion in forearc velocities indicates that there are complexities involved beyond a simple interpretation of sliver motion being driven by oblique convergence. The forearc is reasonably well-fit by rigid rotation about an Euler pole with a rms misfit of residual velocities of 4.9 ± 2.6 mm yr À1. Current motion of the forearc sliver relative to the stable Caribbean plate yields predominantly boundary parallel NW motion with boundary normal extension in the northwestern region averaging $5 mm yr