Estimation Of Topography Research Articles

The Influence of Mantle Structure on Dynamic Topography in Southern Africa

AbstractDue to relatively high terrain and negligible active tectonics, the southern Africa region boasts over 30 independent estimates of dynamic topography. These published estimates display a wide variance due to both the variety of methods used in computation and a lack of constraints on the regional mantle structure. Here we show that a focus on regional mantle structure is important to generate models of lithospheric and mantle dynamics. Global average mantle properties are not representative of a particular region, and it is necessary to generate viscosity profiles specific to a region of interest. We develop a Bayesian inversion using dynamic geoid kernels, existing seismic tomography models, and Slepian functions to invert for a localized radial viscosity profile that best explains the geoid in southern Africa. With an understanding of viscosity uncertainty, we constrain dynamic topography in southern Africa to lie between 1,000 and 2,000 m. Additionally, we model vertical displacements from 112 Global Navigation Satellite System stations across our region to examine the long‐term, long wavelength pattern of present‐day vertical motion, suggesting that a mean of 1.5 mm/yr (1: 0.8–2.0 mm/yr) of vertical motion may be related to ongoing dynamic topography. Our study demonstrates the utility of dynamic geoid kernels in local nonlinear inversions of non‐unique geophysical data. Furthermore, we present evidence that the mantle beneath southern Africa is generating significant dynamic support for and vertical displacement of the lithosphere in this region.

Biometry study of foveal isoplanatic patch variation for adaptive optics retinal imaging.

The change in ocular wavefront aberrations with visual angle determines the isoplanatic patch, defined as the largest field of view over which diffraction-limited retinal imaging can be achieved. Here, we study how the isoplanatic patch at the foveal center varies across 32 schematic eyes, each individualized with optical biometry estimates of corneal and crystalline lens surface topography, assuming a homogeneous refractive index for the crystalline lens. The foveal isoplanatic patches were calculated using real ray tracing through 2, 4, 6 and 8 mm pupil diameters for wavelengths of 400-1200 nm, simulating five adaptive optics (AO) strategies. Three of these strategies, used in flood illumination, point-scanning, and line-scanning ophthalmoscopes, apply the same wavefront correction across the entire field of view, resulting in almost identical isoplanatic patches. Two time-division multiplexing (TDM) strategies are proposed to increase the isoplanatic patch of AO scanning ophthalmoscopes through field-varying wavefront correction. Results revealed substantial variation in isoplanatic patch size across eyes (40-500%), indicating that the field of view in AO ophthalmoscopes should be adjusted for each eye. The median isoplanatic patch size decreases with increasing pupil diameter, coarsely following a power law. No statistically significant correlations were found between isoplanatic patch size and axial length. The foveal isoplanatic patch increases linearly with wavelength, primarily due to its wavelength-dependent definition (wavefront root-mean-squared, RMS <λ/14), rather than aberration chromatism. Additionally, ray tracing reveals that in strongly ametropic eyes, induced aberrations can result in wavefront RMS errors as large as λ/3 for an 8-mm pupil, with implications for wavefront sensing, open-loop ophthalmic AO, spectacle prescription and refractive surgery.

Coupled Fire–Atmosphere Simulations of the Split, Croatia, Wildfire

Abstract The Weather Research and Forecasting Spread FIRE (WRF-SFIRE) model was used to simulate the evolution of the fire perimeter for the Split wildfire in Croatia, evaluated in our previous detailed reconstruction of the event. A four-domain configuration was applied, with the innermost domain having a resolution of 500 m with a planetary boundary layer parameterization. The coupled fire grid was further downscaled to 33.3-m resolution, using the best available estimates of topography and fuel load data. WRF-SFIRE simulated slightly smaller fire areas in the four observed fire perimeter stages of the Split wildfire event, especially the northwestward spread in the last stage along the Adriatic coast. Additional experiments with multiple ignitions in the WRF-SFIRE model, representing the spotting mechanism, successfully simulated the northwestward spread in the event. The study also included sensitivity experiments that turned off the heat and moisture flux coupling between WRF and the fire grid. It was found that the simulated fire area is larger than the model run that included heat and moisture feedback from the fire grid to the atmosphere, which is contrary to some previous studies. The large wind speed and convergence along the fire line in our feedback-on simulation, which constrain the fire when there are feedback processes, were likely the cause of the smaller fire area. Last, areas of pine that burnt in the Split fire event favored crown fire spread, which is not represented in the current modeling framework of WRF-SFIRE.

