ALBERT

Beschreibung: By introducing two types of zenith troposphere delay (ZTD) products in precise point positioning (PPP), we developed the ZTD-corrected PPP and the ZTD-constrained PPP, both of them reduced the PPP convergence time. Both enhanced PPP methods are examined by global empirical ZTD models and regional ZTD corrections. For global ZTD models, we verified that ZTD-corrected PPP will deviate the positioning results, while ZTD-constrained PPP could produce unbiased estimations. Therefore, the latter is utilized to study the performance of global ZTD models (ITG, GPT2w, GZTD and UNB3m). After numerous experiments, we found that the performance of ZTD models was positively related to the real ZTD accuracy, and we proposed a universal tropospheric stochastic model 2SQR(9rms) which denotes double the square of nine times ZTD rms, to constrain ZTD in PPP. The proposed model subsequently was validated by real-time static and kinematic ZTD-constrained PPP on the premise that the ZTD rms on every station was known. Compared with traditional PPP, in static PPP, the number of improved stations is increased by 15.5 per cent (ITG), 14.4 per cent (GPT2w), 11.1 per cent (GZTD) and 8.3 per cent (UNB3m). For kinematic PPP, PPP constrained by ITG model still had the best performance, the number of improved stations is increased by 14.4 per cent, after 30 min of initialization time, 13.4 cm east, 13.4 cm north and 11.7 cm up positioning accuracy was obtained, compared with 15.3 cm east, 15.3 cm north and 14.3 cm up accuracy by traditional PPP. In addition, experiments using regional ZTD corrections to enhance real-time PPP showed that both ZTD-corrected PPP and ZTD-constrained PPP can notably reduce the convergence time on the vertical component (within 15 cm).

Beschreibung: We study fluctuations in the degree-2 zonal spherical harmonic coefficient of the Earth's gravity potential, C 20 , over the period 2003–2015. This coefficient is related to the Earth's oblateness and studying its temporal variations, C 20 , can be used to monitor large-scale mass movements between high and low latitude regions. We examine C 20 inferred from six different sources, including satellite laser ranging (SLR), GRACE and global geophysical fluids models. We further include estimates that we derive from measured variations in the length-of-day (LOD), from the inversion of global crustal displacements as measured by GPS, as well as from the combination of GRACE and the output of an ocean model as described by Sun et al. We apply a sequence of trend and seasonal moving average filters to the different time-series in order to decompose them into an interannual, a seasonal and an intraseasonal component. We then perform a comparison analysis for each component, and we further estimate the noise level contained in the different series using an extended version of the three-cornered-hat method. For the seasonal component, we generally obtain a very good agreement between the different sources, and except for the LOD-derived series, we find that over 90 per cent of the variance in the seasonal components can be explained by the sum of an annual and semiannual oscillation of constant amplitudes and phases, indicating that the seasonal pattern is stable over the considered time period. High consistency between the different estimates is also observed for the intraseasonal component, except for the solution from GRACE, which is known to be affected by a strong tide-like alias with a period of about 161 d. Estimated interannual components from the different sources are generally in agreement with each other, although estimates from GRACE and LOD present some discrepancies. Slight deviations are further observed for the estimate from the geophysical models, likely to be related to the omission of polar ice and groundwater changes in the model combination we use. On the other hand, these processes do not seem to play an important role at seasonal and shorter timescales, as the sum of modelled atmospheric, oceanic and hydrological effects effectively explains the observed C 20 variations at those scales. We generally obtain very good results for the solution from SLR, and we confirm that this well-established technique accurately tracks changes in C 20 . Good agreement is further observed for the estimate from the GPS inversion, showing that this indirect method is successful in capturing fluctuations in C 20 on scales ranging from intra- to interannual. Obtaining accurate estimates from LOD, however, remains a challenging task and more reliable models of atmospheric wind fields are needed in order to obtain high-quality C 20 , in particular at the seasonal scale. The combination of GRACE data and the output of an ocean model appears to be a promising approach, particularly since corresponding C 20 is not affected by tide-like aliases, and generally gives better results than the solution from GRACE, which still seems to be of rather poor quality.

Beschreibung: Archaeomagnetic field models cover longer timescales than historical models and may therefore resolve the motion of geomagnetic features on the core–mantle boundary (CMB) in a more meaningful statistical sense. Here we perform a detailed appraisal of archaeomagnetic field models to infer some aspects of the physics of the outer core. We characterize and compare the identification and tracking of reversed flux patches (RFPs) in order to assess the RFPs robustness. We find similar behaviour within a family of models but differences among different families, suggesting that modelling strategy is more influential than data set. Similarities involve recurrent positions of RFPs, but no preferred direction of motion is found. The tracking of normal flux patches shows similar qualitative behaviour confirming that RFPs identification and tracking is not strongly biased by their relative weakness. We also compare the tracking of RFPs with that of the historical field model gufm1 and with seismic anomalies of the lowermost mantle to explore the possibility that RFPs have preferred locations prescribed by lower mantle lateral heterogeneity. The archaeomagnetic field model that most resembles the historical field is interpreted in terms of core dynamics and core–mantle thermal interactions. This model exhibits correlation between RFPs and low seismic shear velocity in co-latitude and a shift in longitude. These results shed light on core processes, in particular we infer toroidal field lines with azimuthal orientation below the CMB and large fluid upwelling structures with a width of about 80° (Africa) and 110° (Pacific) at the top of the core. Finally, similar preferred locations of RFPs in the past 9 and 3 kyr of the same archaeomagnetic field model suggest that a 3 kyr period is sufficiently long to reliably detect mantle control on core dynamics. This allows estimating an upper bound of 220–310 km for the magnetic boundary layer thickness below the CMB.

Beschreibung: Post-depositional reductive diagenesis usually results in partial or entire cleansing of the pristine palaeomagnetic signal, therefore, its intensity is important to be assessed for sediments that are in the purpose of retrieving palaeomagnetic information. Grain size, rock magnetic and geochemical studies on the entire core, along with scanning electron microscope observations and X-ray diffraction analyses for representative samples were carried out on a Holocene sediment core retrieved from the deep water part of Huguangyan maar lake (HGY), southeast China. The pristine magnetic mineral assemblage of the studied core is domianted by superparamagnetic (SP) and stable single domain titanomagnetite, and high coercivity minerals are not detectable. Based on down-core variations of the average grain size ( M Z ), total organic carbon (TOC), detrital elements (Al, Ti, Fe and Mn) and the concentration and mineralogy of magnetic minerals, the studied core could be divided into three subsections. The uppermost subsection is the least affected by diagenesis, with detrital titanomagnetite as the dominant magnetic mineral. This is owing to low TOC contents, but high detrital input generated by weak Asian summer monsoon intensity during the late Holocene. The intermediate subsection shows down-core progressively enhanced dissolution of detrital titanomagnetite, and concomitant formation of authigenic pyrite and siderite, which indicates down-core progressively enhanced diagenesis generated by down-core progressive increasing TOC content, but decreasing detrital input as the result of down-core progressively strengthened Asian summer monsoon intensity. The pristine magnetic mineral assemblage has been profoundly modified in the lowermost subsection. At certain positions of the lowermost subsection, detrital titanomagnetite has been even completely dissolved via diagenesis, giving place to authigenic pyrite and siderite. High TOC content, but low detrital input generated from strong Asian summer monsoon intensity during the early Holocene are accountable for intensive diagenesis in the lowermost subsection. Complete erasing of detrital magnetic input signal at certain positions of the lowermost subsection, and considerable formation of authigenic siderite indicate that palaeomagnetic records of the studied core have been significantly compromised. The studied core has relatively higher TOC content, lower detrital matter content, calmer sedimentary environments, and less DO available at its water–sediment interface than the cores retrieved at relatively shallower water depths, which all contribute to its relatively stronger diagenesis. Progressive thickening of the upper two subsections with increasing water depth is owing to progressive increase in sedimentation rate with increasing water depth, which is the key factor in determining the thickness of each diagenetic subsection of cores from HGY. It would be better that lake sediments for palaeomagnetic investigations collected at a water depth shallower than the depth of its thermocline.

Beschreibung: The geocentre motion is the motion of the centre of mass of the entire Earth, considered an isolated system, in a terrestrial system of reference. We first derive a formula relating the harmonic degree-1 Lagrangian variation of the gravity at a station to both the harmonic degree-1 vertical displacement of the station and the displacement of the whole Earth's centre of mass. The relationship is independent of the nature of the Earth deformation and is valid for any source of deformation. We impose no constraint on the system of reference, except that its origin must initially coincide with the centre of mass of the spherically symmetric Earth model. Next, we consider the geocentre motion caused by surface loading. In a system of reference whose origin is the centre of mass of the solid Earth, we obtain a specific relationship between the gravity variation at the surface, the geocentre displacement and the load Love number $h^{\prime }_1$ , which demands the Earth's structure and rheological behaviour be known. For various networks of real or fictitious stations, we invert synthetic signals of surface gravity variations caused by atmospheric loading to retrieve the degree-1 variation of gravity. We then select six well-distributed stations of the Global Geodynamics Project, which is a world network of superconducting gravimeters, to invert actual gravity data for the degree-1 variations and determine the geocentre displacement between the end of 2004 and the beginning of 2012, assuming it to be due to surface loading. We find annual and semi-annual displacements with amplitude 0.5–2.3 mm.

Beschreibung: The rocks in the crust and the upper mantle of the Earth are believed to exhibit electrical anisotropy to some extent. It is beneficial to further understand and recognize the propagation of the electromagnetic waves in the Earth by investigating the magnetotelluric (which is one of the main geophysical techniques to probe the deep structures in the Earth) responses of the media with anisotropic conductivity structures. In this study, we examine the magnetotelluric fields over an idealized 2-D model consisting of two segments with axially anisotropic conductivity structures overlying a perfect conductor basement by a quasi-static analytic approach. The resulting analytic solution could not only contribute to the electromagnetic induction theory in the anisotropic Earth but also serve as at least an initial standard solution which could be used to validate the reliability and accuracy of the numerical algorithms developed for modelling the magnetotelluric responses of the 2-D media with much more general anisotropic conductivity.

Beschreibung: We review the theory of the Earth's elastic and gravitational response to a surface disk load. The solutions for displacement of the surface and the geoid are developed using expansions of Legendre polynomials, their derivatives and the load Love numbers. We provide a matlab  function called diskload that computes the solutions for both uncompensated and compensated disk loads. In order to numerically implement the Legendre expansions, it is necessary to choose a harmonic degree, n max , at which to truncate the series used to construct the solutions. We present a rule of thumb (ROT) for choosing an appropriate value of n max , describe the consequences of truncating the expansions prematurely and provide a means to judiciously violate the ROT when that becomes a practical necessity.

Beschreibung: Several works have reported that haematite has non-linear initial susceptibility at room temperature, like pyrrhotite or titanomagnetite, but there is no explanation for the observed behaviours yet. This study sets out to determine which physical property (grain size, foreign cations content and domain walls displacements) controls the initial susceptibility. The performed measurements include microprobe analysis to determine magnetic phases different to haematite; initial susceptibility (300 K); hysteresis loops, SIRM and backfield curves at 77 and 300 K to calculate magnetic parameters and minor loops at 77 K, to analyse initial susceptibility and magnetization behaviours below Morin transition. The magnetic moment study at low temperature is completed with measurements of zero field cooled–field cooled and AC susceptibility in a range from 5 to 300 K. The minor loops show that the non-linearity of initial susceptibility is closely related to Barkhausen jumps. Because of initial magnetic susceptibility is controlled by domain structure it is difficult to establish a mathematical model to separate magnetic subfabrics in haematite-bearing rocks.

Beschreibung: Two types of signals are clearly visible in continuous GPS (cGPS) time-series in Iceland, in particular in the vertical component. The first one is a yearly seasonal cycle, usually sinusoid-like with a minimum in the spring and a maximum in the fall. The second one is a trend of uplift, with higher values the closer the cGPS stations are to the centre of Iceland and ice caps. Here, we study the seasonal cycle signal by deriving its average at 71 GPS sites in Iceland. We estimate the annual and semi-annual components of the cycle in their horizontal and vertical components using a least-squares adjustment. The peak-to-peak amplitude of the cycle of the vertical component at the studied sites ranges from 4 mm near the coastline up to 27 mm at the centre of the Vatnajökull, the largest ice cap in Iceland. The minimum of the seasonal cycle occurs earlier in low lying areas than in the central part of Iceland, consistent with snow load having a large influence on seasonal deformation. Modelling shows that the seasonal cycle is well explained by accounting for elastically induced surface displacements due to snow, atmosphere, reservoir lake and ocean variations. Model displacement fields are derived considering surface loads on a multilayered isotropic spherical Earth. Through forward and inverse modelling, we were able to reproduce a priori information on the average seasonal cycle of known loads (atmosphere, snow in non-glaciated areas and lake reservoir) and get an estimation of other loads (glacier mass balance and ocean). The seasonal glacier mass balance cycle in glaciated areas and snow load in non-glaciated areas are the main contributions to the seasonal deformation. For these loads, induced seasonal vertical displacements range from a few millimetres far from the loads in Iceland, to more than 20 mm at their centres. Lake reservoir load also has to be taken into account on local scale as it can generate up to 20 mm of vertical deformation. Atmosphere load and ocean load are observable and generate vertical displacements in the order of a few millimetres. Inversion results also shows that the Iceland crust is less rigid than the world average. Interannual deviation from the GPS seasonal cycle can occur and are caused by unusual weather conditions over extended period of time.

Beschreibung: A series of important geological events occurred in the Tibetan Plateau area during the Jurassic, such as the collision of the Lhasa and Qiangtang Terranes, the closure of the Meso-Tethyan Ocean, the opening of the Neo-Tethyan Ocean and the cessation of the mega-monsoon. The ~3000 m thick Jurassic sedimentary sequence in the Qiangtang Basin on the central Tibetan Plateau, which is called the Yanshiping (YSP) Group, recorded these geological events. However, the chronology of the sequence is surprisingly poorly constrained. Here, we perform a detailed palaeomagnetic analysis on the ~1060 m thick middle and upper portions of the YSP Group (the Xiali and Suowa Formations) in the YSP section of the eastern Qiangtang Basin. Three bivalve zones at stratigraphic intervals of ~40–140, 640–800 and 940–1040 m are identified, which yield a Bathonian–Callovian age for the Lower Xiali Fm., a Callovian–Oxfordian age for the Lower Suowa Fm. and an Oxfordian–Kimmeridgian age for the Upper Suowa Fm. A total of 544 oriented palaeomagnetic samples were collected from the section. By combining thermal and alternating field demagnetizations, clear characteristic remanent magnetization (ChRM) directions are isolated for most of the samples. The robust ChRM directions pass fold and reversals tests, which support the primary nature of the ChRMs and yield a palaeopole at 76.8°N/297.2°E (dp = 2.2°, dm = 3.7°). A total of 27 normal and 26 reversed polarity zones were successfully recorded in the section. Combined with fossil age constraints, results suggest that the section is plausibly composed of a Callovian-Early Kimmeridgian age sedimentary sequence.

Beschreibung: Seismic waves produced by fault ruptures give rise to gravity perturbations. So far, these perturbations have either been modelled as permanent coseismic gravity change in a half-space or spherical Earth model, or as full time-domain model in infinite space. In this paper, we present the explicit solution of gravity perturbations in time domain produced by a double-couple buried in a homogeneous half-space. This result is especially suited to study gravity perturbations up to a few hundreds of kilometres from the epicentre. It facilitates detailed parametric studies of gravity perturbations from fault rupture, and predicts gravity perturbations of real earthquakes with greatly improved accuracy. The results may serve to develop first designs of gravity-assisted earthquake early-warning systems, made possible by a new generation of ultrasensitive gravity gradiometers, which is currently under development.

Beschreibung: Data on the evolution of Earth's magnetic field intensity are important for understanding the geodynamo and planetary evolution. However, the paleomagnetic record in rocks may be adversely affected by many physical processes, which must be taken into account when analysing the palaeointensity database. This is especially important in the light of an ongoing debate regarding core thermal conductivity values, and how these relate to the Precambrian geodynamo. Here, we demonstrate that several data sets in the Precambrian palaeointensity database overestimate the true paleofield strength due to the presence of non-ideal carriers of palaeointensity signals and/or viscous re-magnetizations. When the palaeointensity overestimates are removed, the Precambrian database does not indicate a robust change in geomagnetic field intensity during the Mesoproterozoic. These findings call into question the recent claim that the solid inner core formed in the Mesoproterozoic, hence constraining the thermal conductivity in the core to ‘moderate’ values. Instead, our analyses indicate that the presently available palaeointensity data are insufficient in number and quality to constrain the timing of solid inner core formation, or the outstanding problem of core thermal conductivity. Very young or very old inner core ages (and attendant high or low core thermal conductivity values) are consistent with the presently known history of Earth's field strength. More promising available data sets that reflect long-term core structure are geomagnetic reversal rate and field morphology. The latter suggests changes that may reflect differences in Archean to Proterozoic core stratification, whereas the former suggest an interval of geodynamo hyperactivity at ca. 550 Ma.

