Ocean Tides Research Articles

Anelastic response of the Earth's crust underneath the Canary Islands revealed from ocean tide loading observations

SUMMARY We report on the analysis of M2 ocean tide loading (OTL) kinematic GPS vertical displacement and tidal gravity measurements using 26 GPS and four gravimetric sites across the Canary Islands archipelago. In this region, the standard deviation among recent ocean tide models is lower than 0.4 cm in amplitude and 0.3° in phase, which are suitably accurate for displacement modelling. However, for gravity we need to model regional ocean tides to achieve enough accuracy in the loading calculations. Particularly, this study improves the predicted OTL gravity variations when global ocean models are replaced with the regional model CIAM2 which assimilates local tide gauge data. These small ocean tide model errors allow us to use the differences between observed and predicted OTL values to study the elastic and anelastic properties of the solid Earth around the Canary Islands. In the prediction of OTL, we first used the recent elastic STW105 and S362ANI seismic models, obtaining average observed minus predicted residuals of 1.2–1.3 mm for vertical displacement and 3 nm s−2 for gravity. After the STW105 and S362ANI models were adjusted for anelasticity, by considering a constant quality factor Q at periods ranging from 1 s to 12.42 hr, the average misfit between observations and predicted OTL values reduced to 0.7–0.8 mm for vertical displacement and to 1 nm s−2 for gravity. However, the average vertical displacement misfit is made up from site misfits less than 0.5 mm in western islands but for the easternmost islands of Lanzarote and Fuerteventura, they still reach up to nearly 2 mm at some sites, which still exceeds the uncertainty in the GPS observations. It is hypothesized that mantle upwelling underneath the Canary Islands, creating spatial variations in the elastic properties, causes the large residuals observed in the eastern islands. We reduced the shear modulus by up to 35 per cent in the upper mantle layer of 24.4–220 km depth. This produced residual observed minus model differences of about 0.7 mm for the sites on Lanzarote and Fuerteventura, comparable to the results obtained for the GPS sites across the rest of the archipelago, whose residuals in turn were also slightly reduced through the VS velocity and shear modulus reductions (by 0.2 mm on average).

Assessment of sub-daily ocean tide loading errors and mitigation of their propagation in multi-GNSS position time series

Assessment of sub-daily ocean tide loading errors and mitigation of their propagation in multi-GNSS position time series

Modeling impacts of saltwater intrusion on methane and nitrous oxide emissions in tidal forested wetlands.

Emissions of methane (CH4 ) and nitrous oxide (N2 O) from soils to the atmosphere can offset the benefits of carbon sequestration for climate change mitigation. While past study has suggested that both CH4 and N2 O emissions from tidal freshwater forested wetlands (TFFW) are generally low, the impacts of coastal droughts and drought-induced saltwater intrusion on CH4 and N2 O emissions remain unclear. In this study, a process-driven biogeochemistry model, Tidal Freshwater Wetland DeNitrification-DeComposition (TFW-DNDC), was applied to examine the responses of CH4 and N2 O emissions to episodic drought-induced saltwater intrusion in TFFW along the Waccamaw River and Savannah River, USA. These sites encompass landscape gradients of both surface and porewater salinity as influenced by Atlantic Ocean tides superimposed on periodic droughts. Surprisingly, CH4 and N2 O emission responsiveness to coastal droughts and drought-induced saltwater intrusion varied greatly between river systems and among local geomorphologic settings. This reflected the complexity of wetland CH4 and N2 O emissions and suggests that simple linkages to salinity may not always be relevant, as non-linear relationships dominated our simulations. Along the Savannah River, N2 O emissions in the moderate-oligohaline tidal forest site tended to increase dramatically under the drought condition, while CH4 emission decreased. For the Waccamaw River, emissions of both CH4 and N2 O in the moderate-oligohaline tidal forest site tended to decrease under the drought condition, but the capacity of the moderate-oligohaline tidal forest to serve as a carbon sink was substantially reduced due to significant declines in net primary productivity and soil organic carbon sequestration rates as salinity killed the dominant freshwater vegetation. These changes in fluxes of CH4 and N2 O reflect crucial synergistic effects of soil salinity and water level on C and N dynamics in TFFW due to drought-induced seawater intrusion.

