ALBERT

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉Spectral induced polarization (SIP) laboratory measurements on water‐saturated rocks show a strong correlation between the electrical polarization strength and the inner surface area of rocks. We investigate the influence of inner surface roughness on the SIP response by simulating the frequency‐dependent complex conductivity of micro‐scale rock models. Starting with smooth grain models, we introduce surface roughness using two different approaches: increasing the surface roughness in a fractal‐like manner, and creating random surface structures, resulting in more natural‐looking surfaces. We find that surface roughness has two distinct effects on the SIP response: (a) a shift in the position and magnitude of the primary relaxation frequency to lower frequencies and lower magnitudes, respectively, and (b) the formation of secondary polarizations above the polarization frequency of the primary polarization. We also compare the relaxation time and normalized chargeability obtained by Debye decomposition and the imaginary conductivity at 1 Hz of our models with mechanistic models and empirical relations. We point out the congruences and offer explanations for the discrepancies between our models and the empirical observations. We conclude that the results of our study are applicable to real rocks and that the SIP method has the potential to detect inner surface roughness. However, the SIP method it not able to discriminate between signals from rough particles and a distribution of smooth particles.〈/p〉

Description: Plain Language Summary: The spectral induced polarization method measures the electrical conductivity of the subsurface at depths of investigation ranging from a few dm to several 100 m. The potential of the method to detect mineral resources, contamination, microbial activity, etc. makes it a promising tool for today's environmental challenges. However, due to the complexity of the underlying physical and chemical processes the interpretation of SIP measurements often remains qualitative. Mathematical and mechanistic models used to describe the processes are usually simplified with respect to geometry, chemistry, and physical properties. In this study, we use computer simulations to investigate the surfaces of spherical particles (grains). After simulating the SIP response of models with smooth and rough surfaces, we analyze the influence of surface roughness on the SIP response. We find that surface roughness causes distinct changes in the SIP response compared to the smooth grain. However, we also find that this introduces an ambiguity in the interpretation of SIP data.〈/p〉

Description: Key Points: 〈list list-type="bullet"〉 〈list-item〉 〈p xml:lang="en"〉Inner surface roughness of rocks has substantial impact on spectral induced polarization response〈/p〉〈/list-item〉 〈list-item〉 〈p xml:lang="en"〉Surface roughness shifts the primary polarization peak to lower frequencies〈/p〉〈/list-item〉 〈list-item〉 〈p xml:lang="en"〉Surface roughness causes additional polarizations above the primary peak frequency〈/p〉〈/list-item〉 〈/list〉 〈/p〉

Description: https://doi.org/10.5281/ZENODO.7049722

Description: Publishing data with a DOI (Digital Object Identifier) is becoming increasingly popular, and many researchers are being asked by journals to publish the data underlying the research results described in the article. Where does the data come from and how do I cite it correctly? […] This publication is a translation of the FID GEO article in 〈a href="http://doi.org/10.23689/fidgeo-5688"〉GMIT 90.〈/a〉

Description: research

Description: The ORCID is an open alphanumeric identifier used to uniquely identify researchers.In an increasingly connected scientific landscape, the growing number of published research outputs (data, text, software) and new open publishing models pose major challenges for scientific publication management. The repositories available to FID GEO, GEO-LEOe-docs (SUB Goettingen) and GFZ Data Services (GFZ), have long been working with PIDs such as the Digital Object Identifier (DOI ) and the Open Researcher and Contributor iD (ORCID iD ), and recommend the use of the ORCID. This publication is a translation of the FID GEO article in 〈a href="http://doi.org/10.23689/fidgeo-5832"〉GMIT 87.〈/a〉

Description: DFG, SUB Göttingen, GFZ Potsdam

Description: research

Description: [...] ROR stands for Research Organization Registry and is an open registry for research organizations. The registry is developed and operated as a joint initiative of the California Digital Library, Crossref and DataCite. ROR already references more than 103,000 organizations, ranging from universities and extramural research institutions to government agencies, research funders and government departments. [...] This publication is a translation of the FID GEO article in 〈a href="http://doi.org/10.23689/fidgeo-5750"〉GMIT 91.〈/a〉

Description: research

Description: Oceanic circulation and mass‐field variability play important roles in exciting Earth's wobbles and length‐of‐day changes (ΔΛ), on time scales from days to several years. Modern descriptions of these effects employ oceanic angular momentum (OAM) series from numerical forward models or ocean state estimates, but nothing is known about how ocean reanalyses with sequential data assimilation (DA) would fare in that context. Here, we compute daily OAM series from three 1/4° global ocean reanalyses that are based on the same hydrodynamic core and input data (e.g., altimetry, Argo) but different DA schemes. Comparisons are carried out (a) among the reanalyses, (b) with an established ocean state estimate, and (c) with Earth rotation data, all focusing on the period 2006–2015. The reanalyses generally provide credible OAM estimates across a range of frequencies, although differences in amplitude spectra indicate a sensitivity to the adopted DA scheme. For periods less than 120 days, the reanalysis‐based OAM series explain ∼40%–50% and ∼30%–40% of the atmosphere‐corrected equatorial and axial geodetic excitation, similar to what is achieved with the state estimate. We find mixed performance of the reanalyses in seasonal excitation budgets, with some questionable mean ocean mass changes affecting the annual cycle in ΔΛ. Modeled excitations at interannual frequencies are more uncertain compared to OAM series from the state estimate and show hints of DA artifacts in one case. If users are to choose any of the tested reanalyses for rotation research, our study points to the Ocean Reanalysis System 5 as the most sensible choice.

Description: Key Points: We evaluate three ocean reanalyses for their skill in explaining Earth rotation variations on different time scales from 2006 to 2015. For periods 〈120 days, reanalyses explain 40%–50% of atmosphere‐reduced polar motion excitation variance, similar to an ocean state estimate. Reanalyses show mixed skill in seasonal excitation budgets and, in one case, hints of data assimilation artifacts at interannual periods.

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: https://resources.marine.copernicus.eu/product-detail/GLOBAL_REANALYSIS_PHY_001_031/INFORMATION

Description: https://isdc.gfz-potsdam.de/ggfc-oceans/oam/

Description: https://www.ncei.noaa.gov/access/metadata/landing-page/bin/iso?id=gov.noaa.ngdc.mgg.dem:316

Description: https://podaac-tools.jpl.nasa.gov/drive/files/GeodeticsGravity/tellus/L3/mascon/RL06/JPL/v02/CRI/netcdf

Description: https://keof.jpl.nasa.gov/combinations/

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉Benthic oxygen dynamics and the exchange of oxygen and other solutes across the sediment‐water interface play a key role for the oxygen budget of many limnic and shallow marine systems. The sediment‐water fluxes are largely determined by two factors: sediment biogeochemistry and the thickness of the diffusive boundary layer that is determined by near‐bottom turbulence. Here, we present a fully coupled benthic‐pelagic modeling system that takes these processes and their interaction into account, focusing especially on the modulation of the sediment‐water fluxes by the effects of near‐bottom turbulence and stratification. We discuss the special numerical methods required to guarantee positivity and mass conservation across the sediment‐water interface in the presence of rapid element transformation, and apply this modeling system to a number of idealized scenarios. Our process‐oriented simulations show that near‐bottom turbulence provides a crucial control on the sediment‐water fluxes, the oxygen penetration depth, and the re‐oxidation of reduced compounds diffusing upward from the deeper benthic layers especially on time scales of a few days, characterizing oceanic tides, internal seiching motions in lakes, and mesoscale atmospheric variability. Our results also show that the response of benthic‐pelagic fluxes to rapid changes in the forcing conditions (e.g., storm events) can only be understood with a fully coupled modeling approach.〈/p〉

Description: Plain Language Summary: Oxygen is one of the most relevant ecosystem parameters in marine systems and in lakes. In shallow systems, the overall oxygen budget is often controlled by the sedimentary oxygen demand, and it is therefore crucial to account for the exchange of oxygen and other solutes between the sediment and the water column. In this, context, a thin (millimeter‐scale) “diffusive sublayer,” located immediately above the sediment surface, is known to play an especially important role as it may form a bottleneck for the solute exchange. As the thickness of this sublayer is controlled by hydrodynamic processes, sediment‐water fluxes are affected by complex feedbacks between physical and biogeochemical processes. Here, we describe a fully coupled numerical modeling system that takes these feedback mechanisms into account, including advanced numerical methods guaranteeing that the total mass of all solutes is conserved (even if these are transformed) and that their concentrations do not become negative. Using a series of idealized examples, it is shown that near‐bottom hydrodynamic processes have an important impact on the sediment‐water fluxes, the depth to which oxygen penetrates into the upper sediment layers, and the re‐oxidation of reduced chemical compounds in the sediments. These feedbacks are particularly important for processes with time scales of a few days, like ocean tides, internal oscillations in lakes, and short‐term atmospheric disturbances (e.g., storm events).〈/p〉

Description: Key Points: First fully coupled benthic‐pelagic modeling system accounting for the effects of near‐bottom turbulence on sediment‐water solute fluxes. Hydrodynamic effects control benthic biogeochemistry and fluxes especially on time scales of a few days, and during extreme events. New numerical methods that guarantee mass conservation and positivity across the sediment‐water interface.

Description: Leibniz Association

Description: German Research Foundation

Description: Swiss National Science Foundation

Description: UK Natural Environment Research Council

Description: Department for Environment Food and Rural Affairs

Description: https://gotm.net/

Description: https://github.com/fabm-model

Description: https://doi.org/10.5281/zenodo.7950383

Description: https://doi.org/10.5281/zenodo.7950866

Description: Methane seeps are typified by the formation of authigenic carbonates, many of which exhibit corrosion surfaces and secondary porosity believed to be caused by microbial carbonate dissolution. Aerobic methane oxidation and sulfur oxidation are two processes capable of inducing carbonate corrosion at methane seeps. Although the potential of aerobic methanotrophy to dissolve carbonate was confirmed in laboratory experiments, this process has not been studied in the environment to date. Here, we report on a carbonate corrosion experiment carried out in the REGAB Pockmark, Gabon‐Congo‐Angola passive margin, in which marble cubes were deployed for 2.5 years at two sites (CAB‐B and CAB‐C) with apparent active methane seepage and one site (CAB‐D) without methane seepage. Marble cubes exposed to active seepage (experiment CAB‐C) were found to be affected by a new type of microbioerosion. Based on 16〈italic toggle="no"〉S r〈/italic〉RNA gene analysis, the biofilms adhering to the bioeroded marble mostly consisted of aerobic methanotrophic bacteria, predominantly belonging to the uncultured Hyd24‐01 clade. The presence of abundant 〈sup〉13〈/sup〉C‐depleted lipid biomarkers including fatty acids (〈italic toggle="no"〉n〈/italic〉‐C〈sub〉16:1ω8c〈/sub〉, 〈italic toggle="no"〉n〈/italic〉‐C〈sub〉18:1ω8c〈/sub〉, 〈italic toggle="no"〉n〈/italic〉‐C〈sub〉16:1ω5t〈/sub〉), various 4‐mono‐ and 4,4‐dimethyl sterols, and diplopterol agrees with the dominance of aerobic methanotrophs in the CAB‐C biofilms. Among the lipids of aerobic methanotrophs, the uncommon 4α‐methylcholest‐8(14)‐en‐3β,25‐diol is interpreted to be a specific biomarker for the Hyd24‐01 clade. The combination of textural, genetic, and organic geochemical evidence suggests that aerobic methanotrophs are the main drivers of carbonate dissolution observed in the CAB‐C experiment at the REGAB pockmark.〈/p〉

Description: Deutscher Akademischer Austauschdienst http://dx.doi.org/10.13039/501100001655

Description: thesis

Description: The distance to failure of the upper crustal rock in the prevalent stress field is of importance to better understand fault reactivation by natural and induced processes as well as to plan and manage georeservoirs. In particular, the contemporary stress state is one of the key ingredients for this assessment. To provide a continuous description of the 3D absolute stress state geomechanical‐numerical models are used. However, stress magnitude data for model calibration are sparse and incomplete and thus, the resulting model uncertainties are large. In order to reduce the uncertainties, we incorporate additional constraints on stress magnitudes to check the plausibility of different data‐based stress states. We use formation integrity tests, borehole breakouts, drilling induced fractures, and observations of seismicity and distinct seismological quiescence. This information is weighted according to its confidence and the agreement with the different modeled stress states is assessed. The information is introduced to a Bayesian approach to estimate weights of the modeled stress states and thereby identify their plausibility. A case study in southern Germany shows the ability of the approach to identify from a wide range of stress states a small number of plausible ones and reject implausible stress states. This significantly reduces the number of stress states and thus lowers the model uncertainties.

Description: Plain Language Summary: The upper crust of the Earth (upper few kilometers) is subject to a kind of pressure, referred to as stress. When the stress becomes larger than the strength of the rock, the rock breaks. Sometimes this can be measured or even felt as a seismic event. It happens naturally but may also happen due to human activity. To prevent such induced seismic event, it is important to know the stress state. But there is only few information on the magnitude of the stress so we need computer models to predict the stress state. These models are often not very precise since there is only few information on the stress magnitudes and in addition they are often contradicting. We use all stress magnitude information individually to model various stress states. Then we look at other information that is related to the stress state but does not provide stress magnitude information on its own. We compare this information with the modeled stress states to find out whether a stress state agrees with the additional information or not. This allows us to identify a few realistic stress state models out of a wide range of possible ones. This reduces the uncertainties of the stress predictions.

Description: Key Points: Bayesian approach to uncertainty quantification and reduction of 3D geomechanical‐numerical models of the undisturbed stress state. Additional constraints on the stress state by formation integrity tests, borehole breakouts, and drilling induced tensile fractures. Constraints on the differential stress by observation of seismicity or distinct seismological quiescence.

Description: Bundesministerium für Bildung, Wissenschaft und Forschung http://dx.doi.org/10.13039/501100013699

Description: Helmholtz Association http://dx.doi.org/10.13039/501100009318

Description: Helmholtz Centre Potsdam ‐ Deutsches GeoForschungsZentrum GFZ

Description: RI Fabrice Cotton

Description: Federal Company for radioactive waste disposal

Description: https://github.com/MorZieg/FAST_Calibration

Description: Spaceborne impact ionization mass spectrometers, such as the Cosmic Dust Analyzer on board the past Cassini spacecraft or the SUrface Dust Analyzer being built for NASA's upcoming Europa Clipper mission, are of crucial importance for the exploration of icy moons in the Solar System, such as Saturn's moon Enceladus or Jupiter's moon Europa. For the interpretation of data produced by these instruments, analogue experiments on Earth are essential. To date, thousands of laboratory mass spectra have been recorded with an analogue experiment for impact ionization mass spectrometers. Simulation of mass spectra of ice grains in space is achieved by a Laser Induced Liquid Beam Ion Desorption (LILBID) approach. The desorbed cations or anions are analyzed in a time‐of‐flight mass spectrometer. The amount of unstructured raw data is increasingly challenging to sort, process, interpret and compare with data from space. Thus far this has been achieved manually for individual mass spectra because no database containing the recorded reference spectra was available. Here we describe the development of a comprehensive, extendable database containing cation and anion mass spectra from the laboratory LILBID facility. The database is based on a Relational Database Management System with a web server interface and enables filtering of the laboratory data using a wide range of parameters. The mass spectra can be compared not only with data from past and future space missions but also mass spectral data generated by other, terrestrial, techniques. The validated and approved subset of the database is available for general public (https://lilbid-db.planet.fu-berlin.de).

Description: Plain Language Summary: Thousands of laboratory mass spectra, each with an individual set of experimental parameters, have been recorded so far using a facility situated at Freie Universität Berlin, Germany. The mass spectra help analyze and interpret data returned from spacecraft in the vicinity of icy moons in the Solar System. The unstructured laboratory data is increasingly challenging to sort and compare to the data from space. We developed an extendable database containing the laboratory data. The database is available for general public and allows filtering the stored data for a wide range of experimental parameters and, in turn, significantly improves analysis of data not only from past space missions but also future missions in particular.

Description: Key Points: We describe the development of a comprehensive spectral database containing laboratory analogue data for spaceborne mass spectrometers. The database is based on a Relational Database Management System with a web interface and accessible for community use. Filtering the laboratory data using a wide range of experimental parameters allows a straightforward analysis of returned flight data.

Description: EC, Horizon 2020 Framework Programme (H2020) http://dx.doi.org/10.13039/100010661

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: https://doi.org/10.5281/zenodo.6863855

Description: https://sbnarchive.psi.edu/pds3/cassini/cda/COCDA_0007.tar.gz

Description: poster

Description: conference

Description: It is generally agreed that the resolution of a regular quadrilateral mesh is the side length of quadrilateral cells. There is less agreement on the resolution of triangular meshes, exacerbated by the fact that the numbers of edges or cells on triangular meshes are approximately three or two times larger than that of vertices. However, the geometrical resolution of triangular meshes, that is, maximum wavenumbers or smallest wavelengths that can be represented on such meshes, is a well defined quantity, known from solid state physics. These wavenumbers are related to a smallest common mesh cell (primitive unit cell), and the set of mesh translations that map it into itself. They do not depend on whether discrete degrees of freedom are placed on vertices, cells or edges. For equilateral triangles the smallest wavelength equals twice the triangle height. Resolutions of quadrilateral and triangular meshes approximately agree if they have the same numbers of vertices.

Description: Plain Language Summary: Some models used in climate studies are formulated on triangular computational meshes. We discuss how to determine the smallest scales that are resolved on such meshes. They are referred to as a mesh resolution. The notion of mesh resolution is commonly used to relate climate model results simulated on different meshes.

Description: Key Points: Geometrical resolution of an equilateral triangular mesh is defined by the height of its triangles. Quadrilateral and triangular meshes with the same number of vertices have approximately the same resolution.

Description: Collaborative Research Centre

Description: German Research Foundation

Description: Efficient compositional models are required to simulate underground gas storage in porous formations where, for example, gas quality (such as purity) and loss of gas due to dissolution are of interest. We first extend the concept of vertical equilibrium (VE) to compositional flow, and derive a compositional VE model by vertical integration. Second, we present a hybrid model that couples the efficient compositional VE model to a compositional full‐dimensional model. Subdomains, where the compositional VE model is valid, are identified during simulation based on a VE criterion that compares the vertical profiles of relative permeability at equilibrium to the ones simulated by the full‐dimensional model. We demonstrate the applicability of the hybrid model by simulating hydrogen storage in a radially symmetric, heterogeneous porous aquifer. The hybrid model shows excellent adaptivity over space and time for different permeability values in the heterogeneous region, and compares well to the full‐dimensional model while being computationally efficient, resulting in a runtime of roughly one‐third of the full‐dimensional model. Based on the results, we assume that for larger simulation scales, the efficiency of this new model will increase even more.