Beschreibung: Continuous gravimetric observations have been made with three successive generations of superconducting gravimeter over 20 yr at Syowa Station ( $39.6\deg$ E, $69.0\deg$ S), East Antarctica. The third-generation instrument, OSG#058, was installed in January 2010 and was calibrated by an absolute gravimeter during January and February, 2010. The estimated scale factor was –73.823 ± 0.053 μGal V –1 (1 μGal = 10 –8 m s –2 ). The first 5 yr of OSG#058 data from 2010 January 7 to 2015 January 10 were decomposed into tidal waves (M3 to Ssa) and other non-tidal components by applying the Bayesian tidal analysis program BAYTAP. Long-term non-tidal gravity residuals, which were obtained by subtracting annual and 18.6 year tidal waves and the predicted gravity response to the Earth's variable rotation, showed significant correlation with the accumulated snow depth measured at Syowa Station. The greatest correlation occurred when the gravity variations lagged the accumulated snow depth by 21 d. To estimate the gravitational effect of the accumulated snow mass, we inferred a conversion factor of 3.13 ± 0.08 μGal m –1 from this relation. The accumulated snow depth at Syowa Station was found to represent an extensive terrestrial water storage (the snow accumulation) around Syowa Station, which was estimated from the Gravity Recovery and Climate Experiment satellite gravity data. The snow accumulation around Syowa Station was detectable by the superconducting gravimeter.

Beschreibung: The new release AIUB-RL02 of monthly gravity models from GRACE GPS and K-Band range-rate data is based on reprocessed satellite orbits referring to the reference frame IGb08. The release is consistent with the IERS2010 conventions. Improvements with respect to its predecessor AIUB-RL01 include the use of reprocessed (RL02) GRACE observations, new atmosphere and ocean dealiasing products (RL05), an upgraded ocean tide model (EOT11A), and the interpolation of shallow ocean tides (admittances). The stochastic parametrization of AIUB-RL02 was adapted to include daily accelerometer scale factors, which drastically reduces spurious signal at the 161 d period in C 20 and at other low degree and order gravity field coefficients. Moreover, the correlation between the noise in the monthly gravity models and solar activity is considerably reduced in the new release. The signal and the noise content of the new AIUB-RL02 monthly gravity fields are studied and calibrated errors are derived from their non-secular and non-seasonal variability. The short-period time-variable signal over the oceans, mostly representing noise, is reduced by 50 per cent with respect to AIUB-RL01. Compared to the official GFZ-RL05a and CSR-RL05 monthly models, the AIUB-RL02 stands out by its low noise at high degrees, a fact emerging from the estimation of seasonal variations for selected river basins and of mass trends in polar regions. Two versions of the monthly AIUB-RL02 gravity models, with spherical harmonics resolution of degree and order 60 and 90, respectively, are available for the time period from March 2003 to March 2014 at the International Center for Global Earth Models or from ftp://ftp.unibe.ch/aiub/GRAVITY/GRACE (last accessed 22 March 2016).

Beschreibung: Red Clay underlying the loess-palaeosol sequences on the Chinese Loess Plateau is an eolian deposit. There is a controversy over whether magnetic susceptibility ( ) variations in Red Clay sequence can be used as an indicator of summer palaeomonsoon intensity. This study investigates the magnetic mineralogy, magnetic concentration and magnetic grain size distribution of Jiaxian Red Clay with multimagnetic methods. Our results indicate that the magnetic properties of Jiaxian Red Clay are similar to those of the Quaternary loess-palaeosol sequences, and ultrafine ferrimagnetic grains produced during pedogenesis are responsible for an increase in susceptibility, therefore the enhancement mechanism of Red Clay is similar to that of the overlying loess-palaeosol sequences. This paper explores variations in the Red Clay sequence through spatial and temporal analysis. The susceptibility variation of six sites along a NNE to SSW transect correlate to palaeoclimatic cycles, so can be used to trace the summer palaeomonsoon intensity from a spatial perspective. However, a simple loess-derived calibration function cannot be used to quantitative reconstruct the palaeomonsoon intensity variations thought time. An adjusted calibration function for palaeosols from Red Clay sequence needs to be developed, so that can be used to quantitative reconstruct palaeomonsoon intensity. Further study is necessary to develop such a transfer function.

Beschreibung: Palaeointensity experiments were carried out to a sample collection from two sections of basalt lava flow sequences of Pliocene age in north central Iceland (Chron C2An) to further refine the knowledge of the behaviour of the palaeomagnetic field. Selection of samples was mainly based on their stability of remanence to thermal demagnetization as well as good reversibility in variations of magnetic susceptibility and saturation magnetization with temperature, which would indicate the presence of magnetite as a product of deuteric oxidation of titanomagnetite. Among 167 lava flows from two sections, 44 flows were selected for the Königsberger–Thellier–Thellier experiment in vacuum. In spite of careful pre-selection of samples, an Arai plot with two linear segments, or a concave-up appearance, was often encountered during the experiments. This non-ideal behaviour was probably caused by an irreversible change in the domain state of the magnetic grains of the pseudo-single-domain (PSD) range. This is assumed because an ideal linear plot was obtained in the second run of the palaeointensity experiment in which a laboratory thermoremanence acquired after the final step of the first run was used as a natural remanence. This experiment was conducted on six selected samples, and no clear difference between the magnetic grains of the experimented and pristine sister samples was found by scanning electron microscope and hysteresis measurements, that is, no occurrence of notable chemical/mineralogical alteration, suggesting that no change in the grain size distribution had occurred. Hence, the two-segment Arai plot was not caused by the reversible multidomain/PSD effect in which the curvature of the Arai plot is dependent on the grain size. Considering that the irreversible change in domain state must have affected data points at not only high temperatures but also low temperatures, f v ≥ 0.5 was adopted as one of the acceptance criteria where f v is a vectorially defined fraction of the linear segment. A measure of curvature k ' was also used to check the linearity of the selected linear segment. It was avoided, however, to reject the result out of hand by the large curvature k of the entire data points because it might still include a linear segment with a large fraction. Combining with the results of Shaw's experiments, 52 palaeointensities were obtained out of 192 specimens, or 11 flow means were obtained out of the 44 lava flows. Most of the palaeointensities were from the upper part of the lava section (Chron C2An.1n) and ranged between 30 and 66 μT. Including two results from the bottom part of the lava section, the mean virtual dipole moment for 2.5–3.5 Ma is 6.3 ± 1.4  x  10 22 Am 2 ( N  = 11), which is ~19 per cent smaller than the present-day dipole moment.

Beschreibung: We document two kinds of traveling ionospheric disturbances, namely, CTIDs (Co-tsunami-Traveling-Ionospheric-disturbances) and ATIDs (Ahead-of-Tsunami-Traveling-Ionospheric-disturbances) related to the Tohoku-Oki tsunami of 2011 March 11. They are referred to the disturbances that remain behind and ahead of the principal tsunami wave front, respectively. We first note their presence in a numerical experiment performed using a simulation code coupling the tsunami, atmosphere and ionosphere. This code uses the tsunami wavefield as an input and simulates acoustic-gravity waves (AGWs) in the atmosphere and TIDs, in the form of total electron content (TEC) disturbance, in the ionosphere. The simulated TEC reveals the excitation of CTIDs (at about 2 TECU) and ATIDs (at about 1 TECU), representing up to 5 per cent disturbance over the ambient electron density, and they arise from the dissipation of AGWs in the thermosphere. A novel outcome is that during the tsunami passage between ~6° and 12° of epicentral distance, strong ATIDs arrive ~20–60 min ahead of the tsunami wave front covering ~3°–10° of distance from the tsunami location. Simulation results are compared with the far-field observations using GNSS satellites and confirm that ATIDs are the first detected TEC maximum, occurring 20–60 min ahead of the tsunami arrival. Our simulation also confirms the presence of largest TEC maximum representing CTIDs, 10–20 min after the first tsunami wave. ATIDs reported in this study have characteristics that can be potentially used for the early warning of the tsunami.

Beschreibung: Apparent acceleration in Gravity Recovery and Climate Experiment (GRACE) Antarctic ice mass time-series may reflect both ice discharge and surface mass balance contributions. However, a recent study suggests there is also contamination from errors in atmospheric pressure de-aliasing fields [European Center for Medium-Range Weather Forecast (ECMWF) operational products] used during GRACE data processing. To further examine this question, we compare GRACE atmospheric pressure de-aliasing (GAA) fields with in situ surface pressure data from coastal and inland stations. Differences between the two are likely due to GAA errors, and provide a measure of error in GRACE solutions. Time-series of differences at individual weather stations are fit to four presumed error components: annual sinusoids, a linear trend, an acceleration term and jumps at times of known ECMWF model changes. Using data from inland stations, we estimate that atmospheric pressure error causes an acceleration error of about +7.0 Gt yr –2 , which is large relative to prior GRACE estimates of Antarctic ice mass acceleration in the range of –12 to –14 Gt yr –2 . We also estimate apparent acceleration rates from other barometric pressure (reanalysis) fields, including ERA-Interim, MERRA and NCEP/DOE. When integrated over East Antarctica, the four mass acceleration estimates (from GAA and the three reanalysis fields) vary considerably (by ~2–16 Gt yr –2 ). This shows the need for further effort to improve atmospheric mass estimates in this region of sparse in situ observations, in order to use GRACE observations to measure ice mass acceleration and related sea level change.

Beschreibung: Understanding the contribution of biogenic magnetic particles into sedimentary assemblages is a current challenge in palaeomagnetism. It has been demonstrated recently that magnetic particles produced through biologically controlled mineralization processes, such as magnetosomes from magnetotactic bacteria, contribute to the recording of natural remanent magnetization in marine and lacustrian sediments. Contributions from other, biologically induced, mineralization types, which are known from multiple laboratory experiments to include magnetic minerals, remain largely unknown. Here, we report magnetic properties of iron minerals formed by a community of iron- and manganese-reducing bacteria isolated from a natural groundwater deposit during a 2 yr long incubation experiment. The main iron phases of the biomineralized mass are lepidocrocite, goethite and magnetite, each of which has environmental significance. Unlike the majority of the previous studies that reported superparamagnetic grain size, and thus no remanence carrying capacity of biologically induced magnetite, hysteresis and first-order reversal curves measurements in our study have not detected significant superparamagnetic contribution. The biomineralized mass, instead, contains a mixture of single-domain to pseudo-single-domain and multidomain magnetite particles that are capable of carrying a stable chemical remanent magnetization. Isothermal remanent magnetization acquisition parameters and first-order reversal curves signatures of the biomineralized samples deviate from previously proposed criteria for the distinction of extracellular (biologically induced) magnetic particles in mixtures. Given its potential significance as a carrier of natural remanent magnetization, environmental requirements, distribution in nature and the efficiency in the geomagnetic field recording by biologically induced mineralization need comprehensive investigation.

Beschreibung: The long-wavelength gravity field contains information about processes in the sublithospheric mantle. As satellite-derived gravity models now provide the long to medium-wavelength gravity field at unprecedented accuracy, techniques used to process gravity data need to be updated. We show that when determining these long-wavelengths, the treatment of topographic-isostatic effect (TIE) and isostatic effects (IE) is a likely source of error. We constructed a global isostatic model and calculated global TIE and IE. These calculations were done for ground stations as well as stations at satellite height. We considered both gravity and gravity gradients. Using these results, we determined how much of the gravity signal comes from distant sources. We find that a significant long-wavelength bias is introduced if far-field effects on the topographic effect are neglected. However, due to isostatic compensation far-field effects of the topographic effect are to a large degree compensated by the far-field IE. This means that far-field effects can be reduced effectively by always considering topographic masses together with their compensating isostatic masses. We show that to correctly represent the ultra-long wavelengths, a global background model should be used. This is demonstrated both globally and for a continental-scale case area in North America. In the case of regional modelling, where the ultra-long wavelengths are not of prime importance, gravity gradients can be used to help minimize correction errors caused by far-field effects.

Beschreibung: Constraining laterally varying structures in planetary interiors is important for understanding both the composition and the internal dynamics of a planet. Recognizing that seismic imaging technique is currently only viable for studying the Earth's interior structures, methods that can be supported by advanced space geodetic techniques may become alternatives to ‘image’ the interiors of other planets. The method of tidal tomography is one possibility, and it relies on high precision measurement of the response of a planet to its body tide. However, it is essential to develop an efficient analytical tool that computes the dependence of tidal response to 3-D interior structures. In this paper, we present a complete formulation of such an analytical tool, which calculates to high accuracy the tidal response of a terrestrial planet with lateral heterogeneities in its elastic and density structures. We treat the lateral heterogeneities as small perturbations and derive the governing equations based on the perturbation theory. In a spherical harmonic representation, equations at each order of perturbation are reduced into multiple matrix equations at harmonics that are allowed by mode couplings, and the total response equals the sum of all those single-harmonic responses, which can be solved semi-analytically. We test our perturbation method by applying it to the Moon with a harmonic degree-1 mantle structure for which the perturbation solutions of the tidal response are compared with those from a fully numerical method. The remarkable agreement between results from these two methods validates the perturbation method. As an example, we then use the perturbation method to evaluate the impact of lunar crustal thickness variations on tidal response of the Moon. We find that lunar crust produces much smaller degree-3 tidal responses than a relatively weak degree-1 structure in the deep lunar mantle. Our calculations show that degree-3 tidal response measurements may hold key constraints on possible degree-1 mantle structure of the Moon, as suggested from previous modelling results.

Beschreibung: Seismic data are primarily used in studies of the Earth's lithospheric structure including the Moho geometry. In regions, where seismic data are sparse or completely absent, gravimetric or combined gravimetric-seismic methods could be applied to determine the Moho depth. In this study, we derive and present generalized expressions for solving the Vening Meinesz–Moritz's (VMM) inverse problem of isostasy for a Moho depth determination from gravity and vertical gravity-gradient data. By solving the (non-linear) Fredholm's integral equation of the first kind, the linearized observation equations, which functionally relate the (given) gravity/gravity-gradient data to the (unknown) Moho depth, are derived in the spectral domain. The VMM gravimetric results are validated by using available seismic and gravimetric Moho models. Our results show that the VMM Moho solutions obtained by solving the VMM problem for gravity and gravity-gradient data are almost the same. This finding indicates that in global applications, using the global gravity/gravity-gradient data coverage, the spherical harmonic expressions for the gravimetric forward and inverse modelling yield (theoretically) the same results. Globally, these gravimetric solutions have also a relatively good agreement with the CRUST1.0 and GEMMA GOCE models in terms of their rms Moho differences (4.7 km and 4.1 km, respectively).

Beschreibung: Monthly means of the magnetic field measurements at ground observatories are a key data source for studying temporal changes of the core magnetic field. However, when they are calculated in the usual way, contributions of external (magnetospheric and ionospheric) origin may remain, which make them less favourable for studying the field generated by dynamo action in the core. We remove external field predictions, including a new way of characterizing the magnetospheric ring current, from the data and then calculate revised monthly means using robust methods. The geomagnetic secular variation (SV) is calculated as the first annual differences of these monthly means, which also removes the static crustal field. SV time-series based on revised monthly means are much less scattered than those calculated from ordinary monthly means, and their variances and correlations between components are smaller. On the annual to decadal timescale, the SV is generated primarily by advection in the fluid outer core. We demonstrate the utility of the revised monthly means by calculating models of the core surface advective flow between 1997 and 2013 directly from the SV data. One set of models assumes flow that is constant over three months; such models exhibit large and rapid temporal variations. For models of this type, less complex flows achieve the same fit to the SV derived from revised monthly means than those from ordinary monthly means. However, those obtained from ordinary monthly means are able to follow excursions in SV that are likely to be external field contamination rather than core signals. Having established that we can find models that fit the data adequately, we then assess how much temporal variability is required. Previous studies have suggested that the flow is consistent with torsional oscillations (TO), solid body-like oscillations of fluid on concentric cylinders with axes aligned along the Earth's rotation axis. TO have been proposed to explain decadal timescale changes in the length-of-day. We invert for flow models where the only temporal changes are consistent with TO, but such models have an unacceptably large data misfit. However, if we relax the TO constraint to allow a little more temporal variability, we can fit the data as well as with flows assumed constant over three months, demonstrating that rapid SV changes can be reproduced by rather small flow changes. Although the flow itself changes slowly, its time derivative can be locally (temporally and spatially) large, in particular when and where core surface secular acceleration peaks. Spherical harmonic expansion coefficients of the flows are not well resolved, and many of them are strongly correlated. Averaging functions, a measure of our ability to determine the flow at a given location from the data distribution available, are poor approximations to the ideal, even when centred on points of the core surface below areas of high observatory density. Both resolution and averaging functions are noticeably worse for the toroidal flow component, which dominates the flow, than the poloidal flow component, except around the magnetic equator where averaging functions for both components are poor.

