Water quality changes of surface waters can be used to assess human impact intensity, but of importance is to consider also impacts of climate change/variability and naturally occurring changes of environmental quality. In Latvia, during the recent decades a major reduction of anthropogenic pressure has happened due to restructuring of economy and industrial production, resulting in major decrease of loading of many groups of pollutants. However, trends and driving factors for other groups of substances have not been much studied. Long term (1980–2012) results of hydrochemical monitoring, performed in rivers of Latvia, are analysed during this study in connection with long-term sets of hydrological and heliophysical data, using standard statistical approaches. Our results indicate that variation of some hydrochemical values, for example COD and total Fe, show clearly visible decadal oscillated character, while variation of some other values, for example phosphate P and total P, show some individual characteristics of decadal oscillations. These results indicate the presence of a large scale, geochemical and geophysical significant process: multiannual pulse of catchment, driven by variation of solar irradiance through complex interactions between global atmospheric circulation, groundwater and surface waterbodies. The process described in our study is significant from a geochemical point of view and must be taken into account in prediction of water quality and quantity. Impacts of natural processes should be considered in the planning of environmental policy.
The runoff of dissolved substances by riverine waters is an important factor in evaluation of not only water quality, but more important, of the character of the processes in waterbody basins (Henriksen et al., 1998). Rivers act as conveyers transporting dissolved substances from continents to seas. Thus, the integral dimensions of adverse processes in water body basins and the increased geochemical mobility of water ingredients can be identified by studying surface-water-driven geochemical flows (Stumm and Morgan, 1996). Increased flows of dissolved substances can affect the living environment of aquatic biota, especially in the receiving water bodies, such as the Gulf of Rīga and the Baltic Sea, in the case of Latvia.
Chemical composition of surface waters is determined by numerous factors and it depends, at first, on natural conditions (physico-geographical character of the region, climate, hydrology, hydrogeology etc.), but also on the level and extent of human impact. Considering many different sources of water ingredients, the study of their contribution to the aquatic chemistry is at utmost importance in evaluation of water quality (Hounslow, 1995). Studies of aquatic chemistry have also been used to reveal the dimension of the human activity-related processes in water body basins which have changed water composition (both at global, regional and local level), to mention industrial production (wastewater and solid waste dumping) and agricultural activities (non-point pollution sources), as well as impacts of hydro-technical constructions and changes of land-use patterns. Also, impacts of climate change of water composition and flows of substances might be considered as major actuality (Whitehead et al., 2009).
Commonly, temporal changes of geochemical flows of dissolved substances, have been done in terms of linear trend analysis, even though the character of variance of these flows can be much more complex (Libiseller and Grimvall, 2002). At the same time river discharge regimes, surface water composition and flows of dissolved substances have not only linear trend of changes, but well expressed periodicity (Pekarova et al., 2003). Oscillations of the discharge intensity, with main frequency of about 20 and 13 years in rivers of Latvia, were reported by Klavins et al. (2002a). Periodic oscillations of groundwater levels and levels of deeper underground aquifers in Latvia, with a main frequency of about 11 years, were reported by Levina et al. (1995). So, it is very credible, that similar oscillational patterns are characteristic also for fluxes of dissolved substances with riverine waters.
The results for the 35-year-period of systematic observations of the concentrations of analysed substances allow estimating more or less credible values of flows of these substances in river waters. The aim of this study is to investigate the oscillation patterns in the temporal variation of flows of dissolved substances in surface waters in Latvia.
Latvia is a country situated on the eastern coast of the Baltic Sea in North
Eastern Europe. Chemical composition of surface waters of Latvia is
influenced by geographical, geological, and geochemical processes. At first,
possible impacts on geographical factors should be mentioned: the surface is
characterised by flat surface topography (57 % of Latvia's territory is
located below 100 m above the sea level, and only 2.5 % of the area
reaches 200 m). The surficial geology has been formed by Quaternary glacial
and ancient sea sediments and by recent exogenous processes. The most widely
distributed parent soil materials are moraine loam and sands which are
resistant to weathering processes. More than one half of Latvian soils are
humic podsols with comparatively high content of organic matter (Klavins et
al., 2002b). Latvia has a comparatively dense network of rivers. The climatic
conditions of Latvia are dominated by transport of cyclonic air masses from
the Atlantic Ocean, leading to comparatively high humidity, uneven
distribution of atmospheric precipitation through the year, mild winters and
moist summers (Klavins et al., 2003). The study area covered the major river
basins of Latvia (Lielupe and Daugava) which supply the dominant part of the
water runoff (24 km
Daugava, also known as Western Dvina is the largest river of Latvia and also
one of the biggest rivers in the Baltic region. It is 1005 km long and its
basin size is 87 900 km
Lielupe is the second biggest river of Latvia, with strongly different
catchment characteristics. Length of Lielupe is 119 km, basin area 17 600 km
Monitoring points under analysis.
Hydrochemical and hydrological data used in this study were collected by the
Centre for Geology, Hydrometeorology and Environment within the National
monitoring program. Monthly records of river discharge from gauging stations
on Daugava (Jēkabpils) and Lielupe (Kalnciems) were used. Sampling and
analytical methods followed the standard methods during all the study period
(Tsirkunov et al., 1992). From 1980 till 1993, water ingredients were
analysed according to standard methods used in USSR, as summarised in
Standard Methods (1973). After 1993 standard methods described in APHA (1992) were used. During 1992–1994, inter-calibration programmes, both at
national and international levels, have demonstrated general acceptability
of the used methods (Stålnacke, 1996). Average concentrations of
chemical substances (Fetot, TOC, COD, P-PO
Flux of total iron and nutrients were estimated, with the help of available
raw monitoring data, while flux of organic matter, taking into account
limits of data availability, was evaluated using COD values for the period
1980–2003 and TOC values for the period 2003–2012. TOC values for the first
23 years of study were calculated from COD concentrations using formula
proposed by Klavins et al. (2012).
As it can be seen in Fig. 2, the character of the variation of discharge values in the studied rivers obviously can be described as oscillated. Signs of oscillations with 20, 11, 8–9, 3–5 and 2 years long periods can be observed, that generally corresponds to the results of earlier studies (Glazacheva, 1988; Klavins et al., 2002a). Increase of water color, known also as brownification, was frequently reported in many lakes and rivers of the Northern hemisphere (Driscoll et al., 2003; Frey and Smith, 2005; Monteith et al., 2007; Erlandsson et al., 2008; Haaland et al., 2010), including Latvian rivers as well (Klavins et al., 2012). Our data show that despite low data quality and short period of observations (1980–2003), statistically significant upward trends of color values can be found in both rivers (Table 3). However, a closer look at obtained curves indicates also an oscillated character of this process in both cases (Fig. 3). The drivers behind brownification are heavily debated and have been ascribed to hydrological factors (Hongve et al., 2004; Erlandsson et al., 2008), as well as an increase in temperature, changes in land-use and reduced acid deposition (Correll et al., 2001; Freeman et al., 2001; Monteith et al., 2007).
Descriptive characteristics of data sets under analysis.
Although it has been difficult to find out a single factor causing brownification, there is an implicit consensus that the direct factor causing this process is the increased concentrations of terrestrially derived organic matter (OM) in the water (Kritzberg and Ekström, 2012).
Variation of discharge in studied rivers (24 monthly moving average).
Variation of water color in studied rivers.
As we see in Fig. 4, temporal variations of TOC flow in both studied
rivers are very similar. Mann-Kendall test results do not show any positive
trend in fluxes of organic carbon in the rivers under analysis. While the
oscillated character of data variance is evident, decadal oscillations, as
well as oscillations with shorter period (3–6 years) can be observed in both
rivers. The character of variation is similar to the character of variation
of discharge data, while statistically significant correlation (
However, it is necessary to admit that the increase in water color cannot be explained only by changes in the quantity and quality of OM, iron also must be mentioned as a potential source of color (Canfield et al., 1984; Kritzberg and Ekström, 2012; Ekström, 2013). Considering low solubility of iron compounds it becomes evident that iron compounds in surface waters exist in form of complexes with OM. Dissolved OM is known to bind to the surface of ferrihydrite and stabilize against aggregation by giving the particle surface a net negative charge (Tipping, 1981), so the flux of organic matter and the flux of iron are strongly connected amongst themselves.
Variation of TOC export in studied rivers.
Variation of total iron export in studied rivers and total solar irradiation (TSI).
Variation of total iron flux in Lielupe indicates clearly visible decadal
oscilations, while period of oscillations in Daugava is shorter (3–6, mainly
4 years). Observed oscillations are highly similar to the cyclic signal of
changes in total solar irradiance (TSI) with a periodicity of
Linear correlations between hydro-chemical values and river discharge.
Differences in both curves can be explained by the difference of iron sources
in the studied areas. In the catchment of Daugava as the most significant
source of iron, the leaching from soils must be mentioned, as indicated also
by a very strong correlation with features of covariation between fluxes of
TOC and Total Iron (
Variation of ammonia nitrogen export in studied rivers.
Variation of phosphate phosphorous export in studied rivers.
Results of Mann Kendall test.
The results of Mann Kendall test (Table 3) show that the decrease of nutrient export amounts can be observed in both rivers. These tendencies can be explained by significant reduction of nutrient loads that took place in last decade of XX century, after the collapse of the Soviet economics (Stalnacke, 1996). However, taking into account the oscillational character of nutrient flows, downward trends in Daugava also can be explained with the minimum in cycle with approximately 20 years long period. In this case, upward trends in nutrient export values can be predicted for Daugava catchment, irrespectively to the reduction of loads, in the next few years. The reduction of phosphorous and ammonia nitrogen exports in Lielupe appears to be more reliable, representing self-purification process of over fertilised agricultural area after sufficient cuts in fertiliser loads.
Oscillational patterns that can be defined as multiannual pulse of catchment are obviously characteristic of export values variation of the dissolved substances from the studied catchments. The oscillational character of geochemical flows can be explained by impacts of multiple ranges of cyclic natural events driven by changes in solar irradiation, e.g. oscillation of precipitation and river discharge, cyclic oscillations of groundwater aquafers, changes in evaporation intensity. Cyclic character of geochemical fluxes must be taken into account researching geochemical processes, e.g. brownification of waters and self-purification of soils. 30 years long sets of geochemical data are too short for a fully credible linear trend analysis.
We have failed to explain upward trends of color in Daugava and Lielupe rivers by increase of export of OM. Our results show that iron must be mentioned as the very significant factor affecting water color in the studied rivers. The character of brownification in both rivers is oscillational and shows significant connection with solar drivers.
Successful recovery of over-fertilized agricultural area is evident in the case of the Lielupe River. Longer observation period is needed to evaluate self-purification of the Daugava River catchment area. Upward trends of nutrient export values from Daugava catchment are possible in next few years, irrespectively to reduction of loads.
Datasets can be found in the Supplement.