Carbon dioxide dynamics in an agricultural headwater stream driven by hydrology and primary production

Abstract. Headwater streams are known to be hotspots for carbon dioxide
(CO2) emissions to the atmosphere and are hence important components in
landscape carbon balances. However, surprisingly little is known about
stream CO2 dynamics and emissions in agricultural settings, a land use
type that globally covers ca. 40 % of the continental area. Here we present
hourly measured in situ stream CO2 concentration data from a 11.3 km2 temperate agricultural headwater catchment covering more than 1
year (in total 339 d excluding periods of ice and snow cover). The stream
CO2 concentrations during the entire study period were generally high
(median 3.44 mg C L−1, corresponding to partial pressures (pCO2) of
4778 µatm) but were also highly variable (IQR = 3.26 mg C L−1). The CO2 concentration dynamics covered a variety of
different timescales from seasonal to hourly, with an interplay of
hydrological and biological controls. The hydrological control was strong
(although with both positive and negative influences dependent on
season), and CO2 concentrations changed rapidly in response to rainfall
and snowmelt events. However, during growing-season base flow and receding
flow conditions, aquatic primary production seemed to control the stream
CO2 dynamics, resulting in elevated diel patterns. During the dry summer
period, rapid rewetting following precipitation events generated high
CO2 pulses exceeding the overall median level of stream CO2 (up
to 3 times higher) observed during the whole study period. This finding
highlights the importance of stream intermittency and its effect on stream
CO2 dynamics. Given the observed high levels of CO2 and its
temporally variable nature, agricultural streams clearly need more attention
in order to understand and incorporate these considerable dynamics in large-scale extrapolations.



Introduction
Fluvial systems (streams and rivers) are estimated to dominate the inland water CO 2 source globally, surpassing CO 2 emissions from lakes and reservoirs by a factor of 6 ( Raymond et al., 2013). However, this estimate relies on a number of assumptions, and the scarcity of empirical data makes it uncertain. One of the critical gaps in the global upscaling is the lack of direct measurements from agriculture-dominated areas (Osborne et al., 2010). Globally, agricultural land covers about 40 % of the total continental area (Ramankutty et al., 2008), but there are few studies specifically focusing on the magnitude and dynamics of CO 2 emissions from agricultural streams. The few studies that do exist have shown that agricultural stream CO 2 concentrations are generally high and up to 5 times greater than those in streams draining forested areas which are more extensively studied (Borges et al., 2018;Bodmer et al., 2016;Wallin et al., 2018). For example, Bodmer et al. (2016) measured partial pressure of CO 2 (pCO 2 ) in German and Polish streams and examined differences between forested and agricultural catchments. They found that pCO 2 was generally 2-3 times higher in agricultural streams compared to streams draining forested areas. Similarly, Borges et al. (2018) found high CO 2 concentra-tions in streams and rivers dominated by agriculture in the river system Meuse, Belgium. They linked the higher pCO 2 in agricultural streams to elevated levels of dissolved organic carbon (DOC), particulate organic carbon (POC) and inorganic nitrogen. On the other hand, Deirmendjian et al. (2019) showed that there was no difference in pCO 2 between forest and cropland streams in south-west France despite higher pCO 2 in forest groundwater compared to cropland groundwater. They explained the similar stream pCO 2 by more efficient gas exchange in the forest streams compared to the low-gradient cropland streams.
There are numerous factors influencing CO 2 patterns in stream systems, and site-specific controls often dominate. Hence, large-scale generalizations are difficult to make (Crawford et al., 2017). Based on high-frequency data, CO 2 concentrations in streams draining nutrient-poor forest and peatlands, as well as tropical forests, are often found related to variations in stream discharge but with site-specific response patterns, with CO 2 found to be either positively or negatively related to stream discharge (Crawford et al., 2017;Dinsmore et al., 2013;Johnson et al., 2007). These response patterns have often been connected to the catchment characteristics and changes in hydrological pathways, which in turn control the dominant source areas (both from a vertical and lateral point of view) of CO 2 in the catchment soils Leith et al., 2015;Dinsmore and Billett, 2008). In contrast, other catchments lack a strong hydrological control and instead display clear diel cycles in stream CO 2 concentration, indicating a metabolic control (Crawford et al., 2017). Here the interplay of photosynthesis and respiration (in stream or terrestrial) could result in large day-to night-time differences in stream CO 2 .
These recent findings concerning dynamics and controls on stream CO 2 concentrations have been possible due to the development of cost-effective CO 2 sensors (e.g. Johnson et al., 2010;Bastviken et al., 2015) which have enabled continuous data collection covering relevant timescales (< hourly resolution). However, very little information about stream CO 2 dynamics exists from agricultural areas, a land use type that is heavily managed by humans, including hydrological drainage, nutrient additions, soil cultivation, etc. As a consequence, CO 2 patterns in agricultural streams could potentially be very different than in other land use types with amplified diel CO 2 dynamics due to high metabolism and/or quicker response to hydrological events due to effective drainage systems.
In addition to the concentration gradient between the stream water and the air above, gas exchange is also highly dependent on the physical conditions at the air-water interface. For stream systems, the gas transfer velocity (often the variable given to describe the efficiency of the air-water gas exchange) is related to a combination of hydrological and morphological conditions of the stream channel, often including slope, velocity, and water depth (Raymond et al., 2012;Wallin et al., 2011). All these variables are proxies for describing the turbulence of the stream water, which controls the gas exchange but is rarely directly measured . Agricultural areas are often located in flat landscapes, resulting in drainage systems that have a low gradient and are slow-flowing (Rhoads et al., 2003;Hughes et al., 2010), conditions that prevent effective air-water gas exchange (Hall and Ulseth, 2019). However, whether the elevated pCO 2 observed in agricultural streams is an effect of land-use-specific hydro-morphological stream conditions preventing efficient gas exchange or an effect of high internal (aquatic) or external (terrestrial) CO 2 production is currently unknown.
Although recent studies have identified agricultural streams as high-pCO 2 systems, there are still large knowledge gaps to be filled in order to improve our understanding concerning the influence of these waterbodies in landscape C cycling. Here we present high-resolution (hourly) CO 2 concentration measurements in a Swedish agricultural headwater stream for more than a year (in total 339 d excluding periods of ice and snow cover). The study aimed to (1) quantify CO 2 concentration levels in an agricultural stream and explore its temporal dynamics and (2) identify the main drivers causing temporal variability in stream CO 2 concentration and how they might vary with season.

