Three-dimensional Hydrodynamic Model Research Articles

Realistic 3D hydrodynamics simulations find significant turbulent entrainment in massive stars

ABSTRACT Our understanding of stellar structure and evolution coming from one-dimensional (1D) stellar models is limited by uncertainties related to multidimensional processes taking place in stellar interiors. 1D models, however, can now be tested and improved with the help of detailed three-dimensional (3D) hydrodynamics models, which can reproduce complex multidimensional processes over short time-scales, thanks to the recent advances in computing resources. Among these processes, turbulent entrainment leading to mixing across convective boundaries is one of the least understood and most impactful. Here, we present the results from a set of hydrodynamics simulations of the neon-burning shell in a massive star, and interpret them in the framework of the turbulent entrainment law from geophysics. Our simulations differ from previous studies in their unprecedented degree of realism in reproducing the stellar environment. Importantly, the strong entrainment found in the simulations highlights the major flaws of the current implementation of convective boundary mixing in 1D stellar models. This study therefore calls for major revisions of how convective boundaries are modelled in 1D, and in particular the implementation of entrainment in these models. This will have important implications for supernova theory, nucleosynthesis, neutron stars, and black holes physics.

Modelling the Dispersion of Seafloor Massive Sulphide Mining Plumes in the Mid Atlantic Ridge Around the Azores

It is increasingly recognised that deep-sea mining of seafloor massive sulphides (SMS) could become an important source of mineral resources. These operations will remove the targeted substrate and produce potentially toxic plumes from in situ seabed excavation and from the return water pumped back down to the seafloor. However, the spatial extent of the impact of deep-sea mining is still uncertain because few field experiments and models of plume dispersion have been conducted. In this study, we used three-dimensional hydrodynamic models of the Azores region together with a theoretical commercial mining operation of polymetallic SMS to simulate the potential dispersal of plumes originating from different phases of mining operations, and to assess the magnitude of potential impacts. Although the model simulations presented here were subject to many caveats, they did reveal some important patterns. The model projected marked differences among sites making generalisations about plume-dispersal patterns in mid-ocean ridges difficult. Nevertheless, the models predicted large horizontal and vertical plume-dispersals above the thresholds adopted. Persistent plumes (temporal frequency >50%, i.e., 6 months out of 12 months) were projected to disperse an average linear distance of 10 to 20 km, cover an area of 17 to 150 km2, and extend more than 800 m in the water column. In fact, the model projected that plumes may disperse beyond the licensed mining areas, reach the flanks and summits of nearby topographic features, and extend into the bathypelagic, mesopelagic, and epipelagic environments. Modelled plume-dispersal overlaps with the predicted distribution of cold-water corals and with existing fishing activities. These potential impacts would be of particular concern in regions such as the Azores, where local populations are highly dependent on the sea for their livelihoods. The findings of this study are an important initial step towards understanding the nature and magnitude of deep-sea mining impacts in space and time.

Impacts of atmospheric stilling and climate warming on cyanobacterial blooms: An individual-based modelling approach

Harmful algal blooms of the freshwater cyanobacteria genus Microcystis are a global problem and are expected to intensify with climate change. In studies of climate change impacts on Microcystis blooms, atmospheric stilling has not been considered. Stilling is expected to occur in some regions of the world with climate warming, and it will affect lake stratification regimes. We tested if stilling could affect water column Microcystis distributions using a novel individual-based model (IBM). Using the IBM coupled to a three-dimensional hydrodynamic model, we assessed responses of colonial Microcystis biomass to wind speed decrease and air temperature increase projected under a future climate. The IBM altered Microcystis colony size using relationships with turbulence from the literature, and included light, temperature, and nutrient effects on Microcystis growth using input data from a shallow urban lake. The model results show that dynamic variations in colony size are critical for accurate prediction of cyanobacterial bloom development and decay. Colony size (mean and variability) increased more than six-fold for a 20% decrease in wind speed compared with a 2 °C increase in air temperature. Our results suggest that atmospheric stilling needs to be included in projections of changes in the frequency, distribution and magnitude of blooms of buoyant, colony-forming cyanobacteria under climate change.