Effects of using the consistent boundary flux method on dynamic topography estimates

SUMMARY Dynamic topography is defined as the deflection of Earth's surface due to the convecting mantle. ASPECT (Advanced Solver for Planetary Evolution, Convection, and Tectonics) is a continually evolving, finite element code that uses modern numerical methods to investigate problems in mantle convection. With ASPECT version 2.0.0 a consistent boundary flux (CBF) algorithm, used to calculate radial stresses at the model boundaries, was implemented into the released version of ASPECT. It has been shown that the CBF algorithm improves the accuracy of dynamic topography calculations by approximately one order of magnitude. We aim to evaluate the influence of the CBF algorithm and explore the geophysical implications of these improved estimates of dynamic topography changes along the East Coast of the United States. We constrain our initial temperature conditions using the tomography models SAVANI, S40RTS and TX2008, and combine them with a corresponding radial viscosity profile (2 for TX2008) and two different boundary conditions for a total of eight experiments. We perform simulations with and without the CBF method, which takes place during post-processing and does not affect the velocity solution. Our dynamic topography calculations are spatially consistent in both approaches, but generally indicate an increase in magnitude using the CBF method (on average ∼15 and ∼76 per cent absolute change in present-day instantaneous and rate of change of dynamic topography, respectively). This enhanced accuracy in dynamic topography calculations can be used to better evaluate the effects of mantle convection on surface processes including vertical land motions, sea level changes, and sedimentation and erosion. We explore results along the US East Coast, where a Pliocene shoreline has been deformed by dynamic topography change. An increased accuracy in estimates of dynamic topography can improve Pleistocene and Pliocene sea level reconstructions, which allow for a better understanding of past sea level changes and ice sheet stability.

Radial surface currents from space: An opportunity for mean dynamic topography estimation?

We analyze the feasibility of so-called radial surface currents as derived from Sentinel-1 Wave mode Doppler shift observations to support the geodetic estimation of the mean dynamic topography (MDT). Broadly defined as subtracting the geoid from the mean sea surface, this separation problem suffers from data inhomogeneities and usually requires strong prior assumptions about the MDT’s smoothness. Recent advancements in calibration and current retrieval methodology for Sentinel-1 data give reason to assess potential gains from this observation type, which can be linked to the gradient of the MDT under geostrophic approximation. Due to current lack of long time-series, we synthesize 10 years of observations from daily surface currents grids and include these in a parametric least-squares adjustment of the MDT. Our results utilizing the synthetic data in place of smoothness constraints are promising, showing 2 cm RMS agreement with a state-of-the-art MDT model.

InSAR-DEM Block Adjustment Model for Upcoming BIOMASS Mission: Considering Atmospheric Effects

The unique P-band synthetic aperture radar (SAR) instrument, BIOMASS, is scheduled for launch in 2024. This satellite will enhance the estimation of subcanopy topography, owing to its strong penetration and fully polarimetric observation capability. In order to conduct global-scale mapping of the subcanopy topography, it is crucial to calibrate systematic errors of different strips through interferometric SAR (InSAR) DEM (digital elevation model) block adjustment. Furthermore, the BIOMASS mission will operate in repeat-pass interferometric mode, facing the atmospheric delay errors introduced by changes in atmospheric conditions. However, the existing block adjustment methods aim to calibrate systematic errors in bistatic mode, which can avoid possible errors from atmospheric effects through interferometry. Therefore, there is still a lack of systematic error calibration methods under the interference of atmospheric effects. To address this issue, we propose a block adjustment model considering atmospheric effects. Our model begins by employing the sub-aperture decomposition technique to form forward-looking and backward-looking interferograms, then multi-resolution weighted correlation analysis based on sub-aperture interferograms (SA-MRWCA) is utilized to detect atmospheric delay errors. Subsequently, the block adjustment model considering atmospheric effects can be established based on the SA-MRWCA. Finally, we use robust Helmert variance component estimation (RHVCE) to build the posterior stochastic model to improve parameter estimation accuracy. Due to the lack of spaceborne P-band data, this paper utilized L-band Advanced Land Observing Satellite (ALOS)-1 PALSAR data, which is also long-wavelength, to emulate systematic error calibration of the BIOMASS mission. We chose climatically diverse inland regions of Asia and the coastal regions of South America to assess the model’s effectiveness. The results show that the proposed block adjustment model considering atmospheric effects improved accuracy by 72.2% in the inland test site, with root mean square error (RMSE) decreasing from 10.85 m to 3.02 m. Moreover, the accuracy in the coastal test site improved by 80.2%, with RMSE decreasing from 16.19 m to 3.22 m.

