Search

Physical geodesy applies potential theory to study the Earth’s gravitational field in space outside and up to a few km inside the Earth’s mass. Among various tools offered by this theory, boundary-value problems are particularly popular for the transformation or continuation of gravitational field parameters across space. Traditional problems, formulated and solved as early as in the nineteenth century, have been gradually supplemented with new problems, as new observational methods and data are available. In most cases, the emphasis is on formulating a functional relationship involving two functions in 3-D space; the values of one function are searched but unobservable; the values of the other function are observable but with errors. Such mathematical models (observation equations) are referred to as deterministic. Since observed data burdened with observational errors are used for their solutions, the relevant stochastic models must be formulated to provide uncertainties of the estimated parameters against which their quality can be evaluated. This article discusses the boundary-value problems of potential theory formulated for gravitational data currently or in the foreseeable future used by physical geodesy. Their solutions in the form of integral formulas and integral equations are reviewed, practical estimators applicable to numerical solutions of the deterministic models are formulated, and their related stochastic models are introduced. Deterministic and stochastic models represent a complete solution to problems in physical geodesy providing estimates of unknown parameters and their error variances (mean squared errors). On the other hand, analyses of error covariances can reveal problems related to the observed data and/or the design of the mathematical models. Numerical experiments demonstrate the applicability of stochastic models in practice. [ABSTRACT FROM AUTHOR]

Relict vegetation of water-splashed petrifying rocks with calcareous tufa formation dominated by the fern Adiantum capillus-veneris and wetland bryophytes (phytosociological class Adiantetea) was investigated in Georgia, Caucasus biodiversity hotspot. The study brings the first phytosociological data on this scarce community in the Caucasus based on a novel dataset of vegetation plot records. A classification analysis revealed two main vegetation communities. They are represented by the Caucasian community recorded in the Mtkvari River valley inside the Tbilisi city (E Georgia) and the Colchic community developed under the humid, warm-temperate climate of Colchis (W Georgia), one of the key refugia of Tertiary flora across W Eurasia. To compare newly recorded vegetation types and so far described associations of the Adiantetea class, we used a broader dataset of the analogous communities reported from Europe and surroundings. It associated the Caucasian community with the pan-Mediterranean association Eucladio-Adiantetum. It was characterized by subhalophytes (e.g. Samolus valerandi), Eucladium verticillatum as a dominant bryophyte, and the absence of relict or endemic vascular plant species. The Colchic community represented a previously undescribed community. Therefore, we designated a new association Saxifrago cymbalariae-Adiantetum capilli-veneris delineated by the characteristic species of the Colchic-Caucasian territory (e.g. Hedera colchica, Saxifraga cymbalaria) and Palustriella commutata as the dominant bryophyte. Both associations belong to the Mediterranean-Atlantic alliance Adiantion. Described communities require conservation attention for their rarity, refugial character and presence of relict and endemic species. Modifications of the hydrological regime and construction activities are among their most important potential threats. [ABSTRACT FROM AUTHOR]

Physical heights is one of the most important topics in physical geodesy. Their original concept, introduced in the 19-th century, defined physical heights as lengths of plumblines of the Earth's gravity field between the geoid and points of interest. There are orthometric heights of surface points, that have been traditionally estimated by spirit levelling and measured gravity; however, the knowledge of the density distribution of topographic masses (masses between the geoid and Earth's surface) is required that significantly affects their determinability. This was also the main reason why a new type of physical heights was proposed in the mid of the 20-th century. Normal heights approximate orthometric heights in a sense that the Earth's gravity field is replaced by the normal gravity field, an analytic model based on the theory of an equipotential ellipsoid. This height system has been introduced since that time in different countries in Europe and beyond. Contrary to the classical height system based on orthometric heights, its counterpart based on normal heights may have slightly different definitions. Moreover, normal heights are often defined as heights of points above the quasigeoid. This contribution reviews alternative definitions of normal heights and respective height systems. It is argued that both orthometric and normal heights refer to the geoid. In the case physical heights are estimated by satellite positioning, normal heights must be computed through the height anomaly estimated at each point of interest, whether it is below, at or above the Earth's surface. On the contrary, orthometric heights of all points along the same plumbline, be it below, at or above the Earth's surface, are estimated by introducing one value of the geoid height. Normal heights of surface points can be estimated by spirit levelling easier than orthometric heights as no topographic mass density hypothesis is required; however, one has to keep in mind the gravity field approximation used both for their definition and realization. [ABSTRACT FROM AUTHOR]

The Caucasus harbours unique forest vegetation so far only little studied using the Braun-Blanquet approach. This study is mostly based on a dataset (N = 110) of original phytosociological relevés of oak-hornbeam and ravine forests in the Eastern Greater Caucasus, Georgia. Their unsupervised classification produced seven communities. Five belong to oak-hornbeam forests (order Lathyro-Carpinetalia caucasicae). Of the zonal Caucasian alliance Crataego-Carpinion, the association Corno australis-Carpinetum inhabits valleys of the Greater Caucasus, and Clinopodio umbrosi-Carpinetum is confined to the warm Eastern Greater Caucasus promontories. The association Astrantio maximae-Carpinetum of the alliance Astrantio-Carpinion represents distinctive Caucasian montaneoak-hornbeam forests. The other two communities, documented by a few relevés, were described at the community level only. Within ravine forests (order Aceretalia pseudoplatani), we introduce a new Caucasian alliance Pachyphragmo macrophyllae-Tilion begoniifoliae with two associations. Valeriano tiliifoliae-Ulmetum glabrae comprises Caucasian montane ravine forests whereas Hedero pastuchovii-Aceretum velutini inhabits the foothills of the Eastern Greater Caucasian. To provide a broader context of the recognized communities, an expanded dataset (N = 231) of original relevés and previously published relevés of Georgian deciduous forests was analysed. It indicated a major turnover in species composition following biogeographical patterns presumably driven by macroclimate and vegetation history. [ABSTRACT FROM AUTHOR]

