Horizontal patterns of water temperature and salinity in an estuarine tidal channel: Ria de Aveiro

This work presents results from two complementary and interconnected approaches to study water temperature and salinity patterns in an estuarine tidal channel. This channel is one of the four main branches of the Ria de Aveiro, a shallow lagoon located in the Northwest coast of the Iberian Peninsula. Longitudinal and cross-sectional fields of water temperature and salinity were determined by spatial interpolation of field measurements. A numerical model (Mohid) was used in a 2D depth-integrated mode in order to compute water temperature and salinity patterns. The main purpose of this work was to determine the horizontal patterns of water temperature and salinity in the study area, evaluating the effects of the main forcing factors. The field results were depth-integrated and compared to numerical model results. These results obtained using extreme tidal and river runoff forcing, are also presented. The field results reveal that, when the river flow is weak, the tidal intrusion is the main forcing mechanism, generating saline and thermal fronts which migrate with the neap/spring tidal cycle. When the river flow increases, the influence of the freshwater extends almost as far as the mouth of the lagoon and vertical stratification is established. Results of numerical modelling reveal that the implemented model reproduces quite well the observed horizontal patterns. The model was also used to study the hydrology of the study area under extreme forcing conditions. When the model is forced with a low river flow (1 m³ s⁻¹) the results confirm that the hydrology is tidally dominated. When the model is forced with a high river flow (1,000 m³ s⁻¹) the hydrology is dominated by freshwater, as would be expected in such an area.     

Model simulations of tidal changes in a coastal lagoon,the Ria de Aveiro (Portugal)

Tidal changes in the Ria de Aveiro lagoon, Portugal, show an average increase of 0.245 m in M₂ amplitude and an average 17.4° decrease in M₂ phase, over 16 years (1987–2004). This paper investigates the causes of these changes using: (a) an analytical model and (b) hydrodynamic (vertically integrated) model of the Ria de Aveiro lagoon.     

A numerical system to study the transport properties in the Ria de Aveiro lagoon

This work presents results from numerical modelling studies on the hydrodynamics and sediments and passive particles transport properties in Ria de Aveiro, a shallow lagoon located on the Portuguese Northwest Atlantic coast. The hydrodynamic of the lagoon was systematically studied, from both the Eulerian and Lagragian point of view, in order to understand the overall circulation in the lagoon, characterize the dynamics of its different channels and assess the transport of suspended cohesive sediments.Responsible Editor: Hans Burchard     

Tidal flow asymmetry in the diurnal regime: bed-load transport and morphologic changes around the Red River Delta

The relation between tidal flow asymmetry and net transport of sediment in the semidiurnal regime has been extensively described. This study reveals that in the diurnal regime, the direction of long-term net bed-load transport and resulting morphologic changes is partly determined by the phase-angle relationship of O₁, K₁, and M₂. Simple analytical formulations of time-averaged bed-load transport were derived which separate the relative contributions of tidal asymmetry from that of residual flow with tidal stirring. In this particular case, the Red River Delta in Vietnam, transports related to tidal asymmetry are larger than those induced by the monsoon currents, and are an order of magnitude larger than those associated with topographic residual flow around the delta. Tide-induced morphologic changes dominate at water depths between 10 and 25 m, of which the patterns of erosion and deposition overlap with observed bathymetric changes. Additional observed changes that occur in more shallow water cannot be explained by tidal asymmetry and are probably related to wave action and to deposition from the buoyant river plume.Responsible Editor: Jens Kappenberg     

SST observations of upwelling induced by tidal straining in the Rhine ROFI

In this paper SST imagery and a three-dimensional numerical model of a river plume were employed to detect upwelling induced by tidal straining in the Rhine ROFI (region of fresh water influence). Previous studies have shown that the Rhine ROFI in the North Sea exhibits strong cross-shore density gradients that compete with tidal and wind mixing to establish stratification. During neap periods with low mixing energy an area measuring 30 km offshore by 100 km alongshore becomes stratified. When the ROFI is stratified strong cross-shore currents are observed, with surface currents rotating anti-cyclonically and bottom currents rotating cyclonically. The cross-shore currents interact with the cross-shore density gradients to produce a semi-diurnal cycle of stratification. Due to continuity requirements imposed by the proximity of the coast, the offshore-directed surface currents and onshore-directed bottom currents should lead to coastal upwelling.     

An idealized model study of flocculation on sediment trapping in an estuarine turbidity maximum

Flocculation has an important impact on particle trapping in estuarine turbidity maximum (ETM) through associated increases in particle settling velocity. To quantify the importance of the flocculation processes, a size-resolved flocculation model is implemented into an ocean circulation model to simulate fine-grained particle trapping in an ETM. The model resolves the particle size from robust small flocs, about 30 μm, to very large flocs, over 1000 μm. An idealized two-dimensional model study is performed to simulate along-channel variations of suspended sediment concentrations driven by gravitational circulation and tidal currents. The results indicate that the flocculation processes play a key role in generating strong tidal asymmetrical variations in suspended sediment concentration and particle trapping. Comparison with observations suggests that the flocculation model produces realistic characteristics of an ETM.     

Tidal currents,energy flux and bottom boundary layer thickness in the Clyde Sea and North Channel of the Irish Sea

