Ice processes on medium‐sized north‐temperate lakes

Lake ice supports a range of socio‐economic and cultural activities including transportation and winter recreational actives. The influence of weather patterns on ice‐cover dynamics of temperate lakes requires further understanding for determining how changes in ice composition will impact ice safety and the range of ecosystem services provided by seasonal ice cover. An investigation of lake ice formation and decay for three lakes in Central Ontario, Canada, took place over the course of two winters, 2015–2016 and 2016–2017, through the use of outdoor digital cameras, a Shallow Water Ice Profiler (upward‐looking sonar), and weekly field measurements. Temperature fluctuations across 0°C promoted substantial early season white ice growth, with lesser amounts of black ice forming later in the season. Ice thickening processes observed were mainly through meltwater, or midwinter rain, refreezing on the ice surface. Snow redistribution was limited, with frequent melt events limiting the duration of fresh snow on the ice, leading to a fairly uniform distribution of white ice across the lakes in 2015–2016 (standard deviations week to week ranging from 3 to 5 cm), but with slightly more variability in 2016–2017 when more snow accumulated over the season (5 to 11 cm). White ice dominated the end‐of‐season ice composition for both seasons representing more than 70% of the total ice thickness, which is a stark contrast to Arctic lake ice that is composed mainly of black ice. This research has provided the first detailed lake ice processes and conditions from medium‐sized north‐temperate lakes and provided important information on temperate region lake ice characteristics that will enhance the understanding of the response of temperate lake ice to climate and provide insight on potential changes to more northern ice regimes under continued climate warming.     

Ice‐jam flood risk assessment and mapping

In northern regions, river ice‐ jam flooding can be more severe than open‐water flooding causing property and infrastructure damages, loss of human life and adverse impacts on aquatic ecosystems. Very little has been performed to assess the risk induced by ice‐related floods because most risk assessments are limited to open‐water floods. The specific objective of this study is to incorporate ice‐jam numerical modelling tools (e.g. RIVICE, Monte‐Carlo simulation) into flood hazard and risk assessment along the Peace River at the Town of Peace River (TPR) in Alberta, Canada. Adequate historical data for different ice‐jam and open‐water flooding events were available for this study site and were useful in developing ice‐affected stage‐frequency curves. These curves were then applied to calibrate a numerical hydraulic model, which simulated different ice jams and flood scenarios along the Peace River at the TPR. A Monte‐Carlo analysis was then carried out to acquire an ensemble of water level profiles to determine the 1 : 100‐year and 1 : 200‐year annual exceedance probability flood stages for the TPR. These flood stages were then used to map flood hazard and vulnerability of the TPR. Finally, the flood risk for a 200‐year return period was calculated to be an average of $32/m²/a ($/m²/a corresponds to a unit of annual expected damages or risk). Copyright © 2016 John Wiley & Sons, Ltd.     

Effects of changing climate on ice cover in three morphometrically different lakes

A one‐dimensional hydrodynamic lake model (DYRESM‐WQ‐I) is employed to simulate ice cover and water temperatures over the period 1911–2014. The effects of climate changes (air temperature and wind speed) on ice cover (ice‐on, ice‐off, ice cover duration, and maximum ice thickness) are modeled and compared for the three different morphometry lakes: Fish Lake, Lake Wingra, and Lake Mendota, located in Madison, Wisconsin, USA. It is found that the ice cover period has decreased due to later ice‐on dates and earlier ice‐off dates, and the annual maximum ice cover thickness has decreased for the three lakes during the last century. Based upon simulated perturbations of daily mean air temperatures across the range of ?10°C to +10°C of historical values, Fish Lake has the most occurrences of no ice cover and Lake Wingra still remains ice covered under extreme conditions (+10°C). Overall, shallower lakes with larger surface areas appear more resilient to ice cover changes caused by climate changes.     

Hydrogeomorphic processes of thermokarst lakes with grounded‐ice and floating‐ice regimes on the Arctic coastal plain,Alaska

