@article{SarrazinKumarBasuetal.2022, author = {Sarrazin, Fanny J. and Kumar, Rohini and Basu, Nandita B. and Musolff, Andreas and Weber, Michael and Van Meter, Kimberly J. and Attinger, Sabine}, title = {Characterizing catchment-scale nitrogen legacies and constraining their uncertainties}, series = {Water resources research}, volume = {58}, journal = {Water resources research}, number = {4}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0043-1397}, doi = {10.1029/2021WR031587}, pages = {32}, year = {2022}, abstract = {Improving nitrogen (N) status in European water bodies is a pressing issue. N levels depend not only on current but also past N inputs to the landscape, that have accumulated through time in legacy stores (e.g., soil, groundwater). Catchment-scale N models, that are commonly used to investigate in-stream N levels, rarely examine the magnitude and dynamics of legacy components. This study aims to gain a better understanding of the long-term fate of the N inputs and its uncertainties, using a legacy-driven N model (ELEMeNT) in Germany's largest national river basin (Weser; 38,450 km(2)) over the period 1960-2015. We estimate the nine model parameters based on a progressive constraining strategy, to assess the value of different observational data sets. We demonstrate that beyond in-stream N loading, soil N content and in-stream N concentration allow to reduce the equifinality in model parameterizations. We find that more than 50\% of the N surplus denitrifies (1480-2210 kg ha(-1)) and the stream export amounts to around 18\% (410-640 kg ha(-1)), leaving behind as much as around 230-780 kg ha(-1) of N in the (soil) source zone and 10-105 kg ha(-1) in the subsurface. A sensitivity analysis reveals the importance of different factors affecting the residual uncertainties in simulated N legacies, namely hydrologic travel time, denitrification rates, a coefficient characterizing the protection of organic N in source zone and N surplus input. Our study calls for proper consideration of uncertainties in N legacy characterization, and discusses possible avenues to further reduce the equifinality in water quality modeling.}, language = {en} } @article{JoziNajafabadiHaberlandRybergetal.2021, author = {Jozi Najafabadi, Azam and Haberland, Christian and Ryberg, Trond and Verwater, Vincent F. and Breton, Eline le and Handy, Mark R. and Weber, Michael}, title = {Relocation of earthquakes in the southern and eastern Alps (Austria, Italy) recorded by the dense, temporary SWATH-D network using a Markov chain Monte Carlo inversion}, series = {Solid earth : SE ; an interaktive open access journal of the European Geosciences Union}, volume = {12}, journal = {Solid earth : SE ; an interaktive open access journal of the European Geosciences Union}, number = {5}, publisher = {Copernicus}, address = {G{\"o}ttingen}, organization = {AlpArray Working Grp, AlpArray SWATH-D Working Grp}, issn = {1869-9529}, doi = {10.5194/se-12-1087-2021}, pages = {1087 -- 1109}, year = {2021}, abstract = {In this study, we analyzed a large seismological dataset from temporary and permanent networks in the southern and eastern Alps to establish high-precision hypocenters and 1-D V-P and V-P/V-S models. The waveform data of a subset of local earthquakes with magnitudes in the range of 1-4.2 M-L were recorded by the dense, temporary SWATH-D network and selected stations of the AlpArray network between September 2017 and the end of 2018. The first arrival times of P and S waves of earthquakes are determined by a semi-automatic procedure. We applied a Markov chain Monte Carlo inversion method to simultaneously calculate robust hypocenters, a 1-D velocity model, and station corrections without prior assumptions, such as initial velocity models or earthquake locations. A further advantage of this method is the derivation of the model parameter uncertainties and noise levels of the data. The precision estimates of the localization procedure is checked by inverting a synthetic travel time dataset from a complex 3-D velocity model and by using the real stations and earthquakes geometry. The location accuracy is further investigated by a quarry blast test. The average uncertainties of the locations of the earthquakes are below 500m in their epicenter and similar to 1.7 km in depth. The earthquake distribution reveals seismicity in the upper crust (0-20 km), which is characterized by pronounced clusters along the Alpine frontal thrust, e.g., the Friuli-Venetia (FV) region, the Giudicarie-Lessini (GL) and Schio-Vicenza domains, the Austroalpine nappes, and the Inntal area. Some seismicity also occurs along the Periadriatic Fault. The general pattern of seismicity reflects head-on convergence of the Adriatic indenter with the Alpine orogenic crust. The seismicity in the FV and GL regions is deeper than the modeled frontal thrusts, which we interpret as indication for southward propagation of the southern Alpine deformation front (blind thrusts).}, language = {en} } @article{JoziNajafabadiHaberlandLeBretonetal.2022, author = {Jozi Najafabadi, Azam and Haberland, Christian and Le Breton, Eline and Handy, Mark R. and Verwater, Vincent F. and Heit, Benjamin and Weber, Michael}, title = {Constraints on crustal structure in the vicinity of the adriatic indenter (European Alps) from Vp and Vp/Vs local earthquake tomography}, series = {Journal of geophysical research : Solid earth}, volume = {127}, journal = {Journal of geophysical research : Solid earth}, number = {2}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9313}, doi = {10.1029/2021JB023160}, pages = {22}, year = {2022}, abstract = {In this study, 3-D models of P-wave velocity (Vp) and P-wave and S-wave ratio (Vp/Vs) of the crust and upper mantle in the Eastern and eastern Southern Alps (northern Italy and southern Austria) were calculated using local earthquake tomography (LET). The data set includes high-quality arrival times from well-constrained hypocenters observed by the dense, temporary seismic networks of the AlpArray AASN and SWATH-D. The resolution of the LET was checked by synthetic tests and analysis of the model resolution matrix. The small inter-station spacing (average of similar to 15 km within the SWATH-D network) allowed us to image crustal structure at unprecedented resolution across a key part of the Alps. The derived P velocity model revealed a highly heterogeneous crustal structure in the target area. One of the main findings is that the lower crust is thickened, forming a bulge at 30-50 km depth just south of and beneath the Periadriatic Fault and the Tauern Window. This indicates that the lower crust decoupled both from its mantle substratum as well as from its upper crust. The Moho, taken to be the iso-velocity contour of Vp = 7.25 km/s, agrees with the Moho depth from previous studies in the European and Adriatic forelands. It is shallower on the Adriatic side than on the European side. This is interpreted to indicate that the European Plate is subducted beneath the Adriatic Plate in the Eastern and eastern Southern Alps.}, language = {en} } @article{EiblRosskopfSciottoetal.2022, author = {Eibl, Eva P. S. and Rosskopf, Martina and Sciotto, Mariangela and Currenti, Gilda and Di Grazia, Giuseppe and Jousset, Philippe and Kr{\"u}ger, Frank and Weber, Michael}, title = {Performance of a rotational sensor to decipher volcano seismic signals on Etna, Italy}, series = {Journal of geophysical research : Solid earth}, volume = {127}, journal = {Journal of geophysical research : Solid earth}, number = {6}, publisher = {Wiley}, address = {Hoboken, NJ}, issn = {0148-0227}, doi = {10.1029/2021JB023617}, pages = {22}, year = {2022}, abstract = {Volcano-seismic signals such as long-period events and tremor are important indicators for volcanic activity and unrest. However, their wavefield is complex and characterization and location using traditional seismological instrumentation is often difficult. In 2019 we recorded the full seismic wavefield using a newly developed 3C rotational sensor co-located with a 3C traditional seismometer on Etna, Italy. We compare the performance of the rotational sensor, the seismometer and the Istituto Nazionale di Geofisica e Vulcanologia-Osservatorio Etneo (INGV-OE) seismic network with respect to the analysis of complex volcano-seismic signals. We create event catalogs for volcano-tectonic (VT) and long-period (LP) events combining a STA/LTA algorithm and cross-correlations. The event detection based on the rotational sensor is as reliable as the seismometer-based detection. The LP events are dominated by SH-type waves. Derived SH phase velocities range from 500 to 1,000 m/s for LP events and 300-400 m/s for volcanic tremor. SH-waves compose the tremor during weak volcanic activity and SH- and SV-waves during sustained strombolian activity. We derive back azimuths using (a) horizontal rotational components and (b) vertical rotation rate and transverse acceleration. The estimated back azimuths are consistent with the INGV-OE event location for (a) VT events with an epicentral distance larger than 3 km and some closer events, (b) LP events and tremor in the main crater area. Measuring the full wavefield we can reliably analyze the back azimuths, phase velocities and wavefield composition for VT, LP events and tremor in regions that are difficult to access such as volcanoes.}, language = {en} } @article{MallonnPoppenhaegerGranzeretal.2022, author = {Mallonn, Matthias and Poppenh{\"a}ger, Katja and Granzer, Thomas and Weber, Michael and Strassmeier, Klaus G.}, title = {Detection capability of ground-based meter-sized telescopes for shallow exoplanet transits}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {657}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {0004-6361}, doi = {10.1051/0004-6361/202140599}, pages = {10}, year = {2022}, abstract = {Meter-sized ground-based telescopes are frequently used today for the follow-up of extrasolar planet candidates. While the transit signal of a Jupiter-sized object can typically be detected to a high level of confidence with small telescope apertures as well, the shallow transit dips of planets with the size of Neptune and smaller are more challenging to reveal. We employ new observational data to illustrate the photometric follow-up capabilities of meter-sized telescopes for shallow exoplanet transits. We describe in detail the capability of distinguishing the photometric signal of an exoplanet transit from an underlying trend in the light curve. The transit depths of the six targets we observed, Kepler-94b, Kepler-63b, K2-100b, K2-138b, K2-138c, and K2-138e, range from 3.9 ppt down to 0.3 ppt. For five targets of this sample, we provide the first ground-based photometric follow-up. The timing of three targets is precisely known from previous observations, and the timing of the other three targets is uncertain and we aim to constrain it. We detect or rule out the transit features significantly in single observations for the targets that show transits of 1.3 ppt or deeper. The shallower transit depths of two targets of 0.6 and 0.8 ppt were detected tentatively in single light curves, and were detected significantly by repeated observations. Only for the target of the shallowest transit depth of 0.3 ppt were we unable to draw a significant conclusion despite combining five individual light curves. An injection-recovery test on our real data shows that we detect transits of 1.3 ppt depth significantly in single light curves if the transit is fully covered, including out-of-transit data toward both sides, in some cases down to 0.7 ppt depth. For Kepler-94b, Kepler-63b, and K2-100b, we were able to verify the ephemeris. In the case of K2-138c with a 0.6 ppt deep transit, we were able to refine it, and in the case of K2-138e, we ruled out the transit in the time interval of more than ±1.5 σ of its current literature ephemeris.}, language = {en} } @article{DahmStillerMechieetal.2020, author = {Dahm, Torsten and Stiller, Manfred and Mechie, James and Heimann, Sebastian and Hensch, Martin and Woith, Heiko and Schmidt, Bernd and Gabriel, Gerald and Weber, Michael}, title = {Seismological and geophysical signatures of the deep crustal magma systems of the cenozoic volcanic fields Beneath the Eifel, Germany}, series = {Geochemistry, geophysics, geosystems}, volume = {21}, journal = {Geochemistry, geophysics, geosystems}, number = {9}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1525-2027}, doi = {10.1029/2020GC009062}, pages = {21}, year = {2020}, abstract = {The Quaternary volcanic fields of the Eifel (Rhineland-Palatinate, Germany) had their last eruptions less than 13,000 years ago. Recently, deep low-frequency (DLF) earthquakes were detected beneath one of the volcanic fields showing evidence of ongoing magmatic activity in the lower crust and upper mantle. In this work, seismic wide- and steep-angle experiments from 1978/1979 and 1987/1988 are compiled, partially reprocessed and interpreted, together with other data to better determine the location, size, shape, and state of magmatic reservoirs in the Eifel region near the crust-mantle boundary. We discuss seismic evidence for a low-velocity gradient layer from 30-36 km depth, which has developed over a large region under all Quaternary volcanic fields of the Rhenish Massif and can be explained by the presence of partial melts. We show that the DLF earthquakes connect the postulated upper mantle reservoir with the upper crust at a depth of about 8 km, directly below one of the youngest phonolitic volcanic centers in the Eifel, where CO(2)originating from the mantle is massively outgassing. A bright spot in the West Eifel between 6 and 10 km depth represents a Tertiary magma reservoir and is seen as a model for a differentiated reservoir beneath the young phonolitic center today. We find that the distribution of volcanic fields is controlled by the Variscan lithospheric structures and terrane boundaries as a whole, which is reflected by an offset of the Moho depth, a wedge-shaped transparent zone in the lower crust and the system of thrusts over about 120 km length.}, language = {en} } @article{ZangStephanssonStenbergetal.2017, author = {Zang, Arno and Stephansson, Ove and Stenberg, Leif and Plenkers, Katrin and von Specht, Sebastian and Milkereit, Claus and Schill, Eva and Kwiatek, Grzegorz and Dresen, Georg and Zimmermann, G{\"u}nter and Dahm, Torsten and Weber, Michael}, title = {Hydraulic fracture monitoring in hard rock at 410 m depth with an advanced fluid-injection protocol and extensive sensor array}, series = {Geophysical journal international}, volume = {208}, journal = {Geophysical journal international}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, pages = {790 -- 813}, year = {2017}, abstract = {In this paper, an underground experiment at the Aspo Hard Rock Laboratory (HRL) is described. Main goal is optimizing geothermal heat exchange in crystalline rock mass at depth by multistage hydraulic fracturing with minimal impact on the environment, that is, seismic events. For this, three arrays with acoustic emission, microseismicity and electromagnetic sensors are installed mapping hydraulic fracture initiation and growth. Fractures are driven by three different water injection schemes (continuous, progressive and pulse pressurization). After a brief review of hydraulic fracture operations in crystalline rock mass at mine scale, the site geology and the stress conditions at Aspo HRL are described. Then, the continuous, single-flow rate and alternative, multiple-flow rate fracture breakdown tests in a horizontal borehole at depth level 410 m are described together with the monitoring networks and sensitivity. Monitoring results include the primary catalogue of acoustic emission hypocentres obtained from four hydraulic fractures with the in situ trigger and localizing network. The continuous versus alternative water injection schemes are discussed in terms of the fracture breakdown pressure, the fracture pattern from impression packer result and the monitoring at the arrays. An example of multistage hydraulic fracturing with several phases of opening and closing of fracture walls is evaluated using data from acoustic emissions, seismic broad-band recordings and electromagnetic signal response. Based on our limited amount of in situ tests (six) and evaluation of three tests in Avro granodiorite, in the multiple-flow rate test with progressively increasing target pressure, the acoustic emission activity starts at a later stage in the fracturing process compared to the conventional fracturing case with continuous water injection. In tendency, also the total number and magnitude of acoustic events are found to be smaller in the progressive treatment with frequent phases of depressurization.}, language = {en} } @article{FeldMechieHuebscheretal.2017, author = {Feld, Christian and Mechie, James and H{\"u}bscher, Christian and Hall, Jeremy and Nicolaides, Stelios and Gurbuz, Cemil and Bauer, Klaus and Louden, Keith and Weber, Michael}, title = {Crustal structure of the eratosthenes seamount, cyprus and S. Turkey from an amphibian wide-angle seismic profile}, series = {Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth}, volume = {700}, journal = {Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0040-1951}, doi = {10.1016/j.tecto.2017.02.003}, pages = {32 -- 59}, year = {2017}, abstract = {In 2010, project CoCoCo (incipient COntinent-COntinent COllision) recorded a 650 km long amphibian N-S wide-angle seismic profile, extending from the Anatolian plateau across southern Turkey and Cyprus to just south of the Eratosthenes Seamount (ESM). The aim of the project is to reveal the impact of the transition from subduction to continent-continent collision of Africa with Anatolia. Arrival picking, finite-differences ray-tracing and inversion of the offshore and on-offshore data produced a tomographic model across southern Cyprus, the accretionary wedge and the ESM. The main features of this model are (1) crustal P-velocities predominantly lower than 6.5 km/s beneath the ESM, (2) crustal thickness between 28 and 37 km, (3) an upper crustal reflection at 5 km depth beneath the ESM, (4) the likely presence of oceanic crust south of the ESM and a transform margin north of it and (5) a 12 km thick ophiolite sequence on Cyprus. Land shots on Turkey, also recorded on Cyprus, gravity data and geological and previous seismic investigations allow to derive a layered velocity model beneath Anatolia and the northern part of Cyprus. The main features of this model are (1) Moho depths of 38-45 km beneath the Anatolian plateau, (2) an upper and lower crust with large lateral changes in velocity and thickness, (3) a north-dipping subducting plate below Cyprus with a steepening of the dip-angle of the plate at about 45 km depth. Thus, the wide-angle seismic and gravity data provide detailed insights into the 2-D geometry and velocity structures associated with the Cyprus Arc collision zone. Finally, integrated analysis of the geophysics and geology allows a comprehensive interpretation of the crustal structure related to the collision process.}, language = {en} } @article{HendriyanaBauerMuksinetal.2018, author = {Hendriyana, Andri and Bauer, Klaus and Muksin, Umar and Weber, Michael}, title = {AIC-based diffraction stacking for local earthquake locations at the Sumatran Fault (Indonesia)}, series = {Geophysical journal international}, volume = {213}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggy045}, pages = {952 -- 962}, year = {2018}, abstract = {We present a new workflow for the localization of seismic events which is based on a diffraction stacking approach. In order to address the effects from complex source radiation patterns, we suggest to compute diffraction stacking from a characteristic function (CF) instead of stacking the original waveform data. A new CF, which is called in the following mAIC (modified from Akaike Information Criterion) is proposed. We demonstrate that both P- and S-wave onsets can be detected accurately. To avoid cross-talk between P and S waves due to inaccurate velocity models, we separate the P and S waves from the mAIC function by making use of polarization attributes. Then, the final image function is represented by the largest eigenvalue as a result of the covariance analysis between P-and S-image functions. Results from synthetic experiments show that the proposed diffraction stacking provides reliable results. The workflow of the diffraction stacking method was finally applied to local earthquake data from Sumatra, Indonesia. Recordings from a temporary network of 42 stations deployed for nine months around the Tarutung pull-apart basin were analysed. The seismic event locations resulting from the diffraction stacking method align along a segment of the Sumatran Fault. A more complex distribution of seismicity is imaged within and around the Tarutung basin. Two lineaments striking N-S were found in the centre of the Tarutung basin which support independent results from structural geology.}, language = {en} } @article{JoussetReinschRybergetal.2018, author = {Jousset, Philippe and Reinsch, Thomas and Ryberg, Trond and Blanck, Hanna and Clarke, Andy and Aghayev, Rufat and Hersir, Gylfi P. and Henninges, Jan and Weber, Michael and Krawczyk, Charlotte M.}, title = {Dynamic strain determination using fibre-optic cables allows imaging of seismological and structural features}, series = {Nature Communications}, volume = {9}, journal = {Nature Communications}, publisher = {Nature Publ. Group}, address = {London}, issn = {2041-1723}, doi = {10.1038/s41467-018-04860-y}, pages = {11}, year = {2018}, abstract = {Natural hazard prediction and efficient crust exploration require dense seismic observations both in time and space. Seismological techniques provide ground-motion data, whose accuracy depends on sensor characteristics and spatial distribution. Here we demonstrate that dynamic strain determination is possible with conventional fibre-optic cables deployed for telecommunication. Extending recently distributed acoustic sensing (DAS) studies, we present high resolution spatially un-aliased broadband strain data. We recorded seismic signals from natural and man-made sources with 4-m spacing along a 15-km-long fibre-optic cable layout on Reykjanes Peninsula, SW-Iceland. We identify with unprecedented resolution structural features such as normal faults and volcanic dykes in the Reykjanes Oblique Rift, allowing us to infer new dynamic fault processes. Conventional seismometer recordings, acquired simultaneously, validate the spectral amplitude DAS response between 0.1 and 100 Hz bandwidth. We suggest that the networks of fibre-optic telecommunication lines worldwide could be used as seismometers opening a new window for Earth hazard assessment and exploration.}, language = {en} }