Residual M2 and S2 ocean tide signals in complex coastal zones identified by X-Track reprocessed altimetry data

While the reliability of reprocessed altimetry data from X-Track has been examined for sea level studies, it has not previously been used to study ocean tides in tidally complex coastal zones. Using a new dataset of short-term tide gauges, the aim of this study is to demonstrate the improvements in X-Track reprocessed data for studying ocean tides for the first time in tidally complex regions in northwest Australia. We show that X-Track can be used to evaluate and potentially correct ocean tide models in challenging coastal areas. Our study comprises three components: (1) evaluate the tidal constants from X-Track and RADS datasets using nearby (<12 km) short-term tide gauge data, (2) assess three ocean tide models (FES2014, DTU10, TPXO9) using the X-Track altimetry dataset and (3) extract residual M2 and S2 tidal constants from sea level anomalies (SLAs) computed using the optimal model (FES2014) for ocean tide correction and then use these to predict the residual tidal height to use as an additional correction for the recomputation of the SLAs.Our results show that X-Track data provides users with more valid observations close to coastlines (as close as 5 km) so that the derived M2 and S2 tidal constants show better agreement with nearby tide gauges compared to RADS for both Jason-3 (J3) and GFO missions. In Van Diemen Gulf, where the footprints were 7 km or less from the coast, RADS did not have any valid observations, whereas tidal constants of M2 computed from X-Track J3 and GFO data agreed with the tide gauges within ∼6% of the M2 tidal height. In contrast, X-Track J3 and GFO extracted S2 tidal constants agreed with nearby tide gauges within ∼11–24% (GFO extreme of 47%) of the S2 tidal height. The FES2014-detided SLA time series from X-Track were then analysed to detect the remaining signals of M2 and S2 tidal constituents, with co-range maps suggesting the unmodelled amplitudes for M2 reaches 40 cm and S2 up to 30 cm in Collier Bay. In Van Diemen Gulf, M2 contributes more to the residual tidal signal than S2 with an amplitude of 40 cm compared to 10 cm respectively, in contrast to Joseph Bonaparte Gulf with ∼10 cm for both. When the SLAs, are corrected for the residual tide heights of M2 and S2, the standard deviation of the SLA residuals is decreased by as much as ∼70% in Collier Bay indicating that X-Track altimetry can identify M2 and S2 tidal residuals in complex coastal regions.

A High Precision Tidal Model Based on Multisource Tidal Harmonic Constants for Offshore China

As the accuracy of analytical instruments improve, the accuracy of ocean tide load displacement (OTLD) in geodesy and ocean dynamics is becoming increasingly demanding. The method of integrating through the load Green’s function using the ocean tide model is the main means of OTLD correction. However, although different tidal models have high accuracy in open ocean regions, they are less precise and less practical in the offshore areas of China, making it difficult to correct OTLDs accurately in complicated coastline areas of China. In this work, we evaluate the accuracy of 13 different tidal models in offshore China, and construct a new tidal model, new.chinasea.2022, based on the combination of high-precision tidal harmonic parameters. The OTLD correction efficiency of new.chinasea.2022 and conventional tidal models was assessed by wRMS analysis of GPS coordinate time series from the Crustal Movement Observation Network of China (CMONOC). It turns out that (1) the newly released tidal model performs better than the old model in offshore China, and there are differences in the performance of different models for each tidal constituent in different seas; (2) The OTLD correction of new.chinasea.2022 performed better in the offshore region of China, improving the correction of FES2004 by 1%-12% and osu.chinasea.2010 by 1%-6.5%.

Evolution of Global Ocean Tide Levels Since the Last Glacial Maximum