In geophysical and geodetic studies, gravity inversion is typically performed such that observed gravity values are first continued downward onto a regular (planar, spherical or spheroidal) surface by solving an inverse integral transform, which originates from a classical solution to the first boundary-value problem in potential theory. A typical example is continuing gravity observed at the topographic surface down to the mean sea level (geoid). Nowadays, gravity-dedicated satellite missions and aerial gravimetry provide gravity data above the topographic surface in addition to classical terrestrial gravity observations. For specific purposes (such as gravity data combination and validation, or quasigeoid determination), satellite and aerial gravity observations have to be continued to the irregular topographic surface. In this study, we address this issue by formulating a functional model for a spectral downward continuation of selected gravitational field quantities to an irregular topographic surface. Moreover, we generalize this functional model to allow for transformation between different types of gravitational field quantities. In particular, we derive spectral weights for estimation of the disturbing potential or disturbing/anomalous gravity at the Earth's surface by combining the first-, second- and third-order radial gradients of the disturbing potential (disturbing gradients). The correctness of the developed combined spectral estimator is verified in a closed-loop test based on synthetic satellite disturbing gradients. The combined spectral estimator is applied to simulated satellite disturbing gradients polluted by a realistic Gaussian noise. Results of the numerical experiments show that the combined spectral estimator puts the highest importance on the least polluted disturbing gradient, while the contribution of the least accurate disturbing gradient is negligible. An important advantage of this spectral combination method is that no matrix inversion with numerical instabilities requiring regularization is needed. [ABSTRACT FROM AUTHOR]

The static Earth’s gravitational field has traditionally been described in geodesy and geophysics by the gravitational potential (geopotential for short), a scalar function of 3-D position. Although not directly observable, geopotential functionals such as its first- and second-order gradients are routinely measured by ground, airborne and/or satellite sensors. In geodesy, these observables are often used for recovery of the static geopotential at some simple reference surface approximating the actual Earth’s surface. A generalized mathematical model is represented by a surface integral equation which originates in solving Dirichlet’s boundary-value problem of the potential theory defined for the harmonic geopotential, spheroidal boundary and globally distributed gradient data. The mathematical model can be used for combining various geopotential gradients without necessity of their re-sampling or prior continuation in space. The model extends the apparatus of integral equations which results from solving boundary-value problems of the potential theory to all geopotential gradients observed by current ground, airborne and satellite sensors. Differences between spherical and spheroidal formulations of integral kernel functions of Green’s kind are investigated. Estimated differences reach relative values at the level of 3% which demonstrates the significance of spheroidal approximation for flattened bodies such as the Earth. The observation model can be used for combined inversion of currently available geopotential gradients while exploring their spectral and stochastic characteristics. The model would be even more relevant to gravitational field modelling of other bodies in space with more pronounced spheroidal geometry than that of the Earth. [ABSTRACT FROM AUTHOR]

The downward continuation (DWC) of the gravity anomalies from the Earth's surface to the geoid is still probably the most problematic step in the precise geoid determination. It is this step that motivates the quasi-geoid users to opt for Molodenskij's rather than Stokes's theory. The reason for this is that the DWC is perceived as suffering from two major flaws: first, a physically meaningful DWC technique requires the knowledge of the irregular topographical density; second, the Poisson DWC, which is the only physically meaningful technique we know, presents itself mathematically in the form of Fredholm integral equation of the 1st kind. As Fredholm integral equations are often numerically ill-conditioned, this makes some people believe that the DWC problem is physically ill-posed. According to a revered French mathematician Hadamard, the DWC problem is physically well-posed and as such gives always a finite and unique solution. The necessity of knowing the topographical density is, of course, a real problem but one that is being solved with an ever increasing accuracy; so sooner or later it will allow us to determine the geoid with the centimetre accuracy. [ABSTRACT FROM AUTHOR]

New spherical integral formulas between components of the second- and third-order gravitational tensors are formulated in this article. First, we review the nomenclature and basic properties of the second- and third-order gravitational tensors. Initial points of mathematical derivations, i.e., the second- and third-order differential operators defined in the spherical local North-oriented reference frame and the analytical solutions of the gradiometric boundary-value problem, are also summarized. Secondly, we apply the third-order differential operators to the analytical solutions of the gradiometric boundary-value problem which gives 30 new integral formulas transforming (1) vertical-vertical, (2) vertical-horizontal and (3) horizontal-horizontal second-order gravitational tensor components onto their third-order counterparts. Using spherical polar coordinates related sub-integral kernels can efficiently be decomposed into azimuthal and isotropic parts. Both spectral and closed forms of the isotropic kernels are provided and their limits are investigated. Thirdly, numerical experiments are performed to test the consistency of the new integral transforms and to investigate properties of the sub-integral kernels. The new mathematical apparatus is valid for any harmonic potential field and may be exploited, e.g., when gravitational/magnetic second- and third-order tensor components become available in the future. The new integral formulas also extend the well-known Meissl diagram and enrich the theoretical apparatus of geodesy. [ABSTRACT FROM AUTHOR]

Values of scalar, vector and second-order tensor parameters of the Earth's gravitational field have been collected by various sensors in geodesy and geophysics. Such observables have been widely exploited in different parametrization methods for the gravitational field modelling. Moreover, theoretical aspects of these quantities have extensively been studied and well understood. On the other hand, new sensors for observing gravitational curvatures, i.e., components of the third-order gravitational tensor, are currently under development. As the gravitational curvatures represent new types of observables, their exploitation for modelling of the Earth's gravitational field is a subject of this study. Firstly, the gravitational curvature tensor is decomposed into six parts which are expanded in terms of third-order tensor spherical harmonics. Secondly, gravitational curvature boundary-value problems defined for four combinations of the gravitational curvatures are formulated and solved in spectral and spatial domains. Thirdly, properties of the corresponding sub-integral kernels are investigated. The presented mathematical formulations reveal some important properties of the gravitational curvatures and extend the so-called Meissl scheme, i.e., an important theoretical framework that relates various parameters of the Earth's gravitational field. [ABSTRACT FROM AUTHOR]

In this study we assume that a gravitational curvature tensor, i.e. a tensor of third-order directional derivatives of the Earth's gravitational potential, is observable at satellite altitudes. Such a tensor is composed of ten different components, i.e. gravitational curvatures, which may be combined into vertical-vertical-vertical, vertical-vertical-horizontal, vertical-horizontal-horizontal and horizontal-horizontal-horizontal gravitational curvatures. Firstly, we study spectral properties of the gravitational curvatures. Secondly, we derive new quadrature formulas for the spherical harmonic analysis of the four gravitational curvatures and provide their corresponding analytical error models. Thirdly, requirements for an instrument that would eventually observe gravitational curvatures by differential accelerometry are investigated. The results reveal that measuring third-order directional derivatives of the gravitational potential imposes very high requirements on the accuracy of deployed accelerometers which are beyond the limits of currently available sensors. For example, for orbital parameters and performance similar to those of the GOCE mission, observing third-order directional derivatives requires accelerometers with the noise level of $${\sim}10^{-17}\,\hbox {m}\,\hbox {s}^{-2}$$ Hz $$^{-1/2}$$ . [ABSTRACT FROM AUTHOR]

