More than 90% of the energy trapped on Earth by increasingly abundant greenhouse gases is absorbed by the ocean. Monitoring the resulting ocean warming remains a challenging sampling problem. To complement existing point measurements, we introduce a method that infers basin-scale deep-ocean temperature changes from the travel times of sound waves that are generated by repeating earthquakes. A first implementation of this seismic ocean thermometry constrains temperature anomalies averaged across a 3000-kilometer-long section in the equatorial East Indian Ocean with a standard error of 0.0060 kelvin. Between 2005 and 2016, we find temperature fluctuations on time scales of 12 months, 6 months, and ~10 days, and we infer a decadal warming trend that substantially exceeds previous estimates.

The ocean is warming in response to accumulating greenhouse gases in the atmosphere, and its heat capacity dominates the climate system's thermal inertia. While the warming has been most pronounced in the surface ocean, heat transfer to the deep ocean importantly slows the climate change experienced at the surface (e.g.,

Over the course of recent decades, N-body simulations have become a standard tool for quantifying the gravitational perturbations that ensue in planet-forming disks. Within the context of such simulations, massive noncentral bodies are routinely classified into "big" and "small" particles, where big objects interact with all other objects self-consistently, while small bodies interact with big bodies but not with each other. Importantly, this grouping translates to an approximation scheme where the orbital evolution of small bodies is dictated entirely by the dynamics of the big bodies, yielding considerable computational advantages with little added cost in terms of astrophysical accuracy. Here we point out, however, that this scheme can also yield spurious dynamical behavior where, even in the absence of big bodies within a simulation, indirect coupling among small bodies can lead to excitation of the constituent "non-interacting" orbits. We demonstrate this self-stirring by carrying out a sequence of numerical experiments, and confirm that this effect is largely independent of the time-step or the integration algorithm employed. Furthermore, adopting the growth of angular momentum deficit as a proxy for dynamical excitation, we explore its dependence on time, the cumulative mass of the system, as well as the total number of particles present in the simulation. Finally, we examine the degree of such indirect excitation within the context of conventional terrestrial planet formation calculations, and conclude that although some level of caution may be warranted, this effect plays a negligible role in driving the simulated dynamical evolution.

Many volcanoes emit a significant portion of the gas they transport to the atmosphere during continual passive degassing rather than during eruptions. To maintain a high gas and thermal flux without erupting magma, the flow field in the volcanic conduit must be approximately balanced with gas-rich, buoyant magma ascending and degassed, heavy magma descending. In vertical conduits, this exchange flow takes the form of core–annular flow, where the gas-rich magma forms a core enclosed by an annulus of degassed magma. The flow dynamics of core–annular flow have been studied extensively in fluid dynamics, but mostly for constant material properties. Our study aims to advance our understanding of how core–annular flow responds to volatile exsolution – a simple, yet ubiquitous disruption in volcanic conduits, which alters both the density and the viscosity of the core fluid. By deriving an evolution equation for the core–annular interface based on a generalized exchange-flow condition using a lubrication approximation, we find that the response of the system to volatile exsolution depends on the conduit flow regime. The same nucleation event can lead to a flow adjustment only in the upper, only in the lower or in both portions of the volcanic conduit. Our results emphasize that the thermodynamic evolution of magma properties and volcanic conduit flow are intricately linked, which may help understand the observed variability of eruptive behaviour at persistently degassing volcanoes.