AbstractThis study addresses the evolution of global tidal dynamics since the Last Glacial Maximum focusing on the extraction of tidal levels that are vital for the interpretation of geologic sea‐level markers. For this purpose, we employ a truly‐global barotropic ocean tide model which considers the non‐local effect of Self‐Attraction and Loading. A comparison to a global tide gauge data set for modern conditions yields agreement levels of 65%–70%. As the chosen model is data‐unconstrained, and the considered dissipation mechanisms are well understood, it does not have to be re‐tuned for altered paleoceanographic conditions. In agreement with prior studies, we find that changes in bathymetry during glaciation and deglaciation do exert critical control over the modeling results with minor impact by ocean stratification and sea ice friction. Simulations of 4 major partial tides are repeated in time steps of 0.5–1 ka and augmented by 4 additional partial tides estimated via linear admittance. These are then used to derive time series from which the tidal levels are determined and provided as a global data set conforming to the HOLSEA format. The modeling results indicate a strengthened tidal resonance by M2, but also by O1, under glacial conditions, in accordance with prior studies. Especially, a number of prominent changes in local resonance conditions are identified, that impact the tidal levels up to several meters difference. Among other regions, resonant features are predicted for the North Atlantic, the South China Sea, and the Arctic Ocean.

Tidal Evolution and Predictable Tide‐Only Inundation Along the East Coast of the United States

AbstractSea level rise (SLR) increases the frequency of floods caused by forces such as spring tides or storm surges. However, with the continuous evolution of global ocean tides, future flooding may be driven only by astronomical tides. Therefore, we introduce an inaction method (IM) that can accurately extract the harmonic constant of tidal components from the time‐frequency spectrum of the tide, and obtain the precise long‐term trend and spatial evolution characteristics of the tides at 18 sites along the East Coast of the United States (ECUS). Except for some sites in the middle of the East Coast, M2 shows an average upward trend of more than 51 mm/century, which is most obvious at three sites in North Carolina, with an average upward trend of more than 80 mm/century. We simultaneously combine the IM with harmonic analysis into an IMHA method to estimate the historical and future “Tide‐only” inundation (TOI) frequency and quantify the effect of tidal evolution on the future TOI frequency. The results show that four stations in the Gulf of Maine have experienced relatively frequent TOI events in the past, with the highest frequency at the Eastport station. The most significant increasing and decreasing effects of tidal evolution on the future TOI frequency under the SLR scenario are at the Wilmington and Lewes stations, respectively. At the spatial scale, there are increases in the northeast and southeast and a decrease in the middle of the ECUS, which is consistent with the spatial characteristics of the long‐term tidal trend.

Sea Tide Influence on Ice Flow of David Drygalski’s Ice Tongue Inferred from Geodetic GNSS Observations and SAR Offset Tracking Analysis

David Glacier and Drygalski Ice Tongue are massive glaciers in Victoria Land, Antarctica. The ice from the East Antarctic Ice Sheet is drained through the former, and then discharged into the western Ross Sea through the latter. David Drygalski is the largest outlet glacier in Northern Victoria Land, floating kilometers out to sea. The floating and grounded part of the David Glacier are the main focus of this article. During the XXI Italian Antarctic Expedition (2005–2006), within the framework of the National Antarctic Research Programme (PNRA), two GNSS stations were installed at different points: the first close to the grounding line of David Glacier, and the second approximately 40 km downstream of the first one. Simultaneous data logging was performed by both GNSS stations for 24 days. In the latest data processing, the kinematic PPP technique was adopted to evaluate the dominant diurnal components and the very small semi-diurnal variations in ice motion induced by the ocean tide and the mean ice flow rates of both GNSS stations. Comparison of the GNSS time series with predicted ocean tide calculated from harmonic coefficients of the nearest tide gauge stations, installed at Cape Roberts and Mario Zucchelli Station, highlight different local response of the glacier to ocean tide, with a minor amplitude of vertical motion at a point partially anchored at the bedrock close to the grounding line. During low tide, the velocity of the ice flow reaches its daily maximum, in accordance with the direction of seawater outflow from the fjord into the ocean, while the greatest daily tidal excursion generates an increase in the horizontal ice flow velocity. With the aim to extend the analysis in spatial terms, five COSMO-SkyMED Stripmap scenes were processed. The comparison of the co-registered offset tracking rates, obtained from SAR images, with the GNSS estimation shows good agreement.

Rapid solidification of Earth’s magma ocean limits early lunar recession