Beschreibung: Geophysical techniques are widely used to monitor volcanic unrest. A number of studies have also demonstrated that hydrological processes can produce or trigger geophysical signals. Hydrologically induced gravity signals have previously been recorded by specifically designed gravity surveys as well as, inadvertently, by volcano monitoring studies. Water table corrections of microgravity surveys are commonplace. However, the fluctuations of the water table beneath survey locations are often poorly known, and such a correction fails to account for changes in water-mass storage in the unsaturated zone. Here, we combine 2-D axis-symmetrical numerical fluid-flow models with an axis-symmetric, distributed-mass, gravity calculation to model gravity changes in response to fluctuating hydrological recharge. Flow simulations are based on tropical volcanic settings where high surface permeabilities promote thick unsaturated zones. Our study highlights that mass storage (saturation) changes within the unsaturated zone beneath a survey point can generate recordable gravity changes. We show that for a tropical climate, recharge variations can generate gravity variations of over 150 μGal; although, we demonstrate that for the scenarios investigated here, the probability of recording such large signals is low. Our modelling results indicate that microgravity survey corrections based on water table elevation may result in errors of up to 100 μGal. The effect of inter-annual recharge fluctuations dominate over seasonal cycles which makes prediction and correction of the hydrological contribution more difficult. Spatial hydrogeological heterogeneity can also impact on the accuracy of relative gravity surveys, and can even result in the introduction of additional survey errors. The loading fluctuations associated with saturation variations in the unsaturated zone may also have implications for other geophysical monitoring techniques, such as geodetic monitoring of ground deformation.

Beschreibung: Several attempts have been made to obtain a radiographic image inside volcanoes using cosmic-ray muons (muography). Muography is expected to resolve highly heterogeneous density profiles near the surface of volcanoes. However, several prior works have failed to make clear observations due to contamination by background noise. The background contamination leads to an overestimation of the muon flux and consequently a significant underestimation of the density in the target mountains. To investigate the origin of the background noise, we performed a Monte Carlo simulation. The main components of the background noise in muography are found to be low-energy protons, electrons and muons in case of detectors without particle identification and with energy thresholds below 1 GeV. This result was confirmed by comparisons with actual observations of nuclear emulsions. This result will be useful for detector design in future works, and in addition some previous works of muography should be reviewed from the view point of background contamination.

Beschreibung: Essential to understanding sea level change and its causes during the last interglacial (LIG) is the quantification of uncertainties. In order to estimate the uncertainties, we develop a statistical framework for the comparison of palaeoclimatic sea level index points and GIA model predictions. For the investigation of uncertainties, as well as to generate better model predictions, we implement a massive ensemble approach by applying a data assimilation scheme based on particle filter methods. The different runs are distinguished through varying ice sheet reconstructions based on oxygen-isotope curves and different parameter selections within the GIA model. This framework has several advantages over earlier work, such as the ability to examine either the contribution of individual observations to the results or the probability of specific input parameters. This exploration of input parameters and data leads to a larger range of estimates than previously published work. We illustrate how the assumptions that enter into the statistical analysis, such as the existence of outliers in the observational database or the initial ice volume history, can introduce large variations to the estimate of the maximum highstand. Thus, caution is required to avoid overinterpreting results. We conclude that there are reasonable doubts whether the data sets previously used in statistical analyses are able to tightly constrain the value of maximum highstand during the LIG.

Beschreibung: A 3-D density model of the crust and upper mantle beneath the Karoo basin is presented here. The model is constrained using potential field, borehole and seismic data. Uplift of the basin by the end of the Cretaceous has resulted in an unusually high plateau (〉1000 m) covering a large portion of South Africa. Isostatic studies show the topography is largely compensated by changes in Moho depths (~35 km on-craton and 〉45 km off-craton) and changes in lithospheric mantle densities between the Kaapvaal Craton and surrounding regions (~50 kg m –3 increase from on- to off-craton). This density contrast is determined by inverted satellite gravity and gravity gradient data. The highest topography along the edge of the plateau (〉1200 m) and a strong Bouguer gravity low over Lesotho, however, can only be explained by a buoyant asthenosphere with a density decrease of around 40 kg m –3 .

Beschreibung: A method to estimate the rotation change in the orientation of the centre-of-figure (CF) frame caused by earthquakes is proposed for the first time. This method involves using the point dislocation theory based on a spherical, non-rotating, perfectly elastic and isotropic (SNREI) Earth. The rotation change in the orientation is related solely to the toroidal displacements of degree one induced by the vertical dip slip dislocation, and the spheroidal displacements induced by an earthquake have no contribution. The effects of two recent large earthquakes, the 2004 Sumatra and the 2011 Tohoku-Oki, are studied. Results showed that the Sumatra and Tohoku-Oki earthquakes both caused the CF frame to rotate by at least tens of μas (micro-arc-second). Although the visible co-seismic displacements are identified and removed from the coordinate time-series, the rotation change due to the unidentified ones and errors in removal is non-negligible. Therefore, the rotation change in the orientation of the CF frame due to seismic deformation should be taken into account in the future in reference frame and geodesy applications.

Beschreibung: We derive scaling relationships for planetary dynamos based on a balance between energy production and Joule dissipation, and between the curl of the buoyancy and Coriolis forces. These scaling relationships are deduced for the particular case of dynamos driven by helical waves, but are shown to have a much broader applicability. They are consistent with the evidence of the numerical dynamos, yielding predictions consistent with published empirical scaling laws and also with the observed transition from dipolar to multipolar dynamos. A direct comparison with the observational evidence for the planets is hampered by the fact that we do not know what sets the smallest scale of the motion in the planets. Nevertheless, we use our scaling relationships to show that the traditional assumption that the Elsasser number is of order unity is inconsistent with the observation that the gas-giant dynamos are dipolar dynamos, as is the more recent suggestion that the strength of the dipole is independent of rotation rate and controlled by the buoyancy flux alone. On the other hand, we show that the observational data is consistent with the hypothesis that a dipolar dynamo saturates at the lowest permissible magnetic energy compatible with a given buoyancy flux.

Beschreibung: We present a detailed palaeomagnetic study from 35 sites on Holocene lava flows of the Tongariro Volcanic Centre, central North Island, New Zealand. Prior to the study the eruption ages of these flows were constrained to within a few thousand years by recently published high-precision 40 Ar/ 39 Ar geochronological data and tephrostratigraphic controls. Correlation of flow mean palaeomagnetic directions with a recently published continuous sediment record from Lake Mavora, Fiordland, allows us to reduce the age uncertainty to 300–500 yr in some cases. Our refined ages significantly improve the chronology of Holocene effusive eruptions of the volcanoes of the Tongariro Volcanic Centre. For instance, differences in the palaeomagnetic directions recorded by lavas from the voluminous Iwikau and Rangataua members suggest that individual effusive periods lasted up to thousands of years and that these bursts have been irregularly spaced over time. While over the last few millennia the effusive eruptive activity from Mt Ruapehu has been relatively quiet, the very young age (200–500 BP) of a Red Crater sourced flow suggests that effusive activity around Mt Tongariro lasted into the past few centuries. This adds an important hazard context to the historical record, which has otherwise comprised frequent relatively small, tephra producing, explosive eruptions without the production of lava flows.

Beschreibung: In this study, we propose an approach for determining the geopotential difference using high-frequency-stability microwave links between satellite and ground station based on Doppler cancellation system. Suppose a satellite and a ground station are equipped with precise optical-atomic clocks (OACs) and oscillators. The ground oscillator emits a signal with frequency f a towards the satellite and the satellite receiver (connected with the satellite oscillator) receives this signal with frequency f b which contains the gravitational frequency shift effect and other signals and noises. After receiving this signal, the satellite oscillator transmits and emits, respectively, two signals with frequencies f b and f c towards the ground station. Via Doppler cancellation technique, the geopotential difference between the satellite and the ground station can be determined based on gravitational frequency shift equation by a combination of these three frequencies. For arbitrary two stations on ground, based on similar procedures as described above, we may determine the geopotential difference between these two stations via a satellite. Our analysis shows that the accuracy can reach 1 m 2 s – 2 based on the clocks’ inaccuracy of about 10 –17 (s s –1 ) level. Since OACs with instability around 10 –18 in several hours and inaccuracy around 10 –18 level have been generated in laboratory, the proposed approach may have prospective applications in geoscience, and especially, based on this approach a unified world height system could be realized with one-centimetre level accuracy in the near future.

Beschreibung: The inversion of magnetotelluric data into subsurface electrical conductivity poses an ill-posed problem. Smoothing constraints are widely employed to estimate a regularized solution. Here, we present an alternative inversion scheme that estimates a sparse representation of the model in a wavelet basis. The objective of the inversion is to determine the few non-zero wavelet coefficients which are required to fit the data. This approach falls into the class of sparsity constrained inversion schemes and minimizes the combination of the data misfit in a least-squares 2 sense and of a model coefficient norm in an 1 sense ( 2 - 1 minimization). The 1 coefficient norm renders the solution sparse in a suitable representation such as the multiresolution wavelet basis, but does not impose explicit structural penalties on the model as it is the case for 2 regularization. The presented numerical algorithm solves the mixed 2 - 1 norm minimization problem for the nonlinear magnetotelluric inverse problem. We demonstrate the feasibility of our algorithm on synthetic 2-D MT data as well as on a real data example. We found that sparse models can be estimated by inversion and that the spatial distribution of non-vanishing coefficients indicates regions in the model which are resolved.

Beschreibung: Earth's liquid core hosts a diverse set of waves with periods ranging from days to thousands of years. One class of waves with periods of several decades is known to arise from an interplay between magnetic, Archimedes and Coriolis forces. These so-called MAC waves are thought to be relevant for interpreting historical fluctuations in the geomagnetic field. In this study, we show that MAC waves provide a good description of time-dependent zonal flow at the top of the core. The same collection of waves also offers a simple explanation for observed fluctuations in the dipole field. Both of these predictions require a stratified layer at the top of the core with a thickness of 130–140 km and a buoyancy frequency comparable to Earth's rotation rate. We extend these predictions to include changes in the length of day (LOD) and find that MAC waves can account for about half of the observed fluctuation at decadal periods. Larger fluctuations are possible when electromagnetic stresses couple MAC waves to flow in the interior of the core. In fact, an idealized model for the coupled motion overestimates the LOD fluctuations, probably reflecting limitations in this idealized model. Our results offer support for stable stratification at the top of the core and suggest a common origin for decadal fluctuations in the geomagnetic field and the LOD.

Beschreibung: A sequence of large earthquakes occurred along the North Anatolian fault in the 20th century. These earthquakes, including the 1999 Izmit/Düzce earthquakes, generally propagated westward towards the Marmara Sea, defining the Main Marmara fault as a potential seismic gap. It is important to conduct a detailed assessment of the seismic hazards along the main Marmara fault because the megacity Istanbul lies only approximately 10 km north of the eastern segment of the Main Marmara fault, which is referred to as the Princes’ Islands Fault segment (PIF). Here, we study the locking status of this fault segment to evaluate the seismic hazard potential. For the first time, combined ascending and descending Interferometric Synthetic Aperture Radar and Global Positioning System observations were used to investigate the crustal deformation associated with the PIF. After careful corrections of the estimated ground velocity, a deformation pattern relating to fault locking near the Princes’ Islands was identified. The modeling results revealed that the slip rate and locking depth of the fault segment show a clear trade-off, which were estimated as 18.9 ± 7.2 mm yr –1 and 12.1 ± 7.0 km, respectively. With a moment accumulation rate of 1.7 ± 0.4  x  10 17 Nm yr –1 (proportional to the product of slip rate and locking depth), our results imply a build-up of a geodetic moment on the PIF and therefore a potential for earthquake hazards in the vicinity of the Istanbul megacity.

Beschreibung: We present a distributed slip model for the 1999 M w 6.3 Chamoli earthquake of north India using interferometric synthetic aperture radar (InSAR) data from both ascending and descending orbits and Bayesian estimation of confidence levels and trade-offs of the model geometry parameters. The results of fault-slip inversion in an elastic half-space show that the earthquake ruptured a $9 _{ - 2.2}^{\circ + 3.4}$ northeast-dipping plane with a maximum slip of ~1 m. The fault plane is located at a depth of ~ $15.9_{ - 3.0}^{ + 1.1}$ km and is ~120 km north of the Main Frontal Thrust, implying that the rupture plane was on the northernmost detachment near the mid-crustal ramp of the Main Himalayan Thrust. The InSAR-determined moment is 3.35 x 10 18 Nm with a shear modulus of 30 GPa, equivalent to M w 6.3, which is smaller than the seismic moment estimates of M w 6.4–6.6. Possible reasons for this discrepancy include the trade-off between moment and depth, uncertainties in seismic moment tensor components for shallow dip-slip earthquakes and the role of earth structure models in the inversions. The released seismic energy from recent earthquakes in the Garhwal region is far less than the accumulated strain energy since the 1803 M s 7.5 earthquake, implying substantial hazard of future great earthquakes.

Beschreibung: Surface nuclear magnetic resonance technique, also called magnetic resonance sounding (MRS), is an emerging geophysical method that can detect the presence and spatial variations of the subsurface water content directly. In this paper, we introduce the MRS central loop geometry, in which the receiver loop is smaller than the transmitter loop and placed in its centre. In addition, using a shielded receiver coil we show how this configuration greatly increases signal-to-noise ratio and improves the resolution of the subsurface layers compared to the typically used coincident loop configuration. We compare sensitivity kernels for different loop configurations and describe advantages of the MRS central loop geometry in terms of superior behaviour of the sensitivity function, increased sensitivity values, reduced noise level of the shielded receiver coil, improved resolution matrix and reduced instrument dead time. With no extra time and effort in the field, central-loop MRS makes it possible to reduce measurement time and to measure data in areas with high anthropogenic noise. The results of our field example agree well with the complementary data, namely airborne electromagnetics, borehole data, and the hydrologic model of the area.

Beschreibung: Frictional heating during earthquake rupture reveals important information on earthquake mechanisms and energy dissipation. The amount of annealing varies widely and is, as yet, poorly constrained. Here we use magnetic susceptibility versus temperature measurements during cycling to increasingly elevated temperatures to constrain the maximum temperature a slip zone has experienced. The case study comprises sheared clay cored from the Japan Trench subduction plate-boundary fault zone (décollement), which accommodated the large slip of the 2011 M w 9.0 Tohoku-oki earthquake. The décollement was cored during the Integrated Ocean Drilling Program (IODP) Expedition 343, the Japan Trench Fast Drilling Project (JFAST). Heating signatures with estimated maximum temperatures ranging from ~300 to over 500 °C are determined close to the multiple slip surfaces within the décollement. Since it is impossible to tie a specific slip surface to a certain earthquake, thermal evidence for the cumulative effect of several earthquakes is unveiled. This as yet preliminary rock magnetic ‘geothermometer’ would be a useful tool to detect seismic heating along faults that experienced medium temperature rise, a range which is difficult to assess with other approaches.

Beschreibung: Feldspars are the most abundant rock-forming minerals in the Earth's crust, but their magnetic properties have not been rigorously studied. This work focuses on the intrinsic magnetic anisotropy of 31 feldspar samples with various chemical compositions. Because feldspar is often twinned or shows exsolution textures, measurements were performed on twinned and exsolved samples as well as single crystals. The anisotropy is controlled by the diamagnetic susceptibility and displays a consistent orientation of principal susceptibility axes; the most negative or minimum susceptibility is parallel to [010], and the maximum (least negative) is close to the crystallographic [001] axis. However, the magnetic anisotropy is weak when compared to other rock-forming minerals, 1.53 x 10 –9 m 3 kg –1 at maximum. Therefore, lower abundance minerals, such as augite, hornblende or biotite, often dominate the bulk paramagnetic anisotropy of a rock. Ferromagnetic anisotropy is not significant in most samples. In the few samples that do show ferromagnetic anisotropy, the principal susceptibility directions of the ferromagnetic subfabric do not display a systematic orientation with respect to the feldspar lattice. These results suggest that palaeointensity estimates of the geomagnetic field made on single crystals of feldspar will not be affected by a systematic orientation of the ferromagnetic inclusions within the feldspar lattice.

Beschreibung: Feldspars are the most abundant rock-forming minerals in the Earth's crust, but their magnetic properties have not been rigorously studied. This work focuses on the intrinsic magnetic anisotropy of 31 feldspar samples with various chemical compositions. Because feldspar is often twinned or shows exsolution textures, measurements were performed on twinned and exsolved samples as well as single crystals. The anisotropy is controlled by the diamagnetic susceptibility and displays a consistent orientation of principal susceptibility axes; the most negative or minimum susceptibility is parallel to [010], and the maximum (least negative) is close to the crystallographic [001] axis. However, the magnetic anisotropy is weak when compared to other rock-forming minerals, 1.53 x 10 –9 m 3 kg –1 at maximum. Therefore, lower abundance minerals, such as augite, hornblende or biotite, often dominate the bulk paramagnetic anisotropy of a rock. Ferromagnetic anisotropy is not significant in most samples. In the few samples that do show ferromagnetic anisotropy, the principal susceptibility directions of the ferromagnetic subfabric do not display a systematic orientation with respect to the feldspar lattice. These results suggest that palaeointensity estimates of the geomagnetic field made on single crystals of feldspar will not be affected by a systematic orientation of the ferromagnetic inclusions within the feldspar lattice.

