The origin of the First Bend of the Yangtze River is key to understanding the birth of the modern Yangtze River. Despite considerable efforts, the timing and mechanism of formation of the First Bend remain highly debated. Inverse river-profile modeling of three tributaries (Chongjiang, Lima, and Gudu) of the Jinsha River, integrated with regional tectonic and geomorphic interpretations, allows the onset of incision at the First Bend to be constrained to 28-20 Ma. The spatio-temporal coincidence of initial river incision and activity of Yulong Thrust Belt in southeastern Tibet highlights thrusting to be fundamental in reshaping the pre-existing stream network at the First Bend. These results enable us to reinterpret a change in sedimentary environment from a braided river to a swamp-like lake in the Jianchuan Basin south of the First Bend, recording the destruction of the hypothesized southwards-flowing paleo-Jinsha and Shuiluo Rivers at ~36-35 Ma by magmatism. During the late Oligoceneearly Miocene, the paleo-Shuiluo River was diverted to the north by focused rock uplift due to thrusting along the Yulong Thrust Belt, which also led to exhumation of the Jianchuan Basin. Diversion of the paleo-Shuiluo River can be explained by capture from a downstream river in the footwall of the Yulong Thrust Belt. Subsequent rapid headward erosion, that was caused by thrusting-induced drop of local base level, is recorded by upstream younging ages for the onset of incision and led to the formation of the First Bend. The combination of new ages for the onset of incision at 28-20 Ma at the First Bend and younger ages upstream indicates northwards expansion of the Jinsha River at a rate of 62 +/- 18 mm/yr. Our results suggest that the origin of the First Bend was likely triggered by thrusting at 28-20 Ma, after which the Yangtze River formed.

1-D site response analysis dominates earthquake engineering practice, while local 2-D/3-D models are often required at sites where the site response is complex. For such sites, the 1-D representation of the soil column can account neither for topographic effects or dipping layers nor for locally generated horizontally propagating surface waves. It then remains a crucial task to identify whether the site response can be modelled sufficiently precisely by 1-D analysis. In this study we develop a method to classify sites according to their 1-D or 2-D/3-D nature. This classification scheme is based on the analysis of surface earthquake recordings and the evaluation of the variability and similarity of the horizontal Fourier spectra. The taxonomy is focused on capturing significant directional dependencies and interevent variabilities indicating a more probable 2-D/3-D structure around the site causing the ground motion to be more variable. While no significant correlation of the 1-D/3-D site index with environmental parameters and site proxies seems to exist, a reduction in the within-site (single-station) variability is found. The reduction is largest (up to 20 per cent) for purely 1-D sites. Although the taxonomy system is developed using surface stations of the KiK-net network in Japan as considerable additional information is available, it can also be applied to any (non-downhole array) site.

Volcanic tremor signals are usually observed before or during volcanic eruptions and must be monitored to evaluate the volcanic activity. A challenge in studying seismic signals of volcanic origin is the coexistence of transient signal swarms and long-lasting volcanic tremor signals. Separating transient events from volcanic tremors can, therefore, contrib-ute to improving upon our understanding of the underlying physical processes. Exploiting the idea of harmonic-percussive separation in musical signal processing, we develop a method to extract the harmonic volcanic tremor signals and to detect tran-sient events from seismic recordings. Based on the similarity properties of spectrogram frames in the time-frequency domain, we decompose the signal into two separate spec-trograms representing repeating (harmonic) and nonrepeating (transient) patterns, which correspond to volcanic tremor signals and earthquake signals, respectively. We reconstruct the harmonic tremor signal in the time domain from the complex spectrogram of the repeating pattern by only considering the phase components for the frequency range in which the tremor amplitude spectrum is significantly contribut-ing to the energy of the signal. The reconstructed signal is, therefore, clean tremor signal without transient events. Furthermore, we derive a characteristic function suitable for the detection of tran-sient events (e.g., earthquakes) by integrating amplitudes of the nonrepeating spectro-gram over frequency at each time frame. Considering transient events like earthquakes, 78% of the events are detected for signal-to-noise ratio = 0.1 in our semisynthetic tests. In addition, we compared the number of detected earthquakes using our method for one month of continuous data recorded during the Holuhraun 2014-2015 eruption in Iceland with the bulletin presented in Agustsdottir et al. (2019). Our single station event detection algorithm identified 84% of the bulletin events. Moreover, we detected a total of 12,619 events, which is more than twice the number of the bulletin events.

Coupled reactive transport simulations are extremely demanding in terms of required computational power, which hampers their application and leads to coarsened and oversimplified domains. The chemical sub-process represents the major bottleneck: its acceleration is an urgent challenge which gathers increasing interdisciplinary interest along with pressing requirements for subsurface utilization such as spent nuclear fuel storage, geothermal energy and CO2 storage. In this context we developed POET (POtsdam rEactive Transport), a research parallel reactive transport simulator integrating algorithmic improvements which decisively speed up coupled simulations. In particular, POET is designed with a master/worker architecture, which ensures computational efficiency in both multicore and cluster compute environments. POET does not rely on contiguous grid partitions for the parallelization of chemistry but forms work packages composed of grid cells distant from each other. Such scattering prevents particularly expensive geochemical simulations, usually concentrated in the vicinity of a reactive front, from generating load imbalance between the available CPUs (central processing units), as is often the case with classical partitions. Furthermore, POET leverages an original implementation of the distributed hash table (DHT) mechanism to cache the results of geochemical simulations for further reuse in subsequent time steps during the coupled simulation. The caching is hence particularly advantageous for initially chemically homogeneous simulations and for smooth reaction fronts. We tune the rounding employed in the DHT on a 2D benchmark to validate the caching approach, and we evaluate the performance gain of POET's master/worker architecture and the DHT speedup on a 3D benchmark comprising around 650 000 grid elements. The runtime for 200 coupling iterations, corresponding to 960 simulation days, reduced from about 24 h on 11 workers to 29 min on 719 workers. Activating the DHT reduces the runtime further to 2 h and 8 min respectively. Only with these kinds of reduced hardware requirements and computational costs is it possible to realistically perform the longterm complex reactive transport simulations, as well as perform the uncertainty analyses required by pressing societal challenges connected with subsurface utilization.

Purpose Root growth, respiration, water uptake as well as root exudation induce biogeochemical patterns in the rhizosphere that can change dynamically over time. Our aim is to develop a method that provides complementary information on 3D root system architecture and biogeochemical gradients around the roots needed for the quantitative description of rhizosphere processes. Methods We captured for the first time the root system architecture of maize plants grown in rectangular rhizotrons in 3D using neutron computed laminography (NCL). Simultaneously, we measured pH and oxygen concentration using fluorescent optodes and the 2D soil water distribution by means of neutron radiography. We co-registered the 3D laminography data with the 2D oxygen and pH maps to analyze the sensor signal as a function of the distance between the roots and the optode. Results The 3D root system architecture was successfully segmented from the laminographic data. We found that exudation of roots in up to 2 mm distance to the pH optode induced patterns of local acidification or alkalization. Over time, oxygen gradients in the rhizosphere emerged for roots up to a distance of 7.5 mm. Conclusion Neutron computed laminography allows for a three-dimensional investigation of root systems grown in laterally extended rhizotrons as the ones designed for 2D optode imaging studies. The 3D information on root position within the rhizotrons derived by NCL explained measured 2D oxygen and pH distribution. The presented new combination of 3D and 2D imaging methods facilitates systematical investigations of a wide range of dynamic processes in the rhizosphere.