scholarly journals Evaluation of ocean circulation models in the computation of the mean dynamic topography for geodetic applications. Case study in the Greek seas

2019 ◽  
Vol 9 (1) ◽  
pp. 154-173
Author(s):  
I. Mintourakis ◽  
G. Panou ◽  
D. Paradissis

Abstract Precise knowledge of the oceanic Mean Dynamic Topography (MDT) is crucial for a number of geodetic applications, such as vertical datum unification and marine geoid modelling. The lack of gravity surveys over many regions of the Greek seas and the incapacity of the space borne gradiometry/gravity missions to resolve the small and medium wavelengths of the geoid led to the investigation of the oceanographic approach for computing the MDT. We compute two new regional MDT surfaces after averaging, for given epochs, the periodic gridded solutions of the Dynamic Ocean Topography (DOT) provided by two ocean circulation models. These newly developed regional MDT surfaces are compared to three state-of-theart models, which represent the oceanographic, the geodetic and the mixed oceanographic/geodetic approaches in the implementation of the MDT, respectively. Based on these comparisons, we discuss the differences between the three approaches for the case study area and we present some valuable findings regarding the computation of the regional MDT. Furthermore, in order to have an estimate of the precision of the oceanographic approach, we apply extensive evaluation tests on the ability of the two regional ocean circulation models to track the sea level variations by comparing their solutions to tide gauge records and satellite altimetry Sea Level Anomalies (SLA) data. The overall findings support the claim that, for the computation of the MDT surface due to the lack of geodetic data and to limitations of the Global Geopotential Models (GGMs) in the case study area, the oceanographic approach is preferable over the geodetic or the mixed oceano-graphic/geodetic approaches.

2016 ◽  
Vol 29 (13) ◽  
pp. 4801-4816 ◽  
Author(s):  
Christopher G. Piecuch ◽  
Sönke Dangendorf ◽  
Rui M. Ponte ◽  
Marta Marcos

Abstract Understanding the relationship between coastal sea level and the variable ocean circulation is crucial for interpreting tide gauge records and projecting sea level rise. In this study, annual sea level records (adjusted for the inverted barometer effect) from tide gauges along the North American northeast coast over 1980–2010 are compared to a set of data-assimilating ocean reanalysis products as well as a global barotropic model solution forced with wind stress and barometric pressure. Correspondence between models and data depends strongly on model and location. At sites north of Cape Hatteras, the barotropic model shows as much (if not more) skill than ocean reanalyses, explaining about 50% of the variance in the adjusted annual tide gauge sea level records. Additional numerical experiments show that annual sea level changes along this coast from the barotropic model are driven by local wind stress over the continental shelf and slope. This result is interpreted in the light of a simple dynamic framework, wherein bottom friction balances surface wind stress in the alongshore direction and geostrophy holds in the across-shore direction. Results highlight the importance of barotropic dynamics on coastal sea level changes on interannual and decadal time scales; they also have implications for diagnosing the uncertainties in current ocean reanalyses, using tide gauge records to infer past changes in ocean circulation, and identifying the physical mechanisms responsible for projected future regional sea level rise.


2009 ◽  
Vol 5 (2) ◽  
pp. 1109-1132 ◽  
Author(s):  
W. Llovel ◽  
A. Cazenave ◽  
P. Rogel ◽  
A. Lombard ◽  
M. Bergé-Nguyen

Abstract. A two-dimensional reconstruction of past sea level is proposed at yearly interval over the period 1950–2003 using tide gauge records at 99 selected sites and 44-year long (1960–2003) 2°×2° gridded dynamic heights from the OPA/NEMO global ocean circulation model with data assimilation. An Empirical Orthogonal Function decomposition of the reconstructed sea level over 1950–2003 displays leading modes that reflect two main components: a long-term (multi-decadal) but regionally variable signal and interannual fluctuations dominated by the signature of El Nino-Southern Oscillation. Tests show that spatial trend patterns of the 54-year long reconstructed sea level (1950–2003) significantly depend on the length of the gridded OPA/NEMO time series used to compute spatial covariance signal used for the reconstruction (i.e., the length of the gridded OPA/NEMO time series). On the other hand, the interannual variability is well reconstructed, even with ~10-year long of the OPA/NEMO model or satellite altimetry-based sea level grids. The robustness of the results is assessed, leaving out successively each of the 99 tide gauges when reconstructing the sea level signal and then comparing observed and reconstructed time series at the non contributing tide gauge site. The reconstruction performs well at most tide gauges, especially at interannual frequency.


