LAGEOS Sensitivity to Ocean Tides

Satellite Laser Ranging (SLR) to LAGEOS has a remarkable contribution to high-precise geodesy and geodynamics through deriving and validating various global geophysical models. This paper validates ocean tide models based on the analysis of satellite altimetry data, coastal tide gauges, and hydrodynamic data, i.e., CSR3.0, TOPEX4.0, CSR4.0A, FES2004, GOT00.2, and the CSRC Schwiderski model. LAGEOS orbits and SLR observation residuals from solutions based on different ocean tide models are compared and examined. It is found that LAGEOS orbits are sensitive to tidal waves larger than 5 mm. The analysis of the aliasing periods of LAGEOS orbits and tidal waves reveals that, in particular, the tidal constituent S₂ is not well established in the recent ocean tide models. Some of the models introduce spurious peaks to empirical orbit parameters, which can be associated with S₂, S_a, and K₂ tidal constituents, and, as a consequence, can be propagated to fundamental parameters derived from LAGEOS observations.     

An estimate of the suspended particulate matter (SPM) transport in the southern North Sea using SeaWiFS images,in situ measurements and numerical model results

A study is presented where satellite images (SeaWiFS), in situ measurements (tidal cycle and snapshot) and a 2D hydrodynamic numerical model have been combined to calculate the long term SPM (Suspended Particulate Matter) transport through the Dover Strait and in the southern North Sea. The total amount of SPM supplied to the North Sea through the Dover Strait is estimated to be 31.74×10⁶ t. The satellite images provide synoptic views of SPM concentration distribution but do not take away the uncertainty of SPM transport calculation. This is due to the fact that SPM concentration varies as a function of tide, wind, spring-neap tidal cycles and seasons. The short term variations (tidal, spring-neap tidal cycle) have not been found in the satellite images, however seasonal variations are clearly visible. Furthermore the SPM concentration in the satellite images is generally lower than in the in situ measurements. The representativness of SPM concentration maps derived from satellites for calculating long term transports has therefore been investigated by comparing the SPM concentration variability from the in situ measurements with those of the remote sensing data. The most important constraints of satellite images are related to the fact that satellite data is evidence of clear sky conditions, whereas in situ measurements from a vessel can be carried out also during rougher meteorological conditions and that due to the too low time resolution of the satellite images the SPM concentration peaks are often missed. It is underlined that SPM concentration measurements should be carried out during at least one tidal cycle in high turbidity areas to obtain representative values of SPM concentration.     

Coastal tides in West Greenland derived from tide gauge records

Based on pressure tide-gauge observations, sea-level records are derived for ten sites along the coast of West Greenland. The ocean tidal signal is extracted by a harmonic tidal analysis. The accuracy of the determined tidal constants is discussed in detail. The tides account for 85% of the observed sea-level standard deviation. The tide gauge records reveal significant shallow-water tidal effects, in particular compound and overtide amplitudes reaching 5 cm. The propagation of the tidal waves into the fjords depends strongly on local conditions and is in some cases accompanied by an amplification of the tidal amplitudes. The observed tidal signals are compared to the predictions of the global ocean tide model FES2004. At the outer coast, a good agreement is found. Inside the fjords, however, the model performs worse and tide gauge observations may still be indispensable when accurate tidal signals are required.     

Numerical Modelling in Geo-Electromagnetics: Advances and Challenges

During the last decade, tremendous advances have been observed in the broad field of numerical modelling for geo-electromagnetic applications. This trend received support due to increasing industrial needs, mainly caused by hydrocarbon and ore exploration industry. On the other hand, the increasing reliability and accuracy of data acquisition techniques further spurs this development. In this review, we will focus on advances and challenges in numerical modelling in geo-electromagnetics. We review recent developments in the discrete solution of the 3-D induction problem in the time and frequency domains. Particularly, advantages and disadvantages of the common numerical techniques for solving partial differential equations such as the Finite Difference and Finite Element methods will be considered.     