Description: Key Points: A compositional vertical equilibrium model is coupled to its full‐dimensional counterpart. A criterion is developed to adaptively identify and assign regions where the vertical equilibrium model is applicable during simulation. A test case of hydrogen storage in a heterogeneous porous aquifer demonstrates efficiency and accuracy of the hybrid model.

Description: Deutsche Forschungsgemeinschaft (DFG) http://dx.doi.org/10.13039/501100001659

Description: https://git.iws.uni-stuttgart.de/dumux-pub/Becker2021b.git

Description: Basal melting of marine‐terminating glaciers, through its impact on the forces that control the flow of the glaciers, is one of the major factors determining sea level rise in a world of global warming. Detailed quantitative understanding of dynamic and thermodynamic processes in melt‐water plumes underneath the ice‐ocean interface is essential for calculating the subglacial melt rate. The aim of this study is therefore to develop a numerical model of high spatial and process resolution to consistently reproduce the transports of heat and salt from the ambient water across the plume into the glacial ice. Based on boundary layer relations for momentum and tracers, stationary analytical solutions for the vertical structure of subglacial non‐rotational plumes are derived, including entrainment at the plume base. These solutions are used to develop and test convergent numerical formulations for the momentum and tracer fluxes across the ice‐ocean interface. After implementation of these formulations into a water‐column model coupled to a second‐moment turbulence closure model, simulations of a transient rotational subglacial plume are performed. The simulated entrainment rate of ambient water entering the plume at its base is compared to existing entrainment parameterizations based on bulk properties of the plume. A sensitivity study with variations of interfacial slope, interfacial roughness and ambient water temperature reveals substantial performance differences between these bulk formulations. An existing entrainment parameterization based on the Froude number and the Ekman number proves to have the highest predictive skill. Recalibration to subglacial plumes using a variable drag coefficient further improves its performance.

Description: Plain Language Summary: In a world of global warming, the melting of glaciers terminating as floating ice tongues into the oceans of Arctic and Antarctic regions allows those glaciers to flow faster and hence to make a considerable contribution to global mean sea‐level rise. Underneath the ice‐ocean interface, turbulent currents of the order of 10 m thickness (so‐called plumes) develop that transport the melt water from the grounding line where the glacier enters the ocean toward the calving front that marks the seaward end of the glacier. At its base, ambient relatively warm and salty ocean water is mixed into the plumes and is vertically transported toward the ice‐ocean interface, where the melting is increased due to the additional heat supply. Understanding these processes is essential for their incorporation into computer models for the prediction of such melt processes. In this study, an accurate simulation model for the water column is constructed that is able to consistently reproduce these processes. The algorithms developed here are proven to provide reliable results also for models with only a few grid points across the plume and can therefore be implemented into climate models with surface‐following coordinates to more accurately simulate future scenarios of sea level rise.

Description: Key Points: A vertically resolving model with second‐moment turbulence closure has been constructed for subglacial plumes. Convergent numerical formulations for the ocean‐to‐ice fluxes of momentum, freshwater and heat have been derived from an analytical model. Model results are consistent with bulk parameterizations for the entrainment of ambient water.

Description: Bundesministerium für Bildung und Forschung (BMBF) http://dx.doi.org/10.13039/501100002347

Description: https://doi.org/10.5281/zenodo.6203838

Description: Hydrogen is a promising alternative to carbon based energy carriers and may be stored in large quantities in subsurface storage deposits. This work assesses the impact of static (density and phase equilibria) and dynamic (viscosity and diffusion coefficients) properties on the pressure field during the injection and extraction of hydrogen in the porous subsurface. In a first step, we derive transport properties for water, hydrogen and their mixture using the Perturbed‐Chain Statistical Associating Fluid Theory equation of state in combination with an entropy scaling approach and compare model predictions to alternative models from the literature. Our model compares excellently to experimental transport coefficients and models from literature with a higher number of adjustable parameters, such as GERG2008, and shows a clear improvement over empirical correlations for transport coefficients of hydrogen. In a second step, we determine the effect of further model reduction by comparing our against a much simpler model applying empirical transport coefficients from the literature. For this purpose, hydrogen is periodically injected into and extracted out of a dome‐shaped porous aquifer under a caprock. Our results show that density and viscosity of hydrogen have the highest impact on the pressure field, and that a thermodynamic model like the new model presented here is essential for modeling the storage aquifer, while keeping the number of coefficients at a minimum. In diffusion‐dominated settings such as the diffusion of hydrogen through the caprock, our developed diffusion coefficients show a much improved dependence on temperature and pressure, leading to a more accurate approximation of the diffusive fluxes.

Description: Key Points: We model the phase behavior of pure hydrogen and the binary hydrogen‐water mixture using the Perturbed‐Chain Statistical Associating Fluid Theory equation of state. New entropy scaling relations for the transport properties of hydrogen and water and diffusion coefficients of their mixture are derived. The impact of the newly derived fluid properties is analyzed for a scenario of hydrogen storage in a porous aquifer.

Description: Deutsche Forschungsgemeinschaft (DFG) http://dx.doi.org/10.13039/501100001659

Description: https://git.iws.uni-stuttgart.de/dumux-pub/sauerborn2020a

Description: Vlasov solvers that operate on a phase‐space grid are highly accurate but also numerically demanding. Coarse velocity space resolutions, which are largely unproblematic in particle‐in‐cell (PIC) simulations, can lead to numerical heating or oscillations in continuum Vlasov methods. To address this issue, we present a new dual Vlasov solver which is based on an established positivity preserving advection scheme for the update of the distribution function and an energy conserving partial differential equation solver for the kinetic update of mean velocity and temperature. The solvers work together via moment fitting during which the maximum entropy part of the distribution function is replaced by the solution from the partial differential equation solver. This numerical scheme makes continuum Vlasov methods competitive with PIC methods concerning computational cost and enables us to model large scale reconnection in Earth's magnetosphere with a fully kinetic continuum method. The simulation results agree well with measurements by the MMS spacecraft.

Description: Key Points: A moment fitting continuum Vlasov solver is presented that preserves positivity of the distribution function and conserves total energy. The method behaves well at low velocity space resolutions, making it competitive with PIC methods concerning computational cost. There is good agreement of the simulations with measurements of magnetic reconnection by the MMS spacecraft.

Description: Helmholtz Association (亥姆霍兹联合会致力) http://dx.doi.org/10.13039/501100009318

Description: https://vlasov.tp1.ruhr-uni-bochum.de/data/paper-JGR-2021

Description: In recent years, Spatial Markov Models have gained popularity in simulating solute transport in heterogeneous formations. They describe the transition times of particles between equidistant observation planes by statistical distributions, assuming correlation of the transit times of individual particles between subsequent steps. By this, the approach naturally captures preasymptotic solute dispersion. In this study, we analyze Spatial Markov Models assuming bivariate log‐normal distributions of the particle slowness (i.e., the inverse velocity) in subsequent transitions. The model is fully parameterized by the mean Eulerian velocity, the variance of the log‐slowness, and the correlation coefficient of log‐slowness in subsequent steps. We derive closed‐form expressions for distance‐dependent ensemble dispersion, which is defined in terms of the second‐central moments of the solute breakthrough curves. We relate the coefficients to the properties of the underlying log‐hydraulic conductivity field assuming second‐order stationarity. The results are consistent with linear stochastic theory in the limit of small log‐conductivity variances, while the approach naturally extends to high‐variance cases. We demonstrate the validity of the approach by comparison to three‐dimensional particle‐tracking simulations of advective transport in heterogeneous media with isotropic, exponential correlation structure for log‐conductivity variances up to five. This study contributes to relating solute dispersion to metrics of the porous‐medium structure in cases of strong heterogeneity.

Description: Key Points: We derive closed‐form expressions of ensemble dispersion in the spatial‐Markov framework of solute transport. The expressions are consistent with linear theory in the limit of small log‐conductivity variances, but extend to high‐variance cases. Comparison to particle‐tracking simulations of advective transport in 3‐D heterogeneous domains show excellent agreement.

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: https://doi.org/10.5281/zenodo.6554308

Description: https://www.hsl.rl.ac.uk/catalogue/hsl_mi20.html

Description: Atmospheric delay corrections for satellite altimetry measurements are essential for deriving highly accurate sea surface heights and reliable global mean sea level (GMSL) trend estimates. A commonly used method to correct for ionospheric path delays are the usage of GNSS‐based Global Ionospheric Maps (GIM). The different orbit heights of GNSS and altimeter satellites require an adaption of GIM corrections to account for free electrons in the Earth plasmasphere. This study shows that the widely used scaling approach based on the International Reference Ionosphere (IRI) is not able to accurately scale the GIM models. The impact of neglecting the plasmaspheric part of the atmosphere strongly correlates with the solar activity of about 11 years. This manifests itself as trend errors in global GMSL. For the Jason period (2002–2021) a trend error of 0.17 mm/year can be shown, which is even larger for smaller periods (e.g., 1.0 mm/year for Jason‐1 lifetime). The application of an additional constant scaling factor of 0.886 can reduce the trend differences to below 0.05 mm/year.

Description: Plain Language Summary: Global mean sea level (GMSL) rise is an important indicator for climate change. To precisely measure this quantity that is only in the order of about 3 mm/year, satellite altimeters are used. Their observations have to be corrected for influences in the Earth atmosphere. This study shows deficiencies in one commonly used correction data set. These corrections, based on observations from the Global Navigation Satellite Systems (GNSS) are not accounting for the higher part of the atmosphere, the plasmasphere. Neglecting this influence derives systematic errors with a 11 years cycle that impacts the estimation of GMSL by up to 1 mm/year, depending on the period under investigation. It is recommended to apply an additional scaling of the available corrections in order to reduce the trend error to below 0.05 mm/year.

Description: Key Points: Global Ionospheric Map corrections in altimetry Sensor Geophysical Data Records are not fully scaled to account for plasmaspheric electron content. Neglecting the plasmaspheric effect leads to trends of up to 1 mm/year in Global mean sea level estimates. The additional application of a scale factor improves the consistency in trend with respect to dual‐frequency satellite altimetry data.

Description: https://podaac.jpl.nasa.gov/

Description: https://www.aviso.altimetry.fr/

Description: https://openadb.dgfi.tum.de/en/products/vertical-total-electron-content/

Description: The subduction of serpentinized slabs is the dominant process to transport “water” into Earth's mantle, and plays a pivotal role for subduction dynamics. Antigorite, the most abundant serpentine mineral in subduction settings, may imprint a seismic signature on serpentinized slabs, making them seismically distinguishable from the dry, non‐serpentinized ones. However, the complete single‐crystal elasticity of antigorite has not been experimentally constrained at high pressures, hindering the use of seismological approaches to detect serpentinization in subducting slabs. Here, we report the full elastic stiffness tensor of antigorite by single‐crystal Brillouin spectroscopy and X‐ray diffraction up to 7.71(5) GPa. We use our results to model seismic properties of antigorite‐bearing rocks and show that their seismological detectability depends on the geometrical relation between seismic wave paths and foliation of serpentinized rocks. In particular, we demonstrate that seismic shear anisotropy shows low sensitivity to serpentinization for a range of relevant geometries.

Description: Plain Language Summary: The subduction of serpentinized slabs plays a key role in the deep recycling of water into the Earth's interior. Antigorite is the main serpentine mineral in subducting slabs, and the most important carrier of water. Antigorite‐bearing rocks are predicted to have a distinct seismic signature, potentially allowing them to be detected with seismological approaches. However, our current knowledge on seismic properties of antigorite‐bearing rocks is limited, mostly hampered by a lack of experimental constraints on single‐crystal elasticity of antigorite at relevant pressures. In this study, state‐of‐the‐art techniques were employed to produce the first experimental description of the complete high‐pressure elasticity of antigorite single crystals. Our experimental data set was implemented in the modeling of seismic properties of antigorite‐bearing rocks at pressures relevant for subduction. Our results were used to discuss the relation between seismic wave path and shear wave anisotropy in serpentinized slabs, and challenge the use of shear wave splitting as a proxy for serpentinization in slabs.

Description: Key Points: Single‐crystal elasticity of antigorite at high pressures is determined by Brillouin spectroscopy and X‐ray diffraction experiments. Seismic signature of serpentinized slabs is constrained in a relevant composition‐pressure space. Serpentinization in slabs may be undetectable through shear wave anisotropy.

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: EC Horizon 2020 Framework Programme http://dx.doi.org/10.13039/100010661

Description: https://doi.org/10.6084/m9.figshare.20348748

Description: https://doi.org/10.6084/m9.figshare.20348781

Description: Geological maps are complex to produce through intensive and expensive field studies. Comparisons of geophysical data with geological conditions are difficult and often only qualitatively possible. The following work therefore examines an automated procedure to better reconcile this information. For this purpose, the terracing method, and a cluster analysis of potential field (gravity and magnetic field) and petrophysical data from the Karasjok and Ligurian Sea regions are used to interpret this geophysical measurements in a geological way. Two different tectonic regions were selected: (1) The Karasjok region is located in Northern Norway, where the Karasjok Greenstone Belt (KGB) dominates geological settings, consisting of abundant ultramafic intrusions, komatiites, gabbroic intrusions, amphibolites and migmatites. (2) The Ligurian-Provençal basin, part of the Western Mediterranean Sea, which is located between the French-Italian coastline and the island of Corsica. Geologically the area is characterised by the spreading zone in the Western Mediterranean. The high-resolution Airborne Gravity Gradient Survey and aeromagnetic datasets of the Karasjok region cover an area of 20 km x 30 km with a data resolution of 50 m. The dataset of the Ligurian basin cover a much larger area with the resolution of 5 km. Data constraints come from former LOBSTER and LISA campaigns and a study in the research group at CAU Kiel, new compilation of the AlpArray Gravity Research Group (AAGRG), besides data of the ICGEM Potsdam (disturbance) and the GOCE mission. By aid of the terracing algorithm, the boundaries of the anomalies are to be sharpened and regions with constant field amplitude were generated. For this purpose, a shape index-based algorithm was applied, which uses the shape index calculated at each field point to grade the function. Through an iterative process and the variation of parameters, the terracing result is refined. The resulting data sets are then further processed using a cluster analysis method. Here, the k-mean algorithm for domain classification is used to divide the geophysical measurement data into groups (cluster) of similar properties. The number of clusters k is specified and the data points are assigned to the respective clusters through an iterative process. Using the data of the datasets mentioned above the results of this applications are successfully compared with the corresponding geological maps of the two areas.

Description: Deutsche Forschungsgemeinschaft (DFG), SPP "MB-4D"

Description: https://doi.org/10.5880/fidgeo.2020.045 (Zahorec et al., 2021) via GFZ Data Services

Description: poster

Description: Aeromagnetic surveys help us to learn about geology. To achieve good coverage, surveys need to be merged. However, conventional methods introduce long-wavelength bias and cannot handle the individual survey quality. We develop a new approach to process large aeromagnetic surveys with an equivalent layer approach and combine them with satellite data. To facilitate the usage of large data sets, we divide the study area into blocks and treat each block individually. We adjust the block size according to the resolution of the equivalent source model. Within each block we solve for equivalent sources using an iterative linear inversion with Tikhonov regularization. We apply a multi-resolution strategy by iteratively decreasing the dipole spacing, dipole depth and block size. In each step, the resolution is applied to the residual of the previous steps. This ensures both a good representation of the large and small-scale structures as well as reasonable computational costs. Advantages of the blockwise inversion are the handling with large data sets due to splitting up the study area and neglecting influences of sources above a certain distance. This reduces computational costs and still fits the data well in comparison with an unblocked inversion. Some structures cannot be resolved well with just one dipole layer, so the multi-resolution strategy enables to have a better fit by separating regional and local sources. For the final compilation, we replace the long wavelengths part of the aeromagnetic data with satellite data to spherical harmonic degree 110. We demonstrate our new approach with a newly compiled large data base for Greenland.

Description: poster

Description: Recent years have shown an increased interest in Polar research and in particular in understanding tectonics and seismic hazard in the Arctic. To understand the seismic activity in the European Arctic, the seismic bulletins should be as complete as possible. We present a new seismic event bulletin for the European Arctic (70° – 90° N, -15° – 75° E), for the 24-year long period 1990 – 2013. The poster will show in detail the merging of the different sources taken in account for the compilation, the homogenization of the data and the relocation of the seismic events. With respect to the ISC bulletin for this region, the new bulletin contains 5,932 new seismic events and 54,630 new seismic onset readings from stations mostly located at regional distances. The gains are distributed over the entire study region, with the most significant contributions across the Svalbard Archipelago, along the Knipovich and northern Mohns Ridges, as well as northern Fennoscandia.

Description: Norwegian Research Council Grant 233973/H30

Description: Russian Foundation for Basic Research Grant 14-05-93080

Description: Russian Foundation for Basic Research Grant 18-05-70018

Description: https://doi.org/10.31905/TYLLQY8T

Description: poster

Description: During the last few years, the use of Ground-Penetrating-Radar (GPR) multi-channel antenna arrays in the Archaeological Prospection increased dramatically. The main advantage of this type of survey is a much faster data acquisition combined with a dense profile spacing. However, most of the common multi-channel arrays consist of antennae with a spacing of not smaller than 8 cm. The aim of our test survey was to evaluate how an even denser spacing of 4 cm that is provided by the IDS Stream-C GPR device at a centre frequency of 600 MHz can improve the detection of small archaeological features. As a test site, we chose the Great Bath in Kempten-Cambodunum. This first capital of the Roman province Rhaetia never has been overbuilt in the following centuries and even today it is used as a grassland declared as an Archaeological Park. Already in 1911, the Great Bath was excavated and beside the walls of different building phases, also small features like a multitude of hypocaust pillars were unearthed. Hence, this building structure depicts an ideal test site and a 40x40m grid covering the main part of the construction was chosen for the application of the antenna array. As a comparison, the same grid was contemporaneously surveyed with a single antenna IDS Duo device (600 and 200 MHz) and a 50 cm profile spacing. Regarding the walls of the Roman bath, the two surveys show comparable results: both datasets represent the stone constructions of 50 – 90 cm width quite well. Furthermore, the depth slices of both devices provide a differentiation of the single building phases. Nevertheless, the resolution for the multi-channel antenna array is of course much higher due to the denser profile spacing. Huger differences occur for the hypocaust pillars of 25 cm lateral length. These features can be mapped in detail with the 4 cm profile spacing of the IDS Stream-C system. Whereas the IDS Duo can only resolve some of the hypocausts, a multitude of them gets visible between 70 and 110 cm depth in the Stream-C data. As a conclusion, it can be stated that standard archaeological remains like stone walls, for sure, can be surveyed with single antenna GPR devices in a common profile spacing of 50 cm. However, in case of the existence of faint archaeological features the application of ultra-dense antenna arrays like the IDS Stream-C is advisable to get a comprehensive overview of a site without the necessity to excavate them.