A new mathematical model for evaluation of the third-order (disturbing) gravitational tensor is formulated in this article. Firstly, we construct corresponding differential operators for the components of the third-order (disturbing) gravitational tensor in a spherical local north-oriented frame. We show that the differential operators may efficiently be decomposed into an azimuthal and an isotropic part. The differential operators are even more simplified for a certain class of isotropic kernels. Secondly, the differential operators are applied to the well-known integrals of Newton, Abel-Poisson, Pizzetti and Hotine. In this way, 40 new integral formulas are derived. The new integral formulas allow for evaluation of the components of the third-order (disturbing) gravitational tensor from density distribution, disturbing gravitational potential, gravity anomalies and gravity disturbances. Thirdly, we investigate the behaviour of the corresponding integral kernels in the spatial domain. The new mathematical formulas extend the theoretical apparatus of geodesy, i.e. the well-known Meissl scheme, and reveal important properties of the third-order gravitational tensor. They may be exploited in geophysical studies, continuation of gravitational field quantities and analysing the gradiometric-geodynamic boundary value problem. [ABSTRACT FROM AUTHOR]

With the advent of geodetic satellite missions mapping almost globally the Earth's gravitational field, new methods and theoretical approaches have been developed and investigated to fully exploit the potential of their new observables. Besides estimating values of numerical coefficients in harmonic series of the gravitational potential, new applications emerged such as data validation and combination. In this contribution, new integral transformations are presented which transform principal components of the terrestrial deflection of the vertical onto disturbing satellite-to-satellite tracking and gradiometric data at altitude. Using spherical approximation, necessary integral kernel functions are derived in both spectral and closed forms. The behaviour of isotropic kernel functions is studied and the new integral transformations are tested in a closed-loop simulation using synthetic terrestrial and satellite data synthesized from a global gravitational model. New integral transformations can be used for data validation and combination purposes. [ABSTRACT FROM AUTHOR]

Integral transformations of gravitational gradients onto a Gravity Recovery And Climate Experiment (GRACE) type of observable are derived in this article. The gravitational gradients represent components of the gravitational tensor in the local north-oriented frame. The GRACE type of observable corresponds to a difference between two gravitational vectors as projected onto the line of sight between the two GRACE satellites. In total, three integral transformations relating vertical-vertical, vertical-horizontal and horizontal-horizontal gravitational gradients with the GRACE type of observable are provided. Spectral and closed forms of corresponding isotropic kernels are derived for each transformation. Special cases show that the integral transformations are general and relate gravitational gradients to many other quantities of the gravitational field, such as the gravitational vector, and its radial and tangential components. Correctness of the mathematical derivations is validated in a closed-loop simulation using synthetic data. [ABSTRACT FROM AUTHOR]

Due to the ESA's satellite mission GOCE launched in March 2009, gravitational gradients sampled along the orbital trajectory approximately 250 km above the Earth's surface have become available. Since 2010, gravitational gradients have routinely been applied in geodesy for the derivation of global Earth's gravitational models provided in terms of fully normalized coefficients in a spherical harmonic series representation of the Earth's gravitational potential. However, in geophysics, gravitational gradients observed by spaceborne instruments have still been applied relatively seldom. This contribution describes their possible geophysical applications in structural studies where gravitational gradients observed at satellite altitudes are compared with those derived by a spectral forward modeling technique using available global models of selected Earth's mass components as input data. In particular, GOCE gravitational gradients are interpreted in terms of a superposition principle of gravitation as combined gravitational effects generated by a homogeneous reference ellipsoid of revolution, mean topographic and ice mass density distributions, depth-dependent mass density contrasts within bathymetry and lateral mass density anomalies with sediments and crustal layers. Respective gravitational effects are one by one removed from gravitational gradients observed at approximately 250 km elevation above ground. Removing respective gravitational gradients from observed gravitational gradients gradually reveals problematic geographic areas with model deficiencies. For the full interpretation of observed gravitational gradients, deficiencies of CRUST2.0 must be corrected and effects of deeper laying mass anomalies not included in the study considered. These findings are confirmed by parameters describing spectral properties of the gravitational gradients. The methodology can be applied for validating Earth's gravitational models and for constraining crustal models in the development phase. [ABSTRACT FROM AUTHOR]

We derive expressions for computing the gravitational field (potential and its radial derivative) generated by an arbitrary homogeneous or laterally varying density contrast layer with a variable depth and thickness based on methods for a spherical harmonic analysis and synthesis of gravity field. The newly derived expressions are utilised in the gravimetric forward modelling of major known density structures within the Earth's crust (excluding the ocean density contrast) beneath the geoid surface. The gravitational field quantities due to the sediments and crust components density contrasts, shown in numerical examples, are computed using the 2 × 2 arc-deg discrete data from the global crustal model CRUST2.0. These density contrasts are defined relative to the adopted value of the reference crustal density of 2670 kgm. All computations are realised globally on a 1 × 1 arc-deg geographical grid at the Earth's surface. The maxima of the gravitational signal due to the sediments density contrast are mainly along continental shelf regions with the largest sedimentary deposits. The corresponding maxima due to the consolidated crust components density contrast are over areas of the largest continental crustal thickness with variable geological structure. [ABSTRACT FROM AUTHOR]

We developed and applied a novel numerical scheme for a gravimetric forward modelling of the Earth's crustal density structures based entirely on methods for a spherical analysis and synthesis of the gravitational field. This numerical scheme utilises expressions for the gravitational potentials and their radial derivatives generated by the homogeneous or laterally varying mass density layers with a variable height/depth and thickness given in terms of spherical harmonics. We used these expressions to compute globally the complete crust-corrected Earth's gravity field and its contribution generated by the Earth's crust. The gravimetric forward modelling of large known mass density structures within the Earth's crust is realised by using global models of the Earth's gravity field (EGM2008), topography/bathymetry (DTM2006.0), continental ice-thickness (ICE-5G), and crustal density structures (CRUST2.0). The crust-corrected gravity field is obtained after modelling and subtracting the gravitational contribution of the Earth's crust from the EGM2008 gravity data. These refined gravity data mainly comprise information on the Moho interface and mantle lithosphere. Numerical results also reveal that the gravitational contribution of the Earth's crust varies globally from 1,843 to 12,010 mGal. This gravitational signal is strongly correlated with the crustal thickness with its maxima in mountainous regions (Himalayas, Tibetan Plateau and Andes) with the presence of large isostatic compensation. The corresponding minima over the open oceans are due to the thin and heavier oceanic crust. [ABSTRACT FROM AUTHOR]

