Representation Of Bathymetry Research Articles

A new high-resolution Coastal Ice-Ocean Prediction System for the East Coast of Canada

The Coastal Ice Ocean Prediction System for the East Coast of Canada (CIOPS-E) was developed and implemented operationally at Environment and Climate Change Canada (ECCC) to support a variety of critical marine applications. These include support for ice services, search and rescue, environmental emergency response and maritime safety. CIOPS-E uses a 1/36° horizontal grid (~ 2 km) to simulate sea ice and ocean conditions over the northwest Atlantic Ocean and the Gulf of St. Lawrence (GSL). Forcing at lateral open boundaries is taken from ECCC’s data assimilative Regional Ice-Ocean Prediction System (RIOPS). A spectral nudging method is applied offshore to keep mesoscale features consistent with RIOPS. Over the continental shelf and GSL, the CIOPS-E solution is free to evolve according to the model dynamics. Overall, CIOPS-E significantly improves the representation of tidal and sub-tidal water levels compared to ECCC’s lower resolution systems: RIOPS (~ 6 km) and the Regional Marine Prediction System – GSL (RMPS-GSL, 5 km). Improvements in the GSL are due to the higher resolution and a better representation of bathymetry, boundary forcing and dynamics in the upper St. Lawrence Estuary. Sea surface temperatures show persistent summertime cold bias, larger in CIOPS-E than in RIOPS, as the latter is constrained by observations. The seasonal cycle of sea ice extent and volume, unconstrained in CIOPS-E, compares well with observational estimates, RIOPS and RMPS-GSL. A greater number of fine-scale features are found in CIOPS-E with narrow leads and more intense ice convergence zones, compared to both RIOPS and RMPS-GSL.

Impact of bathymetry on Indian Ocean circulation in a nested regional ocean model

The Regional Indian Ocean model based on Modular Ocean Model (MOM4p1) was used to understand the importance of a realistic representation of bathymetry on Ocean General Circulation. The model has 1/4° uniform horizontal resolution and is forced with Coordinated Ocean-Ice Reference Experiments (CORE-II) inter-annual forcing with two simulations named BLND (realistic bathymetry) and OM3 (smoothed bathymetry), which only differ in the representation of bathymetry for the years 1992–2005. We also used recent reanalysis products from ORAS5 and SODA3 and ADCP observation to compare the subsurface currents. We show that by the inclusion of realistic bathymetry, there is a significant improvement in the upper ocean salinity, temperature, and currents, particularly near the coast. The salinity and temperature of the upper ocean are very close to the observed value near the coast. The bias in the salinity and temperature was reduced to half in BLND simulation compared to OM3, which led to a more realistic East India Coastal Current (EICC). We show the first evidence of a basin-wide cyclonic gyre over the Bay of Bengal at 1000 m depth during spring, which is just opposite to that of a basin-wide anti-cyclonic gyre at the surface. We found the presence of poleward EICC during spring at 1000 m and 2000 m depth, which is opposite to that of the surface. The presence of this deeper EICC structure is completely absent during fall. We show the presence of a boundary current along the coast of Andaman and Nicobar Island at a depth of 2000 m. The observed Wyrtki Jet (WJ) magnitude and spatial structure are most realistically reproduced in BLND simulation as compared to OM3 simulations. Both ORAS5 and SODA reanalysis products underestimate the WJ magnitude. The presence of the Maldives Islands is responsible for the westward extent of Equatorial Under Current (EUC). The presence of Maldives also creates wakes on the leeward side in the EUC zonal current. During fall, EUC is better defined in the eastern Equatorial Indian Ocean and lies at a depth of between 50 and 100 m, unlike its spring counterpart, in which its core is located slightly deeper, between 100 and 150 m depth. During peak summer months, June–July, a strong eastward zonal jet is present at 1000 m depth, similar to Wyrtki Jet (WJ). Inter-monsoon Jets, i.e., spring and fall jets, are also seen but are in the opposite direction, i.e., westward, unlike eastward in WJ.

Improved Gulf Stream separation through Brinkman penalization