The early evolution of the Earth–Moon system prescribes the tidal environment of the Hadean Earth and holds the key to the formation mechanism of the Moon and its thermal evolution. Estimating its early state by backtracking from the present, however, suffers from substantial uncertainties associated with ocean tides. Tidal evolution during the solidification of Earth’s magma ocean, on the other hand, has the potential to provide robust constraints on the Earth–Moon system before the appearance of a water ocean. Here we show that energy dissipation in a solidifying magma ocean results in considerably more limited lunar recession than previously thought, and that the Moon was probably still at the distance of ∼7–9 Earth radii at the end of solidification. This limited early recession aggravates the often overlooked difficulty of modeling tidal dissipation in Earth’s first billion years, but it also offers a new possibility of resolving the lunar inclination problem by allowing the operation of multiple excitation mechanisms.

Tidal analysis and prediction based on the Fourier basis pursuit spectrum

The high-precision prediction of ocean tides is important for coastal management. Traditional harmonic analysis (HA) regards tides as the cumulative motion of multiple trigonometric functions, ignoring the time-varying characteristics of harmonic constants. Although the time-frequency transformation method can reconstruct the time-varying characteristics of tides, there are problems with tidal component coupling and edge effects. Therefore, we propose a method of reconstructing tidal components based on Fourier basis pursuit (FBP) spectral bandpass filtering (FBPBPF) by combining the Fourier dictionary and the basis pursuit algorithm and apply this method to four tide stations in the Gulf of Maine and Hong Kong station to try to solve the above two problems. The results show that the FBP spectrum can effectively identify and separate the main coupling tidal components. In the Gulf of Maine, FBPBPF can use shorter tide data to obtain better or equivalent prediction results than HA. Especially within 30 days, the prediction accuracy of FBPBPF at the four sites can be improved by 53.5% on average compared with that of HA. At Hong Kong station, FBPBPF has higher prediction accuracy than HA and Inaction Method (IM) in 4 months, but its long-term accuracy is gradually reduced (including all stations). We believe that the frequency resolution is the key factor leading to the divergence of long-term prediction errors.

A superconducting gravimeter on the island of Heligoland for the high-accuracy determination of regional ocean tide loading signals of the North Sea

SUMMARY The superconducting gravimeter GWR iGrav 047 has been installed on the small offshore island of Heligoland in the North Sea approximately at sea level with the overall aim of high-accuracy determination of regional tidal and non-tidal ocean loading signals. For validation, a second gravimeter (gPhoneX 152) has been setup within a gravity gradiometer approach to observe temporal gravity variations in parallel on the upper land of Heligoland. This study covers the determination of regional ocean tide loading (OTL) parameters based on the two continuous gravimetric time-series after elimination of the height-dependent gravity component by empirical transfer functions between the local sea level from a nearby tide gauge and local attraction effects. After reduction of all gravity recordings to sea level, both gravimeters provide very similar height-independent OTL parameters for the eight major diurnal and semidiurnal waves with estimated amplitudes between 0.3 nm s−2 (Q1) and 11 nm s−2 (M2) and RMSE of 0.1–0.2 nm s−2 for 2 yr of iGrav 047 observations and a factor of 2 worse for 1.5 yr of gPhoneX 152 observations. The mean absolute OTL amplitude differences are 0.3 nm s−2 between iGrav 047 and gPhoneX 152, 0.4 nm s−2 between iGrav 047 and the ocean tide model FES2014b and 0.7 nm s−2 between gPhoneX 152 and FES2014b which is in good agreement with the uncertainty estimations. As by-product of this study, OTL vertical displacements are estimated from the height-independent OTL gravity results from iGrav 047 applying proportionality factors ${\rm d}h/{\rm d}g$ for the eight major waves. These height-to-gravity ratios and the corresponding phase shifts are derived from FES2014b. The OTL vertical displacements from iGrav 047 are estimated with amplitudes between 0.4 mm (Q1) and 5.1 mm (M2) and RMSE of 0.1–0.7 mm. These OTL amplitudes agree with FES2014b within 0.0 (M2) and 0.8 mm (K1) with a mean difference of 0.3 mm only. The OTL amplitudes from almost 5 yr of GNSS observations show deviations of up to 6 mm (M2) compared to vertical displacements from both iGrav 047 and FES2014b, which suggests systematic effects included in the estimation of OTL vertical displacements from GNSS. With the demonstrated accuracy, height-independent sensitivity in terms of gravity and vertical displacements along with the high temporal resolution and the even better performance with length of time-series, iGrav 047 delivers the best observational signal for OTL which is representative for a large part of the North Sea.