Beschreibung: We use a Bayesian formalism combined with a grid node discretization for the linear inversion of gravimetric data in terms of 3-D density distribution. The forward modelling and the inversion method are derived from seismological inversion techniques in order to facilitate joint inversion or interpretation of density and seismic velocity models. The Bayesian formulation introduces covariance matrices on model parameters to regularize the ill-posed problem and reduce the non-uniqueness of the solution. This formalism favours smooth solutions and allows us to specify a spatial correlation length and to perform inversions at multiple scales. We also extract resolution parameters from the resolution matrix to discuss how well our density models are resolved. This method is applied to the inversion of data from the volcanic island of Basse-Terre in Guadeloupe, Lesser Antilles. A series of synthetic tests are performed to investigate advantages and limitations of the methodology in this context. This study results in the first 3-D density models of the island of Basse-Terre for which we identify: (i) a southward decrease of densities parallel to the migration of volcanic activity within the island, (ii) three dense anomalies beneath Petite Plaine Valley, Beaugendre Valley and the Grande-Découverte-Carmichaël-Soufrière Complex that may reflect the trace of former major volcanic feeding systems, (iii) shallow low-density anomalies in the southern part of Basse-Terre, especially around La Soufrière active volcano, Piton de Bouillante edifice and along the western coast, reflecting the presence of hydrothermal systems and fractured and altered rocks.

Beschreibung: We compute the gravimetric factor at the Chandler wobble (CW) frequency using time-series from superconducting gravimeters (SG) longer than a decade. We first individually process the polar motion and data at each individual gravity station to estimate the gravimetric factor amplitude and phase, then we make a global analysis by applying a stacking method to different subsets of up to seven SG stations. The stacking is an efficient way of getting rid of local effects and improving the signal-to-noise ratio of the combined data sets. Using the stacking method, we find a gravimetric factor amplitude and phase of 1.118 ± 0.016 and –0.45 ± 0.66 deg, respectively, which is smaller in amplitude than expected. The sources of error are then carefully considered. For both local and global analyses, the uncertainties on our results are reliably constrained by computing the standard deviation of the estimates of the gravimetric factor amplitude and phase for increasing length of the time-series. Constraints on the CW anelastic dissipation can be set since any departure of the gravimetric factor from its elastic value may provide some insights into the dissipative processes that occur at the CW period. In particular, assuming given rheological models for the Earth's mantle enables us to make the link between the gravimetric factor phase and the CW quality factor.

Beschreibung: Palaeomagnetic data obtained from archaeological materials are used for reconstructions of the Earth's magnetic field of the past millennia. While many studies tested the reliability of this recorder for palaeointensity only a few studies did this for direction. The study presents an archaeomagnetic and rock magnetic investigation applied to an experimental pottery kiln, which was operated in 2003 to produce stone ware. This kind of high-quality pottery needs a temperature of at least 1160 °C. Shortly before heating of the kiln direct absolute measurements of the absolute geomagnetic field vector have been carried out close to it. After cooling of the kiln 24 oriented palaeomagnetic samples have been taken. Although Curie temperatures are about 580 °C, that is the typical temperature for magnetite, thermal as well as alternating field demagnetisations reveal also a considerable amount of hematite as magnetic carrier. This mixture of magnetite and hematite is dominated by pseudo-single domain grains. Demagnetisation removed in some cases weak secondary components, but in most cases the specimens carried a single component thermoremanent magnetisation. The mean characteristic remanent magnetisation direction agrees on 95 per cent confidence level with the directly measured field direction. Archaeointensity was obtained from five specimens with the Thellier–Coe method and with the multiple-specimen palaeointensity domain-state corrected method. Six of these specimens also provided a result of the Dekkers–Böhnel method, which overestimated the archaeointensity by about 9 per cent compared to the direct value, while after correction for fraction the value agrees very well. For the multiple-specimen palaeointensity domain-state corrected method only fractions between 25 and 75 per cent have been used and specimens showing alteration have been excluded. Above 450 °C many specimens showed alteration of the magnetic grains. Because median destructive temperatures were often above this value in most cases the fraction was less than 50 per cent. Nevertheless the obtained intensity (48.48 ± 0.24 μ) is on 95 per cent confidence level in agreement with the direct observation. Behaviour of the specimens during the Thellier-experiments was not ideal because of narrow unblocking temperature spectra and alteration. Nevertheless, the obtained mean archaeointensity is also in agreement with the direct field observation. Here the relative palaeointensity error is about 6 per cent and very high compared the multiple-specimen palaeointensity domain-state corrected method. The investigation demonstrates that a pottery kiln can provide a very precise estimate of the ancient geomagnetic field vector.

Beschreibung: Analysing independent 1-yr data sets of 10 European superconducting gravimeters (SG) reveals statistically significant temporal variations of M2 tidal parameters. Both common short-term (〈2 yr) and long-term (〉2 yr) features are identified in all SG time-series but one. The averaged variations of the amplitude factor are about 0.2. The path of load vector variations equivalent to the temporal changes of tidal parameters suggests the presence of an 8.85 yr modulation (lunar perigee). The tidal waves having the potential to modulate M2 with this period belong to the 3rd degree constituents. Their amplitude factors turn out to be much closer to body tide model predictions than that of the main 2nd degree M2, which indicates ocean loading for 3rd degree waves to be less prominent than for 2nd degree waves within the M2 group. These two different responses to the loading suggest that the observed modulation is more due to insufficient frequency resolution of limited time-series rather than to time variable loading. Presently, SG gravity time-series are still too short to prove if time variable loading processes are involved too as in case of the annual M2 modulation known to appear for analysis intervals of less than 1 yr. Whatever the variations are caused by, they provide the upper accuracy limit for earth model validation and permit estimating the temporal stability of SG scale factors and assessing the quality of gravity time-series.

Beschreibung: A sequence of large earthquakes occurred along the North Anatolian fault in the 20th century. These earthquakes, including the 1999 Izmit/Düzce earthquakes, generally propagated westward towards the Marmara Sea, defining the Main Marmara fault as a potential seismic gap. It is important to conduct a detailed assessment of the seismic hazards along the main Marmara fault because the megacity Istanbul lies only approximately 10 km north of the eastern segment of the Main Marmara fault, which is referred to as the Princes’ Islands Fault segment (PIF). Here, we study the locking status of this fault segment to evaluate the seismic hazard potential. For the first time, combined ascending and descending Interferometric Synthetic Aperture Radar and Global Positioning System observations were used to investigate the crustal deformation associated with the PIF. After careful corrections of the estimated ground velocity, a deformation pattern relating to fault locking near the Princes’ Islands was identified. The modeling results revealed that the slip rate and locking depth of the fault segment show a clear trade-off, which were estimated as 18.9 ± 7.2 mm yr –1 and 12.1 ± 7.0 km, respectively. With a moment accumulation rate of 1.7 ± 0.4  x  10 17 Nm yr –1 (proportional to the product of slip rate and locking depth), our results imply a build-up of a geodetic moment on the PIF and therefore a potential for earthquake hazards in the vicinity of the Istanbul megacity.

Beschreibung: Deformation analysis in general and strain analysis in particular using permanent GPS networks require proper analysis of time-series in which all functional effects are taken into consideration and all stochastic effects are captured using an appropriate noise model. This contribution addresses both issues when considering the strain parameters of a GPS network. Estimates of spatial correlation, time correlated noise, and multivariate power spectrum for daily position time-series of the Southern California Integrated GPS Network (SCIGN) stations collected between 1996 and 2011 are obtained. Significant signals with periods of 13.63 d and those related to the GPS draconitic year are identified in these time-series. We aim to assess the effect of a realistic noise model of the series on the uncertainties of the strain parameters including displacements, normal and shear strains, and rotations. For the SCIGN network considered, the following results are highlighted. Contrary to the common belief, the uncertainties of the displacements parameters become smaller when taking a realistic noise model into account. This however was not the case when assessing the noise characteristics of the normal and shear strain, and rotation parameters. The uncertainties increase nearly by a factor of two, in agreement to what is expected. Some of the significant deformation parameters of the white noise model become less significant in case of the realistic noise model.

Beschreibung: We have extended backwards from 2001 to 1979 the current release 05 (RL05) of the Gravity Recovery and Climate Experiment (GRACE) Atmospheric and Oceanic De-aliasing Level-1B (AOD1B) product and studied the impact of this and a previous release 04 (RL04) of the AOD1B product on precise orbits of five altimetry satellites (ERS-1, ERS-2, TOPEX/Poseidon, Envisat and Jason-1) for the time span 1991–2012, as compared to the case when no AOD1B product is used. We have found that using AOD1B RL05 product reduces root mean square (RMS) fits of satellite laser ranging (SLR) observations by about 1.0–6.4 per cent, 2-d arc overlaps in radial, cross-track and along-track directions by about 1.3–12.0, 0.3–10.0 and 2.0–10.0 per cent, respectively, for various satellites tested, as compared to the case without AOD1B product. Using AOD1B RL05 product instead of RL04 one reduces SLR RMS fits by 0.1–0.7 per cent, 2-d arc overlaps in radial, cross-track and along-track directions by 0.1–0.6, 0.1–1.3 and 0.2–1.2 per cent, respectively, for the satellite orbits tested. The multi-mission crossover analysis shows that the application of an AOD1B product reduces the scatter of radial errors by 0.4–2.8 per cent for the satellite missions studied. At the regions with the most pronounced changes the use of the AOD1B products improves the consistency between the sea level as measured by the TOPEX and ERS-2 missions and by the Jason-1 and Envisat missions by 5 to 10 per cent (globally by about 2 per cent). The results of our study show that extended AOD1B RL05 product performs better than the AOD1B RL04 and improves orbits of altimetry satellites and consistency of sea level products.

Beschreibung: In this work, we carry out numerical simulations of the seismo-electromagnetic fields associated with a fault in a porous medium by considering the electrokinetic effect. In addition to porous materials, the adopted layered models comprise solid materials in which the electrokinetic effect is inoperative. First, sensitivity study is performed for the evanescent and direct radiation electromagnetic (EM) waves generated by a double couple point source embedded in a porous half-space below a solid half-space. Results suggest that both the evanescent and direct radiation EM waves are sensitive to some medium properties, for example porosity, salinity, fluid viscosity, and conductivity of solid layer. Then, adopting an eight-layer half-space model, we simulate the seismic and EM wavefields generated by the rupture process of a finite fault. It is shown that the electrokinetic effect is able to generate observable corupture EM signals, but the observability depends on some factors such as the epicentral distance, properties of the medium where the fault is located, and local activity levels of EM noise. Time synchronization coseismic EM signals are recorded when the receiver is close to the ground water level but located in a solid medium. In addition to the post-seismic electric field, our results also show the post-seismic magnetic field which has not been identified in previous simulation studies on the electrokinetic effect. The generation of the post-seismic magnetic field probably requires a sufficiently strong medium heterogeneity or fluid-pressure gradient.

Beschreibung: A new analytical method for the computation of a truncated series of solid spherical harmonic coefficients (HCs) from data on a spheroid (i.e. an oblate ellipsoid of revolution) is derived, using a transformation between surface and solid spherical HCs. A two-step procedure is derived to extend this transformation beyond degree and order (d/o) 520. The method is compared to the Hotine–Jekeli transformation in a numerical study based on the EGM2008 global gravity model. Both methods are shown to achieve submicrometre precision in terms of height anomalies for a model to d/o 2239. However, both methods result in spherical harmonic models that are different by up to 7.6 mm in height anomalies and 2.5 mGal in gravity disturbances due to the different coordinate system used. While the Hotine–Jekeli transformation requires the use of an ellipsoidal coordinate system, the new method uses only spherical polar coordinates. The Hotine–Jekeli transformation is numerically more efficient, but the new method can more easily be extended to cases where (a linear combination of) normal derivatives of the function under consideration are given on the surface of the spheroid. It therefore provides a solution to many types of ellipsoidal boundary-value problems in the spectral domain.

Beschreibung: Fault ruptures in the Earth's crust generate both elastic and electromagnetic (EM) waves. If the corresponding EM signals can be observed, then earthquakes could be detected before the first seismic waves arrive. In this study, I consider the piezomagnetic effect as a mechanism that converts elastic waves to EM energy, and I derive analytical formulas for the conversion process. The situation considered in this study is a whole-space model, in which elastic and EM properties are uniform and isotropic. In this situation, the governing equations of the elastic and EM fields, combined with the piezomagnetic constitutive law, can be solved analytically in the time domain by ignoring the displacement current term. Using the derived formulas, numerical examples are investigated, and the corresponding characteristics of the expected magnetic signals are resolved. I show that temporal variations in the magnetic field depend strongly on the electrical conductivity of the medium, meaning that precise detection of signals generated by the piezomagnetic effect is generally difficult. Expected amplitudes of piezomagnetic signals are estimated to be no larger than 0.3 nT for earthquakes with a moment magnitude of ≥7.0 at a source distance of 25 km; however, this conclusion may not extend to the detection of real earthquakes, because piezomagnetic stress sensitivity is currently poorly constrained.

Beschreibung: We explore Earth's elastic deformation response to ocean tidal loading (OTL) using kinematic Global Positioning System (GPS) observations and forward-modelled predictions across South America. Harmonic coefficients are extracted from up to 14 yr of GPS-inferred receiver locations, which we estimate at 5 min intervals using precise point positioning. We compare the observed OTL-induced surface displacements against predictions derived from spherically symmetric, non-rotating, elastic and isotropic (SNREI) Earth models. We also compare sets of modelled predictions directly for various ocean-tide and Earth-model combinations. The vector differences between predicted displacements computed using separate ocean-tide models reveal uniform-displacement components common to all stations in the South America network. Removal of the network-mean OTL-induced displacements from each site substantially reduces the vector differences between observed and predicted displacements. We focus on the dominant astronomical tidal harmonics from three distinct frequency bands: semidiurnal (M 2 ), diurnal (O 1 ) and fortnightly (M f ). In each band, the observed OTL-induced surface displacements strongly resemble the modelled displacement-response patterns, and the residuals agree to about 0.3 mm or better. Even with the submillimetre correspondence between observations and predictions, we detect regional-scale spatial coherency in the final set of residuals, most notably for the M 2 harmonic. The spatial coherency appears relatively insensitive to the specific choice of ocean-tide or SNREI-Earth model. Varying the load model or 1-D elastic structure yields predicted OTL-induced displacement differences of order 0.1 mm or less for the network. Furthermore, estimates of the observational uncertainty place the noise level below the magnitude of the residual displacements for most stations, supporting our interpretation that random errors cannot account for the entire misfit. Therefore, the spatially coherent residuals may reveal deficiencies in the a priori SNREI Earth models. In particular, the residuals may indicate sensitivity to regional deviations from standard globally averaged Earth structure due to the presence of the South American craton.

Beschreibung: This paper presents a methodology to sample equivalence domain (ED) in nonlinear partial differential equation (PDE)-constrained inverse problems. For this purpose, we first applied state-of-the-art stochastic optimization algorithm called Covariance Matrix Adaptation Evolution Strategy (CMAES) to identify low-misfit regions of the model space. These regions were then randomly sampled to create an ensemble of equivalent models and quantify uncertainty. CMAES is aimed at exploring model space globally and is robust on very ill-conditioned problems. We show that the number of iterations required to converge grows at a moderate rate with respect to number of unknowns and the algorithm is embarrassingly parallel. We formulated the problem by using the generalized Gaussian distribution. This enabled us to seamlessly use arbitrary norms for residual and regularization terms. We show that various regularization norms facilitate studying different classes of equivalent solutions. We further show how performance of the standard Metropolis–Hastings Markov chain Monte Carlo algorithm can be substantially improved by using information CMAES provides. This methodology was tested by using individual and joint inversions of magneotelluric, controlled-source electromagnetic (EM) and global EM induction data.

Beschreibung: Planar faults are widely adopted during inversions to determine slip distributions and fault geometries using geodetic observations; however, little research has been conducted with respect to curved faults. We attribute this to the lack of an appropriate parameterized modelling method. In this paper, we present a curved-fault modelling method (CFMM) that describes a curved fault according to specific parameters, and we also develop a corresponding hybrid iterative inversion algorithm (HIIA) to perform inversions for parametric curved-fault geometries and slips. The results of the strike-component and dip-component synthetic tests show that a complex S-shaped fault surface and a circular slip distribution are successfully recovered, indicating the strong performance of the CFMM and HIIA methods. In addition, we describe and verify a scenario for determining the number of necessary geometrical parameters for the HIIA and examine the case study of the Wenchuan earthquake, which occurred on a complex listric fault surface. During the iteration process of the HIIA, both the fault geometry and slip distribution of the Beichuan and Pengguan faults converge to optimal values, indicating a Beichuan fault (BCF) model with a continuous listric shape and gradual steepening from the southwest to the northeast, which is highly consistent with geological survey results. Both the synthetic and real-world case studies show that the HIIA and the CMFF are superior to the conventional fault modelling method based on rectangular planes and that these models have the potential for use in more integrated research involving inversion studies, such as joint slip/curved-fault-geometry inversions that take into account data resolving power.

