一个三维斜压水动力模型的建立及其在日本琵琶湖中的应用

水温的空间分布对很多湖泊水文现象关系甚大,本文建立了日本琵琶湖的三维斜压水动力数值模式。模拟计算结果表明：１）在不同风速驱动下,湖的上、下两层可形成不同的环流型;２）除形成环流处外,其它部分流速、流向在上、下层表现出明显的垂直切变;３）闭合式温度场储存位能,即使在无风和初始静止流场情况下,也可驱动湖流并形成若干环流特征;４）风对水的驱动作用是形成北湖北部稳定环流的决定性因子,并且本身可构造出温度场     

山体遮挡对滇池风生流的影响初探

用二维风生流数值模型模拟滇池湖流运动。滇池在主导风向西南风作用下,假定湖面风场是均匀的,数值模拟的湖流流态与实测湖流结果相差很大。而考虑山体遮挡影响,根据实测湖流期间现有的风情资料,在湖面上构造一非均匀风场,数值模拟结果与实测值基本一致。山体遮挡对滇池风生流的影响是不容忽视的。建议进一步进行湖流和湖面风向、风速监测,并建立过山气流数学模型,深入研究山体遮挡对湖泊风生流的影响。     

Climate warming alters thermal stability but not stratification phenology in a small north‐temperate lake

Recent climate change represents one of the most serious anthropogenic threats to lake ecosystems in Canada. As meteorological and hydrological conditions are altered by climate change, so too are physical, chemical and biological properties of lakes. The ability to quantify the impact of climate change on the physical properties of lakes represents an integral step in estimating future chemical and biological change. To that end, we have used the dynamic reservoir simulation model, a one‐dimensional vertical heat transfer and mixing model, to hindcast and compare lake temperature‐depth profiles against 30 years of long‐term monitoring data in Harp Lake, Ontario. These temperature profiles were used to calculate annual (June–September) thermal stability values from 1979 to 2009. Comparisons between measured and modelled lake water temperature and thermal stability over three decades showed strong correlation (r² > 0.9). However, despite significant increases in modelled thermal stability over the 30 year record, we found no significant change in the timing of the onset, breakdown or the duration of thermal stratification. Our data suggest that increased air temperature and decreased wind are the primary drivers of enhanced stability in Harp Lake since 1979. The high‐predictive ability of the Harp Lake dynamic reservoir simulation model suggests that its use as a tool in future lake management projects is appropriate. Copyright © 2013 John Wiley & Sons, Ltd.     

Field isotope investigations and numerical simulation of lake hydrodynamics

A three‐dimensional, numerical free‐surface‐flow simulation model is developed to investigate hydrodynamics of a lake and the predictive capabilities of the numerical model are validated by comparing them with field results obtained from Lake Naini, Uttar Pradesh, using environmental isotope (δ¹⁸O, δD and ³H) techniques. This has been carried out in order to understand the hydrodynamics of a lake (i.e. circulations, mixing and transport of water within the lake). Copyright © 2002 John Wiley & Sons, Ltd.     

洪泽湖风生流对风场和水位变化响应特征

深入认识大型湖泊在不同风速、风向和水位下三维风生流结构特征对于湖泊污染控制、生态恢复及资源的开发利用具有重要意义.本文在构建笛卡尔坐标系下洪泽湖三维水动力模型的基础上,利用2次全湖30个点位流场监测数据验证了模型精度.基于1975-2020年长系列风场观测数据,确定了洪泽湖典型风速风向.在此基础上,模拟了16种不同风向,13种不同风速和20种不同水位工况条件下洪泽湖三维风生流结构.结果表明：水动力模型可以较好地刻画洪泽湖三维湖流变化特征.洪泽湖风生流结构随风向变化呈现出较大空间差异.风生流流速随着风速的升高呈加速上升趋势,其中表层水体流速上升幅度远高于其他水层.在2.4 m/s东风驱动下,溧河洼、成子湖和南部湖区垂向平均流速随着水位上升呈先升高后降低的趋势,3个湖区分别在12.7、12.4和12.2 m水位下流速达到最大值.     

Mathematische Modelle zur Simulation von Zirkulationsströmungen in einem See

The three-dimensional, finite-difference model of R.L. Street (Stanford University) has been adapted at the simulation of wind-induced currents in the northern basin of the Lake of Lugano. The model makes use of complete Navier-Stokes equations and a dynamical algorithm to calculate the eddy viscosity. The model predictions are compared with data collected in the measure campaign of 1979. A further application of the model is the cold inflow simulation in a stratified basin with inclusion of turbulent vertical exchange processes. A two-dimensional version of the model, which assumes an hydrostatic pressure distribution, is utilized to simulate the thermal sommer development of the Lake of Lugano, during 40 days. The model includes the thermical interactions and the wind influence. The results obtained show the possible application of the three-dimensional model to the currents simulation in lake basins. The two-dimensional model may be utilized to predict the development of a basin subject to a thermal load.      

