We deployed a dense geodetic and seismological network in the Atacama seismic gap in Chile. We derive a microseismicity catalog of >30,000 events, time series from 70 GNSS stations, and apply a transdimensional Bayesian inversion to estimate interplate locking degree. We identify two highly locked regions of different sizes whose geometries appear to control seismicity patterns. Interface seismicity concentrates beneath the coastline just downdip of the highest locking. A region of lower interplate locking around 27.5ºS coincides with higher seismicity levels, a high number of repeating earthquakes and events extending further towards the trench. Having shown numerous signs of aseismic deformation (slow-slip events and earthquake swarms), this area is situated where the Copiapó Ridge is subducted. While these findings suggest that the structure of the downgoing oceanic plate prescribes patterns of interplate locking and seismicity, we note that the Taltal Ridge further north lacks a similar signature.

The spatial variability of the regional crustal stress in northern Chile is resolved. We infer a margin-parallel compressive crustal stress regime along the coastal region, similar to crustal stress observations in Cascadia and Japan. The Andean Precordillera shows a distinct stress field associated with a strike-slip faulting regime. These results are constrained by over a decade of observations, for which earthquake catalogs report thousands of events in the continental crust. We present focal mechanisms for 817 of these crustal earthquakes, including mechanism qualities. The best mechanisms were grouped and inverted to infer the stress-field variability. We interpret the margin-parallel compression to be caused by the concave shape of the margin and the locking of the plate interface. The inferred strike-slip regime in the Andes agrees with previous studies and has been proposed to be mostly caused by local stresses imposed by a thicker crust.

Double seismic zones (DSZs), parallel planes of intermediate-depth earthquakes inside oceanic slabs, have been observed in a number of subduction zones and may be a ubiquitous feature of downgoing oceanic plates. Focal mechanism observations from DSZ earthquakes sample the intraslab stress field at two distinct depth levels within the downgoing lithosphere. A pattern of downdip compressive over downdip extensive events was early on interpreted to indicate an unbending-dominated intraslab stress field. In the present study, we show that the intraslab stress field in the depth range of DSZs is much more variable than previously thought. Compiling DSZ locations and mechanisms from literature, we observe that the “classical’ pattern of compressive over extensive events is only observed at about half of the DSZ locations around the globe. The occurrence of extensional mechanisms across both planes accounts for most other regions. To obtain an independent estimate of the bending state of slabs at intermediate depths, we compute (un)bending estimates from slab geometries taken from the slab2 compilation of slab surface depths. We find no clear global prevalence of slab unbending at intermediate depths, and the occurrence of DSZ seismicity does not appear to be limited to regions of slab unbending. Focal mechanism observations are frequently inconsistent with (un)bending estimates from slab geometries, which may imply that bending stresses are not always prevalent, and that other stress types such as in-plane tension due to slab pull or shallow compression due to friction along the plate interface may also play an important role.