Thermokarst lakes cover > 20% of the landscape throughout much of the Alaskan Arctic Coastal Plain (ACP) with shallow lakes freezing solid (grounded ice) and deeper lakes maintaining perennial liquid water (floating ice). Thus, lake depth relative to maximum ice thickness (1·5–2·0 m) represents an important threshold that impacts permafrost, aquatic habitat, and potentially geomorphic and hydrologic behaviour. We studied coupled hydrogeomorphic processes of 13 lakes representing a depth gradient across this threshold of maximum ice thickness by analysing remotely sensed, water quality, and climatic data over a 35‐year period. Shoreline erosion rates due to permafrost degradation ranged from < 0·2 m/year in very shallow lakes (0·4 m) up to 1·8 m/year in the deepest lakes (2·6 m). This pattern of thermokarst expansion masked detection of lake hydrologic change using remotely sensed imagery except for the shallowest lakes with stable shorelines. Changes in the surface area of these shallow lakes tracked interannual variation in precipitation minus evaporation (P ? E_L) with periods of full and nearly dry basins. Shorter‐term (2004–2008) specific conductance data indicated a drying pattern across lakes of all depths consistent with the long‐term record for only shallow lakes. Our analysis suggests that grounded‐ice lakes are ice‐free on average 37 days longer than floating‐ice lakes resulting in a longer period of evaporative loss and more frequent negative P ? E_L. These results suggest divergent hydrogeomorphic responses to a changing Arctic climate depending on the threshold created by water depth relative to maximum ice thickness in ACP lakes. Copyright © 2011 John Wiley & Sons, Ltd.     

Modelling inter-annual and spatial variability of ice cover in a temperate lake with complex morphology

The formation of ice cover on lakes alters heat and energy transfer with the water column. The fraction of surface area covered by ice and the timing of ice-on and ice-off therefore affects hydrodynamics and the seasonal development of stratification and related ecosystem processes. Multi-year model simulations of temperate lake ecosystems that freeze partially or completely therefore require simulation of the formation and duration of ice cover. Here we present a multi-year hydrodynamic simulation of an alpine lake with complex morphology (Lower Lake Constance, LLC) using the three-dimensional (3D) model Aquatic Ecosystem Model (AEM3D) over a period of 9 years. LLC is subdivided into three basins (Gnadensee, Zeller See and Rheinsee) which differ in depth, morphological features, hydrodynamic conditions and ice cover phenology and thickness. Model results were validated with field observations and additional information on ice cover derived from a citizen science approach using information from social media. The model reproduced the occurrence of thin ice as well as its inter-annual variability and differentiated the frequency and extent of ice cover between the three sub-basins. It captured that full ice cover occurs almost each winter in Gnadensee, but only rarely in Zeller See and Rheinsee. The results indicate that the 3D model AEM3D is suitable for simulating long-term dynamics of thin ice cover in lakes with complex morphology and inter-annual changes in spatially heterogeneous ice cover.     

Investigation of the snow‐cover dynamics in the Upper Euphrates Basin of Turkey using remotely sensed snow‐cover products and hydrometeorological data

The Euphrates and Tigris rivers serve as the most important water resources in the Middle East. Precipitation in this region falls mostly in the form of snow over the higher elevations of the Euphrates Basin and remains on the ground for nearly half of the year. This snow‐covered area (SCA) is a key element of the hydrological cycle, and monitoring the SCA is crucial for making accurate forecasts of snowmelt discharge, especially for energy production, flood control, irrigation, and reservoir‐operation optimization in the Upper Euphrates (Karasu) Basin. Remote sensing allows the detection of the spatio‐temporal patterns of snow cover across large areas in inaccessible terrain, such as the eastern part of Turkey, which is highly mountainous. In this study, a seasonal evaluation of the snow cover from 2000 to 2009 was performed using 8‐day snow‐cover products (MOD10C2) and the daily snow‐water equivalent (SWE) product. The values of SWE products were obtained using an assimilation process based on the Helsinki University of Technology model using equal area Special Sensor Microwave Imager (SSM/I) Earth‐gridded advanced microwave scanning radiometer—EOS daily brightness‐temperature values. In the Karasu Basin, the SCA percentage for the winter period is 80–90%. The relationship between the SCA and the runoff during the spring period is analysed for the period from 2004 to 2009. An inverse linear relationship between the normalized SCA and the normalized runoff values was obtained (r = 0·74). On the basis of the monthly mean temperature, total precipitation and snow depth observed at meteorological stations in the basin, the decrease in the peak discharges, and early occurrences of the peak discharges in 2008 and 2009 are due to the increase in the mean temperature and the decrease in the precipitation in April. Copyright © 2011 John Wiley & Sons, Ltd.     