Ocean Tides near Hawaii from Satellite Altimeter Data. Part III

Abstract The Chebyshev polynomial fitting (CPF) method has been proved to be effective to construct reliable cotidal charts for the eight major tidal constituents (M2, S2, K1, O1, N2, K2, P1, and Q1) and six minor tidal constituents (2N2, J1, L2, Mu2, Nu2, and T2) near Hawaii in Part I and Part II, respectively. In this paper, this method is extended to estimate the harmonic constants of four long-period tidal constituents (Mf, Mm, Sa, and Ssa). The harmonic constants obtained by this method were compared with those from the TPXO9, Finite Element Solutions 2014 (FES2014), and Empirical Ocean Tide 20 (EOT20) models, using benchmark data from satellite altimeters and eight tide gauges. The accuracies of the Mf and Mm constituents derived from the CPF method are comparable to those from the models, but the accuracies of the Sa and Ssa constituents are significantly higher than those from the FES2014 and EOT20 models. The results indicate that the CPF method is also effective for estimating harmonic constants of long-period tidal constituents. Furthermore, since the CPF method relies only on satellite altimeter data, it is an easier-to-use method than these ocean tide models.

Terdiurnal Radiational Tides

Abstract Terdiurnal atmospheric tides induce an S3 radiational ocean tide, similar to radiational tides S1 and S2 in the diurnal and semidiurnal bands. Although of small amplitude, the terdiurnal tide has some intriguing properties. The tide has an unusually pronounced seasonal variation, manifested by annual sidelines here denoted R3 and T3, which causes the tide to nearly vanish during times near an equinox. Forcing is generally largest in the winter hemisphere. Complicating matters, the two sideline frequencies coincide with those of nonlinear compound tides SK3 and SP3. Whether radiational tides or nonlinear tides (or both) are appearing at any given tide gauge can usually be determined by the relative amplitudes and phase differences of the two sidelines. The amplitudes of R3 and T3 are generally comparable; the amplitudes of SK3 and SP3 are not. Proper identification can lead to a small improvement in tidal prediction, but more importantly can lead to improved physical interpretation. An example from recent measurements under the Ross Ice Shelf bears on the role of nonlinearity in interactions between the ocean tide and the floating ice shelf.

High‐Tide Flooding Along the China Coastline: Past and Future

AbstractCoastal flooding has become a major issue for low‐lying coastal cities in China, and a lot of research has focused on assessing flood risk from storms and associated extreme sea levels. High‐tide flooding (HTF), however, which leads to minor inundation and occurs more frequently as sea level continues to rise has not been assessed comprehensively. Here, we analyze HTF along the China coastline using tide gauge records. We show that the frequency of HTF has increased (Xiamen doubled HTF frequency every 11.5–37.4 years, with the median being 17.5 years). As a result, the cumulative loss ratios of HTF are higher than those of major (or extreme) flooding events in several locations. To gain insights into the processes driving HTF changes, we decompose still water levels that occurred during HTF into five components, including non‐linear trend, interannual to decadal mean sea level (MSL) variability, seasonal MSL cycle, tidal anomaly, and nontidal residuals. It is evident that due to sea‐level rise (SLR) fewer components need to combine to raise the water levels above HTF thresholds. We show that the South China Sea coast already experiences HTF purely driven by high spring tides, and will also see the fastest future increase in the number of tide‐only HTF events. In general, China will experience more HTF days under all warming scenarios as SLRs. This long‐term trend will be modulated by the nodal cycle of ocean tides leading to more rapid increases in HTF in the 2030s and 2050s.

Estimation of the Manning’s n coefficient in multi-constituent tidal models by assimilating satellite observations with the adjoint data assimilation