Tidal inlet response to sediment infilling of the associated bay and possible implications of human activities: the Marennes-Oléron Bay and the Maumusson Inlet,France

Tidal inlet characteristics are controlled by wave energy, tidal range, tidal prism, sediment supply and direction and rates of sand delivered to the inlet. This paper deals with the relations between inlet and lagoon evolutions, linked by the tidal prism. Our study is focused on the Maumusson Inlet and the Marennes-Oléron Bay (first oyster farming area in Europe), located on the western coast of France. The tidal range (2–6 m) and wave climate (mean height: 1.5 m) place this tidal inlet system in the mixed energy (tide, waves), tide-dominated category. The availability of high-resolution bathymetric data since 1824 permits to characterise and quantify accurately morphological changes of both the inlet and the tidal bay. Since 1824, sediment filling of the tidal bay has led to a 20% decrease in its water volume, and a 35% reduction of the inlet throat section. Furthermore, the bay is subjected to a very high anthropic pressure, mainly related to oyster farming. Thus, both natural and human-related processes seem relevant to explain high sedimentation rates. Current measurements, hydrodynamic modelling and cross-sectional area of the inlet throat are used in order to quantify tidal prism changes since 1824. Both flood and ebb tidal prism decreased by 35%. Decrease in the Marennes-Oléron Bay water volume is inferred to be responsible for a part of tidal prism decrease at the inlet. Tidal prisms decrease may also be explained by an increase in frictional resistance to tidal wave propagation, due to a general shoaling and oyster farms in the bay. A conceptual model is proposed, taking into account natural and human-related sedimentation processes, and explaining tidal inlet response to tidal bay evolutions.     

On the stability relationships between tidal asymmetry and morphologies of tidal basins and estuaries

Simple stability relationships are practically useful to provide a rapid assessment of coastal and estuarine landforms in response to human interventions and long‐term climate change. In this contribution, we review a variety of simple stability relationships which are based on the analysis of tidal asymmetry (TA). Most of the existing TA‐based stability relationships are derived using the one‐dimensional tidal flow equations assuming a certain regular shape of the tidal channel cross‐sections. To facilitate analytical solutions, specific assumptions inevitably need to be made, for example by linearizing the friction term and dropping some negligible terms in the tidal flow equations. We find that three major types of TA‐based stability relationships have been proposed between three non‐dimensional channel geometric ratios (represented by the ratio of channel widths, ratio of wet surface areas and ratio of storage volumes) and the tide‐related parameter a/h (i.e. the ratio between tidal amplitude and mean water depth). Based on established geometric relations, we use these non‐dimensional ratios to restate the existing relationships so that they are directly comparable. Available datasets are further extended to examine the utility of these TA‐based relationships. Although a certain agreement is shown for these relationships, we also observe a large scatter of data points which are collected in different types of landscape, hydrodynamic and sedimentological settings over the world. We discuss in detail the potential reasons for this large scatter and subsequently elaborate on the limited applicability of the various TA‐based stability relationships for practical use. We highlight the need to delve further into what constitutes equilibrium and what is needed to develop more robust measures to determine the morphological state of these systems. Copyright © 2018 John Wiley & Sons, Ltd.     

Yardsticks for industrial tomography

Seismic tomography has been developed and applied for decades in seismological applications and for basic research purposes. During the last decade, large-scale applications in the oil and gas industry became standard as tomostatics and velocity modelling for pre-stack depth migration. In this paper, I take a snapshot of some current industrial applications, quantifying practical aspects by yardsticks such as data and model size and I try to draw a road map for the current decade.     

Interpretation of tensor gravity data using an adaptive tilt angle method