Snow cover variability and snowmelt in a high‐altitude ungauged catchment

Snow variability is an integrated indicator of climate change, and it has important impacts on runoff regimes and water availability in high‐altitude catchments. Remote sensing techniques can make it possible to quantitatively detect the snow cover changes and associated hydrological effects in those poorly gauged regions. In this study, the spatial–temporal variations of snow cover and snow melting time in the Tuotuo River basin, which is the headwater of the Yangtze River, were evaluated based on satellite information from the Moderate Resolution Imaging Spectroradiometer snow cover product, and the snow melting equivalent and its contribution to the total runoff and baseflow were estimated by using degree–day model. The results showed that the snow cover percentage and the tendency of snow cover variability increased with rising altitude. From 2000 to 2012, warmer and wetter climate change resulted in an increase of the snow cover area. Since the 1960s, the start time for snow melt has become earlier by 0.9–3 days/10a and the end time of snow melt has become later by 0.6–2.3 days/10a. Under the control of snow cover and snow melting time, the equivalent of snow melting runoff in the Tuotuo River basin has been fluctuating. The average contributions of snowmelt to baseflow and total runoff were 19.6% and 6.8%, respectively. Findings from this study will serve as a reference for future research in areas where observational data are deficient and for planning of future water management strategies for the source region of the Yangtze River. Copyright © 2015 John Wiley & Sons, Ltd.     

Testing a model for predicting the timing and location of shallow landslide initiation in soil‐mantled landscapes

The growing availability of digital topographic data and the increased reliability of precipitation forecasts invite modelling efforts to predict the timing and location of shallow landslides in hilly and mountainous areas in order to reduce risk to an ever‐expanding human population. Here, we exploit a rare data set to develop and test such a model. In a 1·7 km² catchment a near‐annual aerial photographic coverage records just three single storm events over a 45 year period that produced multiple landslides. Such data enable us to test model performance by running the entire rainfall time series and determine whether just those three storms are correctly detected. To do this, we link a dynamic and spatially distributed shallow subsurface runoff model (similar to TOPMODEL) to an in?nite slope model to predict the spatial distribution of shallow landsliding. The spatial distribution of soil depth, a strong control on local landsliding, is predicted from a process‐based model. Because of its common availability, daily rainfall data were used to drive the model. Topographic data were derived from digitized 1 : 24 000 US Geological Survey contour maps. Analysis of the landslides shows that 97 occurred in 1955, 37 in 1982 and ?ve in 1998, although the heaviest rainfall was in 1982. Furthermore, intensity–duration analysis of available daily and hourly rainfall from the closest raingauges does not discriminate those three storms from others that did not generate failures. We explore the question of whether a mechanistic modelling approach is better able to identify landslide‐producing storms. Landslide and soil production parameters were ?xed from studies elsewhere. Four hydrologic parameters characterizing the saturated hydraulic conductivity of the soil and underlying bedrock and its decline with depth were ?rst calibrated on the 1955 landslide record. Success was characterized as the most number of actual landslides predicted with the least amount of total area predicted to be unstable. Because landslide area was consistently overpredicted, a threshold catchment area of predicted slope instability was used to de?ne whether a rainstorm was a signi?cant landslide producer. Many combinations of the four hydrological parameters performed equally well for the 1955 event, but only one combination successfully identi?ed the 1982 storm as the only landslide‐producing storm during the period 1980–86. Application of this parameter combination to the entire 45 year record successfully identi?ed the three events, but also predicted that two other landslide‐producing events should have occurred. This performance is signi?cantly better than the empirical intensity–duration threshold approach, but requires considerable calibration effort. Overprediction of instability, both for storms that produced landslides and for non‐producing storms, appears to arise from at least four causes: (1) coarse rainfall data time scale and inability to document short rainfall bursts and predict pressure wave response; (2) absence of local rainfall data; (3) legacy effect of previous landslides; and (4) inaccurate topographic and soil property data. Greater resolution of spatial and rainfall data, as well as topographic data, coupled with systematic documentation of landslides to create time series to test models, should lead to signi?cant improvements in shallow landslides forecasting. Copyright © 2003 John Wiley & Sons, Ltd.     