Possible Implications of Relatively High Levels of Initial 60Fe in Iron Meteorites for the Noncarbonaceous–Carbonaceous Meteorite Dichotomy and Solar Nebula Formation

Cook et al. found that iron meteorites have an initial abundance ratio of the short-lived isotope 60Fe to the stable isotope 56Fe of 60Fe/56Fe ∼ (6.4 ± 2.0) × 10−7. This appears to require the injection of live 60Fe from a Type II supernova (SN II) into the presolar molecular cloud core, as the observed ratio is over a factor of 10 times higher than would be expected to be found in the ambient interstellar medium (ISM) as a result of galactic chemical evolution. The supernova triggering and injection scenario offers a ready explanation for an elevated initial 60Fe level, and in addition provides a physical mechanism for explaining the noncarbonaceous–carbonaceous (NC–CC) dichotomy of meteorites. The NC–CC scenario hypothesizes the solar nebula first accreted material that was enriched in supernova-derived nuclides, and then later accreted material depleted in supernova-derived nuclides. While the NC–CC dichotomy refers to stable nuclides, not short-lived isotopes like 60Fe, the SN II triggering hypothesis provides an explanation for the otherwise unexplained change in nuclides being accreted by the solar nebula. Three-dimensional hydrodynamical models of SN II shock-triggered collapse show that after triggering collapse of the presolar cloud core, the shock front sweeps away the local ISM while accelerating the resulting protostar/disk to a speed of several kilometers per second, sufficient for the protostar/disk system to encounter within ∼1 Myr the more distant regions of a giant molecular cloud complex that might be expected to have a depleted inventory of supernova-derived nuclides.

CFD simulation and model predictive control of the pipeline transportation system in a trailing suction hopper dredger

Abstract The phenomenon of blocking pipe would appear for a trailing suction hopper dredger (TSHD) rake arm pipe in the construction process. Thus, based on the theory of particle flow mechanics, a transient three-dimensional two-phase hydrodynamic and sediment mixture model was established in this paper. Besides, different particle sizes of sediment concentration and velocity in the piping in the rake arm of the law were analyzed to solve the problem of blocking pipe and the prevention of dredging to provide theoretical support. Then, a model predictive controller was designed to regulate the mudflow in the pipeline by controlling the mud pump speed and compared with the proportional-integral (PI) controller. According to the results, the particle size of sediment affects concentration distribution in the pipeline under the same construction conditions for the TSHD. Besides, the larger the particle size of the sediment, the more significant the difference in the sediment concentration distribution in the pipeline. Similarly, the flow velocity is another influencing factor behind the change in concentration distribution in the pipeline, the increase in the flow velocity will improve the uniformity of concentration distribution in the pipeline, and the critical velocity will maintain the concentration distribution within a reasonable range. In terms of transportation control, the designed predictive model controller is capable of reducing overshoot and shorten the time required for adjustment. To sum up, the research result not only provides a valuable reference for the theoretical analysis and a controller design of the pipeline transportation system in the TSHD but also verifies the feasibility of the computational fluid dynamics model used in the study of pipeline transportation mechanisms.

Hydrodynamic Drivers of Nutrient and Phytoplankton Dynamics in a Subtropical Reservoir