2014 ◽  
Vol 57 (4) ◽  
Author(s):  
Giorgio Spada ◽  
Marco Olivieri ◽  
Gaia Galassi

<p>Observations from the global array of tide gauges show that global sea-level has been rising at an average rate of 1.5-2 mm/yr during the last ~150 years [Douglas 1991, Spada and Galassi 2012]. Although a global sea-level acceleration was initially ruled out [Douglas 1992], subsequent studies [Douglas 1997, Church and White 2006, Jevrejeva et al. 2008, Church and White 2011] have coherently proposed values of ~1 mm/year/century [Olivieri and Spada 2013]. More complex non-linear trends and abrupt sea-level variations have now also been recognized. Globally, these could manifest a regime shift between the late Holocene and the current rhythms of sea-level rise [Gehrels and Woodworth 2013], while locally they result from ocean circulation anomalies, steric effects and wind stress [Bromirski et al. 2011, Merrifield 2011]. Although isostatic readjustment affects the local rates of secular sea-level change [Milne and Mitrovica 1998, Peltier 2004], a possible impact on regional acceleration has been so far discounted [Douglas 1992, Jevrejeva et al. 2008, Woodworth et al. 2009] since the process evolves on a millennium time scale [Turcotte and Schubert 2002]. Here we report a previously unnoticed anomaly in the long-term sea-level acceleration of the Baltic Sea tide gauge records, and we explain it by the classical post-glacial rebound theory and numerical modeling of glacial isostasy. Contrary to previous assumptions, our findings demonstrate that isostatic compensation plays a role in the regional secular sea-level acceleration.</p>


2017 ◽  
Vol 7 (1) ◽  
Author(s):  
H. Bâki Iz ◽  
C. K. Shum ◽  
C. Zhang ◽  
C. Y. Kuo

AbstractThis study demonstrates that relative sea level trends calculated from long-term tide gauge records can be used to estimate relative vertical crustal velocities in a region with high accuracy. A comparison of the weighted averages of the relative sea level trends estimated at six tide gauge stations in two clusters along the Eastern coast of United States, in Florida and in Maryland, reveals a statistically significant regional vertical crustal motion of Maryland with respect to Florida with a subsidence rate of −1.15±0.15 mm/yr identified predominantly due to the ongoing glacial isostatic adjustment process. The estimate is a consilience value to validate vertical crustal velocities calculated from GPS time series as well as towards constraining predictive GIA models in these regions.


2021 ◽  
Vol 11 (11) ◽  
pp. 5286
Author(s):  
Yihao Wu ◽  
Jia Huang ◽  
Hongkai Shi ◽  
Xiufeng He

Mean dynamic topography (MDT) is crucial for research in oceanography and climatology. The optimal interpolation method (OIM) is applied to MDT modeling, where the error variance–covariance information of the observations is established. The global geopotential model (GGM) derived from GOCE (Gravity Field and Steady-State Ocean Circulation Explorer) gravity data and the mean sea surface model derived from satellite altimetry data are combined to construct MDT. Numerical experiments in the Kuroshio over Japan show that the use of recently released GOCE-derived GGM derives a better MDT compared to the previous models. The MDT solution computed based on the sixth-generation model illustrates a lower level of root mean square error (77.0 mm) compared with the ocean reanalysis data, which is 2.4 mm (5.4 mm) smaller than that derived from the fifth-generation (fourth-generation) model. This illustrates that the accumulation of GOCE data and updated data preprocessing methods can be beneficial for MDT recovery. Moreover, the results show that the OIM outperforms the Gaussian filtering approach, where the geostrophic velocity derived from the OIM method has a smaller misfit against the buoy data, by a magnitude of 10 mm/s (17 mm/s) when the zonal (meridional) component is validated. This is mainly due to the error information of input data being used in the optimal interpolation method, which may obtain more reasonable weights of observations than the Gaussian filtering method.