鲫(Carassius auratus)对浅水湖泊水体浊度、营养水平和浮游生物生物量的影响——基于中宇宙模拟实验

鲫(Carassius auratus)是我国各类淡水水体的优势鱼类之一.作为底栖杂食性鱼类,一方面,鲫可以通过排泄和扰动沉积物影响湖泊营养和光照水平,通过"上行效应"促进浮游植物生长;另一方面,鲫也可以捕食浮游动物,通过"下行控制"影响藻类生长以及营养盐循环.对于浅水湖泊,两种途径对于生态系统影响的相对重要性仍有待研究.本研究设计了一个两因素户外中宇宙实验,通过在沉积物表面添加隔网的方式,比较两种情况下(能、否接触沉积物),鲫对水体浊度、营养盐和浮游生物生物量的影响.实验在16个大型钢化玻璃桶(400 L)中进行,持续36 d(2019年8月6日—9月11日).研究结果表明:1)在能接触沉积物的条件下,鲫显著促进了沉积物再悬浮,表现为水体的总悬浮物(TSS)和无机悬浮物(ISS)浓度大幅升高;水体的光衰减系数(Kd)增加,总氮(TN)和总磷(TP)浓度明显升高; 2)在不能接触沉积物的条件下,鲫对水体悬浮物(TSS和ISS)浓度和Kd的影响不明显,但是显著降低了水体TN和TP浓度; 3)在两种情况下,鲫对浮游植物叶绿素a浓度以及浮游动物生物量的影响均不显著.本研究表明鲫只有在能够接触沉积物的条件下,才会显著提高水体浊度和营养水平.因此,在缺乏沉水植被的浅水湖泊中,鲫扰动沉积物产生的"上行效应"可能是其对生态系统产生负面影响的主要途径.     

我国大、中草型浅水型湖泊渔业综合开发技术进展

长江中下游湖大多属浅水草型湖泊,长期以来渔业利用处于薄后状态,1986－1990年国家在江苏鬲湖,湖北保定湖和安徽花园湖等地组织了大,中型浅水草型的湖泊的综合开发技术,包括水大面鱼病害控制,围养区层鱼捕捞技术,水体农业模式在内的科技攻关,本文以这五年取得的成果为主,概述了我国大,中型浅水草型湖泊渔业综合开发技术方面的进展。     

Trapping of aeolian sediments and build‐up of the ice cover of a dry‐based Antarctic lake

Among the perennially frozen lakes of the Dry Valleys of South Victoria Land (Antarctica), some are dry‐based, i.e. frozen to the bottom. One of these is studied here by a multiparametric investigation (isotopic composition in δD and δ¹⁸O, ions, gas and ice texture analyses). A sediment layer about 10 cm thick appearing at a depth of 3·5 m is also studied by grain size, X‐ray diffraction and scanning electron microscope analyses. The information retrieved indicates that this ice‐block lake results from a build‐up in two steps and explains how aeolian sediments were included as a layer into the ice. Copyright © 2002 John Wiley & Sons, Ltd.     

