
==== Front
Environ Monit Assess
Environ Monit Assess
Environmental Monitoring and Assessment
0167-6369
1573-2959
Springer International Publishing Cham

39230772
13048
10.1007/s10661-024-13048-5
Research
Anthropogenic and climatic impacts on historic sediment, carbon, and phosphorus accumulation rates using 210Pbex and 137Cs in a sub-watershed linked to Zarivar Lake, Iran
http://orcid.org/0000-0001-5599-4071
Khodadadi Maral khodadm@miamioh.edu

12
Gibbs Max 3
Swales Andrew 3
Toloza Arsenio 4
Blake William H. 5
1 https://ror.org/05nbqxr67 grid.259956.4 0000 0001 2195 6763 Department of Geology and Environmental Earth Science, Miami University, Oxford, OH USA
2 https://ror.org/05cebxq10 0000 0004 7433 9111 Nuclear Agriculture Research School, Nuclear Science and Technology Research Institute (NSTRI), Karaj, Iran
3 https://ror.org/04hxcaz34 grid.419676.b 0000 0000 9252 5808 National Institute of Water and Atmospheric Research (NIWA), Hamilton, New Zealand
4 https://ror.org/02zt1gg83 grid.420221.7 0000 0004 0403 8399 Soil and Water Management & Crop Nutrition Section and Laboratory, Department of Nuclear Sciences and Applications, Joint FAO/IAEA Division, International Atomic Energy Agency, Vienna, Austria
5 https://ror.org/008n7pv89 grid.11201.33 0000 0001 2219 0747 School of Geography, University of Plymouth, Drake Circus, Plymouth, PL4 8AA UK
4 9 2024
4 9 2024
2024
196 10 88729 3 2024
23 8 2024
© This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
To estimate a watershed’s response to climate change, it is crucial to understand how human activities and climatic extremes have interacted over time. Over the last century, the Zarivar Lake watershed, Iran, has been subjected to various anthropogenic activates, including deforestation and inappropriate land-management practices alongside the implementation of conservation measures like check dams. To understand the effects of these changes on the magnitude of sediment, organic carbon (OC), and phosphorus supplies in a small sub-watershed connected to the lake over the last century, a lake sediment core was dated using 210Pbex and 137Cs as geochronometers. The average mass accumulation rate (MAR), organic carbon accumulation rates (OCAR), and particulate phosphorus accumulation rates (PPAR) of the sediment core were determined to be 6498 ± 2475, 205 ± 85, and 8.9 ± 3.3 g m−2 year−1, respectively. Between the late 1970s and early 1980s, accumulation rates were significantly higher than their averages at 7940 ± 3120, 220 ± 60, and 12.0 ± 2.8 g m−2 year−1 respectively. During this period, the watershed underwent extensive deforestation (12%) on steep slopes, coinciding with higher mean annual precipitations (more than double). Conversely, after 2009, when check dams were installed in the sub-watershed, the sediment load to the lake became negligible. The results of this research indicate that anthropogenic activities had a pronounced effect on MAR, OCAR, and PPAR, causing them to fluctuate from negligible amounts to values twice the averages over the last century, amplified by climatic factors. These results imply that implementing climate-smart watershed management strategies, such as constructing additional check dams and terraces, reinforcing restrictions on deforestation, and minimum tillage practices, can facilitate protection of lacustrine ecosystems under accelerating climate change conditions.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s10661-024-13048-5.

Keywords

Sediment chronology
Particulate P
Geochemical indices
Land-use change
δ13C
Inland lakes
http://dx.doi.org/10.13039/501100004493 International Atomic Energy Agency 19167 Khodadadi Maral issue-copyright-statement© Springer Nature Switzerland AG 2024
==== Body
pmcIntroduction

The increased influx of fine sediments and their associated chemical constituents into reservoirs can indeed lead to adverse effects on water quality, aquatic life (Kemp et al., 2011), reservoir storage capacity (Kondolf et al., 2014), and primary production (Harmon et al., 2014). Climatic extremes (Kortelainen et al., 2013) and anthropogenic activities, particularly land-use change (LUC; e.g., Foster et al., 2003; Valero-Garcés et al., 2006; Anderson et al., 2013; Lin et al., 2021), are known to influence the amounts of fine sediments entering reservoirs. However, understanding the effect of human–environment interactions at the basin scale over historical times can be highly complex (Dearing et al., 2008). This complexity can be attributed to the limited availability of precise data on past climatic events and human activities, as well as the intricate interplay between human and climate influences. In this context, sediment cores can provide high-resolution temporal records of landscape response, aiding understanding how anthropogenic and climatic influences propagate through space and time (Canuel et al., 2017; Foster et al., 2003). Such knowledge can be applied to understand the current state of ecosystems (Zan et al., 2012) and predict the future ecosystem responses to climate change (Canuel et al., 2017), facilitating the implementation of effective conservation strategies and climate-smart agriculture practices to assure food security. In other words, reconstructing past environmental conditions enables informed future management decisions to protect and restore freshwater ecosystems in the face of increasing anthropogenic pressures and climate change.

Elevated dissolved organic carbon (OC) loadings to lakes cause heterotrophic conditions, where the rate of ecosystem respiration surpasses primary productivity, leading to a net efflux of CO2 into the atmosphere (Cole et al., 2007). In addition, lakes tend to accumulate large amounts of OC in the sediment, thus playing a vital role in the carbon cycle (Anderson et al., 2013; Kortelainen et al., 2004; Tranvik et al., 2009). The amount of accumulated OC in the sediment signifies carbon sequestration in short-to-long-term time scales (Downing et al., 2008). Organic carbon accumulation rate (OCAR) in lakes represents the difference between inputs from aquatic primary production and catchment fluxes and losses through respiration and outflow (Tranvik et al., 2009). Overall, understanding how OCAR may have changed over time, particularly in response to changes in land use and climate, could provide valuable insights into how future global change processes might affect OCAR in lakes (Anderson et al., 2013; Kortelainen et al., 2013; Lin et al., 2021) and, consequently, their storage capacity (Kondolf et al., 2014). Besides, excessively high concentrations of nutrients like nitrogen (N) and phosphorus (P) stimulate algal growth and are the most common cause of eutrophication in freshwater reservoirs (Sterner, 2008). P, in particular, is a limiting nutrient for primary production in freshwater ecosystems, playing a crucial role in influencing the trophic level of water bodies (Sterner, 2008). Additionally, in shallow eutrophic lakes, sediments, both internal sediment P and suspended particulate matter P, are the predominant sources of P (Radbourne et al., 2019), with particulate phosphorus (PP) constituting the largest portion of total phosphorus (TP; Yang et al., 2020). PP in lakes can originate from terrestrial sources or the decomposition of aquatic plants (Liu et al., 2020). The sediment-associated PP buried in deposited sediments can be released as soluble phosphate into the water column under anoxic conditions and contribute to eutrophication in water bodies (Heathcote et al., 2015). These factors underscore the importance of identifying procedures that impact the transport and fate of OC and PP to maintain the ecological health of water bodies. The fallout radionuclide (FRN), including unsupported lead-210 excess (210Pbex) and cesium-137 (137Cs), can be used as geochronometers to date sediment cores and estimate changes in mass accumulation rates (MAR) over time (Appleby, 2008). The reduction rate of 210Pbex activity with depth establishes an age-depth relationship in the core (IAEA, 2003). While 210Pbex has been proven useful to reconstruct geochronology in sediments deposited over the last 100 years, due to its natural radioactive decay rate in a sediment column, an independent marker such as 137Cs, is needed to verify calculations resulting from 210Pbex activity (Foucher et al., 2021; Smith, 2001). More precisely, 210Pbex serves as a tool to date the entire sediment column, whereas 137Cs is considered to be an event marker that is only employed for particular dates to validate the 210Pbex dating curve. These “date stamps” are indicated by distinct peaks, e.g., that of 1963 indicating the highest fallout on the global scale and 1986 corresponding with Chernobyl incident (Northern Hemisphere). This FRN dating technique proved to be a particularly useful tool for comprehending climate change impacts (Chen et al., 2014), to determine the impacts of historical events during the previous century (Harmon et al., 2014; Sanders et al., 2006), and to indicate OC and other nutrients’ burial rates in the last 100 years (Anderson et al., 2013; Breithaupt et al., 2014; Sanders et al., 2014, 2017).

Zarivar Lake, Kurdistan Province, Iran, is a freshwater reservoir and a habitat for diverse plant and animal life (IMA, 2008), including threatened and endangered species, e.g., Oxyura leucocephala (Maroufi & Ehyaei, 2016). The soil erosion rate estimated by the 137Cs method in cultivated lands situated in one of Zarivar Lake’s sub-watersheds was found to be five times higher than that of the forested areas in the same watershed, at 26 Mg ha⁻1 year⁻1 (Khodadadi et al., 2024) compared to 5 Mg ha−1 year−1 (Khodadadi et al., 2023) and less than twice the country’s average of 15 Mg ha−1 year−1. One main challenge the lake encounters is the continuous input of sediments and associated nutrient/contaminants via the river network of the watershed (Sharifinia et al., 2013). Torabi Kachoosangi et al. (2020) reported the average sedimentation rate of 5.6 mm year−1 based on 210Pbex activity profile in one core retrieved from Zarivar Lake. The lake has experienced hypereutrophic conditions over the last few decades, which can accelerate the primary productivity of the lake (Ebrahimi Mohammadi et al., 2012; Hamidian & Hassanzadeh, 2013; Sharifinia et al., 2013), posing threats to aquatic life and both local inhabitants. Moreover, over the last century, the Zarivar Lake watershed has been subjected to various anthropogenic activities, including deforestation and inappropriate land-management practices alongside the implementation of conservation measures like check dams. Nevertheless, agricultural practices such as minimum tillage, contour farming, and the construction of terraces and check dams can serve as adaptive management measures to mitigate the impacts of climate change. Zarivar Lake was selected as the study site due to its high soil erosion and sedimentation rates, significant history of land-use changes, and the construction of numerous check dams in 2008–2009, providing an opportunity to analyze their hydrologic impacts. Thus, the unique characteristics of the watershed make it an ideal site for studying anthropogenic impacts under changing climate conditions. Furthermore, the region is expected to undergo substantial climate changes in the future, which is likely to influence soil erosion rates, with modifications expected in the erosive power of rainfall, soil moisture, plant biomass, litter cover, soil erodibility, and land use, especially in semi-arid areas (Nunes & Nearing, 2011). These changes, in turn, will likely impact soil erosion, the amounts of sediment entering the lakes, nutrient delivery, and internal lake production.

Zarivar Lake watershed has been the subject of a number of studies, including Ebrahimi Mohammadi et al. (2012), Hamidian and Hassanzadeh (2013), Sharifinia et al. (2013), Torabi Kachoosangi et al. (2020), and Khodadadi et al., (2023, 2024). However, to date, there has been no investigation into the accumulation rates of OC and PP in the lake, nor has the role of anthropogenic activities in upstream watersheds on MAR and nutrient burial in sediments been explored over time. In addition, despite the significant impacts of LUC on sediment and nutrient supply, as indicated by Anderson et al. (2013) and Lin et al. (2021), this effect has only been examined at a limited number of sites worldwide and thus needs further investigation. Thus, the objectives of the study were to (i) quantify land-use/land-cover (LULC) changes using aerial photos and satellite data between 1955 and 2008; (ii) estimate the sediment, organic carbon, and particulate phosphorus accumulation rates using 210Pbex; and (iii) understand anthropogenic and climatic impacts on sediment, OC, and PP accumulation rates using LULC changes, mean annual precipitation, δ13C, and elemental indices over the past decades in one of sub-watersheds of Zarivar Lake.

Materials and methods

Study area

Zarivar Lake (35° 33′ 15″ N and 46° 7′ 25″ E) is a natural shallow reservoir situated approximately 15 km east of Iran-Iraq border (Fig. 1). The lake, which is inhabited by aquatic plants within its shallow regions, e.g., Typha spp. (i.e., cattails or common bulrush) and Phragmites spp. (common reed), has a surface area of ca. 18.3 km2 and a maximum depth greater than 6 m (IMA, 2008). Most of the rivers in the lake watershed are seasonal (non-permanent). A dike was constructed at the outlet of the watershed in 1995; therefore, the lake only drains at high water levels. Additionally, in 2004, the Ghezelchesoo River was partly diverted to an artificial seasonal tributary which flows into the lake (IMA, 2008). Subsequent to the construction of the diversion dam redirecting the Ghezelchesoo River, the lake’s watershed area expanded from 107.2 to about 296.4 km2 (Fig. 1). The mean annual precipitation (MAP) and temperature (MAT) are at 991 mm and 12.8 °C, respectively (Iran Meteorological Organization).Fig. 1 The location of Zarivar lake watershed, Z3 watershed, and the sediment cores in marshland (this study) and in lake (Torabi Kachoosangi et al., 2020) in Iran. The coordinates of the cores are also provided