The advantage of a smooth representation of bathymetry in terrain-following σ-coordinate ocean models is compromised by the need to avoid numerical errors on steep slopes associated with pressure gradient discretization or spurious diapycnal diffusion. Geopotential z-coordinate models avoid these errors, but greatly underrepresent the interaction of flow with a topographic slope, especially when the bathymetry is underresolved. Hybrid coordinate models are also deficient because it is difficult to find a satisfactory compromise between z and σ coordinates. More general vertical coordinates (not just combinations of z and σ) can also be used, in particular for ocean interior, but without solving the problems associated with the representation of bathymetry. With volume penalization, we do not seek a compromise, but rather a correction to the usual coordinate systems that realistically recovers continuous and steep bathymetry. The Brinkman volume penalization method studied here is a modified version of the one introduced in Debreu et al.(2020) that simplifies the numerical implementation of the penalization, increases robustness and improves its computational performance for realistic long-term simulations, while preserving accuracy. We apply this penalization method to the Gulf Stream separation problem that has puzzled modelers for decades. The method improves the representation of the flow-topography interaction and achieves realistic separation of the Gulf Stream at resolutions as coarse as 1/8°. In addition, it provides a tool to separate the effect of eddy activity and topographic slope when changing grid resolution. This has never before been possible because at coarse resolution none of the usual coordinate systems can properly represent a steep continental slope. Our results show that realistic bathymetry is more important than eddy activity in ensuring realistic Gulf Stream separation, even though many recent studies tend to focus on the eddy activity. A steep slope can exert a stabilizing influence that promotes a strong mean slope current with strong inertia that helps it separate from the coast at the topographic curvature of Cape Hatteras. We anticipate that a successful topographic slope correction will be very valuable for climate models, as their current resolution is far from sufficient to represent western boundary currents (WBCs) using traditional coordinate systems. Our results suggest that a climate model with a 1/4° resolution using volume penalization — and perhaps also some parameterization of the eddy-mean flow interaction to energize the WBCs — could represent ocean circulation much more realistically than a model at the same resolution, but without volume penalization.

Systematic Comparison of Tsunami Simulations on the Chilean Coast Based on Different Numerical Approaches

Tsunami inundation estimates are of crucial importance to hazard and risk assessments. In the context of tsunami forecast, numerical simulations are becoming more feasible with the growth of computational power. Uncertainties regarding source determination within the first minutes after a tsunami generation might be a major concern in the issuing of an appropriate warning on the coast. However, it is also crucial to investigate differences emerging from the chosen algorithms for the tsunami simulations due to a dependency of the outcomes on the suitable model settings. In this study, we compare the tsunami inundation in three cities in central Chile (Coquimbo, Viña del Mar, and Valparaíso) using three different models (TsunAWI, Tsunami-HySEA, COMCOT) while varying the parameters such as bottom friction. TsunAWI operates on triangular meshes with variable resolution, whereas the other two codes use nested grids for the coastal area. As initial conditions of the experiments, three seismic sources (2010 Mw 8.8 Maule, 2015 Mw 8.3 Coquimbo, and 1730 Mw 9.1 Valparaíso) are considered for the experiments. Inundation areas are determined with high-resolution topo-bathymetric datasets based on specific wetting and drying implementations of the numerical models. We compare each model’s results and sensitivities with respect to parameters such as bottom friction and bathymetry representation in the varying mesh geometries. The outcomes show consistent estimates for the nearshore wave amplitude of the leading wave crest based on identical seismic source models within the codes. However, with respect to inundation, we show high sensitivity to Manning values where a non-linear behaviour is difficult to predict. Differences between the relative decrease in inundation areas and the Manning n-range (0.015–0.060) are high (11–65%), with a strong dependency on the characterization of the local topo-bathymery in the Coquimbo and Valparaíso areas. Since simulations carried out with such models are used to generate hazard estimates and warning products in an early tsunami warning context, it is crucial to investigate differences that emerge from the chosen algorithms for the tsunami simulations.

Regional scale hydrodynamic modeling of the river-floodplain-reservoir continuum