Evaluating the potential of a blended passive microwave‐interactive multi‐sensor product for improved mapping of snow cover and estimations of snow water equivalent

The retrieval of Snow Water Equivalent (SWE) from remote sensing satellites continues to be a very challenging problem. In this paper, we evaluate the accuracy of a new SWE product derived from the blending of a passive microwave SWE product based on the Advanced Microwave Sounding Unit (AMSU) with a multi‐sensor snow cover extent product based on the Interactive Multi‐sensor Snow and Ice Mapping System (IMS). The microwave measurements have the ability to penetrate the snow pack, and thus, the retrieval of SWE is best accomplished using the AMSU. On the other hand, the IMS maps snow cover more reliably due to the use of multiple satellite and ground observations. The evolution of global snow cover from the blended, the AMSU and the IMS products was examined during the 2006 snow season. Despite the overall good inter‐product agreement, it was shown that the retrievals of snow cover extent in the blended product are improved when using IMS, with implications for improved microwave retrievals of SWE. In a separate investigation, the skill of the microwave SWE product was also examined for its ability to correctly estimate SWE globally and regionally. Qualitative evaluation of global SWE retrievals suggested dependence on land surface temperature: the lower the temperature, the higher the SWE retrieved. This temperature bias was attributed in part to temperature effects on those snow properties that impact microwave response. Therefore, algorithm modifications are needed with more dynamical adjustments to account for changing snow cover. Quantitative evaluation over Slovakia in central Europe, for a limited period in 2006, showed reasonably good performance for SWE less than 100 mm. Sensitivity to deeper snow decreased significantly. Copyright © 2007 John Wiley & Sons, Ltd.     

Ensemble modelling of ice cover for a reservoir affected by pumped‐storage operation and climate change

Ensemble modelling was used to assess the robustness of projected impacts of pumped‐storage (PS) operation and climate change on reservoir ice cover. To this end, three one‐dimensional and a two‐dimensional laterally averaged hydrodynamic model were set up. For the latter, the strength of the impacts with increasing distance from the dam was also investigated. Climate change effects were simulated by forcing the models with 150 years of synthetic meteorological time series created with a weather generator based on available air temperature scenarios for Switzerland. Future climate by the end of the 21^st century was projected to shorten the ice‐covered period by ~2 months and decrease ice thicknesses by ~13 cm. Under current climate conditions, the ice cover would already be affected by extended PS operation. For example, the average probability of ice coverage on a specific day was projected to decrease by ~13% for current climate and could further be reduced from ~45% to ~10% for future climate. Overall, the results of all models were consistent. Although the number of winters without ice cover was projected to increase for all one‐dimensional models, studying individual segments of the two‐dimensional model showed that the impact was pronounced for segments close to the PS intake/outlet. In summary, the reservoir's ice cover is expected to partially vanish with higher probability of open water conditions closer to the PS intake/outlet.     

Estimation of winter leaf area index and sky view fraction for snow modelling in boreal coniferous forests: consequences on snow mass and energy balance

Leaf area index (LAI) and canopy coverage are important parameters when modelling snow process in coniferous forests, controlling interception and transmitting radiation. Estimates of LAI and sky view factor show large variability depending on the estimation method used, and it is not clear how this is reflected in the calculated snow processes beneath the canopy. In this study, the winter LAI and sky view fraction were estimated using different optical and biomass‐based approximations in several boreal coniferous forest stands in Fennoscandia with different stand density, age and site latitude. The biomass‐based estimate of LAI derived from forest inventory data was close to the values derived from the optical measurements at most sites, suggesting that forest inventory data can be used as input to snow hydrological modelling. Heterogeneity of tree species and site fertility, as well as edge effects between different forest compartments, caused differences in the LAI estimates at some sites. A snow energy and mass balance model (SNOWPACK) was applied to detect how the differences in the estimated values of the winter LAI and sky view fraction were reflected in simulated snow processes. In the simulations, an increase in LAI and a decrease in sky view fraction changed the snow surface energy balance by decreasing shortwave radiation input and increasing longwave radiation input. Changes in LAI and sky view fraction affected directly snow accumulation through altered throughfall fraction and indirectly snowmelt through the changed surface energy balance. Changes in LAI and sky view fraction had a greater impact on mean incoming radiation beneath the canopy than on other energy fluxes. Snowmelt was affected more than snow accumulation. The effect of canopy parameters on evaporation loss from intercepted snow was comparable with the effect of variation in governing meteorological variables such as precipitation intensity and air temperature. Copyright © 2012 John Wiley & Sons, Ltd.     