The present study was conducted in an upstream sub-watershed connected to the lake in the western part of the watershed, namely, the “Z3” sub-watershed (Fig. 1), encompassing a total area of approximately 2.97 km2. The selection of this watershed was motivated by findings from Ebrahimi Mohammadi et al. (2012), which highlighted that the northern and western tributaries contributed significantly higher amounts of suspended sediment and nutrients to Zarivar Lake. The watershed has an average altitude of 1518 m a.s.l. (varying between 1296 and 1884 m a.s.l.) and an average slope of 31% (ranging from to 0.1 to 82%). While the parent materials in ridges consist of metamorphic rocks of complexly folded black shale and minor metamorphosed limestone and sandstones, in the lowlands, Quaternary sediments, mainly lacustrine deposits, can be found (Geological Survey and Mineral Exploration of Iran). The sub-watershed’s soils fall into the Typic Haploxeralfs and Typic Haploxerepts categories based on Soil Taxonomy (2014). In the WRB classification, they are classified as Haplic Luvisols and Haplic Cambisols (WRB, 2014). The topsoil textures are predominantly silty loam and loam with significant variation in stone content (IMA, 2008). The landscape in this region is characterized by mountains covered by oak forests, hills, and plains which are used mainly for agriculture.

The sub-watershed encompasses three distinct land-use categories, comprising deciduous oak forests and irrigated and dry-farming agricultural lands (Fig. 1). Within the irrigated zones, common crops include alfalfa, tobacco, and strawberries. Conversely, dry-farming lands feature cereal crops such as wheat, barley, lentils, and vineyards. Notably, vineyards are predominantly situated on sloped terrains. It is important to mention that, after the harvest of cereal crops, the soil surface typically remains exposed, making it vulnerable to erosion during heavy rainfall events. Generally, different soil erosion features, including sheet, rill, stream bank, and gully (only one stabilized gully was found), have been observed in the sub-watershed. However, conservation measures have been taken to reduce sediment runoff in the sub-watershed, such as terraces and check dams in the sub-watershed in 2008–2009. The sub-watershed is an ungauged watershed without streamflow data.

Sediment core sampling and laboratory analyses