Mixing processes in lakes: mechanisms and ecological relevance

In Lake Baldegg, Switzerland (surface area 5.3 km², maximum depth 66 m) the analysis of data from moored instrument systems (atmospheric boundary layer, lake temperature distribution, bottom currents) was correlated to the long-term development of vertical mixing as seen from profiles of natural isotopes (radon-222, tritium and helium-3) and chemical species. The investigation shows: 1. Vertical mixing coefficients below 25 m are small. Consequently the vertical concentration distribution of sediment emanating species in the deep hypolimnion is controlled by the bottom topography. 2. Renewal of deep hypolimnic water is significant even during stratification. 3. Weakly damped internal waves characterize the internal dynamics during stratification. 4. Horizontal bottom currents play an important role in the hypolimnion mixing and can be correlated to internal waves during stratification.     

Modelling mixing and circulation in subglacial Lake Vostok,Antarctica

Lake Vostok, isolated from direct exchange with the atmosphere by about 4 km of ice for millions of years, provides a unique environment. This inaccessibility raises the importance of numerical models to investigate the physical conditions within the lake. Using a three-dimensional numerical model and the best available geometry, we test different parameter settings to define a standard model configuration suitable for studying flow in this subglacial lake. From our model runs we find a baroclinic circulation within the lake that splits into three different parts: Along a topographic ridge in the northern part of Lake Vostok, bottom water masses are transported eastward, diverging away from the ridge. In the lake’s surface layer, the flow in these two vertical overturning cells has opposite directions. In the southern part of the lake, where freezing occurs across about 3,500 km², two opposing gyres split the water column vertically. The general flow is stronger in the southern basin with horizontal velocities in the order of 1 mm/s. The strongest upwelling, found in the eastern part of this basin, is about 25 μm/s. We estimate the lower limit of the overturning timescale to be about 2.5 years vertically and 8.6 years horizontally. The basal mass loss of ice from the ice sheet floating on the lake is 5.6 mm/year (equivalent to a fresh water flux of 2.78 m³/s, or a basal ice loss of 0.09 km³/year). This imbalance indicates either a constant growth of the lake or its continuous (or periodical) discharge into a subglacial drainage system.     

青藏高原拉昂错热力分层和混合层深度变化特征观测

湖泊热力结构不仅影响湖泊内部生态环境,而且与区域气象和气候系统相互影响,但目前对湖泊垂直温度的观测研究仍非常匮乏.本研究基于青藏高原拉昂错连续的湖温和气象观测,分析了小时尺度和日尺度热力分层规律和混合层深度的变化特征.结果表明:拉昂错为冷多次完全混合型湖泊;湖表温度8月达到最大值,湖面敞水区和沿岸的湖表温度季节震荡相同,但存在一定的空间差异;与空气温度相比,湖表温度变幅小,降温更慢,峰值滞后约1个月.同时发现,拉昂错混合层深度变化较为规律,辐射和风速是影响湖泊混合层深度的主要外界因子.探明真实的高原湖泊分层规律,有利于提高湖泊模拟能力,为优化湖泊参数化方案提供参考.     

Sensitivity of the multiannual thermal dynamics of a deep pre‐alpine lake to climatic change

The study of the multiannual thermal dynamics of Lake Iseo, a deep lake in the Italian pre‐alpine area, is presented. Interflow was found to be the dominant river entrance mode, suggesting future susceptibility of the lake thermal structure to the overall effects of climate change expected in the upstream alpine watershed. A lake model employed the results of a long‐term hydrologic model to simulate the effects of a climate change scenario on the lake's thermal evolution for the period 2012–2050. The model predicts an overall average increase in the lake water temperature of 0.012 °C/year and a reinforced Schmidt thermal stability of the water column in the winter up to 800 J/m². Both these effects may further hinder the deep circulation process, which is vital for the oxygenation of deep water. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical analysis of the Black Sea currents and mesoscale eddies in 2006 and 2011