In geophysical studies investigating the lithosphere structure, the gravitational field generated by the ocean density contrast (i.e., bathymetry-generated gravitational field) represents a significant amount of the signal to be modelled and subsequently removed from the Earth's gravity field. The ocean density contrast is typically calculated as the difference between the mean density values of the Earth's crust and seawater. The approximation of the actual seawater density distribution by its mean value yields relative errors up to about 2% in computed quantities of the gravitational field. To reduce these errors, a more realistic model of the seawater density distribution is utilized based on the analysis of existing oceanographic data of salinity, temperature, and pressure (depth). We study the accuracy of the bathymetry-generated gravitational field quantities formulated for a depth-dependent model of the seawater density distribution. This density distribution approximates the seawater density variations due to an increasing pressure with depth, whereas smaller lateral density variations caused by salinity, temperature, and other oceanographic factors are not taken into consideration. The error analysis reveals that the approximation of the seawater density by the depth-dependent density model reduces the maximum errors to less than 0.6%. The corresponding depth-averaged errors are below 0.1%. The depth-dependent seawater density model is further facilitated in expressions for computing the bathymetry-generated gravitational field quantities by means of the spherical bathymetric (ocean bottom depth) functions. The numerical realization reveals large differences in the results obtained with and without consideration of the depth-dependent seawater density distribution. The maxima of absolute differences reach 201 m/s and 16.5 mGal in computed values of the potential and attraction, respectively. The application of the depth-dependent seawater density model thus significantly improves the accuracy in the forward modelling of the bathymetric gravitational field quantities. [ABSTRACT FROM AUTHOR]

Spherical harmonic transformation is of practical interest in geodesy for transformation of globally distributed quantities such as gravity between space and frequency domains. The increasing spatial resolution of the latest and forthcoming gravitational models pose true computational challenges for classical algorithms since serious numerical instabilities arise during the computation of the respective base functions of the spherical harmonic expansion. A possible solution is the evaluation of the associated Legendre functions in the Fourier domain where numerical instabilities can be circumvented by an independent frequency-wise scaling of numerical coefficients into a numerically suitable double precision range. It is then rather straightforward to commit global fast data transformation into the Fourier domain and to evaluate subsequently spherical harmonic coefficients. For the inverse, the computation of respective Fourier coefficients from a given spherical harmonic model is performed as an inverse Fast Fourier Transform into globally distributed data points. The two-step formulation turns out to be stable even for very high resolutions as well as efficient when using state-of-the-art shared memory/multi-core architectures. In principle, any functional of the geopotential can be computed in this way. To give an example for the overall performance of the algorithm, we transformed an equiangular 1 arcmin grid of terrain elevation data corresponding to spherical harmonic degree and order 10800. [ABSTRACT FROM AUTHOR]

spectral representation of the topographic corrections to gravity field quantities is formulated in terms of spherical height functions. When computing the far-zone contributions to the topographic corrections, various types of the truncation coefficients are applied to a spectral representation of Newton's integral. In this study we utilise Molodensky's truncation coefficients in deriving the expressions for the far-zone contributions to topographic corrections. The expressions for computing the far-zone gravity field contributions corrected for the effect of topography are then obtained by combining the expressions for the far-zone contributions to the gravity field quantities and to the respective topographic corrections, both expressed in terms of Molodensky's truncation coefficients. The numerical examples of the far-zone contributions to the topographic corrections and to the topography-corrected gravity field quantities are given over the study area situated in the Canadian Rocky Mountains with adjacent planes. Coefficients of the global elevation and geopotential models are used as the input data. [ABSTRACT FROM AUTHOR]

We improve the precision and computation speed of the fully-normalized associated Legendre functions (fnALFs) for ultra-high degrees and orders of spherical harmonic transforms. We take advantage of their numerical behaviour of and propose two new methods for solving an underflow/overflow problem in their calculation. We specifically discuss the application of the two methods in the fixed-order increasing-degree recursion computation technique. The first method uses successive ratios of fnALFs and the second method, called the Midway method, starts iteration from tiny initial values, which are still in the range of the IEEE double-precision environment, rather than from sectorial fnALFs. The underflow/overflow problem in the successive ratio method is handled by using a logarithm-based method and the extended range arithmetic. We validate both methods using numerical tests and compare their results with the X-number method in terms of precision, stability, and speed. The results show that the relative precision of the proposed methods is better than 10−9 for the maximum degree of 100000, compared to results derived by the high precision Wolfram's Mathematica software. Average CPU times required for evaluation of fnALFs over different latitudes demonstrate that the two proposed methods are faster by about 10–30% and 20–90% with respect to the X-number method for the maximum degree in the range of 50–65000. [ABSTRACT FROM AUTHOR]

An Correction to this paper has been published: https://doi.org/10.1007/s12224-021-09391-x [ABSTRACT FROM AUTHOR]

The right-hand side of Eq. (35) in the original article equals one. Thus, Eq. (35) is incorrect, and we suggest the following correct form instead: [ABSTRACT FROM AUTHOR]

River gravel-bar habitats are highly endangered. They are still well-preserved in the Caucasus, but developing conservation strategies is burdened by the lack of data from this region. We studied vegetation and habitat types on gravel bars of 22 rivers in Georgia, including successional stages from open early-successional herbaceous vegetation to scrub. We distinguished five vegetation types based on vegetation physiognomy and β-flexible clustering of species composition, and described them as phytosociological vegetation units: Early-successional herbaceous vegetation at higher elevations was described as the new association Epilobietum colchici and that at lower elevations as the Petrorhagia saxifraga–Crepis foetida community. The grassland dominated by Calamagrostis pseudophragmites and scrub vegetation were assigned to the associations previously described from Central Europe (Tussilagini farfarae-Calamagrostietum pseudophragmitae, Salici purpureae-Myricarietum germanicae and Salici incanae-Hippophaëtum rhamnoidis). We established diagnostic plant species for each type using the fidelity calculation and related these types to environmental variables. We further compared them with the previously published data on gravel-bar vegetation from the Russian part of the Caucasus and with European systems of habitat classification. This study demonstrates that vegetation and habitat types occurring in Georgia correspond to those recognized earlier in Europe, and can be easily linked to the European systems of habitat classification. Unlike in other parts of Europe, these habitats are still well-preserved on rivers with natural hydrological dynamics, but they are threatened by plans of dam building and other river regulations. Our study provides baseline data for developing conservation strategies for the Caucasian gravel-bar habitats. [ABSTRACT FROM AUTHOR]