2018 ◽  
Vol 41 (6) ◽  
pp. 517-545 ◽  
Author(s):  
Ole Baltazar Andersen ◽  
Karina Nielsen ◽  
Per Knudsen ◽  
Chris W. Hughes ◽  
Rory Bingham ◽  
...  

2021 ◽  
Author(s):  
Mika Rantanen ◽  
Jani Särkkä ◽  
Jani Räihä ◽  
Matti Kämäräinen ◽  
Kirsti Jylhä

&lt;p&gt;Extremely high sea levels on the Finnish coast are typically caused by close passages of extratropical cyclones (ETCs), which raise the sea level with their associated extreme winds and lower air pressure. For coastal infrastructure, such as nuclear power plants, it is crucial to study physically possible sea level heights associated with ETCs. Such sea levels are not straightforward to determine from observational datasets only, because tide gauge records&amp;#160; cover about 100 years and do not necessarily capture the most extreme cases having return periods longer than 100 years.&lt;/p&gt;&lt;p&gt;In this study, a method for generating an ensemble of synthetic low-pressure systems is being developed to investigate the extreme sea level heights on the Finnish coast of Baltic sea. As input parameters for the method, the point of origin, velocity of the center of the cyclone and depth of the pressure anomaly need to be given. Based on the input parameters, the method forms an idealized low-pressure system using a two-dimensional Gaussian function. In order to find extreme, but still reasonable values for the input parameters, cyclone tracks from ERA5 reanalysis data will be analysed.&lt;/p&gt;&lt;p&gt;The ensemble of synthetic low pressure systems (i.e. the wind and pressure data) is used as an input to a numerical sea level model. As a result, we have an ensemble of simulated sea levels, from which we can determine the properties of the ETCs that induce the highest sea levels on a given location on the coast. The preliminary simulation results show that this method works well, forming a basis for studies on extreme sea levels.&amp;#160;&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;


2021 ◽  
Author(s):  
Lojze Žust ◽  
Matjaž Ličer ◽  
Anja Fettich ◽  
Matej Kristan

&lt;p&gt;Interactions between atmospheric forcing, topographic constraints to air and water flow, and resonant character of the basin make sea level modeling in Adriatic a challenging problem. In this study we present an ensemble deep-neural-network-based sea level forecasting method HIDRA, which outperforms our setup of the general ocean circulation model ensemble (NEMO v3.6) for all forecast lead times and at a minuscule fraction of the numerical cost (order of 2 &amp;#215; 10&lt;sup&gt;-6&lt;/sup&gt;). HIDRA exhibits larger bias but lower RMSE than our setup of NEMO over most of the residual sea level bins. It introduces a trainable atmospheric spatial encoder and employs fusion of atmospheric and sea level features into a self-contained network which enables discriminative feature learning. HIDRA architecture building blocks are experimentally analyzed in detail and compared to alternative approaches. Results show the importance of sea level input for forecast lead times below 24 h and the importance of atmospheric input for longer lead times. The best performance is achieved by considering the input as the total sea level, split into disjoint sets of tidal and residual signals. This enables HIDRA to optimize the prediction fidelity with respect to atmospheric forcing while compensating for the errors in the tidal model. HIDRA is trained and analysed on a ten-year (2006-2016) timeseries of atmospheric surface fields from a single member of ECMWF atmospheric ensemble. In the testing phase, both HIDRA and NEMO ensemble systems are forced by the ECMWF atmospheric ensemble. Their performance is evaluated on a one-year (2019) hourly time series from tide gauge in Koper (Slovenia). Spectral and continuous wavelet analysis of the forecasts at the semi-diurnal frequency (12 h)&lt;sup&gt;-1&lt;/sup&gt; and at the ground-state basin seiche frequency (21.5 h)&lt;sup&gt;-1&lt;/sup&gt; is performed. The energy at the basin seiche in the HIDRA forecast is close to the observed, while our setup of NEMO underestimates it. Analyses of the January 2015 and November 2019 storm surges indicate that HIDRA has learned to mimic the timing and amplitude of basin seiches.&lt;/p&gt;


Sign in / Sign up

Export Citation Format

Share Document