Two prognostic experiments taking into account real atmospheric forcing for 2006 and 2011 were carried out based on the eddy-resolving numerical model with a horizontal resolution of 1.6 km for the Black Sea. The main dynamic features such as the Rim Current, the Sevastopol, and Batumi anticyclones are reproduced in both experiments. The model results are confirmed via observation data. We accomplished the analysis of simulated circulation and energetics. The results demonstrate that both the vertical viscosity and vertical diffusion along with the energy inflow from the wind have been the main contributors to the annual and seasonal budgets of kinetic and potential energies of the Black Sea circulation. It is shown that two regimes of the Black Sea general circulation are implemented depending on a magnitude of wind contribution to the kinetic energy in winter. Intensive mesoscale eddy formation was observed along the Anatolian, Caucasian, and Crimean coasts. The analysis of the Black Sea circulation and eddy energetics allowed us to conclude that the generation and development of the mesoscale coastal eddies is associated with the barotropic instability in case of intensive coastal currents and is associated with both the barotropic and baroclinic instability in case of weak coastal currents.     

The baroclinic response to wind in a small two-basin lake

Field data and two linear layered models were used to examine the baroclinic response to wind in a small elongated two-basin lake, Amisk Lake (Alberta, Canada). For the first vertical baroclinic mode, wind-forced horizontal modes were simulated using a dynamic two-layer variable cross-section (TVC) model. The first horizontal mode, H1, was found to dominate the exchange between the two basins of Amisk Lake. The highest velocities associated with H1 occurred in the constricted channel connecting the two basins. This high velocity led to strong damping which brought H1 in near-resonance with the diurnal wind despite a difference in periods. Along with H1, the second horizontal mode, H2, was detected at a thermistor mooring. The TVC showed that H2 resulted from the coupling between along-thalweg wind distribution and the bathymetry of Amisk Lake. The damping for H2 was found to be weaker than for H1, likely because of weaker bottom drag in the connecting channel. The evolution of vertical H1 modes were simulated using a multi-layered box model. In response to wind pulses, the first vertical mode V1H1 initially dominated over higher vertical modes causing two-layer exchange. Following the faster decay of V1H1, higher vertical modes shifted the exchange to three and more layers. Our study shows the importance of the coupling between wind stress distribution and lake bathymetry in exciting horizontal modes, reveals the effect of damping on resonance with wind, and explains the evolution of exchange associated with vertical modes.     

Upwelling and surface cold patches in the Yellow Sea in summer: Effects of tidal mixing on the vertical circulation

A three-dimensional, prognostic, wave–tide–circulation coupled numerical model is developed to study the effects of tidal mixing on the summertime vertical circulation in the Yellow Sea (YS). The distribution and mechanisms of upwelling are investigated by numerical means. Validated by historical tide gauge data, satellite sea surface temperature (SST) data, and cruise observation data, the model shows satisfactory performances in reproducing the dominant tidal system and three-dimensional sea temperature structure. Model results suggest that strong tidal mixing plays an important role in the formation of the vertical circulation in the YS. The Yellow Sea Cold Water Mass (YSCWM) is fringed by typical tidal mixing fronts (TMFs), which separate the cold, stratified water at the offshore side from the warm, well-mixed, shallow water at the other side. Considerable baroclinic gradient across the TMF makes the frontal zone the spot where the most active vertical circulation occurs; a secondary circulation is triggered with a distinct upwelling branch occurring mainly on the mixed side of the front. The numerical model produces systematic upwelling belts surrounding the YSCWM, and the upwelling is essentially induced by the TMF over sloping topography. The relative importance of tidal mixing and wind forcing for upwelling is further examined in numerical experiments. The southerly wind enhances the upwelling off the western coasts, but its overall influences in the whole YS are less important than tidal mixing. As shown by both satellite data and numerical modeling, the summertime SST field in the YS is featured by the stable existence of several site-selective surface cold patches (SCPs), most of which scatter in the waters off convex coastlines. One of the SCPs is found off Subei Bank, and the others are located off the eastern tip of Shandong Peninsula and off the three tips of Korean Peninsula. Two processes give rise to the SCP: on the one hand, TMF-induced upwelling supplies cold water from the deep layer; on the other hand, tidal mixing itself can stir the bottom water upward and homogenize the water column vertically. In the waters around the tips of peninsula in the YS, the tidal currents are extraordinarily strong, which provides a possible explanation for the site-selectivity of the SCPs.     

风场对藻类在太湖中迁移影响的动力学研究

在Ｗｅｂｓｔｅｒ等人工作的基础上,以三维湖流为背景,考虑了波浪及藻类自身浮力的影响,建立了一个太湖梅梁湾三维藻类迁多模型,以研究在不同几场作用下藻类在湖泊中的迁移过程。模拟结果表明,不同风场对于藻类在湖泊中的水平及垂直分布影响很大,并且存在着一临界风速,其范围在２－３ｍ／ｓ之间,当几速小于临床风速时,水面可以拟看作水动力学光滑,没有波浪产生,在水表面藻类顺着风向迅速的向迎风岸的边漂移,形成藻类大量     