River floodplains and reservoirs interact throughout a basin drainage network, defining a coupled human-water system with multiple feedbacks. Recent modeling developments have aimed to improve the representation of such processes at regional to continental scales. However, most large-scale hydrological models adopt simplified lumped reservoir schemes, where an offline routine is run with inflows estimated by the model, with limited consideration of the complementarity between floodplains and reservoirs on altering the hydrological regime at regional scale. This paper presents a novel approach that fully couples river-floodplain-reservoir hydrodynamic and hydrological models, significantly improving the representation of reservoir dynamics and operation in the river-floodplain-reservoir continuum at large scale and across multiple dam cascades. The model is applied to the Paraná River Basin with explicit simulation of 31 large dams and river hydraulic variables at basin scale. Three types of reservoir bathymetry representation are compared, from lumped to distributed methods, combined with three reservoir operation schemes and varying degrees of input data requirement within two parameterization scenarios (global and regional setups). The operation schemes were more relevant than the reservoir bathymetry representation to estimate downstream flows and water levels. While the data-driven operation scheme, based on linear regressions between observed water levels and dam outflows, provided the best estimates of both active storage and discharges, the more generic operation reasonably estimated discharges and peak attenuation, albeit not as accurately for active storage. The global parameterization of reservoir operation resulted in poorer performance compared to the regional-based one, but it satisfactorily modeled discharge and peak attenuation. Regarding the reservoir bathymetry representation, a basin scale comparison of the lumped and distributed schemes indicated the inability of the former to represent backwater effects. This was further corroborated by validating the longitudinal water level profile of Itaipu dam with ICESat satellite altimetry data. Finally, the model was used to show the complementarity between floodplains and reservoirs on attenuating floods at regional scale. Large scale models should move beyond offline coupling strategies, and include regional-based, data-driven reservoir operation schemes together with a distributed representation of reservoir bathymetry into river-floodplain hydraulic schemes. This will largely improve the estimation of river discharges, water levels and flood storage, and thus the model ability to represent the regional scale river-floodplain-reservoir continuum.

Computational Models of Trajectory Investigation of Marine Geophysical Fields and Its Implementation for Solving Problems of Map-Aided Navigation

Comprehensive investigations of geophysical fields (GPF) of the ocean are among the top priority problems in underwater robotics. Problems reside in developing automated systems capable of real-time operation for gathering, accumulating, and processing diverse geophysical information. The data’s total volume is applied for long-term or real-time monitoring of marine areas, environment or objects surveillance, mapping certain regions, or anomalies in geographical coordinates. One of the main elements of such systems is the computational models sensitive to properties of geophysical fields, features of search routes of underwater vehicles, and particular aspects of missions performed during trajectory measurements and features of navigational support tasks. The work considers computational models of comprehensive interpretation of the results of trajectory measurements of geophysical fields using an autonomous underwater vehicle (AUV), and estimation of accuracy of map-aided navigation on the reconstructed map of geophysical fields. Algorithms and software consider distinctive features of representation of bathymetry, magnetometry, and gravimetry data visualized in 2D and 3D. Algorithm of map-aided navigation by the reconstructed fields is discussed. The final assessments of the considered models’ accuracy consider averaged errors of measurements, mapping, and inertial navigation. Obtained assessments are based on theoretical investigations, results of model experiments, and experimental trials of AUV’s systems under actual operating conditions.

A Numerical Study of Hypoxia in Chesapeake Bay Using an Unstructured Grid Model: Validation and Sensitivity to Bathymetry Representation

AbstractA three‐dimensional unstructured‐grid hydrodynamic and water quality model (Semi‐implicit Cross‐scale Hydroscience Intergrated System Model‐Integrated Compartment Model) is applied successfully for Chesapeake Bay. The model is validated with observations of salinity, chlorophyll‐a, dissolved oxygen, nutrients, and phytoplankton productions from the year 1991 to 1995 for the mainstem and some major tributaries, based on multiple model skill scores. Model experiments are conducted to test the importance of having (1) an accurate representation of bathymetry to correctly predict hypoxia and other processes and (2) a high‐resolution model grid for tributaries to correctly simulate water quality variables. Comparison with the model experiment results with bathymetry smoothing indicates that bathymetry smoothing, as commonly used for many systems, changes the stratification and lateral circulation pattern, resulting in more salt intrusion into shallow water regions, and an increase in the freshwater age. Consequently, a model with bathymetry smoothing can lead to an unrealistic prediction of the distribution of hypoxia and phytoplankton production. Local grid refinement shows significant improvement of model simulations on local stratification and water quality variables. Overall, the use of high‐resolution unstructured grid model leads to a faithful representation of the complex geometry, and thus a seamless cross‐scale capability for simulating water quality processes in the Bay including tributaries and tidal creeks.