Chlorophyll-a (chla) is an important parameter to assess water quality in lakes and reservoirs, since it is a proxy for phytoplankton biomass and primary production. The increasing availability of data with high spatial and temporal resolution allows assessing short-term dynamics and small-scale variations of chla within larger water bodies. In freshwater reservoirs, the nutrient concentrations and the physical conditions that control phytoplankton growth vary along their longitudinal extend. Here, we analyze how the flow paths of the inflowing river water into density-stratified reservoirs affect the vertical and longitudinal distribution of nutrients and chla. We combine spatially resolved and high-frequency measurements of chla from satellite remote sensing and in-situ sensors, with numerical simulations using a three-dimensional hydrodynamic model to assess the influence of density currents on chla dynamics along a subtropical drinking water reservoir in the south of Brazil. Chla did not have pronounced seasonal dynamics (4.9 ± 1.2 µg L−1, at the continuous measurement station); instead, spatial variability along the reservoir was more pronounced (4.4 ± 2.1 µg L−1, all monitored points within the euphotic zone). Most of the nutrients from the inflowing river were consumed in the upstream region, and phytoplankton in the lacustrine zone depended on internal loading. Temporal variability was observed only in the upstream region, and phytoplankton produced in that area was transported downstream by density currents, resulting in large concentrations of chla below the euphotic zone. The results of a hydrodynamic model that simulates the present state are in good agreement with the observations. Two simulated scenarios, where the density current patterns were altered, resulted only in slight variations in density currents, indicating that the influence of the main inflow was of minor relevance in chla concentrations in downstream regions of the reservoir. Our results highlight the importance of two-dimensional hydrodynamic processes in regulating phytoplankton dynamics in reservoirs.

Potential Impact of Late 1980s Regime Shift on the Collapse of Walleye Pollock Catch in the Western East/Japan Sea

Walleye pollock (Gadus chalcogramma) caught in the Korean fishing area dramatically decreased in the late 1980s. To investigate the potential impact of the late 1980s climate regime shift on the collapse of the pollock catch, we developed a three-dimensional hydrodynamic model with data assimilation and a particle tracking model. Data-assimilated reanalysis showed that sea surface temperature increased by approximately 2°C in the spawning area of pollock in the late 1980s. The suitable spawning area in the East Korean Bay decreased due to warming in the late 1980s. Spawned eggs of walleye pollock were tracked using a particle tracking model for 30 days in January and February during 1983–1992. The number of individuals transported to the nursery within the Korean fishing area from the spawning area was reduced by 74% in the late 1980s. The intensified East Korean Warm Current (EKWC) could be responsible for the decreased number of individuals transported to the southern area in the late 1980s. Warming in the Korean fishing area could also cause a decrease in pollock. These oceanic changes might be linked to climate regime shifts in the late 1980s. The warming regime with positive Arctic Oscillation and weakened monsoon intensified the northward flow of the EKWC and accelerated the warming of the spawning and fishing areas in the late 1980s.

The calculation of flushing time for the upper Pasur River Estuary, Bangladesh

Estimation of estuarine flushing time, a time required to transport of pollutants or any other properties from estuaries to the coast, is very important for its resource management. In this study, we estimated flushing time (T) of the upper Pasur River Estuary (UPRE) for understanding the water quality condition in the dry and wet seasons. High-resolution salinity data were collected from the PRE at high water in the dry and wet seasons in 2014 and 2019. Flushing time was calculated using the freshwater fraction method (FFM) as well as e-folding flushing time scales was estimated empirically using the salinity (non-reactive conservative) and monthly river discharge data. System flushing during the dry season was thirteen times weaker than flushing during the wet season owing to decreasing river discharge by nearly 94%. In addition, the daily exchange volume was decreased by eight times during the dry season than during the wet season. As a consequence, the conditions of the UPRE are more dynamic during the wet season due to receiving huge amount of river discharge. During the dry season, only the e-folding time scales showed higher values in the salinity maximum zone (salt plug area). This implied that the e-folding time scale is an empirical approach and was able to encompass the tidal dispersion process whereas the FFM was unable to include that process. As the PRE is a macrotidal estuary, the tide assists to flush dissolved substances from the UPRE to the coast during the dry season having negligible river discharge. In addition, there was no significant variation in water quality parameter between the salt plug area (SP) and downstream of salt plug area (DSP) during the dry season. In order to have more accuracy, a three-dimensional hydrodynamic model would be useful to compute estuarine time scales precisely.

Effects of flood discharge on the water quality of a drinking water reservoir in China – Characteristics and management strategies