The regional gravimetric geoid solved using boundary-value problems of the potential theory is usually determined in two computational steps: (1) downward continuing ground gravity data onto the geoid using inverse Poisson's integral equation in a mass-free space and (2) evaluating geoidal heights by applying Stokes integral to downward continued gravity. In this contribution, the two integration steps are combined in one step and the so-called one-step integration method in spherical approximation is implemented to compute the regional gravimetric geoid model. Advantages of using the one-step integration method instead of the two integration steps include less computational cost, more stable numerical computation and better utilization of input ground gravity data (reduced in each integration step to avoid edge effects). A discrete form of the one-step integral equation is used to convert mean values of ground gravity anomalies into mean values of geoidal heights. To evaluate mean values of the integral kernel in the vicinity of the computation point, a fast and numerically accurate analytical formula is proposed using planar approximation. The proposed formula is tested to determine the regional gravimetric geoid of the Auvergne test area, France. Results show a good agreement of the estimated geoid with geoidal heights estimated at GNSS-levelling reference points, with the standard deviation for the difference of 3.3 cm. Considering the uncertainty of geoidal heights derived at the GNSS/levelling reference points, one can conclude the geoid models computed by the one-step and two-step integration methods have negligible differences. Thus, the one-step method can be recommended for regional geoid modelling with its methodological and numerical advantages. [ABSTRACT FROM AUTHOR]

The uncertainties of the geoidal heights estimated from ground gravity data caused by their spatial distribution and noise are investigated in this study. To test these effects, the geoidal heights are estimated from synthetic ground gravity data using the Stokes- Helmert approach. Five different magnitudes of the random noise in ground gravity data and three types of their spatial distribution are considered in the study, namely grid, semigrid and random. The noise propagation is estimated for the two major computational steps of the Stokes-Helmert approach, i.e., the downward continuation of ground gravity and Stokes's integration. Numerical results show that in order to achieve the 1-cm geoid, the ground gravity data should be distributed on the grid or semi-grid with the average angular distance less than 2′. If they are randomly distributed (scattered gravity points), the 1-cm geoid cannot be estimated if the average angular distance between scattered gravity points is larger than 1′. Besides, the noise of the gravity data for the tree types of their spatial distribution should be below 1 mGal to estimate the 1-cm geoid. The advantage of interpolating scattered gravity points onto the regular grid, rather than using them directly, is also investigated in this study. Numerical test shows that it is always worth interpolating the scattered points to the regular grid except if the scattered gravity points are sparser than 5′. [ABSTRACT FROM AUTHOR]

The disjunct occurrence of Pedicularis exaltata in the White Carpathians (Czech Republic), isolated by more than 500 km from the nearest populations in the Eastern Carpathians, has been considered one of the mysteries of the Western Carpathian flora. We used molecular methods (AFLP, ITS and cpDNA sequencing) to reconstruct a possible scenario of the evolution of P. exaltata and its closely related congener P. hacquetii and to evaluate their differentiation. We paid particular attention to the origin of the isolated population in the White Carpathians. We also analysed the vegetation composition at the sampling sites to characterize the habitat preferences of the Pedicularis species and compare different sites. The pattern of molecular variability does not support a species boundary between P. exaltata and P. hacquetii. These assumed species should be merged into one bearing the name Pedicularis hacquetii following the priority rule, as also proposed by several early morphology-based taxonomic studies. The cpDNA variability pattern supports an evolutionary scenario involving a distribution centre (refugium) in the Eastern Carpathians from which the species expanded to the more westerly parts of its extant geographical range. Low population differentiation in AFLPs, high gene diversity and high DW index in isolated populations indicate that this expansion was contiguous rather than based on long-distance dispersal. Under such a scenario, the White Carpathian population is considered one of isolated relict populations scattered throughout the Carpathians and the East European Plain. A complete phylogeographic reconstruction of P. hacquetii s.l., however, requires ancient DNA analysis of herbarium specimens of nowadays extinct populations of the East European Plain. Extant habitats of P. hacquetii s.l. consist of species-rich dry-mesic to mesic (sub)montane grasslands and subalpine tall-herb growths, which may be relics of vegetation widespread in the late Glacial/early Holocene, when the contiguous expansion of the species possibly proceeded. The analysis of species with high fidelity to P. hacquetii s.l. throughout its range indicated that middle altitudes of south-western White Carpathians belong to the regions most suitable for the taxon within the Czech Republic and Slovakia. We found no support for the late introduction/long-distance dispersal scenarios and consider the relic survival scenario most parsimonious for the White Carpathians. [ABSTRACT FROM AUTHOR]

The launch of gravity-dedicated satellite missions at the beginning of the new millennium led to an accuracy improvement of global Earth gravity field models (GGMs). One of these missions was the Gravity field and steady-state Ocean Circulation Explorer (GOCE) launched in 2009. As the first European Space Agency's Earth Explorer Mission, the satellite carried a novel instrument, a 3-D gradiometer, which allowed measurement of the second-order directional derivatives of the gravitational potential (gravitational gradients) with a uniform quality and a near-global coverage. The main mission goal was to determine the static Earth's gravity field with the ambitious precision of 1-2 cm in terms of geoid heights and 1 mGal in terms of gravity anomalies for spatial resolution of 100 km (half wavelength at the equator). More than three years of the outstanding measurements resulted in three levels of data products (Level 0, Level 1b and Level 2), six releases of GGMs, and several global grids of gravitational gradients. The grids, which represent a step between gravitational gradients measured directly along the GOCE orbit and those represented by GGMs, found their usage mainly in geophysical applications. In this contribution, we validate the official Level 2 product GRD_SPW_2 using height anomalies over two test areas located in central and northern Europe (Czechia/Slovakia and Norway). A mathematical model based on the least-squares spectral weighting is employed with corresponding spectral weights estimated for validation of gravitational gradient grids. This model continues gravitational gradients from the mean orbital altitude of GOCE down to the irregular Earth's surface (not to a sphere) and transforms them to height anomalies in one computational step. Analytical downward continuation errors of the model are estimated using a closed-loop test. Prior to the comparison of height anomalies estimated from gravitational gradients with their reference values derived from Global Navigation Satellite Systems (GNSS)/levelling over the two test areas, the gravitational gradients and reference data are corrected for all systematic effects such as the tide system conversion. Moreover, the high-frequency part of the gravitational signal is estimated and subtracted from reference data as it is attenuated in the gravitational gradients measured by GOCE. A relative improvement between the release 6 and release 2 gradient grids reaches 48 % in terms of height anomalies in Czechia/Slovakia. The relative improvement in Norway is even more significant and reaches 55 % . The release 6 of the official Level 2 product GRD_SPW_2 gained the absolute accuracy with the standard deviation of 8.7 cm over Czechia/Slovakia and 9.3 cm over Norway. Article Highlights: A mathematical model based on the least-squares spectral weighting for validation of the official Level 2 product GRD_SPW_2 is presented Possibilities of the spectral combination method for the downward continuation of gravitational gradients and their transformation into height anomalies are investigated A review of quality of the official Level 2 product GRD_SPW_2 with respect to GNSS/levelling points over two test areas (Czechia/Slovakia and Norway) is provided [ABSTRACT FROM AUTHOR]

