Modelling phosphorus dynamics in multi-branch river systems: A study of the Black River, Lake Simcoe, Ontario, Canada

High rates of nutrient loading from agricultural and urban development have resulted in surface water eutrophication and groundwater contamination in regions of Ontario. In Lake Simcoe (Ontario, Canada), anthropogenic nutrient contributions have contributed to increased algal growth, low hypolimnetic oxygen concentrations, and impaired fish reproduction. An ambitious programme has been initiated to reduce phosphorus loads to the lake, aiming to achieve at least a 40% reduction in phosphorus loads by 2045. Achievement of this target necessitates effective remediation strategies, which will rely upon an improved understanding of controls on nutrient export from tributaries of Lake Simcoe as well as improved understanding of the importance of phosphorus cycling within the lake. In this paper, we describe a new model structure for the integrated dynamic and process-based model INCA-P, which allows fully-distributed applications, suited to branched river networks. We demonstrate application of this model to the Black River, a tributary of Lake Simcoe, and use INCA-P to simulate the fluxes of P entering the lake system, apportion phosphorus among different sources in the catchment, and explore future scenarios of land-use change and nutrient management to identify high priority sites for implementation of watershed best management practises.

Immune factors and fatty acid composition in human milk from river/lake, coastal and inland regions of China

Breast milk fatty acid composition may be affected by maternal diet during gestation and lactation. The influence of dietary and breast milk fatty acids on breast milk immune factors is poorly defined. We determined the fatty acid composition and immune factor concentrations of breast milk from women residing in river & lake, coastal, and inland regions of China, which differ in their consumption of lean fish and oily fish. Breast milk samples were collected on days 3 to 5 (colostrum), 14 and 28 post-partum and analysed for soluble CD14 (sCD14), transforming growth factor (TGF)-β1, TGF-β2, secretory immunoglobulin A (sIgA) and fatty acids. The fatty acid composition of breast milk differed between regions and with time post-partum. The concentrations of all four immune factors in breast milk decreased over time, with sCD14, sIgA and TGF-β1 being highest in colostrum in the river & lake region. Breast milk DHA and arachidonic acid (AA) were positively associated, and γ-linolenic acid and EPA negatively associated, with the concentrations of each of the four immune factors. In conclusion, breast milk fatty acids and immune factors differ between regions in China characterised by different patterns of fish consumption and change during the course of lactation. A higher breast milk DHA and AA concentration is associated with higher concentrations of immune factors in breast milk, suggesting a role for these fatty acids in promoting gastrointestinal and immune maturation of the infant.

Understanding lake and catchment history as a tool for integrated lake management

Sustainable lake management for nutrient-enriched lakes must be underpinned by an understanding of both the functioning of the lake, and the origins of changes in nutrient loading from the catchment. To date, limnologists have tended to focus on studying the impact of nutrient enrichment on the lake biota, and the dynamics of nutrient cycling between the water column, biota and sediments within the lake. Relatively less attention has been paid to understanding the specific origins of nutrient loading from the catchment and nutrient transport pathways linking the lake to its catchment. As such, when devising catchment management strategies to reduce nutrient loading on enriched lakes, assumptions have been made regarding the relative significance of non-point versus point sources in the catchment. These are not always supported by research conducted on catchment nutrient dynamics in other fields of freshwater science. Studies on nutrient enrichment in lakes need to take account of the history of catchment use and management specific to each lake in order to devise targeted and sustainable management strategies to reduce nutrient loading to enriched lakes. Here a modelling approach which allows quantification of the relative contribution of nutrients from each specific point and non-point catchment source over the course of catchment history is presented. The approach has been applied to three contrasting catchments in the U.K. for the period 1931 to present. These are the catchment of Slapton Ley in south Devon, the River Esk in Cumbria and the Deben Estuary in Suffolk. Each catchment showed marked variations in the nature and intensity of land use and management. The model output quantifies the relative importance of point source versus non-point livestock and land use sources in each of the catchments, and demonstrates the necessity for an understanding of site-specific catchment history in devising suitable management strategies for the reduction of nutrient loading on enriched lakes.

Pediastrum species as potential indicators of lake-level change in tropical South America