Stormflow runoff is an important non-point source of pollution in drinking water reservoirs. Rationally managed flood discharge processes at estuaries can reduce the high concentration of pollutants carried by runoff, and thus their impacts on water plant operations. In this study, the physical and chemical water quality parameters upstream from a dam were measured for a flood discharge process in the Jinpen Reservoir of Northwest China. The results showed that the time needed for the flood to reach the reservoir was ∼9 h after flooding began, and this lag effect meant that water in the metalimnion and hypolimnion were vented in advance. Consequently, the undercurrent intruded into the hypolimnion. The water temperature increased by 1.83 °C, and the Schmidt stability index decreased from 3291.37 J m−2 to 2496.32 J m−2. Flood discharge can effectively reduce the volume of pollutants; however, it cannot completely prevent the deterioration of water quality in the main reservoir. The turbidity of the reservoir still exceeded 300 NTU after the flood discharge. When the outflow discharge decreased, the critical height of aspiration also decreased, and dissolved pollutants could neither be vented nor precipitated quickly, not allowing their concentration within the reservoir to decrease further. A three-dimensional hydrodynamic model was successfully used to simulate flood processes and determine an optimal flood discharge plan. We found that the elevations of the outlet and undercurrent layer are essential to determine the efficiency of flood discharge, while appropriate timing and outflow discharge volumes may further enhance the results. Therefore, using a hydrodynamic model to predict the position of an undercurrent layer, and opening spillways near the elevation of the undercurrent layer can improve the efficacy of flood discharge. Our findings provide novel insights that may be used to improve the operation and management of source water reservoirs.

Assessing the potential impacts of the Canal Istanbul on the physical oceanography of the Turkish Straits System

The Turkish Straits System (TSS) is a region of intense mixing between the Black Sea and the Mediterranean Sea exerting crucial control on the water & substance exchange through the Bosphorus and Dardanelles straits between these two sea basins. A waterway project, namely the Canal Istanbul, was proposed to be built within this system parallel to the Bosphorus. Even tough environmental impacts of the project have been of concern, a few studies that address these scientific and public concerns under the light of scientific methods are available. In order to fill this gap in the literature, a three-dimensional hydrodynamic model with realistic forcing has been set up for long-term integration. The study includes simulating of two scenarios for the adjacent basins; (i) existing configuration of the TSS, and (ii) the proposed extension with the recent planned route of Canal Istanbul that was added to the validated model. The results reveals that the canal project would increase the exchange flow between the Black Sea and the Mediterranean Sea in both directions while causing slight changes in temperature & salinity along the TSS. Findings of this study can be used in further research that would analyze the impacts of the proposed Canal Istanbul on the water quality and ecosystem of the TSS.

Modeling of streamflow in a 30 km long reach spanning 5 years using OpenFOAM 5.x

Abstract. Developing accurate and efficient modeling techniques for streamflow at the tens-of-kilometers spatial scale and multi-year temporal scale is critical for evaluating and predicting the impact of climate- and human-induced discharge variations on river hydrodynamics. However, achieving such a goal is challenging because of limited surveys of streambed hydraulic roughness, uncertain boundary condition specifications, and high computational costs. We demonstrate that accurate and efficient three-dimensional (3-D) hydrodynamic modeling of natural rivers at 30 km and 5-year scales is feasible using the following three techniques within OpenFOAM, an open-source computational fluid dynamics platform: (1) generating a distributed hydraulic roughness field for the streambed by integrating water-stage observation data, a rough wall theory, and a local roughness optimization and adjustment strategy; (2) prescribing the boundary condition for the inflow and outflow by integrating precomputed results of a one-dimensional (1-D) hydraulic model with the 3-D model; and (3) reducing computational time using multiple parallel runs constrained by 1-D inflow and outflow boundary conditions. Streamflow modeling for a 30 km long reach in the Columbia River (CR) over 58 months can be achieved in less than 6 d using 1.1 million CPU hours. The mean error between the modeled and the observed water stages for our simulated CR reach ranges from −16 to 9 cm (equivalent to approximately ±7 % relative to the average water depth) at seven locations during most of the years between 2011 and 2019. We can reproduce the velocity distribution measured by the acoustic Doppler current profiler (ADCP). The correlation coefficients of the depth-averaged velocity between the model and ADCP measurements are in the range between 0.71 and 0.83 at 75 % of the survey cross sections. With the validated model, we further show that the relative importance of dynamic pressure versus hydrostatic pressure varies with discharge variations and topography heterogeneity. Given the model's high accuracy and computational efficiency, the model framework provides a generic approach to evaluate and predict the impacts of climate- and human-induced discharge variations on river hydrodynamics at tens-of-kilometers and decadal scales.