Beschreibung: Recent high-resolution aeromagnetic surveys in South Norway have revealed numerous remanent anomalies over Mesoproterozoic metamorphic rocks. Studies on the nature of the minerals that are the remanent carriers has led to discoveries of titanohematite samples with unusual magnetic properties caused by nanoscale exsolution lamellae with their related lamellar magnetism. Here we focus on a rock unit dominated by quartz-plagioclase-biotite granulite containing titanohematite grains with a strong lattice-preferred orientation parallel to regional foliation. When samples with their natural remanent magnetization (NRM), acquired nearly 1 billion years ago, are cooled to 10 K and hysteresis loops measured, these loops show bi-modal exchange bias caused by the magnetism induced within the ilmenite by antiferromagnetic coupling with the adjacent lamellar NRM. By contrast when the samples are cooled in a strong magnetic field (1.5 Tesla), this results in unimodal lamellar magnetism, and, below the T N of ilmenite it adopts a consistent negative orientation, giving rise to unimodal negative exchange bias of 〉500 mT. The results presented here cover the chemical and magnetic properties, Mossbauer results and transmission electron microscopy of the titanohematite and ilmenite lamellae. Initial magnetic experiments indicated the shifts found in the exchange-bias experiments were directly related to the orientation of the sample to the applied field and the initial state of the NRM. In most samples with these unusual magnetic properties, ilmenite lamellae could not be seen in an optical or a scanning electron microscope. However magnetic experiments gave proof of the presence of ilmenite, later confirmed by Mössbauer spectroscopy. Several attempts were made to identify ilmenite in TEM studies, finally successful in showing ilmenite lamellae parallel to (001) of hematite with thicknesses ~1.2 to 1.7 nm and aspect ratios 7–13. Here we compare new TEM images and the magnetic behaviour of these samples to the MOD2 samples that previously showed extraordinary exchange bias properties, and investigate further the nature of these magnetic minerals.

Beschreibung: Globally gridded estimates of monthly-mean anomalies of terrestrial water storage (TWS) are estimated from the most recent GRACE release 05a of GFZ Potsdam in order to provide non-geodetic users a convenient access to state-of-the-art GRACE monitoring data. We use an ensemble of five global land model simulations with different physics and different atmospheric forcing to obtain reliable gridded scaling factors required to correct for spatial leakage introduced during data processing. To allow for the application of this data-set for large-scale monitoring tasks, model validation efforts, and subsequently also data assimilation experiments, globally gridded estimates of TWS uncertainties that include (i) measurement, (ii) leakage and (iii) re-scaling errors are provided as well. The results are generally consistent with the gridded data provided by Tellus, but deviate in some basins which are largely affected by the uncertainties of the model information required for re-scaling, where the approach based on the median of a small ensemble of global land models introduced in this paper leads to more robust results.

Beschreibung: Rock-magnetic and geochemical characteristics of three Vertisol profiles with different degree of textural differentiation have been studied. Thermomagnetic analyses, thermal demagnetization of laboratory remanences and acquisition of isothermal remanence curves are applied for identification of iron oxide mineralogy. The main magnetic minerals in Vertisols are ferrihydrite, single-domain magnetite, maghemite and hematite. Variations in magnetic susceptibility, anhysteretic remanent magnetization, isothermal remanent magnetization, as well as different ratios (Xarm/X, ARM/SIRM, S-ratio) along depth are studied. Concentration of magnetic minerals in Vertisols is low, influenced by the intense reductomorphic processes. The lowest magnetic susceptibility is found in the most texturally differentiated soil. However, rock-magnetic data suggest the presence of small, but well defined fraction of single domain-like magnetite with relatively wide grain-size distribution found in those parts of the profiles, which are subjected to most intense and frequent seasonal changes in oxidation-reduction conditions. It is suggested that this fraction is formed as a result of transformations of ferrihydrite under repeated cycles of anaerobic/aerobic conditions. Based on geochemical data, CALMAG weathering index was calculated for the three Vertisols. Using the established relation between CALMAG and mean annual precipitation (MAP), palaeo-MAP was evaluated for the studied profiles. The obtained MAP estimations fall in the range 1000–1200 mm and are much higher compared to contemporary precipitation in the area (MAP in the interval 540–770 mm). This finding confirms the relict character of Vertisols on Bulgarian territory and gives more information about the palaeoclimate during the initial stages of Vertisol formation.

Beschreibung: Regional recovery of the disturbing gravitational potential in the area of Central Europe from satellite gravitational gradients data is discussed in this contribution. The disturbing gravitational potential is obtained by inverting surface integral formulas which transform the disturbing gravitational potential onto disturbing gravitational gradients in the spherical local north-oriented frame. Two numerical approaches that solve the inverse problem are considered. In the first approach, the integral formulas are rigorously decomposed into two parts, that is, the effects of the gradient data within near and distant zones. While the effect of the near zone data is sought as an inverse problem, the effect of the distant zone data is synthesized from the global gravitational model GGM05S using spectral weights given by truncation error coefficients up to the degree 150. In the second approach, a reference gravitational field up to the degree 180 is applied to reduce and smooth measured gravitational gradients. In both cases we recovered the disturbing gravitational potential from each of the four well-measured gravitational gradients of the GOCE satellite separately as well as from their combination. Obtained results are compared with the EGM2008, DIR-r2, TIM-r2 and SPW-r2 global gravitational models. The best fit was achieved for EGM2008 and the second approach combining all four well-measured gravitational gradients with rms of 1.231 m 2  s –2 .

Beschreibung: We present efficient Fourier-domain algorithms for modelling gravity effects due to topographic masses. The well-known Parker's formula originally based on the standard fast Fourier transform (FFT) algorithm is modified by applying the Gauss–FFT method instead. Numerical precision of the forward and inverse Fourier transforms embedded in Parker's formula and its extended forms are significantly improved by the Gauss–FFT method. The topographic model is composed of two major aspects, the geometry and the density. Versatile geometric representations, including the mass line model, the mass prism model, the polyhedron model and smoother topographic models interpolated from discrete data sets using high-order splines or pre-defined by analytical functions, in combination with density distributions that vary both laterally and vertically in rather arbitrary ways following exponential or general polynomial functions, now can be treated in a consistent framework by applying the Gauss–FFT method. The method presented has been numerically checked by space-domain analytical and hybrid analytical/numerical solutions already established in the literature. Synthetic and real model tests show that both the Gauss–FFT method and the standard FFT method run much faster than space-domain solutions, with the Gauss–FFT method being superior in numerical accuracy. When truncation errors are negligible, the Gauss–FFT method can provide forward results almost identical to space-domain analytical or semi-numerical solutions in much less time.

Beschreibung: The gravity gradient tensor (GGT) has been increasingly used in practical applications, but the advantages and the disadvantages of the analysis of GGT components versus the analysis of the vertical component of the gravity field are still debated. We analyse the performance of joint inversion of GGT components versus separate inversion of the gravity field alone, or of one tensor component. We perform our analysis by inspection of the Picard Plot, a Singular Value Decomposition tool, and analyse both synthetic data and gradiometer measurements carried out at the Vredefort structure, South Africa. We show that the main factors controlling the reliability of the inversion are algebraic ambiguity (the difference between the number of unknowns and the number of available data points) and signal-to-noise ratio. Provided that algebraic ambiguity is kept low and the noise level is small enough so that a sufficient number of SVD components can be included in the regularized solution, we find that: (i) the choice of tensor components involved in the inversion is not crucial to the overall reliability of the reconstructions; (ii) GGT inversion can yield the same resolution as inversion with a denser distribution of gravity data points, but with the advantage of using fewer measurement stations.

Beschreibung: This paper compares GRACE (Gravity Recovery and Climate Experiment) and ICESat (Ice, Cloud and land Elevation Satellite) observations to confirm whether the observed gravity increase in the Tibetan Plateau (TP) was primarily caused by lake storage gain, and comprehensively analyses the changing pattern of lake level over 2003–2009. An improved automated method was used to obtain lake-level changes and the underestimation of lake water storage was considered due to lake area expansion and lake density. The result demonstrates that GRACE recorded a mass gain (16.43 ± 1.65/11.79 ± 1.25 gt a –1 ) in the total/inner TP, of which lake storage increase accounts for (8.78 ± 0.75/7.53 ± 0.56 gt a –1 ) based on ICESat. The northwestern residual may be stored in new lakes and soil moisture as a result of net precipitation gain. According to the character of the lake-level changes, we divide the TP into four subregions. Generally, the changing pattern of lake level concurs with the distribution of precipitation, which is increasing in the inner TP and decreasing in the upstream area of the Indus and Brahmaputra Rivers. An excess of rainfall in the northeastern TP in the summer of 2005 and 2009 caused a simultaneous large increase in water level in many lakes. The correlation of lake changes with precipitation demonstrates that precipitation rather than glacial melt is the main cause of lake-level change in most places. Nonetheless, the meltwater is a considerable supplement for lakes near glaciers such as Selin Co and Nam Co, which partly explains why GRACE indicates a much weaker signal in this region.

Beschreibung: This work presents MARE2DEM, a freely available code for 2-D anisotropic inversion of magnetotelluric (MT) data and frequency-domain controlled-source electromagnetic (CSEM) data from onshore and offshore surveys. MARE2DEM parametrizes the inverse model using a grid of arbitrarily shaped polygons, where unstructured triangular or quadrilateral grids are typically used due to their ease of construction. Unstructured grids provide significantly more geometric flexibility and parameter efficiency than the structured rectangular grids commonly used by most other inversion codes. Transmitter and receiver components located on topographic slopes can be tilted parallel to the boundary so that the simulated electromagnetic fields accurately reproduce the real survey geometry. The forward solution is implemented with a goal-oriented adaptive finite-element method that automatically generates and refines unstructured triangular element grids that conform to the inversion parameter grid, ensuring accurate responses as the model conductivity changes. This dual-grid approach is significantly more efficient than the conventional use of a single grid for both the forward and inverse meshes since the more detailed finite-element meshes required for accurate responses do not increase the memory requirements of the inverse problem. Forward solutions are computed in parallel with a highly efficient scaling by partitioning the data into smaller independent modeling tasks consisting of subsets of the input frequencies, transmitters and receivers. Non-linear inversion is carried out with a new Occam inversion approach that requires fewer forward calls. Dense matrix operations are optimized for memory and parallel scalability using the ScaLAPACK parallel library. Free parameters can be bounded using a new non-linear transformation that leaves the transformed parameters nearly the same as the original parameters within the bounds, thereby reducing non-linear smoothing effects. Data balancing normalization weights for the joint inversion of two or more data sets encourages the inversion to fit each data type equally well. A synthetic joint inversion of marine CSEM and MT data illustrates the algorithm's performance and parallel scaling on up to 480 processing cores. CSEM inversion of data from the Middle America Trench offshore Nicaragua demonstrates a real world application. The source code and MATLAB interface tools are freely available at http://mare2dem.ucsd.edu .

Beschreibung: Traditional processing of Global Navigation Satellite System (GNSS) data using dedicated scientific software has provided the highest levels of positional accuracy, and has been used extensively in geophysical deformation studies. To achieve these accuracies a significant level of understanding and training is required, limiting their availability to the general scientific community. Various online GNSS processing services, now freely available, address some of these difficulties and allow users to easily process their own GNSS data and potentially obtain high quality results. Previous research into these services has focused on Continually Operating Reference Station (CORS) GNSS data. Less research exists on the results achievable with these services using large campaign GNSS data sets, which are inherently noisier than CORS data. Even less research exists on the quality of velocity fields derived from campaign GNSS data processed through online precise point positioning services. Particularly, whether they are suitable for geodynamic and deformation studies where precise and reliable velocities are needed. In this research, we process a very large campaign GPS data set (spanning 10 yr) with the online Jet Propulsion Laboratory Automated Precise Positioning Service. This data set is taken from a GNSS network specifically designed and surveyed to measure deformation through the central North Island of New Zealand. This includes regional CORS stations. We then use these coordinates to derive a horizontal and vertical velocity field. This is the first time that a large campaign GPS data set has been processed solely using an online service and the solutions used to determine a horizontal and vertical velocity field. We compared this velocity field to that of another well utilized GNSS scientific software package. The results show a good agreement between the CORS positions and campaign station velocities obtained from the two approaches. We discuss the implications of these results for how future GNSS campaign field surveys might be conducted and how their data might be processed.

Beschreibung: The spectral complex conductivity of a water-bearing sand during interaction with carbon dioxide (CO 2 ) is influenced by multiple, simultaneous processes. These processes include partial saturation due to the replacement of conductive pore water with CO 2 and chemical interaction of the reactive CO 2 with the bulk fluid and the grain-water interface. We present a laboratory study on the spectral induced polarization of water-bearing sands during exposure to and flow-through by CO 2 . Conductivity spectra were measured successfully at pressures up to 30 MPa and 80 °C during active flow and at steady-state conditions concentrating on the frequency range between 0.0014 and 100 Hz. The frequency range between 0.1 and 100 Hz turned out to be most indicative for potential monitoring applications. The presented data show that the impact of CO 2 on the electrolytic conductivity may be covered by a model for pore-water conductivity, which depends on salinity, pressure and temperature and has been derived from earlier investigations of the pore-water phase. The new data covering the three-phase system CO 2 -brine-sand further show that chemical interaction causes a reduction of surface conductivity by almost 20 per cent, which could be related to the low pH-value in the acidic environment due to CO 2 dissolution and the dissociation of carbonic acid. The quantification of the total CO 2 effect may be used as a correction during monitoring of a sequestration in terms of saturation. We show that this leads to a correct reconstruction of fluid saturation from electrical measurements. In addition, an indicator for changes of the inner surface area, which is related to mineral dissolution or precipitation processes, can be computed from the imaginary part of conductivity. The low frequency range between 0.0014 and 0.1 Hz shows additional characteristics, which deviate from the behaviour at higher frequencies. A Debye decomposition approach is applied to isolate the feature dominating the data at low frequencies. We conclude from our study that electrical conductivity is not only a highly sensitive indicator for CO 2 saturation in pore space. When it is measured in its full spectral and complex form it contains additional information on the chemical state of the system, which holds the potential of getting access to both saturation and interface properties with one monitoring method.

Beschreibung: We have developed a method to estimate the geometry, location and densities of anomalies coming from 2-D gravity data based on compact gravity inversion technique. Compact gravity inversion is simple, fast and user friendly but severely depends on the number of model parameters, that is, by increasing the model parameters, the anomalies tend to concentrate near the surface. To overcome this ambiguity new weighting functions based on density contrast, depth, and compactness models have been introduced. Variable compactness factors have been defined here to get either a sharp or a smooth model based on the depth of the source or existence of prior information. Depth weighting derived from one station of gravity data whereas the effect of gravity data is 2-D and 3-D. To compensate this limitation an innovating weighting function namely kernel function has been introduced which multiplies with weight and compactness matrixes to yield a general model weighting function. The method is tested using three different sets of synthetic examples: a body at various depths (20, 40, 80 and 140 m), two bodies at the same depth but various distances to estimate lateral resolution and three bodies with negative and positive density contrast in different depths. The method is also applied to three real gravity data of Woodlawn massive sulphide body, sulphides mineralization of British Colombia and iron ore body of Missouri. The method produces solutions consistent with the known geologic attributes of the gravity sources, illustrating its potential practicality.

Beschreibung: Our study attempts to identify a characteristic magnetic signature of overbank sediments exhibiting anthropogenically induced magnetic enhancement and thereby to distinguish them from unenhanced sediments with weak magnetic background values, using a novel approach based on data mining methods, thus providing a mean of rapid pollution determination. Data were obtained from 539 bulk samples from vertical profiles through overbank sediment, collected on seven rivers in the eastern Czech Republic and three rivers in northwest England. k -Means clustering and hierarchical clustering methods, paired group (UPGMA) and Ward's method, were used to divide the samples to natural groups according to their attributes. Interparametric ratios: SIRM/; SIRM/ARM; and S –0.1T were chosen as attributes for analyses making the resultant model more widely applicable as magnetic concentration values can differ by two orders. Division into three clusters appeared to be optimal and corresponded to inherent clusters in the data scatter. Clustering managed to separate samples with relatively weak anthropogenically induced enhancement, relatively strong anthropogenically induced enhancement and samples lacking enhancement. To describe the clusters explicitly and thus obtain a discrete magnetic signature, classification rules (JRip method) and decision trees (J4.8 and Simple Cart methods) were used. Samples lacking anthropogenic enhancement typically exhibited an S –0.1T 〈 c. 0.5, SIRM/ARM 〈 c. 150 and SIRM/ 〈 c. 6000 A m –1 . Samples with magnetic enhancement all exhibited an S –0.1T 〉 0.5. Samples with relatively stronger anthropogenic enhancement were unequivocally distinguished from the samples with weaker enhancement by an SIRM/ARM 〉 c. 150. Samples with SIRM/ARM in a range c. 126–150 were classified as relatively strongly enhanced when their SIRM/ 〉 18 000 A m –1 and relatively less enhanced when their SIRM/ 〈 18 000 A m –1 . An additional rule was arbitrary added to exclude samples with fd% 〉 6 per cent from anthropogenically enhanced clusters as samples with natural magnetic enhancement. The characteristics of the clusters resulted mainly from the relationship between SIRM/ARM and the S –0.1T , and SIRM/ and the S –0.1T . Both SIRM/ARM and SIRM/ increase with increasing S –0.1T values reflecting a greater level of anthropogenic magnetic particles. Overall, data mining methods demonstrated good potential for utilization in environmental magnetism.