Description: poster

Description: The Floodrisk project takes a muti- and interdisciplinary look on the effects of the rise in mine water level in abandoned coal mine regions in Germany. Such effects are heterogeneous ground uplift, stress changes due to the pore pressure changes and the reactivation potential of faults. One of the most directly measurable effects is the induced seismicity. It is known from previous studies that the flooding of old mines can lead to a renewed increase in induced microseismicity in these regions. We focused on the observation of the eastern Ruhr area and investigate in detail the relationship between mine water rise and induced seismicity in the Haus Aden dewatering area. For this purpose, we operate a network of up to 30 short period seismic stations in the region of the former "Bergwerk Ost" colliery, which had the highest seismicity rate in the Ruhr area during active mining. Continuous monitoring of seismicity and mine water levels is available for this region from the active mining phase, through the post-mining phase to flooding. Since the beginning of the flooding, more than 20000 onsets were picked and over 1700 induced events were localised in a magnitude range from -0.7 up to 2.6 MLv. For some larger events, focal mechanisms could be determined. The spatial distribution of hypocentres is divided into two areas, with few events in the central study area and over 95% of earthquakes in its eastern part. Many of these events are spatially clustered and some show quite high waveform similarity. This allows relative localisation to increase the accuracy of the location. Comparing the old galleries,which today serve as the main underground waterways, with the localisations from the relative localisation, strong correlations can be seen. The measured temporal trend of the mine water level, after pumps were shut down in mid-2019, shows a strong correlation with the temporal evolution of the observed micro seismicity. In the first months after the pumps are switched off, the water levels at the observation points rise only slowly and isolated microseismic events occur again. In November 2019, the rise in water levels doubled and at the same time, the strongest induced event in the measurement period was recorded with a magnitude of 2.6 MLv. In the years 2020, 2021 66 and 58 events 〉= MLv 1 were observed, respectively. In contrast to this number only 2- 9 events 〉= MLv 1 per year were observed in the post-minig phase before flooding.

Description: Bundesministerium für Bildung und Forschung

Description: poster

Description: In the last few years, several Unmanned Aerial Vehicle (UAV) based magnetometer systems have been tested for archaeological prospection. Due to their higher sensitivity, scalar magnetometers have been preferred as test sensors. However, total field magnetometer are vulnerable to disturbances, especially those generated by the UAV itself. Therefore, most UAV scalar magnetometer systems use the method of increasing the distance between magnetic sensors and the UAV to reduce interference. But freely suspended sensors tend to swing on ropes under the UAV and can produce data that are strongly influenced by heading errors. For our test, we therefore chose the UAV-fixed, compact setup of the SENSYS MagDrone R4, which is equipped with five three-axis FGM3D/75 fluxgate sensors at 50 cm spacing and is covering a swath width of 2.5 m. The 200Hz sampling rate of the R4 allows easy filtering of interference generated by the UAV and external disturbances like power lines or infrastructure. Magnetograms with a spatial resolution of up to 0.20 m per pixel were produced from the data. At Ganacker, we chose the former infrastructure core of the World War II German Air Force airfield as a test site. A wide range of archaeological structures and features with high magnetic contrast were expected on this area. The test site is currently an open agricultural area with a quite flat terrain. Hence, the R4 could be operated at a fixed flight height of just one metre above the surface of the terrain that is controlled actively by a radar sensor. An area of around 110 hectares were prospected within only four days. The MagDrone R4 system thus offers an outstanding survey area progress that cannot even reached by common vehicle-moved multi-sensor arrays. Here, we present the first results of this test survey by comparing the magnetograms, historical and current geodata. Most of the expected archaeological features and several unknown ones were detected by the R4 system. Our results show that the R4 system is well suited for mapping large archaeological structures with high magnetisations. In the future, we want to compare the R4 data with data from a ground-based fluxgate magnetometer. We also want to test whether the system is suitable for detecting archaeological features that have lower magnetic susceptibility and remanence contrasts with the surrounding soil.

Description: poster

Description: Knowledge sharing in academia has been considered indispensable and is becoming a priority in most European funding schemes. Although we are already quite familiar with the different possibilities to publish our results in open-access journals, open science means way more than that. Open science aims at opening up research processes and granting access to research outputs to researchers, professionals and amateur scientists. There are different ways to ensure the storage and reusability of our data, making it available to other scientists. Furthermore, most of the scientific disciplines migrate their analyses to open-source environments (e.g., R, Phyton). However, tons of code produced remain stored in our personal computers either because we do not know the appropriate tools to share them with our colleagues or because we believe that it is not well structured. In this short course, you would learn how to establish links between publications, data, software and methods. Hence, we will discuss with our experts: i) the options to share our data and code with other peers, ii) obtain some tips to better organize our scripts, and iii) uncover potential barriers to sharing research and discuss possible solutions.

Description: DFG, SUB Göttingen, GFZ Potsdam

Description: presentation

Description: The rise of Open Science practices is impacting the entire scientific publishing culture. The transition to Open Access for text publications goes hand in hand with the growing demand to make scientific data and software available to the general public. The FAIR data principles play a key role in this, designed to make research and the underlying data easily findable, accessible, interoperable, and reusable for humans and machines. Geosciences data are as diverse as their content. They range from large real-time data streams of international observing networks to small data sets produced by individual researchers at their laboratories. Consequently, there need to be different strategies for data management and publication in which research data repositories can be important partners for the researchers. The Specialised Information Service for Geoscience (FID GEO) is a DFG-funded project that is promoting a holistic approach of Open Science that includes scholarly literature, data, samples, and scientific software equally and pushes for their interlinkage. FID GEO has become an important player for connecting researchers, data repositories, information infrastructures, German geoscientific societies, and publishers. FID GEO actively provides data and text publishing services through its affiliated repositories GFZ Data Services and GEO-LEOe-docs, as well as on-demand digitization of printed geoscience literature and maps. In addition, FID GEO aims to inform the German-based geoscience community about all aspects of Open Science and FAIR data by bringing the discussions to the individual disciplines through various communication channels.

Description: DFG, SUB Göttingen, GFZ Potsdam

Description: http://doi.org/10.5194/egusphere-egu22-13354

Description: presentation

Description: The rise of Open Science practices is impacting the entire scientific publishing culture. The transition to Open Access for text publications goes hand in hand with the growing demand to make scientific data and software available to the general public. The FAIR data principles play a key role in this, designed to make research and the underlying data easily findable, accessible, interoperable, and reusable for humans and machines. Geosciences data are as diverse as their content. They range from large real-time data streams of international observing networks to small data sets produced by individual researchers at their laboratories. Consequently, there need to be different strategies for data management and publication in which research data repositories can be important partners for the researchers. The Specialised Information Service for Geoscience (FID GEO) is a DFG-funded project that is promoting a holistic approach of Open Science that includes scholarly literature, data, samples, and scientific software equally and pushes for their interlinkage. FID GEO has become an important player for connecting researchers, data repositories, information infrastructures, German geoscientific societies, and publishers. FID GEO actively provides data and text publishing services through its affiliated repositories GFZ Data Services and GEO-LEOe-docs, as well as on-demand digitization of printed geoscience literature and maps. In addition, FID GEO aims to inform the German-based geoscience community about all aspects of Open Science and FAIR data by bringing the discussions to the individual disciplines through various communication channels.

Description: DFG, SUB Göttingen, GFZ Potsdam

Description: presentation

Description: FZ Data Services is a domain repository for geosciences data comprising the Earth, Space and Environmental Sciences. It assigns digital object identifier (DOI) to data and scientific software since 2004. Hosted at the GFZ German Research Centre for Geosciences (GFZ), the repository has a focus on the curation of long-tail data by domain scientists on one hand, but also provides DOI minting services for several global monitoring networks/observatories in geodesy and geophysics and collaborative projects. Furthermore, as Allocating Agent for the International Generic Sample Number (IGSN), the globally unique persistent identifier for physical samples, GFZ is providing IGSN minting services for physical samples. GFZ Data Services increases the interoperability of long-tail data through (1) the provision of comprehensive domain-specific data description via standardised and machine-readable metadata with controlled domain vocabularies. Metadata is (2) complemeted with comprehensive and standardised technical data descriptions or reports; and (3) by embedding the research data in wider context by providing cross-references through Persistent Identifiers (DOI, IGSN, ORCID, Fundref) to related research products (text, data, software) and people or institutions involved. In addition to the task as a research data publisher, GFZ Data Services is the central node for research data management at the GFZ with information on metadata, data formats, the data publication workflow, FAQ, links to different versions of our metadata editor and downloadable data description templates. Specific data publication guidance is complemented by more general information on data management, like a data management roadmap for PhD students.

Description: GFZ Potsdam

Description: presentation

Description: The shift towards Open Science practices is increasingly demanded by science policy. The transition to Open Access for text publications goes hand in hand with a growing demand to make data, scientific software and samples, freely and FAIRly available to the general public. A persistent problem here is the clear and permanent accessibility and re-usability of scientific publications. This development affects both the scientific publication culture as well as the information infrastructures and poses major challenges to the German- based geosciences community. The specialized information service for geosciences (FID GEO) is a service funded by the German Research Foundation (DFG) and supports the cultural change towards Open Access publications. Hereby, FID GEO pursues a holistic approach to Open Science, including scientific literature, data, samples, and scientific software, and aims to promote their interconnection. FID GEO actively provides data and text publishing services through the affiliated repositories GFZ Data Services and GEO-LEOe-docs, as well as an on-demand digitization service of printed geoscientific literature and maps. The focus here is on the services and information systems that ensure permanently available and reliably citable publications of writings and data. Specifically, the service for text publications is provided in the FID GEOs own subject repository GEO-LEOe- docs. The affiliated research data repository GFZ Data Services is available for the publication of research data and scientific software from the earth and environmental sciences. In addition, FID GEO aims to comprehensively inform the German-based geoscientific community about Open Science and FAIR data by bringing the discussions to the individual disciplines through various communication channels. To strengthen the open information culture in the geosciences, FID GEO collaborates with strategic (inter)national initiatives such as NFDI4Earth, COPDESS and OneGeochemistry.

Description: DFG, SUB Göttingen, GFZ Potsdam

Description: poster

Description: The Berlin Declaration from 2003 was the starting point for Open Access to scholarly publications. Today, however, we speak about Open Science that reaches far beyond Open Access and represents collaborative, transparent and accessible research that includes all kinds of research results: scholarly literature, research data, software, samples, instruments, etc. In addition, efforts such as the FAIR Principles and the Enabling FAIR Data Commitment Statement, combined with increasing demands for machine accessibility to data, have raised user expectations towards the capabilities of research data repositories and datacentres. These repositories are often key partners supporting researchers in fulfilling the new requirements. This presentation will draw the line from major statements and requirements of Open Science, delineate the role of research data repositories as well as major research infrastructures, like the fNFDI4Earth or EPOS (European Plate Observing System) are additional players in making research data accessible in harmonised form.

Description: DFG, SUB Göttingen, GFZ Potsdam

Description: presentation

Description: How open is our research? A discussion on the state of Open Science in the field of mineralogy and petrology. Vortrag bei der Arbeitsgruppe Mineralogie+Petrologie der Uni Potsdam, 29.04.2022 A brief look at current trends toward an open science with a focused look at the state of the art in mineralogy and petrology. Followed by a discussion on how research in mineralogy and petrology can become more open.

Description: DFG, GFZ Potsdam

Description: presentation

Description: We present the first‐generation full‐waveform tomographic model (SinoScope 1.0) for the crust‐mantle structure beneath China and adjacent regions. The three‐component seismograms from 410 earthquakes recorded at 2,427 stations are employed in iterative gradient‐based inversions for three successively broadened period bands of 70–120 s, 50–120 s, and 30–120 s. Synthetic seismograms were computed using GPU‐accelerated spectral‐element simulations of seismic wave propagation in 3‐D anelastic models, and Fréchet derivatives were calculated based on an adjoint‐state method facilitated by a checkpointing algorithm. The inversion involved 352 iterations, which required 18,600 wavefield simulations. SinoScope 1.0 is described in terms of isotropic P‐wave (VP), horizontally and vertically polarized S‐wave velocities (VSH and VSV), and mass density (ρ), which are independently constrained with the same data set coupled with a stochastic L‐BFGS quasi‐Newton optimization scheme. It systematically reduced differences between observed and synthetic full‐length seismograms. We performed a detailed resolution analysis by repairing input random parametric perturbations, indicating that resolution lengths can approach the half propagated wavelength within the well‐covered areas. SinoScope 1.0 reveals strong lateral heterogeneities in the lithosphere, and features correlate well with geological observations, such as sedimentary basins, Holocene volcanoes, Tibetan Plateau, Philippine Sea Plate, and various tectonic units. The asthenosphere lies below the lithosphere beneath East and Southeast Asia, bounded by subduction trenches and cratonic blocks. Furthermore, we observe an enhanced image of well‐known slabs along strongly curved subduction zones, which do not exist in the initial model.