Swimming behavior of emigrating Chinook Salmon smolts

Swimming behavior of Chinook Salmon (Oncorhynchus tshawytscha) smolts affects transit time, route selection and survival in complex aquatic ecosystems. Behavior quantified at the river reach and junction scale is of particular importance for route selection and predator avoidance, though few studies have developed field-based approaches for quantifying swimming behavior of juvenile migratory fishes at this fine spatial scale. Two-dimensional acoustic fish telemetry at a river junction was combined with a three-dimensional hydrodynamic model to estimate in situ emigration swimming behavior of federally-threatened juvenile Chinook salmon smolts. Fish velocity over ground was estimated from telemetry, while the hydrodynamic model supplied simultaneous, colocated water velocities, with swimming velocity defined by the vector difference of the two velocities. Resulting swimming speeds were centered around 2 body lengths/second, and included distinct behaviors of positive rheotaxis, negative rheotaxis, lateral swimming, and passive transport. Lateral movement increased during the day, and positive rheotaxis increased in response to local hydrodynamic velocities. Swim velocity estimates were sensitive to the combination of vertical shear in water velocities and vertical distribution of fish.

A Numerical Tool for Dissolved Oxygen Simulation in the Western Mississippi Sound

Bazgirkhoob, H.; Linhoss, A., and Armandei, M., 2022. A numerical tool for dissolved oxygen simulation in the western Mississippi Sound. Journal of Coastal Research, 38(4), 699–711. Coconut Creek (Florida), ISSN 0749-0208. Dissolved oxygen is an indicator of aquatic health. Hydrodynamic–water quality models facilitate the understanding of flow dynamics and dissolved oxygen in estuarine systems. In this study, dissolved oxygen concentration, seasonality, and vertical stratification were studied in the western Mississippi Sound, Northern Gulf of Mexico. Visual Environmental Fluid Dynamic Code and the water quality model CE-QUAL-ICM were linked together to develop a three-dimensional hydrodynamic and water quality model of the western Mississippi Sound. The model was verified and validated for nutrients and chlorophyl a. Then, the modeled results for dissolved oxygen were compared to measured data. Mean observed dissolved oxygen ranged between 6.29 and 8.48 mg/L in the top and bottom layers at each site. The root mean squared error between measured and modeled dissolved oxygen ranged between 1.30 and 2.87 mg/L. Results indicate that the model was able to capture measured seasonal trends and vertical stratification. The model simulated more realistic results in shallower areas compared to deeper areas. Overall, the model simulated slightly higher dissolved oxygen than what was measured, although 27% of the stations showed no statistical difference between the modeled and measured dissolved oxygen. These results are useful for interpreting the value of this model for future ecological studies. The model can be used as a tool to provide better understanding of the dissolved oxygen dynamics due to multiple scenarios such as hurricanes and storms.

The Impacts of Freshwater Input and Surface Wind Velocity on the Strength and Extent of a Large High Latitude River Plume

Arctic Ocean physical and biogeochemical properties are strongly influenced by freshwater input from land and through the Bering Strait, where the mean currents transport water northward from the Bering Sea. The Yukon River is one of the largest rivers in North America and the Arctic, contributing large quantities of freshwater and terrigenous material to the coastal ocean in the northern Bering Sea. However, a detailed analysis of the coastal hydrodynamics at the outflow of the river has not been conducted in this remote but regionally important river. A three-dimensional hydrodynamic model was built to represent the lower Yukon River and coastal ocean for the ice-free months in 7 years. On average, a large anticyclonic eddy persisted at the main outflow of the Yukon that recirculates water back toward the coast where the currents converge to form a mean northward transport along the delta. Interannual spatial variance in salinity was relatively small, while there was substantial variance in u and v current velocity. u velocity spatial variance was correlated to the volume of freshwater discharge across years, while v velocity spatial variance was correlated to the N–S wind velocity. During strong wind events, plume structure was substantially altered: southerly winds deepened the plume and enhanced northward transport, while northerly winds shoaled and strengthened the pycnocline, and reversed the flow toward the south. The variability in plume dispersion on short time scales due to wind forcing has implications for where terrigenous material is processed in and settles out of the water column.