The significance of coastal bathymetry representation for modelling the tidal response to mean sea level rise in the German Bight

Abstract. Due to climate change an accelerated mean sea level rise is expected. One key question for the development of adaptation measures is how mean sea level rise affects tidal dynamics in shelf seas such as the North Sea. Owing to its low-lying coastal areas, the German Bight (located in the southeast of the North Sea) will be especially affected. Numerical hydrodynamic models help to understand how mean sea level rise changes tidal dynamics. Models cannot adequately represent all processes in overall detail. One limiting factor is the resolution of the model grid. In this study we investigate which role the representation of the coastal bathymetry plays when analysing the response of tidal dynamics to mean sea level rise. Using a shelf model including the whole North Sea and a high-resolution hydrodynamic model of the German Bight we investigate the changes in M2 amplitude due to a mean sea level rise of 0.8 and 10 m. The shelf model and the German Bight Model react in different ways. In the simulations with a mean sea level rise of 0.8 m the M2 amplitude in the shelf model generally increases in the region of the German Bight. In contrast, the M2 amplitude in the German Bight Model increases only in some coastal areas and decreases in the northern part of the German Bight. In the simulations with a mean sea level rise of 10 m the M2 amplitude increases in both models with largely similar spatial patterns. In two case studies we adjust the German Bight Model in order to more closely resemble the shelf model. We find that a different resolution of the bathymetry results in different energy dissipation changes in response to mean sea level rise. Our results show that the resolution of the bathymetry especially in flat intertidal areas plays a crucial role for modelling the impact of mean sea level rise.

Ocean phosphorus inventory: large uncertainties in future projections on millennial timescales and their consequences for ocean deoxygenation

Abstract. Previous studies have suggested that enhanced weathering and benthic phosphorus (P) fluxes, triggered by climate warming, can increase the oceanic P inventory on millennial timescales, promoting ocean productivity and deoxygenation. In this study, we assessed the major uncertainties in projected P inventories and their imprint on ocean deoxygenation using an Earth system model of intermediate complexity for the same business-as-usual carbon dioxide (CO2) emission scenario until the year 2300 and subsequent linear decline to zero emissions until the year 3000. Our set of model experiments under the same climate scenarios but differing in their biogeochemical P parameterizations suggest a large spread in the simulated oceanic P inventory due to uncertainties in (1) assumptions for weathering parameters, (2) the representation of bathymetry on slopes and shelves in the model bathymetry, (3) the parametrization of benthic P fluxes and (4) the representation of sediment P inventories. Considering the weathering parameters closest to the present day, a limited P reservoir and prescribed anthropogenic P fluxes, we find a +30 % increase in the total global ocean P inventory by the year 5000 relative to pre-industrial levels, caused by global warming. Weathering, benthic and anthropogenic fluxes of P contributed +25 %, +3 % and +2 %, respectively. The total range of oceanic P inventory changes across all model simulations varied between +2 % and +60 %. Suboxic volumes were up to 5 times larger than in a model simulation with a constant oceanic P inventory. Considerably large amounts of the additional P left the ocean surface unused by phytoplankton via physical transport processes as preformed P. In the model, nitrogen fixation was not able to adjust the oceanic nitrogen inventory to the increasing P levels or to compensate for the nitrogen loss due to increased denitrification. This is because low temperatures and iron limitation inhibited the uptake of the extra P and growth by nitrogen fixers in polar and lower-latitude regions. We suggest that uncertainties in P weathering, nitrogen fixation and benthic P feedbacks need to be reduced to achieve more reliable projections of oceanic deoxygenation on millennial timescales.

Robot-Assisted Measurement for Hydrologic Understanding in Data Sparse Regions