Beschreibung: The construction of geomagnetic, archaeomagnetic or palaeomagnetic field models requires some prior knowledge about the actual field, which can be gathered from the statistical properties of the field over a variety of length-scales and timescales. However, available geomagnetic data on centennial to millennial periods are too sparse to infer directly these statistical properties. We thus use high-resolution numerical simulations of the geodynamo to test a method for estimating the temporal power spectra (or equivalently the autocovariance functions) of the individual Gauss coefficients that describe the geomagnetic field outside the Earth's fluid outer core. Based on the spectral analysis of our simulations, we argue that a prior for the observational geomagnetic field over decennial to millennial periods can be constructed from the statistics of the field during the short satellite era. The method rests on the assumption that time-series of spherical harmonic coefficients can be considered as realizations of stationary and differentiable stochastic processes, namely order 2 autoregressive (AR2) processes. In the framework of these processes, the statistics of Gauss coefficients are well constrained by their variance and one or two timescales. We find that the time spectra in the dynamo simulations of all Gauss coefficients but the axial dipole are well approximated by the spectra of AR2 processes characterized by only one timescale. The process parameters can simply be deduced from instantaneous estimates of the spatial power spectra of the magnetic field and of its first time derivative. Some deviations of the Gauss coefficients statistics from this minimal model are also discussed. Characterizing the axial dipole clearly requires a more sophisticated AR2 process, with a second distinct timescale.

Beschreibung: This paper resurrects a version of Poisson's Partial Differential Equation (PDE) associated with the gravitational field at the Earth's surface and illustrates how the PDE possesses a capability to extract the mass density of Earth's topography from land-based gravity data. Herein, first we propound a theorem which mathematically introduces this version of Poisson's PDE adapted for the Earth's surface and then we use this PDE to develop a method of approximating the terrain mass density. Also, we carry out a real case study showing how the proposed approach is able to be applied to a set of land-based gravity data. In the case study, the method is summarized by an algorithm and applied to a set of gravity stations located along a part of the north coast of the Persian Gulf in the south of Iran. The results were numerically validated via rock-samplings as well as a geological map. Also, the method was compared with two conventional methods of mass density reduction. The numerical experiments indicate that the Poisson PDE at the Earth's surface has the capability to extract the mass density from land-based gravity data and is able to provide an alternative and somewhat more precise method of estimating the terrain mass density.

Beschreibung: The Luonan Basin, located in the transitional zone between temperate and subtropical China, is an important locality for human evolution during the early to middle Pleistocene. The loess deposits in the Luonan Basin contain numerous in situ lithic artefacts; the deposits also constitute suitable material for dating the artefacts and are potentially useful for reconstructing the climatic fluctuations which is important for studying the adaptation and occupation of the area by early humans. We carried out a combined rock magnetic and geochemical investigation of a loess sequence from the Liuwan Palaeolithic site in the Luonan Basin. The results indicate a mixture of magnetic minerals, including magnetite/maghemite and hematite/goethite. Magnetic susceptibility was used as a palaeoclimate proxy on the Chinese Loess Plateau; however, its application to the Luonan Basin may be problematic because the provenance of the loess parent material, as well as the depositional environment, differs from that of the Chinese Loess Plateau. We found that rock magnetic parameters related to the grain size of magnetic minerals, such as SIRM/ and ARM /SIRM, are better palaeoclimatic indicators than magnetic susceptibility. Overall, the magnetic results, together with the results of bulk grain-size and chemical index of alteration, indicate that the interglacial environment of early humans in Luonan Basin was warmer and more humid than the coeval environment of the Chinese Loess Plateau.

Beschreibung: The Moho surface can be determined according to isostatic theories, and among them, the recent Vening Meinesz-Moritz (VMM) theory of isostasy has been successfully applied for this purpose. In this paper, this method is studied from a theoretical prospective and its connection to the Airy–Heiskanen (AH) and Vening Meinesz original theories are presented. Jeffrey's inverse solution to isostasy is developed according to the recent developments of the VMM method and both are compared in similar situations. It is shown that they are generalizations of the AH model in a global and continuous domain. In the VMM spherical harmonic solution for Moho depth, the mean Moho depth contributes only to the zero-degree term of the series, while in Jeffrey's solution it contributes to all frequencies. In addition, the VMM spherical harmonic series is improved further so that the mean Moho can contribute to all frequencies of the solution. This modification makes the VMM global solution superior to the Jeffrey one, but in a global scale, the difference between both solutions is less than 3 km. Both solutions are asymptotically convergent and we present two methods to obtain smooth solutions for Moho from them.

Beschreibung: We develop the theory for the magnetic to magnetic and kinetic to magnetic energy transfer between different spherical harmonic degrees due to the interaction of fluid flow and radial magnetic field at the top of the Earth's core. We show that non-zero secular variation of the total magnetic energy could be significant and may provide evidence for the existence of stretching secular variation, which suggests the existence of radial motions at the top of the Earth's core—whole core convection or MAC waves. However, the uncertainties of the small scales of the geomagnetic field prevent a definite conclusion. Combining core field and flow models we calculate the detailed magnetic to magnetic and kinetic to magnetic energy transfer matrices. The magnetic to magnetic energy transfer shows a complex behaviour with local and non-local transfers. The spectra of magnetic to magnetic energy transfers show clear maxima and minima, suggesting an energy cascade. The kinetic to magnetic energy transfers, which are much weaker due to the weak poloidal flow, are either local or non-local between degree one and higher degrees. The patterns observed in the matrices resemble energy transfer patterns that are typically found in 3-D MHD numerical simulations.

Beschreibung: In a pioneering study, Wahr & Bergen developed the widely adopted, pseudo-normal mode framework for predicting the impact of anelastic effects on the Earth's body tides. Lau et al. have recently derived an extended normal mode treatment of the problem (as well as a minor variant of the theory known as the direct solution method) that makes full use of theoretical developments in free oscillation seismology spanning the last quarter century and that avoids a series of assumptions and approximations adopted in the traditional theory for predicting anelastic effects. There are two noteworthy differences between these two theories: (1) the traditional theory only considers perturbations to the eigenmodes of an elastic Earth, whereas the new theory augments this set of modes to include the relaxation modes that arise in anelastic behaviour; and (2) the traditional theory approximates the complex perturbation to the tidal Love number as a scaled version of the complex perturbation to the elastic moduli, whereas the new theory computes the full complex perturbation to each eigenmode. In this study, we highlight the above differences using a series of synthetic calculations, and demonstrate that the traditional theory can introduce significant error in predictions of the complex perturbation to the Love numbers due to anelasticity and the related predictions of tidal lag angles. For the simplified Earth models we adopt, the computed lag angles differ by ~20 per cent. The assumptions in the traditional theory have important implications for previous studies that use model predictions to correct observables for body tide signals or that analyse observations of body tide deformation to infer mantle anelastic structure. Finally, we also highlight the fundamental difference between apparent attenuation (i.e. attenuation inferred from observations or predicted using the above theories) and intrinsic attenuation (i.e. the material property investigated through experiments), where both are often expressed in terms of lag angles or Q –1 . In particular, we demonstrate the potentially significant (factor of two or more) bias introduced in estimates of Q –1 and its frequency dependence in studies that have treated Q –1 determined from tidal phase lags or measured experimentally as being equal. The observed or theoretically predicted lag angle (or apparent Q –1 ) differs from the intrinsic, material property due to inertia, self-gravity and effects associated with the energy budget. By accounting for these differences we derive, for a special case, an expression that accurately maps apparent attenuation predicted using the extended normal mode formalism of Lau et al. into intrinsic attenuation. The theory allows for more generalized mappings which may be used to robustly connect observations and predictions of tidal lag angles to results from laboratory experiments of mantle materials.

Beschreibung: The magnetic susceptibility measured in alternating field can in general be resolved into a component that is in-phase with the applied field and a component that is out-of-phase. While in non-conductive diamagnetic, paramagnetic and many ferromagnetic materials the phase is effectively zero, in some ferromagnetic minerals, such as pyrrhotite, hematite, titanomagnetite or small magnetically viscous grains of magnetite, it is clearly non-zero. The anisotropy of out-of-phase susceptibility (opAMS) can then be used as a tool for the direct determination of the magnetic subfabrics of the minerals with non-zero phase. The error in determination of out-of-phase susceptibility non-linearly increases with decreasing phase angle. This may result in imprecise determination of the opAMS in specimens with very low phase angle. The degree of opAMS is higher than that of ipAMS, which may in contrast result in slightly increasing precision n the opAMS determination. It is highly recommended to inspect the results of the statistical tests of each specimen and to exclude the specimens whose opAMS is determined with insufficient precision from further processing. In rocks, whose magnetism is dominated by the mineral with non-zero out-of-phase susceptibility, the principal directions of the opAMS and ipAMS are virtually coaxial, while the degree of opAMS is higher than that of ipAMS. In some specific cases, the opAMS provides us with similar data to those provided by anisotropies of low-field dependent susceptibility and frequency-dependent susceptibility. The advantage of the opAMS compared to the other two anisotropies is its simultaneous measurement with the ipAMS during one measuring process, while the other two anisotropies require the AMS measurements in several fields or at least at two operating frequencies.

Beschreibung: The temperature dependence of magnetic susceptibility ( –T curves) and magnetization ( M–T curves) has been used as a routine rock magnetic tool to characterize the magnetic mineralogy and magnetic granulometry of Chinese loess/palaeosols. However, palaeoclimatic interpretation of these thermomagnetic analyses remains controversial. In the present study, total organic carbon (TOC), thermomagnetic and low-temperature magnetic experiments on Mid-Late Pleistocene loess/palaeosols in Central Asia and the Chinese Loess Plateau (CLP) have been conducted. We found that the M ( T ) cooling curves at room temperature were mostly lower than the corresponding heating curves, whereas for the ( T ) analyses the cooling curves at room temperature were always much higher than the heating curves. Low-temperature magnetic measurements demonstrated that a large amount of superparamagnetic ferrimagnetic particles were produced during the thermal treatment and resulted in the aforementioned differences. This finding further indicated that the use of the M–T curves to estimate the relative content of maghemite in the loess/palaeosols from the CLP was problematic. In addition, a positive correlation exists between the TOC and the frequency-dependent susceptibility ( FD ) in the CLP, suggesting that stronger pedogenesis would result in the simultaneous increase in the content of both maghemite and organic matter. Consequently, the parameters 1 (representing the relative content of pedogenic maghemite), 2 ([ ph – ] + 1 ) and ph (related to the organic matter concentration), which can be calculated from the –T analyses, can potentially be used as new indicators of pedogenesis and palaeoclimate in Central Asia and the CLP.

Beschreibung: We consider a new approach to both the forward and inverse problems in post-seismic deformation. We present a method for forward modelling post-seismic deformation in a self-gravitating, heterogeneous and compressible earth with a variety of linear and nonlinear rheologies. We further demonstrate how the adjoint method can be applied to the inverse problem both to invert for rheological structure and to calculate the sensitivity of a given surface measurement to changes in rheology or time-dependence of the source. Both the forward and inverse aspects are illustrated with several numerical examples implemented in a spherically symmetric earth model.

Beschreibung: This is the second of three papers investigating properties of titanohematite-bearing quartzo-feldspathic rocks that create a significant remanent magnetic anomaly in the Modum District, South Norway. The first paper provided initial magnetic results, mineralogical characterization and evidence for the presence of lamellar magnetism. In this paper, knowledge of lamellar magnetic properties is explored through experiments where ilmenite lamellae were magnetized below 57 K, and interact magnetically along interfaces with the titanohematite host. Samples with known NRM directions were placed in specific orientations in an MPMS then cooled in zero field to 5 K, where hysteresis loops were measured in fields up to 5 Tesla. This assured that results were ultimately related to the natural lamellar magnetism produced during cooling ~1 billion years ago. In a second set of experiments the same oriented samples, were subjected to a +5 Tesla (T) field then field cooled to 5 K before hysteresis experiments. The first experiments consistently produced asymmetric shifted hysteresis loops with two loop separations, one in a positive field and one in a negative field. Without exception, when the NRM was oriented toward the negative field end of the MPMS, the bimodal loop showed a dominant loop separation in a positive field. By contrast, when the NRM was oriented toward the positive field end of the MPMS, the bimodal loop showed a dominant loop separation in a negative field. Both observations are consistent with antiferromagnetic coupling between the hard magnetization of ilmenite and the more easily shifted lamellar magnetism of the hematite. The bimodal nature of the loops indicates that the NRMs are vector sums of natural lamellar moments, which are oriented both positively and negatively, and that these opposite moments control the orientations of ilmenite magnetizations when cooling through 57 K. Here, extreme exchange biases up to 1.68 T were measured. The second set of experiments produced asymmetric shifted hysteresis loops with one opening always in the negative field. These observations indicate that the +5 T field applied at room temperature rotated the hematite lamellar magnetism in a positive direction, so that upon cooling all the ilmenite lamellae acquired negative magnetic moments, thus causing unimodal negatively shifted loops. Here, the largest exchange bias among the unimodal loops was only 0.7 T. These results will be used in paper III to build a better understanding of lamellar magnetism at the atomic layer scale.

Beschreibung: The Seiland Igneous Province (SIP) is the largest complex of mafic and ultramafic intrusions in northern Fennoscandia intruded at ca . 580–560 Ma. The depth extent and the deep structure of the SIP are mainly unknown apart from three profiles modelled by gravity and refraction seismic data. Utilizing 3-D gravity modelling, a complex model of the deep subsurface structure of the SIP has been developed. The structure is presented in a multiprofile model ranging from the surface to the Moho. The mafic/ultramafic rocks of the SIP are modelled with densities of 3100 and 3300 kg m –3 , the surrounding rocks by densities of 2700 and 2900 kg m –3 for upper and lower crust, respectively. This density model explains the pronounced positive Bouguer gravity anomaly of up to 100 mGal above background. Its minimum volume is estimated from the subsurface model to 17 000 km 3 and as such we revise downwards the earlier estimations of 25 000 km 3 . The new subsurface model suggests that most of the SIP has a thickness between 2 and 4 km. An area with roots in an annular pattern is found and two deep-reaching roots have been identified located below the islands of Seiland and Sørøy. The depth of these roots is estimated to approximatively 9 km. The SIP is presently interpreted to be in the Caledonian Kalak Nappe Complex and the roots depth constrains its minimum thickness which is larger than earlier estimated. Furthermore, the rather undisturbed shape of the annular root pattern indicates that the SIP has not been subjected to strong tectonic reworking during the Caledonian orogeny.

Beschreibung: Non-heating palaeointensity methods are a vital tool to explore magnetic field strength variations recorded by thermally sensitive materials of both terrestrial and extraterrestrial origin. One such method is the calibrated pseudo-Thellier method in which a specimen's natural remanent magnetization is alternating field demagnetized and replaced with a laboratory induced anhysteretic remanent magnetization (as an analogue of a thermoremanent magnetization, TRM). Using a set of 56 volcanic specimens given laboratory TRMs in fields of 10–130 μT, we refine the calibration of the pseudo-Thellier method and better define the uncertainty associated with its palaeointensity estimates. Our new calibration, obtained from 32 selected specimens, resolves the issue of non-zero intercept, which is theoretically predicted, but not satisfied by any previous calibration. The range of individual specimen calibration factors, however, is relatively large, but consistent with the variability expected for SD magnetite. We explore a number of rock magnetic parameters in an attempt to identify selection thresholds for reducing the calibration scatter, but fail to find a suitable choice. We infer that our careful selection process, which incorporates more statistics then previous studies, may be largely screening out any strong rock magnetic dependence. Some subtle grain size or mineralogical dependencies, however, remain after selection, but cannot be discerned from the scatter expected for grain size variability of SD magnetite. As a consequence of the variability in the calibration factor, the uncertainty associated with pseudo-Thellier results is much larger than previously indicated. The scatter of the calibration is ~25 per cent of the mean value, which implies that, when combined with the scatter of results typically obtained from a single site, the uncertainty of averaged pseudo-Thellier results will always be 〉25 per cent. As such, pseudo-Thellier results should be complementary to, and cross-validated with results from other methods. Nevertheless, the pseudo-Thellier method remains a valuable tool for obtaining palaeointensity estimates from thermally sensitive terrestrial and extraterrestrial materials and with careful data selection and analysis can yield results that are accurate to within a factor of 4 or better.