Study area
The study was conducted within the 11.3 km 2 Sundbromark (SBM) catchment (59 • 55 N, 17 • 32 E), located 5 km NW of the city of Uppsala, Sweden (Fig. 1b). The catchment is a part of the hydro-meteorological observatory Marsta that was established in the late 1940s (Halldin et al., 1999). The 30-year (1960-1991) mean annual temperature for the area is 5.3 • C (mean January and July temperatures are −4.5 and 16.0 • C, respectively) and with a mean annual precipitation of 535 mm. The length of the growing season is on average ca. 210 d from early April to the end of October (Swedish Meteorological and Hydrological Institute, SMHI). The catchment is dominated by agricultural land (86 %) mainly used for cereal production and pasture, and with minor influence of forest (8 %) and urban areas (6 %). The area is flat with only 28 m elevation difference from 41 m a.s.l. at the highest point to 13 m a.s.l. at the catchment outlet (Table 1). The bedrock consists of gneissic granites, and the soils are dominated by post-glacial clay at lower elevations and with some influence of glacial clay and silt at higher elevations. Although the bedrock does not contain any known carbonates, the soils are alkaline due to glacial-carbonatecontaining deposits, resulting in a stream pH ranging between 7.4 and 8.4 (Table 2), and with high electrical conductivity (EC, 791-1908 µS cm −1 ) (Osterman, 2018). The nutrient and DOC levels of the stream water (Table 1) are at the lower end (within the 25th percentile) of monitored agri-   To explore how representative the SBM catchment is for streams draining agricultural areas in the region, a snapshot sampling survey was performed across 10 streams (denoted region UPP 2 in the study by Audet et al., 2020) of various sizes (catchment area 8.5-740 km 2 ) and agricultural influences (30 %-86 %) distributed within a radius of 10 km from the city centre of Uppsala (Fig. 1a, Table S1 in the Supplement).

