At periods greater than 1000 seconds, Earth's seismic free oscillations have anomalously large amplitude when referenced to the Harvard Centroid Moment Tensor fault mechanism, which is estimated from 300- to 500-second surface waves. By using more realistic rupture models on a steeper fault derived from seismic body and surface waves, we approximated free oscillation amplitudes with a seismic moment (6.5 × 10 ²² Newton·meters) that corresponds to a moment magnitude of 9.15. With a rupture duration of 600 seconds, the fault-rupture models represent seismic observations adequately but underpredict geodetic displacements that argue for slow fault motion beneath the Nicobar and Andaman islands.

On February 6, 2023, two large earthquakes occurred near the Turkish town of Kahramanmaraş. The moment magnitude (Mw) 7.8 mainshock ruptured a 310 km-long segment of the left-lateral East Anatolian Fault, propagating through multiple releasing step-overs. The Mw 7.6 aftershock involved nearby left-lateral strike-slip faults of the East Anatolian Fault Zone, causing a 150 km-long rupture. We use remote-sensing observations to constrain the spatial distribution of coseismic slip for these two events and the February 20 Mw 6.4 aftershock near Antakya. Pixel tracking of optical and synthetic aperture radar data of the Sentinel-2 and Sentinel-1 satellites, respectively, provide near-field surface displacements. High-rate Global Navigation Satellite System data constrain each event separately. Coseismic slip extends from the surface to about 15 km depth with a shallow slip deficit. Most aftershocks cluster at major fault bends, surround the regions of high coseismic slip, or extend outward of the ruptured faults. For the mainshock, rupture propagation stopped southward at the diffuse termination of the East Anatolian fault and tapered off northward into the Pütürge segment, some 20 km south of the 2020 Mw 6.8 Elaziğ earthquake, highlighting a potential seismic gap. These events underscore the high seismic potential of immature fault systems.

We have developed the open-source software Tesseroids, a set of command-line programs to perform forward modeling of gravitational fields in spherical coordinates. The software is implemented in the C programming language and uses tesseroids (spherical prisms) for the discretization of the subsurface mass distribution. The gravitational fields of tesseroids are calculated numerically using the Gauss-Legendre quadrature (GLQ). We have improved upon an adaptive discretization algorithm to guarantee the accuracy of the GLQ integration. Our implementation of adaptive discretization uses a “stack-based” algorithm instead of recursion to achieve more control over execution errors and corner cases. The algorithm is controlled by a scalar value called the distance-size ratio ($D$) that determines the accuracy of the integration as well as the computation time. We have determined optimal values of $D$ for the gravitational potential, gravitational acceleration, and gravity gradient tensor by comparing the computed tesseroids effects with those of a homogenous spherical shell. The values required for a maximum relative error of 0.1% of the shell effects are $D=1$ for the gravitational potential, $D=1.5$ for the gravitational acceleration, and $D=8$ for the gravity gradients. Contrary to previous assumptions, our results show that the potential and its first and second derivatives require different values of $D$ to achieve the same accuracy. These values were incorporated as defaults in the software.