Please be informed that the wrong table was inadvertently published in the following paper [ABSTRACT FROM AUTHOR]

The geoid-to-quasigeoid separation is often computed only approximately as a function of the simple planar Bouguer gravity anomaly and the height of the computation point while disregarding the contributions of terrain geometry and anomalous topographic density as well as the sub-geoid masses. In this study we demonstrate that these contributions are significant and, therefore, should be taken into consideration when investigating the relation between the normal and orthometric heights particularly in the mountainous, polar and geologically complex regions. These contributions are evaluated by applying the spectral expressions for gravimetric forward modelling and using the EIGEN-6C4 gravity model, the Earth2014 datasets of terrain, ice thickness and inland bathymetry and the CRUST1.0 sediment and (consolidated) crustal density data. Since the global crustal density models currently available (e.g. CRUST1.0) have a limited accuracy and resolution, the comparison of individual density contributions is-for consistency-realized with a limited spectral resolution up to a spherical harmonic degree 360 (or 180). The results reveal that the topographic contribution globally varies between $$-$$ 0.33 and 0.57 m, with maxima in Himalaya and Tibet. The contribution of ice considerably modifies the geoid-to-quasigeoid separation over large parts of Antarctica and Greenland, where it reaches $$\sim $$ 0.2 m. The contributions of sediments and bedrock are less pronounced, with the values typically varying only within a few centimetres. These results, however, have still possibly large uncertainties due to the lack of information on the actual sediment and bedrock density. The contribution of lakes is mostly negligible; its maxima over the Laurentian Great Lakes and the Baikal Lake reach only several millimetres. The contribution of the sub-geoid masses is significant. It is everywhere negative and reaches extreme values of $$-$$ 4.43 m. According to our estimates, the geoid-to-quasigeoid separation globally varies within $$-$$ 4.19 and 0.26 m while the corresponding values computed according to a classical definition are only negative and reach extreme values of $$-$$ 3.5 m. A comparison of these results reveals that inaccuracies caused by disregarding the terrain geometry and mass density heterogeneities distributed within the topography and below the geoid surface can reach $$\pm $$ 2 m or more in the mountainous regions. [ABSTRACT FROM AUTHOR]

Integral transforms of the disturbing gravitational potential derived from satellite altimetry onto satellite gradiometric data are formulated, investigated and applied in this article. First, corresponding differential operators, that relate the disturbing gravitational potential to the six components of the disturbing gradiometric tensor in the spherical local north-oriented frame, are applied to the spherical Abel-Poisson integral equation. This yields six new integral equations for which respective kernel functions are given in both spectral and spatial forms. Second, truncation error formulas for each of the integral transforms are provided in the spectral form. Also expressions for the corresponding truncation error coefficients are derived. Third, practical estimators for evaluation of the disturbing gravitational gradients are formulated, and their correctness and expected accuracy are investigated. Finally, the practical estimators are applied for validation of a sample of the gradiometric data provided by the GOCE satellite mission. Obtained results demonstrate applicability of the new apparatus as an alternative validation method of the satellite gravitational gradients. [ABSTRACT FROM AUTHOR]

In the last few decades, satellites have acquired various potential data sets hundreds of kilometers above the Earth's surface. Conventionally, these global magnetic and gravitational data sets are approximated by using spherical harmonics that allow straightforward work with both fields outside the Earth's mass. In this article, we present an alternative approach for working with potential data in mass-free space given over a regular coordinate grid on a spherical surface. The algorithm is based on an iterative scheme and the Poisson integral equation for the sphere. With help from the Fourier transform, global potential (magnetic or gravitational) data can efficiently be continued from a mean orbital sphere down to a reference surface without using the spherical harmonics. This is illustrated both with simulated magnetic field data and with real data from the satellite gradiometry mission GOCE. In the case of simulated magnetic data and the downward continuation for 450 km, we have achieved a root mean square at the level of 0.05 nT, while it was <1 E (eotvos) for real GOCE data continued for 250 km. The crucial point is to apply the algorithm twice as a large part of noise can be removed from the input data. [ABSTRACT FROM AUTHOR]

We compute globally the consolidated crust-stripped gravity disturbances/anomalies. These refined gravity field quantities are obtained from the EGM2008 gravity data after applying the topographic and crust density contrasts stripping corrections computed using the global topography/bathymetry model DTM2006.0, the global continental ice-thickness data ICE-5G, and the global crustal model CRUST2.0. All crust components density contrasts are defined relative to the reference crustal density of 2,670 kg/m. We demonstrate that the consolidated crust-stripped gravity data have the strongest correlation with the crustal thickness. Therefore, they are the most suitable gravity data type for the recovery of the Moho density interface by means of the gravimetric modelling or inversion. The consolidated crust-stripped gravity data and the CRUST2.0 crust-thickness data are used to estimate the global average value of the crust-mantle density contrast. This is done by minimising the correlation between these refined gravity and crust-thickness data by adding the crust-mantle density contrast to the original reference crustal density of 2,670 kg/m. The estimated values of 485 kg/m (for the refined gravity disturbances) and 481 kg/m (for the refined gravity anomalies) very closely agree with the value of the crust-mantle density contrast of 480 kg/m, which is adopted in the definition of the Preliminary Reference Earth Model (PREM). This agreement is more likely due to the fact that our results of the gravimetric forward modelling are significantly constrained by the CRUST2.0 model density structure and crust-thickness data derived purely based on methods of seismic refraction. [ABSTRACT FROM AUTHOR]