Impacts of regional warming on long‐term hypolimnetic anoxia and dissolved oxygen concentration in a deep lake

Although few reports have described long‐term continuous anoxia in aquatic systems, Lake Ikeda in Japan experienced such conditions in the hypolimnion from 1990 to 2010. The present study aimed to assess temporal fluctuations in the lake's thermal stability from 1978 to 2011 to understand the influence of regional climate change on hypolimnetic anoxia in this lake. Because complete vertical mixing, which supplies dissolved oxygen (DO) to the hypolimnion, potentially occurs on February, we calculated the Schmidt stability index (S) in February and compared it with hypolimnetic DO dynamics. Vertical water temperature profiles were calculated using a one‐dimensional model, and calculated temperatures and meteorological data were used to analyse annual fluctuations in water temperatures, thermocline depth, meteorological variables and S. We estimated that mean annual air and volume‐weighted water temperatures increased by 0.028 and 0.033 °C year^?1, respectively, from 1978 to 2011. Between 1986 and 1990, S and water temperature increased abruptly, probably due to a large upwards trend in air temperature (+0.239 °C year^?1). We hypothesize that a mixing regime that lacked overturn took effect at this time and that this regime lasted until 2011, when S was particularly small. These results demonstrate that abrupt climate warming in the late 1980s likely triggered the termination of complete mixing and caused the 21‐year period of successive anoxia in Lake Ikeda. We conclude that the lake response to a rapid shift in regional climate conditions was a key factor in changing the hypolimnetic water environment and that thermal stability in winter is a critical environmental factor controlling the mixing regime and anoxic conditions in deep lakes. Copyright © 2014 John Wiley & Sons, Ltd.     

Circulation and thermal structure in Lake Huron and Georgian Bay: Application of a nested-grid hydrodynamic model

A nested-grid hydrodynamic modelling system is used to study circulation and temperature distributions in Lake Huron (LH) and adjacent areas. This nested system is based on the three-dimensional, primitive-equation z-level ocean model. The nested system consists of two sub-components: a coarse-resolution outer model covering LH and Georgian Bay (GB) with a horizontal resolution of roughly 2.5 km, and the fine-resolution inner model covering eastern LH and northwestern GB with a horizontal resolution of roughly 900 m. Both the outer and inner models have 30 z-levels in the vertical. To assess the model performance, we simulate the three-dimensional circulation and temperature distributions of LH and GB in 1974–1975 and compare the model results with observations made in the lake. We demonstrate that outer model of the nested system simulates reasonably well the large-scale circulation and seasonal evolution of thermal stratifications in LH and GB, and the inner model produces reasonably well the three-dimensional flow and thermal structure over the coastal boundary layer close to the eastern shore of the lake.     

The SOMS model and its application to Lake Neuchâtel

A three-dimensional numerical circulation model (SOMS) based on primitive equations is described. The algorithm, by which Coriolis and vertical diffusion terms are treated implicitly while mass is still conserved exactly (algebraically), is discussed in detail. The model is applied to Lake Neuchâtel (Switzerland), to determine the general circulation under influence of the most prevailing wind.     

Numerical study of three-dimensional shelf circulation on the Scotian Shelf using a shelf circulation model

A numerical shelf circulation model was developed for the Scotian Shelf, using a nested-grid setup consisting of a three-dimensional baroclinic inner model embedded inside a two-dimensional barotropic outer model. The shelf circulation model is based on the Princeton Ocean Model and driven by three-hourly atmospheric forcing provided by a numerical weather forecast model and by tidal forcing specified at the inner model's open boundaries based on pre-calculated tidal harmonic constants. The outer model simulates the depth-mean circulation forced by wind and atmospheric pressure fields over the northwest Atlantic Ocean with a horizontal resolution of 1/12°. The inner model simulates the three-dimensional circulation over the Gulf of St. Lawrence, the Scotian Shelf, and the adjacent slope with a horizontal resolution of 1/16°. The performance of the shelf circulation model is assessed by comparing model results with oceanographic observations made along the Atlantic coast of Nova Scotia and in the vicinity of Sable Island (on the Scotian Shelf) during two periods: October 2000–March 2001 and April–June 2002. Analysis of model results on Sable Island Bank indicates that tidal currents account for as much as ∼80% of the total variance of near-bottom currents, and currents driven by local winds account for ∼30% of the variance of the non-tidal near-bottom currents. Shelf waves generated remotely by winds and propagating into the region also play an important role in the near-bottom circulation on the bank.