Field sampling and analysis
The measurements were conducted from 26 September 2017 to 12 December 2018 (in total 339 d of measurements excluding periods of ice and snow cover). Stream CO 2 concentration was monitored using an eosGP sensor (Eosense, Dartmouth, Canada). The sensor was covered by copper tape in order to avoid biofouling. Sensor accuracy is <1 % of the calibrated range (0 %-2 % CO 2 ) +1 % of the reading, corresponding to a maximum error of ca. 0.3 mg C L −1 based on the maximum CO 2 measured in the current study. The CO 2 sensor was calibrated against known gas standards (400, 1000, 5000, and 20 000 ppm) before and after deployment. No significant drift (exceeding the above-given uncertainty) in the instrument was observed during the period. Volume fraction outputs from the sensor were corrected for variations in temperature and pressure (atmospheric and water depth) using the method described in Johnson et al. (2010) and expressed in milligrammes of carbon per litre.
Water level, water temperature, and EC were measured together with CO 2 concentration at a V -notch weir (Fig. S1 in the Supplement). Water level was measured using a pressure transducer (1400, MJK Automation, Sweden) mounted in a stilling well representing the stream water level at the Vnotch weir. Discharge was calculated from a stage-discharge rating curve based on a series of manual measurements and according to a rating curve presented in Holmqvist (1998). Water temperature and EC were monitored using a thermocouple (Type T) and a CS547A-L conductivity sensor (Campbell, UK), respectively. The sensors (except for the pressure transducer) were deployed under the water surface attached to a wooden rod in the centre of the stream just upstream of the weir. All sensors were connected to a CR1000X data logger (Campbell, UK) measuring at a 1 min interval and storing average values at a temporal resolution of 30 (in 2017) or 60 min (in 2018).
Stable isotopic analysis of the dissolved inorganic carbon (DIC) (δ 13 C-DIC) was performed on six occasions during the falling limb of the snowmelt discharge peak in 2018 in order to explore the temporal variability in DIC source. At each sampling occasion a sample for analysis of δ 13 C-DIC was taken in a 60 mL glass vial completely filled with stream water and closed airtight with a rubber septum below the water surface. In order to preserve the sample, 1 mL of highly concentrated ZnCl 2 solution was injected in each sample (with subsequent release of 1 mL of sample in order to keep atmospheric pressure) directly after sample collection. Samples were kept cold and dark until analysis. Prior to analysis, 2 mL of sample was injected into 12 mL septum-sealed pre-combusted glass vials (Labco Limited) pre-filled with He gas and pre-injected with 0.1 mL of concentrated phosphoric acid in order to convert all DIC species to CO 2 (g) (Campeau et al., 2017a). The samples were analysed using an isotope ratio mass spectrometer (DELTA V Plus, Thermo Fisher Scientific, Bremen, Germany) GasBench II (Thermo Fisher Scientific, Bremen, Germany) measuring the CO 2 in the headspace. Each sample was analysed seven times (sample volume; 100 µL per sample), and the first two injections for each sample were discarded to avoid memory effects. The mean of the other five samples was taken to give the final result. The δ 13 C-DIC values are given in terms of deviation from known carbonate standards in per mille, where R is the (1) Precipitation, air temperature, and incoming shortwave (global) radiation data (Fig. 2) were obtained from the Marsta meteorological observatory located within the catchment ca. 2.5 km from the stream sampling station (Halldin et al., 1999). In the absence of direct measurements of photosynthetically active radiation (PAR), shortwave incoming radiation was used as a proxy for available photosynthetic light.
A spatial sampling campaign for CO 2 concentration, pH, EC, and water temperature was conducted on 21 June 2018 across 10 agricultural streams (including the SBM stream) located in different catchments around the city of Uppsala (Fig. 1a). The sampling was performed between 10.00 and 14.00 during the day. Samples for CO 2 analysis were collected using the headspace method (Hope et al., 2004;Kokic et al., 2015). Briefly, 30 mL of bubble-free water was collected in 60 mL polypropylene syringes and equilibrated with a known volume of ambient air by shaking vigorously for 1 min. The equilibrated headspace (15-20 mL) was recovered and analysed on an ultra-portable greenhouse gas analyser (UGGA) (Los Gatos Research, USA) equipped with a soda lime filter and manual injection port. In situ CO 2 concentration was calculated from the UGGA-determined parts per million values using Henry's law considering stream temperature (Weiss, 1974), atmospheric pressure, the added ambient air, and the water-air volume ratio in the syringe. The pH, EC, and water temperature were measured in situ in the streams with handheld instruments, for pH with a pH110 pH meter (VWR, USA) and for EC and temperature with a HI 99300 (Hanna Instr., USA).

