Search

Seismic moment accumulation rate is a fundamental parameter for assessing seismic hazard. It can be estimated geodetically from either fault‐based modeling, or strain rate‐based calculations, where fault‐based models largely depend on the rheological layering and the number of faults. The strain‐rate method depends on an unknown (Kostrov) thickness used to convert strain rate into moment rate. In Part 1 of this study, we use three published fault‐based models from southern California to establish the value of the Kostrov thickness such that the total moment from the strain‐rate approach, calculated from the fault model‐predicted strain rate, matches the fault‐based approach. Constrained thickness estimates of 7.3, 9.7, and 11.5 km (6.4–13.0 km, including uncertainties) suggest that the 11 km value used in previous studies may be too large and a lower value may be more accurate. In Part 2 we use calibrated values of Kostrov thickness, along with the latest compilation of GNSS velocity data, to partition moment rate into on‐fault and off‐fault moment rate, where off‐fault varies from 32%–43% of the total moment rate. The largest uncertainty is related to the method used to interpolate sparse GNSS data. Lastly, we compare our estimates of total moment rate (mean: 2.13 ± 0.42 × 1019 Nm/yr) with the historical seismic catalog. Results suggest that including uncertainties in Kostrov thickness brings fault‐based geodetic moment rate closer to the seismic moment release (particularly when aseismic afterslip is accounted for), while the (uncertain) values of off‐fault moment rate push geodetic moment rates to be larger than seismic moment rates. Plain Language Summary: One commonly‐used method of evaluating the future earthquake potential of a known fault is to estimate the moment rate accumulating along that fault. In other words, one can estimate how much earthquake energy is accumulating over time, caused by the forces that propel plate boundary motion along fault systems. Two ways to do this are: (a) using detailed fault models of plate boundary motion, which are constrained by satellite‐based measurements of surface motion or (b) using those satellite‐based measurements directly through the estimation of regional strain rate. Here, we investigate the uncertainties present in the latter method by comparing previously published estimates of fault‐based moment rate, with their equivalent strain rate‐based versions. We estimate uncertainties related to the strain‐rate based method, which can lead to an over‐ or under‐estimation of earthquake potential on a given fault system. In addition, we explore the off‐fault moment rate accumulation in southern California which can help identify areas of enhanced moment accumulation and thus increased seismic hazard. Key Points: Previously published fault‐based models can constrain a Kostrov thickness range of 6.4–13.0 km, including uncertaintiesDifferences between strain rate moment estimates and fault‐based estimates suggest 32%–43% off‐fault moment rate in southern CaliforniaThe largest uncertainty in estimating moment rate from strain rate is the degree of smoothing used to estimate strain from geodetic data [ABSTRACT FROM AUTHOR]

A coherent portrayal of global bathymetry requires that depths are inferred between sparsely distributed direct depth measurements. Depths can be interpolated in the gaps using alternate information such as satellite‐derived gravity and a mapping from gravity to depth. We designed and trained a neural network on a collection of 50 million depth soundings to predict bathymetry globally using gravity anomalies. We find the best result is achieved by pre‐filtering depth and gravity in accordance with isostatic admittance theory described in previous predicted depth studies. When training the model, if the training and testing split is a random partition at the same resolution as the data, the training and testing sets will not be independent, and model misfit is underestimated. We solve this problem by partitioning the training and testing set with geographic bins. Our final predicted depth model improves on old predicted depth model RMSE by 16%, from 165 to 138 m. Among constrained grid cells, 80% of the predicted values are within 128 m of the true value. Improvements to this model will continue with additional depth measurements, but predictions at higher spatial resolution, being limited by upward continuation of gravity, should not be attempted with this method. Plain Language Summary: Only a fraction of the seafloor has been mapped by shipboard measurements. In the unmapped regions of the ocean, we must estimate the depth of the seafloor using information from the Earth's gravity field. Models predicting seafloor depth using gravity typically determine the linear relationship between gravity and depth in some regions and use the established relationships to make global predicted depth maps. Here, we describe a new method for predicting depth globally using gravity, decades of shipboard depth measurements, and a neural network regression. Ultimately, our model shows a clear improvement over the reference model. Key Points: We present a new method for global bathymetry prediction using a machine learning algorithmThe new predicted depth model improves on the reference model by all error metrics [ABSTRACT FROM AUTHOR]