The bottom friction is critical for the dissipation of the global tidal energy. The bottom friction coefficient is traditionally determined using the Manning’s n formulation in tidal models. The Manning’s n coefficient in the Manning’s n formulation is vital for the accurate simulation and prediction of the tide in coastal shallow waters, but it cannot be directly measured and contains large amounts of uncertainties. Based on a two-dimensional multi-constituent tidal model with the adjoint data assimilation, the estimation of the Manning’s n coefficient is investigated by assimilating satellite observations in the Bohai, Yellow and East China Seas with the simulation of four principal tidal constituents M2, S2, K1 and O1. In the twin experiments, the Manning’s n coefficient is assumed to be constant, and it is estimated by assimilating the synthetic observations at the spatial locations of the satellite tracks. Regardless the inclusion of artificial random observational errors associated with synthetic observations, the model performance is improved as evaluated by the independent synthetic observations. The prescribed ‘real’ Manning’s n coefficient is reasonably estimated, indicating that the adjoint data assimilation is an effective method to estimate the Manning’s n coefficient in multi-constituent tidal models. In the practical experiments, the errors between the independent observations at the tidal gauge stations and the corresponding simulated results of the four principal tidal constituents are substantially decreased under both scenarios of the constant and spatially-temporally varying Manning’s n coefficient estimated by assimilating the satellite observations with the adjoint data assimilation. In addition, the estimated spatial and temporal variation trend is robust and not affected by the model settings. The spatially-temporally varying Manning’s n coefficient is negatively correlated with the current speed and shows significant spatial variation in the shallow water areas. This study demonstrates that the Manning’s n coefficient can be reasonably estimated by the adjoint data assimilation, which allows significant improvement in accurate simulation of the ocean tide.

Occurrence and sources of micro-plastics in various water bodies, sediments, and fishes in Ansan, South Korea

In this study, the Pearson correlation coefficients were determined to derive correlations between micro-plastics (MPs) in carp and river crabs. MPs were detected for various water sources, including four rivers and four main waterways, sediments, and fish, using Fourier transform infrared spectrometry (FTIR), microscopic analysis, and image mapping. Carp and river crabs had coefficients of 0.888 and 0.724, respectively, which showed a high positive correlation. In water samples, the MPs detected in rivers were higher than those in the main waterway. However, in sediment samples, the MPs detected in the main waterway were higher than those in the rivers. It is believed that MPs are carried toward shore by ocean tide. The size of most of the sediment MPs was 20-49µm, representing 64.1% of the entire population. The plastics detected in this study were polyethylene terephthalate (PET), polypropylene (PP), and polyethylene (PE), which originate from synthetic fibers, scrubs, and packing material. MP pollution by non-point pollution sources was investigated, with the abundance of MPs increasing by 2 to 3 times between the dry and wet seasons in water and sediment, respectively. It was determined that the inflow of MPs into rivers could have been due to non-point source pollutants from household items, roads, plants, and soil around the water sources.

Grounding Zone of Amery Ice Shelf, Antarctica, From Differential Synthetic‐Aperture Radar Interferometry

AbstractWe employ a time series of Sentinel‐1 differential radar interferometry data from 2018 to detect the variability in grounding line position of the Fisher, Mellor, and Lambert glaciers, which drain about 47 billion tons of ice per year from East Antarctica. We observe kilometer‐scale tidal migration, two orders of magnitude larger than expected for ice flowing over a hard bed. The migration is not in phase with changes in oceanic tide. In two estuaries underlaid by subglacial channels, we observe two states of migration that switch on and off over time scales of several weeks. The range of vertical motion reveals a water column thickness of 2–20 cm. Such intrusions of seawater over wide grounding zones are not accounted for in physical models. Including them will add vigorous melting of grounded ice that will enhance the sensitivity of glaciers to ocean warming and increase projections of mass loss.

Assessment of a near-polar pair mission for detecting the Earth's temporal gravity field