This article describes the field application of small, low-cost robots for remote surface data collection and an automated workflow to support water balance computations and hydrologic understanding where water availability data is sparse. Current elevation measurement approaches, such as manual surveying and LiDAR, are costly and infrequent, leading to potential inefficiencies for quantifying the dynamic hydrologic storage capacity of the land surface over large areas. Experiments to evaluate a team of two different robots, including an unmanned aerial vehicle (UAV) and an unmanned surface vehicle (USV), to collect hydrologic surface data utilizing sonar and visual sensors were conducted at three different field sites within the Arkavathy Basin river network located near Bangalore in Karnataka, South India. Visual sensors were used on the UAV to capture high resolution imagery for topographic characterization, and sonar sensors were deployed on the USV to capture bathymetric readings; the data streams were fused in an automated workflow to determine the storage capacity of agricultural reservoirs (also known as ``tanks'') at the three field sites. This study suggests: (i) this robot-assisted methodology is low-cost and suitable for novice users, and (ii) storage capacity data collected at previously unmapped locations revealed strong power-type relationships between surface area, stage, and storage volume, which can be incorporated into modeling of landscape-scale hydrology. This methodology is of importance to water researchers and practitioners because it produces local, high-resolution representations of bathymetry and topography and enables water balance computations at small-watershed scales, which offer insight into the present-day dynamics of a strongly human impacted watershed.

Groundwater fluxes in a shallow seasonal wetland pond: The effect of bathymetric uncertainty on predicted water and solute balances

The successful management of groundwater dependent shallow seasonal wetlands requires a sound understanding of groundwater fluxes. However, such fluxes are hard to quantify. Water volume and solute mass balance models can be used in order to derive an estimate of groundwater fluxes within such systems. This approach is particularly attractive, as it can be undertaken using measurable environmental variables, such as; rainfall, evaporation, pond level and salinity. Groundwater fluxes estimated from such an approach are subject to uncertainty in the measured variables as well as in the process representation and in parameters within the model. However, the shallow nature of seasonal wetland ponds means water volume and surface area can change rapidly and non-linearly with depth, requiring an accurate representation of the wetland pond bathymetry. Unfortunately, detailed bathymetry is rarely available and simplifying assumptions regarding the bathymetry have to be made. However, the implications of these assumptions are typically not quantified. We systematically quantify the uncertainty implications for eight different representations of wetland bathymetry for a shallow seasonal wetland pond in South Australia. The predictive uncertainty estimation methods provided in the Model-Independent Parameter Estimation and Uncertainty Analysis software (PEST) are used to quantify the effect of bathymetric uncertainty on the modelled fluxes. We demonstrate that bathymetry can be successfully represented within the model in a simple parametric form using a cubic Bezier curve, allowing an assessment of bathymetric uncertainty due to measurement error and survey detail on the derived groundwater fluxes compared with the fixed bathymetry models. Findings show that different bathymetry conceptualisations can result in very different mass balance components and hence process conceptualisations, despite equally good fits to observed data, potentially leading to poor management decisions for the wetlands. Model predictive uncertainty increases with the crudity of the bathymetry representation, however, approximations that capture the general shape of the wetland pond such as a power law or Bezier curve show only a small increase in prediction uncertainty compared to the full dGPS surveyed bathymetry, implying these may be sufficient for most modelling purposes.

Impact of remote forcing, model resolution and bathymetry on predictions of currents on the shelf

Impacts of remote forcing, model resolution and bathymetry on current predictions at two moorings located on the shelf of the Monterey Bay area are investigated. We consider three Monterey Bay model configurations which differ in resolution and bathymetry representation, and we specify open boundary conditions for these three configurations from two larger scale models, which have different accuracy in the representation of the remote forcing (in the form of poleward propagating along the coast coastally-trapped Kelvin type waves).Comparisons of correlations between observed and model currents as well as visual comparisons show that the most critical element in reproducing currents on the shelf is accurate representation of the remote forcing. Our results also show that accurate representation of bathymetry is the second most critical factor in reproducing observed currents.

Limiters for spectral propagation velocities in SWAN

