Regional Seafloor Topography by Extended Kalman Filtering of Marine Gravity Data without Ship-Track Information

An iterative Extended Kalman Filter (EKF) approach is proposed to recover a regional set of topographic heights composing an undersea volcanic mount by the successive combination of large numbers of gravity measurements at sea surface using altimetry satellite-derived grids and taking the error uncertainties into account. The integration of the non-linear Newtonian operators versus the radial and angular distances (and its first derivatives) enables the estimation process to accelerate and requires only few iterations, instead of summing Legendre polynomial series or using noise-degraded 2D-FFT decomposition. To show the effectiveness of the EKF approach, we apply it to the real case of the bathymetry around the Great Meteor seamount in the Atlantic Ocean by combining only geoid height/free-air anomaly datasets and using ship-track soundings as reference for validation. Topography of the Great Meteor seamounts structures are well-reconstructed, especially when regional compensation is considered. Best solution gives a RMS equal to 400 m with respect to the single beam depth observations and it is comparable to RMS obtained for ETOPO1 of about 365 m. Larger discrepancies are located in the seamount flanks due to missing high-resolution information for gradients. This approach can improve the knowledge of seafloor topography in regions where few echo-sounder measurements are available.

Lithosphere–asthenosphere interactions beneath northeast China and the origin of its intraplate volcanism

Northeast China hosts one of the largest Cenozoic intraplate volcanic regions in the world. However, the mechanisms that generate the volcanism, its spatial-temporal distribution, and compositional signatures remain highly debated due to the lack of high-resolution images of the mantle’s thermochemical structure. We jointly inverted new surface-wave dispersion data, surface heat flow, geoid height, and elevation data to image the fine-scale thermal and compositional structures beneath northeast China and infer regions of partial melting in the mantle. Our model reveals a complex circulation pattern in the asthenosphere and a highly variable lithospheric structure. Combining predictions from our model with independent geochemical data from recent basaltic volcanism, we demonstrate that the generation, location, and composition of intraplate volcanism in this region are controlled by the interaction between shallow asthenospheric circulation and lithospheric thickness. The modeling approach and correlations between basaltic composition and mantle state identified in our study are globally applicable to assessing mantle conditions over time in other continental regions.

Using Heat Flow Density Values Obtained in the Gulf of Cadiz and Gorringe Bank, Atlantic Ocean

The geothermal heat flow measured at the surface of the Earth is originated by different heat sources located at different depths of the planet. The main sources of heat flow in the crust are associated with radioactive decay of Uranium, Thorium and Potassium, in rocks. In some regions, additional heat sources must be considered such as exothermic chemical reactions. The value of the heat flow coming from deep regions, designated by “heat from the mantle”, must be obtained using indirect methods. In this work, the geoid height was used as indicator of alterations “in heat from the mantle” values, considering that the density decrease in regions with geoid height increase is related to high temperature values in the upper part of the mantle. The region on study is located in the Atlantic Ocean, SW of Cape St. Vincent and Cadiz Gulf. Temperature-depth values were obtained in twelve points of the region considering heat flow by conduction in the vertical direction, using published heat flow and thermal conductivity data. Layered models were made using data obtained in published seismic profiles. Moho depth values were used as lower boundary of the crust and mantle heat flow variations were made according geoid height increases. Ocean depth values between 2.5 and 4.3 km were used. A value of 5°C was used for temperature at the upper boundary (ocean bottom) of the models. Temperature calculus stops when a value of 1350 °C was attained. Lithosphere thickness is obtained considering this temperature value as temperature at the bottom of the lithosphere. Heat flow density values from 36 to 65.8 mW m−2 were used in the work with “heat from the mantle” values from 33 to 35 mw m−2. Curie Point Temperature (600°C) depths from 33 to 36 km were obtained. Lithosphere thickness values about 97 km were obtained in all the models.

The implementation and advancement of a regional geodetic vertical datum

The essential parameter in computing three-dimensional coordinate system is the height or depth of the Earth’s surface. It represents a particular reference surface that recognised as a vertical datum. The vertical datum is alienated into two foremost categories recognised as Mean Sea Level and Lowest Astronomical Tide. Different modifications approach, techniques and software programs are developed to determine vertical datum of a region with respect to geoid surface. This paper presents an effort to review and discuss the implementations and advancement of geodetic vertical datum based on geoid height reference surface. Hence, there are eight countries will be extracted and outlined in this paper consist of the United States of America, Australia, Taiwan, New Zealand, South Korea, Thailand, Philippines and Malaysia. An overview of geodetic vertical datum which implemented in these countries are summarised to support the future development of a regional vertical datum model. Then, the overview will also be utilised and analysed based on the essential elements and parameters for vertical datum model determination which include: data gathering, data input and analysis approach in order to develop a geodetic vertical datum model with good accuracy. These attempt and initiative are vital for the current and future implementation and advancement of geodetic vertical datum in the region of Malaysia across land and marine areas.