Automated Estimation of Sub-Canopy Topography Combined with Single-Baseline Single-Polarization TanDEM-X InSAR and ICESat-2 Data

TanDEM-X bistatic interferometric system successfully generated a high-precision, high-resolution global digital elevation model (DEM). However, in forested areas, two core problems make it difficult to obtain sub-canopy topography: (1) the penetrability of short-wave signals is limited, and the DEM obtained in dense forest areas contains a significant forest signal, that is, the scattering phase center (SPC) height; and (2) the single-baseline and single-polarization TanDEM-X interferometric synthetic aperture radar (InSAR) data cannot provide sufficient observations to make the existing physical model reversible for estimating the real surface phase, whereas the introduction of optical data makes it difficult to ensure data synchronization and availability of cloud-free data. To overcome these problems in accurately estimating sub-canopy topography from TanDEM-X InSAR data, this study proposes a practical method of sub-canopy topography estimation based on the following innovations: (1) An orthogonal polynomial model was established using TanDEM-X interferometric coherence and slope to estimate the SPC height. Interferometric coherence records forest height and dielectric property information from an InSAR perspective and has spatiotemporal consistency with the InSAR-derived DEM. (2) Introduce Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) data to provide more observational information and automatically screen ICESat-2 control points with similar forest and slope conditions in the local area to suppress forest spatial heterogeneity. (3) A weighted least squares criterion was used to solve this model to estimate the SPC height. The results were validated at four test sites using high-precision airborne light detection and ranging (LiDAR) data as a reference. Compared to the InSAR-derived DEM, the accuracy of the sub-canopy topography was improved by nearly 60%, on average. Furthermore, we investigated the necessity of local modeling, confirming the potential of the proposed method for estimating sub-canopy topography by relying only on TanDEM-X and ICESat-2 data.

Imaging of moho topography with conditional generative adversarial network from observed gravity anomalies

Accurate estimation of Moho topography plays a crucial role in understanding Earth’s structure, geodynamic processes, and resource exploration. This study presents a novel approach that utilizes conditional Generative Adversarial Networks (cGAN) to reveal Moho topography based on observed gravity anomalies. Synthetic training datasets of Moho topography were generated using the FFT filtering method due to the scarcity of true datasets. Spherical prism-based forward gravity modeling was employed to evaluate the resulting gravity anomalies. We compared the performance of our developed deep learning algorithm cGAN (conditional Generative Adversarial Networks) with a traditional inversion technique using various synthetic datasets, and a real case study in southern peninsular India, a geologically diverse region comprising ancient continental tectonic blocks. Bott’s inversion scheme was employed as a verification method for the Moho surface estimation using the presented deep learning model. Using spherical prism-based forward gravity modeling, observed gravity anomalies were corrected for multiple factors such as topography, bathymetry, sediments, crustal heterogeneities, and mantle heterogeneities. By removing these effects, we isolated the gravity contribution solely related to pure Moho undulation. The mean Moho depth and density contrast between the crust and mantle were derived from seismic constraints for improving estimation accuracy. The findings demonstrate the potential of the cGAN and spherical prism-based gravity modeling approach in accurately estimating the Moho topography, offering insights into Earth’s subsurface structures and enhancing our understanding of geodynamic processes and resource exploration efforts.

Mantle Structure Beneath the Damara Belt in South‐Central Africa Imaged Using Adaptively Parameterized P‐Wave Tomography