The mid‐ocean ridge system comprises a series of spreading ridges, transform faults, propagating ridges, and other non‐transform offsets. Transform faults remain stable for millions of years leaving long linear scars, or fracture zones, on older seafloor. Propagating ridges migrate in the ridge parallel direction leaving V‐shaped or W‐shaped scars on older seafloor. Vertical gravity gradient maps can now resolve the details of the ridge segmentation. For slow‐ and intermediate‐spreading ridges, there appears to be an offset length threshold above which adjacent ridges do not propagate so remain as stable transform faults. We propose this threshold is due to the yield strength of the lithosphere, and we develop a model framework based on a force balance wherein forces driving propagation must exceed the integrated shear strength of the offset zone. We apply this model framework to 5 major propagating ridges, 55 seesaw propagating ridges, and 69 transform faults. The model correctly predicts the migration of 4 out of 5 major propagating ridges and the stability of transform faults, but the results for seesaw propagators are less accurate. Model predictions for direction of ridge propagation are mixed as well. This model framework simplifies deformation in the shear zone, but can possibly explain why non‐transform deformation is preferred at short offsets. Plain Language Summary: Mid‐ocean ridges are constructive plate boundaries where new crust is created. In map view, the system resembles a stair‐step configuration of alternating spreading ridges and ridge offsets. Some ridges and offsets, typically large ones, remain fixed and maintain their plan‐view shape over many millions of years, while other ridges, usually those bound by shorter offsets, may slowly grow and shrink–such behavior is revealed in maps of the seafloor. The different behavior is possibly due to the material strength of the oceanic crust and upper mantle which, if great, will inhibit ridge growth. To test our hypothesis, we estimate the total material strength at identified ridge offsets and compare this to an estimate of forces contributing to ridge growth. Our estimates can explain why large offsets maintain their shape, and may explain why short offsets do not and allow some segments to grow and shrink. Key Points: Non‐transform offsets at slower‐spreading ridges rarely propagate when offsets exceed 30 km, possibly a result of lithospheric strengthWe develop a model framework that balances material strength at ridge offsets and forces driving ridge propagationGreater strength of the lithosphere as ridge offset increases may limit ridge propagation [ABSTRACT FROM AUTHOR]

