@article{OeztuerkBozzolanHolcombeetal.2022,
  author    = {{\"O}zt{\"u}rk, Ugur and Bozzolan, Elisa and Holcombe, Elizabeth A. and Shukla, Roopam and Pianosi, Francesca and Wagener, Thorsten},
  title     = {How climate change and unplanned urban sprawl bring more landslides},
  series = {Nature : the international weekly journal of science},
  volume    = {608},
  journal   = {Nature : the international weekly journal of science},
  number    = {7922},
  publisher = {Nature portfolio},
  address   = {Berlin},
  issn      = {0028-0836},
  doi       = {10.1038/d41586-022-02141-9},
  pages     = {262 -- 265},
  year      = {2022},
  abstract  = {More settlements will suffer as heavy rains and unregulated construction destabilize slopes in the tropics, models show.},
  language  = {en}
}
@article{SongJieGaoetal.2022,
  author    = {Song, Lina and Jie, Dongmei and Gao, Guizai and Liu, Lidan and Liu, Hongyan and Li, Dehui and Liu, Ying},
  title     = {Application of a topsoil phytolith dataset to quantitative paleoclimate reconstruction in Northeast China},
  series = {Palaeogeography, palaeoclimatology, palaeoecology : an international journal for the geo-sciences},
  volume    = {601},
  journal   = {Palaeogeography, palaeoclimatology, palaeoecology : an international journal for the geo-sciences},
  publisher = {Elsevier},
  address   = {Amsterdam},
  issn      = {0031-0182},
  doi       = {10.1016/j.palaeo.2022.111108},
  pages     = {12},
  year      = {2022},
  abstract  = {Although phytoliths are recognized as an important proxy for paleoenvironmental reconstruction, the quantitative relationship between phytoliths and climate is still debated. In order to provide an improved basis for phytolith-based paleoclimate reconstructions, a representative modern phytolith dataset is essential. Here, we synthesize a modern topsoil phytolith dataset for Northeast China, analyze its climatic significance, and apply it to a fossil phytolith series from the Hani peat core in Northeast China. The dataset comprises 660 topsoil phytolith assemblages from 289 sample sites. We compiled modern meteorological data to assess the quantitative relationship between the phytolith assemblages and climatic variables. Detrended correspondence analysis (DCA) and Redundancy analysis (RDA) were used to determine the dominant climatic variable influencing the phytolith distributions. The results showed that mean annual temperature (MAT) is the dominant variable controlling the spatial distribution of phytoliths, accounting for 8.91\% of the total variance. Transfer function based on inverse deshrinking locally-weighted weighted averaging (LWWA_Inv) was developed for MAT (R-_boot(2) = 0.86, RMSEP = 1.02 degrees C). Applying the LWWA_Inv transfer function to fossil phytolith records from the Hani peat core enables quantitative inferences to be made about Holocene climate changes in Northeast China. Overall, combined with the LWWA_Inv method, the topsoil phytolith dataset of Northeast China can be used for reliable quantitative MAT reconstruction.},
  language  = {en}
}
@article{SteirouGerlitzSunetal.2022,
  author    = {Steirou, Eva and Gerlitz, Lars and Sun, Xun and Apel, Heiko and Agarwal, Ankit and Totz, Sonja Juliana and Merz, Bruno},
  title     = {Towards seasonal forecasting of flood probabilities in Europe using climate and catchment information},
  series = {Scientific reports},
  volume    = {12},
  journal   = {Scientific reports},
  number    = {1},
  publisher = {Nature portfolio},
  address   = {Berlin},
  issn      = {2045-2322},
  doi       = {10.1038/s41598-022-16633-1},
  pages     = {10},
  year      = {2022},
  abstract  = {We investigate whether the distribution of maximum seasonal streamflow is significantly affected by catchment or climate state of the season/month ahead. We fit the Generalized Extreme Value (GEV) distribution to extreme seasonal streamflow for around 600 stations across Europe by conditioning the GEV location and scale parameters on 14 indices, which represent the season-ahead climate or catchment state. The comparison of these climate-informed models with the classical GEV distribution, with time-constant parameters, suggests that there is a substantial potential for seasonal forecasting of flood probabilities. The potential varies between seasons and regions. Overall, the season-ahead catchment wetness shows the highest potential, although climate indices based on large-scale atmospheric circulation, sea surface temperature or sea ice concentration also show some skill for certain regions and seasons. Spatially coherent patterns and a substantial fraction of climate-informed models are promising signs towards early alerts to increase flood preparedness already a season ahead.},
  language  = {en}
}
@article{AstudilloSotomayorJaraMunozMelnicketal.2021,
  author    = {Astudillo-Sotomayor, Luis and Jara Mu{\~n}oz, Julius and Melnick, Daniel and Cort{\´e}s-Aranda, Joaqu{\´i}n and Tassara, Andr{\´e}s and Strecker, Manfred},
  title     = {Fast Holocene slip and localized strain along the Liqui{\~n}e-Ofqui strike-slip fault system, Chile},
  series = {Scientific reports},
  volume    = {11},
  journal   = {Scientific reports},
  number    = {1},
  publisher = {Macmillan Publishers Limited, part of Springer Nature},
  address   = {London},
  issn      = {2045-2322},
  doi       = {10.1038/s41598-021-85036-5},
  pages     = {10},
  year      = {2021},
  abstract  = {In active tectonic settings dominated by strike-slip kinematics, slip partitioning across subparallel faults is a common feature; therefore, assessing the degree of partitioning and strain localization is paramount for seismic hazard assessments. Here, we estimate a slip rate of 18.8 +/- 2.0 mm/year over the past 9.0 +/- 0.1 ka for a single strand of the Liquirie-Ofqui Fault System, which straddles the Main Cordillera in Southern Chile. This Holocene rate accounts for similar to 82\% of the trench-parallel component of oblique plate convergence and is similar to million-year estimates integrated over the entire fault system. Our results imply that strain localizes on a single fault at millennial time scale but over longer time scales strain localization is not sustained. The fast millennial slip rate in the absence of historical Mw> 6.5 earthquakes along the Liquine-Ofqui Fault System implies either a component of aseismic slip or Mw similar to 7 earthquakes involving multi-trace ruptures and > 150-year repeat times. Our results have implications for the understanding of strike-slip fault system dynamics within volcanic arcs and seismic hazard assessments.},
  language  = {en}
}
@article{LopezCominoCescaNiemzetal.2021,
  author    = {L{\´o}pez-Comino, Jos{\´e} {\´A}ngel and Cesca, Simone and Niemz, Peter and Dahm, Torsten and Zang, Arno},
  title     = {Rupture directivity in 3D inferred from acoustic emissions events in a mine-scale hydraulic fracturing experiment},
  series = {Frontiers in Earth Science},
  volume    = {9},
  journal   = {Frontiers in Earth Science},
  publisher = {Frontiers Media},
  address   = {Lausanne},
  issn      = {2296-6463},
  doi       = {10.3389/feart.2021.670757},
  pages     = {9},
  year      = {2021},
  abstract  = {Rupture directivity, implying a predominant earthquake rupture propagation direction, is typically inferred upon the identification of 2D azimuthal patterns of seismic observations for weak to large earthquakes using surface-monitoring networks. However, the recent increase of 3D monitoring networks deployed in the shallow subsurface and underground laboratories toward the monitoring of microseismicity allows to extend the directivity analysis to 3D modeling, beyond the usual range of magnitudes. The high-quality full waveforms recorded for the largest, decimeter-scale acoustic emission (AE) events during a meter-scale hydraulic fracturing experiment in granites at similar to 410 m depth allow us to resolve the apparent durations observed at each AE sensor to analyze 3D-directivity effects. Unilateral and (asymmetric) bilateral ruptures are then characterized by the introduction of a parameter kappa, representing the angle between the directivity vector and the station vector. While the cloud of AE activity indicates the planes of the hydrofractures, the resolved directivity vectors show off-plane orientations, indicating that rupture planes of microfractures on a scale of centimeters have different geometries. Our results reveal a general alignment of the rupture directivity with the orientation of the minimum horizontal stress, implying that not only the slip direction but also the fracture growth produced by the fluid injections is controlled by the local stress conditions.},
  language  = {en}
}
@article{EsfahaniGholamiOhrnberger2020,
  author    = {Esfahani, Reza Dokht Dolatabadi and Gholami, Ali and Ohrnberger, Matthias},
  title     = {An inexact augmented Lagrangian method for nonlinear dispersion-curve inversion using Dix-type global linear approximation},
  series = {Geophysics : a journal of general and applied geophysics},
  volume    = {85},
  journal   = {Geophysics : a journal of general and applied geophysics},
  number    = {3},
  publisher = {GeoScienceWorld},
  address   = {Tulsa, Okla.},
  issn      = {0016-8033},
  doi       = {10.1190/geo2019-0717.1},
  pages     = {EN77 -- EN85},
  year      = {2020},
  abstract  = {Dispersion-curve inversion of Rayleigh waves to infer subsurface shear-wave velocity is a long-standing problem in seismology. Due to nonlinearity and ill-posedness, sophisticated regularization techniques are required to solve the problem for a stable velocity model. We have formulated the problem as a minimization problem with nonlinear operator constraint and then solve it by using an inexact augmented Lagrangian method, taking advantage of the Haney-Tsai Dix-type relation (a global linear approximation of the nonlinear forward operator). This replaces the original regularized nonlinear problem with iterative minimization of a more tractable regularized linear problem followed by a nonlinear update of the phase velocity (data) in which the update can be performed accurately with any forward modeling engine, for example, the finite-element method. The algorithm allows discretizing the medium with thin layers (for the finite-element method) and thus omitting the layer thicknesses from the unknowns and also allows incorporating arbitrary regularizations to shape the desired velocity model. In this research, we use total variation regularization to retrieve the shear-wave velocity model. We use two synthetic and two real data examples to illustrate the performance of the inversion algorithm with total variation regularization. We find that the method is fast and stable, and it converges to the solution of the original nonlinear problem.},
  language  = {en}
}
@article{SpallanzaniKogaCichyetal.2022,
  author    = {Spallanzani, Roberta and Koga, Kenneth T. and Cichy, Sarah B. and Wiedenbeck, Michael and Schmidt, Burkhard C. and Oelze, Marcus and Wilke, Max},
  title     = {Lithium and boron diffusivity and isotopic fractionation in hydrated rhyolitic melts},
  series = {Contributions to mineralogy and petrology},
  volume    = {177},
  journal   = {Contributions to mineralogy and petrology},
  number    = {8},
  publisher = {Springer},
  address   = {New York},
  issn      = {0010-7999},
  doi       = {10.1007/s00410-022-01937-2},
  pages     = {17},
  year      = {2022},
  abstract  = {Lithium and boron are trace components of magmas, released during exsolution of a gas phase during volcanic activity. In this study, we determine the diffusivity and isotopic fractionation of Li and B in hydrous silicate melts. Two glasses were synthesized with the same rhyolitic composition (4.2 wt\% water), having different Li and B contents; these were studied in diffusion-couple experiments that were performed using an internally heated pressure vessel, operated at 300 MPa in the temperature range 700-1250 degrees C for durations from 0 s to 24 h. From this we determined activation energies for Li and B diffusion of 57 +/- 4 kJ/mol and 152 +/- 15 kJ/mol with pre-exponential factors of 1.53 x 10(-7) m(2)/s and 3.80 x 10(-8) m(2)/s, respectively. Lithium isotopic fractionation during diffusion gave beta values between 0.15 and 0.20, whereas B showed no clear isotopic fractionation. Our Li diffusivities and isotopic fractionation results differ somewhat from earlier published values, but overall confirm that Li diffusivity increases with water content. Our results on B diffusion show that similarly to Li, B mobility increases in the presence of water. By applying the Eyring relation, we confirm that B diffusivity is limited by viscous flow in silicate melts. Our results on Li and B diffusion present a new tool for understanding degassing-related processes, offering a potential geospeedometer to measure volcanic ascent rates.},
  language  = {en}
}
@article{KothaWeatherillBindietal.2022,
  author    = {Kotha, Sreeram Reddy and Weatherill, Graeme and Bindi, Dino and Cotton, Fabrice},
  title     = {Near-source magnitude scaling of spectral accelerations},
  series = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering},
  volume    = {20},
  journal   = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering},
  number    = {3},
  publisher = {Springer},
  address   = {Dordrecht},
  issn      = {1570-761X},
  doi       = {10.1007/s10518-021-01308-5},
  pages     = {1343 -- 1370},
  year      = {2022},
  abstract  = {Ground-motion models (GMMs) are often used to predict the random distribution of Spectral accelerations (SAs) at a site due to a nearby earthquake. In probabilistic seismic hazard and risk assessment, large earthquakes occurring close to a site are considered as critical scenarios. GMMs are expected to predict realistic SAs with low within-model uncertainty (sigma(mu)) for such rare scenarios. However, the datasets used to regress GMMs are usually deficient of data from critical scenarios. The (Kotha et al., A Regionally Adaptable Ground-Motion Model for Shallow Crustal Earthquakes in Europe Bulletin of Earthquake Engineering 18:4091-4125, 2020) GMM developed from the Engineering strong motion (ESM) dataset was found to predict decreasing short-period SAs with increasing M-W >= M-h = 6.2, and with large sigma(mu) at near-source distances <= 30km. In this study, we updated the parametrisation of the GMM based on analyses of ESM and the Near source strong motion (NESS) datasets. With M-h = 5.7, we could rectify the M-W scaling issue, while also reducing sigma(mu). at M-W >= M-h. We then evaluated the GMM against NESS data, and found that the SAs from a few large, thrust-faulting events in California, New Zealand, Japan, and Mexico are significantly higher than GMM median predictions. However, recordings from these events were mostly made on soft-soil geology, and contain anisotropic pulse-like effects. A more thorough non-ergodic treatment of NESS was not possible because most sites sampled unique events in very diverse tectonic environments. We provide an updated set of GMM coefficients,sigma(mu), and heteroscedastic variance models; while also cautioning against its application for M-W <= 4 in low-moderate seismicity regions without evaluating the homogeneity of M-W estimates between pan-European ESM and regional datasets.},
  language  = {en}
}
@phdthesis{Koyan2024,
  author    = {Koyan, Philipp},
  title     = {3D attribute analysis and classification to interpret ground-penetrating radar (GPR) data collected across sedimentary environments: Synthetic studies and field examples},
  doi       = {10.25932/publishup-63948},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-639488},
  school      = {Universit{\"a}t Potsdam},
  pages     = {xi, 115, A51},
  year      = {2024},
  abstract  = {Die Untersuchung des oberfl{\"a}chennahen Untergrundes erfolgt heutzutage bei Frage- stellungen aus den Bereichen des Bauwesens, der Arch{\"a}ologie oder der Geologie und Hydrologie oft mittels zerst{\"o}rungsfreier beziehungsweise zerst{\"o}rungsarmer Methoden der angewandten Geophysik. Ein Bereich, der eine immer zentralere Rolle in Forschung und Ingenieurwesen einnimmt, ist die Untersuchung von sediment{\"a}ren Umgebungen, zum Beispiel zur Charakterisierung oberfl{\"a}chennaher Grundwassersysteme. Ein in diesem Kontext h{\"a}ufig eingesetztes Verfahren ist das des Georadars (oftmals GPR - aus dem Englischen ground-penetrating radar). Dabei werden kurze elektromagnetische Impulse von einer Antenne in den Untergrund ausgesendet, welche dort wiederum an Kontrasten der elektromagnetischen Eigenschaften (wie zum Beispiel an der Grundwasseroberfl{\"a}che) reflektiert, gebrochen oder gestreut werden. Eine Empfangsantenne zeichnet diese Signale in Form derer Amplituden und Laufzeiten auf. Eine Analyse dieser aufgezeichneten Signale erm{\"o}glicht Aussagen {\"u}ber den Untergrund, beispielsweise {\"u}ber die Tiefenlage der Grundwasseroberfl{\"a}che oder die Lagerung und Charakteristika oberfl{\"a}chennaher Sedimentschichten. Dank des hohen Aufl{\"o}sungsverm{\"o}gens der GPR-Methode sowie stetiger technologischer Entwicklungen erfolgt heutzutage die Aufzeichnung von GPR- Daten immer h{\"a}ufiger in 3D. Trotz des hohen zeitlichen und technischen Aufwandes f{\"u}r die Datenakquisition und -bearbeitung werden die resultierenden 3D-Datens{\"a}tze, welche den Untergrund hochaufl{\"o}send abbilden, typischerweise von Hand interpretiert. Dies ist in der Regel ein {\"a}ußerst zeitaufwendiger Analyseschritt. Daher werden oft repr{\"a}sentative 2D-Schnitte aus dem 3D-Datensatz gew{\"a}hlt, in denen markante Reflektionsstrukuren markiert werden. Aus diesen Strukturen werden dann sich {\"a}hnelnde Bereiche im Untergrund als so genannte Radar-Fazies zusammengefasst. Die anhand von 2D-Schnitten erlangten Resultate werden dann als repr{\"a}sentativ f{\"u}r die gesamte untersuchte Fl{\"a}che angesehen. In dieser Form durchgef{\"u}hrte Interpretationen sind folglich oft unvollst{\"a}ndig sowie zudem in hohem Maße von der Expertise der Interpretierenden abh{\"a}ngig und daher in der Regel nicht reproduzierbar. Eine vielversprechende Alternative beziehungsweise Erg{\"a}nzung zur manuellen In- terpretation ist die Verwendung von so genannten GPR-Attributen. Dabei werden nicht die aufgezeichneten Daten selbst, sondern daraus abgeleitete Gr{\"o}ßen, welche die markanten Reflexionsstrukturen in 3D charakterisieren, zur Interpretation herangezogen. In dieser Arbeit wird anhand verschiedener Feld- und Modelldatens{\"a}tze untersucht, welche Attribute sich daf{\"u}r insbesondere eignen. Zudem zeigt diese Arbeit, wie ausgew{\"a}hlte Attribute mittels spezieller Bearbeitungs- und Klassifizierungsmethoden zur Erstellung von 3D-Faziesmodellen genutzt werden k{\"o}nnen. Dank der M{\"o}glichkeit der Erstellung so genannter attributbasierter 3D-GPR-Faziesmodelle k{\"o}nnen zuk{\"u}nftige Interpretationen zu gewissen Teilen automatisiert und somit effizienter durchgef{\"u}hrt werden. Weiterhin beschreiben die so erhaltenen Resultate den untersuchten Untergrund in reproduzierbarer Art und Weise sowie umf{\"a}nglicher als es bisher mittels manueller Interpretationsmethoden typischerweise m{\"o}glich war.},
  language  = {en}
}
@article{MerzBassoFischeretal.2022,
  author    = {Merz, Bruno and Basso, Stefano and Fischer, Svenja and Lun, David and Bloeschl, Guenter and Merz, Ralf and Guse, Bjorn and Viglione, Alberto and Vorogushyn, Sergiy and Macdonald, Elena and Wietzke, Luzie and Schumann, Andreas},
  title     = {Understanding heavy tails of flood peak distributions},
  series = {Water resources research},
  volume    = {58},
  journal   = {Water resources research},
  number    = {6},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {0043-1397},
  doi       = {10.1029/2021WR030506},
  pages     = {37},
  year      = {2022},
  abstract  = {Statistical distributions of flood peak discharge often show heavy tail behavior, that is, extreme floods are more likely to occur than would be predicted by commonly used distributions that have exponential asymptotic behavior. This heavy tail behavior may surprise flood managers and citizens, as human intuition tends to expect light tail behavior, and the heaviness of the tails is very difficult to predict, which may lead to unnecessarily high flood damage. Despite its high importance, the literature on the heavy tail behavior of flood distributions is rather fragmented. In this review, we provide a coherent overview of the processes causing heavy flood tails and the implications for science and practice. Specifically, we propose nine hypotheses on the mechanisms causing heavy tails in flood peak distributions related to processes in the atmosphere, the catchment, and the river system. We then discuss to which extent the current knowledge supports or contradicts these hypotheses. We also discuss the statistical conditions for the emergence of heavy tail behavior based on derived distribution theory and relate them to the hypotheses and flood generation mechanisms. We review the degree to which the heaviness of the tails can be predicted from process knowledge and data. Finally, we recommend further research toward testing the hypotheses and improving the prediction of heavy tails.},
  language  = {en}
}