Beschreibung: As a consequence of measuring time variations of the electric and the magnetic field, which are related to current flow and charge distribution, magnetotelluric (MT) data in 2-D and 3-D environments are not only sensitive to the geoelectrical structures below the measuring points but also to any lateral anomalies surrounding the acquisition site. This behaviour complicates the characterization of the electrical resistivity distribution of the subsurface, particularly in complex areas. In this manuscript we assess the main advantages of complementing the standard MT impedance tensor ( Z ) data with interstation horizontal magnetic tensor ( H ) and geomagnetic transfer function ( T ) data in constraining the subsurface in a 3-D environment beneath a MT profile. Our analysis was performed using synthetic responses with added normally distributed and scattered random noise. The sensitivity of each type of data to different resistivity anomalies was evaluated, showing that the degree to which each site and each period is affected by the same anomaly depends on the type of data. A dimensionality analysis, using Z , H and T data, identified the presence of the 3-D anomalies close to the profile, suggesting a 3-D approach for recovering the electrical resistivity values of the subsurface. Finally, the capacity for recovering the geoelectrical structures of the subsurface was evaluated by performing joint inversion using different data combinations, quantifying the differences between the true synthetic model and the models from inversion process. Four main improvements were observed when performing joint inversion of Z , H and T data: (1) superior precision and accuracy at characterizing the electrical resistivity values of the anomalies below and outside the profile; (2) the potential to recover high electrical resistivity anomalies that are poorly recovered using Z data alone; (3) improvement in the characterization of the bottom and lateral boundaries of the anomalies with low electrical resistivity; and (4) superior imaging of the horizontal continuity of structures with low electrical resistivity. These advantages offer new opportunities for the MT method by making the results from a MT profile in a 3-D environment more convincing, supporting the possibility of high-resolution studies in 3-D areas without expending a large amount of economical and computational resources, and also offering better resolution of targets with high electrical resistivity.

Beschreibung: Estimating the relief of the Moho from gravity data is a computationally intensive nonlinear inverse problem. What is more, the modelling must take the Earths curvature into account when the study area is of regional scale or greater. We present a regularized nonlinear gravity inversion method that has a low computational footprint and employs a spherical Earth approximation. To achieve this, we combine the highly efficient Bott's method with smoothness regularization and a discretization of the anomalous Moho into tesseroids (spherical prisms). The computational efficiency of our method is attained by harnessing the fact that all matrices involved are sparse. The inversion results are controlled by three hyperparameters: the regularization parameter, the anomalous Moho density-contrast, and the reference Moho depth. We estimate the regularization parameter using the method of hold-out cross-validation. Additionally, we estimate the density-contrast and the reference depth using knowledge of the Moho depth at certain points. We apply the proposed method to estimate the Moho depth for the South American continent using satellite gravity data and seismological data. The final Moho model is in accordance with previous gravity-derived models and seismological data. The misfit to the gravity and seismological data is worse in the Andes and best in oceanic areas, central Brazil and Patagonia, and along the Atlantic coast. Similarly to previous results, the model suggests a thinner crust of 30–35 km under the Andean foreland basins. Discrepancies with the seismological data are greatest in the Guyana Shield, the central Solimões and Amazonas Basins, the Paraná Basin, and the Borborema province. These differences suggest the existence of crustal or mantle density anomalies that were unaccounted for during gravity data processing.

Beschreibung: Geophysical parameters of the deep Earth's interior can be evaluated through the resonance effects associated with the core and inner-core wobbles on the forced nutations of the Earth's figure axis, as observed by very long baseline interferometry (VLBI), or on the diurnal tidal waves, retrieved from the time-varying surface gravity recorded by superconducting gravimeters (SGs). In this paper, we inverse for the rotational mode parameters from both techniques to retrieve geophysical parameters of the deep Earth. We analyse surface gravity data from 15 SG stations and VLBI delays accumulated over the last 35 yr. We show existing correlations between several basic Earth parameters and then decide to inverse for the rotational modes parameters. We employ a Bayesian inversion based on the Metropolis–Hastings algorithm with a Markov-chain Monte Carlo method. We obtain estimates of the free core nutation resonant period and quality factor that are consistent for both techniques. We also attempt an inversion for the free inner-core nutation (FICN) resonant period from gravity data. The most probable solution gives a period close to the annual prograde term (or S 1 tide). However the 95 per cent confidence interval extends the possible values between roughly 28 and 725 d for gravity, and from 362 to 414 d from nutation data, depending on the prior bounds. The precisions of the estimated long-period nutation and respective small diurnal tidal constituents are hence not accurate enough for a correct determination of the FICN complex frequency.

Beschreibung: We evaluate the benefit of different global geophysical loading products on the internal scatter of GPS position time-series from 240 globally distributed sites. We focus on the non-tidal atmospheric pressure loading predicted from NASA's Modern-Era Retrospective Analysis for Research and Applications (MERRA-NATML) and the European Center for Medium-Range Weather Forecasts operational model (ECMWF-NATML), non-tidal ocean loading predicted from the Ocean Model for Circulation and Tides model (OMCT-NTOL), and the continental water storage loading predicted from the MERRA model (MERRA-CWSL) and the GFZ's Land Surface Discharge Model (LSDM-CWSL), respectively. The result shows that the root mean square (RMS) discrepancy of different CWSL models is larger than that of NATML models in the vertical component due to the varying model parameters and input data choices. We discuss the performance of different loading models and their combination to reduce the weighted RMS of GPS up-coordinates. MERRA-NATML & OMCT-NTOL & MERRA-CWSL reduced the weighted RMS (WRMS) in 96 per cent (JPL solutions) and 86 per cent (SOPAC solutions) of the cases, and ECMWF-NATML & OMCT-NTOL & LSDM-CWSL reduced the WRMS in 95 per cent (JPL solutions) and 88 per cent (SOPAC solutions) of the cases. The result shows that local effects and technical uncertainties in GPS time-series hamper the meaningful comparison between GPS observations and mass loading models. Hence, simply using the RMS of the time-series as the assessment criteria may lead to biased comparison results. Nonetheless, we give a detailed comparison (differences in phase and amplitude at seasonal timescales) for eight representative stations located adjacent to great rivers, lakes and reservoirs. We find that LSDM can provide a complementary model to study the small-scale hydrological loading like loading extremes along river channels. However, such small-scale hydrological loading effects are still instable to be modelled in some regions with its current accuracy. Finally, we discuss the impacts of mass loading corrections on the velocity and noise estimates. The noise reductions have the consistent performance as WRMS reductions for most sites, whereas some sites have their noise increased even though GPS signal WRMS is decreased there, suggesting that our posterior correction is potentially feasible, but not sufficient.

Beschreibung: Although an amount of hydrocarbon has been discovered in the West Korea Bay Basin (WKBB), located in the North Korean offshore area, geophysical investigations associated with these hydrocarbon reservoirs are not permitted because of the current geopolitical situation. Interpretation of satellite-derived potential field data can be alternatively used to image the 3-D density distribution in the sedimentary basin associated with hydrocarbon deposits. We interpreted the TRIDENT satellite-derived gravity field data to provide detailed insights into the spatial distribution of sedimentary density structures in the WKBB. We used 3-D forward density modelling for the interpretation that incorporated constraints from existing geological and geophysical information. The gravity data interpretation and the 3-D forward modelling showed that there are two modelled areas in the central subbasin that are characterized by very low density structures, with a maximum density of about 2000 kg m –3 , indicating some type of hydrocarbon reservoir. One of the anticipated hydrocarbon reservoirs is located in the southern part of the central subbasin with a volume of about 250 km 3 at a depth of about 3000 m in the Cretaceous/Jurassic layer. The other hydrocarbon reservoir should exist in the northern part of the central subbasin, with an average volume of about 300 km 3 at a depth of about 2500 m.

Beschreibung: Inference of globally averaged eustatic sea level (ESL) rise since the Last Glacial Maximum (LGM) highly depends on the interpretation of relative sea level (RSL) observations at Barbados and Bonaparte Gulf, Australia, which are sensitive to the viscosity structure of Earth's mantle. Here we examine the RSL changes at the LGM for Barbados and Bonaparte Gulf ( ${\rm{RSL}}_{\rm{L}}^{{\rm{Bar}}}$ and ${\rm{RSL}}_{\rm{L}}^{{\rm{Bon}}}$ ), differential RSL for both sites ( $\Delta {\rm{RSL}}_{\rm{L}}^{{\rm{Bar}},{\rm{Bon}}})$ and rate of change of degree-two harmonics of Earth's geopotential due to glacial isostatic adjustment (GIA) process (GIA-induced 2 ) to infer the ESL component and viscosity structure of Earth's mantle. Differential RSL, $\Delta {\rm{RSL}}_{\rm{L}}^{{\rm{Bar}},{\rm{Bon}}}$ and GIA-induced 2 are dominantly sensitive to the lower-mantle viscosity, and nearly insensitive to the upper-mantle rheological structure and GIA ice models with an ESL component of about (120–130) m. The comparison between the predicted and observationally derived $\Delta {\rm{RSL}}_{\rm{L}}^{{\rm{Bar}},{\rm{Bon}}}$ indicates the lower-mantle viscosity higher than ~2  x  10 22 Pa s, and the observationally derived GIA-induced 2 of –(6.0–6.5)  x  10 –11 yr –1 indicates two permissible solutions for the lower mantle, ~10 22 and (5–10)  x  10 22 Pa s. That is, the effective lower-mantle viscosity inferred from these two observational constraints is (5–10)  x  10 22 Pa s. The LGM RSL changes at both sites, ${\rm{RSL}}_{\rm{L}}^{{\rm{Bar}}}$ and ${\rm{RSL}}_{\rm{L}}^{{\rm{Bon}}}$ , are also sensitive to the ESL component and upper-mantle viscosity as well as the lower-mantle viscosity. The permissible upper-mantle viscosity increases with decreasing ESL component due to the sensitivity of the LGM sea level at Bonaparte Gulf ( ${\rm{RSL}}_{\rm{L}}^{{\rm{Bon}}}$ ) to the upper-mantle viscosity, and inferred upper-mantle viscosity for adopted lithospheric thicknesses of 65 and 100 km is (1–3)  x  10 20 Pa s for ESL~130 m and (4–10)  x  10 20 Pa s for ESL~125 m. The former solution of (1–3)  x  10 20 Pa s is consistent with the inferences from the postglacial differential RSL changes in the Australian region and also inversion study of far-field sea-level data. The inference of the viscosity structure based on these four observational constraints is, however, relatively insensitive to the viscosity structure of D'' layer.

Beschreibung: We present the theory, algorithms and implementation of a parallel finite-volume algorithm for the solution of the incompressible magnetohydrodynamic (MHD) equations using unstructured grids that are applicable for a wide variety of geometries. Our method implements a mixed Adams-Bashforth/Crank-Nicolson scheme for the nonlinear terms in the MHD equations and we prove that it is stable independent of the time step. To ensure that the solenoidal condition is met for the magnetic field, we use a method whereby a pseudo-pressure is introduced into the induction equation; since we are concerned with incompressible flows, the resulting Poisson equation for the pseudo-pressure is solved alongside the equivalent Poisson problem for the velocity field. We validate our code in a variety of geometries including periodic boxes, spheres, spherical shells, spheroids and ellipsoids; for the finite geometries we implement the so-called ferromagnetic or pseudo-vacuum boundary conditions appropriate for a surrounding medium with infinite magnetic permeability. This implies that the magnetic field must be purely perpendicular to the boundary. We present a number of comparisons against previous results and against analytical solutions, which verify the code's accuracy. This documents the code's reliability as a prelude to its use in more difficult problems. We finally present a new simple drifting solution for thermal convection in a spherical shell that successfully sustains a magnetic field of simple geometry. By dint of its rapid stabilization from the given initial conditions, we deem it suitable as a benchmark against which other self-consistent dynamo codes can be tested.

Beschreibung: We consider an inviscid and perfectly conducting fluid sphere in rapid rotation and permeated by a background magnetic field. Such a system admits normal modes in the form of torsional oscillations, namely azimuthal motions of cylinders coaxial with the rotation axis. We analyse this system for a particular background magnetic field that provides a new closed form normal mode solution. We derive Wentzel–Kramers–Brillouin–Jeffreys (WKBJ) approximations to the normal modes, and focus particularly on the reflections that take place on the rotation axis and at the equator. We propose a procedure to calculate the reflection coefficients and we discuss the analogy of our findings with well-known seismological results. Our analytical results are tested against numerical calculations and show good agreement.

Beschreibung: The 69 Ma Prospector Mountain stock is located in southwestern Yukon in the northern Canadian Cordillera. This massive monzonite-syenogranite stock is thought to be the intrusive volcanic centre for the surrounding coeval Carmacks Group volcanics. Anomalous palaeomagnetic data from these volcanics have provided the only evidence for the commonly posited hypothesis that the Yukon-Tanana terrane (YTT) was part of the far-travelled (1950 ± 600 km northward) ‘Baja BC’ terrane from 70 to 50 Ma. All other geologic evidence and averaged palaeomagnetic data support a northward displacement of ~415 ± 15 km. This study provides a direct test of the Carmacks volcanics’ estimate and examines the possible causes of its anomalous results. Both the stock and volcanics are unmetamorphosed and rest unconformably on metamorphosed basement rocks of the YTT. Palaeomagnetic and mineral magnetic results from 17 of 19 tested sites (218 specimens) in the Prospector stock and its peripheral skarn isolated a stable thermoremanent magnetization (TRM) in magnetite or low-Ti titanomagnetite that was mostly determined on demagnetization between temperatures of 500 and 580 °C. The TRM has a direction of Decl. = 8.3°, Incl. = 82.4° ( N = 17, k = 71.9, α 95 = 4.2°), providing a non-significant northwards translation estimate of 70 ± 880 km for the YTT. The normal-polarity TRM direction at Prospector Mountain provides a highly significant palaeomagnetic reversals test with the reversed-polarity TRM of the 70 Ma Swede Dome stock, another volcanic centre of the Carmacks Group about 190 km to the north. The test affirms to a high probability that both stocks carry primary TRMs and have not been tectonically tilted significantly since emplacement. Combining the palaeopoles for these two stocks with that for the 75 Ma Mount Lorne volcanic centre stock about 210 km south of Prospector Mountain yields a combined northward translation estimate of 330 ± 400 km for the YTT since ~71 Ma. This estimate agrees closely with the 415 ± 15 km estimate from geological constraints on strike-slip motion along the inboard Tintina fault zone since ~70 Ma, and with the palaeomagnetic estimate of 650 ± 450 km northward since ~108 Ma from four mid-Cretaceous batholiths. In contrast, the palaeomagnetic directions from four areas of the ~70 Ma Carmacks volcanics are poorly clustered. Nonetheless their estimates of ~1900 km northward have been used for a quarter of a century to support a far-travelled ‘Baja BC’ tectonic model for the YTT. This paper discusses and concludes that the anomalous far-travelled estimates from the Carmacks volcanics are attributable to the unfortunate additive effects of inadequate averaging of secular variation, dipole offset error, unrecognized primary dip and other possible causes.

Beschreibung: Previous magnetotelluric (MT) studies of the high-temperature Coso geothermal system in California identified a subvertical feature of low resistivity (2–5 Ohm m) and appreciable lateral extent (〉1 km) in the producing zone of the East Flank field. However, these models could not reproduce gross 3-D effects in the recorded data. We perform 3-D full-tensor inversion and retrieve a resistivity model that out-performs previous 2-D and 3-D off-diagonal models in terms of its fit to the complete 3-D MT data set as well as the degree of modelling bias. Inclusion of secondary Z xx and Z yy data components leads to a robust east-dip (60) to the previously identified conductive East Flank reservoir feature, which correlates strongly with recently mapped surface faults, downhole well temperatures, 3-D seismic reflection data, and local microseismicity. We perform synthetic forward modelling to test the best-fit dip of this conductor using the response at a nearby MT station. We interpret the dipping conductor as a fractured and fluidized compartment, which is structurally controlled by an unmapped blind East Flank fault zone.

Beschreibung: The island of Terceira (38.5°N, 27.2°W) in the Azores Archipelago boasts a wealth of well-dated recent volcanic products; it is therefore the ideal location to improve the resolution of the full-vector geomagnetic field record for the Mid-Atlantic Ocean. We investigated 21 (sub-)sites sampled from 10 different cooling units; six of these cooling units are from the Holocene, the other four are approximately coeval with the age of known geomagnetic excursions. After successfully determining the palaeodirections from all but one site, we applied the multimethod palaeointensity approach to all sites. Here, we used three different palaeointensity methods: the IZZI-Thellier protocol, the Multispecimen technique, and the calibrated pseudo-Thellier method. For the Holocene sites four robust palaeointensity estimates were obtained, they mostly agree with data from the nearby island of São Miguel and data from the Iberian Peninsula and West Africa. Two older sites, dated at 25.7 and 46.3 ka, produced low palaeointensities: 24.4 and 24.1 μT, respectively. As both the IZZI-Thellier and pseudo-Thellier methods yielded technically acceptable results for these sites we used them to improve the pseudo-Thellier calibration relation and establish its applicability to lower palaeointensities. Sixty per cent of all sampled cooling units/ages in our study yielded a technically successful palaeointensity estimate; this unusually high success rate is attributed both to our sampling strategy and the implementation of the multimethod approach.