Seamounts are isolated elevations in the seafloor with circular or elliptical plans, comparatively steep slopes, and relatively small summit area (Menard, 1964). The vertical gravity gradient (VGG), which is the curvature of the ocean surface topography derived from satellite altimeter measurements, has been used to map the global distribution of seamounts (Kim & Wessel, 2011, https://doi.org/10.1111/j.1365-246x.2011.05076.x). We used the latest grid of VGG to update and refine the global seamount catalog; we identified 19,325 new seamounts, expanding a previously published catalog having 24,643 seamounts. Seven hundred thirty‐nine well‐surveyed seamounts, having heights ranging from 421 to 2,500 m, were used to estimate the typical radially symmetric seamount morphology. First, an Empirical Orthogonal Function (EOF) analysis was used to demonstrate that these small seamounts have a basal radius that is linearly related to their height—their shapes are scale invariant. Two methods were then used to compute this characteristic base to height ratio: an average Gaussian fit to the stack of all profiles and an individual Gaussian fit for each seamount in the sample. The first method combined the radial normalized height data from all 739 seamounts to form median and median‐absolute deviation. These data were fit by a 2‐parameter Gaussian model that explained 99.82% of the variance. The second method used the Gaussian function to individually model each seamount in the sample and further establish the Gaussian model. Using this characteristic Gaussian shape we show that VGG can be used to estimate the height of small seamounts to an accuracy of ∼270 m. Key Points: We used the latest vertical gravity gradient maps to update and refine a global seamount catalog, finding 19,325 new seamountsSmaller seamounts (<2,500 m tall) having good bathymetry coverage (739) were modeled with a radially symmetric Gaussian functionTwo modeling approaches show that smaller seamounts have a sigma to height ratio of 2.4 which agrees with an earlier study by Smith (1988) [ABSTRACT FROM AUTHOR]

The 2019 Ridgecrest conjugate Mw6.4 and Mw7.1 events resulted in several meters of strike‐slip and dip‐slip along an intricate rupture, extending from the surface down to 15 km. Now with >2 years of post‐rupture observations, we utilize these results to better understand vertical postseismic deformation from the Ridgecrest sequence and illuminate the emerging significance of vertical earthquake cycle deformation data. We determine the cumulative vertical displacement observed by the continuous GNSS network since Ridgecrest, which requires additional time series analyses to adequately resolve vertical deformation compared to the horizontal. Using a Maxwell‐type viscoelastic relaxation model, with a best fit time‐averaged asthenosphere viscosity of 4e17 Pa·s and a laterally heterogeneous lithosphere, we find that viscoelastic relaxation accounts for a majority of the cumulative vertical deformation at Ridgecrest and strongly controls far‐field observations in all north‐east‐up components. The viscoelastic model alone generally underpredicts deformation from GNSS and the remaining nonviscoelastic displacement is most prominent in the horizontal near‐field (−16 to 19 mm), revealing a deformation pattern matching the coseismic observations. This suggests that multiple deformation mechanisms are contributing to Ridgecrest's postseismic displacement, where afterslip likely dominates the near‐field while viscoelastic relaxation controls the far‐field. Similar deformation at individual GNSS stations has been observed for past earthquakes and additionally reveals long‐term transient viscosity over several years. Moreover, the greater temporal and spatial resolution of the GNSS array for Ridgecrest will help resolve the evolution of deformation for the entire network of observations as regional postseismic deformation persists for the next several years. Plain Language Summary: The 2019 Ridgecrest earthquake sequence is one the most well observed seismic events in California's history. We take advantage of the unprecedented amount of satellite observations and previous modeling efforts to better understand the vertical post‐earthquake signal resulting from Ridgecrest. We compare cumulative surface displacement since the earthquakes, derived from Global Navigation Satellite System (GNSS) point observations to a viscoelastic model, which allows us to quantify postseismic deformation due to viscoelastic (i.e., viscoelastic relaxation) and nonviscoelastic (i.e., afterslip and poroelastic rebound) contributions. We find that viscoelastic relaxation accounted for a majority of the currently observed vertical deformation since Ridgecrest and strongly controlled far‐field observations in all north‐east‐up components, emphasizing the importance of utilizing both horizontal and vertical observations when developing earthquake cycle models. The remaining observed deformation insinuates that multiple postseismic deformation mechanisms, and thus nonviscoelastic contributions, are also needed to reproduce vertical postseismic displacement, where afterslip is the most likely mechanism here. We additionally see similarities between Ridgecrest and past large regional earthquakes and emphasize the ability of our deformation model, derived from 3D surface observations, to provide important insight on crustal parameters and the characteristics of the seismic region for past, current, and future events. Key Points: Viscoelastic relaxation accounts for the majority of vertical postseismic displacement >1 year after Ridgecrest and controls the far‐fieldVertical observations provide important constraints on time‐averaged and transient upper mantle viscosityAdditional time series analysis techniques are needed to produce coherent vertical deformation patterns compared to the horizontal [ABSTRACT FROM AUTHOR]

To date, ∼20% of the ocean floor has been surveyed by ships at a spatial resolution of 400 m or better. The remaining 80% has depth predicted from satellite altimeter‐derived gravity measurements at a relatively low resolution. There are many remote ocean areas in the southern hemisphere that will not be completely mapped at 400 m resolution during this decade. This study is focused on the development of synthetic bathymetry to fill the gaps. There are two types of seafloor features that are not typically well resolved by satellite gravity; abyssal hills and small seamounts (<2.5 km tall). We generate synthetic realizations of abyssal hills by combining the measured statistical properties of mapped abyssal hills with regional geology including fossil spreading rate/orientation, rms height from satellite gravity, and sediment thickness. With recent improvements in accuracy and resolution, it is now possible to detect all seamounts taller than about 800 m in satellite‐derived gravity and their location can be determined to an accuracy of better than 1 km. However, the width of the gravity anomaly is much greater than the actual width of the seamount so the seamount predicted from gravity will underestimate the true seamount height and overestimate its base dimension. In this study, we use the amplitude of the vertical gravity gradient (VGG) to estimate the mass of the seamount and then use their characteristic shape, based on well‐surveyed seamounts, to replace the smooth‐predicted seamount with a seamount having a more realistic shape. Plain Language Summary: The floor of the deep ocean remains as the last uncharted frontier in the inner solar system. The deep seawater (>1,000 m) prevent any type of exploration from conventional satellite remote sensing. Echosounders aboard large vessels have mapped about 20% of the seafloor, however, vast areas in the southern hemisphere will not be mapped in our lifetimes. The deep ocean floor has more than 90% of the active volcanoes; hydrothermal circulation of seawater through the crust of the seafloor spreading ridges replenishes the nutrients needed for life on Earth. This study is an effort to fill the large gaps in seafloor coverage by creating a synthetic abyssal hill fabric using geological information such as the age of the seafloor, the spreading rate and direction when it formed, and the thickness of the sediments covering the original topography. In addition, we use the latest satellite‐derived gravity to estimate the locations and shapes of about 20,000 uncharted seamounts. The combination of mapped (20%) and synthetic (80%) topography is useful for modeling ocean circulation and ocean tides although it may give a false impression that 100% of the seafloor has been mapped. Key Points: ∼20% of the ocean floor topography has been surveyed by ships, the remaining 80% is predicted by satellite altimetryWe increased the resolution of the predicted depth using spectral properties of abyssal hills and the characteristic shapes of seamountsWe estimate the height and radius of 19,000 uncharted seamounts [ABSTRACT FROM AUTHOR]

Measuring crustal strain and seismic moment accumulation, is crucial for understanding the growth and distribution of seismic hazards along major fault systems. Here, we develop a methodology to integrate 4.5 years (2015–2019.5) of Sentinel‐1 Interferometric Synthetic Aperture Radar (InSAR) and continuous Global Navigation Satellite System (GNSS) time series to achieve 6 to 12‐day sampling of surface displacements at ∼500 m spatial resolution over the entire San Andreas fault system. Numerous interesting deformation signals are identified with this product (video link: https://www.youtube.com/watch?v=SxNLQKmHWpY). We decompose the line‐of‐sight InSAR displacements into three dimensions by combining the deformation azimuth from a GNSS‐derived interseismic fault model. We then construct strain rate maps using a smoothing interpolator with constraints from elasticity. The resulting deformation field reveals a wide array of crustal deformation processes including, on‐ and off‐fault secular and transient tectonic deformation, creep rates on all the major faults, and vertical signals associated with hydrological processes. The strain rate maps show significant off‐fault components that were not captured by GNSS‐only models. These results are important in assessing the seismic hazard in the region. Plain Language Summary: Seismic hazard models rely on accurate measurements of small motion over large areas on the Earth's crust. Traditional geodetic models based on Global Navigation Satellite System (GNSS) data cannot resolve small scale deformation patterns, mainly due to expensive and limited station deployment. Interferometric Synthetic Aperture Radar (InSAR) has become the emerging tool for mapping surface deformation, with its advantages of low‐cost and full‐coverage. Yet InSAR measurements, compared to GNSS, come with larger biases from the atmospheric noise, especially over length scales greater than 80 km. Here, we combined the two methods to resolve fine spatial scales and achieve high accuracy. Our results are presented as deformation time‐series over the entire San Andreas fault system (video link: https://www.youtube.com/watch?v=SxNLQKmHWpY). From these deformation time series, we have estimated fault creep rates and strain accumulation. One important finding is that there is significant off‐fault strain, though we suspect this is mainly due to hydrological processes. These results will advance our knowledge of the earthquake cycle, strain/moment accumulation, and the associated seismic hazards. Key Points: A practical approach is developed to integrate Sentinel‐1 Interferometric Synthetic Aperture Radar and Global Navigation Satellite System time‐series over the entire San Andreas fault systemThe product is used to estimate fault creep and three components of horizontal crustal strain which shows notable off‐fault portionChallenges remain in separating tectonic and hydrologic sources and whether hydrologic strain will increase seismic hazards [ABSTRACT FROM AUTHOR]

Rheologic variations in the Earth's crust (like elastic plate thickness [EPT] or crustal rigidity) modulate the rate at which seismic moment accumulates for potentially hazardous faults of the San Andreas Fault System (SAFS). To quantify rates of seismic moment accumulation, Global Navigation Satellite Systems, and Interferometric Synthetic Aperture Radar data were used to constrain surface deformation rates of a four‐dimensional viscoelastic deformation model that incorporates rheological variations spanning a 900 km section of the SAFS. Lateral variations in EPT, estimated from surface heat flow and seismic depth to the lithosphere‐asthenosphere boundary, were converted to lateral variations in rigidity and then used to solve for seismic moment accumulation rates on 32 fault segments. We find a cluster of elevated seismic moment rates (11–20 × 1015 Nm year−1 km−1) along the main SAFS trace spanning the historical Mw 7.9 1857 Fort Tejon earthquake rupture length; present‐day seismic moment magnitude on these segments ranges from Mw 7.2–7.6. We also find that the average plate thickness in the Salton Trough is reduced to only 60% of the regional average, which results in a ∼60% decrease in moment accumulation rate along the Imperial fault. Likewise, a 30% increase of average plate thickness results in at least a ∼30% increase in moment rate and even larger increases are identified in regions of complex plate heterogeneity. These results emphasize the importance of considering rheological variations when estimating seismic hazard, suggesting that meaningful changes in seismic moment accumulation are revealed when considering spatial variations in crustal rheology. Plain Language Summary: The earthquake potential of faults (i.e., moment magnitude) is determined by a fault's deep slip rate, the area of the fault, the time since the last major rupture, and the average crustal rigidity surrounding the fault. Here, we explore how earthquake potential is affected by variations in crustal rheology along the San Andreas Fault System (SAFS). We use measurements of surface heat flow and seismic velocity to estimate lateral variations in crustal rigidity and tectonic plate thickness. We then estimate seismic potential rate for 32 segments of the SAFS using surface deformation measurements from Global Navigation Satellite Systems and radar interferometry as constraints. We find high seismic potential rates along several segments of the San Andreas fault that also participated in the Mw 7.9 1857 Fort Tejon earthquake and have not ruptured in over 160 years. We also find that the change in earthquake potential scales almost linearly with plate thickness variations surrounding the fault. For example, the Salton Trough region has a plate thickness that is only 60% of the regional average, which lowers the earthquake potential along the Imperial fault by roughly 60%. These findings suggest that variations in crustal rheology have an important impact on earthquake magnitude forecasts. Key Points: Variable rheology was implemented in a four‐dimensional viscoelastic model to explore earthquake cycle deformation of the San Andreas Fault SystemElastic plate thickness variations are estimated from heat flow and depth to the lithosphere‐asthenosphere boundaryVariations in earthquake magnitude forecasts should be carefully considered when implementing variable crustal parameters [ABSTRACT FROM AUTHOR]

The segmentation of mid‐ocean ridges is a defining characteristic of seafloor spreading, yet some tectonic processes operating at segment boundaries remain poorly understood. Here, we analyze new satellite‐derived vertical gravity gradient (VGG) data, which reveal an abundance of off‐axis seafloor features that are oblique to ridges and transform faults and thus reflect the occurrence of ridge propagation at some segment boundaries. However, unlike many propagating ridges, these features commonly reverse direction leaving W‐shaped signatures in the seafloor which we refer to as "seesaw propagators" (SSPs). Using the VGG, we have documented these globally and find that: (1) SSPs are ubiquitous on seafloor that formed at half spreading rates between 10 and 40 mm year−1 and their total length is about 1/3 that of fracture zones. (2) The lithospheric age offset across SSPs (0–2.5 Ma) is less than transform faults (2–10 Ma), which likely reflects a rheological threshold, whereby only young and weak lithosphere allows for "non‐rigid" SSP behavior. (3) Isostatic modeling of well‐surveyed SSPs confirms that they formed on young and thin (3 km) lithosphere. (4) The directional changes of SSPs in both time and space appear largely uncorrelated and cannot be explained by previous regional‐scale models invoked to explain unidirectional ridge propagation and thus require a different driving force. Key Points: The process of "seesaw" ridge propagation is ubiquitous on seafloor that forms at half spreading rates between 10 and 40 mm yr−1Globally, the total length of seesaw propagators is approximately one‐third that of fracture zonesDriving mechanisms previously proposed to explain unidirectional ridge propagation fail account for the behavior of seesaw propagators [ABSTRACT FROM AUTHOR]

Our understanding of plate boundary deformation has been enhanced by transient signals observed against the backdrop of time‐independent secular motions. We make use of a new analysis of displacement time series from about 1,000 continuous Global Positioning System (GPS) stations in California from 1999 to 2018 to distinguish tectonic and nontectonic transients from secular motion. A primary objective is to define a high‐resolution three‐dimensional reference frame (datum) for California that can be rapidly maintained with geodetic data to accommodate both secular and time‐dependent motions. To this end, we compare the displacements to those predicted by a horizontal secular fault slip model for the region and construct displacement and strain rate fields. Over the past 19 years, California has experienced 19 geodetically detectable earthquakes and widespread postseismic deformation. We observe postseismic strain rate variations as large as 1,000 nstrain/year with moment releases equivalent up to an Mw6.8 earthquake. We find significant secular differences up to 10 mm/year with the fault slip model, from the Mendocino Triple Junction to the southern Cascadia subduction zone, the northern Basin and Range, and the Santa Barbara channel. Secular vertical uplift is observed across the Transverse Ranges, Coastal Ranges, Sierra Nevada, as well as large‐scale postseismic uplift after the 1999 Mw7.1 Hector Mine and 2010 Mw7.2 El Mayor‐Cucapah earthquakes. We also identify areas of vertical land motions due to anthropogenic, natural, and magmatic processes. Finally, we demonstrate the utility of the kinematic datum by improving the accuracy of high‐spatial‐resolution 12‐day repeat‐cycle Sentinel‐1 Interferometric Synthetic Aperture Radar displacement and velocity maps. Key Points: Displacement fields from 19 years of GPS data in California are preferable to velocity fields in characterizing 3‐D transient motionsMapping of secular and transient motions at high temporal and spatial resolution is required for maintaining a 3‐D kinematic reference frameDisplacement fields derived from continuous GPS and InSAR provide 200‐m resolution with millimeter‐level velocity uncertainties [ABSTRACT FROM AUTHOR]

Observations of shallow fault creep reveal increasingly complex time‐dependent slip histories that include quasi‐steady creep and triggered as well as spontaneous accelerated slip events. Here we report a recent slow slip event on the southern San Andreas fault triggered by the 2017 Mw8.2 Chiapas (Mexico) earthquake that occurred 3,000 km away. Geodetic and geologic observations indicate that surface slip on the order of 10 mm occurred on a 40‐km‐long section of the southern San Andreas fault between the Mecca Hills and Bombay Beach, starting minutes after the Chiapas earthquake and continuing for more than a year. Both the magnitude and the depth extent of creep vary along strike. We derive a high‐resolution map of surface displacements by combining Sentinel‐1 Interferometric Synthetic Aperture Radar acquisitions from different lines of sight. Interferometric Synthetic Aperture Radar‐derived displacements are in good agreement with the creepmeter data and field mapping of surface offsets. Inversions of surface displacement data using dislocation models indicate that the highest amplitudes of surface slip are associated with shallow (<1 km) transient slip. We performed 2‐D simulations of shallow creep on a strike‐slip fault obeying rate‐and‐state friction to constrain frictional properties of the top few kilometers of the upper crust that can produce the observed behavior. Key Points: We present geodetic and geologic observations of slow slip on the southern SAF triggered by the 2017 Chiapas (Mexico) earthquakeThe slow slip event produced surface offsets on the order of 5–10 mm, with significant variations along strikeWe interpret the observed complexity in shallow fault slip in the context of rate‐and‐state friction models [ABSTRACT FROM AUTHOR]

The existing bathymetry map of the Arctic is a compilation of ship soundings and digitized contours. Due to the presence of all‐year sea ice, costly operations and political restrictions, dense and full coverage of the Arctic is not possible, leaving huge gaps between the existing surveys. In this paper, we make use of the existing Arctic bathymetry IBCAOv3 and invert Arctic bathymetry from the recent altimetric gravity model DTU17, whose accuracy is improved significantly with revised data processing strategy. The long and short wavelength components are preserved from IBCAOv3. The band‐pass‐filtering function proposed by Smith and Sandwell (1994, https://doi.org/10.1029/94JB00988) is adapted for the Arctic by reducing the cutoff wavelength. The predicted bathymetry is within 100 m on 85.8% of the grid nodes, when compared to the IBCAOv3. The consistency of the prediction is validated with two independent profiles from Healy cruises conducted in 2016 over the Chukchi Cap. A questionable valley in the IBCAOv3 is detected with gravity and at this spot, bathymetry predicted from gravity is consistent with independent multibeam soundings. The gravity‐inverted bathymetry could be combined with ship soundings for the next generation of Arctic bathymetry map. Key Points: The first ever Arctic bathymetry inverted from DTU17 gravity anomalies is presentedGood agreement between the ship soundings and inverted bathymetry is achieved [ABSTRACT FROM AUTHOR]