Snow model sensitivity analysis to understand spatial and temporal snow dynamics in a high‐elevation catchment

In this paper, we addressed a sensitivity analysis of the snow module of the GEOtop2.0 model at point and catchment scale in a small high‐elevation catchment in the Eastern Italian Alps (catchment size: 61 km²). Simulated snow depth and snow water equivalent at the point scale were compared with measured data at four locations from 2009 to 2013. At the catchment scale, simulated snow‐covered area (SCA) was compared with binary snow cover maps derived from moderate‐resolution imaging spectroradiometer (MODIS) and Landsat satellite imagery. Sensitivity analyses were used to assess the effect of different model parameterizations on model performance at both scales and the effect of different thresholds of simulated snow depth on the agreement with MODIS data. Our results at point scale indicated that modifying only the “snow correction factor” resulted in substantial improvements of the snow model and effectively compensated inaccurate winter precipitation by enhancing snow accumulation. SCA inaccuracies at catchment scale during accumulation and melt period were affected little by different snow depth thresholds when using calibrated winter precipitation from point scale. However, inaccuracies were strongly controlled by topographic characteristics and model parameterizations driving snow albedo (“snow ageing coefficient” and “extinction of snow albedo”) during accumulation and melt period. Although highest accuracies (overall accuracy = 1 in 86% of the catchment area) were observed during winter, lower accuracies (overall accuracy < 0.7) occurred during the early accumulation and melt period (in 29% and 23%, respectively), mostly present in areas with grassland and forest, slopes of 20–40°, areas exposed NW or areas with a topographic roughness index of ?0.25 to 0 m. These findings may give recommendations for defining more effective model parameterization strategies and guide future work, in which simulated and MODIS SCA may be combined to generate improved products for SCA monitoring in Alpine catchments.     

River ice cover influence on sediment transportation at present and under projected hydroclimatic conditions

A large number of rivers are frozen annually, and the river ice cover has an influence on the geomorphological processes. These processes in cohesive sediment rivers are not fully understood. Therefore, this paper demonstrates the impact of river ice cover on sediment transport, i.e. turbidity, suspended sediment loads and erosion potential, compared with a river with ice‐free flow conditions. The present sediment transportation conditions during the annual cycle are analysed, and the implications of climate change on wintertime geomorphological processes are estimated. A one‐dimensional hydrodynamic model has been applied to the Kokemäenjoki River in Southwest Finland. The shear stress forces directed to the river bed are simulated with present and projected hydroclimatic conditions. The results of shear stress simulations indicate that a thermally formed smooth ice cover diminishes river bed erosion, compared with an ice‐free river with similar discharges. Based on long‐term field data, the river ice cover reduces turbidity statistically significantly. Furthermore, suspended sediment concentrations measured in ice‐free and ice‐covered river water reveal a diminishing effect of ice cover on riverine sediment load. The hydrodynamic simulations suggest that the influence of rippled ice cover on shear stress is varying. Climate change is projected to increase the winter discharges by 27–77% on average by 2070–2099. Thus, the increasing winter discharges and possible diminishing ice cover periods both increase the erosion potential of the river bed. Hence, the wintertime sediment load of the river is expected to become larger in the future. Copyright © 2015 John Wiley & Sons, Ltd.     