Potential and scientific requirements of optical clock networks for validating satellite gravity missions

<p>The GRACE mission, now continued as the GRACE-FO mission, has provided an unprecedented quantification of large-scale changes in the water cycle.<br>Meanwhile, stationary optical clocks show fractional instabilities below 10<sup>-18</sup> when averaged over an hour, and continue to be improved in terms of precision and accuracy, uptime, and transportability. The frequency of a clock is affected by the gravitational redshift, and thus depends on the local geopotential; a relative frequency change of 10<sup>-18</sup> corresponds to a geoid height change of about 1 cm. This effect could be exploited for sensing temporal geopotential changes via a network of clocks distributed at the Earth's surface. <br>Here, we concentrate on how the measurements of an ensemble of optical clocks connected accross Europe via optical fibre links could be used to validate and complement gravity field solutions from GRACE-FO and potential future gravity missions.<br>Through simulations it is shown how hydrology (water storage) and atmosphere (dry and wet air mass) variations over Europe could be observed with clock comparisons in a future network. We assume different scenarios for clock and GNSS uncertainties, where we deem the latter to be nessecary to separate local height changes from the mass redistribution signals. Our findings suggest that even under conservative assumptions -- a clock error of 10<sup>-18</sup> and vertical height control error of 1.4 mm for daily measurements -- hydrological signals at the annual time scale and atmospheric signals down to the weekly time scale could be observed.<br>However, the requirements to an optical clock network used for validation of GRACE-FO and future gravity missions are higher than that, which is demonstrated along with the according spatial resolutions.</p>

Harmonic Correction for Residual Terrain Modelling (RTM) Technique in Physical Geodesy Applications

<p>       In physical geodesy, the harmonic correction (HC), as one of the main problems when using residual terrain modelling (RTM), has become a research focus of high-frequency gravity field modelling. Over past decades, though various methods have been proposed to handle the HC issues for RTM technique, most of them focused on the HC for RTM gravity anomaly rather than other gravity functionals, such as RTM geoid height and gravity gradient. In practice, the HC for RTM geoid height was generally assumed to be negligible, but a quantification is yet studied. In this study, besides the highlighted HC for gravity anomaly in previous studies, the expressions of HC terms for RTM geoid height are provided in the framework of the classical condensation method under infinite Bouguer plate approximation. The errors involved by various assumption of the classical condensation method, e.g., mass inconsistency between infinite masses in the HC and limited masses in the RTM, and planar assumption of the Earth’s surface, are further studied. Based on the derived formulas, the quantification of HC for RTM geoid height when reference surface is expanded to degree and order of 2,159 is given. Our results showed the significance of HC for RTM geoid height, with values up to ~10 cm, in cm-level and mm-level geoid determination. With integration masses extending up to a sufficient distance, such as 1° from calculation point for the determination of RTM geoid height, the errors due to an infinite Bouguer plate approximation are neglectable small. The validation through comparison with terrestrial measurements proved that the HC terms provided in this study can improve the accuracy of RTM derived geoid height and are expected to be useful for applications of RTM technique in regional and global gravity field modelling.</p>

Continental Water Storage Changes Sensed by GRACE Satellite Gravimetry

Since its launch in March 2002, the Gravity Recovery And Climate Experiment (GRACE) mission has been mapping the time variations of the Earth’s gravity field with a precision of 2–3 cm in terms of geoid height at the surface resolution of 300–400 km. The unprecedented precision of this twin satellite system enables to detect tiny changes of gravity that are due to the water mass variations inside the fluid envelops of our planet. Once they are corrected from known gravitational contributions of the atmosphere and the oceans, the monthly and (bi)weekly GRACE solutions reveal the continental water storage redistributions, and mainly the dominant seasonal cycle in the largest drainage river basins such as Amazon, Congo, Mississippi. The potential differences measured between the twin GRACE satellites represent the sum of integrated surface waters (lakes and rivers), soil moisture, snow, ice and groundwater. Once they are inverted for estimating surface water mass densities, GRACE solutions are also used to establish the long-term mass balance of the ice sheets impacted by global warming, for quantifying the interannual variations of the major aquifers, as well as for surveying the hydrological signatures of intense meteorological events lasting a few days such as tropical hurricanes. This chapter describes GRACE gravity products and the different data processings used for mapping continental water storage variations, it also presents the most remarkable results concerning global continental hydrology and climate changes.