Full Tensor Gravity Gradiometry (FTG) data are routinely used in exploration programmes to evaluate and explore geological complexities hosting hydrocarbon and mineral resources. FTG data are typically used to map a host structure and locate target responses of interest using a myriad of imaging techniques. Identified anomalies of interest are then examined using 2D and 3D forward and inverse modelling methods for depth estimation. However, such methods tend to be time consuming and reliant on an independent constraint for clarification. This paper presents a semi‐automatic method to interpret FTG data using an adaptive tilt angle approach. The present method uses only the three vertical tensor components of the FTG data (T_zx, T_zy and T_zz) with a scale value that is related to the nature of the source (point anomaly or linear anomaly). With this adaptation, it is possible to estimate the location and depth of simple buried gravity sources such as point masses, line masses and vertical and horizontal thin sheets, provided that these sources exist in isolation and that the FTG data have been sufficiently filtered to minimize the influence of noise. Computation times are fast, producing plausible results of single solution depth estimates t hat relate directly to anomalies. For thick sheets, the method can resolve the thickness of these layers assuming the depth to the top is known from drilling or other independent geophysical data. We demonstrate the practical utility of the method using examples of FTG data acquired over the Vinton Salt Dome, Louisiana, USA and basalt flows in the Faeroe‐Shetland Basin, UK. A major benefit of the method is the ability to quickly construct depth maps. Such results are used to produce best estimate initial depth to source maps that can act as initial models for any detailed quantitative modelling exercises using 2D/3D forward/inverse modelling techniques.     

Field measurements of intertidal bar evolution on a high‐energy beach system

Nearshore bars play a pivotal role in coastal behaviour, helping to protect and restore beach systems particularly in post‐storm conditions. Examination of bar behaviour under various forcing conditions is important to help understand the short‐ to medium‐term evolution of sandy beach systems. This study carried out over a nine‐week period examines, the behaviour of three intertidal bars along a high energy sandy beach system in northwest Ireland using high‐frequency topographic surveys and detailed nearshore hydrodynamic modelling. Results show that, in general, there was onshore migration for all the bars during the study period, despite the variability observed between bars, which was driven mostly by wave dominated processes. Under the prevailing conditions migration rates of up to 1.83 m day^?1 and as low as 0.07 m day^?1 were observed. During higher wave energy events the migration rates of the bars decelerated in their onshore route, however, under lower wave energy conditions, they quickly accelerated maintaining their shoreward migration direction. Tidal influence appears to be subordinate in these conditions, being restricted to moderating the localized wave energy at low tides and in maintaining runnel configurations providing accommodation space for advancing slip faces. The study highlights the intricate behavioural patterns of intertidal bar behaviour along a high energy sandy coastline and provides new insights into the relative importance of wave and tidal forcing on bar behaviour over a relatively short time period. Copyright © 2016 John Wiley & Sons, Ltd.     

Tidal dissipation from free drift sea ice in the Barents Sea assessed using GNSS beacon observations

Ocean Dynamics - One of the major challenges facing global hydrodynamic tidal models is the modelling of the interaction between sea ice and tides in high-latitude waters. Recent studies have shown...     

Impact of tidal power dams upon tides and environmental conditions in the Sea of Okhotsk

The transformation of natural tidal sea-level and currents is studied resulting from large-scale tidal power plant (TPP) dams in bays of the Sea of Okhotsk (SO). Some effects due to this transformation are estimated based on predictive modelling and a number of expected changes in amplitudes and phases, and spectral composition of tidal oscillations are described. Changes of morphometric properties of basins change the character of tidal motions even on significant distance from a dam. That is why, it is impossible to estimate this impact as usual boundary-value problems. The problem is solved based on “impedance” conditions on the open boundary of the model area, allowing to take into account the radiation of the additional perturbations induced by both waves reflected from the dam and nonlinear effects inside the area. In general, the transformation effects are proportional to the dam size and depend essentially on the dam location, the creation of which can change dissipative and resonance properties of the bays. The changes in tidal energetics of SO due to the dam construction are also considered to show noticeable reconstruction of horizontal energy fluxes and changes in the energy dissipation. Possible environmental consequences are related mainly to the transformation of tidal currents.     

Finite element modeling of the July 12, 1993 Hokkaido Nansei-Oki tsunami