Delineation of the stream network and catchment characteristics
Catchment area and characteristics were calculated in QGIS 3.8 based on a high-resolution (2 m × 2 m) digital elevation model (DEM) derived from lidar data (GSD elevation data, grid 2+, Swedish Land Survey). Land use distribution within the catchment was derived from the CORINE Land Cover 2018 product (European Environment Agency), and soil and bedrock characteristics were based on digital versions of the Quaternary deposits (1 : 25 000-1 : 100 000) and bedrock (1 : 50 000-1 : 250 000) maps (Swedish Geological Survey).

Data analysis
Out of the total data set (339 d) from the SBM catchment, only data measured at discharge rates >0 L s −1 (i.e excluding standing water or completely dry conditions) were used in the analysis of the stream CO 2 data (Fig. S1). For further evaluation of the control on stream CO 2 concentration, the data set was divided into four periods (autumn, 49 d; snowmelt, 17 d; spring, 91 d; and the dry period, 138 d) according to distinct phases in the hydrograph (Fig. 3, Table S2). The stream CO 2 dynamics observed among the different periods were examined visually, and any hydrological controls on the CO 2 were identified by the presence and direction of CO 2 -discharge hysteresis loops (Evans and Davies, 1998). Similar hysteresis analysis was used to investigate diel patterns in the CO 2 concentration data. Spearman's rank correlation coefficient was used to test for monotonic relationships between the diel amplitude in stream CO 2 concentration and potential drivers. Correlations were considered significant if p<0.05. The software JMP 14.2.0 (SAS Institute Inc., Cary, NC, USA) was used for all statistical calculations.

Results
The mean air temperature and total precipitation for the entire period (26 September 2017-12 December 2018) were 6.8 • C and 704 mm, respectively. The summer and autumn of 2018 were dry with generally low precipitation; the exception was on July 29 with 82 mm of rain within 24 h (Fig. 2). Mean and median stream discharge for the study period were 30.6 and 0.9 L s −1 , respectively, and with a total range from 0 to 668 L s −1 (corresponding to a range from 0 to 5.0 mm d −1 ) (Fig. 3). However, due to a high water table exceeding the range of the pressure transducer, the absolute peak discharge occurring from 5 to 7 April was missed in the measurements. The large skewness between mean and median discharge was an effect of the large number of days without water flow over the weir during summer and autumn 2018, 128 d (38 %) out of the study period. According to frequency analysis, 67 % of the days had a mean daily discharge <5 L s −1 . Despite the few days with discharge >100 L s −1 (7 % of the entire period), those days accounted for 69 % of the accumulated discharge. The majority (84 %) of these high-discharge days occurred during the snowmelt in April.

General CO 2 patterns
The stream CO 2 concentrations during the entire study period (median and mean of 3.44 and 3.94 mg C L −1 , respectively, corresponding to a pCO 2 of 4778 and 5324 µatm) were highly variable (IQR = 3.26 mg C L −1 ) (Fig. 3) and displayed a bimodal distribution with frequency peaks at ∼ 2.7 and ∼ 6.1 mg C L −1 (Fig. S2). The lower peak was associated with the snowmelt and spring period, whereas the higher peak was attributed to the autumn period in 2017 and to rain events during the dry period of summer-autumn 2018. In addition to the bimodal shape, a very distinct peak in frequently measured concentrations was observed at ∼ 1.6 mg C L −1 . This peak was attributed to the minimum concentration values for the diel cycles observed during the spring period.