AbstractMany seismic tomography studies have indicated that the African Large Low Velocity Province (LLVP) extends from the lower mantle beneath southern Africa into the upper mantle beneath eastern Africa; however, it has been questioned whether the LLVP structure may also extend to the north or northwest beneath south‐central Africa. Debates regarding the upper mantle structure beneath the Damara Belt contribute to this uncertainty. Some studies suggest the Damara Belt is underlain by thermally perturbed upper mantle; however, other studies indicate the region is not associated with anomalous structure. Here, we use a comprehensive P‐wave travel‐time data set and an adaptive model parameterization to develop a new tomographic model for the Damara Belt and surrounding regions. Our results show that seismically slow structure beneath the Damara Belt is relegated to depths greater than ∼1,200 km, indicating that the LLVP is not significantly affecting this region. However, further to the northeast, the LLVP structure obliquely rises and crosses the mantle transition zone near the Irumide Belt, where it then extends into the upper mantle. The seismic structure beneath the Damara Belt and neighboring areas in our model correlates well with tectonic observations at the surface, including variations in heat flow, the distribution of geothermal features, the locations of rifts, and estimates of dynamic topography.

Old detrital AFT ages from East Greenland do not require plateau erosion

Accurate estimates of past topography are required to reliably reconstruct past ice sheets to infer palaeoclimate. For this reason, understanding erosion rates across East Greenland is crucial to constrain landscape evolution driven by tectonics and climate-dependent erosion rates. Here we analyse published apatite fission-track (AFT) data to constrain the spatial pattern of AFT bedrock ages across the landscape. We compare these bedrock ages with published detrital distribution to highlight ambiguity in the pattern of erosion. In contrast to earlier work, we regress a simple model of exhumation pace through the bedrock ages such that age can vary as a function of both elevation and position. The resulting iso-age surfaces allow us to determine potential source areas for detrital AFT age distributions. We find that old ages observed in detrital distributions are just as likely to be sourced from low-elevation locations that are far from the coast as from high-elevation locations close to the coast. Additional data from lower temperature systems are thus required to make firm conclusions on landscape evolution in the region and distinguish between the two landscape-forming scenarios.

Integrating Hydrography Observations and Geodetic Data for Enhanced Dynamic Topography Estimation

Dynamic topography (DT) refers to the time-varying component of the sea surface height influenced by factors like ocean currents, temperature, and salinity gradients. Accurate estimation of DT is crucial for comprehending oceanic circulation patterns and their impact on climate. This study introduces two approaches to estimating DT: (1) utilizing satellite altimetry to directly observe sea surface height and (2) considering the steric and non-steric components of sea level anomalies. The steric term is calculated using salinity and temperature data obtained from local buoy data, Argo observations, and the World Ocean Atlas model. The non-steric term is calculated using GRACE Satellite gravimetry data. To estimate the assimilated DT, four methods are utilized, including variance component estimation (VCE), Bayesian theory, Kalman filter, and 3D variational (3DVAR). These methods assimilate the two aforementioned schemes. The validity of the estimated DT is assessed by comparing the calculated sea surface current, derived from the obtained DT, with observations from local current meter stations. The results indicate that the VCE method outperforms other methods in determining the final DT. Furthermore, incorporating the steric and non-steric terms of sea level in determining DT in coastal areas enhances the accuracy of estimating sea surface currents.

Simultaneous Estimation of Subcanopy Topography and Forest Height With Single-Baseline Single-Polarization TanDEM-X Interferometric Data Combined With ICESat-2 Data

Simultaneous Estimation of Subcanopy Topography and Forest Height With Single-Baseline Single-Polarization TanDEM-X Interferometric Data Combined With ICESat-2 Data

Estimation of Sea Surface Topography Using Jason-3 Satellite Altimetry Across the Eastern Indonesian Seas in a period of 2016 to 2021