We analyse spatial and spectral characteristics of various refined gravity data used for modelling and gravimetric interpretation of the crust-mantle interface and the mantle-lithosphere structure. Depending on the purpose of the study, refined gravity data have either a strong or weak correlation with the Moho depths (Moho geometry). The compilation of the refined gravity data is purely based on available information on the crustal density structure obtained from seismic surveys without adopting any isostatic hypothesis. We demonstrate that the crust-stripped relative-to-mantle gravity data have a weak correlation with the CRUST2.0 Moho depths of about 0.02. Since gravitational signals due to the crustal density structure and the Moho geometry are subtracted from gravity field, these refined gravity data comprise mainly the information on the mantle lithosphere and sub-lithospheric mantle. On the other hand, the consolidated crust-stripped gravity data, obtained from the gravity field after applying the crust density contrast stripping corrections, comprise mainly the gravitational signal of the Moho geometry, although they also contain the gravitational signal due to anomalous mass density structures within the mantle. In the absence of global models of the mantle structure, the best possible option of computing refined gravity data, suitable for the recovery/refinement of the Moho interface, is to subtract the complete crust-corrected gravity data from the consolidated crust-stripped gravity data. These refined gravity data, that is, the homogenous crust gravity data, have a strong absolute correlation of about 0.99 with the CRUST2.0 Moho depths due to removing a gravitational signal of inhomogeneous density structures within the crust and mantle. Results of the spectral signal decomposition and the subsequent correlation analysis reveal that the correlation of the homogenous crust gravity data with the Moho depths is larger than 0.9 over the investigated harmonic spectrum up to harmonic degree 90. The crust-stripped relative-to-mantle gravity data correlate substantially with the Moho depths above harmonic degree 50 where the correlation exceeds 0.5. [ABSTRACT FROM AUTHOR]

This paper represents a milestone in the UNB effort to formulate an accurate and self-consistent theory for regional geoid determination. To get the geoid to a sub-centimetre accuracy, we had to formulate the theory in a spherical rather than linear approximation, advance the modelling of the effect of topographic mass density, formulate the solid spherical Bouguer anomaly, develop the probabilistic downward continuation approach, incorporate improved satellite determined global gravitational models and introduce a whole host of smaller improvements. Having adopted Auvergne, an area in France as our testing ground, where the mean standard deviation of observed gravity values is 0.5 mGal, according to the Institute Geographique Nationale (Duquenne in Proceedings of the 1st international symposium of the international gravity field service "gravity field of the earth", International gravity field service meeting, Istanbul, Turkey, 2006), we obtained the standard deviation of the gravity anomalies continued downward to the geoid, as estimated by minimizing the L2 norm of their residuals, to be in average 3 times larger than those on the surface with large spikes underneath the highest topographic points. The standard deviations of resulting geoidal heights range from a few millimetres to just over 6 cm for the highest topographic points in the Alpine region (just short of 2000 m). The mean standard deviations of the geoidal heights for the whole region are only 0.6 cm, which should be considered quite reasonable even if one acknowledges that the area of Auvergne is mostly flat. As one should expect, the main contributing factors to these uncertainties are the Poisson probabilistic downward continuation process, with the maximum standard deviation just short of 6 cm (the average value of 2.5 mm) and the topographic density uncertainties, with the maximum value of 5.6 cm (the average value of 3.0 mm). The comparison of our geoidal heights with the testing geoidal heights, obtained for a set of 75 control points (regularly spaced throughout the region), shows the mean shift of 13 cm which is believed to reflect the displacement of the French vertical datum from the geoid due to sea surface topography. The mean root square error of the misfit is 3.3 cm. This misfit, when we consider the estimated accuracy of our geoid, indicates that the mean standard deviation of the "test geoid" is about 3 cm, which makes it about 5 times less accurate than the Stokes–Helmert computed geoid. [ABSTRACT FROM AUTHOR]

In terrestrial studies, Bouguer gravity data is routinely computed by adopting various numerical schemes, starting from the most basic concept of approximating the actual topography by an infinite Bouguer plate, through adding a planar terrain correction to account for a local/regional terrain geometry, to more advanced schemes that involve the computation of the topographic gravity correction by taking into consideration a gravitational contribution of the whole topography while adopting a spherical (or ellipsoidal) approximation. Moreover, the topographic density information has significantly improved the gravity forward modeling and interpretations, especially in polar regions (by accounting for a density contrast of polar glaciers) and in regions characterized by a complex geological structure. Whereas in geodetic studies (such as a gravimetric geoid modeling) only the gravitational contribution of topographic masses above the geoid is computed and subsequently removed from observed (free-air) gravity data, geophysical studies focusing on interpreting the Earth's inner structure usually require the application of additional stripping gravity corrections that account for known anomalous density structures in order to reveal an unknown (and sought) density structure or density interface. In planetary studies, numerical schemes applied to compile Bouguer gravity maps might differ from terrestrial studies due to two reasons. While in terrestrial studies the topography is defined by physical heights above the geoid surface (i.e., the geoid-referenced topography), in planetary studies the topography is commonly described by geometric heights above the geometric reference surface (i.e., the geometric-referenced topography). Moreover, large parts of a planetary surface have negative heights. This obviously has implications on the computation of the topographic gravity correction and consequently Bouguer gravity data because in this case the application of this correction not only removes the gravitational contribution of a topographic mass surplus, but also compensates for a topographic mass deficit. In this study, we examine numerically possible options of computing the topographic gravity correction and consequently the Bouguer gravity data in planetary applications. In agreement with a theoretical definition of the Bouguer gravity correction, the Bouguer gravity maps compiled based on adopting the geoid-referenced topography are the most relevant. In this case, the application of the topographic gravity correction removes only the gravitational contribution of the topography. Alternative options based on using geometric heights, on the other hand, subtract an additional gravitational signal, spatially closely correlated with the geoidal undulations, that is often attributed to deep mantle density heterogeneities, mantle plumes or other phenomena that are not directly related to a topographic density distribution. [ABSTRACT FROM AUTHOR]

A composite of powders of semi-Heusler ferromagnetic shape memory and pure titanium was successfully prepared by spark plasma sintering at the temperature of 950 °C. Sintering resulted in the formation of small precipitates and intermetallic phases at the heterogeneous interfaces. Various complementary experimental methods were used to fully characterize the microstructure. Imaging methods including transmission and scanning electron microscopy with energy dispersive X-ray spectroscopy revealed a position and chemical composition of individual intermetallic phases and precipitates. The crystalline structure of the phases was examined by a joint refinement of X-ray and neutron diffraction patterns. It was found that Co38Ni33Al29 decomposes into the B2-(Co,Ni)Al matrix and A1-(Co,Ni,Al) particles during sintering, while Al, Co and Ni diffuse into Ti forming an eutectic two phase structure with C9-Ti2(Co,Ni) precipitates. Complicated interface intermetallic structure containing C9-Ti2(Co,Ni), B2-(Co,Ni)Ti and L21-(Co,Ni)(Al,Ti) was completely revealed. In addition, C9-Ti2(Co,Ni) and A1-(Co,Ni,Al) precipitates were investigated by an advanced method of small angle neutron scattering. This study proves that powder metallurgy followed by spark plasma sintering is an appropriate technique to prepare bulk composites from very dissimilar materials. [ABSTRACT FROM AUTHOR]