A fault plane model and a finite element hydrodynamic model are applied to the simulation of the Hokkaido Nansei-Oki tsunami of July 12, 1993. The joint performance of the models is assessed based on the overall ability to reproduce observed tsunami waveforms and to preserve mass and energy during tsunami propagation. While a number of observed characteristics of the waveforms are satisfactorily reproduced (in particular, amplitudes and arrival times at tidal gauges relatively close to the source, and general patterns of energy concentration), others are only marginally so (notably, wave periods at the same gauges, and wave heights along Okushiri); differences between observations and simulations are traceable to both the fault plane and the hydrodynamic models. Nonnegligible losses of energy occur throughout the simulated tsunami propagation. These losses seem to be due to a combination of factors, including numerical damping and possible deficiencies of the shallow water equations in preserving energy.     

Hydromechanics for the formation and development of radial sandbanks (I)

Based upon the long-term observation of field data, a two-dimensional numerical model is applied to simulating the tidal flow covering from the neap tide to spring tide in the radial sandbank area in the southern Yellow Sea. From the development of tidal current ridges under the hydrodynamic action, multi-purpose analysis and study are carried out, which include the propagation process of tidal wave, the distributions of tidal wave energy rate and tidal range, the tidal ellipses and traces. It is shown that the tidal current is the major dynamic factor for the formation and development of the radial sandbanks, and the differences of tidal wave energy rate and current strength determine the distinct plane shapes of ridges and troughs in this region.     

海潮对卫星重力场恢复的影响

本文讨论了海潮对卫星重力测量的影响问题. 首先介绍了海潮对卫星重力测量影响的基本理论；采用FES02和TPXO6海潮模型计算了海潮负荷对卫星重力结果前60阶的影响；并用两个模型之间的差异作为海潮模型精度的估计量，据此计算了海潮模型误差对卫星重力结果的影响. 与GRACE恢复的重力场精度的比较说明：海潮对重力场40阶以下的影响都超过了目前重力场恢复精度；尽管由于卫星测高技术的发展，海潮模型的精度有了很大的提高，但目前的全球海潮模型用于GRACE重力场恢复的前12阶的改正还是不够精确. 另外，我们也利用中国东海和南海潮汐资料以及FES02海潮模型讨论了中国近海潮汐效应对GRACE观测的影响. 结果说明该影响与海潮模型的误差相当. 这反映了当前海潮模型的不确定度，因此通过结合全球验潮站资料有望提高海潮对卫星重力测量的改正精度.     

Parameterization of wind turbine impacts on hydrodynamics and sediment transport

Monopile foundations of offshore wind turbines modify the hydrodynamics and sediment transport at local and regional scales. The aim of this work is to assess these modifications and to parameterize them in a regional model. In the present study, this is achieved through a regional circulation model, coupled with a sediment transport module, using two approaches. One approach is to explicitly model the monopiles in the mesh as dry cells, and the other is to parameterize them by adding a drag force term to the momentum and turbulence equations. Idealised cases are run using hydrodynamical conditions and sediment grain sizes typical from the area located off Courseulles-sur-Mer (Normandy, France), where an offshore windfarm is under planning, to assess the capacity of the model to reproduce the effect of the monopile on the environment. Then, the model is applied to a real configuration on an area including the future offshore windfarm of Courseulles-sur-Mer. Four monopiles are represented in the model using both approaches, and modifications of the hydrodynamics and sediment transport are assessed over a tidal cycle. In relation to local hydrodynamic effects, it is observed that currents increase at the side of the monopile and decrease in front of and downstream of the monopile. In relation to sediment transport effect, the results show that resuspension and erosion occur around the monopile in locations where the current speed increases due to the monopile presence, and sediments deposit downstream where the bed shear stress is lower. During the tidal cycle, wakes downstream of the monopile reach the following monopile and modify the velocity magnitude and suspended sediment concentration patterns around the second monopile.     

Uncertainty Analysis of Strong-Motion and Seismic Hazard