The measurement of sea level height relative to a fixed reference point, such as the geoid or mean sea level, is known as sea surface topography. There are a number of factors that affect it, such as undulation, sea surface height, and tides. Measurement of sea surface topography is very important to understand ocean dynamics. One of the factors that influences the dynamics of the sea is the flow across Indonesia. The Indonesian traffic flow is the mass flow of water from the Pacific Ocean to the Indian Ocean. The difference in pressure between the two oceans is what causes the flow. The purpose of this study is to provide precise information about the sea surface topography in the East Indonesian Sea region during the 2016–2021 period, which is required for various marine applications. This approach involves processing and analysing Jason-3 altimetry satellite data, which provides precise sea level elevation measurements. The data is then used to determine the value of the sea surface topography, which shows the spatial and temporal variations in sea level height in the region. The average sea surface topography value every month from 2016 to 2021 ranges from 1.05m to 1.30m, with the highest value occurring in January 2021 at 1.2542m and the lowest value occurring in February 2016 at 1.0498m. For the accuracy of processing results, data for 2016–2021 is used as training data for model estimation and 2021 as test data. The root mean square error value is 0.0762m. Changes in sea surface topography values tend to increase every year, except for 2019. The decline that occurred in 2019 was caused by the ENSO phenomenon, which had an effect of -0.807 on sea surface topography with an inverse correlation. The findings of this study provide important details about sea surface topography in Eastern Indonesian waters, which can be applied to enhance our understanding of ocean dynamics.

A transient coupled general circulation model (CGCM) simulation of the past 3 million years

Abstract. Driven primarily by variations in the earth's axis wobble, tilt, and orbit eccentricity, our planet experienced massive glacial/interglacial reorganizations of climate and atmospheric CO2 concentrations during the Pleistocene (2.58 million years ago (Ma)–11.7 thousand years ago (ka)). Even after decades of research, the underlying climate response mechanisms to these astronomical forcings have not been fully understood. To further quantify the sensitivity of the earth system to orbital-scale forcings, we conducted an unprecedented quasi-continuous coupled general climate model simulation with the Community Earth System Model version 1.2 (CESM1.2, ∼3.75∘ horizontal resolution), which covers the climatic history of the past 3 million years (3 Myr). In addition to the astronomical insolation changes, CESM1.2 is forced by estimates of CO2 and ice-sheet topography which were obtained from a simulation previously conducted with the CLIMBER-2 earth system model of intermediate complexity. Our 3 Ma simulation consists of 42 transient interglacial/glacial simulation chunks, which were partly run in parallel to save computing time. The chunks were subsequently merged, accounting for spin-up and overlap effects to yield a quasi-continuous trajectory. The computer model data were compared against a plethora of paleo-proxy data and large-scale climate reconstructions. For the period from the Mid-Pleistocene Transition (MPT, ∼1 Ma) to the late Pleistocene we find good agreement between simulated and reconstructed temperatures in terms of phase and amplitude (−5.7 ∘C temperature difference between Last Glacial Maximum and Holocene). For the earlier part (3–1 Ma), differences in orbital-scale variability occur between model simulation and the reconstructions, indicating potential biases in the applied CO2 forcing. Our model-proxy data comparison also extends to the westerlies, which show unexpectedly large variance on precessional timescales, and hydroclimate variables in major monsoon regions. Eccentricity-modulated precessional variability is also responsible for the simulated changes in the amplitude and flavors of the El Niño–Southern Oscillation. We further identify two major modes of planetary energy transport, which played a crucial role in Pleistocene climate variability: the first obliquity and CO2-driven mode is linked to changes in the Equator-to-pole temperature gradient; the second mode regulates the interhemispheric heat imbalance in unison with the eccentricity-modulated precession cycle. During the MPT, a pronounced qualitative shift occurs in the second mode of planetary energy transport: the post-MPT eccentricity-paced variability synchronizes with the CO2 forced signal. This synchronized feature is coherent with changes in global atmospheric and ocean circulations, which might contribute to an intensification of glacial cycle feedbacks and amplitudes. Comparison of this paleo-simulation with greenhouse warming simulations reveals that for an RCP8.5 greenhouse gas emission scenario, the projected global mean surface temperature changes over the next 7 decades would be comparable to the late Pleistocene glacial-interglacial range; but the anthropogenic warming rate will exceed any previous ones by a factor of ∼100.

A mountain ridge model for quantifying oblique mountain wave propagation and distribution