In geodesy, the geoid and the quasigeoid are used as a reference surface for heights. Despite some similarities between these two concepts, the differences between the geoid and the quasigeoid (i.e. the geoid-to-quasigeoid correction) have to be taken into consideration in some specific applications which require a high accuracy. Over the world's oceans and marginal seas, the quasigeoid and the geoid are identical. Over the continents, however, the geoid-to-quasigeoid correction could reach up to several metres especially in the mountainous, polar and geologically complex regions. Various methods have been developed and applied to compute this correction regionally in the spatial domain using detailed gravity, terrain and crustal density data. These methods utilize the gravimetric forward modelling of the topographic density structure and the direct/inverse solutions to the boundary-value problems in physical geodesy. In this article, we provide a brief summary of existing theoretical and numerical studies on the geoid-to-quasigeoid correction. We then compare these methods with the newly developed procedure and discuss some numerical and practical aspects of computing this correction. In global applications, the geoid-to-quasigeoid correction can conveniently be computed in the spectral domain. For this purpose, we derive and present also the spectral expressions for computing this correction based on applying methods for a spherical harmonic analysis and synthesis of global gravity, terrain and crustal structure models. We argue that the newly developed procedure for the regional gravity-to-potential conversion, applied for computing the geoid-to-quasigeoid correction in the spatial domain, is numerically more stable than the existing inverse models which utilize the gravity downward continuation. Moreover, compared to existing spectral expressions, our definition in the spectral domain takes not only the terrain geometry but also the mass density heterogeneities within the whole Earth into consideration. In this way, the geoid-to-quasigeoid correction and the respective geoid model could be determined more accurately. [ABSTRACT FROM AUTHOR]

Downward continuation of the potential data (gravity or magnetic) helps to interpret contributing sources by transforming the physical signal to their neighborhood. This paper applies three continuation approaches (Landweber's iteration, the direct method and the inverse-matrix approach with Tikhonov's regularization), based on the Poisson integral equation, and four integration schemes to the second vertical derivative of the anomalous potential $$T_{\rm zz}$$ obtained from GOCE data. In the experiments, $$T_{\rm zz}$$ was downward continued for 250 km in the area of Central Europe with special attention to edge effects. For the integration schemes, which use prior information outside the region of interest to reduce edge effects, the best agreement with TIM-r4 was achieved by the inverse-matrix approach with Tikhonov's regularization (RMS = 1.15 eotvos). On the contrary, without prior information the iterative approach performed with RMS = 1.17 eotvos. [ABSTRACT FROM AUTHOR]

The global crustal model CRUST1.0 (refined using additional global datasets of the solid topography, polar ice sheets and geoid) is used in this study to estimate the average densities of major crustal structures. We further use this refined model to compile the gravity field quantities generated by the Earth's crustal structures and to investigate their spatial and spectral characteristics and their correlation with the crustal geometry in context of the gravimetric Moho determination. The analysis shows that the average crustal density is 2,830 kg/m, while it decreases to 2,490 kg/m when including the seawater. The average density of the oceanic crust (without the seawater) is 2,860 kg/m, and the average continental crustal density (including the continental shelves) is 2,790 kg/m. The correlation analysis reveals that the gravity field corrected for major known anomalous crustal density structures has a maximum (absolute) correlation with the Moho geometry. The Moho signature in these gravity data is seen mainly at the long-to-medium wavelengths. At higher frequencies, the Moho signature is weakening due to a noise in gravity data, which is mainly attributed to crustal model uncertainties. The Moho determination thus requires a combination of gravity and seismic data. In global studies, gravimetric methods can help improving seismic results, because (1) large parts of the world are not yet sufficiently covered by seismic surveys and (2) global gravity models have a relatively high accuracy and resolution. In regional and local studies, the gravimetric Moho determination requires either a detailed crustal density model or seismic data (for a combined gravity and seismic data inversion). We also demonstrate that the Earth's long-wavelength gravity spectrum comprises not only the gravitational signal of deep mantle heterogeneities (including the core-mantle boundary zone), but also shallow crustal structures. Consequently, the application of spectral filtering in the gravimetric Moho determination will remove not only the gravitational signal of (unknown) mantle heterogeneities, but also the Moho signature at the long-wavelength gravity spectrum. [ABSTRACT FROM AUTHOR]

New integral formulas for upward/downward continuation of gravitational gradients onto gravitational gradients are derived in this article. They provide more options for continuation of gravitational gradient combinations and extend available mathematical apparatus formulated for this purpose up to now. The starting point represents the analytical solution of the spherical gradiometric boundary value problem in the spatial domain. Applying corresponding differential operators on the analytical solution of the spherical gradiometric boundary value problem, a total of 18 integral formulas are provided. Spatial and spectral forms of isotropic kernels are given and their behaviour for parameters of a GOCE-like satellite is investigated. Correctness of the new integral formulas and the isotropic kernels is tested in a closed-loop simulation. The derived integral formulas and the isotropic kernels form a theoretical basis for validation purposes and geophysical applications of satellite gradiometric data as provided currently by the GOCE mission. They also extend the well-known Meissl scheme. [ABSTRACT FROM AUTHOR]

The details of the mechanism of persistent luminescence were probed by investigating the trap level structure of SrMgSiO:Eu,R materials (R: Y, La-Lu, excluding Pm and Eu) with thermoluminescence (TL) measurements and Density Functional Theory (DFT) calculations. The TL results indicated that the shallowest traps for each SrMgSiO:Eu,R material above room temperature were always ca. 0.7 eV corresponding to a strong TL maximum at ca. 90 °C. This main trap energy was only slightly modified by the different co-dopants, which, in contrast, had a significant effect on the depths of the deeper traps. The combined results of the trap level energies obtained from the experimental data and DFT calculations suggest that the main trap responsible for the persistent luminescence of the SrMgSiO:Eu,R materials is created by charge compensation lattice defects, identified tentatively as oxygen vacancies, induced by the R co-dopants. [ABSTRACT FROM AUTHOR]