Description: Plain Language Summary: Subduction (a geological process where the oceanic lithosphere descends into the Earth's mantle at convergent boundaries) has been the dominant plate‐tectonic process in the broad Asian region since the Mesozoic (252–66 million years ago). The nature of cold subducting slabs and hot mantle flows can record and affect the tectonic evolution of the overriding lithospheric plates. In this study, we conduct the full waveform inversion on a large data set to image the crust‐mantle structure of this region. The computationally demanding simulations were performed on two of the world's fastest supercomputing facilities. Our new model delivers seismic illumination of the region at unprecedented resolution and exhibits sharper and more detailed shear wave velocity structure in the lithosphere with greatly improved correlations with surface tectonic units compared to previous tomographic models. The narrow low‐velocity layer (generally referred to as asthenosphere) below the lithosphere is present beneath East and Southeast Asia, bounded by subduction trenches and cratonic blocks. The continuous and intense subduction processes are responsible for high‐velocity anomalous bodies in the mantle and the formation of the asthenosphere mentioned above.

Description: Key Points: We construct a new full‐waveform tomographic model of the broad Asian region for 30–120 s period via adjoint and spectral‐element methods. The resolution analysis shows reasonably good resolution in the frequency band of interest and limited trade‐offs between model parameters. Our model shed new light on the subsurface behavior of cold subducting slabs & hot mantle flows and their relation to the overriding plates.

Description: China Scholarship Council http://dx.doi.org/10.13039/501100004543

Description: Swiss National Supercomputing Center

Description: European Unions Horizon 2020

Description: The Collaborative Seismic Earth Model

Description: National Natural Science Foundation of China

Description: Korea Meteorological Administration Research and Development Program

Description: Korean Government

Description: Ministry of Education

Description: https://doi.org/10.5281/zenodo.6597380

Description: Abrupt fluid emissions from shallow marine sediments pose a threat to seafloor installations like wind farms and offshore cables. Quantifying such fluid emissions and linking pockmarks, the seafloor manifestations of fluid escape, to flow in the sub‐seafloor remains notoriously difficult due to an incomplete understanding of the underlying physical processes. Here, using a compositional multi‐phase flow model, we test plausible gas sources for pockmarks in the south‐eastern North Sea, which recent observations suggest have formed in response to major storms. We find that the mobilization of pre‐existing gas pockets is unlikely because free gas, due to its high compressibility, damps the propagation of storm‐induced pressure changes deeper into the subsurface. Rather, our results point to spontaneous appearance of a free gas phase via storm‐induced gas exsolution from pore fluids. This mechanism is primarily driven by the pressure‐sensitivity of gas solubility, and the appearance of free gas is largely confined to sediments in the vicinity of the seafloor. We show that in highly permeable sediments containing gas‐rich pore fluids, wave‐induced pressure changes result in the appearance of a persistent gas phase. This suggests that seafloor fluid escape structures are not always proxies for overpressured shallow gas and that periodic seafloor pressure changes can induce persistent free gas phase to spontaneously appear.

Description: Plain Language Summary: Thousands of pockmarks, circular depressions in the seafloor, were reported in North Sea, presumably formed in response to wave motions during major storms. It has been hypothesized that these pockmarks formed as pre‐existing shallow free‐gas pockets were mobilized by pressure changes of the waves. However, mechanisms that could have mobilized free‐gas are not yet constrained. Moreover, large scale free‐gas accumulations have not been reported in this region, and therefore, commonly invoked mechanisms like tensile failure and breaching of capillary seals are hard to justify as they rely on the presence of pre‐existing gas pockets. Here, through modeling studies, we tackle the question of the source of the observed free‐gas. Our study consists of two parts: First, assuming that some hitherto unknown shallow free‐gas pocket is indeed present, we test whether storm‐induced pressure changes could breach capillary seals. We find that free‐gas damps pressure changes due to its high compressibility, making the mobilization of pre‐existing gas unlikely. In the second part, we propose an alternative mechanism where free‐gas spontaneously appears due to exsolution from pore‐fluids. We test the feasibility of this mechanism and show how periodic pressure changes can lead to a persistent gas phase, that could explain the elusive gas source linked to these pockmarks.

Description: Key Points: Storm‐induced pressure changes can lead to spontaneous appearance of free gas phase near the seafloor. This process is driven by pressure‐sensitive phase instabilities. This mechanism could help explain elusive gas sources in recently observed pockmarks in the North Sea.

Description: Aker BP (AkerBP) http://dx.doi.org/10.13039/100016998

Description: Deutsche Forschungsgemeinschaft, DFG

Description: This GGR Bibliographic Review is a survey of approximately 5200 geoanalytical publications for the year 2021. Selected articles, numbering over 340, containing measurement results for relevant geological and environmental reference materials are listed with individual summaries of target analytes, relevant reference materials and producers. A brief summary of a selection of these publications is included that highlights notable developments in geoanalytical studies, newly developed or characterised RMs, and new datasets of established reference materials that have been re‐analysed using improved or state‐of‐the‐art measurement techniques.

Description: Key Points: Literature review of 5200 geoanalytical publications for the year 2021. 344 selected articles with summaries of target analytes, relevant reference materials and producers. Selected publications include data obtained by new analytical developments and improved analytical protocols for established RMs, and identifies recently developed RMs for specific scientific topics.

Description: http://georem.mpch-mainz.gwdg.de/

Description: Geometric 3D models are a very efficient tool to visualize geological units and structural features that have been presented before just in two dimensions on maps or cross-sections. Most of the information of 3D models is presented as 3D views, virtual wells or horizontal or vertical cross-sections. However, are there further options to transfer as much as possible of the complex information of a 3D model in an adequate way to the user? Is it useful and promising to analyse 3D objects like surfaces or volumes in GIS software? In our investigation we performed a GIS based analysis of an existing geotechnical-geological 3D model of periglacial sediments. The two steps were multiple raster calculations to create geotechnical maps and a digital analysis of surface parameters based on geomorphological techniques and statistics. The investigation area is located in southern Lower Saxony and covers the city of Goettingen and surrounding regions within the valley of the river Leine. The valley is filled by unconsolidated, periglacial sediments of Quaternary age with a variable thickness ranging from 1 to 70 m. The analysed 3D model was constructed with GoCAD in a former project (Nix et al. 2009). The model is based on a heterogenous dataset comprising well data, thematic maps, and outcrop descriptions. Finally, the surfaces and volumes of the following units were modelled, with a special focus on their different geotechnical properties: (1) anthropogenic material, (2) floodplain and slope deposits, (3) freshwater limestone, peat and organic clay, (4) loess, displaced loess, and loess loam, (5) fluvial gravel, (6) outwash fan material, (7) solifluction material, (8) mixed, hetereougenous fillings of subrosion sinks and (9) the surface of the underlying hardrocks. Each top and bottom surface of the Gocad volumes was exported as raster file with additional information stored in an associated attribute table. In ArcGIS various geoprocessing tools were used to calculate and analyse these rasters and to develope thematic geotechnical and geological maps. The geomorphological analysis was subdivided in several steps. Firstly, the surfaces were described visually, concerning their outline, shape and distribution, as well as superficial structures like distinct edges, holes, channels. Secondarily, descriptive statistic parameters of thickness, area and elevation of each surface were calculated. Thirdly, geoprocessing tools of the Spatial Analyst were performed on each surface. Finally, several surfaces were combined to analyse them together, calculating ratios and overlay combinations. Seven thematical geoengineering and geological maps were created, each of them presenting one portion of the three-dimensional dataset: Map of the stratigraphy and depth of the Quaternary base, Map of the thickness of the Quaternary sediments, Distribution map of model units 1 m below ground level, Distribution map of model units 2 m below ground level, Maps of types of different foundation soils, Distribution map of sediments with low loading capacity and Map of distribution and quality of the wells. While the map creation focused on the geotechnical aspects of the model, the applied geomorphological analysis revealed various parameters and values that are related to the geological formation of the model units. Despite the complex dataset represented by the analysed 3D model, thematical information could be transfered into 2D as thematic maps. Some geological characteristics and parameters of the model units were extracted by the descriptive and GIS-based analysis. References Nix, T., Wagner, B., Lange, T. , Fritz, J., Sauter, M. (2009): 3D-Baugrundmodell der quartären Sedimente des Leinetals bei Göttingen. – 17. Tagung für Ingenieurgeologie, S. 223-227, Zittau

Description: poster

Description: The mineralogy, chemical composition, and physical properties of cratonic mantle eclogites with oceanic crustal protoliths can be modified by secondary processes involving interaction with fluids and melts, generated in various slab lithologies upon subduction (auto‐metasomatism) or mantle metasomatism after emplacement into the cratonic lithosphere. Here we combine new and published data to isolate these signatures and evaluate their effects on the chemical and physical properties of eclogite. Mantle metasomatism involving kimberlite‐like, ultramafic carbonated melts (UM carbonated melts) is ubiquitous though not pervasive, and affected between ~20% and 40% of the eclogite population at the various localities investigated here, predominantly at ~60–150 km depth, overlapping cratonic midlithospheric seismic discontinuities. Its hallmarks include lower jadeite component in clinopyroxene and grossular component in garnet, an increase in bulk‐rock MgO ± SiO2, and decrease in FeO and Al2O3 contents, and LREE‐enrichment accompanied by higher Sr, Pb, Th, U, and in part Zr and Nb, as well as lower Li, Cu ± Zn. This is mediated by addition of a high‐temperature pyroxene from a UM carbonated melt, followed by redistribution of this component into garnet and clinopyroxene. As clinopyroxene‐garnet trace‐element distribution coefficients increase with decreasing garnet grossular component, clinopyroxene is the main carrier of the metasomatic signatures. UM carbonated melt‐metasomatism at 〉130–150 km has destroyed the diamond inventory at some localities. These mineralogical and chemical changes contribute to low densities, with implications for eclogite gravitational stability, but negligible changes in shear‐wave velocities, and, if accompanied by H2O‐enrichment, will enhance electrical conductivities compared to unenriched eclogites.

Description: Plain Language Summary: Oceanic crust formed at spreading ridges is recycled in subduction zones and undergoes metamorphism to eclogite. Some of this material is captured in the overlying lithospheric mantle, where it is exhumed by passing magmas. Having formed in spreading ridges, these eclogites have proven invaluable archives for the onset of plate tectonics, for the construction of cratons during subduction/collision, as probes of the convecting mantle from which their precursors formed, and as generators of heterogeneity upon recycling into Earth's convecting mantle. During subduction and until exhumation, interaction with fluids and melts (called metasomatism) can change the mineralogy, chemical composition, and physical properties of mantle eclogites, complicating their interpretation, but a comprehensive study of these effects is lacking so far. We investigated mantle eclogites from ancient continents (cratons) around the globe in order to define hallmarks of metasomatism by subduction‐related fluids and small‐volume ultramafic carbonated mantle melts. We find that the latter is pervasive and occurs predominantly at midlithospheric depths where seismic discontinuities are detected, typically causing diamond destruction and a reduction in density. This has consequences for their gravitational stability and for the interpretation of shearwave velocities in cratons.

Description: Key Points: Exploration of metasomatic effects during subduction of ancient oceanic crust and after its emplacement into cratonic lithospheric mantle. Metasomatism by kimberlite‐like ultramafic melt affected between 20% and 40% of mantle eclogite suites worldwide, mostly at 2–5 GPa. Metasomatism lowers FeO, hence density in eclogite; no significant effect on shearwave velocities.

Description: German Research Foundation http://dx.doi.org/10.13039/501100001659

Description: National Research Foundation (NRF) http://dx.doi.org/10.13039/501100001321

Description: Wilhelm and Else Heraeus Foundation http://dx.doi.org/10.13039/501100011618

Description: Deutsche Forschungsgemeinschaft (DFG, INST

Description: research

Description: Pyrroloquinoline quinone (PQQ) is an important cofactor of calcium‐ and lanthanide‐dependent alcohol dehydrogenases, and has been known for over 30 years. Crystal structures of Ca–MDH enzymes (MDH is methanol dehydrogenase) have been known for some time; however, crystal structures of PQQ with biorelevant metal ions have been lacking in the literature for decades. We report here the first crystal structure analysis of a Ca–PQQ complex outside the protein environment, namely, poly[[undecaaquabis(μ‐4,5‐dioxo‐4,5‐dihydro‐1H‐pyrrolo[2,3‐f]quinoline‐2,7,9‐tricarboxylato)tricalcium(II)] dihydrate], {[Ca3(C14H3N2O8)2(H2O)11]·2H2O}n. The complex crystallized as Ca3PQQ2·13H2O with Ca2+ in three different positions and PQQ3−, including an extensive hydrogen‐bond network. Similarities and differences to the recently reported structure with biorelevant europium (Eu2PQQ2) are discussed.

Description: Pyrroloquinoline quinone (PQQ) is an important cofactor of calcium‐ and lanthanide‐dependent alcohol dehydrogenases. The crystal structure of a Ca–PQQ complex (Ca3PQQ2·13H2O) is reported for the first time outside a protein environment. image

Description: research

Description: In Structural Geology, many projects start with intensive field-based data acquisition campaigns, which might be performed in quite different types of natural or artificial outcrops. For some years, this field work has been substantially influenced and transformed by various close-range sensing techniques that allow the field geologist to create a digital outcrop model (DOM) and to take along plenty of geometrical and spectral information about the outcropping rocks. In general, DOMs can be utilized for outcrop visualization, documentation, manual outcrop analysis (“point-picking”), extraction of spectral data and/or semi-automatic extraction of geometric data. Within a structural investigation DOMs might be deployed for fold analysis, fault analysis, extraction of fracture networks, fracture roughness estimation, detection of neotectonic activities or digitization of geological features for 3D-models of various scales resulting in a large number of analyzing techniques. Latter might be carried out on point clouds or meshes (with or without spectral information) and may differ in pre-processing and processing steps as well as in software solution. Therefore, the analyzing structural geologist faces various tools, data formats, file types, operations and outcomes. Our investigation focus on the compilation of useful, transparent, sustainable and comparable workflows or “pipelines”, which can be executed by open-source/open-access solutions.

Description: poster

Description: Quantitative fold structure analyses at different scales are essential for deducing deformation mechanisms and the reconstruction of the deformation history of orogens. However, not only the field surveying of fold structures, especially in view of their quantification in three dimensions with the classical tools as measuring tape, grid mapping with measuring tapes, geological compass, field book and camera is a time consuming and laborious job, but also the construction of a georeferenced 3D-model of fold structures based on classical data. Another crucial aspect of the classical field surveying of folds is the limitation by poor outcrop conditions. Reasons might be restricted or no accessibility due to high outcrop walls, water or fences, limited visibility because of vegetation, difficult measurability due to very smooth walls or complexity as a result of irregular outcrop walls or distant outcrops. Furthermore, inappropriate oriented outcrop surfaces in respect to the fold geometry can make a survey even worse. Over the past years modern 3D surveying techniques like terrestrial Lidar and digital photogrammetry became progressively affordable for geological field work and now start to complement or replace traditional methods. We started to utilize these techniques on fold structure surveying and to apply quantitative fold structure analysis on different outcrop settings in Central Germany. Different workflows were developed and tested to optimize data conversion, handling and representation. We applied a laser scanner and a single lens reflex camera, complemented by a differential GPS device and laser tachymeter. Data conversion, correction and analysis were done by means of different free as well as commercial software packages. To test different outcrop situations, different quarries, salt mines and steep cliffs, exposing from single fold to complex folds in limestone, greywacke, cherts, rock salt or potassium salt, were selected. As a result, exact 3D point clouds of all exposed folds could be generated by the use of both techniques. The resultant point clouds are suited as excellent visualisation objects as well as base for accurate geometrical measurements in the range of mm or cm of single and complex folds. In addition, the point clouds serve as input dataset for the construction of detailed geological 3D models comprising punctual, linear and plane fold elements. In summary, terrestrial Lidar and digital photogrammetry are excellent field techniques to survey and document exposed folds in the range of few meters to tens of meters, especially under poor outcrop conditions. Different fold sections can now easily be correlated in 3D space to construct complete fold structures with their 3D-fold geometry. Certain fold elements, e.g. axial planes, can be reconstructed much faster and much more accurate compared to the classical approach. The only limiting factor are the very large datasets and the processing power. As next steps, we seek (1) to involve drones for completing data sets from inaccessible areas and perspectives and (2) to incorporate 3D-microfabric analysis data in the fold models as e.g. grain shape an crystallographic preferred orientations to better understand stain paths and deformation mechanisms.

Description: poster

Description: Under certain conditions, ocean surface gravity waves (SGW) interact with the seafloor underneath to trigger relatively faint but measurable seismic waves known as ocean microseisms. Cyclonic storms (e.g. hurricanes, typhoons) wandering over the ocean are major (non-stationary) sources of the former, thus opening the possibility of tracking and studying cyclones by means of their corresponding microseims. For this purpose, we identified storm-related microseisms hidden in the ambient seismic wavefield via array processing. Polarization beamforming, a robust and well-known technique is implemented. The analyses hinge on surface waves (Love and Rayleigh) which, in contrast to P-waves, are stronger but only constrain direction of arrival (without source remoteness). We use a few land-based virtual seismic arrays surrounding the North Atlantic to investigate the signatures of major hurricanes in the microseismic band (0.05-0.16 Hz), in a joint attempt to continuously triangulate their tracks. Our findings show that storm microseisms are intermittently excited with modulated amplitude at localized oceanic regions, particularly over the shallow continental shelves and slopes, having maximum amplitudes virtually independent of storm category. In most cases no detection was possible over deep oceanic regions, nor at distant arrays. Additionally, the rear quadrants and trailing swells of the cyclone provide the optimum SGW spectrum for the generation of microseisms, often shifted more than 500 km off the "eye". As a result of the aforementioned and added to the strong attenuation of storm microseisms, the inversion of tracks or physical properties of storms using a few far-field arrays is discontinuous in most cases, being reliable only if benchmark atmospheric and/or oceanic data is available for comparison. Even if challenging due to the complexity of the coupled phenomena responsible for microseisms, the inversion of site properties, such as bathymetric parameters (e.g. depth, seabed geomorphology), near- bottom geology or SGW spectrum might be possible if storms are treated as natural sources in time-lapse ambient noise investigations. This will likely require near-field (land and underwater) observations using optimal arrays or dense, widespread sensor networks. Improved detection and understanding of ocean microseisms carries a great potential to contribute to mechanically coupled atmosphere-ocean-earth models.

Description: Universität Hamburg

Description: poster

Description: Observations in polar regions show that sea ice deformations are often narrow linear features. These long bands of deformations are referred to as Linear Kinematic Features (LKFs). Viscous‐plastic sea ice models have the capability to simulate LKFs and more generally sea ice deformations. Moreover, viscous‐plastic models simulate a larger number and more refined LKFs as the spatial resolution is increased. Besides grid spacing, other aspects of a numerical implementation, such as the placement of velocities and the associated degrees of freedom, may impact the formation of simulated LKFs. To explore these effects this study compares numerical solutions of sea ice models with different velocity staggering in a benchmark problem. Discretizations based on A‐,B‐, and C‐grid systems on quadrilateral meshes have similar resolution properties as an approximation with an A‐grid staggering on triangular grids (with the same total number of vertices). CD‐grid approximations with a given grid spacing have properties, specifically the number and length of simulated LKFs, that are qualitatively similar to approximations on conventional Arakawa A‐grid, B‐grid, and C‐grid approaches with half the grid spacing or less, making the CD‐discretization more efficient with respect to grid resolution. One reason for this behavior is the fact that the CD‐grid approach has a higher number of degrees of freedom to discretize the velocity field. The higher effective resolution of the CD‐discretization makes it an attractive alternative to conventional discretizations.

Description: Plain Language Summary: Sea ice in the Arctic and Antarctic Oceans plays an important role in the exchange of heat and freshwater between the atmosphere and the ocean and hence in the climate in general. Satellite observations of polar regions show that the ice drift sometimes produces long features that are either cracks (leads) and zones of thicker sea ice (pressure ridges). This phenomenon is called deformation. It is mathematically described by the non‐uniform way in which the ice moves. For numerical models of sea ice motion it is difficult to represent this deformation accurately. Details of the numerics may affect the way these models simulate leads and ridges, their number and length. Specifically, we find by comparing different numerical models, that the way the model variables are ordered on a computational grid to solve the mathematical equations of sea ice motion has an effect of how many deformation features can be represented on a grid with a given spacing between grid points. A new discretization (ordering of model variables) turns out to resolve more details of the approximated field than traditional methods.

Description: Key Points: The placement of the sea ice velocity has a mayor influence on the number of simulated linear kinematic features (LKFs). The CD‐grid resolves twice as many LKFs compared to A, B, C‐grids. A, B, C‐grids on quadrilateral meshes resolve a similar number of LKFs as A‐grids on triangular meshes (with the same total number of nodes).

Description: Porosity is one of the key properties of fluvial sediments. It is defined as the ratio of pore volume to total volume. In river science, porosity is often assumed to be spatially constant, which might be a gross simplification of reality. Ignoring the spatial variations in porosity can cause errors in morphological, ecological, hydrological, hydrogeological and sedimentological applications. Although detailed information about spatial porosity variations can be obtained from porosity measurements at field sites, such information has never been collected where these variations might be important. In this study, field porosity measurements were carried out to quantify the magnitude of the spatial porosity variation for four different sedimentological environments of a braided river: a confluence, a tributary delta, a braid bar and a secondary channel. A nuclear density gauge was used for the measurement of porosity. The nuclear density gauge proved to be a time‐saving and labour‐saving technique that produces accurate porosity values with a root mean square error of 0.03. The four sedimentological environments showed significant differences in porosity, with mean porosity being lower for confluence and bar than for delta and secondary channel. Semi‐variogram analysis showed the absence of any spatial correlation in porosity for distances beyond 4 m. This shows that distance cannot be used as a parameter for porosity extrapolation in a fluvial system unless the extrapolation distance is less than 4 m. At least eight measurements of porosity are required to obtain a reliable estimate of mean porosity in a sedimentary environment, i.e. with uncertainty 〈0.03. Although grain size characteristics were found to have a significant impact on porosity, the relationships between these parameters and porosity were not very strong in this study. The unique porosity dataset, presented in this article, provides a valuable source of information for researchers and river managers.

Description: Deutsche Forschungsgemeinschaft

Description: Mixed sand‐ and gravel‐bed rivers record erosion, transport, and fining signals in their bedload size distributions. Thus, grain‐size data are imperative for studying these processes. However, collecting hundreds to thousands of pebble measurements in steep and dynamic high‐mountain river settings remains challenging. Using the recently published digital grain‐sizing algorithm PebbleCounts, we were able to survey seven large (≥ 1,000 m2) channel cross‐sections and measure thousands to tens‐of‐thousands of grains per survey along a 100‐km stretch of the trunk stream of the Toro Basin in Northwest Argentina. The study region traverses a steep topographic and environmental gradient on the eastern margin of the Central Andean Plateau. Careful counting and validation allows us to identify measurement errors and constrain percentile uncertainties using large sample sizes. In the coarse ≥2.5 cm fraction of bedload, only the uppermost size percentiles (≥95th) vary significantly downstream, whereas the 50th and 84th percentiles show less variability. We note a relation between increases in these upper percentiles and along‐channel junctions with large, steep tributaries. This signal is strongly influenced by lithology and geologic structures, and mixed with local hillslope input. In steep catchments like the Toro Basin, we suggest nonlinear relationships between geomorphic metrics and grain size, whereby the steepest parts of the landscape exert primary control on the upper grain‐size percentiles. Thus, average or median metrics that do not apply weights or thresholds to steeper topography may be less predictive of grain‐size distributions in such settings.

Description: Plain Language Summary: Rock fragments on hillsides are transported to rivers, eventually becoming pebbles, sand, and mud as they are carried downstream by flowing water. The initial size of the pebbles, the way the size changes downstream, and the overprinting of the sizes with new pebbles from other hills and tributaries all form a complex process that can be difficult to disentangle. Yet studying the size of the pebbles at a given stream location or in a sedimentary deposit can provide insights into the conditions of their transport in terms of local upstream patterns of erosion, tectonics, and climate. We show that just looking at the size of the large pebbles on a riverbed can be used to infer the sources of material, but, since there are fewer large pebbles, they require more measurements to quantify. This necessitates new methods for pebble measurement using modern image‐processing tools.

Description: Key Points: Complex grain‐size distributions in dynamic mountain rivers can be computed via thousands of measurements from PebbleCounts. Many measurements allow robust estimation of higher percentiles and we observe the most significant changes in the ≥95th percentile. Downstream grain‐size variation is nonlinearly related to variations in topographic steepness and lithology.

Description: Deutsches Zentrum für Luft‐ und Raumfahrt (DLR) http://dx.doi.org/10.13039/501100002946

Description: Deutsche Forschungsgemeinschaft (DFG) http://dx.doi.org/10.13039/501100001659

Description: Bundesministerium für Bildung und Forschung (BMBF) http://dx.doi.org/10.13039/501100002347

Description: Brandenburger Staatsministerium für Wissenschaft, Forschung und Kultur (MWFK) http://dx.doi.org/10.13039/501100004581

Description: https://zenodo.org/record/5089789

Description: research

Description: A reanalysis is a physically consistent set of optimally merged simulated model states and historical observational data, using data assimilation. High computational costs for modeled processes and assimilation algorithms has led to Earth system specific reanalysis products for the atmosphere, the ocean and the land separately. Recent developments include the advanced uncertainty quantification and the generation of biogeochemical reanalysis for land and ocean. Here, we review atmospheric and oceanic reanalyzes, and more in detail biogeochemical ocean and terrestrial reanalyzes. In particular, we identify land surface, hydrologic and carbon cycle reanalyzes which are nowadays produced in targeted projects for very specific purposes. Although a future joint reanalysis of land surface, hydrologic, and carbon processes represents an analysis of important ecosystem variables, biotic ecosystem variables are assimilated only to a very limited extent. Continuous data sets of ecosystem variables are needed to explore biotic‐abiotic interactions and the response of ecosystems to global change. Based on the review of existing achievements, we identify five major steps required to develop terrestrial ecosystem reanalysis to deliver continuous data streams on ecosystem dynamics.

Description: Plain Language Summary: A reanalysis is a unique set of continuous variables produced by optimally merging a numerical model and observed data. The data are merged with the model using available uncertainty estimates to generate the best possible estimate of the target variables. The framework for generating a reanalysis consists of the model, the data, and the model‐data‐fusion algorithm. The very specific requirements of reanalysis frameworks have led to the development of Earth‐compartment specific reanalysis for the atmosphere, the ocean and land. Here, we review atmospheric and oceanic reanalyzes, and in more detail biogeochemical ocean and terrestrial reanalyzes. In particular, we identify land surface, hydrologic, and carbon cycle reanalyzes which are nowadays produced in targeted projects for very specific purposes. Based on a review of existing achievements, we identify five major steps required to develop reanalysis for terrestrial ecosystem to shed more light on biotic and abiotic interactions. In the future, terrestrial ecosystem reanalysis will deliver continuous data streams on the state and the development of terrestrial ecosystems.

Description: Key Points: Reanalyzes provide decades‐long model‐data‐driven harmonized and continuous data sets for new scientific discoveries. Novel global scale reanalyzes quantify the biogeochemical ocean cycle, terrestrial carbon cycle, land surface, and hydrologic processes. New observation technology and modeling capabilities allow in the near future production of advanced terrestrial ecosystem reanalysis.

Description: European Union's Horizon 2020 research and innovation programme

Description: Deutsche Forschungsgemeinschaft

Description: U.S. Department of Energy

Description: Emory University's Halle Institute for Global Research and the Halle Foundation Collaborative Research

Description: NSF

Description: NASA

Description: Natural Environment Research Council

Description: European Union'’s Horizon 2020 research and innovation programme

Description: NSERC Discovery program, the Ocean Frontier Institute, and MEOPAR

Description: Research Foundation Flanders (FWO)

Description: Helmholtz Association

Description: NASA Terrestrial Ecosystems

Description: Molecular hydrogen (H2) released during serpentinization of mantle rocks is one of the main fuels for chemosynthetic life. Processes of H2 production at slow-spreading mid-ocean ridges (MORs) have received much attention in the past. Less well understood is serpentinization at passive continental margins where different rock types are involved (lherzolite instead of harzburgite/dunite at MORs) and the alteration temperatures tend to be lower (〈200°C vs. 〉200°C). To help closing this knowledge gap we investigated drill core samples from the West Iberia margin. Lherzolitic compositions and spinel geochemistry indicate that the exhumed peridotites resemble sub-continental lithospheric mantle. The rocks are strongly serpentinized, mainly consist of serpentine with little magnetite, and are generally brucite-free. Serpentine can be uncommonly Fe- rich, with XMg Mg/(Mg + Fe) 〈 0.8, and shows distinct compositional trends toward a cronstedtite endmember. Bulk rock and silicate fraction Fe(III)/ Fe ratios are 0.6–0.92 and 0.58–0.8, respectively; our data show that 2/3 of the ferric Fe is accounted for by Fe(III)- serpentine. Mass balance and thermodynamic calculations suggest that the sample’s initial serpentinization produced ∼120 to 〉300 mmol H2 per kg rock. The cold, late-stage weathering of the serpentinites at the seafloor caused additional H2 formation. These results suggest that the H2 generation potential evolves during the transition from continental break-up to ultraslow and, eventually, slow MOR spreading. Metamorphic phase assemblages systematically vary between these settings, which has consequences for H2 yields during serpentinization. At magma-poor rifted margins and ultraslow- spreading MORs, serpentine hosts most Fe(III). Hydrogen yields of 120 to 〉300 mmol and 50–150 mmol H2 per kg rock, respectively, may be expected at temperatures of 〈200°C. At slow-spreading MORs, in contrast, serpentinization may produce 200–350 mmol H2, most of which is related to magnetite formation at 〉200°C. Since, in comparison to slow-spreading MORs, geothermal gradients at magma-poor margins and ultraslow-spreading MORs are lower, larger volumes of low-temperature serpentinite should form in these settings. Serpentinization of lherzolitic rocks at magma-poor margins should produce particularly high amounts of H2 under conditions within the habitable zone. Magma-poor margins may hence be more relevant environments for hydrogenotrophic microbial life than previously thought.

Description: research

Description: https://doi.org/10 .1594/PANGAEA.920208

Description: research

Description: Few data exist that provide insight into pro- cesses affecting the long-term carbon cycle at shallow fore- arc depths. To better understand the mobilization of C in sediments and crust of the subducting slab, we investigated carbonate materials that originate from the subduction chan- nel at the Mariana forearc (〈 20 km) and were recovered during International Ocean Discovery Program Expedition 366. Calcium carbonates occur as vein precipitates within metavolcanic and metasedimentary clasts. The clasts repre- sent portions of the subducting lithosphere, including ocean island basalt, that were altered at lower blueschist facies conditions and were subsequently transported to the fore- arc seafloor by serpentinite mud volcanism. Euhedral arag- onite and calcite and the lack of deformation within the veins suggest carbonate formation in a stress-free environ- ment after peak metamorphism affected their hosts. Inter- growth with barite and marked negative Ce anomalies in car- bonate attest the precipitation within a generally oxic envi- ronment, that is an environment not controlled by serpen- tinization. Strontium and O isotopic compositions in car- bonate (87Sr/86Sr = 0.7052 to 0.7054, δ18OVSMOW = 20 to 24 ‰) imply precipitation from slab-derived fluids at tem- peratures between ∼ 130 and 300 ◦C. These temperature es- timates are consistent with the presence of blueschist facies phases such as lawsonite coexisting with the carbonates in some veins. Incorporated C is inorganic (δ13CVPDB = −1 ‰ to +4 ‰) and likely derived from the decarbonation of cal- careous sediment and/or oceanic crust. These findings pro- vide evidence for the mobilization of C in the downgoing slab at depths of 〈 20 km. Our study shows for the first time in detail that a portion of this C forms carbonate precipitates in the subduction channel of an active convergent margin. This process may be an important asset in understanding the deep carbon cycle since it highlights that some C is lost from the subducting lithosphere before reaching greater depths.

Description: research

Description: Changes in penguin populations on the Antarctic Peninsula have been linked to several environmental factors, but the potentially devastating impact of volcanic activity has not been considered. Here we use detailed biogeochemical analyses to track past penguin colony change over the last 8,500 years on Ardley Island, home to one of the Antarctic Peninsula’s largest breeding populations of gentoo penguins. The first sustained penguin colony was established on Ardley Island c. 6,700 years ago, pre-dating sub-fossil evidence of Peninsula- wide occupation by c. 1,000 years. The colony experienced five population maxima during the Holocene. Overall, we find no consistent relationships with local-regional atmospheric and ocean temperatures or sea-ice conditions, although the colony population maximum, c. 4,000–3,000 years ago, corresponds with regionally elevated temperatures. Instead, at least three of the five phases of penguin colony expansion were abruptly ended by large eruptions from the Deception Island volcano, resulting in near-complete local extinction of the colony, with, on average, 400–800 years required for sustainable recovery.

Description: https://static-content.springer.com/esm/art%3A10.1038%2Fncomms14914/MediaObjects/41467_2017_BFncomms14914_MOESM3080_ESM.xlsx

Description: https://static-content.springer.com/esm/art%3A10.1038%2Fncomms14914/MediaObjects/41467_2017_BFncomms14914_MOESM3081_ESM.xlsx

Description: https://static-content.springer.com/esm/art%3A10.1038%2Fncomms14914/MediaObjects/41467_2017_BFncomms14914_MOESM3082_ESM.xlsx

Description: https://static-content.springer.com/esm/art%3A10.1038%2Fncomms14914/MediaObjects/41467_2017_BFncomms14914_MOESM3083_ESM.xlsx

Description: https://static-content.springer.com/esm/art%3A10.1038%2Fncomms14914/MediaObjects/41467_2017_BFncomms14914_MOESM3084_ESM.xlsx

Description: https://static-content.springer.com/esm/art%3A10.1038%2Fncomms14914/MediaObjects/41467_2017_BFncomms14914_MOESM3085_ESM.xlsx

Description: https://static-content.springer.com/esm/art%3A10.1038%2Fncomms14914/MediaObjects/41467_2017_BFncomms14914_MOESM3086_ESM.xlsx

Description: research

Description: Determination of earth pressures is one of the fundamental tasks in geotechnical engineering. Although many different methods have been utilized to present passive earth pressure coefficients, the influence of non-associated plasticity on the passive earth pressure problem has not been discussed intensively. In this study, finite-element limit analysis and displacement finite-element analysis are applied for frictional materials. Results are compared with selected data from literature in terms of passive earth pressure coefficients, shape of failure mechanism and robustness of the numerical simulation. The results of this study show that passive earth pressure coefficients determined with an associated flow rule are comparable to the Sokolovski solution. However, comparison with a non-associated flow rule reveals that passive earth pressure coefficients are significantly over predicted when following an associated flow rule. Moreover, this study reveals that computational costs for determination of passive earth pressure are considerably larger following a non-associated flow rule. Additionally, the study shows that numerical instabilities arise and failure surfaces become non-unique. It is shown that this problem may be overcome by applying the approach suggested by Davis (Soil Mech 341–354, 1968).

Description: Ruhr-Universität Bochum (1007)

Description: The pure cross-anisotropy is understood as a special scaling of strain (or stress). The scaled tensor is used as an argument in the elastic stiffness (or compliance). Such anisotropy can be overlaid on the top of any elastic stiffness, in particular on one obtained from an elastic potential with its own stress-induced anisotropy. This superposition does not violate the Second Law. The method can be also applied to other functions like plastic potentials or yield surfaces, wherever some cross-anisotropy is desired. The pure cross-anisotropy is described by the sedimentation vector and at most two constants. Scaling with more than two purely anisotropic constants is shown impossible. The formulation was compared with experiments and alternative approaches. Static and dynamic calibration of the pure anisotropy is also discussed. Graphic representation of stiffness with the popular response envelopes requires some enhancement for anisotropy. Several examples are presented. All derivations and examples were accomplished using the algebra program Mathematica.

Description: Karlsruher Institut für Technologie (KIT) (4220)

Description: The near-surface geology of northern Germany is characterized by glacial deposits, deformed by rising Permian and Upper Triassic salt structures. Ground motions potentially associated with salt tectonic processes are very slow and are superimposed by signals of e.g. hydrological and anthropogenic sources. To measure them requires the detection of motion rates in the range of a few millimeters per year with sufficient spatial coverage. For large areas little is known about the rates and the characteristics of ground motions, even though they directly affect anthropogenic infrastructure and could have an impact on the future use of the underground for storage purposes or the exploitation of geothermal energy. To measure ground motion, we use radar interferometric time series data provided by the German Aerospace Center and the Federal Institute for Geosciences and Natural Resources' Ground motion service. These data are based on Synthetic Aperture Radar images acquired by ESA's ERS and Sentinel satellites. Time-series analyses are possible for temporally stable backscattering objects (persistent scatterers) on the ground. Generally, this results in spatially dense observations over built-up areas and sparse observations over rural areas. We use a set of geostatistical methods to analyze these time series data. We detect signals of large-scale surface-deforming processes such as the subsidence of the marshes. We see small-scale signals like the swelling of Permian anhydrite at the Segeberger "Kalkberg". And we can observe subsidence processes over the historic town of Lübeck. Our work extends the area of application of the PS-InSAR technique from areas with high motion rates to regions with particular low motion rates. We discuss methods that can be used to link ERS data to the Sentinel-1 data, in particular, to separate long-term motion processes from short-term effects. We are working on techniques that shall help to decompose different signal sources. Finally, we aim to prepare a set of tools, that can be used by the community.

Description: poster

Description: The objective of the current study is to utilize an innovative method called “change probabilities” for describing fracture roughness. In order to detect and visualize anisotropy of rock joint surfaces, the roughness of one-dimensional profiles taken in different directions is quantified. The central quantifiers, change probabilities, are based on counting monotonic changes in discretizations of a profile. These probabilities, which usually vary with the scale, can be reinterpreted as scale-dependent Hurst exponents. For a large class of Gaussian stochastic processes, change probabilities are shown to be directly related to the classical Hurst exponent, which generalizes a relationship known for fractional Brownian motion. While related to this classical roughness measure, the proposed method is more generally applicable, therefore increasing the flexibility of modeling and investigating surface profiles. In particular, it allows a quick and efficient visualization and detection of roughness anisotropy and scale dependence of roughness.

Description: Karlsruher Institut für Technologie http://dx.doi.org/10.13039/100009133

Description: Bundesministerium für Bildung und Forschung http://dx.doi.org/10.13039/501100002347

Description: After the publication of the type-profiles for the estimation of the joint roughness coefficient (JRC) a discussion evolved about how to adequately use these traces. Based on the chart numerous researchers assembled mathematical correlations with various parameters seeking objectivity in the determination of JRC. Within these works differences concerning the database and the mathematical implementations exist. Consequently, each correlation, although predominantly the same parameters are used, leads to different JRC values. In theory, for any arbitrary profile, irrespective of the particular calculation approach, the same JRC should result. This is a requisite because of the referencing of all correlations to the 10 type-profiles. However, it is shown in this study that in most cases equal or even satisfactorily similar results are not obtained. The discrepancies are vast when non-standard profiles are evaluated, in this case, more than 40,000 traces from six different rock surfaces that cover a broad range of roughness categories. The simple intuitive parameter Z2 served as an agent for the statistical methods because of its broad use and consequently good comparability. On the part of the fractal approaches, three definitions were used. However, JRC inferred from fractal correlations are very much dependent on the particular calculation routine. In fact, the theory of fractals is overly complex for the sparse and low-resolution type-profiles. In summary, fractal approaches do not produce safer or more reliable estimates of roughness compared to simple statistical means and using Z2 perfectly suffices to determine the class of JRC.

Description: Projekt DEAL

Description: In this study, we carried out reactivation potential analysis of discontinuities revealed from four exploration boreholes penetrating heavily faulted and folded Upper Carboniferous rock strata of the Ruhr region. We performed this study based on the notion that slip is controlled by the ratio of shear to effective normal stresses acting on a pre-existing plane of weakness in the prevailing stress field configuration. The results of this analysis were supported by indicators of localized fluid flow, both on micro- and macro-scales, which confirm relationship between secondary permeability and in situ stress state in the Ruhr region. Findings from this study, in conjunction with results of destructive laboratory testing, indicate that the steep NW–SE- and NNE–SSW-striking planar discontinuities are likely to be either close to the critical state or critically stressed in the in situ stress configuration in the Ruhr region. These planar structures, as evidenced by indicators of localized permeability, are the main fluid pathways in the studied region. The NE–SW-striking discontinuities, on the other hand, are most likely to be closed and hydraulically inactive in the prevailing stress state. Based on results gained from this study, implications for utilization of deep geothermal energy in the region were discussed.

Description: Bundesministerium für Bildung und Forschung http://dx.doi.org/10.13039/501100002347

Description: Fraunhofer-Einrichtung für Energieinfrastruktur und Geothermie IEG (1050)

Description: The coexistence of a wide variety of subsurface uses in urban areas requires increasingly demanding geological prediction capacities for characterizing the geological heterogeneities at a small-scale. In particular, detrital systems are characterized by the presence of highly varying sediment mixtures which control the non-constant spatial distribution of properties, therefore presenting a crucial aspect for understanding the small-scale spatial variability of physical properties. The proposed methodology uses the lithological descriptions from drilled boreholes and implements sequential indicator simulation to simulate the cumulative frequencies of each lithological class in the whole sediment mixture. The resulting distributions are expressed by a set of voxel models, referred to as Di models. This solution is able to predict the relative amounts of each grain fraction on a cell-by-cell basis and therefore also derive a virtual grain size distribution. Its implementation allows the modeler to flexibly choose both the grain fractions to be modeled and the precision in the relative quantification. The concept of information entropy is adapted as a measure of the disorder state of the clasts mixture, resulting in the concept of “Model Lithological Uniformity,” proposed as a measure of the degree of detrital homogeneity. Moreover, the “Most Uniform Lithological Model” is presented as a distribution of the most prevailing lithologies. This method was tested in the city of Munich (Germany) using a dataset of over 20,000 boreholes, providing a significant step forward in capturing the spatial heterogeneity of detrital systems and addressing model scenarios for applications requiring variable relative amounts of grain fractions.

Description: Bayerisches Staatsministerium für Umwelt und Verbraucherschutz http://dx.doi.org/10.13039/501100010219

Description: Technische Universität München (1025)

Description: Over the last years, installations of wind turbines (WTs) increased worldwide. Owing to negative effects on humans, WTs are often installed in areas with low population density. Because of low anthropogenic noise, these areas are also well suited for sites of seismological stations. As a consequence, WTs are often installed in the same areas as seismological stations. By comparing the noise in recorded data before and after installation of WTs, seismologists noticed a substantial worsening of station quality leading to conflicts between the operators of WTs and earthquake services. In this study, we compare different techniques to reduce or eliminate the disturbing signal from WTs at seismological stations. For this purpose, we selected a seismological station that shows a significant correlation between the power spectral density and the hourly windspeed measurements. Usually, spectral filtering is used to suppress noise in seismic data processing. However, this approach is not effective when noise and signal have overlapping frequency bands which is the case for WT noise. As a first method, we applied the continuous wavelet transform (CWT) on our data to obtain a time-scale representation. From this representation, we estimated a noise threshold function (Langston & Mousavi, 2019) either from noise before the theoretical P-arrival (pre-noise) or using a noise signal from the past with similar ground velocity conditions at the surrounding WTs. Therefore, we installed low cost seismometers at the surrounding WTs to find similar signals at each WT. From these similar signals, we obtain a noise model at the seismological station, which is used to estimate the threshold function. As a second method, we used a denoising autoencoder (DAE) that learns mapping functions to distinguish between noise and signal (Zhu et al., 2019). In our tests, the threshold function performs well when the event is visible in the raw or spectral filtered data, but it fails when WT noise dominates and the event is hidden. In these cases, the DAE removes the WT noise from the data. However, the DAE must be trained with typical noise samples and high signal-to-noise ratio events to distinguish between signal and interfering noise. Using the threshold function and pre-noise can be applied immediately on real-time data and has a low computational cost. Using a noise model from our prerecorded database at the seismological station does not improve the result and it is more time consuming to find similar ground velocity conditions at the surrounding WTs.

Description: poster

Description: As part of the German continental seismic reflection program (Deutsches Kontinentales Reflexionsseismisches Programm, DEKORP), three large seismic traverses (with the sub-profiles: DEKORP'84-2S and '86-2N; DEKORP'88-9N; DEKORP'90-3A and '90-3B) were measured in the state of Hesse in Germany. The main research topic of DEKORP were deep seismic studies to investigate the lithospheric structure beneath Germany. Thus, for acquisition, strong sources were used to image in these depths, resulting in an excellent S/N ratio, but the main focus was not on the uppermost kilometres. From today's perspective, however, this depth range is of great interest for a wide range of possible technical applications (including medium-deep and deep geothermal projects). The DEKORP profiles cover approx. 450 km in the state of Hesse and mostly cross areas where only insufficient geological data exist (i.e. only few deep boreholes). In order to close or reduce these knowledge gaps, these DEKORP lines were reprocessed in 2019/20. The focus of the reprocessing was on improving the resolution / mapping of geological structures down to a depth of 6 km (approx. 3 s TWT) to describe the prolongation of faults and geological structures in more detail than in previous studies. Nevertheless, deeper structures were also reinterpreted and compared to previous interpretations. The results were directly incorporated into the new geological 3D model of the state of Hesse, developed by the Technical University of Darmstadt (Hessen3D 2.0, BMWi-FKZ: 0325944). In order to achieve these goals and in view of the fact that today's processing methods have improved considerably compared to the 1990‘s, a state-of-the-art reprocessing was applied for all DEKORP profiles traversing the state of Hesse. In comparison to the original processing, additional processing steps like CRS instead of CDP stacking, turning-ray tomography and prestack depth migration were carried out. We present exemplary results of the reprocessing as well as initial geological reinterpretations for the profile DEKORP'88-9N.

Description: poster

Description: The work presented here is part of ongoing studies in the AlpArray initiative and the priority program "MB-4D" regarding the modelling of the lithosphere in the Ligurian Sea (northwestern Mediterranean Sea). It will be based on constraining data from LOBSTER and LISA campaigns of past GEOMAR projects and a study in our research group at CAU Kiel. Our motivation is the combination and interdisciplinary interpretation of independent information from geology, tectonics, geophysics, and petrology. The existing gravity fields, especially the new compilation of the AlpArray Gravity Research Group (AAGRG) is considered as database (high resolution Free Air- and Bouguer anomalies) and the isostatic residual field, besides data of the ICGEM Potsdam (disturbance) and the ESA GOCE gradients for gravity and data for the magnetic field anomaly. The gravity and magnetic fields are analyzed using Euler deconvolution with regularization (R. Pašteka, Comenius University Bratislava) and application of curvature analysis we use both, the fields themselves and their gradients. Besides the calculation of the so-called "3rd derivative" of the gravity potential, we also investigate a possible use of the invariants of the gravity field based on gradient data and compare and correlate the results with structural and tectonic maps in the area of the Ligurian Sea and the adjacent French and Italian mainland. The findings from these comparisons will later be used to initiate the compilation of 3D density and susceptibility models for the studied region.

Description: poster

Description: FloodRisk is an interdisciplinary project focusing on the effects of mine water level rise in bandoned coal mine regions in Germany. Such effects are heterogeneous ground uplift, stress changes due to the change in pore pressure and the reactivation of potential faults. One of the most directly measurable effects is certainly the induced micro seismicity. It is known from previous studies that the flooding of old mines can lead to a renewed increase level in induced micro seismicity in these regions. In this study the relationship between mine water rise, fluid-induced stress changes and induced seismicity in the Haus Aden dewatering area in the eastern Ruhr area (Germany) will be investigated in more detail. For this purpose, we operate a network of currently 21 short period seismic stations in the region of the former "Bergwerk Ost" colliery, which had the highest seismicity rate in the Ruhr area during active underground coal mining. This network is still to be expanded to cover the entire water drainage area, about 30 Raspberry Shake sensors are waiting for the possibility of installation. Nevertheless, the existing network registered almost 1000 induced micro seismic events in a magnitude range from -0.7 up to 2.6 MLv. Many of these events are spatially clustered and some show quite high waveform similarity. This allows relative localisation and can increase the accuracy of the location. The depth location of the earthquakes, within the limits of localisation accuracy, agrees very well with the distribution of seismicity at the time of active mining. The spatial distribution so far seems to be limited by a large inactive transverse fault in the west. It needs to be clarified what influence this fault has on the propagation of mine water in the underground. The measured temporal trend of the mine water level, after pumps were shut down in mid-2019, shows a strong correlation with the temporal evolution of the observed micro seismicity. In the first months after the pumps are switched off, the water levels at the observation points rise only slowly and isolated microseismic events occur again. In November 2019, the rise in water levels doubled and at the same time, the strongest induced event in the measurement period was recorded with a magnitude of 2.6 MLv . In the following months, the seismicity rate ranged from 8 to 34 events above 0.5 MLv per month, some of which were felt.

Description: Bundesministerium für Bildung und Forschung

Description: poster

Description: We introduce an approach for 3D joint interpretation of potential fields and its derivatives under the condition of constraining data and information. The interactive 3D gravity and magnetic application IGMAS (Interactive Gravity and Magnetic Application System) has been around for more than 30 years, initially developed on a mainframe and then transferred to the first DOS PCs, before it was adapted to Linux in the ’90s and finally implemented as a cross-platform Java application with GUI. Since 2019 IGMAS+ is maintained and developed in the Helmholtz Centre Potsdam – GFZ German Research Centre by the staff of Section 4.5 – Basin Modelling and ID2 – eScience Centre. The core of IGMAS+ applies an analytical solution of the volume integral for the gravity and magnetic effect of a homogeneous body. It is based on the reduction of the three-folded integral to an integral over the bounding polyhedrons that are formed by triangles. Later the algorithm has been extended to cover all elements of the gravity tensor as well and the optimized storage enables fast leastsquares inversion of densities and changes to the model geometry and this flexibility makes geometry changes easy. Because of the triangular model structure of model interfaces, IGMAS can handle complex structures (multi- Z surfaces) like the overhangs of salt domes and variable densities due to voxelization. To account for the curvature of the Earth, we use spherical geometries. Therefore IGMAS+ is capable to handle models from big-scale to regional and small-scale models (meters) used in Applied Geophysics.

Description: poster

Description: A new crustal model of the Volga-Uralian subcraton was built. The compilation of the model was subdivided in two steps: (1) inverse gravity modeling followed by (2) thorough forward gravity modeling. For inverse gravity modeling GOCE gravity gradients were used. The effect of the Earth sphericity was taken into account by using tesseroids. Density contrasts between crust and mantle were varied laterally according to the tectonic units present in the region. The model is constrained by the available seismic data including receiver function studies, and deep reflection and refraction profiles. The Moho discontinuity obtained during the gravity inversion was consequently modified, and complemented by the sedimentary cover, upper crust, lower crust, and lithospheric mantle layers in the process of forward gravity modeling. Obtained model showed crustal thickness variation from 34 to more than 55 km in some areas. The thinnest crust with the thickness below 40 km appeared on the Pericaspian basin with the thickest sedimentary column. A relatively thin crust was found along the central Russia rift system, while the thickest crust is located underneath Ural Mountains as well as in the center of the Volga-Uralian subcraton. In both areas the crustal thickness exceeds 50 km. At the same time, the gravity misfit of ca. 95 mGal between the measured Bouguer gravity anomaly and forward calculated gravity field was revealed in the central area of the Volga-Uralian subcraton. This misfit was interpreted and modeled as high-density lower crust which can possibly represent an underplated material. In the end, the new crustal model of Volga-Uralian subcraton respects the gravity and seismic constraints, and reflects the main geological features of the region. This model will be used for further geothermal analysis of the area.

Description: poster

Description: poster

Description: This presentation summarizes input data, procedures and results of probabilistic seismic hazard assessment (PSHA) in Bangladesh in the framework of the project ‘Geo information for Urban Planning and Adaptation to Climate Change’. It is a cooperation of the Geological Survey of Bangladesh (GSB) and the Federal Institute for Geosciences and Natural Resources (BGR) of Germany. The main aim of the project is to provide city planners with “Ground Suitability Maps”, which display different geo-factors. Seismic hazard is one of the geo-factors that contributes to these maps. For the derivation of “Ground Suitability Maps”, the influence of the local underground conditions will be taken into account additionally. A major part of Bangladesh is located in earthquake prone regions due to active tectonics. The Indian plate moves north-eastward towards the Eurasian plate at a velocity of about 6 cm/year. This motion leads to thrusting to the north (Himalaya) and to subduction to the east together with strike-slip mechanism. The thrusting and subduction processes have caused large historical earthquakes even inside Bangladesh (e.g. 1885 Bengal Earthquake M7 and 1918 Srimangal Earthquake M7.6). Therefore, it is crucial to assess seismic hazard in urban planning in Bangladesh. The input databases were compiled from the literature, reviewed and evaluated in this study. These are earthquake catalogs, the distribution of active faults and ground motion prediction equations. The most consistent and reliable databases were selected to be used in PSHA. The data of the earthquake catalog were declustered to eliminate the duplicated events, aftershocks and foreshocks. The spatial distribution of areal seismic sources was characterized using the distributions of earthquakes in the catalog and active faults. The completeness analysis of the earthquake catalog was performed and the Gutenberg-Richter magnitude recurrence distribution was derived for each seismic source. The results of PSHA are presented in the form of peak ground acceleration (PGA) maps with 10% exceedance probability in 50 years. As usual in regional PSHA, the results were compiled assuming bedrock as underground condition (so-called engineering bedrock with shear velocity of Vs30≥760 m/s). The northern and eastern parts of Bangladesh show the highest seismic hazard with PGA around 0.4 g with 10% exceedance probability in 50 years. This observation was expected because of the active tectonics in these parts.

Description: poster

Description: As part of the EPOS-Norway infrastructure project, NORSAR received funding from the Research Council of Norway for a new regional seismic array on Bjørnøya (Bear Island) in the European Arctic. After along planning phase, a six-element broadband array was installed by NORSAR staff in August 2019 and has been providing data to NORSAR in near real-time since then. Due to several logistical and administrative constraints the 6-element array has an aperture of only 300 m. All sites are equipped with Kinemetrics MBB-2 sensors and Earth Data EDR-209 digitizers that are installed in near-surface vaults. Data are automatically copied to the Norwegian node of the European Integrated Data Archive (EIDA) and are openly available. Due to environmental restrictions less than the planned 9 array sites could be installed on Bjørnøya and the non-used instruments are now available to extend the broadband station Hornsund (HSPB), Southern Spitsbergen, to another small aperture broadband array, also with 6 sites. The array installation had to be postponed because of the ongoing pandemic and is now planned for the Arctic summer 2021.

Description: poster

Description: Seismotectonic regions are a basic input in seismic hazard assessment. Several seismotectonic zonations for Germany have been proposed in the past. We have developed a new regionalization based on the definition in the German Nuclear Safety Standard: “A seismotectonic unit is a region for which uniformity is assumed regarding seismic activity, geological structure and development and, in particular, regarding neotectonic conditions”. Our new concept aims for transparent implementation of geological criteria, which we initially analyze separately from seismicity. We strive for a better documentation and justification of the geological elements used to delimit seismotectonic regions, based on an analysis of the geological evolution in six time slices from the Permian (300 Ma) to the Present. The time slices are separated by marked changes in the tectonic regime and associated with the development of new fault systems or reactivation of existing ones. The present-day fault network comprises faults from all time slices. For each time slice, a subset of active faults has been extracted based on geological evidence for fault activity at that time. Uncertainties of these age assignments are documented. The fault subsets delimit regions of different strain intensity. The superposition of strain intensity distributions across all time slices identifies regions affected by polyphase deformation and regions nearly undeformed over geological time, potentially indicating areas of increased or reduced present-day seismic hazard. Our new zonation consists of fewer regions than earlier ones. The geological zonation correlates well with recent seismicity in areas of Cenozoic rifting and reasonably well with less frequent earthquakes in a belt affected by Mesozoic extension and contraction. However, a few stronger earthquake cluster in regions of low geological strain. The most prominent earthquake clusters (Swabian Jura, Vogtland / NW Bohemia) also defy a simple correlation with known geological structures.

Description: poster

Description: poster

Description: As a result of the joint project ‘Subsurface Potentials for Storage and Economic Use in the North German Basin’ (German acronym: TUNB) the Geological Surveys of Northern Germany and the Federal Institute for Geosciences and Natural Resources (BGR) finalized a coherent geological 3D-model of the deep subsurface of the North German Basin in early 2021. The model consists of 13 major base surfaces from Oligocene to Zechstein, fault surfaces and hull surfaces of salt diapirs. In the northwestern part it is based on the datasets of the Tectonic Atlas of NW-Germany (GTA) along with well and seismic data from the hydrocarbon industry. Additionally to modelling the onshore part of Schleswig-Holstein and Hamburg and insuring cross-border consistency to the neighboring federal states, the Geological Survey of Schleswig-Holstein (LLUR) reconstructed a 3D-large-scale velocity model based on previous work from Jaritz et al. (1991). Their velocity approach was developed within the GTA-project and is based on sonic-log and check-shot-velocities. It assumes a linear velocity increase, which is calculated from specific global gradients for different major lithostratigraphic layers and laterally varying starting-velocities. To validate the constructed 3D-velocity-model, its velocities were compared to velocities measured at boreholes by oil and gas companies. In general, a good agreement was found between modeled and measured data (deviation 〈 5%), in which the average velocities seemed to resemble the check-shot data more accurately than the interval velocities the sonic-log-measurements. In distinct locations, the velocity model was used to convert the newly constructed TUNB-horizons from the depth- to the time-domain in order to compare them to seismic sections. Whereas overall a good agreement between horizons and seismic reflectors was found, differences were identified especially in structural complex areas. Whether these can be attributed to earlier interpretations from the GTA, the modelling of the horizons or insufficiencies in the velocity-model has yet to be determined. A follow-up project to the TUNB-project is anticipated to start in early summer 2021. The goal of the project is to derive a consistent velocity model over large parts of the North German Basin. Main challenges will be integrating available borehole and seismic data into existing velocity-modelling approaches with a special focus on establishing cross border consistency to eastern federal states.

Description: poster

Description: About 25% of the Earth’s mid-ocean ridges spread at ultraslow rates of less than 20 mm/yr. However, most of these ultraslow spreading ridges are located in geographically remote areas, which hamper investigation. Consequently, how the crust forms and ages at such spreading centres, which traditional models predict to be magma-starved and cold, remains poorly understood. One of the most accessible ultra-slow spreading centres is the Mid Cayman Spreading Centre (MCSC), in the Caribbean Sea, with spreading rates of ~15-17 mm/yr. CAYSEIS project was proposed to survey the Cayman Trough area in order to obtain new data that constraints the nature of the crust, tectonic structures, lithologies outcropping and hydrothermal processes taking place in this area. Understanding the sub-seabed geophysical structure of the MCSC is key to understanding not only the lithologies and structures exposed at the seabed, but more fundamentally, how they are related at depth and what role hydrothermal fluid flow plays in the geodynamics of ultraslow spreading. CAYSEIS was a joint and multidisciplinary programme of German, British and US American top tier scientists designed for the obtaining of a new high-quality dataset, including 3D Wide-Angle Seismic (WAS), magnetic, gravimetric and seismological data. During the CAYMAN project, we took leverage of the CAYSEIS dataset to invert a 3D tomographic model of the Cayman Trough lithosphere using the Tomo3D code (Meléndez et al., 2015; 2019). This is one of the first times that the Tomo3D code is used for 3D inversion of real datasets. Thus, we are checking our results comparing them with tomographic inversions of 2D lines and testing the different parameters to obtain the more accurate and higher resolution model as possible. The results of this experiment will show not only the lithospheric structure along and across the MSCS, including the exhumed Ocean Core Complexes in the surrounding areas, but the 3D lithospheric configuration of the region which is important to understand the crustal formation processes and the evolution of ultra-slow spreading settings.

Description: Poster

Description: poster

Description: The ongoing activity of Mount Etna and the proximity to the nearby population requires constant monitoring. Infrasound recordings play an important role in volcanic observation because explosive activity near or above ground as well as shallow tremor processes are easier to identify with airborne sound waves than with seismic waves that are significantly scattered and refracted in the volcano edifice. However, infrasound signals are often blurred by noise, in case of Mount Etna, mostly wind induced noise. manual distinction of noisy data from real volcanogenic signals brings along a considerable effort and requires expert knowledge. At Mount Etna five summit craters are currently known with fluctuating levels of activity. This leads to a wide variety of infrasound signal patterns interfered by changing noise levels. In order to distinguish waveforms of noise from signals of volcanic origin we apply unsupervised pattern recognition techniques. We show that by extracting features from the amplitude spectrum different infrasound regimes can be distinguished with Self-Organizing maps (SOMs). This technique provides an option to color-code the results for an intuitive interpretation and allows even for a more detailed recognition of transitional activity regimes. We create a reference data set from multiple months of infrasound waveforms to include as many activity regimes as possible to train the SOM. This enables a fast classification of new data.

Description: poster

Description: DFG FOR 2825

Description: poster

Description: poster

Description: Due to a complex interplay between the Earth and overlying ice sheets, a large variety of subglacial landforms developed. One example is the in the North German Basin widely spread phenomenon of tunnel valleys. An observed correlation to underlying salt structures is often explained mechanically. We focus on an alternative hypothesis based on thermodynamic processes: As salt better conducts heat than the surrounding rocks, the geothermal heat flux is augmented above salt structures. This leads to melting processes at the interface between the Earth and the ice sheet. The subglacial rivers finally erode the tunnel valleys. To test this hypothesis, we model related hydrothermal processes by means of a finite-difference open-source code (SHEMAT-Suit). The model accounts for heat conduction, groundwater flows, processes in the glaciothermal system such as the motion and spatiotemporal temperature evolution within the ice, and finally the coupling of both at the subglacial interface to account for the feedback mechanisms. Glaciothermal system and coupling processes are incorporated based on an idealized 1D model for the ice cover. We present a scaling analysis to discuss dominant processes. Our results show that a purely conductive subsurface (complete absence of groundwater flow) leads to a very moderate increase of the geothermal heat flux above salt structures. This implies a slight increase of the melting rates, which by itself is not enough to trigger tunnel valley erosion. Additional hydrothermal flows e.g. through fault zones may increase the subglacial melting rates. In this contribution, we will present results from a case study in the Southern North Sea. A 2D seismic section includes two tunnel valleys above salt structures. To model the state prior to erosion and sedimentation during and after the Quaternary glaciations, the Quaternary strata is replaced by strata with the same physical properties and thicknesses than the Paleocene to Miocene strata. Simulation runs with SHEMAT-Suite calculated the subsurface temperature distribution and the geothermal heat flux distribution at the subglacial interface. This allows assessing the subglacial melting rates along with the temperature profile within the ice cover for a number of glaciation scenarios. Current results show that thermodynamic processes reinforce the formation of tunnel valleys together with e.g. mechanical weakening by faulting.

Description: poster

Description: A combination of noninvasive geophysical magnetic gradiometery and electrical resistivity tomography (ERT) was employed to locate the remains of an old church form the 9th century in Neuss-Norf, Germany. The gradiometery survey was carried out along 27 parallel profiles oriented approximately E–W with a profile spacing of 1 m and a sampling spatial interval of 0.5 m along each profile in order to detect the distributions of the possible buried walls of the church and other related archaeological remains. The vertical gradient of the magnetic field, with a fixed distance of 1.04 m between the sensors, was measured. The lower sensor was fixed at a height of 0.32 m from the ground surface. The magnetic data were transferred to the frequency domain using Fast Fourier transform then reduced to the magnetic pole. The analytic signal and power spectrum techniques were applied to the obtained magnetic data. Moreover, ERT measurements were performed based on the results of the magnetic survey along 12 profiles utilizing the Wenner and Dipole-Dipole arrays with 0.5 m electrode spacing. The ERT data from both arrays were merged into one dataset to form a non-conventional mixed array. The ERT data were inverted into 2D resistivity models using robust (blocky) inversion technique, and then a 3D resistivity prospective was created. The combined interpretation of the magnetic and ERT showed that the archaeological structures are close to the ground surface with a maximum depth of up to 2 m. We successfully detected anomalous zones that could be associated with the walls of at least one ancient church-building in addition to several possible archaeological structures in the survey area. A considerable agreement between the results of both methods was observed. Highly magnetic sources that could be associated with metallic objects within tomb-like structures were detected. An archaeological map of the possible location of the old church and the assumed surrounding tombs and features was constructed. Finally, some promising places were suggested in order to start an archaeological excavation in the site based on the findings of our research.

Description: poster

Description: The German continental seismic reflection program DEKORP (DEutsches KOntinentales Reflexionsseismisches Programm) was carried out in the years between 1984 and 1999. The aim of DEKORP was to investigate the deep crustal structure of Germany with high-resolution near-vertical incidence (mostly vibro)seismic acquisition, supplemented by wide-angle seismic and other target-oriented piggy-back experiments, all complemented by optimized methods of data processing and interpretation. The DEKORP project was an equivalent to many other deep-seismic programs world-wide such as COCORP, BIRPS, LITHOPROBE, ECORS, CROP, BELCORP, IBERSEIS and others. The resulting DEKORP database consists of approximately 40 crustal-scale 2D-seismic reflection lines covering a total of ca. 4 700 km and one 3D-seismic survey covering ca. 400 km², recorded in close connection with the German Continental Deep Drilling Program (KTB). Nowadays, re-recording of these seismic traverses in the same extent and quality would often not be possible anymore due to increased acquisition costs and tightened permission requirements. Therefore these datasets provide unique and deep insights into the subsurface below Germany covering the earth’s crust from the surface to the upper mantle. Currently, many of the original raw data are still stored on old storage media and in formats, which can only be read by special devices, programs and experts. To prevent the final loss of this valuable geoscientific treasure an initiative at GFZ transcripts all relevant DEKORP data to modern formats and media. Over the last few years the demand for DEKORP data continuously increased. Several academic institutions and commercial companies reprocess and/or reinterpret these data, which lead to significant improvements in the quality of the results. Fields of applications are geothermal development, hazard analysis, hydrocarbon/shale gas exploration, underground gas storage, tunnel construction, disposal of nuclear waste and more. To simplify the data access for the scientific as well as for the commercial geo-community, a well-structured provision and utilisation concept is being developed. The concept includes so-called data publications with DOIs, a defined license model and automised retrieval for each of the surveys providing raw data, processed data, meta data, related links and more. The plan aims to have all relevant DEKORP datasets compiled and prepared for access via web interface till 2022.

Description: poster

Description: NEXD is an open source software package for the simulation of seismic waves in complex geological media. This includes elastic, viscoelastic, porous and fractured media with complex geometries. For the computation of the wave fields, the nodal discontinuous Galerkin approach (NDG) is used. The NDG approach combines unstructured tetrahedral meshes with an element-wise, high-order spatial interpolation of the wave field based on Lagrange polynomials. NEXD offers capabilities for modeling wave propagation in one-, two- and three-dimensional settings of very different spatial scale with little logistical overhead. It allows the import of external triangular (2D) and tetrahedral (3D) meshes provided by independent meshing software and can be run in a parallel computing environment. The computation of adjoint wavefields and an interface for the computation of waveform sensitivity kernels are offered. The method is verified by means of symmetry tests and the method of exact solutions. The capabilities of NEXD are demonstrated through, for example, a 2D synthetic survey of a geological carbon storage site. The most recent developments have been the inclusion of porous media in 2D and the inversion capabilities to the latest release versions of the 2D and 3D codes as well as the release of the 1D code. NEXD is available on GitHub: https://github.com/seismology-RUB.

Description: poster

Description: The ocean worlds of our Solar System, like Saturn's moon Enceladus and Jupiter's moon Europa are covered with ice. Recently, these icy moons gained further scientific interest, as they are attributed some potential to sustain or host extraterrestrial life in a subglacial ocean. The investigation of these moons will also help to understand the evolution of the Solar System. The in-situ exploration of these moons requires novel technological solutions as well as intelligent data acquisition and interpretation tools. In 2020, the DLR Space Administration started the TRIPLE project (Technologies for Rapid Ice Penetration and subglacial Lake Exploration) which develops an integrated concept for a melting probe that launches an autonomous underwater vehicle (nanoAUV) into a scientifically interesting water reservoir and an AstroBioLab for in-situ analysis. These three components build up the TRIPLE system. As part of a second project stage, it is envisioned to build the TRIPLE system and test it in Antarctica in 2026. In this contribution, we are going to present the general concept of TRIPLE with a focus on the geophysically most relevant aspects. To navigate the melting probe through the ice, a forefield reconnaissance system (TRIPLE-FRS) based on combined radar and sonar techniques is designed. This will include radar antennas directly integrated into the melting head combined with a pulse amplifier and a piezoelectric acoustic transducer just behind the melting head. In addition, an in-situ permittivity sensor will be implemented to account for the ice structure dependent propagation speed of electromagnetic waves. With this system, obstacles as well as the ice-water interface at the bottom of the icy shell could be detected. To deliver key parameters such as transit time and overall energy requirement, a virtual test bed for strategic mission planning is currently under development. This consists of the Ice Data Hub that combines available data from Earth or any other planetary body – measured or taken from the literature – and allows display, interpretation and export of data, as well as trajectory models for the melting probe. We develop high-fidelity thermal contact models for the phase change as well as macroscopic trajectory models that consider the thermodynamic melting process and the convective loss of heat via the melt-water flow.

Description: Bundesministerium für Wirtschaft und Energie (FKZ: 50NA1908, 50RK2050, 50RK2051, 50RK2052, 50RK2053)

Description: poster

Description: Deutsche Forschungsgemeinschaft

Description: poster

Description: Modelling the propagation of seismic waves in porous media gets more and more popular in the seismological community since it is an important but challenging task in the field of computational seismology. The fluid content of, for example, reservoir rocks or soils, and the interaction between the fluid and the rock or between different immiscible fluids has to be taken into account to accurately describe seismic wave propagation through such porous media. Often, numerical models are based on the elastic wave equation and some might include artificially introduced attenuation. This simplifies the problem but only approximates the true physics involved. Hence, the results are also simplified and could lack accuracy or miss phenomena in some applications. The aim of the conducted work was the consistent derivation of a theory for seismic wave propagation in porous media saturated by two immiscible fluids and the accompanying numerical solution for the derived wave equation. The theory is based on Biot's theory of poroelasticity. Starting from the basic conservation equations (energy, momentum, etc.) and generally accepted laws, the theory was derived using a macroscopic approach which demands that the wavelength is significantly larger than the size of the heterogeneities in the medium due to the size of the grains and pores or due to effects on the mesoscopic scale. This condition is usually fulfilled for seismic waves since the typical wavelength of seismic waves is in the order of 10 m to 10 km. Fluid flow is described by a Darcy type flow law and interactions between the fluids by means of capillary pressure curve models. In addition, consistent boundary conditions on interfaces between poroelastic media and elastic or acoustic media are derived from this poroelastic theory itself. The nodal discontinuous Galerkin method is used for the numerical modelling. The poroelastic solver is integrated into the 1D and 2D codes of the larger software package NEXD that uses the nodal discontinuous Galerkin method to solve wave equations. The implementation has been verified using symmetry tests and the method of exact solutions. This work has potential for applications in various scientific fields like, for example, exploration and monitoring of hydrocarbon or geothermal reservoirs as well as CO2 storage sites.

Description: poster

Description: i-EDDA-Research School (BGR Spandau); i-EDDA = Innovative Exploration Drilling and Data Acquisition, https://www.iedda.eu/ ; Datenmanagement Workshop: 05.11.2021

Description: DFG, GFZ Potsdam

Description: presentation

Description: We present context-aware benchmarking and performance engineering of a mature TByte-scale air quality database system which was created by the Tropospheric Ozone Assessment Report (TOAR) and contains one of the world’s largest collections of near-surface air quality measurements. A special feature of our data service https://join.fz-juelich.de is on-demand processing of several air quality metrics directly from the TOAR database. As a service that is used by more than 350 users of the international air quality research community, our web service must be easily accessible and functionally flexible, while delivering good performance. The current on-demand calculations of air quality metrics outside the database together with the necessary transfer of large volume raw data are identified as the major performance bottleneck. In this study, we therefore explore and benchmark in-database approaches for the statistical processing, which results in performance enhancements of up to 32%.

Description: H2020 European Research Council https://doi.org/10.13039/100010663

Description: H2020 European Research Council https://doi.org/10.13039/100010663

Description: H2020 European Research Council https://doi.org/10.13039/100010663

Description: Forschungszentrum Jülich GmbH (4205)

Description: https://join.fz-juelich.de/

Description: Increasing ice loss of the Antarctic Ice Sheet (AIS) due to global climate change affects the orientation of the Earth’s spin axis with respect to an Earth-fixed reference system (polar motion). Here the contribution of the decreasing AIS to the excitation of polar motion is quantified from precise time variable gravity field observations of the Gravity Recovery and Climate Experiment (GRACE) and from measurements of the changing ice sheet elevation from altimeter satellites. While the GRACE gravity field models need to be reduced by noise and leakage effects from neighboring subsystems, the ice volume changes observed by satellite altimetry have to be converted into ice mass changes. In this study we investigate how much individual gravimetry and altimetry solutions differ from each other. We show that due to combination of individual solutions systematic and random errors of the data processing can be reduced and the robustness of the geodetic derived AIS polar motion excitations can be increased. We investigate the interannual variability of the Antarctic polar motion excitation functions by means of piecewise linear trends. We find that the long-term behavior of the three ice sheet subregions: EAIS (East Antarctic Ice Sheet), WAIS (West Antarctic Ice Sheet) and APIS (Antarctic Peninsula Ice Sheet) is quite different. While APIS polar motion excitations show no significant interannual variations during the study period 2003-2015, the trend of the WAIS and EAIS polar motion excitations increased in 2006 and again in 2009 while it started slightly to decline in 2013. AIS mass changes explain about 45% of the observed magnitude of the polar motion vector (excluding glacial isosatic adjustment). They cause the pole position vector to drift along 59〈sup〉◦〈/sup〉 East longitude with an amplitude of 2.7 mas/yr. Thus the contribution of the AIS has to be considered to close the budget of the geophysical excitation functions of polar motion.

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: European Space Agency (FR)

Description: Technische Universität München (1025)

Description: Even though micropolar theories are widely applied for engineering applications such as the design of metamaterials, applications in the study of the Earth’s interior still remain limited and in particular in seismology. This is due to the lack of understanding of the required elastic material parameters present in the theory as well as the eigenfrequency ωr which is not observed in seismic data. By showing that the general dynamic equations of the Timoshenko’s beam is a particular case of the micropolar theory we are able to connect micropolar elastic parameters to physically measurable quantities. We then present an alternative micropolar model that, based on the same physical basis as the original model, circumvents the problem of the original eigenfrequency ωr laking in seismological data. We finally validate our model with a seismic experiment and show it is relevant to explain observed seismic dispersion curves.

Description: Westfälische Wilhelms-Universität Münster (1056)

Description: On the Qinghai-Tibetan Plateau, the ecological crisis appears in every physical and social aspect, and typically affects the rangelands on which most local people depend for their livelihood. Until today, however, changes in land use and vegetation cover on this Plateau have not yet been adequately and practically treated in studying the framework of sustainability. This thesis deals first with the basic ecological influences and interactions before reviewing their implications for development. The limitations of an ecological perspective is recognized, but this stems mainly from the misapplication of perspective and should not detract from its value. Ecology is undoubtedly the most basic and pervasive of the many determinants of the pastoral production systems of the high-frigid land. To ignore ecology or to fail to take it into account is to court disaster. The amount of literature reflects the importance of the subject but is not reviewed here in any detail; rather the aim is simply to set the scene for later discussions.

Description: thesis

Description: DFG, SUB Göttingen

Description: Unter Berücksichtigung der Forschungssituation und der Resultate der wissenschaftstheoretischen Vororientierung stellt sich diese Arbeit zum einen die Aufgabe, die bei SOULAVIE mit den Bezeichnungen "géographie physique des plantes“, "géographie physique des animaux" und "geographie physique de l'homme et de la femme" verbundenen Grundannahmen, Fragestellungen, Verfahrensweisen und Ergebnisse aus seinen naturhistorischen Schriften zu erschließen, sie zu systematisieren und als Gedankensystem zu rekonstruieren. Zum anderen fragt sie nach dem Spannungsverhältnis von Konventionalität und Neuartigkeit, in dem diese Forschungsansätze SOULAVIES zu den naturwissenschaftlichen Bemühungen seiner Epoche stehen. Da die wesentlichen Aussagen eines einzelnen Forschers nur aus ihrer Einbettung in die Wissenschaft ihrer Zeit verständlich werden, erfordert auch die Bearbeitung der erstgenannten Problemstellung einen Rekurs auf den zeitgenössischen wissenschaftlichen Hintergrund. Es versteht sich von selbst, dass angesichts dieses doppelten Rückbezuges auch Licht auf die ideengeschichtlichen Wurzeln der geographischen Konzeptionen SOULAVIES fällt.

Description: thesis

Description: DFG, SUB Göttingen

Description: Der vorliegende Band 52 der Abhandlungen „Anthropogeographie“ des Instituts für Geographische Wissenschaften der Freien Universität Berlin erfasst mit 54 Beiträgen den wissenschaftlichen Forschungsstand der IGU-Kommission “Urban Development and Urban Life”. Die vorgelegten Fassungen stellen die diskutierte und revidierte Form der Beiträge dar, die während der Jahrestagung der Kommission vom 15. bis 20. August 1994 in Berlin gehalten wurden. Die Tagung selbst stand unter dem Thema “Urban Transition and Quality of Life” und forderte damit in den jeweils zwei Einführungs- und zahlreichen weiterführenden Referaten der drei Arbeitsthemen 1) “Restructuring Urban Systems in Central and Eastern Europe“, 2) “Unemployment and Informal Economy in Cities” und 3) “Urban Identity, City Image and Urban Marketing” zu Bewertungen des Forschungsstandes bzw. kritischen Stellungnahmen durch aktuelle Forschungen heraus.