Abstract. Following the current understanding of gravity waves (GWs) and especially mountain waves (MWs), they have a high potential for horizontal propagation from their source. This horizontal propagation and therefore the transport of energy is usually not well represented in MW parameterizations of numerical weather prediction and general circulation models. In this study, we present a mountain wave model (MWM) for the quantification of horizontal propagation of orographic gravity waves. This model determines MW source locations from topography data and estimates MW parameters from a fit of idealized Gaussian-shaped mountains to the elevation. Propagation and refraction of these MWs in the atmosphere are modeled using the Gravity-wave Regional Or Global Ray Tracer (GROGRAT). Ray tracing of each MW individually allows for an estimation of momentum transport due to both vertical and horizontal propagation. The MWM is a capable tool for the analysis of MW propagation and global MW activity and can support the understanding of observations and improvement of MW parameterizations in GCMs. This study presents the model itself and gives validations of MW-induced temperature perturbations to ECMWF Integrated Forecast System (IFS) numerical weather prediction data and estimations of GW momentum flux (GWMF) compared to HIgh Resolution Dynamics Limb Sounder (HIRDLS) satellite observations. The MWM is capable of reproducing the general features and amplitudes of both of these data sets and, in addition, is used to explain some observational features by investigating MW parameters along their trajectories.

Event-related potential and oscillatory cortical activities of artistic methodology in information visualization design in human–computer interface

Data visualization is essential for interface design in human–computer interaction. With the advances in information technology, artists and designers have used creative thoughts and artistic methodologies in information design and data visualization. We design a method that minimizes readers’ time to memorize data or information through event-related potential (ERP), electroencephalography (EEG) topography, brain connectivity, and volume source estimation analyses. We developed four information visualization interfaces by importing career data from the World Bank Data: integrated generative artistic design (IGAD), separated generative artistic design (SGAD), integrated line chart design (ILCD), and separated line chart design (SLCD). We then compared subjective evaluations and brain electrophysiological biofeedback on user experience (UX) with the aesthetics and usability of user interface design. We enroled 52 healthy subjects (26 men and 26 women) for psychological or physiological differences. The line chart wireframe design and initial animation arousal were found to be more acceptable than the generative style through P200 and N200 analyses. Meanwhile, late positive complex (LPC) and subjective evaluations showed that SGAD is the most beautiful interface for information visualization design. LPC showed generative art as more appealing than the line chart in the left and middle brain regions from the frontal-to-central lobes. We found that decreased delta, theta, and alpha and increased gamma lead to significant memory enhancement in ILCD compared to IGAD. Men and women were found to have different relative power activations in delta, theta, alpha, and gamma. The sex-based difference in visual working memory in theta and gamma connectivity in frontal-to-central brain regions occurred in the ILCD interface design. The difference between delta activation on line chart design and theta and gamma activated in the frontal-parietal-occipital brain in SGAD helped reveal how different designs influence the same data result through connectivity analysis. Volume estimation revealed that the activations vary between brain regions across the one-minute analysis. The subjective evaluation produced different results from LPC objective evaluation aesthetic judgement but the same results in SGAD interface design. We found that memory performance decreases from SLCD, ILCD, SGAD, to IGAD. We also found that most people quickly memorized incrementing and decrementing data trends, which can help future human–computer interface and information visualization design.

Underlying topography and forest height estimation from SAR tomography based on a nonparametric spectrum estimation method with low sidelobes

ABSTRACT The underlying topography and forest height play an indispensable role in many fields, including geomorphology, civil engineering construction, forest investigation, and the modeling of natural disasters. As a new microwave remote sensing technology with three-dimensional imaging capability, synthetic aperture radar (SAR) tomography (TomoSAR) has already been proven to be an important tool for underlying topography and forest height estimation. Many spectrum estimation methods have now been proposed for TomoSAR. However, most of the commonly used methods are susceptible to noise and inevitably produce sidelobes, resulting in a reduced accuracy for the inversion of forest structural parameters. In this paper, to solve this problem, a nonparametric spectrum estimation method with low sidelobes – the G-Pisarenko method – is introduced. This method performs a logarithmic operation on the covariance matrix to obtain the main scattering characteristics of the objects of interest while suppressing the noise as much as possible. The effectiveness of the proposed method is demonstrated by the use of both simulated data and P-band airborne SAR data from a tropical forest region in Gabon, Africa. The results show that the proposed method can reduce the sidelobes and improve the estimation accuracy for the underlying topography and forest height.