A high-resolution three-dimensional model of the Clyde Sea and the adjacent North Channel of the Irish Sea is used to compute the major diurnal and semidiurnal tides in the region, the associated energy fluxes and thickness of the bottom boundary layer. Initially, the accuracy of the model is assessed by performing a detailed comparison of computed tidal elevations and currents in the region, against an extensive database that exists for the M₂, S₂, N₂, K₁ and O₁ tides. Subsequently, the model is used to compute the tidal energy flux vectors in the region. These show that the major energy flux is confined to the North Channel region, with little energy flux into the Clyde Sea. Comparison with the observed energy flux in the North Channel shows that its across-channel distribution and its magnitude are particularly sensitive to the phase difference between elevation and current. Consequently, small changes in the computed values of these parameters due to slight changes of the order of the uncertainty in the open-boundary values to the model, can significantly influence the computed energy flux. The thickness of the bottom boundary layer in the region is computed using a number of formulations. Depending upon the definition adopted, the empirical coefficient C used to determine its thickness varies over the range 0.1 to 0.3, in good agreement with values found in the literature. In the North Channel, the boundary layer thickness occupies the whole water depth, and hence tidal turbulence produced at the sea bed keeps the region well mixed. In the Clyde Sea, the boundary layer thickness is a small fraction of the depth, and hence the region stratifies.Responsible Editor: Phil Dyke     

Complete Coulomb stress changes induced by the Ms7.6 earthquake in Lancang-Gengma,Yunnan and triggering of aftershocks by dynamic and static stress

The spatiotemporal evolution patterns of complete Coulomb stress changes caused by 1988 Ms7.6 earthquake in Lancang-Gengma,Yunnan,are calculated and studied.And the triggering problems of Ms7.2 Gengma shock occurring 13 minutes after the main shock and of Ms5.0―6.9 aftershocks within 24 days after the main shock are discussed.The results show that the spatial distribution patterns of complete Coulomb stress changes of the Ms7.6 main shock are strongly asymmetric.The areas of positive dynamic and static Coulomb stress are both coincident well with the strong aftershocks' locations.The Ms7.2 Gengma shock and most of strong aftershocks are subjected to the triggering effect of dynamic and static Coulomb stresses induced by the Ms7.6 Lancang earthquake.     

云南武定M6.5地震动态库仑破裂应力变化的数值模拟及其与余震活动的关系

在离散波数法（DWN）基础上,计算了武定M6.5地震断层破裂在周围介质中产生的位移时程（位移理论地震图）和动态位移场;进行弹性动力学转换后,求得应变时程和动态应变场;最后得到了武定M6.5地震所产生的动态库仑破裂应力变化量和动态库仑破裂应力变化场,进而研究其与后续余震的关系. 结果表明,动态应力最大峰值和静态应力的正区均呈非对称性分布,两者的分布特征与余震的分布特征基本一致. 在动态应力峰值为正的确定区应力值超过了0.1 MPa触发阈值, 在静态应力值为正的确定区应力值超过了0.01 MPa触发阈值. 这说明动态应力和静态应力均有助于余震的发生.      

Numerical simulation of dynamic Coulomb stress changes induced by M6.5 earthquake in Wuding, Yunnan and its relationship with aftershocks

Based on the discrete wavenumber method, we calculate the fields of dynamic Coulomb rupture stress changes and static stress changes caused by M6.5 earthquake in Wuding, and study their relationship with the subsequent after- shocks. The results show that the spatial distribution patterns of the positive region of dynamic stress peak value and static stress peak value are similarly asymmetric, which are basically identical with distribution features of aftershock. The dynamic stress peak value and the stati...     

Numerical simulation of dynamic Coulomb stress changes induced by M 6.5 earthquake in Wuding, Yunnan and its relationship with aftershocks

Based on the discrete wavenumber method, we calculate the fields of dynamic Coulomb rupture stress changes and static stress changes caused by M6.5 earthquake in Wuding, and study their relationship with the subsequent aftershocks. The results show that the spatial distribution patterns of the positive region of dynamic stress peak value and static stress peak value are similarly asymmetric, which are basically identical with distribution features of aftershock. The dynamic stress peak value and the static stress in the positive region are more than 0.1 MPa and 0.01 MPa of the triggering threshold, respectively, which indicates that the dynamic and static stresses are helpful for the occurrence of aftershock. This suggests that both influences of dynamic and static stresses should be considered other than only either of them when studying aftershock triggering in near field. Foundation item: National Natural Science Foundation of China (40364001) and Joint Seismological Science Foundation of China (605014).     

Analysis of flood propagation changes in the Kienstock—Bratislava reach of the Danube River/Analyse des changements de propagation des crues dans le tronçon Kienstock—Bratislava du Fleuve Danube

Abstract

In the first part of this study, a flood wave transformation analysis for the largest historical floods in the Danube River reach Kienstock–Bratislava was carried out. For the simulation of the historical (1899 and 1954) flood propagation, the nonlinear river model NLN-Danube (calibrated on the recent river reach conditions) was used. It was shown that the simulated peak discharges were not changed significantly when compared to their historical counterparts. However, the simulated hydrographs exhibit a significant acceleration of the flood wave movement at discharges of between 5000 and 9000 m³ s^-1. In the second part, the travel time-water level relationships between Kienstock and Bratislava were analysed on a dataset of the flood peak water levels for the period 1991–2002. An empirical regression routing scheme for the Danube short-term water level forecast at Bratislava station was derived. This is based on the measured water level at Kienstock gauging station.