A modelling framework to track phosphorus sources of the drinking water intakes in a large eutrophic lake

Excess phosphorus (P) in freshwater lakes is a global challenge due to its potential cause of harmful algal blooms threating drinking water safety. However, quantifying the P sources for a specific site in a large lake is extremely challenging due to the complex interaction between internal and external P loading to surface water. To address this challenge, this study developed a modelling framework to track P sources of two drinking water intakes in a large shallow lake (Lake Taihu) in China. The framework proposed a new-developed index (PSCI, P source contribution index) to quantify the contributions of various P sources for the drinking water intakes. PSCI was derived by a three-dimensional hydrodynamic and water quality model that describes the complex processes of P sources, sinks and transportation in both horizontal and vertical directions within the lake. Application of the framework to these two drinking water intakes (Jinshu and Shangshan) achieved a surprising finding that internal P loading from sediment was a significant P source with a contribution as high as 47.1% and 30.4%, respectively. Central Lake Taihu had a large contribution of 49.5% and 68.3% for Jinshu and Shangshan, and inflow river (Wangyu River) had a contribution of < 5%. The high contribution of internal P loading was mainly due to the sediment P accumulation for years and large external P loading. Our study highlighted the important role of internal P loading affecting the P concentration of drinking water intake in a large eutrophic lake, and demonstrated the high value of the modelling framework in quantifying the P sources for a specific site (e.g., drinking water intake) in a lake.

The GAPS Programme at TNG

Context. Because of its proximity to an active K-type star, the hot Jupiter WASP-80b has been identified as a possible excellent target for detecting and measuring He I absorption in the upper atmosphere. Aims. Our aim was to look for, and eventually measure and model, metastable He I atmospheric absorption. Methods. We observed four primary transits of WASP-80b in the optical and near-infrared using the HARPS-N and GIANO-B high-resolution spectrographs attached to the Telescopio Nazionale Galileo telescope, focusing the analysis on the He I triplet. We further employed a three-dimensional hydrodynamic aeronomy model to understand the observational results. Results. We did not find any signature of planetary absorption at the position of the He I triplet with an upper limit of 0.7% (i.e. 1.11 planetary radii; 95% confidence level). We re-estimated the high-energy stellar emission, which we combined with a stellar photospheric model, to generate the input for the hydrodynamic modelling. We determined that, assuming a solar He to H abundance ratio, He I absorption should have been detected. Considering a stellar wind 25 times weaker than solar, we could reproduce the non-detection only by assuming a He to H abundance ratio about 16 times smaller than solar. Instead, considering a stellar wind ten times stronger than solar, we could reproduce the non-detection only with a He to H abundance ratio about ten times smaller than solar. We attempted to understand this result by collecting all past He I measurements and looking for correlations with high-energy stellar emission and planetary gravity, but without success. Conclusions. WASP-80b is not the only planet with an estimated sub-solar He to H abundance ratio, which suggests the presence of efficient physical mechanisms (e.g. phase separation, magnetic fields) capable of significantly modifying the He to H content in the upper atmosphere of hot Jupiters. The planetary macroscopic properties and the shape of the stellar spectral energy distribution are not sufficient for predicting the presence or absence of detectable metastable He in a planetary atmosphere, since the He abundance also appears to play a major role.

Coastal upwelling by wind-driven forcing in the Caspian Sea: A numerical analysis