SUMMARYTo ensure the global observability of next generation gravimetric mission (NGGM), different agencies have to repeatedly launch satellites to about 89.0° orbit inclination. However, due to the poor isotropy of observation system, only minor improvement in terms of temporal gravity field estimation can be obtained via these repeatedly launched polar pair missions. To ensure the global observability as well as the isotropy of observation system, a near-polar pair mission rather than a polar mission is likely an optimal selection, especially considering that the Gravity Recovery and Climate Experiment (GRACE) Follow-On mission has been already in operation. In this study, for the upcoming NGGMs (for instance, the Chinese NGGM), we design a closed-loop simulation to assess the performance of a near-polar mission at a near-circular orbit with about 500 km altitude for detecting the Earth's temporal gravity field, and the main conclusions are summarized as follows. (1) Based on the statistic results, 85.0° is selected as the optimal orbit inclination for the near-polar mission, which provides 37 per cent noise reduction in terms of cumulative geoid height error in spectral domain, but also 31 per cent noise reduction in terms of mean oceanic root-mean-square (RMS) error in spatial domain when compared to the 89.0° polar mission (89-PM). (2) To figure out the reason of these noise reductions, we also compare the contribution of single error components [including instrument error, atmospheric and oceanic (AO) aliasing error, ocean tide error]. It indicates that the outperformance of the 85.0° near-polar mission (85-NPM) is mainly derived from the relatively smaller non-tidal atmospheric and oceanic mass variation error. (3) Although there are inevitable 5.0° polar gaps in the 85-NPM, the analysis result (including the comparison via in-orbit observations and simulated retrieved solutions) still demonstrates a comparable performance of the 85-NPM in tracking mass variations over the Antarctic, and even an outperformance with 12 per cent noise reduction over the Greenland when compared to the 89-PM. The result confirms the feasibility of implementing a near-polar mission as a stand-alone mission or a complementary observation system for the repeatedly launched polar missions, which offers an alternative option of launching the Chinese NGGM satellites to an 85.0° inclination orbit instead of 89.0°.

Barotropic tides in MPAS-Ocean (E3SM V2): impact of ice shelf cavities

Abstract. Oceanic tides are seldom represented in Earth system models (ESMs) owing to the need for high horizontal resolution to accurately represent the associated barotropic waves close to coasts. This paper presents results of tides implemented in the Model for Prediction Across Scales–Ocean or MPAS-Ocean, which is the ocean component within the U.S. Department of Energy developed Energy Exascale Earth System Model (E3SM). MPAS-Ocean circumvents the limitation of low resolution using unstructured global meshing. We are at this stage simulating the largest semidiurnal (M2, S2, N2) and diurnal (K1, O1) tidal constituents in a single-layer version of MPAS-O. First, we show that the tidal constituents calculated using MPAS-Ocean closely agree with the results of the global tidal prediction model TPXO8 when suitably tuned topographic wave drag and bottom drag coefficients are employed. Thereafter, we present the sensitivity of global tidal evolution due to the presence of Antarctic ice shelf cavities. The effect of ice shelves on the amplitude and phase of tidal constituents are presented. Lower values of complex errors (with respect to TPXO8 results) for the M2 tidal constituents are observed when the ice shelf is added in the simulations, with particularly strong improvement in the Southern Ocean. Our work points towards future research with varying Antarctic ice shelf geometries and sea ice coupling that might lead to better comparison and prediction of tides and thus better prediction of sea-level rise and also the future climate variability.

Multidecadal and climatological surface current simulations for the southwestern Indian Ocean at 1∕50° resolution

Abstract. The Western INDian Ocean Simulation (WINDS) is a regional configuration of the Coastal and Regional Ocean Community Model (CROCO) for the southwestern Indian Ocean. WINDS has a horizontal resolution of 1/50∘ (∼2 km) and spans a latitudinal range of 23.5∘ S–0∘ N and a longitudinal range from the East African coast to 77.5∘ E. We ran two experiments using the WINDS configuration: WINDS-M, a full 28-year multidecadal run (1993–2020); and WINDS-C, a 10-year climatological control run with monthly climatological forcing. WINDS was primarily run for buoyant Lagrangian particle tracking applications, and horizontal surface velocities are output at a temporal resolution of 30 min. Other surface fields are output daily, and the full 3D temperature, salinity, and velocity fields are output every 5 d. We demonstrate that WINDS successfully manages to reproduce surface temperature, salinity, currents, and tides in the southwestern Indian Ocean, and it is therefore appropriate for use in regional marine dispersal studies for buoyant particles or other applications using high-resolution surface ocean properties.