Effect of permafrost thaw on the dynamics of lakes recharged by ice‐jam floods: case study of Yukon Flats,Alaska

Large river floods are a key water source for many lakes in fluvial periglacial settings. Where permeable sediments occur, the distribution of permafrost may play an important role in the routing of floodwaters across a floodplain. This relationship is explored for lakes in the discontinuous permafrost of Yukon Flats, interior Alaska, using an analysis that integrates satellite‐derived gradients in water surface elevation, knowledge of hydrogeology, and hydrologic modelling. We observed gradients in water surface elevation between neighbouring lakes ranging from 0.001 to 0.004. These high gradients, despite a ubiquitous layer of continuous shallow gravel across the flats, are consistent with limited groundwater flow across lake basins resulting from the presence of permafrost. Permafrost impedes the propagation of floodwaters in the shallow subsurface and constrains transmission to ‘fill‐and‐spill’ over topographic depressions (surface sills), as we observed for the Twelvemile‐Buddy Lake pair following a May 2013 ice‐jam flood on the Yukon River. Model results indicate that permafrost table deepening of 1–11 m in gravel, depending on watershed geometry and subsurface properties, could shift important routing of floodwater to lakes from overland flow (fill‐and‐spill) to shallow groundwater flow (‘fill‐and‐seep’). Such a shift is possible in the next several hundred years of ground surface warming and may bring about more synchronous water level changes between neighbouring lakes following large flood events. This relationship offers a potentially useful tool, well suited to remote sensing, for identifying long‐term changes in shallow groundwater flow resulting from thawing of permafrost. Copyright © 2015 John Wiley & Sons, Ltd.     

Effect of sediment resuspension on underwater light field in shallow lakes in the middle and lower reaches of the Yangtze River: A case study in Longgan Lake and Taihu Lake

Based on three continuous in situ underwater light field measurement under different wind waves conditions in Longgan Lake, Meiliang Bay of Taihu Lake in July 2003 and littoral zone near TLLER in July 2004, respectively, the effects of sediment resuspension caused by wind waves on PAR diffuse attenuation, absorption coefficients and euphotic depths are analyzed. In Longgan Lake, PAR diffuse attenuation coefficients during small, middle and large wind waves were 1.74, 2.02 and 2.45 m-1, respectively, and the corresponding PAR spectral diffuse attenuations ranged from 0.98 to 2.97, 1.34 to 3.95 and 1.80 to 5.40 m-1, respectively. In Meiliang Bay, PAR diffuse attenuation coefficients were 2.63, 3.72, 4.37 m-1 during small, middle and large wind waves. PAR diffuse attenuation coefficients increased by 41% and 66% from small to middle, large wind waves, respectively. Absorption coefficients integrated over the range of PAR of CDOM, phytoplankton were 0.26, 0.28 m-1; 0.76, 0.49 m-1, respectively during middle and large wind waves. Absorption coefficients integrated over the range of PAR of non-algal particulate matter and total suspended particulate matter increased from 0.94 to 1.73 m-1, and from 1.70 to 2.22 m-1, respectively during middle and large wind waves. Relative contributions of absorption coefficients of non-algal particulate matter to total absorption coefficient integrated over the range of PAR were 44.14%, 65.05%, respectively, during middle and large wind waves. PAR euphotic depths decreased by 0.40, 0.19, 0.20 m from middle to large wind waves in Longganhu Lake, Meliang Bay and littoral zone near TLLER. Significant correlations were found between transparency, PAR diffuse attenuation coefficients, euphotic depths and total suspended paniculate matter, wind velocity, wave height. Most significant correlations were found between transparency, PAR diffuse attenuation coefficients, euphotic depths and inorganic suspended paniculate matter but low correlations for chlorophyll a, dissolved organic carbon. Increase of total suspended paniculate matter, especially inorganic suspended paniculate matter caused by wind waves was the dominant factor affecting underwater light field in shallow lakes in the middle and lower reaches of the Yangtze River based on observations at three stations.