Beschreibung: Two lava flows from the Ceboruco volcano in west-central Mexico were sampled for palaeomagnetic dating. The younger one was emitted in 1870 and used to validate the method, while the older one known as Ceboruco flow is of unknown age but probably younger than ~1005 AD and older than 1528 AD. Each flow was sampled in at least four sites, in order to unravel between site variations. For the 1870 flow, between site differences were notable and additionally post-cooling block movements were important; therefore, two sites had to be rejected. Three sites from the vent area and one at the tip of the 1870 flow provided well-constrained directions. This is also true for Ceboruco lava flow, and overall mean directions and palaeointensities were then used for palaeomagnetic dating applying the Matlab tool archaeo_dating and the global palaeosecular variation model SHA.DIF.14k. For the 1870 lava flow, the dating resulted in an age ranging between 1755 and 1871 AD (95 per cent probability level), which includes the real emplacement age. In addition, the Ceboruco lava flow was dated between 1000 and 1134 AD, which is close to the large plinian Jala eruption producing the crater of Ceboruco volcano around 1005 AD. This age is older than previously assumed and suggests an emplacement only shortly after the Jala eruption. As this lava flow is considered to be the youngest one of seven post-Jala lava flows, the age also defines a period of inactivity of Ceboruco volcano of about 730–860 yr before the historic 1870 eruption. Future volcanic hazard analysis will have to take into account this result. Our work also shows that multiple sampling of single lava flows is important to obtain a reliable mean direction. Sampling sites have to be carefully selected so that they represent un-tilted parts of the flows. We interpret this to be the case for the Ceboruco lava flow, while three of the six sites of the 1870 lava flow may have been partly or completely affected by movements after thermoremanent magnetization acquisition. Unfortunately, no better sites were found for this flow.

Beschreibung: We develop a high-resolution regional gravity field model by a combination of spherical harmonics, band-limited spherical radial basis functions (SRBFs) and the residual terrain model (RTM) technique. As the main input data set, we employ a dense terrestrial gravity database (3–6 stations km –2 ), which enables gravity field modelling up to very short spatial scales. The approach is based on the remove–compute–restore methodology in which all the parts of the signal that can be modelled are removed prior to the least-squares adjustment in order to smooth the input gravity data. To this end, we utilize degree-2159 spherical harmonic models and the RTM technique using topographic models at 2 arcsec resolution. The residual short-scale gravity signal is modelled via the band-limited Shannon SRBF expanded up to degree 21 600, which corresponds to a spatial resolution of 30 arcsec. The combined model is validated against GNSS/levelling-based height anomalies, independent surface gravity data, deflections of the vertical and terrestrial vertical gravity gradients achieving an accuracy of 2.7 cm, 0.53 mGal, 0.39 arcsec and 279 E in terms of the RMS error, respectively. A key aspect of the combined approach, especially in mountainous areas, is the quality of the RTM. We therefore compare the performance of two RTM techniques within the innermost zone, the tesseroids and the polyhedron. It is shown that the polyhedron-based approach should be preferred in rugged terrain if a high-quality RTM is required. In addition, we deal with the RTM computations at points located below the reference surface of the residual terrain which is known to be a rather delicate issue.

Beschreibung: Archaeomagnetic observations are key to recovering the behaviour of the geomagnetic field over the past few millennia. The corresponding data set presents a highly heterogeneous distribution in both space and time. Furthermore, the data are affected by substantial age and experimental uncertainties. In order to mitigate these detrimental properties, time-dependent global archaeomagnetic field models are usually constructed under spatial and temporal regularization constraints, with the use of bootstrap techniques to account for data uncertainties. The models so obtained are the product of an adjustable trade-off between goodness-of-fit and model complexity. The spatial complexity is penalized by means of a norm reflecting the minimization of Ohmic dissipation within the core. We propose in this study to resort to alternative spatial constraints relying on the statistics of a numerical dynamo simulation with Earth-like features. To that end, we introduce a dynamo norm in an ensemble least-squares iterative framework, the goal of which is to produce single-epoch models of the archaeomagnetic field. We first validate this approach using synthetic data. We next construct a redistributed archaeomagnetic data set between 1200 BC and 2000 AD by binning the data in windows of 40-yr width. Since the dynamo norm is not adjustable, we can legitimately calculate a resolution matrix to quantify the resolving power of the available archaeomagnetic data set. Gauss coefficients are resolved up to spherical harmonic degree 3 for the first thousand years of the interval, to degree 4 for the next thousand years and to degree 5 during the last millennium. These conclusions are based on the distribution and uncertainties that characterize the data set, and do not take into account the possible presence of outliers. Comparison between our model, called AmR, and previously published archaeomagnetic field models confirms the archaeomagnetic resolution analysis: it highlights the dichotomy between data-driven coefficients for which model predictions coincide (within their respective uncertainties), and prior-driven coefficients. This study opens the way to physics-based models of the archaeomagnetic field; future work will be devoted to integrating the framework here introduced into a time-dependent ensemble assimilation scheme.

Beschreibung: We carried out an archaeomagnetic directional study on 38 oriented samples (bricks and baked clays) collected from four archaeological locations at three provinces in China. The ages of our samples, spanning from ~3000 BCE to ~1300 CE, were constrained using a combination of archaeological context, radiocarbon dating and stratigraphic information. Rock magnetic results demonstrate that the main magnetic minerals of the studied samples are magnetite and/or hematite in single domain and superparamagnetic states. A total of 20 new reliable archaeodirectional data from 12 independent sites are obtained after thermal demagnetization experiments. These are the first set of archaeodirectional data in China produced since the 1990s. The published data are largely from the past 2 kyr and data from older time periods are rare. Our new data, especially those from period older than 3 ka, fill many gaps of the presently published dataset and will provide strong constraints on palaeosecular variation of the geomagnetic field in Eastern Asia and on the improvement of global models. Quite a few inflection points in the direction of the geomagnetic field are recorded in Eastern Asia over the past 10 kyr and some of them synchronize with the maximums or minimums of the palaeointensity. The palaeosecular variation rates are very low (based on present data distribution) before 2000 BCE and then start to increase and fluctuate afterward, which is generally consistent with the pattern of palaeointensity variations in this area.

Beschreibung: In this paper we present the potential of a new compact superconducting gravimeter (GWR iGrav) designed for groundwater monitoring. At first, 3 yr of continuous gravity data are evaluated and the performance of the instrument is investigated. With repeated absolute gravity measurements using a Micro-g Lacoste FG5, the calibration factor (–894.8 nm s –2 V –1 ) and the long-term drift of this instrument (45 nm s –2 yr –1 ) are estimated for the first time with a high precision and found to be respectively constant and linear for this particular iGrav. The low noise level performance is found similar to those of previous superconducting gravimeters and leads to gravity residuals coherent with local hydrology. The iGrav is located in a fully instrumented hydrogeophysical observatory on the Durzon karstic basin (Larzac plateau, south of France). Rain gauges and a flux tower (evapo-transpiration measurements) are used to evaluate the groundwater mass balance at the local scale. Water mass balance demonstrates that the karst is only capacitive: all the rainwater is temporarily stored in the matrix and fast transfers to the spring through fractures are insignificant in this area. Moreover, the upper part of the karst around the observatory appears to be representative of slow transfer of the whole catchment. Indeed, slow transfer estimated on the site fully supports the low-flow discharge at the only spring which represents all groundwater outflows from the catchment. In the last part of the paper, reservoir models are used to characterize the water transfer and storage processes. Particular highlights are done on the advantages of continuous gravity data (compared to repeated campaigns) and on the importance of local accurate meteorological data to limit misinterpretation of the gravity observations. The results are complementary with previous studies at the basin scale and show a clear potential for continuous gravity time-series assimilation in hydrological simulations, even on heterogeneous karstic systems.

Beschreibung: A thorough understanding of time-dependent noise in Global Navigation Satellite System (GNSS) position time-series is necessary for computing uncertainties in any signals found in the data. However, estimation of time-correlated noise is a challenging task and is complicated by the difficulty in separating noise from signal, the features of greatest interest in the time-series. In this paper, we investigate how linear trends affect the estimation of noise in daily GNSS position time-series. We use synthetic time-series to study the relationship between linear trends and estimates of time-correlated noise for the six most commonly cited noise models. We find that the effects of added linear trends, or conversely de-trending, vary depending on the noise model. The commonly adopted model of random walk (RW), flicker noise (FN) and white noise (WN) is the most severely affected by de-trending, with estimates of low-amplitude RW most severely biased. FN plus WN is least affected by adding or removing trends. Non-integer power-law noise estimates are also less affected by de-trending, but are very sensitive to the addition of trend when the spectral index is less than one. We derive an analytical relationship between linear trends and the estimated RW variance for the special case of pure RW noise. Overall, we find that to ascertain the correct noise model for GNSS position time-series and to estimate the correct noise parameters, it is important to have independent constraints on the actual trends in the data.

Beschreibung: Seismic waves sensitive to the outermost part of the Earth's liquid core seem to be affected by a stably stratified layer at the core–mantle boundary. Such a layer could have an observable signature in both long-term and short-term variations of the magnetic field of the Earth, which are used to probe the flow at the top of the core. Indeed, with the recent SWARM mission, it seems reasonable to be able to identify waves propagating in the core with period of several months, which may play an important role in the large-scale dynamics. In this paper, we characterize the influence of a stratified layer at the top of the core on deep quasi-geostrophic (Rossby) waves. We compute numerically the quasi-geostrophic eigenmodes of a rapidly rotating spherical shell, with a stably stratified layer near the outer boundary. Two simple models of stratification are taken into account, which are scaled with commonly adopted values of the Brunt–Väisälä frequency in the Earth's core. In the absence of magnetic field, we find that both azimuthal wavelength and frequency of the eigenmodes control their penetration into the stratified layer: the higher the phase speed, the higher the permeability of the stratified layer to the wave motion. We also show that the theory developed by Takehiro & Lister for thermal convection extends to the whole family of Rossby waves in the core. Adding a magnetic field, the penetrative behaviour of the quasi-geostrophic modes (the so-called fast branch) is insensitive to the imposed magnetic field and only weakly sensitive to the precise shape of the stratification. Based on these results, the large-scale and high-frequency modes (1–2 month periods) may be detectable in the geomagnetic data measured at the Earth's surface, especially in the equatorial area where the modes can be trapped.

Beschreibung: Previous formalisms for determining the static perturbation of spherically symmetric self-gravitating elastic Earth models due to displacement dislocations deal with each infinitesimal element of the fault system in its epicentral reference frame. In this work, we overcome this restriction and present novel and compact formulas for obtaining the perturbation due to the whole fault system in an arbitrary and common reference frame. Furthermore, we show that, even in an arbitrary reference frame, it is still possible to discriminate the contributions associated with the polar, bipolar and quadrupolar patterns of the seismic source response, as well as their relation with the along strike, along dip and tensile components of the displacement dislocation. These results allow a better understanding of the relation between the static perturbation and the whole fault system, and find direct applications in geodetic problems, like the modelling of long-wavelength geoid or gravity data from GRACE and GOCE space missions and of the perturbation of the deviatoric inertia tensor of the Earth.

Beschreibung: Sensitive instruments like strainmeters and tiltmeters are necessary for measuring slowly varying low amplitude Earth deformations. Nonetheless, laser and fibre interferometers are particularly suitable for interrogating such instruments due to their extreme precision and accuracy. In this paper, a practical design of a simple pendulum borehole tiltmeter based on laser fibre interferometric displacement sensors is presented. A prototype instrument has been constructed using welded borosilicate with a pendulum length of 0.85 m resulting in a main resonance frequency of 0.6 Hz. By implementing three coplanar extrinsic fibre Fabry-Perot interferometric probes and appropriate signal filtering, our instrument provides tilt measurements that are insensitive to parasitic deformations caused by temperature and pressure variations. This prototype has been installed in an underground facility (Rustrel, France) where results show accurate measurements of Earth strains derived from Earth and ocean tides, local hydrologic effects, as well as local and remote earthquakes. The large dynamic range and the high sensitivity of this tiltmeter render it an invaluable tool for numerous geophysical applications such as transient fault motion, volcanic strain and reservoir monitoring.

Beschreibung: A 3-D magnetotelluric (MT) inversion code using unstructured tetrahedral elements has been developed in order to correct the topographic effect by directly incorporating it into computational grids. The electromagnetic field and response functions get distorted at the observation sites of MT surveys because of the undulating surface topography, and without correcting this distortion, the subsurface structure can be misinterpreted. Of the two methods proposed to correct the topographic effect, the method incorporating topography explicitly in the inversion is applicable to a wider range of surveys. For forward problems, it has been shown that the finite element method using unstructured tetrahedral elements is useful for the incorporation of topography. Therefore, this paper shows the applicability of unstructured tetrahedral elements in MT inversion using the newly developed code. The inversion code is capable of using the impedance tensor, the vertical magnetic transfer function (VMTF), and the phase tensor as observational data, and it estimates the subsurface resistivity values and the distortion tensor of each observation site. The forward part of the code was verified using two test models, one incorporating topographic effect and one without, and the verifications showed that the results were almost the same as those of previous works. The developed inversion code was then applied to synthetic data from a MT survey, and was verified as being able to recover the resistivity structure as well as other inversion codes. Finally, to confirm its applicability to the data affected by topography, inversion was performed using the synthetic data of the model that included two overlapping mountains. In each of the cases using the impedance tensor, the VMTF and the phase tensor, by including the topography in the mesh, the subsurface resistivity was determined more proficiently than in the case using the flat-surface mesh. Although the locations of the anomalies were not accurately estimated by the inversion using distorted impedance tensors due to the slightly undervalued gain, these locations were correctly estimated by using undistorted impedance tensors or adding VMTFs in the data. Therefore, it can be concluded that the inversion using the unstructured tetrahedral element effectively prevents the misinterpretation of subsurface resistivity and recovers subsurface resistivity proficiently by representing the topography in the computational mesh.

Beschreibung: Separating the contribution of different hematite coercivity grains to the magnetic fabric is a standing problem in rock magnetism because of the common occurrence of thermochemical alterations when measuring the anisotropy of thermal remanence. A technique that eliminates this bias is presented, which is useful when there is a need to separate the fabric of detrital from pigmentary hematite, for example. The method is based on stepwise thermal demagnetization of saturation isothermal remanent magnetizations (IRMs) applied orthogonally on three sister specimens, allowing calculation of the anisotropy tensor from the three components of each demagnetized IRM vector, avoiding the necessity of having to apply IRMs to thermochemically altered specimens. Vector subtraction allows determining the anisotropy tensor for specific unblocking-temperature ranges. The anisotropies of the pigmentary, specular and total hematite of the Mauch Chunk Formation red beds of Pennsylvania have been measured from an oriented block sample and results are compared to previous anisotropy measurements performed using the high-field anisotropy of isothermal remanence technique (hf-AIR), which measures total undifferentiated hematite. Experiments were conducted using non-saturating 1 T and fully saturating 5.5 T fields: both experimental sets seem capable of measuring the orientation of the specularite anisotropy principal axes, but 5.5 T are needed to capture the orientation of the higher coercivity pigmentary grains. The magnitudes of the principal axes, instead, are only faithfully measured using 5.5 T fields and yield somewhat higher anisotropies than those measured by hf-AIR. The fundamental requirement for this technique is homogeneous material among the three sister specimens, which is a significant limitation; homogeneity tests allow assessment of applicability of the method and reliability of the results.

Beschreibung: This study investigates pressure effects on the magnetic properties of non-interacting single-domain (SD) magnetite. Using a high-pressure cell specially designed for a Magnetic Property Measurement System, magnetic hysteresis measurements were conducted under high pressures of up to 1 GPa on natural plagioclase crystals containing much acicular SD magnetite. Coercivity and saturation magnetization were nearly constant with pressure, while saturation remanent magnetization and coercivity of remanence decreased with pressure at moderate rates of –8 per cent GPa –1 and –18 per cent GPa –1 , respectively. These results suggest that temperature effects govern the magnetic behaviour of acicular SD magnetite grains in the middle and lower crusts.

Beschreibung: Data from the Gravity Recovery and Climate Experiment (GRACE) satellite mission can be used to estimate the mass change rate for separate drainage systems (DSs) of the Greenland Ice Sheet (GrIS). One approach to do so is by inversion of the level-2 spherical harmonic data to surface mass changes in predefined regions, or mascons. However, the inversion can be numerically unstable for some individual DSs. This occurs mainly for DSs with a small mass change signal that are located in the interior region of Greenland. In this study, we present a modified mascon inversion approach with an improved implementation of the constraint equations to obtain better estimates for individual DSs. We use separate constraints for mass change variability in the coastal zone, where run-off takes place, and for the ice sheet interior above 2000 m, where mass changes are smaller. A multi-objective optimization approach is used to find optimal prior variances for these two areas based on a simulation model. Correlations between adjacent DSs are suppressed when our optimized prior variances are used, while the mass balance estimates for the combination of the DSs that make up the GrIS above 2000 m are not affected significantly. The resulting mass balance estimates for some DSs in the interior are significantly improved compared to an inversion with a single constraint, as determined by a comparison with mass balance estimates from surface mass balance modelling and discharge measurements. The rate of mass change of the GrIS for the period of January 2003 to December 2012 is found to be –266.1 ± 17.2 Gt yr –1 in the coastal zone and areas below 2000 m, and +8.2 ± 8.6 Gt yr –1 in the interior region.