In 2016, three sediment cores (in excess of 50 cm depth) were collected from the marshland, away from the shoreline (Fig. 1) where the water depth reaches around 3 m in wet seasons, using a purpose-built gravity corer comprising stainless-steel coring tube and a transparent plexiglass tube with a 9.4 cm diameter. This setup facilitated the visual assessment of changes in sediment texture. For this study, one of the three cores (core#2 in Fig. 1) was sliced into 1-cm-depth increments, except for the upper section of the core (i.e., 0–7 cm), where dense plant residuals made it impossible. Sediment from the top 7 cm was separated from plant material and composited for further analysis. To minimize contamination and disturbance during sample collection and preparation, we followed Gibbs (2014) sampling protocols. Water samples above the sediment columns were collected and placed in dark bottles for total dissolved phosphorus (TDP; ascorbic acid blue method, Murphy & Riley, 1962) analyses. All samples were kept in an ice box (− 4 °C) during transport to avoid changes in chemical properties.

Increment sediment samples were oven-dried at 60 °C, disaggregated, and gently ground using a mortar and pestle. The samples were analyzed for 137Cs, 210Pb, and 226Ra by gamma spectrometry using N-type HPGe detector at the Soil Science Unit at the IAEA Seibersdorf Laboratory. Prior to analyzing 210Pb, samples were sealed for 1 month to achieve equilibrium between 226Ra and its daughter 214Pb. 137Cs, total 210Pb, and 214Pb activities were determined from the net peak areas of gamma rays at 661.6, 46.5, and 352 keV, respectively. The 210Pbex mass activity was calculated by subtracting 226Ra from the total 210Pb measurement. Depending on the FRN activities, counting times ranged between 12 and 24 h, to acquire a precision of 5% to a maximum of 20% at the 95% level of confidence.

The sources of buried OC can be identified using carbon stable isotope signatures, i.e., δ13C (e.g., Sanders et al., 2014). OC and δ13C were measured in 27 of 45 increment sediment samples, which were randomly selected from various depths within the core. Initially, inorganic carbonates were removed using 10% HCl. The OC content in each sample was determined as loss on ignition (LOI) by combustion of dried soil at 500 °C for 3 h, which was then multiplied by 0.47 (Perie & Ouimet, 2008) to estimate total organic carbon. The δ13C isotopic signature of organic carbon was measured in acidified samples using a Delta V Plus continuous flow isotope ratio mass spectrometer (CF-IRMS) linked to a Flash 2000 elemental analyzer on a MAS 200 R autosampler (Thermo-Fisher Scientific, Bremen, Germany) at Environmental Stable Isotope Laboratory, the National Institute of Water and Atmospheric Research, Wellington, New Zealand. To make δ13C values with different ages comparable, the Suess effect—a decrease in δ13C signature of atmospheric carbon dioxide (CO2) due to burning fossil fuels with δ13C depleted CO2 emissions (commencing in the 1700s)—correction was required (Verburg, 2007). This correction was performed using Verburg’s (2007) sixth-order polynomial equation to calculate isotopic depletion. Then, the absolute δ13C value (− 8.65‰) of present-day CO2 (the year 2016) was added to compute the change in δ13C isotopic values for each year of the core (Gibbs, 2014).

Sediment increment samples were also analyzed for major and minor elements using wavelength dispersive X-ray fluorescence (WD-XRF) with the OMNIAN application on an Axios Max instrument (Malvern PANalytical, Malvern, UK) at the University of Plymouth, Plymouth, UK. The procedure followed the method described by Wynants et al. (2021). Prior to pelleting, samples were homogenized by milling for 20 min at 300 rpm to reduce shadowing effects and prevent the preferential analysis of finer particles (Willis et al., 2011). Instrument calibration and drift were evaluated using internal quality control protocols with a multi-element glass sample provided by the manufacturer. Triplicates of randomly selected samples were prepared to assess the repeatability of the method. The instrument measurements were validated with stream sediment certified reference material (GBW07318, LGC, Middlesex, UK). Geochemical constituents of the samples measured were P, Na, Mg, Al, Si, S, Cl, K, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Ga, Br, Rb, Sr, Y, Zr, Nb, Sn, I, Cs, Ba, La, Ce, Nd, Hf, Ta, Hg, Pb, Th, and U.

Core dating and MAR, OCAR, and PPAR calculation

The constant rate supply (CRS) model was applied to establish the age-depth relationship according to the 210Pbex depth profile (Appleby & Oldfield, 1978). The CRS model is considered to be a very practical dating approach and is widely used and validated by numerous independent approaches (Appleby, 2008; Foucher et al., 2021). The CRS model presumes that the input of 210Pbex to the sediment is constant throughout time (Krishnaswamy et al., 1971), which is a widely accepted assumption (e.g., Appleby & Oldfield, 1992) and that MAR can be modeled using the 210Pbex profile shape and radioactive decay (Appleby & Oldfield, 1978). We calculated CRS ages and MARs for each increment depth using the equations provided by Appleby (2001). MAR was reported as g cm−2 year−1 or kg m−2 year−1 and converted to sediment accumulation rate (SAR; mm year−1) using corresponding bulk density at each interval. Wet and dry bulk densities were calculated for sediment increments using the core radius, increment depth, and weight of each increment. For each increment depth, the mid-point of that layer was used in calculations.

OCAR and PPAR were computed by simply multiplying the total sediment mass accumulation rate at each increment depth by the respective proportion of OC or PP and expressed as g m−2 year−1 (e.g., see Waters et al., 2019). Subsequently, to evaluate significant differences between estimated accumulation rates in different sections of the core, we used a Kruskal–Wallis test followed by Dunn’s post-hoc test with Bonferroni correction in R. Significance was assessed at p < 0.05 level.

Geochemical and geoaccumulation indices

Using elemental compositions, the geochemical indices were calculated to assess sediment maturity, chemical weathering, and depositional environment in different sections of the core. Two indices were applied to deduce the level of chemical maturity of the sediment column: the index of chemical maturity (ICV; see Table 1) and the SiO2/Al2O3 ratio. Additionally, the chemical index of alteration (CIA; refer to Table 1) was used to evaluate the extent of chemical weathering in the sediments. To understand the bottom lake conditions, i.e., oxic, sub-oxic, and anoxic, various elemental ratios such as U/Th, V/Cr, and Ni/Co were applied (Armstrong-Altrin et al., 2015; Jones & Manning, 1994). Table 1 The average (range) of geochemical indices for sediment maturity, chemical weathering, and depositional environment in different sections of the Zarivar Lake core

Depositional environment	Geochemical index	Observed mean (range)	Classification	Classification criteria	Reference	
Upper (0–15 cm)	Middle (15–25 cm)	Lower (25–50 cm)	
Sediment maturity	Si2O/Al2O3	2.74 (2.68–2.81)	2.75 (2.65–2.97)	3.05 (2.89–3.34)	Immature	Immature < 10

Mature > 10

	Bakkiaraj et al. (2010)	
ICV*	1.02 (0.94–1.24)	1.02 (0.92–1.17)	0.92 (0.83–1.08)	Immature	Immature < 2

Mature > 2

	Armstrong-Altrin et al. (2015)	
Chemical weathering	CIA**	69.37 (62.93–72.75)	70.04 (64.06–75.21)	73.55 (67.79–76.80)	Moderate to intense chemical weathering	Intense chemical weathering > 70	Zaid (2015)	
Deposition condition	U/Th	0.30 (0.24–0.34)	0.30 (0.24–0.36)	0.34 (0.29–0.40)	Oxic	Anoxic > 1.25

Oxic < 0.75

	Jones and Manning (1994)	
V/Cr	0.99 (0.96–1.00)	0.99 (0.96–1.04)	1.12 (0.98–1.32)	Oxic	Anoxic > 4.5

Oxic < 2

	
Ni/Co	4.64 (4.15–5.02)	4.70 (4.11–5.16)	5.06 (4.48–5.97)	Oxic to sub-oxic	Anoxic > 5 Oxic < 5	
*Index of chemical maturity; ICV = (Fe2O3 + K2O + Na2O + CaO + MgO + MnO + TiO2)/Al2O3 (Cox et al., 1995)

**The chemical index of alteration; CIA = [Al2O3/(Al2O3 + CaO + Na2O + K2O)] * 100

We also utilized the degree of contamination (Dc) and pollution load index (PLI) as key indicators to assess temporal variations in overall contamination levels as a function of depth. Initially, the contamination factor (CF) for each element was determined by dividing its concentration in the core by the corresponding average shale values (AVS; Wedepohl, 1995), which were used as the background value (El-Safa et al., 2022). The degree of contamination (Dc) was calculated as the sum of the concentration factors of pollutants in a given sample (Håkanson, 1980), while the pollution load index (PLI) was determined as the geometric mean of the concentration factors (El-Safa et al., 2022).

Evaluating anthropogenic and climatic controls on accumulation rates over the last century

Quantifying land-use/land-cover changes through time

Given that LU correlates with underlying elements causing changes in accumulation rates (Ferland et al., 2014; Kastowski et al., 2011), the record was paired with LU data. Land-use and land-cover changes from 1955 to 2008 were quantified using a combination of high-resolution and medium-resolution remote sensing imagery. High-resolution black and white aerial photographs, which lack spectral data, are the only source of historical LU maps in many areas (Hudak & Wessman, 1997). Imagery used for the analysis included the following: aerial photographs which were available on 8/25/1955 and 10/11/1975; Landsat 4–5 imagery available on 8/24/1982, 7/6/1984, 7/29/1985, and 8/1/1989; and Indian Remote Sensing Satellite (IRS) panchromatic (PAN) imagery available on 6/30/2008. IRS PAN provides a spectral range of 0.50–0.75 m and a spatial resolution of 5.8 m (Patwary & Parvaiz, 2009). Initially, the aerial and IRS PAN images were georeferenced in the UTM WGS84 system. To classify satellite data as forest or non-forest classes, photo-interpreted LU features along with ground-based information were used. LU classification of both types of photos is usually limited to visual interpretation involving digitizing LU classes and estimating their areas (Hudak & Wessman, 1997). Here, to produce LU maps, visual image interpretation was used in ArcGIS 10.4.4 environment.

In addition, cloud-free scenes of Landsat Thematic Mapper (TM) images (with a spatial resolution of 30 × 30 m), taken during the summer, were chosen to best identify the spectral signatures of the various LU classes. Also, for seasonal consistency, it was attempted to select them on similar dates in different years. Applied satellite data were pre-processed, meaning they had been corrected radiometrically and geometrically. Maximum likelihood classifier (Enderle & Weih, 2005) was used to classify images. Specifically, combined supervised and unsupervised classifications were applied in ERDAS 11 to enhance the accuracy of the LU maps (Muhati et al., 2018).

To enable a more accurate understanding of the relationship between land-use change and accumulation rate, the temporal aspect of deforestation was accounted for by developing a new parameter, LUC factor. Thus, LUC factors were developed for each increment by multiplying the percentage of land-use change in that period by the corresponding average slope gradient (in percentage). This approach allowed for the consideration of the extent of deforestation on slopes, taking into account the temporal sequence of deforestation.

Furthermore, areas within the forest exhibiting sparse vegetation were identified on high-resolution images, such as aerial photographs from 1955 and 1975, as well as IRS PAN image from 2008. To accomplish this, a combination of visual image interpretation and unsupervised classifications was employed using ArcGIS 10.4.4 software.

Quantifying climate changes through time

MAP and MAT data were obtained from daily records of two regional meteorological stations, the local Marivan station and the more distant Sanandaj station (Iran Meteorological Organization). Meteorological data was available from the Marivan station beginning in 1992. Sanandaj station is located around 90 km from the study area and has meteorological data available from 1959 to the present (Table S1). Sanandaj experiences a semi-arid cold climate, while Marivan station has a humid cold climate.

A simple yearly correlation model was built for the annual precipitation of Marivan and Sanandaj stations prior to 1992, which was then applied to produce MAP between 1959 and 1992. A trend analysis of climatic parameters was done using the non-parametric Mann–Kendall test and the non-parametric Sen method in Sanandaj and Marivan meteorological stations (Salmi, 2002). As MAPs were only available since 1959 (depth ~ 33 cm), we assumed MAPs in previous years (n = 4) to be equal to the average MAP (1139.0 mm). Ultimately, Pearson’s correlation and principal component analysis (PCA) were conducted to explore potential relationships between sediment parameters, climatic and land-use change factors, and accumulation rates using R-Studio. Additionally, flood reports were collected from local and national news sources via online platforms; however, they were not used in data interpretation due to concerns about their comprehensiveness.

Upscaling MAR and OCAR to the Zarivar Lake

Although the primary objective of this study was to investigate the anthropogenic and climatic impacts on sediment accumulation rates in a small sub-watershed (i.e., Z3) connected to Zarivar Lake, the accumulation rates from this study were also used in a two-compartment model to improve estimates for in the entire lake and inform whole-large management strategies. The estimation of whole lake MAR was done by combining observations from the core collected for this study in the marshlands with observations from a core collected in the center of the lake, as reported by Torabi Kachoosangi et al. (2020). Core locations are shown in Fig. 1. In this approach, it was assumed that the center-lake core was representative of the non-vegetated lake regions (8.8 km2) and that the marshland core was representative of the marshland areas (8.9 km2). As accumulation rates within marshlands and center lake might differ due to lower velocity and higher deposition rates as well as higher nutrient cycling in marshlands compared to the lake center.

The approximate calculations were conducted for the whole lake, including sedimentation rate, sediment production per unit area of the lake watershed, the amount of sediment stored behind the check dams of the Z3 sub-watershed over an 8-year period, and the OCAR for recent sediments. To achieve these, a method was employed to upscale the burial rates of single sample data to estimate whole lake burial rates. This method involved computing a focusing factor, as described by Engstrom (2005). The focusing factor is defined as the ratio of the observed 210Pbex accumulation rate in the core (F[210Pbsite]) to the known atmospheric flux of 210Pbex in the region (F[210Pbatm]; Rippey & Douglas, 2004), which is approximately 216 Bq m−2 year−1 in the region under investigation (Khodadadi et al., 2023). The sediment production per unit area of the Zarivar Lake watershed was estimated by dividing the annual mass accumulation in both the waterbody and marshland by the total area of the watershed. Additionally, using an average net soil erosion rate and sediment production per unit area of Z3 watershed, we estimated the sediment delivery ratio (SDR) and the volume of sediment trapped behind the check dams over an 8-year period. The former was approximated by dividing the annual sediment production by the erosion amount, while the latter was determined by subtracting the difference between the erosion amount and sediment production over the 8-year period. Similarly, whole lake OCAR was computed based on surface sediments (ca. 10 years of accumulation), following the method described by Hobbs et al. (2013).

Results

Geochronology and accumulation rates

The 137Cs profile indicated two distinct peaks at depths of 13 and 30 cm, corresponding to 1986 and 1963 peaks, respectively (Fig. 2a), while 210Pbex profile represented an extremely non-uniform depth profile. The irregular 210Pbex profile indicated that there may have been significant changes in the sedimentation rate over the past century. For such variable 210Pbex depth profile, the CRS model is normally used for the stratigraphic chronology (Fig. 2b). However, a small discrepancy between the CRS model dates and the 1986 and 1963 depths determined from the 137Cs records was observed. Therefore, an approach, known as composite modeling, was applied to calculate the corrected 210Pbex dates for the Zarivar Lake core (Fig. 2c; for more details on composite models, see Appleby, 2008). This was in accordance with elsewhere, e.g., in Blelham Tarn Lake, UK (Appleby, 2001, 2008) and Żabińskie Lake, Poland (Tylmann et al., 2016). This can be due to various causes such as flood events, sediment slumps, turbidity currents, and major land-use and land-management changes (Appleby, 2001). Here, the most probable reason for this discrepancy was the substantial reduction of sediment supply following the construction of the check dams in the river network since 2008. In general, in the case of an irregular 210Pbex record, the CRS model is normally applied in a piecewise way in various zones of the core in which the chronostratigraphic dates are deployed as reference points, an approach that yields corrected 210Pbex dates (Appleby, 2008; Tylmann et al., 2016). Thus, composite CRS model was applied to calculate the corrected 210Pbex dates and MAR. Based on this method, we sub-divided the record into two sections, i.e., below 13 cm (pre-1986) and above 13 cm (post-1986), and calculated age and MAR for all intervals in each section.Fig. 2 Depth profile of 137Cs (with analytical error bars (a)); 210Pb, 214Pb, and 210Pbex (b); 210Pbex (with analytical error bars (c)); and composite CRS model chronology with the 1986 and 1963 depths indicated by the 137Cs records (d) (the record was sub-divided into two sections at 14 cm, i.e., 1986). Sedimentation rates (g cm−2 year−1) and errors from the CRS model are also shown

Unusually high SARs were observed at the bottom of the core (> 50 cm; Figs. 2b and 3), which could be attributed to reaching the end of the 210Pbex profile, uncertainty of 210Pb measurement, or rainfalls with a 100-year return period. As calculating the ages in depth with certainty in depths > 50 cm proved challenging without applying independent tracers, e.g., 14C (Tylmann et al., 2016), this section was excluded from further calculations.Fig. 3 Estimated SARs (mm year−1) using composite CRS model in the Zarivar Lake. Dashed horizontal lines bracket a period with high SARs. Average SAR (dashed green line), estimated MAPs (solid blue line), and reported floods in the region and country ( ×) through time are shown as well. Note: Flood reports were collected from local and national news sources through online platforms and were not used in the interpretation of the data

Sedimentation rate varied considerably through time (Fig. 3), and a period of significantly higher-than-average SARs was identified. Subsequently, we sub-divided the core into three sections based on SAR, i.e., 0–15 cm (upper), 15–25 cm (middle), and > 25 cm (lower). ANOVA showed a significant difference in MAR across sections (p < 0.01; Fig. 4a). The mean dry bulk density stood at 0.99 ± 0.19 g cm−3, being the highest in the upper 15 cm (1.08 ± 0.26 g cm−3), mid-section at 1.04 ± 0.09 g cm−3, and slightly lower at 0.91 ± 0.18 g cm−3 in the bottom section. The average CRS-derived MAR was 6498 ± 2475 g m−2 year−1 (SAR = 6.9 mm year−1).Fig. 4 The range of MAR (a), OCAR (b), and PPAR (c) in upper (0–15 cm), middle (15–25 cm), lower (25–50 cm) sections of the sediment core. The Kruskal–Wallis test indicated a significant difference among the groups for MAR (χ² (chi-square) = 6.70, p = 0.0350), OCAR (χ² = 10.01, p = 0.0067), and PPAR (χ² = 15.58, p = 0.0004). Pairwise comparisons were conducted using Dunn’s post hoc test with Bonferroni correction. Different letters indicate significant differences at the p-level < 0.05

The mean OC content and PP for all increments were 3.2 ± 0.9% and 1035 ± 108 mg kg−1, respectively (Fig. 5a). OCAR and PPAR estimates were calculated based on MAR derived from 210Pbex-composite CRS model and representing an 80-year time span. Mean OCAR and PPAR in all intervals stood at 205 ± 85 g m−2 year−1 and 8.9 ± 3.3 g m−2 year−1, respectively (Fig. 5b). OCAR was significantly lower in the upper zone of the core (Fig. 4b). The PPAR in the middle section was significantly higher than in the two other sections (Fig. 4c).Fig. 5 Depth profile of various parameters in the Zarivar Lake, including total organic carbon (OC; %), particulate phosphorus (PP; mg kg−1), and δ13C (Suess effect corrected; ‰) contents (a) and depth profiles of mass accumulation rate (MAR; kg m−2 year−1), organic carbon accumulation rate (OCAR; g m−2 year−1), and particulate phosphorus accumulation rate (PPAR; g m−2 year−1) based on the composite 210Pbex-CRS model (b). Brown arrows show the general trend of the parameters with depth. Dashed lines bracket a period with distinct high MARs

Elemental geochemistry

The ICV ranged from 0.83 to 1.24 and the SiO2/Al2O3 ratio varied between 2.65 and 3.34, highlighting a low level of compositional maturity across various sections of the core (Table 1). The CIA values of 69.37, 70.04, and 73.55 for the upper, middle, and lower core sections, respectively, revealed moderate weathering in the two top sections but an intense weathering condition for sediment in the bottom section (Table 1). Zarivar Lake sediments demonstrated low U/Th varying from 0.24 to 0.40 and V/Cr ranging from 0.96 to 1.23, suggesting an oxic depositional environment (Jones & Manning, 1994; Table 1). Ni/Co ratios were at 4.64, 4.70, and 5.06 in the upper, middle, and lower sections, respectively, indicating an oxic depositional environment in the upper parts but a sub-oxic one in the lower zone.

The multi-element pollution indices, i.e., Dc and PLI, were computed based on the concentrations of Al, Cu, Co, Cr, Fe, Mn, Ni, P, V, and Zn (Table S2). Dc values ranged from 11.5 to 14.1, indicating a moderate contamination level (8 < Dc < 16; Håkanson, 1980), while PLI values varied between 1.0 and 1.3 in the samples, placing them in the polluted group (PLI > 1; Håkanson, 1980). Both indices exhibited significant differences across different core sections, with the lowest values observed at the bottom core section, followed by the top section. This suggests that elevated levels of contamination can enter the lake at higher MARs. This is also reflected in the positive significant correlations (at p > 0.05) observed between Dc and PLI and MAR (Fig. S1).

Quantifying land-use/land-cover change over time

Aerial photos and satellite data illustrated that compared to 1955, the net forest cover had declined by 6 and 18% in 1975 and 1984, respectively, decreasing from 75% of the watershed area 1955 to 57% in 1984 (Fig. 6a). It is notable that 12% of forests been cleared between 1975 and 1984 and in addition, based on the digital elevation model (DEM), after 1975, deforestation took place on steeper slopes (< 10% vs 20–30%) (Fig. 6a). In contrast, from 1984 to 2008, major LUCs did not occur except for the quite limited expansion of developed areas.Fig. 6 3D images of temporal sequence land-use classifications and the percentage of forest converted to agriculture in the watershed between 1955 and 2008 (a) and expansion of areas with sparse vegetation in forest in 1975 (pink) and 2008 (red) (b). Note that investigating satellite data for 1984, 1985, 1989, and 2008 revealed no significant changes in land use within the watershed

In addition, the visual interpretation of high-resolution data, i.e., the aerial photo from 1975 and the 2008 IRS PAN imagery, showed that the areas with sparse vegetation had almost doubled from 0.12 to 0.22 km2 during the period between 1975 and 2008, with most of these areas being situated on southern slopes (Fig. 6b). This could be due to both human interferences like grazing and wood harvest (cutting branches from trees) as well as climate change.

Quantifying climate change over time

The trend analysis of climatic parameters using the Mann–Kendall test and Sen’s slope estimates (Salmi, 2002) showed a clear decreasing trend in MAP in the study area (1959–2017; Fig. S2). However, despite an upward trend in average minimum and maximum temperature in the regional data (Sanandaj station), there was no trend detected in the watershed temperature data (Marivan station; Fig. S2). This could be due to the presence of the lake preventing temperature variabilities between 1992 and 2017. Consequently, the impacts of temperature on accumulation rates were considered to be marginal in the region.

Discussion

Interpreting erosional history through changes of sediment geochemistry indices through time

Elemental geochemistry can help assess sediment maturity, chemical weathering, and depositional environments in a core. A matured sediment is composed mainly of stable minerals such as quartz, so it is unlikely to undergo further weathering processes due to the absence of unstable minerals (Pettijohn, 1975). Based on the ICV (> 0.83) and the SiO2/Al2O3 ratio (> 2.65) values, a low compositional maturity was proved in the core (Table 1). High ratios of SiO2/Al2O3 usually correspond with high quartz content, indicating chemically mature sediments (Mahamuda, 2015). Owing to the dominance of expandable clays (smectite group or 2:1) measured by semi-quantitative X-ray diffraction (XRD) of the clay fraction of soils in the region (particles < 2 μm; Khodadadi et al., 2024), such SiO2/Al2O3 ratios were expected in the sediments (particles < 63 μm). The CIA values revealed moderate weathering in the two top sections but an intense weathering condition for the bottom section (Table 1). This suggests older contributions were from forest topsoil and more recent inputs eroded from agriculture or deeper in the soil profile, e.g., incisions.

Anthropogenic and climatic controls on SAR

The average SAR of the lake was estimated to be at 6.90 ± 3.45 mm year−1, varying between 2.02 and 17.73 mm year−1 (Fig. 3). This value was close to the average sedimentation rate of 5.6 mm year−1 documented by Torabi Kachoosangi et al., (2020; Fig. 1), but the small difference can be due to differences in sampling locations within Zarivar Lake. Although sedimentation rates in different lakes can be influenced by many factors, making comparisons challenging, we present findings from several studies conducted in similar environmental settings. For instance, Morellón et al. (2011) studied climate changes and human activities recorded in the sediments of a shallow lake, draining a small catchment of 1.07 km2, in Spain. The catchment’s lowlands are under farming lands dedicated to barley cultivation, while higher elevation areas are mostly occupied by scrublands and oak forests. Two short cores were collected and sub-sampled every 1 cm for 137Cs and 210Pb dating. The sedimentation rate over the past 50 years was estimated at 3.41 mm year−1, which is lower than our estimated value. This difference can be attributed to the larger area and unsustainable management practices in Z3 sub-watershed. In a study, Valero-Garcés et al. (2006) examined human impact on a Mediterranean lake in Spain with a watershed of 900 km2 and a MAP of 538 mm. Soil erosion intensified in the mid-twentieth century due to farm machinery introduction. In the early 1980s, the lake became a protected area with restricted agricultural practices and implemented conservation measures. Two sediment cores were retrieved for 210Pb and 137Cs dating, revealing a mean sedimentation rate of 15 mm year⁻1 over the past half-century. This rate was substantially higher than our 80-year average, which may be due to the extremely large area of the watershed compared to Z3. In another study, Barreiro-Lostres et al. (2015) investigated human-climate interactions on SAR in a Mediterranean mountain lake in Spain, with a basin area of 11 km2 and a MAP of 956 mm. The uplands are dominated by pine trees, while cereal fields are found at lower elevations. Using 210Pb and 137Cs dating, they determined a SAR of 10 mm year⁻1 for the past 60 years. Although higher than our predicted value, it is notable that this sediment comes from a watershed four times larger than our study site, indicating unsustainable management practices in our studied watershed.

To gain a comprehensive understanding of the anthropogenic and climatic impacts on SAR, two specific time periods, 1975–1988 and 2008–2016, characterized by the highest levels of human interference in the watershed, were examined. Between 2008 and 2016, the amount of sediment entering the lake substantially reduced. It can be inferred that check dams managed to store almost all mobilized sediments in the watershed, demonstrating the effectiveness of the measures taken. This is supported by field observations of substantial sediment trapped behind the numerous check dams, which resulted in the majority of them being almost full of sediment in merely 8 years (see Fig. S4). Khodadadi et al. (2022) evaluated the main sediment sources in the sub-watershed before (1995–2009) and after (2009–2016) introducing the check dams to the sub-watershed using geochemical fingerprints. The results showed a significant shift in the main sediment sources from channel banks to uncultivated subsoil following the establishment of the check dams, which verifies their effectiveness in stream bank stabilization, at least for the duration of their infill stage. The efficiency of check dams in decreasing catchment sediment yield has been confirmed elsewhere (e.g., Boix-Fayos et al., 2008; Tiessen et al., 2011). However, their impact is known to be short-lived (Tiessen et al., 2011), and further investigation on their effectiveness is thus recommended in the sub-watershed in the near future, particularly due to the observed dramatic decline in storage capacity.

Between the late 1970s and early 1980s, the MAR was significantly higher than its average at 7.94 ± 3.12 kg m−2 year−1 (Fig. 4a). Scrutinizing LULC changes through time (Fig. 6a) and estimated MAPs (Fig. 3) in the watershed suggested that this may have been the result of combined human activities and climatic events (Figs. 4a and 5a). As previously noted, most of the LUCs (12%) in the watershed occurred between 1975 and 1984. Extensive deforestation after Iran’s revolution (1979) and the war between Iraq and Iran (1980–1988) were likely to have played a role. More precisely, between 1978 and 1980, the area experienced extensive deforestation subsequent to the revolution and probably caused by lack of law reinforcement, coinciding with higher MAPs (more than double). In the period between 1981 and 1983, limited military activities began in mountainous areas (building ditches) and the population started to increase (from around 350 persons in 1976 to more than 750 persons in 1986; Statistical Center of Iran), happening along with the higher MAPs. Together, these factors seem to have caused excessive amounts of sediment to enter the lake.

To elucidate the impact of LUC and MAP on SAR, correlation coefficients were established between accumulation rates and sediment properties, MAP, and land-use change factor up to the depth of 40 cm, which corresponds to the maximum depth of measurement for OC (1955; Fig. 7a). While a significant correlation was observed between SAR and LUC factors (r = 0.53, p < 0.01), there was no significant relationship between SAR and MAP (r = 0.32). This suggests that human activities had more pronounced effects on SAR, while climatic factors amplified anthropogenic impacts. Notable is that in the absence of extensive LUCs but high amounts of estimated MAP, similar to the events occurring in 1968, SAR increased only slightly (Fig. 3). All in all, it can be claimed that human-induced activities have resulted in varying levels of SAR, ranging from values that are twice the average rates to insignificant amounts. Olley and Wasson (2003) showed how small climate variabilities and a range of human activities combined had impacted sediment flux in a watershed in Australia over a period of 180 years, resulting in a non-linear relationship between SAR and rainfall, runoff, and sediment transport capacity. Having reviewed numerous case studies, Foster et al. (2003) recognized the history of LUCs as the primary element influencing sediment flux in the USA. It should be noted that accumulation rate and flux are not the same: the former refers to the amount of materials or nutrients that build up in a specific area over time, while the latter refers to the rate at which materials or nutrients are transferred into or out of a system.Fig. 7 Pearson correlation coefficient (a) and PCA1 and PCA2 (b) of sediment properties, MAP, land-use change factor (percentage of land-use change multiplied by its average slope), and accumulation rates up to 40 cm (n = 27). Non-colored cells in 7a indicate that correlation was insignificant at p < 0.05. Note: Si.Al = SiO2/Al2O3 ratio and LUC.Slope = LUC factor

Although the accumulation rates were highly sensitive to both anthropogenic and climatic factors that occurred in the sub-watershed, some of the changes in the lake seem not to impact accretion rates. For instance, the diversion dam constructed on a river, which drains agricultural lands containing high amounts of sediment and nutrients, did not exert any considerable influence on accumulation rates in 2004 or thereabouts. This may be due to the core location, which was close to the studied sub-watershed in the west of the lake.

Anthropogenic and climatic controls on OCAR and PPAR

The mean OCAR in the entire core was estimated to be 205 ± 85 g m−2 year−1, which is 20-fold greater than the reported global means of 4 to 14 g m−2 year−1 (Cole et al., 2007; Tranvik et al., 2009). However, agriculturally impacted lakes, like Zarivar Lake, are subject to receiving high amounts of nutrients, which in turn stimulates autochthonous primary production, as well as sediments, caused by higher susceptibility to erosion of tilled soils. Consequently, the OCAR in these lakes is likely to be higher than lakes in catchments with less agricultural activity (Downing et al., 2008). There was previously no report on OCARs for Zarivar Lake; thus, we compared our value with a study in a similar environmental setting. Dominik et al. (2023) studied the organic carbon burial rate in a Mediterranean lagoon in Italy. Using four cores and the 210Pb dating technique, they estimated a mean SAR of 3.6 mm year⁻1 and OCAR of 49 g m⁻2 year⁻1. While the SAR is somewhat close to our estimated average, the OCAR is substantially smaller than our predicted value, highlighting unsustainable management in our study area.

The depth profile of the OC indicates that, except for the lowest increment depth (OC = 6.4%), OC fluctuated between 1.9 and 4.3% in the rest of the depth intervals. OCAR at upper, middle, and lower sections were 140 ± 30, 220 ± 60, and 240 ± 100 g m−2 year−1, respectively (Fig. 5b). Torabi Kachoosangi et al. (2020) also reported an increasing OC trend with depth, particularly at depths > 26 cm, in a core in the center of Zarivar Lake, which is consistent with our findings. Several factors can cause such an increasing trend in OC with depth. Firstly, the coefficient of variation (CV) of OCAR in the upper section of the core was significantly lower than in the other zones. LULC changes greatly affect OC export from terrestrial ecosystems (Smith et al., 1999), meaning that OCAR is likely to be relatively constant when a significant LULC change has not occurred (Anderson et al., 2009; Kortelainen et al., 2004). This finding is consistent with our results, as the extent of LUCs was minimal after 1984. Secondly, as U/Th and V/Cr ratios showed that the record was in oxic conditions in all intervals (Table 1), a dramatic reduction in OC in the core can be a reflection of lower OC contents of the sources, i.e., soils of the region. To illustrate, Khodadadi et al. (2023) reported that following deforestation (in 1981) on a hillslope in the sub-watershed, soil loss increased by five times, resulting in a 10-cm net soil loss, and soil organic carbon stock (SOCS) decreased by 30% (14 Mg C ha−1) within a course of 35 years. As such, 8 Mg C ha−1 of SOCS was removed by erosion, and the remaining seemed to be lost via emissions over the period in question. Thus, progressive conversion of oak forests to cropland along with an increase in areas with sparse vegetation (Fig. 7b) typically characterized by a lower OC content (Fig. S3) might have reduced OC content in the topsoils of the region through time. Together, it seemed that a decline in forest area and density, climate change (lower MAPs), and accelerated soil erosion had led to a lower OC in soils of the region since the 1980s. This hypothesis, however, may need further testing to ensure that sub-oxic conditions in lower layers, shown by the other geochemical index, i.e., Ni/Co ratio, have not impacted the organic material decomposition since it is known that organic material in anaerobic sediments tends to decompose more slowly (Dong et al., 2012; Ferland et al., 2014).

Another reason for lower OC content and decreases in OCARs in younger sediments could be the result of a change in OC sources through time, e.g., probably from forests with high OC to mainly subsoil sources with lower OC (Fig. S3). These changes of sediment sources over time were also inferred by a positive significant correlation coefficient between CIA and SiO2/Al2O3 ratio and OCAR (Fig. 7a), implying a change in sediment source over time from more mature (under forest) to younger sources (under cultivation). Furthermore, Lin et al. (2021) also found a notable decline in OCAR after the 1980s. This decrease was hypothesized to have been caused by increased OM decomposition and OC mineralization in the water column induced by increasing temperatures and lake-water thermal stratification. Although here, the temperature change was assumed to be insignificant over recent decades, it might have played a role as well.

While SAR and PPAR showed a significant correlation with LUC factors but an insignificant one with MAP, OCAR displayed non-significant relationships with both parameters (Fig. 7a), indicating that controlling factors of the SAR and PPAR drastically differ from those of OCAR (Fig. 7b). Nonetheless, many studies on lacustrine OC burial rates have also recognized LU practices as the primary contributor to observed patterns (Einola et al., 2011; Ferland et al., 2014; Kastowski et al., 2011; Kortelainen et al., 2004; Sanders et al., 2017). Anderson et al. (2013) emphasized the role of anthropogenic activities on OCAR in lake sediments, highlighting their importance over climatic factors. However, Kortelainen et al. (2013) claimed that climate was a more important element impacting OC in aquatic systems rather than land use. Here, it seems that OCAR is impacted by LUC and MAP factors to some extent (r ~ 0.2); however, other factors may play a role as well, e.g., change in OC sources and their OC content over time.

PPAR followed the same general pattern of the MAR with depth (Fig. 5b); i.e., PPAR was significantly higher in the mid-section, at 12.0 ± 2.8 g m−2 year−1, than in the upper and lower sections, at 6.4 ± 1.8 and 7.4 ± 1.9 g m−2 year−1, respectively (Fig. 5c). Significant correlations were observed between PPAR and MAR (0.99) as well as between PPAR and OCAR (0.52), the first of which considerably exceeds the second correlation coefficient. This implies that the process impacted inorganic particulate matter and PP transport and fate may have been similar, whereas some other processes may have impacted OCAR in the sediment core. However, it is important to note that the investigation of these processes and their changes over time is beyond the scope of this study.

There are distinctive isotope signatures in C3 vs. C4 plants (i.e., cool weather grasses, trees, and shrubs vs warm weather grasses, cereals, and crops) owing to a marked difference in 13C discrimination in their photosynthetic pathways (Farquhar et al., 1989). C3 and C4 plants and algae have a mean δ13C value of − 27‰ (− 32 to − 20‰), − 12‰ (− 17 to − 9‰), and − 22‰ (− 30 to − 20 ‰; France, 1996), respectively. The mean δ13C for all intervals was − 28.53 ± 0.52 ‰, and there was an incremental decrease with depth (Fig. 6a), signifying that C3 plants (e.g., forest) have covered the region and no abrupt ecological shifts (e.g., C3 to C4) have occurred. Sanders et al., (2014, 2017) reported a gradual decrease in δ13C with depth in mangrove sediments, which related to a growth in algal production. The TDP measurements at two points in the lake at 1.21 and 1.02 mg l−1 were in line with reported values by Ebrahimi Mohammadi et al. (2012). Thus, the lake had eutrophic conditions over the last few decades, leading to an acceleration in its primary productivity. This often includes higher rates of algal growth. Therefore, here, also higher values of δ13C in younger sediments may also have been caused by algal growth in response to larger amounts of nutrients entering the lake. Another possible explanation for SOC δ13C variations with depth could be that the vegetation found in the past might have been slightly different from existing flora (Mermut & Acton, 1984). Overall, δ13C imposes some limitations to reliably describing OC dynamics in ecosystems not having experienced a shift in vegetation from C3 to C4 plants (Trumbore, 2009).

Estimating whole-lake accumulation rates

To assess the average whole-lake sedimentation rates, the focusing factor was calculated using the observed 210Pbex accumulation rate in the core to the known atmospheric flux of 210Pbex in the region. The sediment focusing factors were at ca. 1.0 and 1.2 for central (Torabi Kachoosangi et al., 2020) and marshland (this study) cores, respectively. The values were in agreement with values reported by Rippey and Douglas (2004) for lake sediments in the UK and Ireland, with focusing factors varying from 0.37 to 7.77. The whole-lake average sedimentation rate was estimated at approximately 7 mm year−1 and the annual sedimentation amount was 100 Gg year−1. The amount of sediment production per unit area of the Zarivar Lake watershed (296.4 km2) was computed to be approximately 3.4 Mg ha−1 year−1. Further, using the two cores, the average whole lake OCAR for near-surface sediments was estimated at 300 g m−2 year−1. From this, an estimated 5 Gg OC year−1 was accumulated in the lake during the past 10 years.

As for the Z3 sub-watershed, using an average net soil erosion rate of 12.6 Mg ha−1 year−1 (unpublished data), the sub-watershed area, and the sediment production per unit area of the Zarivar Lake watershed (3.4 Mg ha−1 year−1), the SDR of the sub-watershed would be around 27%. The magnitude of the trapped sediments behind the check dams during 8 years was estimated to be 8 Gg. This calculation highlights the effectiveness of check dams at this site in reducing sediment load to the lake. Overall, these estimations can provide valuable insights which could inform lake management strategies.

However, several limitations should be acknowledged in our study. To begin with, while upper increments could be precisely dated using 210Pbex, dating bottom core increments requires 14C measurements. Furthermore, the study could gain insight from investigating the impacts of shifts in sediment sources over time to elucidate variations in MAR, OCAR, and PPAR. Moreover, the study could benefit from additional measurements, such as N, δ15N, and TP, to further illuminate the eutrophic condition of the lake over the past decades (Zan et al., 2012). Additionally, although we calculated rough estimates of entire lake accretions, enhancing the precision of sedimentation estimations across the entire lake watershed would require collecting more cores from both the marshlands and the lake. Lastly, applying our methodology to broader spatial scales can help evaluate the impacts of land-use and climate changes in each sub-watershed and across the entire lake watershed, contributing to a more holistic understanding of the lake ecosystem dynamics.

Conclusion

Over the last century, Zarivar Lake has been subjected to anthropogenic interferences, including intensive LUCs and the introduction of some soil conservation practices. To understand the extent of the effects that these human activities and climate change on the magnitude of sediment, OC, and nutrients accretions through time, a short sediment core was collected within the lake. LUCs were quantified using aerial photos in 1955 and 1975; Landsat 4–5 in 1982, 1984, 1985, and 1989; and IRS PAN in 2008. This study revealed that MAR, OCAR, and PPAR during the late 1970s and early 1980s were significantly higher than average, which is thought to result from massive human interferences cooccurring with extreme climatic events. In contrast, human intervention in the form of check dams proved effective conservation measures, storing almost all sediment in the watershed between 2008 and 2016.

Together, these results illustrate that human activities caused accumulation rates to fluctuate from substantial amounts, with values twice the averages, to insignificant amounts following the installation of check dams, during the studied time window. Differentiating anthropogenic from climatic effects can be immensely complex at the basin scale over historical times, particularly in watersheds subjected to intense human interferences, and was not entirely possible in this study. However, the results suggested that human activities had more pronounced effects on MAR, OCAR, and PPAR while climatic factors amplified anthropogenic impacts. Overall, our study showed that applying a composite CRS model is recommended when the record experiences a hiatus in sedimentation for a while. The findings clearly demonstrated that land-use changes on steep slopes are responsible for the elevated accumulation rates of sediment, OC, and PP in lakes, particularly when these changes coincide with extreme weather events. Furthermore, while land-use and climate changes can accelerate the influx of fine sediments and associated nutrients into inland lakes, watershed conservation practices could effectively reduce the supply of fine particles. Our findings clearly demonstrated that applied conservation measures were able to mitigate the adverse impacts of climate variabilities over the past decades, i.e., the decreasing trend of MAP. However, it is important to note that check dams are considered short-term conservation measures. This means that the lake’s future accumulation rates might fluctuate significantly, particularly if new conversion measures are not placed and LUC, soil erosion, and lake eutrophication are not regulated. Therefore, it is imperative for local policymakers to implement additional check dams in the Z3 sub-watershed. These findings also provide a basis for global policymakers to enforce restrictions on deforestation on steep slopes and monitor the long-term effectiveness of implemented conservation measures.

The applied methodology has proven effective in evaluating the impacts of LUC, climate variability, and the effectiveness of conservation measures in ungauged watersheds linked to reservoirs. However, further research is essential to fully understand the impacts of LUC and the effectiveness of conservation methods in inland lakes under extreme weather events and prolonged climate shifts to ensure the long-term health of water bodies.

Supplementary Information

Below is the link to the electronic supplementary material.Supplementary file1 (DOCX 2.05 MB)

Acknowledgements

This study has been finalized to support the IAEA Coordinated Research Project (CRP) on “Nuclear techniques for a better understanding of the impact of climate change on soil erosion in upland agroecosystems” (D1.50.17). The authors would like to thank Dr. Peter Strauss, Dr. Moncef Benmansour, and Dr. Amy T. Hansen for their valuable comments and feedback in preparing this manuscript.

Author contribution

The authors have made significant contributions to the preparation, writing, and revision of the manuscript. Maral Khodadadi contributed to funding acquisition, data collection and analyses, software utilization, interpretation of data, and drafting of the work. Max Gibbs and William H. Blake were responsible for overseeing all aspects of the work, ensuring that any questions pertaining to the accuracy or integrity of the content were properly addressed and resolved during manuscript revision. Andrew Swales and Arsenio Toloza facilitated the analysis process. All authors have thoroughly reviewed and endorsed the final manuscript.

Funding

This study has been finalized to support the IAEA Coordinated Research Project (CRP) on “Nuclear techniques for a better understanding of the impact of climate change on soil erosion in upland agroecosystems” (D1.50.17).

Data availability

Data availability The datasets used and analyzed during the current study are available from the authors upon reasonable request.

Declarations

Competing interests

The authors declare no competing interests.

Highlights

• The construction of check dams reduced lake sedimentation substantially

• Organic carbon burial rate was 20-fold greater than the global mean for lakes

• Human activity has had a more significant impact on accumulation rates than climate variability

• Climatic factors amplified anthropogenic impacts.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

Anderson NJ D’andrea W Fritz SC Holocene carbon burial by lakes in SW Greenland Global Change Biology 2009 15 11 2590 2598 10.1111/j.1365-2486.2009.01942.x
Anderson, N. J., D’andrea, W., & Fritz, S. C. (2009). Holocene carbon burial by lakes in SW Greenland. Global Change Biology, 15(11), 2590–2598. 10.1111/j.1365-2486.2009.01942.x10.1111/j.1365-2486.2009.01942.x
Anderson NJ Dietz RD Engstrom DR Land-use change, not climate, controls organic carbon burial in lakes Proceedings of the Royal Society b: Biological Sciences 2013 280 1769 2013 1278 10.1098/rspb.2013.1278
Anderson, N. J., Dietz, R. D., & Engstrom, D. R. (2013). Land-use change, not climate, controls organic carbon burial in lakes. Proceedings of the Royal Society b: Biological Sciences, 280(1769), 2013–1278. 10.1098/rspb.2013.127810.1098/rspb.2013.1278
Appleby PG Last WM Smol JP Chronostratigraphic techniques in recent sediments Tracking environmental change using lake sediments volume 1: Basin analysis, coring, and chronological techniques 2001 Kluwer Academic 171 203
Appleby, P. G. (2001). Chronostratigraphic techniques in recent sediments. In W. M. Last & J. P. Smol (Eds.), Tracking environmental change using lake sediments volume 1: Basin analysis, coring, and chronological techniques (pp. 171–203). Kluwer Academic.
Appleby PG Three decades of dating recent sediments by fallout radionuclides: A review Holocene 2008 18 1 83 93 10.1177/0959683607085598
Appleby, P. G. (2008). Three decades of dating recent sediments by fallout radionuclides: A review. Holocene, 18(1), 83–93. 10.1177/095968360708559810.1177/0959683607085598
Appleby PG Oldfield F The calculation of lead-210 dates assuming a constant rate of supply of unsupported 210Pb to the sediment CATENA 1978 5 1 1 8 10.1016/S0341-8162(78)80002-2
Appleby, P. G., & Oldfield, F. (1978). The calculation of lead-210 dates assuming a constant rate of supply of unsupported 210Pb to the sediment. CATENA, 5(1), 1–8. 10.1016/S0341-8162(78)80002-210.1016/S0341-8162(78)80002-2
Appleby, P. G., & Oldfield, F. (1992). Applications of lead-210 to sedimentation studies. In Uranium-series disequilibrium: Applications to Earth, Marine, and Environmental Sciences. 2. ed.
Armstrong-Altrin JS Machain-Castillo ML Rosales-Hoz L Carranza-Edwards A Sanchez-Cabeza JA Ruíz-Fernández AC Provenance and depositional history of continental slope sediments in the Southwestern Gulf of Mexico unraveled by geochemical analysis Continental Shelf Research 2015 95 15 26 10.1016/j.csr.2015.01.003
Armstrong-Altrin, J. S., Machain-Castillo, M. L., Rosales-Hoz, L., Carranza-Edwards, A., Sanchez-Cabeza, J. A., & Ruíz-Fernández, A. C. (2015). Provenance and depositional history of continental slope sediments in the Southwestern Gulf of Mexico unraveled by geochemical analysis. Continental Shelf Research, 95, 15–26. 10.1016/j.csr.2015.01.00310.1016/j.csr.2015.01.003
Bakkiaraj D Nagendra R Nagarajan R Armstrong-Altrin JS Geochemistry of sandstones from the Upper Cretaceous Sillakkudi Formation, Cauvery Basin, southern India: Implication for provenance Journal of the Geological Society of India 2010 76 453 467 10.1007/s12594-010-0128-3
Bakkiaraj, D., Nagendra, R., Nagarajan, R., & Armstrong-Altrin, J. S. (2010). Geochemistry of sandstones from the Upper Cretaceous Sillakkudi Formation, Cauvery Basin, southern India: Implication for provenance. Journal of the Geological Society of India, 76, 453–467.10.1007/s12594-010-0128-3
Barreiro-Lostres F Brown E Moreno A Morellón M Abbott M Hillman A Giralt S Valero-Garcés B Sediment delivery and lake dynamics in a Mediterranean mountain watershed: Human-climate interactions during the last millennium (El Tobar Lake record, Iberian Range, Spain) Science of the Total Environment 2015 533 506 519 10.1016/j.scitotenv.2015.06.123 26186465
Barreiro-Lostres, F., Brown, E., Moreno, A., Morellón, M., Abbott, M., Hillman, A., Giralt, S., & Valero-Garcés, B. (2015). Sediment delivery and lake dynamics in a Mediterranean mountain watershed: Human-climate interactions during the last millennium (El Tobar Lake record, Iberian Range, Spain). Science of the Total Environment, 533, 506–519. 10.1016/j.scitotenv.2015.06.12326186465 10.1016/j.scitotenv.2015.06.123
Boix-Fayos C de Vente J Martínez-Mena M Barberá GG Castillo V The impact of land use change and check-dams on catchment sediment yield Hydrological Processes: An International Journal 2008 22 25 4922 4935 10.1002/hyp.7115
Boix-Fayos, C., de Vente, J., Martínez-Mena, M., Barberá, G. G., & Castillo, V. (2008). The impact of land use change and check-dams on catchment sediment yield. Hydrological Processes: An International Journal, 22(25), 4922–4935.10.1002/hyp.7115
Breithaupt JL Smoak JM Smith TJ III Sanders CJ Temporal variability of carbon and nutrient burial, sediment accretion, and mass accumulation over the past century in a carbonate platform mangrove forest of the Florida Everglades Journal of Geophysical Research: Biogeosciences 2014 119 10 2032 2048 10.1002/2014JG00271
Breithaupt, J. L., Smoak, J. M., Smith, T. J., III., & Sanders, C. J. (2014). Temporal variability of carbon and nutrient burial, sediment accretion, and mass accumulation over the past century in a carbonate platform mangrove forest of the Florida Everglades. Journal of Geophysical Research: Biogeosciences, 119(10), 2032–2048. 10.1002/2014JG0027110.1002/2014JG00271
Canuel, E. A., Brush, G. S., Cronin, T. M., Lockwood, R., & Zimmerman, A. R. (2017). Paleoecology studies in Chesapeake Bay: A model system for understanding interactions between climate, anthropogenic activities and the environment. In Applications of Paleoenvironmental Techniques in Estuarine Studies (495–527). Springer, Dordrecht. 10.1007/978-94-024-0990-1_20
Chen, B., Zheng, Z., Huang, K., Zheng, Y., Zhang, G., Zhang, Q., & Huang, X. (2014). Radionuclide dating of recent sediment and the validation of pollen-environment reconstruction in a small watershed reservoir in Southeastern China. Catena, 115, 29–38. 10.1016/j.catena.2013.11.010
Cole JJ Prairie YT Caraco NF McDowell WH Tranvik LJ Striegl RG Duarte CM Kortelainen P Downing JA Middelburg JJ Melack J Plumbing the global carbon cycle: Integrating inland waters into the terrestrial carbon budget Ecosystems 2007 10 172 185 10.1007/s10021-006-9013-8
Cole, J. J., Prairie, Y. T., Caraco, N. F., McDowell, W. H., Tranvik, L. J., Striegl, R. G., Duarte, C. M., Kortelainen, P., Downing, J. A., Middelburg, J. J., & Melack, J. (2007). Plumbing the global carbon cycle: Integrating inland waters into the terrestrial carbon budget. Ecosystems, 10, 172–185. 10.1007/s10021-006-9013-810.1007/s10021-006-9013-8
Cox R Lowe DR Cullers RL The influence of sediment recycling and basement composition on evolution of mudrock chemistry in the southwestern United States Geochimica Et Cosmochimica Acta 1995 59 14 2919 2940 10.1007/s10021-006-9013-8
Cox, R., Lowe, D. R., & Cullers, R. L. (1995). The influence of sediment recycling and basement composition on evolution of mudrock chemistry in the southwestern United States. Geochimica Et Cosmochimica Acta, 59(14), 2919–2940. 10.1007/s10021-006-9013-810.1007/s10021-006-9013-8
Dearing JA Jones RT Shen J Yang X Boyle JF Foster GC Crook DS Elvin MJD Using multiple archives to understand past and present climate–human–environment interactions: The lake Erhai catchment, Yunnan Province, China Journal of Paleolimnology 2008 40 3 31 10.1007/s10933-007-9182-2
Dearing, J. A., Jones, R. T., Shen, J., Yang, X., Boyle, J. F., Foster, G. C., Crook, D. S., & Elvin, M. J. D. (2008). Using multiple archives to understand past and present climate–human–environment interactions: The lake Erhai catchment, Yunnan Province, China. Journal of Paleolimnology, 40, 3–31. 10.1007/s10933-007-9182-210.1007/s10933-007-9182-2
Dominik J Leoni S Cassin D Guarneri I Bellucci LG Zonta R Eutrophication history and organic carbon burial rate recorded in sediment cores from the Mar Piccolo of Taranto (Italy) Environmental Science and Pollution Research 2023 30 19 56713 56730 10.1007/s11356-023-26317-6 36929258
Dominik, J., Leoni, S., Cassin, D., Guarneri, I., Bellucci, L. G., & Zonta, R. (2023). Eutrophication history and organic carbon burial rate recorded in sediment cores from the Mar Piccolo of Taranto (Italy). Environmental Science and Pollution Research, 30(19), 56713–56730. 10.1007/s11356-023-26317-636929258 10.1007/s11356-023-26317-6
Dong X Anderson NJ Yang X Chen X Shen JI Carbon burial by shallow lakes on the Yangtze floodplain and its relevance to regional carbon sequestration Global Change Biology 2012 18 7 2205 2217 10.1111/j.1365-2486.2012.02697.x
Dong, X., Anderson, N. J., Yang, X., Chen, X., & Shen, J. I. (2012). Carbon burial by shallow lakes on the Yangtze floodplain and its relevance to regional carbon sequestration. Global Change Biology, 18(7), 2205–2217. 10.1111/j.1365-2486.2012.02697.x10.1111/j.1365-2486.2012.02697.x
Downing, J. A., Cole, J. J., Middelburg, J. J., Striegl, R. G., Duarte, C. M., Kortelainen, P., Prairie, Y. T., & Laube, K. A. (2008). Sediment organic carbon burial in agriculturally eutrophic impoundments over the last century. Global Biogeochemical Cycles, 22(1). 10.1029/2006GB002854
EbrahimiMohammadi MS Sadeghi SH Chapi K Analysis of runoff, suspended sediment and nutrient yield from different tributaries to Zarivar lake in event and base flows Journal of Soil Water Resource Conservation 2012 2 1 61 76
EbrahimiMohammadi, M. S., Sadeghi, S. H., & Chapi, K. (2012). Analysis of runoff, suspended sediment and nutrient yield from different tributaries to Zarivar lake in event and base flows. Journal of Soil Water Resource Conservation, 2(1), 61–76. In Persian.
Einola, E., Rantakari, M., Kankaala, P., Kortelainen, P., Ojala, A., Pajunen, H., Mäkelä, S., & Arvola, L. (2011). Carbon pools and fluxes in a chain of five boreal lakes: A dry and wet year comparison. Journal of Geophysical Research: Biogeosciences, 116(G3). 10.1029/2010JG001636
El-Safa MMA Elsayed S Elsherbiny O Elmetwalli AH Gad M Moghanm FS Eid EM Taher MA El-Morsy MH Osman HE Saleh AH Environmental assessment of potentially toxic elements using pollution indices and data-driven modeling in surface sediment of the littoral shelf of the Mediterranean Sea coast and Gamasa estuary, Egypt Journal of Marine Science and Engineering 2022 10 6 816 10.3390/jmse10060816
El-Safa, M. M. A., Elsayed, S., Elsherbiny, O., Elmetwalli, A. H., Gad, M., Moghanm, F. S., Eid, E. M., Taher, M. A., El-Morsy, M. H., Osman, H. E., & Saleh, A. H. (2022). Environmental assessment of potentially toxic elements using pollution indices and data-driven modeling in surface sediment of the littoral shelf of the Mediterranean Sea coast and Gamasa estuary, Egypt. Journal of Marine Science and Engineering, 10(6), 816.10.3390/jmse10060816
Enderle DI Weih RC Jr Integrating supervised and unsupervised classification methods to develop a more accurate land cover classification Journal of the Arkansas Academy of Science 2005 59 1 65 73
Enderle, D. I., & Weih, R. C., Jr. (2005). Integrating supervised and unsupervised classification methods to develop a more accurate land cover classification. Journal of the Arkansas Academy of Science, 59(1), 65–73.
Engstrom DR Long-term changes in iron and phosphorus sedimentation in Vadnais Lake, Minnesota, resulting from ferric chloride addition and hypolimnetic aeration Lake and Reservoir Management 2005 21 1 95 105 10.1080/07438140509354417
Engstrom, D. R. (2005). Long-term changes in iron and phosphorus sedimentation in Vadnais Lake, Minnesota, resulting from ferric chloride addition and hypolimnetic aeration. Lake and Reservoir Management, 21(1), 95–105. 10.1080/0743814050935441710.1080/07438140509354417
Farquhar GD Ehleringer JR Hubick KT Carbon isotope discrimination and photosynthesis Annual Review of Plant Biology 1989 40 1 503 537 10.1146/annurev.pp.40.060189.002443
Farquhar, G. D., Ehleringer, J. R., & Hubick, K. T. (1989). Carbon isotope discrimination and photosynthesis. Annual Review of Plant Biology, 40(1), 503–537.10.1146/annurev.pp.40.060189.002443
Ferland ME Prairie YT Teodoru C del Giorgio PA Linking organic carbon sedimentation, burial efficiency, and long-term accumulation in boreal lakes Journal of Geophysical Research: Biogeosciences 2014 119 5 836 847 10.1002/2013JG002345
Ferland, M. E., Prairie, Y. T., Teodoru, C., & del Giorgio, P. A. (2014). Linking organic carbon sedimentation, burial efficiency, and long-term accumulation in boreal lakes. Journal of Geophysical Research: Biogeosciences, 119(5), 836–847. 10.1002/2013JG00234510.1002/2013JG002345
Foster GC Dearing JA Jones RT Crook DS Siddle DJ Harvey AM James PA Appleby PG Thompson R Nicholson J Loizeau JL Meteorological and land use controls on past and present hydro-geomorphic processes in the pre-alpine environment: An integrated lake–catchment study at the Petit Lac d’Annecy, France Hydrological Processes 2003 17 16 3287 3305 10.1002/hyp.1387
Foster, G. C., Dearing, J. A., Jones, R. T., Crook, D. S., Siddle, D. J., Harvey, A. M., James, P. A., Appleby, P. G., Thompson, R., Nicholson, J., & Loizeau, J. L. (2003). Meteorological and land use controls on past and present hydro-geomorphic processes in the pre-alpine environment: An integrated lake–catchment study at the Petit Lac d’Annecy, France. Hydrological Processes, 17(16), 3287–3305. 10.1002/hyp.138710.1002/hyp.1387
Foucher A Chaboche PA Sabatier P Evrard O A worldwide meta-analysis (1977–2020) of sediment core dating using fallout radionuclides including 137Cs and 210Pbxs Earth System Science Data Discussions 2021 2021 1 61 10.5194/essd-13-4951-2021
Foucher, A., Chaboche, P. A., Sabatier, P., & Evrard, O. (2021). A worldwide meta-analysis (1977–2020) of sediment core dating using fallout radionuclides including 137Cs and 210Pbxs. Earth System Science Data Discussions, 2021, 1–61. 10.5194/essd-13-4951-202110.5194/essd-13-4951-2021
France RL Carbon-13 enrichment in benthic compared to planktonic algae: Foodweb implications Oceanographic Literature Review 1996 5 43 490
France, R. L. (1996). Carbon-13 enrichment in benthic compared to planktonic algae: Foodweb implications. Oceanographic Literature Review, 5(43), 490.
Gibbs, M.M. (2014). Protocols on the use of the CSSI technique to identify and apportion soil sources from land use. NIWA report HAM2013–106 to Joint FAO/IAEA Division of Nuclear Techniques in Food and Agriculture (Contract Number 15491/R1) available at https://www-naweb.iaea.org/nafa/swmn/public/CSSI-technique-protocols-revised-2013.pdf
Håkanson, L. (1980). An ecological risk index for aquatic pollution control. A sedimentological approach. Water research, 14(8), 975–1001.
Hamidian AH Hassanzadeh M Can the PAMOLARE 1-layer model predict eutrophication in hypertrophic lakes? A case study: The Zaribar Lake, Iran International Journal of Aquatic Biology 2013 1 3 100 108 10.22034/ijab.v1i3.58
Hamidian, A. H., & Hassanzadeh, M. (2013). Can the PAMOLARE 1-layer model predict eutrophication in hypertrophic lakes? A case study: The Zaribar Lake, Iran. International Journal of Aquatic Biology, 1(3), 100–108. 10.22034/ijab.v1i3.5810.22034/ijab.v1i3.58
Harmon TS Smoak JM Waters MN Sanders CJ Hydrologic fragmentation-induced eutrophication in Dove Sound, Upper Florida Keys, USA Environmental Earth Sciences 2014 71 4387 4395 10.1007/s12665-013-2832-y
Harmon, T. S., Smoak, J. M., Waters, M. N., & Sanders, C. J. (2014). Hydrologic fragmentation-induced eutrophication in Dove Sound, Upper Florida Keys, USA. Environmental Earth Sciences, 71, 4387–4395. 10.1007/s12665-013-2832-y10.1007/s12665-013-2832-y
Heathcote AJ Ramstack Hobbs JM Anderson NJ Frings P Engstrom DR Downing JA Diatom floristic change and lake paleoproduction as evidence of recent eutrophication in shallow lakes of the midwestern USA Journal of Paleolimnology 2015 53 17 34 10.1007/s10933-014-9804-4
Heathcote, A. J., Ramstack Hobbs, J. M., Anderson, N. J., Frings, P., Engstrom, D. R., & Downing, J. A. (2015). Diatom floristic change and lake paleoproduction as evidence of recent eutrophication in shallow lakes of the midwestern USA. Journal of Paleolimnology, 53, 17–34. 10.1007/s10933-014-9804-410.1007/s10933-014-9804-4
Hobbs WO Engstrom DR Scottler SP Zimmer KD Cotner JB Estimating modern carbon burial rates in lakes using a single sediment sample Limnology and Oceanography: Methods 2013 11 6 316 326 10.4319/lom.2013.11.316
Hobbs, W. O., Engstrom, D. R., Scottler, S. P., Zimmer, K. D., & Cotner, J. B. (2013). Estimating modern carbon burial rates in lakes using a single sediment sample. Limnology and Oceanography: Methods, 11(6), 316–326. 10.4319/lom.2013.11.31610.4319/lom.2013.11.316
Hudak, A.T., & Wessman, C.A. (1997). Textural analysis of aerial photography to characterize large scale land cover change. In Proceedings of the ESRI Users Conference (Vol. 643), San Diego, CA. https://www.esri.com/library/userconf/proc97/PROC97/TO650/PAP643P. Accessed 15 Aug 2020.
IAEA (International Atomic Energy Agency) Collection and preparation of bottom sediment samples for analysis of radionuclides and trace elements IAEA TECDOC No. 2003 1360 136
IAEA (International Atomic Energy Agency). (2003). Collection and preparation of bottom sediment samples for analysis of radionuclides and trace elements. IAEA TECDOC No., 1360, 136.
IMA (Iranian Ministry of Agriculture). (2008). Detailed and executive watershed management studies of Zarivar Lake catchment, Marivan (1996–7), Bureau of statistics and information. (In Persian)
Jones B Manning DA Comparison of geochemical indices used for the interpretation of palaeoredox conditions in ancient mudstones Chemical Geology 1994 111 1–4 111 129 10.1016/0009-2541(94)90085-X
Jones, B., & Manning, D. A. (1994). Comparison of geochemical indices used for the interpretation of palaeoredox conditions in ancient mudstones. Chemical Geology, 111(1–4), 111–129. 10.1016/0009-2541(94)90085-X10.1016/0009-2541(94)90085-X
Kastowski, M., Hinderer, M., & Vecsei, A. (2011). Long‐term carbon burial in European lakes: Analysis and estimate. Global Biogeochemical Cycles, 25(3). 10.1029/2010GB003874.
Kemp P Sear D Collins A Naden P Jones I The impacts of fine sediment on riverine fish Hydrological Processes 2011 25 11 1800 1821 10.1002/hyp.7940
Kemp, P., Sear, D., Collins, A., Naden, P., & Jones, I. (2011). The impacts of fine sediment on riverine fish. Hydrological Processes, 25(11), 1800–1821. 10.1002/hyp.794010.1002/hyp.7940
Khodadadi M Alewell C Mirzaei M Ehssan-Malahat E Asadzadeh F Strauss P Meusburger K Understanding deforestation impacts on soil erosion rates using 137Cs, 239+240Pu, and 210Pbex and soil physicochemical properties in western Iran Journal of Environmental Radioactivity 2023 257 107078 10.1016/j.jenvrad.2022.107078 36423536
Khodadadi, M., Alewell, C., Mirzaei, M., Ehssan-Malahat, E., Asadzadeh, F., Strauss, P., & Meusburger, K. (2023). Understanding deforestation impacts on soil erosion rates using 137Cs, 239+240Pu, and 210Pbex and soil physicochemical properties in western Iran. Journal of Environmental Radioactivity, 257, 107078. 10.1016/j.jenvrad.2022.10707836423536 10.1016/j.jenvrad.2022.107078
Khodadadi M Meusburger K Mirzaei M Strauss P Blake WH Moghaseh E Alewell C Spatial cross-correlation of surface soil physicochemical properties with soil erosion estimated by fallout radionuclides in croplands in a semi-humid region of Iran CATENA 2024 237 107836 10.1016/j.catena.2024.107836
Khodadadi, M., Meusburger, K., Mirzaei, M., Strauss, P., Blake, W. H., Moghaseh, E., & Alewell, C. (2024). Spatial cross-correlation of surface soil physicochemical properties with soil erosion estimated by fallout radionuclides in croplands in a semi-humid region of Iran. CATENA, 237, 107836. 10.1016/j.catena.2024.10783610.1016/j.catena.2024.107836
Khodadadi, M., Blake, W.H., Ajili, F., Ahmadi, E., Toloza, A., & Fuentealba, L.O. (2022). Evaluating the impact of installing check-dams on the sediment sources using geochemical fingerprints in Zarivar Lake watershed, Iran. International Symposium on Managing Land and Water for Climate-Smart Agriculture, pp. 25–29 Jul 2022, ID #225.
Kondolf GM Gao Y Annandale GW Morris GL Jiang E Zhang J Cao Y Carling P Fu K Guo Q Hotchkiss R Sustainable sediment management in reservoirs and regulated rivers: Experiences from five continents Earth's Future 2014 2 5 256 280 10.1002/2013EF000184
Kondolf, G. M., Gao, Y., Annandale, G. W., Morris, G. L., Jiang, E., Zhang, J., Cao, Y., Carling, P., Fu, K., Guo, Q., & Hotchkiss, R. (2014). Sustainable sediment management in reservoirs and regulated rivers: Experiences from five continents. Earth’s Future, 2(5), 256–280. 10.1002/2013EF00018410.1002/2013EF000184
Kortelainen P Pajunen H Rantakari M Saarnisto M A large carbon pool and small sink in boreal Holocene lake sediments Global Change Biology 2004 10 10 1648 1653 10.1111/j.1365-2486.2004.00848.x
Kortelainen, P., Pajunen, H., Rantakari, M., & Saarnisto, M. (2004). A large carbon pool and small sink in boreal Holocene lake sediments. Global Change Biology, 10(10), 1648–1653. 10.1111/j.1365-2486.2004.00848.x10.1111/j.1365-2486.2004.00848.x
Kortelainen P Rantakari M Pajunen H Huttunen JT Mattsson T Juutinen S Larmola T Alm J Silvola J Martikainen PJ Carbon evasion/accumulation ratio in boreal lakes is linked to nitrogen Global Biogeochemical Cycles 2013 27 2 363 374 10.1002/gbc.20036
Kortelainen, P., Rantakari, M., Pajunen, H., Huttunen, J. T., Mattsson, T., Juutinen, S., Larmola, T., Alm, J., Silvola, J., & Martikainen, P. J. (2013). Carbon evasion/accumulation ratio in boreal lakes is linked to nitrogen. Global Biogeochemical Cycles, 27(2), 363–374. 10.1002/gbc.2003610.1002/gbc.20036
Krishnaswamy S Lal D Martin JM Meybeck M Geochronology of lake sediments Earth and Planetary Science Letters 1971 11 1–5 407 414 10.1016/0012-821X(71)90202-0
Krishnaswamy, S., Lal, D., Martin, J. M., & Meybeck, M. (1971). Geochronology of lake sediments. Earth and Planetary Science Letters, 11(1–5), 407–414.10.1016/0012-821X(71)90202-0
Lin Q Liu E Zhang E Nath B Bindler R Liu J Shen J Organic carbon burial in a large, deep alpine lake (southwest China) in response to changes in climate, land use and nutrient supply over the past~ 100 years CATENA 2021 202 105240 10.1016/j.catena.2021.105240
Lin, Q., Liu, E., Zhang, E., Nath, B., Bindler, R., Liu, J., & Shen, J. (2021). Organic carbon burial in a large, deep alpine lake (southwest China) in response to changes in climate, land use and nutrient supply over the past~ 100 years. CATENA, 202, 105240. 10.1016/j.catena.2021.10524010.1016/j.catena.2021.105240
Liu Y Ren Z Qu X Zhang M Yu Y Zhang Y Peng W Microbial community structure and functional properties in permanently and seasonally flooded areas in Poyang Lake Scientific Reports 2020 10 1 4819 10.1038/s41598-020-61569-z 32179796
Liu, Y., Ren, Z., Qu, X., Zhang, M., Yu, Y., Zhang, Y., & Peng, W. (2020). Microbial community structure and functional properties in permanently and seasonally flooded areas in Poyang Lake. Scientific Reports, 10(1), 4819. 10.1038/s41598-020-61569-z32179796 10.1038/s41598-020-61569-z
Mahamuda, A. (2015). The sedimentology of Gambaga massifs of the northeastern Voltaian, Ghana (Doctoral dissertation, University of Ghana).
Maroufi H Ehyaei A Zarivar Lake: The largest freshwater spring in Iran Iran Nature 2016 1 1 65 60
Maroufi, H., & Ehyaei, A. (2016). Zarivar Lake: The largest freshwater spring in Iran. Iran Nature, 1(1), 65–60. In Persian.
Mermut, A. R., & Acton, D. F. (1984). The age of some Holocene soils on the Ear Lake terraces in Saskatchewan. Canadian journal of soil science, 64(2), 163–172.
Morellón M Valero-Garcés B González-Sampériz P Vegas-Vilarrúbia T Rubio E Rieradevall M Delgado-Huertas A Mata P Romero O Engstrom DR López-Vicente M Climate changes and human activities recorded in the sediments of Lake Estanya (NE Spain) during the Medieval Warm Period and Little Ice Age Journal of Paleolimnology 2011 46 423 452 10.1007/s10933-009-9346-3
Morellón, M., Valero-Garcés, B., González-Sampériz, P., Vegas-Vilarrúbia, T., Rubio, E., Rieradevall, M., Delgado-Huertas, A., Mata, P., Romero, O., Engstrom, D. R., & López-Vicente, M. (2011). Climate changes and human activities recorded in the sediments of Lake Estanya (NE Spain) during the Medieval Warm Period and Little Ice Age. Journal of Paleolimnology, 46, 423–452. 10.1007/s10933-009-9346-310.1007/s10933-009-9346-3
Muhati GL Olago D Olaka L Land use and land cover changes in a sub-humid Montane forest in an arid setting: A case study of the Marsabit forest reserve in northern Kenya Global Ecology and Conservation 2018 16 e00512 10.1016/j.gecco.2018.e00512
Muhati, G. L., Olago, D., & Olaka, L. (2018). Land use and land cover changes in a sub-humid Montane forest in an arid setting: A case study of the Marsabit forest reserve in northern Kenya. Global Ecology and Conservation, 16, e00512. 10.1016/j.gecco.2018.e0051210.1016/j.gecco.2018.e00512
Murphy J Riley JP A modified single solution method for the determination of phosphate in natural waters Analytica Chimica Acta 1962 27 31 36 10.1016/S0003-2670(00)88444-5
Murphy, J., & Riley, J. P. (1962). A modified single solution method for the determination of phosphate in natural waters. Analytica Chimica Acta, 27, 31–36.10.1016/S0003-2670(00)88444-5
Nunes, J. P., & Nearing, M. A. (2011). Modelling impacts of climatic change: case studies using the new generation of erosion models. Handbook of erosion modelling, John Wiley & Sons, Ltd, 289–312.
Olley JM Wasson RJ Changes in the flux of sediment in the Upper Murrumbidgee catchment, Southeastern Australia, since European settlement Hydrological Processes 2003 17 16 3307 3320 10.1002/hyp.1388
Olley, J. M., & Wasson, R. J. (2003). Changes in the flux of sediment in the Upper Murrumbidgee catchment, Southeastern Australia, since European settlement. Hydrological Processes, 17(16), 3307–3320. 10.1002/hyp.138810.1002/hyp.1388
Patwary MAA Parvaiz I IRS-LISS-III and PAN data analysis for landslide susceptibility mapping using heuristic approach in active tectonic region of Himalaya Journal of Indian Society Remote Sensing 2009 37 3 493 509 10.1007/s12524-009-0036-4
Patwary, M. A. A., & Parvaiz, I. (2009). IRS-LISS-III and PAN data analysis for landslide susceptibility mapping using heuristic approach in active tectonic region of Himalaya. Journal of Indian Society Remote Sensing, 37(3), 493–509. 10.1007/s12524-009-0036-410.1007/s12524-009-0036-4
Perie, C., & Ouimet, R. (2008). Organic carbon, organic matter and bulk density relationships in boreal forest soils. Canadian journal of soil science, 88(3), 315–325.
Pettijohn, F.J. (1975). Sedimentary rocks (Vol. 3: 628). New York: Harper & Row.
Radbourne AD Elliott JA Maberly SC Ryves DB Anderson NJ The impacts of changing nutrient load and climate on a deep, eutrophic, monomictic lake Freshwater Biology 2019 64 6 1169 1182 10.1111/fwb.13293
Radbourne, A. D., Elliott, J. A., Maberly, S. C., Ryves, D. B., & Anderson, N. J. (2019). The impacts of changing nutrient load and climate on a deep, eutrophic, monomictic lake. Freshwater Biology, 64(6), 1169–1182. 10.1111/fwb.1329310.1111/fwb.13293
Rippey, B., & Douglas, R. W. (2004). Reconstructing regional‐scale lead contamination of the atmosphere (1850–1980) in the United Kingdom and Ireland using lake sediments. Global Biogeochemical Cycles, 18(4). 10.1029/2004GB002305
Salmi, T. (2002). Detecting trends of annual values of atmospheric pollutants by the Mann-Kendall test and Sen’s slope estimates-the Excel template application MAKESENS. Ilmatieteen laitos.
Sanders CJ Santos IR Silva-Filho EV Patchineelam SR Mercury flux to estuarine sediments, derived from Pb-210 and Cs-137 geochronologies (Guaratuba Bay, Brazil) Marine Pollution Bulletin 2006 52 9 1085 1089 10.1016/j.marpolbul.2006.06.004 16901512
Sanders, C. J., Santos, I. R., Silva-Filho, E. V., & Patchineelam, S. R. (2006). Mercury flux to estuarine sediments, derived from Pb-210 and Cs-137 geochronologies (Guaratuba Bay, Brazil). Marine Pollution Bulletin, 52(9), 1085–1089. 10.1016/j.marpolbul.2006.06.00416901512 10.1016/j.marpolbul.2006.06.004
Sanders CJ Eyre BD Santos IR Machado W Luiz-Silva W Smoak JM Breithaupt JL Ketterer ME Sanders L Marotta H Silva-Filho E Elevated rates of organic carbon, nitrogen, and phosphorus accumulation in a highly impacted mangrove wetland Geophysical Research Letters 2014 41 7 2475 2480 10.1002/2014GL059789
Sanders, C. J., Eyre, B. D., Santos, I. R., Machado, W., Luiz-Silva, W., Smoak, J. M., Breithaupt, J. L., Ketterer, M. E., Sanders, L., Marotta, H., & Silva-Filho, E. (2014). Elevated rates of organic carbon, nitrogen, and phosphorus accumulation in a highly impacted mangrove wetland. Geophysical Research Letters, 41(7), 2475–2480. 10.1002/2014GL05978910.1002/2014GL059789
Sanders LM Taffs KH Stokes DJ Sanders CJ Smoak JM Enrich-Prast A Macklin PA Santos IR Marotta H Carbon accumulation in Amazonian floodplain lakes: A significant component of Amazon budgets? Limnology and Oceanography Letters 2017 2 1 29 35 10.1002/lol2.10034
Sanders, L. M., Taffs, K. H., Stokes, D. J., Sanders, C. J., Smoak, J. M., Enrich-Prast, A., Macklin, P. A., Santos, I. R., & Marotta, H. (2017). Carbon accumulation in Amazonian floodplain lakes: A significant component of Amazon budgets? Limnology and Oceanography Letters, 2(1), 29–35. 10.1002/lol2.1003410.1002/lol2.10034
Sharifinia M Ramezanpour Z Imanpour J Mahmoudifard A Rahmani T Water quality assessment of the Zarivar Lake using physico-chemical parameters and NSF-WQI indicator, Kurdistan Province-Iran International Journal of Advanced Biological and Biomedical Research 2013 1 3 302 312
Sharifinia, M., Ramezanpour, Z., Imanpour, J., Mahmoudifard, A., & Rahmani, T. (2013). Water quality assessment of the Zarivar Lake using physico-chemical parameters and NSF-WQI indicator, Kurdistan Province-Iran. International Journal of Advanced Biological and Biomedical Research, 1(3), 302–312.
Smith JN Why should we believe 210Pb sediment geochronologies? Journal of Environmental Radioactivity 2001 55 121 123 10.1016/s0265-931x(00)00152-1 11398372
Smith, J. N. (2001). Why should we believe 210Pb sediment geochronologies? Journal of Environmental Radioactivity, 55, 121–123. 10.1016/s0265-931x(00)00152-111398372 10.1016/s0265-931x(00)00152-1
Smith VH Tilman GD Nekola JC Eutrophication: Impacts of excess nutrient inputs on freshwater, marine, and terrestrial ecosystems Environmental Pollution 1999 100 1–3 179 196 10.1016/S0269-7491(99)00091-3 15093117
Smith, V. H., Tilman, G. D., & Nekola, J. C. (1999). Eutrophication: Impacts of excess nutrient inputs on freshwater, marine, and terrestrial ecosystems. Environmental Pollution, 100(1–3), 179–196. 10.1016/S0269-7491(99)00091-315093117 10.1016/S0269-7491(99)00091-3
Soil Taxonomy (2014). Soil survey staff. Keys to soil taxonomy, 10th ed. United States Department of Agriculture, Natural Resources Conservation Services, 331.
Sterner RW Review paper on the phosphorus limitation paradigm for lakes International Review Hydrobiology 2008 93 4–5 433 445 10.1002/iroh.200811068
Sterner, R. W. (2008). Review paper on the phosphorus limitation paradigm for lakes. International Review Hydrobiology, 93(4–5), 433–445. 10.1002/iroh.20081106810.1002/iroh.200811068
Tiessen KHD Elliott JA Stainton M Yarotski J Flaten DN Lobb DA The effectiveness of small-scale headwater storage dams and reservoirs on stream water quality and quantity in the Canadian Prairies Journal of Soil and Water Conservation 2011 66 3 158 171 10.2489/jswc.66.3.158
Tiessen, K. H. D., Elliott, J. A., Stainton, M., Yarotski, J., Flaten, D. N., & Lobb, D. A. (2011). The effectiveness of small-scale headwater storage dams and reservoirs on stream water quality and quantity in the Canadian Prairies. Journal of Soil and Water Conservation, 66(3), 158–171.10.2489/jswc.66.3.158
TorabiKachoosangi F Karbassi A Sarang A Noori R Sedimentation rate determination and heavy metal pollution assessment in Zariwar Lake, Iran SN Applied Sciences 2020 2 1 10 10.1007/s42452-020-03279-9
TorabiKachoosangi, F., Karbassi, A., Sarang, A., & Noori, R. (2020). Sedimentation rate determination and heavy metal pollution assessment in Zariwar Lake, Iran. SN Applied Sciences, 2, 1–10. 10.1007/s42452-020-03279-910.1007/s42452-020-03279-9
Tranvik LJ Downing JA Cotner JB Loiselle SA Striegl RG Ballatore TJ Dillon P Finlay K Fortino K Knoll LB Kortelainen PL Lakes and reservoirs as regulators of carbon cycling and climate Limnology and Oceanography 2009 54 6part2 2298 2314 10.4319/lo.2009.54.6_part_2.2298
Tranvik, L. J., Downing, J. A., Cotner, J. B., Loiselle, S. A., Striegl, R. G., Ballatore, T. J., Dillon, P., Finlay, K., Fortino, K., Knoll, L. B., & Kortelainen, P. L. (2009). Lakes and reservoirs as regulators of carbon cycling and climate. Limnology and Oceanography, 54(6part2), 2298–2314. 10.4319/lo.2009.54.6_part_2.229810.4319/lo.2009.54.6_part_2.2298
Trumbore S Radiocarbon and soil carbon dynamics Annual Review Earth Planet Science 2009 37 1 47 66 10.1146/annurev.earth.36.031207.124300
Trumbore, S. (2009). Radiocarbon and soil carbon dynamics. Annual Review Earth Planet Science, 37(1), 47–66. 10.1146/annurev.earth.36.031207.12430010.1146/annurev.earth.36.031207.124300
Tylmann W Bonk A Goslar T Wulf S Grosjean M Calibrating 210Pb dating results with varve chronology and independent chronostratigraphic markers: Problems and implications Quaternary Geochronology 2016 32 1 10 10.1016/j.quageo.2015.11.004
Tylmann, W., Bonk, A., Goslar, T., Wulf, S., & Grosjean, M. (2016). Calibrating 210Pb dating results with varve chronology and independent chronostratigraphic markers: Problems and implications. Quaternary Geochronology, 32, 1–10. 10.1016/j.quageo.2015.11.00410.1016/j.quageo.2015.11.004
Valero-Garcés BL González-Sampériz P Navas A Machín J Mata P Delgado-Huertas A Bao R Moreno A Carrión JS Schwalb A González-Barrios A Human impact since medieval times and recent ecological restorationin a Mediterranean Lake: The Laguna Zoñar, Southern Spain Journal of Paleolimnology 2006 35 441 465 10.1007/s10933-005-1995-2
Valero-Garcés, B. L., González-Sampériz, P., Navas, A., Machín, J., Mata, P., Delgado-Huertas, A., Bao, R., Moreno, A., Carrión, J. S., Schwalb, A., & González-Barrios, A. (2006). Human impact since medieval times and recent ecological restorationin a Mediterranean Lake: The Laguna Zoñar, Southern Spain. Journal of Paleolimnology, 35, 441–465. 10.1007/s10933-005-1995-210.1007/s10933-005-1995-2
Verburg P The need to correct for the Suess effect in the application of δ13C in sediment of autotrophic Lake Tanganyika, as a productivity proxy in the Anthropocene Journal of Paleolimnology 2007 37 4 591 602 10.1007/s10933-006-9056-z
Verburg, P. (2007). The need to correct for the Suess effect in the application of δ13C in sediment of autotrophic Lake Tanganyika, as a productivity proxy in the Anthropocene. Journal of Paleolimnology, 37(4), 591–602. 10.1007/s10933-006-9056-z10.1007/s10933-006-9056-z
Waters, M.N., Kenney, W.F., Brenner, M. and Webster, B.C., (2019). Organic carbon sequestration in sediments of subtropical Florida lakes. PloS One, 14(12), e0226273. 10.1371/journal.pone.0226273
Wedepohl KH The composition of the continental crust Geochimica Et Cosmochimica Acta 1995 59 1217 1239 10.1016/0016-7037(95)00038-2
Wedepohl, K. H. (1995). The composition of the continental crust. Geochimica Et Cosmochimica Acta, 59, 1217–1239.10.1016/0016-7037(95)00038-2
Willis JP Turner K Pritchard G XRF in the workplace: A guide to practical XRF spectrometry 2011 PANalytical
Willis, J. P., Turner, K., & Pritchard, G. (2011). XRF in the workplace: A guide to practical XRF spectrometry. PANalytical.
WRB (World Reference Base for soil resources) (2014). International soil classification system for naming soils and creating legends for soil maps. World Soil Resources Report, 106.
Wynants M Munishi L Mtei K Bodé S Patrick A Taylor A Gilvear D Ndakidemi P Blake WH Boeckx P Soil erosion and sediment transport in Tanzania: Part I–sediment source tracing in three neighbouring river catchments Earth Surface Processes and Landforms 2021 46 15 3096 3111 10.1002/esp.5217
Wynants, M., Munishi, L., Mtei, K., Bodé, S., Patrick, A., Taylor, A., Gilvear, D., Ndakidemi, P., Blake, W. H., & Boeckx, P. (2021). Soil erosion and sediment transport in Tanzania: Part I–sediment source tracing in three neighbouring river catchments. Earth Surface Processes and Landforms, 46(15), 3096–3111. 10.1002/esp.521710.1002/esp.5217
Yang P Yang C Yin H Dynamics of phosphorus composition in suspended particulate matter from a turbid eutrophic shallow lake (Lake Chaohu, China): Implications for phosphorus cycling and management Science of the Total Environment 2020 741 140203 10.1016/j.scitotenv.2020.140203 32570068
Yang, P., Yang, C., & Yin, H. (2020). Dynamics of phosphorus composition in suspended particulate matter from a turbid eutrophic shallow lake (Lake Chaohu, China): Implications for phosphorus cycling and management. Science of the Total Environment, 741, 140203. 10.1016/j.scitotenv.2020.14020332570068 10.1016/j.scitotenv.2020.140203
Zaid SM Geochemistry of sandstones from the Pliocene Gabir Formation, north Marsa Alam, Red Sea, Egypt: Implication for provenance, weathering and tectonic setting Journal of African Earth Science 2015 102 1 17 10.1016/j.jafrearsci.2014.10.016
Zaid, S. M. (2015). Geochemistry of sandstones from the Pliocene Gabir Formation, north Marsa Alam, Red Sea, Egypt: Implication for provenance, weathering and tectonic setting. Journal of African Earth Science, 102, 1–17. 10.1016/j.jafrearsci.2014.10.01610.1016/j.jafrearsci.2014.10.016
Zan F Huo S Xi B Zhu C Liao H Zhang J Yeager KM A 100-year sedimentary record of natural and anthropogenic impacts on a shallow eutrophic lake, Lake Chaohu, China Journal of Environmental Monitoring 2012 14 3 804 816 10.1039/c1em10760g 22237732
Zan, F., Huo, S., Xi, B., Zhu, C., Liao, H., Zhang, J., & Yeager, K. M. (2012). A 100-year sedimentary record of natural and anthropogenic impacts on a shallow eutrophic lake, Lake Chaohu, China. Journal of Environmental Monitoring, 14(3), 804–816. 10.1039/c1em10760g22237732 10.1039/c1em10760g