As phase-averaged spectral wave models continue to grow in sophistication, they are applied more frequently throughout the ocean, from the generation of waves in deep water to their dissipation in the nearshore. Mesh spacings are varied within the computational domain, either through the use of nested, structured meshes or a single, unstructured mesh. This approach is economical, but it can cause accuracy errors in regions where the input parameters are under-resolved. For instance, in regions with a coarse representation of bathymetry, refraction can focus excessive wave energy at a single mesh vertex, causing the computed solution to become non-physical. Limiters based on the Courant–Friedrichs–Lewy (CFL) criteria are proposed for the spectral propagation (refraction and frequency shifting) velocities in SWAN. These limiters are not required for model stability, but they improve accuracy by reducing local errors that would otherwise spread throughout the computational domain. As demonstrated on test cases in deep and shallow water, these limiters prevent the excessive directional turning and frequency shifting of wave energy and control the largest errors in under-resolved regions.

On a momentum conservative z-layer unstructured C-grid ocean model with flooding

We present a z-layer unstructured C-grid finite volume hydrostatic model. An efficient and highly scalable implicit technique for solution of the free surface equation is combined with an Eulerian approach for the advection of momentum. We show that an accurate velocity reconstruction procedure is of crucial importance not only for discretization of the Coriolis term, but also for the correct advection of momentum, especially in the multilayer case. Unlike other z-layer models the method presented here ensures that the staircase representation of bathymetry and free surface has no influence on the vertical structure of the flow. The method is therefore guaranteed to be strictly momentum conservative, also in the layers containing the free surface and bed. A number of test cases are presented to show that the model is able to accurately simulate Coriolis dominated flows and flooding and drying processes both in the depth-averaged case and in the presence of multiple z-layers. We use a simulation of 2004 Indian Ocean tsunami to evaluate the ability of the method to simulate fast propagating tsunami waves and detailed inundation processes. Results obtained using two different rupture models are compared to the tide gauge arrival times, satellite altimetry data and the inundation observations in the Banda Aceh area. The comparison is used not only to assess the quality of the underlying rupture models but also to determine the value of the available data sources for such an assessment.

Adaptive volume penalization for ocean modeling

The development of various volume penalization techniques for use in modeling topographical features in the ocean is the focus of this paper. Due to the complicated geometry inherent in ocean boundaries, the stair-step representation used in the majority of current global ocean circulation models causes accuracy and numerical stability problems. Brinkman penalization is the basis for the methods developed here and is a numerical technique used to enforce no-slip boundary conditions through the addition of a term to the governing equations. The second aspect to this proposed approach is that all governing equations are solved on a nonuniform, adaptive grid through the use of the adaptive wavelet collocation method. This method solves the governing equations on temporally and spatially varying meshes, which allows higher effective resolution to be obtained with less computational cost. When penalization methods are coupled with the adaptive wavelet collocation method, the flow near the boundary can be well-resolved. It is especially useful for simulations of boundary currents and tsunamis, where flow near the boundary is important. This paper will give a thorough analysis of these methods applied to the shallow water equations, as well as some preliminary work applying these methods to volume penalization for bathymetry representation for use in either the nonhydrostatic or hydrostatic primitive equations.

Computations of the energy flux to mixing processes via baroclinic wave drag on barotropic tides

The geographical distribution of barotropic to baroclinic transfer of tidal energy by baroclinic wave drag in the abyssal ocean is estimated. Using tidal velocities from a state-of-the-art numerical tidal model, the total loss of barotropic tidal energy in the deep ocean (between 70°S and 70°N and at depths greater than 1000 m) is estimated to be about 0.7 TW ( M 2) corresponding to a mean value of the energy flux ( e) of 2.4×10 −3 W/m 2. The distribution of e is however highly skewed with a median of about 10 −6 W/m 2. Only 10% of the area is responsible for more than 97% of the total energy transfer. To assess the possible influence of the relatively coarse bathymetry representation upon the present estimate, complementary calculations using better resolved sea floor topography are carried out over a control area around the Hawaiian Ridge. There are no major differences between the results achieved using the two different bathymetry databases. Fluxes of about 16 GW or 6×10 −3 W/m 2 are computed in both cases, and the main contributions to the total fluxes originate in the same range of e-values and cover equally large parts of the total area. It is not clear whether the present model is valid at flat or subcritical bottom slopes. However, for the Hawaiian region, only 2% of the total energy flux as calculated in the present study originates in areas of critical and subcritical slopes.