We explored the potential for using Pediastrum (Meyen), a genus of green alga commonly found in palaeoecological studies, as a proxy for lake-level change in tropical South America. The study site, Laguna La Gaiba (LLG) (17°45′S, 57°40′W), is a broad, shallow lake located along the course of the Paraguay River in the Pantanal, a 135,000-km2 tropical wetland located mostly in western Brazil, but extending into eastern Bolivia. Fourteen surface sediment samples were taken from LLG across a range of lake depths (2-5.2 m) and analyzed for Pediastrum. We found seven species, of which P. musteri (Tell et Mataloni), P. argentiniense (Bourr. et Tell), and P. cf. angulosum (Ehrenb.) ex Menegh. were identified as potential indicators of lake level. Results of the modern dataset were applied to 31 fossil Pediastrum assemblages spanning the early Holocene (12.0 kyr BP) to present to infer past lake level changes qualitatively. Early Holocene (12.0-9.8 kyr BP) assemblages do not show a clear signal, though abundance of P. simplex (Meyen) suggests relatively high lake levels. Absence of P. musteri, characteristic of deep, open water, and abundance of macrophyte-associated taxa indicate lake levels were lowest from 9.8 to 3.0 kyr BP. A shift to wetter conditions began at 4.4 kyr BP, indicated by the appearance of P. musteri, though inferred lake levels did not reach modern values until 1.4 kyr BP. The Pediastrum-inferred mid-Holocene lowstand is consistent with lower precipitation, previously inferred using pollen from this site, and is also in agreement with evidence for widespread drought in the South American tropics during the middle Holocene. An inference for steadily increasing lake level from 4.4 kyr BP to present is consistent with diatom-inferred water level rise at Lake Titicaca, and demonstrates coherence with the broad pattern of increasing monsoon strength from the late Holocene until present in tropical South America.

Late Quaternary highstands at Lake Chilwa, Malawi: frequency, timing and possible forcing mechanisms in the last 44ka

The extensive shoreline deposits of Lake Chilwa, southern Malawi, a shallow water body today covering 600 km2 of a basin of 7500 km2, are investigated for their record of late Quaternary highstands. OSL dating, applied to 36 samples from five sediment cores from the northern and western marginal sand ridges, reveal a highstand record spanning 44 ka. Using two different grouping methods, highstand phases are identified at 43.7–33.3 ka, 26.2–21.0 ka and 17.9–12.0 ka (total error method) or 38.4–35.5 ka, 24.3–22.3 ka, 16.2–15.1 ka and 13.5–12.7 ka (Finite Mixture Model age components) with two further discrete events recorded at 11.01 ± 0.76 ka and 8.52 ± 0.56 ka. Highstands are comparable to the timing of wet phases from other basins in East and southern Africa, demonstrating wet conditions in the region before the LGM, which was dry, and a wet Lateglacial, which commenced earlier in the southern compared to northern hemisphere in East Africa. We find no evidence that wet phases are insolation driven, but analysis of the dataset and GCM modelling experiments suggest that Heinrich events may be associated with enhanced monsoon activity in East Africa in both timing and as a possible causal mechanism.

Tropical climates at the Last Glacial Maximum: a new synthesis of terrestrial palaeoclimate data. I. Vegetation, lake-levels and geochemistry.