In this article we present the modelling of uncertainty in strong-motion studies for engineering applications, particularly for the assessment of earthquake hazard. We examine and quantify the sources of uncertainty in the basic variables involved in ground motion estimation equations, including those associated with the seismological parameters, which we derive from a considerable number of strong-motion records. Models derived from regression analysis result in ground motion equations with uncertain parameters, which are directly related to the selected basic variables thus providing an uncertainty measure for the derivative variable. These uncertainties are exemplified and quantified. An alternative approach is presented which is based on theoretical modelling defining a functional relationship on a set of independent basic variables. Uncertainty in the derivative variable is then readily obtained when the uncertainties of the basic variables have been defined. In order to simplify the presentation, only the case of shallow strike-slip earthquakes is presented. We conclude that the uncertainty is approximately the same as given by the residuals typical for regression modelling. This implies that uncertainty in ground motion modelling cannot be reduced below certain limits, which is in accordance with findings reported in the literature. Finally we discuss the implications of the presented methodology in hazard analyses, which is sensitive to the truncation of the internal error term, commonly given as an integral part of ground motion estimation equations. The presented methodology does not suffer from this shortcoming; it does not require truncation of the error term and yields realistic hazard estimates. This revised version was published online in July 2006 with corrections to the Cover Date.     

Scales and structure of frontal adjustment and freshwater export in a region of freshwater influence

Sea surface temperature satellite imagery and a regional hydrodynamic model are used to investigate the variability and structure of the Liverpool Bay thermohaline front. A statistically based water mass classification technique is used to locate the front in both data sets. The front moves between 5 and 35 km in response to spring–neap changes in tidal mixing, an adjustment that is much greater than at other shelf-sea fronts. Superimposed on top of this fortnightly cycle are semi-diurnal movements of 5–10 km driven by flood and ebb tidal currents. Seasonal variability in the freshwater discharge and the density difference between buoyant inflow and more saline Irish Sea water give rise to two different dynamical regimes. During winter, when cold inflow reduces the buoyancy of the plume, a bottom-advected front develops. Over the summer, when warm river water provides additional buoyancy, a surface-advected plume detaches from the bottom and propagates much larger distances across the bay. Decoupled from near-bed processes, the position of the surface front is more variable. Fortnightly stratification and re-mixing over large areas of Liverpool Bay is a potentially important mechanism by which freshwater, and its nutrient and pollutant loads, are exported from the coastal plume system. Based on length scales estimated from model and satellite data, the erosion of post-neap stratification is estimated to be responsible for exporting approximately 19% of the fresh estuarine discharge annually entering the system. Although the baroclinic residual circulation makes a more significant contribution to freshwater fluxes, the episodic nature of the spring–neap cycle may have important implications for biogeochemical cycles within the bay.     

Coupled 1D–3D hydrodynamic modelling,with application to the Pearl River Delta

Within the hydrodynamic modelling community, it is common practice to apply different modelling systems for coastal waters and river systems. Whereas for coastal waters 3D finite difference or finite element grids are commonly used, river systems are generally modelled using 1D networks. Each of these systems is tailored towards specific applications. Three-dimensional coastal water models are designed to model the horizontal and vertical variability in coastal waters and are less well suited for representing the complex geometry and cross-sectional areas of river networks. On the other hand, 1D river network models are designed to accurately represent complex river network geometries and complex structures like weirs, barrages and dams. A disadvantage, however, is that they are unable to resolve complex spatial flow variability. In real life, however, coastal oceans and rivers interact. In deltaic estuaries, both tidal intrusion of seawater into the upstream river network and river discharge into open waters play a role. This is frequently approached by modelling the systems independently, with off-line coupling of the lateral boundary forcing. This implies that the river and the coastal model run sequentially, providing lateral discharge (1D) and water level (3D) forcing to each other without the possibility of direct feedback or interaction between these processes. An additional disadvantage is that due to the time aggregation usually applied to exchanged quantities, mass conservation is difficult to ensure. In this paper, we propose an approach that couples a 3D hydrodynamic modelling system for coastal waters (Delft3D) with a 1D modelling system for river hydraulics (SOBEK) online. This implies that contrary to off-line coupling, the hydrodynamic quantities are exchanged between the 1D and 3D domains during runtime to resolve the real-time exchange and interaction between the coastal waters and river network. This allows for accurate and mass conserving modelling of complex coastal waters and river network systems, whilst the advantages of both systems are maintained and used in an optimal and computationally efficient way. The coupled 1D–3D system is used to model the flows in the Pearl River Delta (Guangdong, China), which are determined by the interaction of the upstream network of the Pearl River and the open waters of the South China Sea. The highly complex upstream river network is modelled in 1D, simulating river discharges for the dry and wet monsoon periods. The 3D coastal model simulates the flow due to the external (ocean) periodic tidal forcing, the salinity distribution for both dry and wet seasons, as well as residual water levels (sea level anomalies) originating from the South China Sea. The model is calibrated and its performance extensively assessed against field measurements, resulting in a mean root mean square (RMS) error of below 6% for water levels over the entire Pearl River Delta. The model also represents both the discharge distribution over the river network and salinity transport processes with good accuracy, resolving the discharge distribution over the main branches of the river network within 5% of reported annual mean values and RMS errors for salinity in the range of 2 ppt (dry season) to 5 ppt (wet season).     

