Simulation of long-term thermal characteristics of three Estonian lakes

A one-dimensional surface energy-balance lake model, coupled to a thermodynamic model of lake ice, is used to simulate variations in the temperature of and evaporation from three Estonian lakes: Karujärv, Viljandi and Kirjaku. The model is driven by daily climate data, derived by cubic-spline interpolation from monthly mean data, and was run for periods of 8 years (Kirjaku) up to 30 years (Viljandi). Simulated surface water temperature is in good agreement with observations: mean differences between simulated and observed temperatures are from −0.8°C to +0.1°C. The simulated duration of snow and ice cover is comparable with observed. However, the model generally underpredicts ice thickness and overpredicts snow depth. Sensitivity analyses suggest that the model results are robust across a wide range (0.1–2.0 m−1) of lake extinction coefficient: surface temperature differs by less than 0.5°C between extreme values of the extinction coefficient. The model results are more sensitive to snow and ice albedos. However, changing the snow (0.2–0.9) and ice (0.15–0.55) albedos within realistic ranges does not improve the simulations of snow depth and ice thickness. The underestimation of ice thickness is correlated with the overestimation of snow cover, since a thick snow layer insulates the ice and limits ice formation. The overestimation of snow cover results from the assumption that all the simulated winter precipitation occurs as snow, a direct consequence of using daily climate data derived by interpolation from mean monthly data.

Modeling the photo-oxidation of dissolved organic matter by ultraviolet radiation in freshwater lakes: implications for mercury bioavailability

Uncertainties in projected ultraviolet (UV) radiation may lead to future increases in UV irradiation of freshwater lakes. Because dissolved organic carbon (DOC) is the main binding phase for mercury (Hg) in freshwater lakes, an increase in DOC photo-oxidation may affect Hg speciation and bioavailability. We quantified the effect of DOC concentration on the rate of abiotic DOC photo-oxidation for five lakes (DOC = 3.27–12.3 mg L−1) in Kejimkujik National Park, Canada. Samples were irradiated with UV-A or UV-B radiation over a 72-h period. UV-B radiation was found to be 2.36 times more efficient at photo-oxidizing DOC than UV-A, with energy-normalized rates of dissolved inorganic carbon (DIC) production ranging from 3.8 × 10−5 to 1.1 × 10−4 mg L−1 J−1 for UV-A, and from 6.0 × 10−5 to 3.1 × 10−4 mg L−1 J−1 for UV-B. Energy normalized rates of DIC production were positively correlated with DOC concentrations. Diffuse integrated attenuation coefficients were quantified in situ (UV-A Kd = 0.056–0.180 J cm−1; UV-B Kd = 0.015–0.165 J cm−1) and a quantitative depth-integrated model for yearly DIC photo-production in each lake was developed. The model predicts that, UV-A produces between 3.2 and 100 times more DIC (1521–2851 mg m−2 year−1) than UV-B radiation (29.17–746.7 mg m−2 year−1). Future increases in UV radiation may increase DIC production and increase Hg bioavailability in low DOC lakes to a greater extent than in high DOC lakes.

Commensalism in an agroecosystem: hydraulic redistribution by deep-rooted legumes improves survival of a droughted shallow-rooted legume companion

We investigated commensalism of water use among annual shallow-rooted and perennial deep-rooted pasture legumes by examining the effect of hydraulic lift by Cullen pallidum (N.T.Burb.) J.W.Grimes and Medicago sativa on growth, survival and nutrient uptake of Trifolium subterraneum L. A vertically split-root design allowed separate control of soil water in top and bottom soil. Thirty-five days after watering ceased in the top tube, but soil remained at field capacity in the bottom tube, an increase in shallow soil water content by hydraulic lift was 5.6 and 5.9 g kg−1 soil overnight for C. pallidum and M. sativa, respectively. Trifolium subterraneum in this treatment maintained higher leaf water potentials (with M. sativa) or exhibited a slower decline (with C. pallidum) than without companion perennial plants; and shoot biomass of T. subterraneum was 56% (with C. pallidum) and 67% (with M. sativa) of that when both top and bottom tubes were at field capacity. Uptake of rubidium (a potassium analog) and phosphorus by T. subterraneum was not facilitated by hydraulic lift. Interestingly, phosphorus content was threefold greater, and shoot biomass 1.5–3.3-fold greater when T. subterraneum was interplanted with C. pallidum compared with M. sativa, although dry weight of C. pallidum was much greater than that of M. sativa. This study showed that interplanting with deep-rooted perennial legumes has benefited the survival of T. subterraneum.

Global climatology of surface water temperatures of large lakes by remote sensing

Lake surface water temperatures (LSWTs) of 246 globally distributed large lakes were derived from Along-Track Scanning Radiometers (ATSR) for the period 1991–2011. The climatological cycles of mean LSWT derived from these data quantify on a global scale the responses of large lakes' surface temperatures to the annual cycle of forcing by solar radiation and the ambient meteorological conditions. LSWT cycles reflect the twice annual peak in net solar radiation for lakes between 1°S to 12°N. For lakes without a lake-mean seasonal ice cover, LSWT extremes exceed air temperatures by 0.5–1.7 °C for maximum and 0.7–1.9 °C for minimum temperature. The summer maximum LSWTs of lakes from 25°S to 35°N show a linear decrease with increasing altitude; −3.76 ± 0.17 °C km−1 (inline image = 0.95), marginally lower than the corresponding air temperature decrease with altitude −4.15 ± 0.24 °C km−1 (inline image = 0.95). Lake altitude of tropical lakes account for 0.78–0.83 (inline image) of the variation in the March to June LSWT–air temperature differences, with differences decreasing by 1.9 °C as the altitude increases from 500 to 1800 m above sea level (a.s.l.) We define an ‘open water phase’ as the length of time the lake-mean LSWT remains above 4 °C. There is a strong global correlation between the start and end of the lake-mean open water phase and the spring and fall 0 °C air temperature transition days, (inline image = 0.74 and 0.80, respectively), allowing for a good estimation of timing and length of the open water phase of lakes without LSWT observations. Lake depth, lake altitude and distance from coast further explain some of the inter-lake variation in the start and end of the open water phase.

A moving point approach to model shallow ice sheets: a study case with radially-symmetrical ice sheets

Predicting the evolution of ice sheets requires numerical models able to accurately track the migration of ice sheet continental margins or grounding lines. We introduce a physically based moving point approach for the flow of ice sheets based on the conservation of local masses. This allows the ice sheet margins to be tracked explicitly and the waiting time behaviours to be modelled efficiently. A finite difference moving point scheme is derived and applied in a simplified context (continental radially-symmetrical shallow ice approximation). The scheme, which is inexpensive, is validated by comparing the results with moving-margin exact solutions and steady states. In both cases the scheme is able to track the position of the ice sheet margin with high precision.

A moving-point approach to model shallow ice sheets: a study case with radially symmetrical ice sheets

Predicting the evolution of ice sheets requires numerical models able to accurately track the migration of ice sheet continental margins or grounding lines. We introduce a physically based moving-point approach for the flow of ice sheets based on the conservation of local masses. This allows the ice sheet margins to be tracked explicitly. Our approach is also well suited to capture waiting-time behaviour efficiently. A finite-difference moving-point scheme is derived and applied in a simplified context (continental radially symmetrical shallow ice approximation). The scheme, which is inexpensive, is verified by comparing the results with steady states obtained from an analytic solution and with exact moving-margin transient solutions. In both cases the scheme is able to track the position of the ice sheet margin with high accuracy.