Description: 〈html〉 〈body〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0004.pdf"〉Titelseite〈/a〉〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0005.pdf"〉Titelseite〈/a〉〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0006.pdf"〉Vorwort〈/a〉〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0007.pdf"〉Foreword〈/a〉〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0008.pdf"〉Inhaltsverzeichnis〈/a〉〈/li〉 〈/body〉 〈/html〉 〈br〉 〈html〉 〈body〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0009.pdf"〉I. Introduction〈/a〉〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0010.pdf"〉Reflections Between Urban and University Structures and Planning〈/a〉〈br〉(Dietz, K., Braun, G. O.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0011.pdf"〉Berlin - Metropolis in Transition: The State of Urban Development in 1994〈/a〉〈br〉(Ellger, C.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0014.pdf"〉Planning for the Future of Berlin〈/a〉〈br〉(Branoner, W.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0016.pdf"〉Strategic Planning in Berlin〈/a〉〈br〉(Braun, G. O.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0019.pdf"〉Change in the Urban Structure in the Five New Länder: The Example of Eisenach〈/a〉〈br〉(Cassel, M., Schwaderer, G.)〈/li〉 〈/body〉 〈/html〉 〈br〉 〈html〉 〈body〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0021.pdf"〉II. Restructuring of Inter- and Intra-Urban Systems〈/a〉〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0022.pdf"〉Questions on Hierarchical Change in Urban Systems in Eastern and Central Europe〈/a〉〈br〉(Rey, V.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0026.pdf"〉Urban and Urban System Development in Central and Eastern Europe〈/a〉〈br〉(Grimm, F.-D., Taege, G.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0030.pdf"〉The Impact of Recent Socio-Political and Economic Changes on the Development of Slovenian Towns〈/a〉〈br〉(Pak, M.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0033.pdf"〉An Urban Evolution Model Applied to Romania's Towns〈/a〉〈br〉(Ianos, I.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0038.pdf"〉Warsaw: Development Problems and Strategies〈/a〉〈br〉(Korcelli, P., Potrykowska, A., Weclawowicz, G.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0039.pdf"〉The Urban System and Emerging Structure: An Application of Gibb's Measure to the Case of India〈/a〉〈br〉(Mookherjee, D.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0041.pdf"〉Reorganisation of the Urban System in Sri Lanka〈/a〉〈br〉(Wanasinghe, D. S.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0046.pdf"〉The Formation of an European Urban System〈/a〉〈br〉(Rozenblat, C., Pumain, D.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0051.pdf"〉On the Impact of Internationalization Process on the Finnish Settlement System〈/a〉〈br〉(Palomäki, M.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0055.pdf"〉Global Perspective of Seoul as a World City in the Region of Pacific Rims〈/a〉〈br〉(Kim, I.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0057.pdf"〉Urban Systems Dynamics: Evidence for the Toronto Urban System: 1930-1991〈/a〉〈br〉(Preston, R. E.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0060.pdf"〉The Growth of Toronto: A Market-Share Approach〈/a〉〈br〉(Simmons, J.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0064.pdf"〉Industrial Restructuring and Urban Development: An Examination in Metropolitan Detroit〈/a〉〈br〉(Sinclair, R.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0067.pdf"〉Counter-Urbanization in Perspective: Images and Reality in Settlement System Change〈/a〉〈br〉(Davies, W.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0069.pdf"〉Urban Region as a Geoeconomic and Geopolitical Unit in "A Europe of Regions"〈/a〉〈br〉(Vartiainen, P.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0070.pdf"〉The Ethnic and Social Division of a World City: The Case of London〈/a〉〈br〉(Petsimeris, P.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0075.pdf"〉Transformation in the South African Apartheid City – Residential Desegregation in Pietersburg: A Case Study〈/a〉〈br〉(Donaldson, S. E., Kotze, N. J.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0078.pdf"〉Intra-Urban Distribution of Ethnic Minorities in Munich and Duisburg〈/a〉〈br〉(Yamamoto, K.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0083.pdf"〉New Urban Poverty in the Basque Country and Navarra〈/a〉〈br〉(Ferrer, M., d’Entremont, A., Ciscar, I., Saracibar, M. A.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0085.pdf"〉Polarization and the Dutch Welfare State. The Case of Amsterdam〈/a〉〈br〉(Ostendorf, W., Musterd, S.)〈/li〉 〈/body〉 〈/html〉 〈br〉 〈html〉 〈body〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0088.pdf"〉III. Informal Economy〈/a〉〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0089.pdf"〉The Labour Market and the Role of Informality in Urban Mexico〈/a〉〈br〉(Aguilar, A. G.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0092.pdf"〉The Informal Economy in Chinese Cities〈/a〉〈br〉(Qi, S. D.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0093.pdf"〉The Urban Informal Sector: A Third World Experience〈/a〉〈br〉(Geeta Reddy, A., Rajkumar, G.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0094.pdf"〉Informal Urban Economy in the Lake Chad Basin〈/a〉〈br〉(Simeu-Kamdem, M.)〈/li〉 〈/body〉 〈/html〉 〈br〉 〈html〉 〈body〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0095.pdf"〉IV. Recent Urban Processes〈/a〉〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0096.pdf"〉City Challenge: Involving the Community in UK Urban Policy? 〈/a〉〈br〉(Lewis, J. R.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0098.pdf"〉Intra-Metropolitan Relocation of Work Places: The Case of Stockholm〈/a〉〈br〉(Mahieu, R.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0102.pdf"〉Problems of Mobility and Traffic in the Ruhr〈/a〉〈br〉(Lötscher, L., Fleisgarten, S.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0103.pdf"〉Moscow 1993: The Emergence of Housing Market and the Problems of Intra-Urban Segregation〈/a〉〈br〉(Vendina, O. I.)〈/li〉 〈/body〉 〈/html〉 〈br〉 〈html〉 〈body〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0107.pdf"〉V. Urban Marketing, City Management, Urban Identity〈/a〉〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0108.pdf"〉Urban Marketing: A Review〈/a〉〈br〉(Borchert, J. G.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0109.pdf"〉Urban Identity, City Image and Urban Marketing〈/a〉〈br〉(Dematteis, G.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0110.pdf"〉The Hidden Faces of the City〈/a〉〈br〉(Racine, J.-B.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0111.pdf"〉The City of Bergen. Image and Marketing〈/a〉〈br〉(Sjøholt, P.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0113.pdf"〉Structural Underpinnings of an Olympic City Bid〈/a〉〈br〉(Davies, R. J.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0116.pdf"〉Double Cities: Identity and Marketing of a New Urban Product〈/a〉〈br〉(Buursink, J.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0118.pdf"〉Corporate Space and Emerging Spatial Order in Japan〈/a〉〈br〉(Fujita, N.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0125.pdf"〉The Development of International Tourism in Tallinn, Estonia〈/a〉〈br〉(Kuus, M.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0126.pdf"〉Conflict, Consent, Cooperation: Comprehensive Planning in Germany Beyond Market and State〈/a〉〈br〉(Helbrecht, I.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0127.pdf"〉Urban Strategies: Mega Events. A Copenhagen Perspective〈/a〉〈br〉(Thor Andersen, H., Wichmann Matthiessen, C.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0129.pdf"〉From Spaces of Planning to Places of Resistance〈/a〉〈br〉(Haarni, T.)〈/li〉 〈/body〉 〈/html〉 〈br〉 〈html〉 〈body〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0132.pdf"〉VI. Re-Cycling Urban Landscape〈/a〉〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0133.pdf"〉The Role of Gentrification in the Changing Ecology of Income: Evidence From Canadian Cities and Implications for Further Research〈/a〉〈br〉(Bourne, L. S.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0134.pdf"〉Gentrification and the Youth Movements of the 1960s〈/a〉〈br〉(Ley, D.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0135.pdf"〉Recycling Urban Landscapes - Beyond the Power〈/a〉〈br〉(Andersson, H.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0136.pdf"〉Why was the Central City Population Stabilized? The Case of Copenhagen〈/a〉〈br〉(Illeris, S.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0138.pdf"〉Toronto's Underground City: Excavating the Terms of Access〈/a〉〈br〉(Hopkins, J.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0140.pdf"〉The Interpretation of Sea in Urban Planning and Everyday Life: The Case of Helsinki〈/a〉〈br〉(Karvinen, M.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0141.pdf"〉Urban Destinies - What are the Trends?〈/a〉〈br〉(Wärneryd, O.)〈/li〉 〈/body〉 〈/html〉 〈br〉 〈html〉 〈body〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0142.pdf"〉VII. Urban Modelling〈/a〉〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0143.pdf"〉DUVA – A Concept for a Metadata Driven Statistical Production and Information System〈/a〉〈br〉(Kopp, N.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0145.pdf"〉A New Era for Urban Modelling?〈/a〉〈br〉(Pumain, D.)〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0146.pdf"〉Competing Order Parameters in a Self-Organizing City〈/a〉〈br〉(Portugali, J., Benenson, I.)〈/li〉 〈/body〉 〈/html〉 〈br〉 〈html〉 〈body〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0150.pdf"〉Authors and Addresses〈/a〉〈/li〉 〈li〉〈a href="https://gdz.sub.uni-goettingen.de/download/pdf/PPN1030505985/LOG_0151.pdf"〉[Werbung]〈/a〉〈/li〉 〈/body〉 〈/html〉

Description: conference

Description: DFG, SUB Göttingen

Description: Since the 14th century, moderate seismic activity with 14 earthquakes of magnitude MW≥5.0 occurred in Western Europe in a region extending from the Lower Rhine Graben (LRG) to the southern North Sea. In this paper, we investigate how well this seismic activity could reflect that of the future. The observed earthquake activity in the LRG is continuous and concentrates on the Quaternary normal faults delimiting the LRG, which are also the source of large surface rupturing Holocene and Late Pleistocene earthquakes. The estimated magnitudes of these past earthquakes range from 6.3±0.3 to 7.0±0.3 while their average recurrence on individual faults varies from ten thousand to a few ten thousand years, which makes foreseeing future activity over the long-term possible. Three of the largest historical earthquakes with MW≥5.5 occurred outside the LRG. Late Quaternary activity along the fault zones suspected to be the source of two of these earthquakes, i.e. the 1580 Strait of Dover and 1692 northern Belgian Ardennes earthquakes, is very elusive if it exists. Hence, similar earthquakes would be very infrequent at these locations suggesting that the seismicity outside of the LRG would be episodic and clustered on some faults during periods of a few hundreds of years interrupted by long periods of inactivity typically lasting for some tens to hundreds of thousand years. Seismic moment release estimation and its comparison between recent geological and historical seismicity periods lead us to suggest that the high seismicity level observed between AD 1350 and AD 1700 west of the LRG would be uncommon.

Description: research

Description: In the weekly newspaper of Osnabrück (Germany) of November 3, 1770, a report about a local earthquake was published. Pastor Buck described ground motion effects in the manor ‘Haus Horst’, 1.5 km away from the village of Alfhausen: ‘roof tiles rattled, a chimney fell down, inside the house the top of a stove was overturned, abraded chalk trickled down in all rooms; in the nearby villages, people felt the shaking, and especially the churches suffered noticeably’. The epicentral intensity was estimated to VII (MSK) by Ahorner et al. (1970), but later modified to VI (EMS) by Meier & Grünthal (1992) considering Buck’s report in detail. Since this event is the only documented earthquake in this region, a reliable characterization of its parameters is important. Our re-examination reveals that some reported effects are quite contradictory. Contrary to Buck’s statement, no documents of damages on churches or costs of repairs could be found in the parish registers. As a result, the event appears to be a tectonic earthquake with an epicenter at Alfhausen / Haus Horst. Applying intensity-attenuation relationships, a revised value of the epicentral intensity of I_0≤V (EMS-98) with a focal depth of z≥2 km was derived. A cavity collapse due to leaching processes as a cause of the effects can be ruled out here. However, several details given in the primary source turned out to be unrealistic or at least exaggerated. The tectonic earthquake on September 3, 1770 near Alfhausen should be classified therefore as uncertain or even doubtful.

Description: research

Description: Central Europe is an intraplate domain which is characterized by low to moderate seismicity with records of larger seismic events occurring in historical and recent times. These records of seismicity are restricted to just over one thousand years. This does not reflect the long seismic cycles in Central Europe which are expected to be in the order of tens of thousands of years. Therefore, we have developed a paleoseismic database (PalSeisDB) that documents the records of paleoseismic evidence (trenches, soft-sediment deformation, mass movements, etc.) and extends the earthquake record to at least one seismic cycle. It is intended to serve as one important basis for future seismic hazard assessments. In the compilation of PalSeisDB, paleoseismic evidence features are documented at 129 different locations in the area of Germany and adjacent regions.

Description: research

Description: Current procedures to collect macroseismic data in Germany are diverse and scattered. At least 10 institutions collect macroseismic data by internet. Several institutes have a long tradition in seismology and have collected macroseismic data using paper forms for many decades. In addition, the responsibilities for geoscientific issues in Germany are a matter of the federal states and several of them have a state earthquake service. The only institution that automatically calculates and maps intensities online in near real time is Erdbebenstation Bensberg in cooperation with the Royal Observatory of Belgium. Baden-Württemberg uses a short form internet questionnaire at the moment. 5 state earthquake services (Bayern, Hessen, Niedersachsen, Rheinland-Pfalz, Sachsen) have implemented the standard German earthquake questionnaire (Kaiser 2014) which is adapted from the standard questionnaire developed by the ESC Working Group on Internet Macroseismology published by Musson & Cecić (2012). Most institutions express their strong need to implement standard procedures for automatic intensity assignment and a standard format for the exchange of questionnaire responses. References Kaiser, D. (2014): Der neue einheitliche Erdbeben-Fragebogen. Mitteilungen / Deutsche Geophysikalische Gesellschaft, 2/2014, 29-33. Musson, R. M. W. & Cecić, I. (2012): Intensity and Intensity Scales. In: New Manual of Seismological Observatory Practice 2 (NMSOP-2).- Deutsches GeoForschungsZentrum GFZ, 1-41; Potsdam. doi:10.2312/GFZ.NMSOP-2_ch12

Description: lecture

Description: We develop empirical relationships between the surface wave magnitude MS and macroseismic data, i.e. the epicentral intensity I0, isoseismal radii R(I) of different intensities I and the focal depth h. The basis of this study is formed by carefully selected instrumental parts (since 1900) of 2 earthquake catalogues: Kárník 1996 (Europe and the Mediterranean), and Shebalin et al. 1998 (Central and Eastern Europe). We use the orthogonal regression because we presume that all parameters are in error and because it has the advantage to provide a reversible regression equation. From Shebalin et al.1998 catalogue we obtain MS = 0.65 I0 + 1.90 log(h) – 1.62 with equivalent error δMS = ±0.21. In order to establish a relationship between MS and isoseismal radii we apply a theoretically based model which takes into account both exponential decay and geometrical spreading. From Shebalin et al. 1998catalogue we find MS = 0.673 I + 2.44 log (S(I)) + 0.00163 S(I) – 2.48 with δMS = ±0.28. Here I is the macroseismic intensity (I = 3…9) of the isoseismal in the focal distance S(I) [km]. Kárník 1996 gives isoseismal radii for I = 3 and 5. We obtain: MS = 0.808 I + 2.84 log (S(I)) + 0.00190 S(I) – 3.71 with δMS = ±0.65. These equations make possible reliably estimates of MS . We recommend them for application. The use of high quality data only as input in the regression analysis provides reliable relationships to estimate magnitudes. The magnitude estimation of a historical earthquake from the epicentral intensity gives reliable results only if the focal depth is known well enough. The relationship using isoseismal radii is of greater practical importance as it allows more reliable magnitude estimations of historical earthquakes. We observe regional variations in the relationships which need further investigation.

Description: poster

Description: The earthquake catalogue of southwestern Germany for the last millennium now contains about 30,000 digital macroseismic intensity data points (IDPs). Intensity assessments are based mainly on primary sources using the European Macroseismic Scale 1998 (EMS). The article describes a guideline for best practice of conventional macroseismic evaluation in application to historical and modern-time earthquakes in SW-Germany. Suitability of various diagnostics for intensity assessment is discussed. Assumptions to estimate damage grades and vulnerability classes of buildings are presented. Data restrictions and treatment of special cases are outlined. Further topics are quantification of uncertainties and IDP quality as well as substitutes for intensity. An essential task is to bridge the gap between information from historical sources and seismological needs for use in the earthquake catalogue, thus all issues have a focus on historical earthquakes. Questions of completeness, subjectivity, transparency, and interdisciplinary work are addressed also. Special emphasis is given to a well balanced use of the EMS scale throughout all time periods leading to consistent assessments in the catalogue.

Description: research

Description: Considering the importance of migration research and the dearth in literature on migration, the following study has been done, with the intention of examining the complex interrelation between internal migration and the development process in Bangladesh. Thereby, emphasis has been given to the investigation and analysis of the causes of rural-urban migration and its spatial as well as socioeconomic consequences, with examples from Dhaka. To construct a framework for the present study, some consideration has been given to both theoretical as well as empirical studies on internal migration. The discussion of general theories which follows in the next section, has mainly been done to assess the validity of these in the analysis of internal migration in the Third World countries.

Description: research

Description: DFG, SUB Göttingen