Modeling of dam break wave propagation in a partially ice-covered channel

During the last four decades, several numerical formulations and specialized software have been developed in response to studies about dam break (DB) wave propagation and its hydraulic and environmental impacts on downstream hydraulic structures and valleys. These methods cannot, however, be used to predict wave propagation within partially covered channels or reservoirs located upstream of hydraulic structures. In fact, such problems require the modelling of the complex transition from a free surface flow into a pressurized one. Because rivers or channels partially covered with ice sheets are typical examples commonly met in winter in northern climates, it is vitally important to assess ice-cover effects on the DB wave propagation and develop appropriate tools to predict resulting hydrodynamic loads on hydraulic structures downstream. This paper proposes an original numerical formulation to model wave propagation and hydrodynamic pressure in partially covered channels. The proposed formulation uses one-dimensional St. Venant equations to simulate open-water flow and water hammer equations to simulate pressure flow within the partially covered channel. To illustrate the use of the hydrodynamic pressures obtained, a case study is presented where a channel cover and a dam located downstream are modelled using finite elements to investigate their dynamic structural response.     

A model of sand transport in Treporti channel: northern Venice lagoon

The origin of the sands in the Venice lagoon has been the subject of an extensive field survey in parallel with numerical modelling. Four transects along Treporti and Burano canals were conducted from which 33 bottom sediment samples were collected. These samples were analysed for grain size and sorting to examine any trends in the granulometry of these sediments that might shed light on transport paths. The modelling study consists of three parts: the sediment transport model sedtrans96 was used with a finite-element hydrodynamic model (Shyfem) and an empirical wave model (US Army Corps of Engineering) to simulate sand transport in the Treporti canal. A type of link box model was created where finite elements of the hydrodynamic model have been combined to macro-boxes on which the water and sediment flux over the sections, and a mass balance has been computed. Several grain size classes were simulated; the distributions before and after the simulation were examined. Idealised wind and tidal values were initially used to force 12 h simulations to test the sediment transport sensitivity. Finally, a full-year simulation (1987) has been carried out using measured tidal and wind data. Only a part of Venice lagoon was covered by the simulation: a major channel (Treporti) running from Lido inlet towards the northern lagoon. The total sand transport through all of the sections was computed for 1 year. Sediment mass balance was determined, and the resulting trends of erosion and deposition were computed. There were no trends in the median grain diameter and sorting of bottom samples from the Treporti canal; all sands were fine (120 μm, one outlier of 300 μm was removed). The absence of a trend in grain size suggests that there is no significant import of sand to the lagoon through the Lido inlet. The results from the simulations seem therefore to confirm the hypothesis of reworking of sand within the lagoon. The computed erosion is some centimeters per year diagnostic of channel scouring and enlargement with time. The Treporti canal is subject to strong current velocities of around 1 m/s, which hold fine sand in suspension and thus prevent sedimentation.