In this study, wind-driven coastal upwelling in the Caspian Sea was investigated using a developed three-dimensional hydrodynamic numerical model based on the Princeton Ocean Model (POM). The model was forced with wind fields and atmospheric fluxes from the ECMWF database and it considers freshwater inflows from the Volga, Kura and Ural Rivers. This model was implemented for 10 years (2008–2018). Findings indicated that the upwelling in the Caspian Sea was due to effects of wind and bottom topography, often occurring from May to September. In June and July, in the eastern part of the middle and sometimes southern basins, up to 3°C water temperature difference occurs between coastal and offshore areas. The vertical temperature gradient in the middle basin was larger than that in the southern basin. Upwelling in August in the eastern coasts of the middle basin within 25 km of coast from the depth of 15 m to the surface was shown, which was due to the effects of wind and bottom topography. In the middle basin, the highest vertical velocities caused by upwelling in June, July and August were 12, 13.82, and 10.36 m/day, respectively.

Three-dimensional hydrodynamic model of the Rance estuary (France) influenced by the world’s second largest tidal power plant

The Rance estuary is a small steep-sided ria, located in the Brittany coast of northern France, with a maximum spring tidal range of 13.5m and an average river discharge of 7m3/s. Taking advantage of this significant tidal range, the Rance tidal power station (RTPS) was built in the 1960s as the world’s first and largest tidal power plant, with peak output capacity of 240 Megawatts. It is currently the second world’s largest tidal power installation after the Sihwa-Lake tidal barrage. The RTPS has two active parts: a barrage of 6 sluice gates and a structure of 24 turbines. Despite a well-known effect of the plant on damping estuarine water levels, little attention has been given to currents vertical distribution and the plant's impact on the dynamics of freshwater-saltwater interface. Therefore, a three-dimensional model of the Rance estuary was developed. Moreover, currents and salinity measurements were carried out to validate the numerical model. Simulated and measured currents showed that (i)the RTPS induces an acceleration of flood currents directly upstream of the sluice gates and (ii)ebb currents are strengthened by the narrowness of the Saint-Hubert-Port. Finally, salinity analyses assessed the dynamics of the freshwater-saltwater interface which is pushed further upstream during summer.

Evaluation of the Impact of Meteorological Factors on the Stratification of Structure in Lake Biwa, Japan

Hypoxia in Lake Biwa, Japan remains a serious water environmental problem. One of the causes of hypoxia in the lake is the formation of a thermocline, which is largely affected by meteorological factors, such as (1) air temperature, (2) wind speed, and (3) precipitation. However, the effects of these three meteorological factors on the formation of the thermocline have not been clarified quantitatively. In this study, applying a three-dimensional hydrodynamic model to Lake Biwa, the effects of each of the three meteorological elements on the formation of the thermocline was quantitatively analyzed to clarify the governing factors of meteorological conditions in the formation of anoxic oxygen. Sensitivity analysis of the stratification structure in Lake Biwa was performed by changing the three meteorological factors of (1) air temperature, (2) wind speed, and (3) precipitation. As a result, the change in wind speed gives the greatest effect on the stratification structure, the change in air temperature makes the difference in the stratification structure from the surface layer to the vicinity of the thermocline, and the change in precipitation affects it less than the others.

Multi-column modelling of lake Geneva for climate applications

The interaction between large inland water bodies and the atmosphere impacts the evolution of regional weather and climate, which in turn affects the lake dynamics, thermodynamics, ice-formation, and, therefore, ecosystems. Over the last decades, various approaches have been used to model lake thermodynamics and dynamics in standalone mode or coupled to numerical atmospheric models. We assess a turbulence-closure k-epsilon multi-column lake model in standalone mode as a computationally-efficient alternative to a full three-dimensional hydrodynamic model in the case of Lake Geneva. While it struggles to reproduce some short-term features, the multi-column model reasonably reproduces the seasonal mean of the thermal horizontal and vertical structures governing heat and mass exchanges between the lake surface and the lower atmosphere (stratified period, thermocline depth, stability of the water column). As it requires typically two orders of magnitude less computational ressources, it may allow a two-way coupling with a RCM on timescales or spatial resolutions where full 3D lake models are too demanding.