Control of interaction force in constant-height contact mode atomic force microscopy

Contact mode is a versatile and widely used technique for imaging samples using the Atomic Force Microscope (AFM). When contact mode imaging is performed in constant-height mode, it enables linear and faster response but leads to uncontrolled tip-sample forces. Here, a control strategy based on magnetic actuation is proposed to achieve high-bandwidth control of the tip-sample forces in constant-height contact mode AFM. A magnetic particle attached to the AFM probe is actuated by an external solenoid and employed for force regulation. A quasi-static model has been proposed and employed to develop the control strategy. Likewise, the contact natural-frequency, which decides the limit of achievable speed, has been shown to be significantly higher relative to the free probe and to be relatively insensitive to the particle size. Subsequently, a setup is developed to validate the control strategy and demonstrate reduction of tip-sample force variation by over a factor of 12 compared to conventional constant-height mode operation. Likewise, in comparison with conventional contact mode AFM, an improvement of linearity by over a factor of 9 and improvement in response speed by a factor of 100 have been demonstrated while imaging hard samples. The system has been shown to image topography at speeds of 2.44 frames per second while regulating the interaction force. Finally, the stiffness of a sample has also been characterized using the developed system and simultaneous estimation of topography has also been demonstrated. They are shown to agree well with theoretical expectations.

Earth's Isostatic and Dynamic Topography—A Critical Perspective

AbstractEarth's topography arises from the linear superposition of isostatic and dynamic contributions. The isostatic contribution reflects the distribution of thickness and density of the crust overlying a static, non‐convecting mantle. We argue that isostatic topography should be limited to the crust, thereby delimiting all sources for dynamic topography below the Moho. Dynamic topography is the component of the topography produced by normal stresses acting on the Moho that deflect the isostatic topography away from crustal isostatic equilibrium largely as a consequence of mantle flow dynamics. These normal stresses arise from pressure variations and vertical gradients of the radial flow in the convecting mantle. The best estimate of dynamic topography is from the residual topography, which is the difference between observed topography and crustal isostatic topography. Dynamic and residual topography are the same. It is clear that thermal anomalies horizontally advected by plate motions would not exist if the mantle were not convecting, therefore their contribution to topography is inherently dynamic in origin. The global integral of dynamic topography that encompasses all non‐crustal buoyancy sources is demonstrated to be equal to zero. It follows that mantle convection cannot change the mean radius or mean elevation of the Earth. Since changes in ocean basin volume driven by changes in mean depth of the oceans are inherently part of dynamic topography, thereby requiring that continental elevations must also change, such that the global integral of these perturbations must also be equal to zero. This constraint has important implications for global long‐term sea level and the stratigraphic record, among other features of the Earth system impacted by changes in Earth's dynamic topography.

Marine Gravimetry and Its Improvements to Seafloor Topography Estimation in the Southwestern Coastal Area of the Baltic Sea

Marine gravimetry provides high-quality gravity measurements, particularly in coastal areas. After the update of new sensors in GFZ’s air-marine gravimeter Chekan-AM, gravimetry measurements showed a significant improvement from the first new campaign DENEB2017 with an accuracy of 0.3/2=0.21 mGal @ 1 km along the tracks, which is at the highest accuracy level of marine gravimetry. Then, these measurements were used to assess gravity data derived from satellite altimetry (about 3 mGal) and a new finding is that a bias of −1.5 mGal exists in the study area. Additionally, ship soundings were used to assess existing seafloor topography models. We found that the accuracy of SRTM model and SIO model is at a level of 2 m, while the accuracy of the regional model EMODnet reaches the lever of sub-meters. Furthermore, a bias of 0.7 m exists and jumps above 5 m in the SRTM model near the coast of Sweden. Finally, new combined gravity anomalies with sounding data are used to reveal the fine structure of ocean topography. Our estimated seafloor topography model is more accurate than existing digital elevation data sets such as EMODnet, SRTM and SIO models and, furthermore, shows some more detailed structure of seafloor topography. The marine gravimetry and sounding measurements as well as the estimated seafloor topography are crucial for future geoid determination, 3D-navigation and resource exploration in the Baltic Sea.