Palaeodata in synthesis form are needed as benchmarks for the Palaeoclimate Modelling Intercomparison Project (PMIP). Advances since the last synthesis of terrestrial palaeodata from the last glacial maximum (LGM) call for a new evaluation, especially of data from the tropics. Here pollen, plant-macrofossil, lake-level, noble gas (from groundwater) and δ18O (from speleothems) data are compiled for 18±2 ka (14C), 32 °N–33 °S. The reliability of the data was evaluated using explicit criteria and some types of data were re-analysed using consistent methods in order to derive a set of mutually consistent palaeoclimate estimates of mean temperature of the coldest month (MTCO), mean annual temperature (MAT), plant available moisture (PAM) and runoff (P-E). Cold-month temperature (MAT) anomalies from plant data range from −1 to −2 K near sea level in Indonesia and the S Pacific, through −6 to −8 K at many high-elevation sites to −8 to −15 K in S China and the SE USA. MAT anomalies from groundwater or speleothems seem more uniform (−4 to −6 K), but the data are as yet sparse; a clear divergence between MAT and cold-month estimates from the same region is seen only in the SE USA, where cold-air advection is expected to have enhanced cooling in winter. Regression of all cold-month anomalies against site elevation yielded an estimated average cooling of −2.5 to −3 K at modern sea level, increasing to ≈−6 K by 3000 m. However, Neotropical sites showed larger than the average sea-level cooling (−5 to −6 K) and a non-significant elevation effect, whereas W and S Pacific sites showed much less sea-level cooling (−1 K) and a stronger elevation effect. These findings support the inference that tropical sea-surface temperatures (SSTs) were lower than the CLIMAP estimates, but they limit the plausible average tropical sea-surface cooling, and they support the existence of CLIMAP-like geographic patterns in SST anomalies. Trends of PAM and lake levels indicate wet LGM conditions in the W USA, and at the highest elevations, with generally dry conditions elsewhere. These results suggest a colder-than-present ocean surface producing a weaker hydrological cycle, more arid continents, and arguably steeper-than-present terrestrial lapse rates. Such linkages are supported by recent observations on freezing-level height and tropical SSTs; moreover, simulations of “greenhouse” and LGM climates point to several possible feedback processes by which low-level temperature anomalies might be amplified aloft.

Simulated climate and biomes of Africa during the Late Quaternary: comparison with pollen and lake status data

New compilations of African pollen and lake data are compared with climate (CCM1, NCAR, Boulder) and vegetation (BIOME 1.2, GSG, Lund) simulations for the last glacial maximum (LGM) and early to mid-Holocene (EMH). The simulated LGM climate was ca 4°C colder and drier than present, with maximum reduction in precipitation in semi-arid regions. Biome simulations show lowering of montane vegetation belts and expansion of southern xerophytic associations, but no change in the distribution of deserts and tropical rain forests. The lakes show LGM conditions similar or drier than present throughout northern and tropical Africa. Pollen data indicate lowering of montane vegetation belts, the stability of the Sahara, and a reduction of rain forest. The paleoenvironmental data are consistent with the simulated changes in temperature and moisture budgets, although they suggest the climate model underestimates equatorial aridity. EMH simulations show temperatures slightly less than present and increased monsoonal precipitation in the eastern Sahara and East Africa. Biome simulations show an upward shift of montane vegetation belts, fragmentation of xerophytic vegetation in southern Africa, and a major northward shift of the southern margin of the eastern Sahara. The lakes indicate conditions wetter than present across northern Africa. Pollen data show an upward shift of the montane forests, the northward shift of the southern margin of the Sahara, and a major extension of tropical rain forest. The lake and pollen data confirm monsoon expansion in eastern Africa, but the climate model fails to simulate the wet conditions in western Africa.

Evaluation of modelled regional water balance using lake status data: a comparison of 6ka simulations with the NCAR CCM

Variations in lake area and depth reflect climatically induced changes in the water balance of overflowing as well as closed lakes. A new global data base of lake status has been assembled, and is used to compare two simulations for 6 ka (6000 yr ago) made with successive R15 versions of the NCAR Community Climate Model (CCM). Simulated water balance was expressed as anomalies of annual precipitation minus evaporation (P-E); observed water balance as anomalies of lake status. Comparisons were made visually, by comparing regional averages, and by a statistic that compares the signs of simulated P-E anomalies (smoothly interpolated to the lake sites) with the status anomalies. Both CCM0 and CCM1 showed enhanced Northern-Hemisphere monsoons at 6 ka. Both underestimated the effect, but CCM1 fitted the spatial patterns better. In the northern mid- and high-latitudes the two versions differed more, and fitted the data less satisfactorily. CCM1 performed better than CCM0 in North America and central Eurasia, but not in Europe. Both models (especially CCM0) simulated excessive aridity in interior Eurasia. The models were systematically wrong in the southern mid-latitudes. Problems may have been caused by inadequate treatment of changes in sea-surface conditions in both models. Palaeolake status data will continue to provide a benchmark for the evaluation of modelling improvements.