Controls on stream CO 2 concentration
The autumn period started dry with low discharge (<3 L s −1 ) for the initial month of measurements. The CO 2 concentrations were at the same time highly dynamic but unrelated to variations in discharge. The CO 2 concentration reached the maximum for the autumn (10.89 mg C L −1 , which was also the maximum for the entire study period) during late October followed by a decline in CO 2 to ca. 2 mg C L −1 in early November. During November and December four main rain events were identified which all displayed an increasing stream CO 2 concentration with increasing discharge. In three of these events a positive clock-wise hysteresis loop was observed (Fig. 4) where the CO 2 concentration reached its maximum before the discharge did. At the last event during autumn 2017, the relationship between CO 2 concentra- Figure 6. Time series of (a) stream CO 2 concentration (black) and discharge (red) and (b) water temperature (black) and shortwave incoming radiation (SR, red) covering the period April-July 2018. Note the reverse axis for shortwave incoming radiation.
tion and discharge was close to linear, but still positive. During the snowmelt period the hydrograph was characterized by a diel cycle with melting during daytime, resulting in daily discharge peaks which were suppressed during night-time freezing. In contrast to the autumn events, the daily discharge peaks were negatively related to the stream CO 2 concentration and with an anticlockwise hysteresis loop where the minimum CO 2 concentration was reached before the highest discharge of the event (Fig. 5). After the snowmelt discharge peak, the spring and early summer periods (late April to early July) were dry with limited precipitation and with a steady decline in runoff (Fig. 3). During this period the CO 2 concentration displayed a pronounced diel cycle with daily maximum and minimum CO 2 concentrations reached during early mornings (06:00) and late afternoons (18:00), respectively (Fig. 6). The medium amplitude of the diel CO 2 cycle for this period was 2.03 mg C L −1 , corresponding to pCO 2 = 2974 µatm (IQR = 1.23 mg C L −1 , corresponding to pCO 2 = 2212 µatm), and with the size of the diel CO 2 concentration amplitude being related to both the daily mean water temperature and the shortwave radiation (Fig. 7). The diel pattern displayed a clear negative anti-clockwise CO 2stream water temperature hysteresis loop, where the median CO 2 concentration could differ by up to 75 % between daytime and night-time despite being measured at the same stream water temperature (Fig. 8). From early July the stream dried out and hence no runoff over the V -notch weir was generated. During this period the CO 2 sensor mostly recorded an atmospheric signal. However, for five rain events during the summer and early autumn, runoff was generated, which allowed stream CO 2 determination for shorter periods (Fig. 9). During these runoff events (<2 d long), high CO 2 concentration pulses were recorded (up to 11 mg C L −1 ). At all events CO 2 was recorded for a longer period than the discharge as the small dam above the V -notch weir was still water-filled for some time after runoff over the weir ceased. Also, common for all events was that the stream CO 2 concentration continued to increase although the discharge peak had passed. During 29 July a heavy rainstorm occurred with 82 mm of precipitation in 24 h. Although more than 15 % of the long-term annual mean precipitation fell during 1 d, low discharge was generated (maximum discharge 6.1 L s −1 ) due to high evapotranspiration and dry soils (Figs. 3 and 9). However, the rainstorm event resulted in close to the highest stream CO 2 concentration (10.81 mg C L −1 ) being observed during the studied period. As soon as the stream was more permanently refilled in early December and with discharge generated over the weir, the stream CO 2 concentration was back to similarly high levels (typically 5-8 mg C L −1 ) as observed in the autumn of 2017.

Sources of DIC
The δ 13 C-DIC data collected during the falling limb of the spring discharge peak (discharge range 130-9.6 L s −1 ) ranged from −13.8 ‰ to −12.2 ‰. This narrow range suggests a relatively constant source of inorganic C during the spring period. Although there was a tendency towards more negative δ 13 C-DIC values at higher discharge, no significant relationship was found (Fig. 10). δ 13 C-DIC was also unrelated to the stream CO 2 concentration (data not shown).  . Stream CO 2 concentration (black) and discharge (red) for the dry period (July-September 2018). Periods when the CO 2 sensor was above the water table capturing an atmospheric signal (i.e. with concentrations <0.5 mg C L −1 ) are highlighted by the lower box.

Spatial representativeness
The 10 streams manually sampled around Uppsala displayed a wide range in CO 2 concentrations (1.8-4.6 mg C L −1 ) on the day of sampling (21 June 2018) and with the SBM stream (site 3 in Table S1) being close to the overall median (SBM, 2.7 mg C L −1 ; overall median, 3.0 mg C L −1 ) (Table S1). Furthermore, the CO 2 concentration manually sampled at SBM was close to the sensor-recorded CO 2 (2.59 mg C L −1 ) at the hour of sampling. The SBM stream was also close to the spatial median DOC concentration but slightly elevated in NO 3 and PO 4 . The CO 2 concentration was on a spatial scale related to pH but unrelated to catchment area or land use distribution within the catchment. Furthermore, the CO 2 concentration was on a spatial scale unrelated to mean stream concentrations of DOC, PO 4 , and NO 3 , although these variables were sampled during a different period than the CO 2 .

Discussion
In order to produce large-scale estimates of the exchange of greenhouse gases (GHGs) between inland surface waters and the atmosphere, a basic requirement is to know the aqueous concentrations of the gases of interest and how they might vary over time. Headwater streams have been identified as "hotspots" for CO 2 emissions (Raymond et al., 2013;Wallin et al., 2018), but there are limited data capturing the temporal resolution, specifically from streams draining agricultural regions, making large-scale generalizations uncertain. Due to effective drainage, high-nutrient conditions and often high sunlight exposure (due to limited tree cover), agricultural streams could potentially be very different in their CO 2 dynamics compared with streams draining other environments. Here we continuously measured stream CO 2 concentration in a headwater catchment dominated by agricultural land use (86 %) covering more than 1 year of the snowfree period. In line with findings from similar studies from other environments (arctic tundra, boreal forest, temperate peatlands, alpine areas) (e.g. Rocher-Ros et al., 2019;Riml et al., 2019;Crawford et al., 2017;Peter et al., 2014;Dinsmore et al., 2013), we found a mixture of controls on stream CO 2 operating at different timescales, generating a highly dynamic stream CO 2 concentration pattern. These timescales cover seasonal patterns to diel cycles, or even shorter scales associated with discharge events. Both the magnitude of CO 2 concentrations and their associated temporal dynamics were found to be high in the current agricultural stream when compared with the literature. The mean CO 2 concentration (3.94 mg C L −1 corresponding to a pCO 2 of 5324 µatm) is at the high end when compared with other high-frequency CO 2 data sets covering low-order (less than third stream order) catchments draining multiple environments, including arctic tundra, boreal forest, hemi-boreal forest, temperate forest, temperate peatlands, and alpine areas (typically ranging from ca. 0.2 to 6 mg C L −1 ) (Crawford et al., 2017;Natchimuthu et al., 2017;Peter et al., 2014;Dinsmore et al., 2013). Still, CO 2 concentrations in SBM do not seem to be exceptionally high compared to snapshot-based data from other agricultural streams.
The spatial variability seen in this study, although only based on snapshot samples, and previous studies indicates that CO 2 concentrations in agricultural streams are comparably high (Borges et al., 2018;Bodmer et al., 2016;Sand-Jensen and Staehr, 2012). In addition, the observed temporal dynamics presented here are, to our knowledge, among the most pronounced in the literature, although the number of high-frequency stream CO 2 data sets is limited, for example, the rapid decrease in stream CO 2 during the autumn of 2017, the strong diel cycle (diel amplitude up to almost 5.0 mg C L −1 ) during the spring-early summer period, or the rapid and high CO 2 pulses (up to 11.0 mg C L −1 ) occurring in accordance with rain events during the dry late summerautumn period. These high CO 2 dynamics clearly illustrate the need for continuous high-frequency CO 2 concentration measurements in streams in general, and in agricultural streams more specifically. Without such high-frequency data, representative estimates of agricultural stream CO 2 will be associated with high uncertainty. Although based on measurements from a single stream, these findings in turn indicate that current large-scale stream CO 2 emission estimates (e.g. Raymond et al., 2013;Humborg et al., 2010), which are largely based on snapshot concentration data with low (or no) resolution in time, might be specifically uncertain for agricultural regions.
According to our continuous data the highly dynamic pattern in stream CO 2 concentration is driven by a complex interplay of hydrology and biology. The high autumn concentrations observed in both 2017 and 2018 are likely an effect of high respiration of organic matter in the stream channel and/or in the adjacent soil water (Fig. 3d). This is supported by efficient aquatic microbial DOC degradation (<800 µg C L −1 d −1 ) observed during the autumn period across the 10 streams (agricultural land use, 30 %-86 %) included in the spatial sampling campaign (Peacock et al., unpublished data). This could be compared with organic C degradation rates determined in boreal forest and mire streams displaying typically lower rates (<300 µg C L −1 d −1 ; Berggren et al., 2009). The positive CO 2 -discharge relationships indicated that event flow pathways were in contact with soils with higher concentrations of CO 2 compared to flow pathways during base flow (Evans and Davies, 1998;Seibert et al., 2009). Also, the clockwise shape of the hysteresis loop suggests that there is a buildup of CO 2 in the catchment that is flushed out during rain events (Fig. 4). The CO 2 pool seems to be limited as the CO 2 concentration drops before the maximum discharge peak occurs, and vertical patterns in the CO 2 soil profile control the stream CO 2 depending on dominating flow paths (Evans and Davies, 1998;Öquist et al., 2009). This could explain the fact that the stream CO 2 increase did not reach any source limitation for rain events of lower magnitude (Fig. 4d). Similar positive CO 2 concentration-discharge patterns have been observed across different low-order streams (e.g. Crawford et al., 2017;Dinsmore et al., 2013), but the absolute patterns are often concluded to be highly site-specific and even event-specific. Here we suggest, by exploring the hysteresis loops, that such positive relationships are influenced by the size of the available catchment CO 2 pool or the hydrological connectivity to it. In a highly drained low-elevation agricultural landscape where much of the stream runoff is generated through drainage pipes (Castellano et al., 2019), the extent and spatial distribution of these terrestrial source areas and connections between groundwater and surface water are central for the CO 2 patterns observed in the stream. Strong hydrological control has been found for DOC in agricultural streams in the USA and France, where high-discharge events flush allochthonous DOC, via subsurface drainage pipes, into streams (Morel et al., 2009;Royer and David, 2005). In contrast to the seasonally variable CO 2 -discharge response patterns observed in the current study, Morel et al. (2009) suggested that stream DOC is non-limited and would continue to rise until the maximum discharge peak is reached. Whether this discrepancy in source limitation between CO 2 and DOC (although based on different studies and environments) indicates differences in the source areas of the different carbon components requires further investigation.
In contrast to the patterns observed during the autumn, during the snowmelt period the stream CO 2 was diluted when discharge increased following a diel pattern (Fig. 5). The melting and freezing between daytime and night-time suggest that meltwater from the surface snowpack during daytime to a larger extent reached the stream without picking up an elevated CO 2 signal. Similar dilution patterns in conjunction with snowmelt have been observed in catchments of various land use but specifically in peatland catchments with limited forest cover (e.g. Wallin et al., 2013). The similarity between this agricultural catchment and open peatlands could potentially be the effect of an efficient melting of the snowpack. Both unforested peatlands and agricultural fields are open areas subject to direct sunlight and wind and rain exposure, while the soil under the snow remains frozen. As a result, a large share of the meltwater will never infiltrate the soil but instead reaches the surface drainage system as overland flow (Laudon et al., 2007). This is further accompanied by the low hydraulic conductivity of clay soils, which dominate the catchment of the current study. Although we did not capture the 2-3 d of peak spring flood (due to a water level out of the range of the pressure transducer), it was evident that the stream CO 2 concentration was diluted from ca. 6.0 to ca. 2.0 mg C L −1 during these days, something that is further supported by the similar drop in EC during the peak spring flood from ca. 900 to ca. 150 µS cm −1 . However, as soon as the discharge peak passed, the stream CO 2 concen-tration recovered rapidly to the pre-peak levels, suggesting a shift to hydrological pathways that mobilize a high CO 2 pool, again supported by the concurrent increase in EC. April and May 2018 were characterized by warm and clear weather with an average 4.2 • C higher air temperature and 255 more sun hours than the 30-year mean (1961( -1990. Altogether, this stimulates a kick-start of the aquatic primary production upon snowmelt, which likely explains the steady decline in CO 2 that occurred during late April-early May. During the spring and early summer, a strong diel pattern in CO 2 concentration further developed, likely driven by aquatic primary production consuming CO 2 during daytime. Such diel CO 2 patterns are commonly observed in stream CO 2 time series at base flow or during receding flow conditions (e.g. Riml et al., 2019;Peter et al., 2014) and are especially pronounced in amplitude in nutrient-rich streams or in streams without canopy shading (Alberts et al., 2017;Crawford et al., 2017;Rocher-Ros et al., 2019). Initial evaluation of the δ 13 C-DIC data collected during the spring period suggests a relatively steady mixture of geogenic and biogenic DIC although somehow related to variations in discharge (Fig. 10). However, given the suppressed stream CO 2 during the spring period, together with the strong diel cycle caused by aquatic primary production, fractionation of a strict biogenic DIC pool (with a δ 13 C-DIC from −28 ‰ to −20 ‰) could theoretically push the δ 13 C-DIC towards the less negative values observed in the current study (from −13.8 ‰ to −12.2 ‰) (Campeau et al., 2017b). Combined studies on aquatic metabolism, C dynamics, and stable isotopic composition is recommended to disentangle the dynamic CO 2 source patterns in this type of agricultural system.
The spring and early summer of 2018 were generally dry, leading to the stream channel drying out during long periods. The rapid rewetting periods (<2 d) that occurred following larger precipitation events resulted in high CO 2 pulses (3-11 mg C L −1 ), generally exceeding the overall median level of stream CO 2 (3.44 mg C L −1 ) observed during the study period. The intermittent nature of streams, with distinct drying and rewetting episodes, is known to generate high CO 2 concentration pulses and subsequent emissions (Marcé et al., 2019). Such rapid pulses are generally suggested to be a result of intense respiration in the stream bed sediments upon rewetting, or due to a rapid mobilization of terrestrial C, both organic (DOC) and inorganic (CO 2 ) in connection to precipitation events. However, the high CO 2 pulses upon rewetting have mostly been found in areas that display pronounced dry and wet seasons, e.g. Mediterranean areas or Australia (e.g. Gómez-Gener et al., 2015;Looman et al., 2017). Here we show that such stream intermittency can also cause high and rapid CO 2 pulses in a Swedish agricultural setting, highlighting the need for expanding the geographical coverage of studies that investigate stream intermittency in relation to GHG dynamics and emissions. Areas that display stream intermittency will likely also increase in the future given the predicted changes in temperature and precipitation patterns.
An obvious tool in this work is the use of continuous sensorbased measurements which allows the capture of the episodic and unpredictable nature of these phenomena.

Conclusions
It is evident from the current study that the stream CO 2 dynamics in an agricultural headwater catchment are highly variable across a variety of different timescales and with an interplay of hydrological and biological controls. The hydrological control was strong (although with both positive and negative influences dependent on season) and rapid in response to rainfall and snowmelt events. However, during growing-season base flow and receding flow conditions, the aquatic primary production seems to control the stream CO 2 dynamics, which in turn sets the basis for atmospheric emissions. During the dry summer period, rapid rewetting following precipitation events generated high CO 2 pulses exceeding the overall median level of stream CO 2 (up to 3× higher). This finding thus highlights the importance of stream intermittency in agricultural areas and its effect on stream CO 2 dynamics. Given the observed high levels of CO 2 and its temporally variable nature, agricultural streams clearly need more attention in order to understand and incorporate these considerable dynamics in large-scale extrapolations.
Author contributions. MBW and MW conceived the idea and designed the study. MBW funded and instrumented the catchment and analysed the data. MW conducted the GIS analysis. JA, MP, and ES provided ideas and data. MBW wrote the manuscript with great support from all co-authors.