A Multi-Disciplinary Approach to Understand Hydrologic and Geochemical Processes at Koiliaris Critical Zone Observatory

: Koiliaris CZO is a European Critical Zone Observatory (CZO) typical of the Mediterranean karstic geomorphology, which represents watersheds a ﬀ ected by humans over the centuries. This study aims to provide information that underpins the hydrologic and geochemical processes functioning at Koiliaris CZO. Linking geomorphologic and tectonic analysis improved the delineation of a karstic area which extends outside of the Koiliaris watershed and identiﬁed how structural elements inﬂuence the regional hydrology. The ﬂuctuation in the river ﬂow represents processes occurring in the karst and the periodic signal is related to Earth tide stressing of the karstic reservoirs. The conceptualization of a two-reservoir, well-mixed karstic system is conﬁrmed by both the geomorphologic and tidal analysis. The hydrologic response is fast and it is manifested especially during extreme events where 70% of the precipitation becomes surface runo ﬀ , creating major ﬂood events. The di ﬀ erent sampling sites in the Koiliaris CZO were geochemically clustered and the quantiﬁcation of the weathering ﬂuxes showed that 25 mm / 1000 years and 39 mm / 1000 years of carbonate were removed by chemical weathering for the Keramianos ephemeral river and the springs, respectively. These studies illustrate the importance of critical zone science and transdisciplinary studies on water and soil management.


Introduction
The main anthropogenic drivers causing ecosystem change and potential tipping points are climate change, land-use change, and above-and below-ground biodiversity loss [1]. In many aspects, the carrying capacity of the planet has been exceeded with potential adverse consequences [2,3]. The degree of utilization of ecosystem resources is such that almost two thirds of ecosystem services have been degraded in just fifty years [4]. In addition, humans have to mitigate and adapt the impacts of climate change that would exert a significant pressure on all ecosystems around the globe and the services they provide [5].
One of the strategies the European Commission is using to address these global environmental challenges is through long-term ecosystem research (LTER) and the development of a pan-European research infrastructure. The European LTER network conducts "integrated long-term research on Table 1. Research Areas and Key Scientific Questions Addressed at the Koiliaris CZO and research citations related to them.

Stream and Ground Water Chemistry [16-19]
• How anthropogenic practices affect hydrogeochemical pathways between the terrestrial and aquatic environments in a Mediterranean watershed • What are the impacts of agricultural land use practices on water quality? • What are the impacts of livestock grazing on water quality? • Identify hydrogeochemical pathways and quantify chemical fluxes between the terrestrial and aquatic environments

Soil Degradation and Formation [11,20-31]
• How grazing affects nutrient (C, N, P) cycling in degraded soils of Koiliaris CZO highlands. What are the impacts of livestock grazing on soil and water quality? • What are the processes controlling aggregation and nutrient sequestration in soils? • How agricultural practices affect aggregation and nutrient sequestration in soils? • What are the factors and processes affecting soil (aggregate) formation and degradation? Rates of formation and degradation? • What is the areal extend of soil at threat to degradation?

Hydrologic and Geochemical Modeling of Complex Terrain [16,17,20,32-38]
• Model parameterization to improve hydrologic and geochemical predictability • Scaling up methodology -what are the uncertainties of parameter estimation at different scales.

Sustainable Management of Water and Soil Resources [1,12,39-47]
• How is climate change going to affect water availability and management? • Soil sustainability requires long-term social change. What is the best governance paradigm for such a change? • Considering the fact that sustainability demands a mixture of governmental policies, market initiatives and civil-society interventions what are the more appropriate tools and mechanisms for more efficient public participation?

High Frequency Environmental Monitoring [15,48]
• Development of intelligent monitoring system with high sampling frequency and how this is improving the data availability for environmental decision making? • How is high frequency sampling going to improve the data availability flash flood forecasting?

Site Description
The Koiliaris CZO is situated in the north-western part of Crete, near the city of Chania, Greece and extends within the latitude 35 • Figure 1). The climate of the area is characterized as Mediterranean, semi-arid. The total area of the watershed is 132 km 2 and the altitude ranges from 0 to 2120 m above mean sea level (amsl). Slope range from 43% at high altitude to 1-2% at lower altitude. The main land uses are intensively grazed shrubland and pasture (67.3%); olive, citrus groves, vines, and vegetables (32.1%); and mixed forest (0.6%). The length of the drainage network is 44.8 km, consisting of the intermittent tributary of Keramianos (13.8 km), two ephemeral streams providing surface runoff, located in the eastern part of the basin and feeding the Anavreti tributary (total length is 27.2 km), and, finally, the karstic springs of Stylos (permanent flow) (S1) and the karstic spring of Anavreti (intermittent flow) (S5), located just upstream of the Agios Georgios hydrometric station (R1), which merge with the rest of the streams to form the main segment of the Koiliaris River ( Figure 1). These springs are fed by an extended area of karst which is located outside the basin boundaries and occupies an estimated area of 80 km 2 . The extended karstic area has been approximated following the geomorphologic characteristics and the faults formation and direction [15,16]. One of the aims of this work is to improve the delineation of the extended karst. The Keramianos tributary drains a small sub-catchment in the eastern part of the basin, enters a karstic gorge and joins the remaining of the river and springs network. The Keramianos tributary loses most of its water in the two faults that crosscut the gorge and generates flash floods when the precipitation in the sub-basin exceeds 120 mm [17].

Geological Setting
The geology of the Koiliaris CZO is composed of alpine nappes and post-alpine rocks. Specifically, the predominant geologic rock types are: limestones, dolomites, marbles, and recrystallized limestones with cherts of the Plattenkalk; Tripolis; and Trypali series; marls of Neogene deposits [49]; and schists (Arna unit). Quaternary alluvial deposits crop out in the northern part of the Koiliaris basin ( Figure 1) [16].
The karstic system of the White Mountains ( Figure 1) is composed of an autochthonous geotectonic unit, the Plattenkalk nappe and two allochthonous units, the Western Crete unit (Trypali limestones and phyllites) and the Arna unit (metamorphic schists-denoted as crystalline schists in Figure 1). the alpine nappes, the Plattenkalk nappe stratigraphy from the lower to the upper strata consists of (1) karstified un-bedded dolomites, (2) dolomitic marbles, (3) thin-bedded calcitic marbles with chert layers and nodules, and (4) metamarls and limeschists. The Trypali nappe consists of recrystallized brecciated un-bedded limestones thrusted over Plattenkalke nappe and it is highly karstified. In the western part of the basin, Western Crete Phyllites-Quartzites [49] cover the Trypali limestones. Alluvial sediments found near the river corridor and neogene deposits at the lower elevations of the basin comprise the post alpine rocks.
The major tectonic features in Crete comprised of primary structures (Late Oligocene-Early Miocene) of mega-folds with hinge of E-W trending and thrusting of nappes southward [49]. The primary structures are mainly the result of thrusting of "inner" nappes (Trypali, Western Crete Phyllites-Quartzites, etc. units) onto the autochthonous of Plattenkalk unit. Later tectonic contacts (Miocene) compose of extensional detachments (dipping at north and south with 10 • -30 • ) and normal faults (Middle-Upper Miocene-present) of higher angle which possibly dissect the detachment thrust faults [49][50][51] (Supplementary Figure S1).

Data Collection
The White Mountains, located in the south of the watershed, provide the main supply of water for the basin. Meteorological data are acquired by two meteorological stations within the basin, installed and managed by the Laboratory of Hydrogeochemical Engineering and Remediation of Soil The hydrologic monitoring of Koiliaris CZO has also been carried out since 2004, with the installation of the hydrometric station at Agios Georgios (R1), strategically located downstream of the confluence of Keramianos tributary with the main, permanent river that drains the Stylos (S1) and Anavreti (S5) springs, recording water level data every 10 min (Water level & Temperature logger, Level TROLL 500 Data Logger, In-Situ Inc., Fort Collins, CO, USA). H.E.R.S. Lab has also installed two hydrometric stations ( Figure 1) along the Keramianos tributary, one at the gorge entrance and one at the gorge exit,  Figure 1. Schematic of the regional geology, the drainage network and the monitoring network of Koiliaris CZO. Figure 1. Schematic of the regional geology, the drainage network and the monitoring network of Koiliaris CZO.

Data Collection
There are 9 sampling locations within the basin where monthly water sampling of surface and groundwater is taking place. The sites are: Agios Georgios stream (R1) located at the junction of Stylos karst springs and Keramianos, Katochori spring (S2) located near the entrance of Keramianos gorge and Katochori temporary river (R2), Stylos spring (S1) located at the center of Stylos village, Dug well (G2) located at Stylos village, Anavreti spring (S5) located near Nio Chorio village, Zourpos spring (S4) located at Kalives village, Armenoi spring (S3) located at Armenoi village and finally, Kampoi well (G3) located near Kampoi village (Figure 1). Field and laboratory measurements are conducted once a month. The field studies include the measurement of pH, electrical conductivity, water temperature and dissolved oxygen.

Geologic and Geomorphologic Analysis
A geomorphological and geological analysis was conducted to better define the area of the extended karst and understand how structural elements influence the regional hydrology. The southern boundary of the extended karst was located using the south-north divide of the White Mountains [16]. However, recently obtained cave and sinkhole locations [52] by various speleological expeditions necessitated the reexamination of the delineation the extended karst. The Digital Elevation Model (Shuttle Radar Topography Mission, SRTM, data with a pixel size 90 × 90) was used to generate hillshade (Supplementary Figure S1) and aspect maps [53]. The geological maps of the area were used for the digitization of the sedimentary rock dipping direction and the construction of the geological cross sections [54][55][56]. Geologic cross sections were created in order to identify major morphological and tectonic elements such as anticlines and mega folds and combine them with normal faulting in order to reassess the extended karst of Koiliaris CZO.

Water Level Frequency Analysis
Since the beginning of the study at Koiliaris CZO, the diurnal fluctuation of flow was observed using the high frequency flow measurements at the Agios Georgios gauging station (R1). Figure 2 illustrates this sub-daily variation of water level that consists of a persistent double-peaked daily cycle, especially distinguishable during low flow periods when only the Stylos springs feed the river. The amplitude of this periodic signal (2-4 cm) is orders of magnitude lower than typical flood events (50-250 cm, Supplementary Figure S3) triggered by rain events and it an order of magnitude lower than the overall long-term variability (average stage 2007-2017 is 60 cm) observed in streamflow. Moraetis and co-workers [15] has shown that this fluctuation in the river flow represents processes occurring in the karst. In addition, they used stage and pressure measurement within the karst to postulate that the "transfer" of the atmospheric barometric pressure through the sinkholes within the karst is a potential reason for the diurnal fluctuation and causes the karst to act as a "bladder pump". This could be one reason to support the hypothesis that the water in the karstic aquifer is well-mixed.
Matlab code of Pawlowicz et al. [70], keeping constituents that have a sufficiently large signal to noise ratio (SNR > 2). T_TIDE performs harmonic analysis and identifies tidal signals having at least 95% CL by using nonlinear bootstrapped error estimates (fitting a sum of sines and cosines with known frequencies to the data, using a multi-linear least-square regression). The coefficients that minimize the misfit are the amplitudes and phases for all selected tidal constituents [70]. The analysis was done for each period and also for both ER & LR sub periods with the original water level data.

Hydrologic Analysis
The objective of the hydrologic analysis was to confirm the hydrologic budget and hydrologic pathways using the new delineation of the extended karst and gain better understanding of the hydrologic pathways and response of the karst during extreme events. The Karst-SWAT version of the soil and water assessment tool (SWAT) was used to simulate the hydrologic and geochemical response of the basin. It was chosen because of its ability to capture the water cycle of multiple catchments of different climates, according to numerous studies and, particularly, its good performance in the Mediterranean region [16,32,39,71,72]. The SWAT model is a deterministic and continuous time basin scale model which simulates the hydrology, sediment yield and water quality of ungauged watersheds [73] and estimates the impact of agricultural management practices on agricultural yields. The main model inputs consist of a digital elevation model (DEM), a land use map, a soil map and daily meteorological data of precipitation and temperature. The watershed is first subdivided into sub-catchments and each sub-catchment into hydrologic response units (HRUs), which are a unique combination of soil type, land use and land slope. The model output provides, Diurnal fluctuations of river flow have been observed in many locations worldwide and several mechanisms may explain them [57]; (a) the diel fluctuation in mobilized water caused by freezing and thawing in cold regions; (b) "early afternoon rainfall" events that are frequently observed under tropical climates; (c) temperature-controlled variations in streambed hydraulic conductivity; (d) the diel evapotranspiration and snow melt signals which causes oscillations in groundwater table; (e) anthropogenic water abstraction for drinking and irrigation demand; (f) fluctuations due to tidal processes.
In the case of the Koiliaris CZO, it is clear that the two first phenomena cannot explain the diurnal signal due to the Mediterranean climate of Crete. In addition, water temperature within the karstic aquifer is fairly constant around 13 • C and there are no transmission losses in the riverbed from the springs to the sea. Therefore, we do not expect any losses or large flow fluctuations due to evaporation which in any case would not produce a double-peaked signal in the flow. The effect of water abstractions can also be excluded since the pattern persists throughout the year.
The only plausible explanation left is that of tidal effects which in fact, form a whole family of processes caused by the gravitational pull of the Moon and the Sun [58]. The cyclic interaction between these celestial bodies and the rotating Earth result in the periodic motion of oceans, as well as the atmosphere and the crust. These motions, in turn, produce mechanical perturbations at the Earth's surface that may affect continental water bodies in various ways and at different timescales, from a few hours up to several months [58][59][60]. For instance, tidal loading effects on groundwater reservoirs have been recognized and studied in the early developments of modern hydrogeology [61][62][63], and even have been used to characterize reservoir physical properties like permeability and storage [64][65][66]. Indeed, in confined aquifers, the dynamic stressing caused mainly by body tides (also named Earth tides) induce rock deformation, and due to the poroelastic coupling between the porous material and water, tidal signals are frequently observed in groundwater levels. Among the very many tidal constituents that may be measured in geophysical signals [67], the most significant in hydrology are the sinusoidal forcing with periods around 24 and 12 h, corresponding to the diurnal and semi-diurnal main astronomical phenomena.
Tidal effects in streamflow are most often encountered in coastal areas where ocean tides bring about a significant oscillation of the base level at a basin's outlet. These so-called "tidal rivers" undergo double-peaked daily tidal fluctuations with amplitudes ranging from a few centimeters to more than 10 m [68]. As the Agios Georgios gauging station is located at an elevation of 15 amsl and given that the maximum tidal range in Crete is less than 0.1 m, the sea tides cannot directly influence streamflow in Koiliaris. Consequently, the remaining possible mechanism is the stressing of the karst's porous matrix by Earth tides and/or ocean tidal loading that creates an oscillating hydraulic gradient and thus an oscillating spring flow. Such a mechanism to explain double-peaked daily fluctuations in stream hydrographs has very rarely been documented, as argued by Briciu [69]. In fact, Briciu [69] reported many occurrences of diurnal and especially semi-diurnal components in river flow data across the globe, and noticed that watersheds formed by limestones and dolomites were particularly prone to this effect. As the dominant Plattenkalk formation in Koiliaris is mainly composed by these rocks, it strengthens the hypothesis of a body tide dominated process.
To test and verify this hypothesis, a frequency and tidal analysis was performed using the stage measurements at Koiliaris CZO. The data were collected by using vented pressure transducers which give a direct value of the water level with the atmospheric pressure removed. Selected periods of the streamflow hydrograph at the Agios Georgios station are first analyzed using Continuous Wavelet transform (CWT). Prior to the CWT analysis, the data trend was estimated and removed by applying a 144-term (no of measurements per day) centered moving average. For CWT, the classic Morlet mother wavelet is used with 16 voices per octave. The tidal analysis is performed with the T_TIDE Matlab code of Pawlowicz et al. [70], keeping constituents that have a sufficiently large signal to noise ratio (SNR > 2). T_TIDE performs harmonic analysis and identifies tidal signals having at least 95% CL by using nonlinear bootstrapped error estimates (fitting a sum of sines and cosines with known frequencies to the data, using a multi-linear least-square regression). The coefficients that minimize the misfit are the amplitudes and phases for all selected tidal constituents [70]. The analysis was done for each period and also for both ER & LR sub periods with the original water level data.

Hydrologic Analysis
The objective of the hydrologic analysis was to confirm the hydrologic budget and hydrologic pathways using the new delineation of the extended karst and gain better understanding of the hydrologic pathways and response of the karst during extreme events. The Karst-SWAT version of the soil and water assessment tool (SWAT) was used to simulate the hydrologic and geochemical response of the basin. It was chosen because of its ability to capture the water cycle of multiple catchments of different climates, according to numerous studies and, particularly, its good performance in the Mediterranean region [16,32,39,71,72]. The SWAT model is a deterministic and continuous time basin scale model which simulates the hydrology, sediment yield and water quality of ungauged watersheds [73] and estimates the impact of agricultural management practices on agricultural yields. The main model inputs consist of a digital elevation model (DEM), a land use map, a soil map and daily meteorological data of precipitation and temperature. The watershed is first subdivided into sub-catchments and each sub-catchment into hydrologic response units (HRUs), which are a unique combination of soil type, land use and land slope. The model output provides, among others, the surface runoff for each HRU, but also the deep aquifer recharge, which is useful for the modeling of the karstic springs, and is another reason for the SWAT model selection.
Moraetis et al. [15] concluded that the karstic system of the area consists of two underground reservoirs: the upper reservoir with a faster response and the lower reservoir with a slower response. Based on this information, Nikolaidis et al. [16] augmented the SWAT model by adding a karstic component to account for the karstic springs flow. The SWAT model was calibrated in a way that allows fast infiltration of water that falls on the "karstic" HRUs into the deep groundwater. The model uses the deep groundwater flow from SWAT as input into a two-reservoir model and analytically solves hydrologic mass balance equations on a daily step. The model output is the daily total flow of the karstic springs (Stylos and Anavreti).
This coupled model allows for the determination of every hydrologic parameter for each HRU, and also enables the estimation of the hydrologic pathways during average and extreme conditions. An additional analysis was conducted to assess the flow of Anavreti and Stylos separately. Based on observations, the intermittent Anvreti spring discharges only if the groundwater stage exceeds a certain level, therefore this spring usually only runs during the months of December, January, February, and March and rarely during April. To define the point for the Anavreti spring to flow, the observed monthly flows of the Stylos springs for the 1973-2004 period were correlated to the total modeled springs flow (simulation of Stylos and Anavreti) for the above-mentioned months. The correlation was satisfactory (r 2 = 0.56, Q Stylos = 0.48 × Q spr ) and provided an estimation of the Stylos springs discharge, based on the modeled karstic flow, as well as an estimation of the Anavreti spring discharge

Geochemical Analysis
The objective of the geochemical analysis was to assess the geochemical time series of 9 surface and groundwater stations and determine the geochemical fluxes and weathering rates as well as the processes affecting the geochemistry. To achieve this objective, a principal components analysis (PCA) was conducted using the time series of the physico-chemical parameters in order to categorize and cluster the sampling sites as well as determine the parameters that differentiate their clustering. The carbonate weathering rates were calculated in Keramianos ephemeral river and in the springs of Stylos and Armenoi of Koiliaris CZO. Carbonate weathering rates (F carb ) were calculated by multiplying the concentration of carbonate-derived chemicals (Ca, Mg, and half of HCO 3 − ) with the water runoff value as follows [74]: where F carb is expressed in tn/km 2 /year or tn/km 2 /month, and M CaCO3 and M MgCO3 refer to the molar mass of CaCO 3 andMgCO 3 and are expressed in g/mol. The carbonate weathering rates were related with runoff, since runoff exerts a strong control on carbonate weathering rates [75] and were compared with the rates obtained from various catchments in Jura Mountains in the Alps.

Geomorphology and Delineation of the Extended Karst
The composite geomorphologic map in (Supplementary Figure S1) presents critical geomorphologic and structural elements of the extended karst of Koiliaris River. It presents the aspect map (overlain by hillshade with 60% transparency) along with structural (north and south dipping direction arrows) and lithology data (extent of Mesozoic Limestones blue dashed line). In addition, the location of caves and sinkholes (black dots in Figure 1 and yellow dots in Figure S1) are depicted.
The stratigraphic column of White Mountains is composed of gray-black dolomites, recrystallized limestones, limestones with silex and limeschists (Plattenkalk Unit). The area includes significant meso-to mega-scale folds [76,77].
Several publications [49,50,78,79] demonstrate the morphologic and tectonic conditions of Crete which are summarized in four phases: i.
Eocene to Early Oligocene east to west compression created the anticlines striking NNE-SSW ii.
Late Oligocene north-south compression created the anticline striking E-W iii.
Early to Middle Miocene crustal extension created extensive detachment faults and Middle to Upper Miocene grabens filled with shallow marine sediments iv.
Upper Miocene-Quaternary created normal faulting with three sets of striking directions faults E-W (or ENE-WSW), NE-SW and NW-SE.
Manutsoglu and co-workers [80] identified a mega-anticline underlining the Samaria Gorge (shown in Supplementary Figure S1 with number 1) which was created during the Eocene to early Oligocene east-west compression. In this study, we identified two other mega-anticlines that regulate the hydrogeology in the northern part of White Mountains. Anticline number 2 is situated in the boundary between north and south aspect of White Mountain flanks with a strike direction NE-SW. The hinge of this anticline can be considered as the divide of White Mountains between north and south aspects flanks. Point 2B (Supplementary Figures S1 and S2b) shows the estimated position of the flank aspect change. We postulate that this anticline was created during the Late Oligocene north-south compression. Along this anticline, 739 caves and sinkholes are situated in the gray-black dolomite (Supplementary Figures S1 and S2b and black dots in Figure 1). The caves are mainly spread in elevation between 1700 and 2200 m where the outcrop of gray-black dolomites is located (Supplementary Figure S1). The presence of the extensive system of karstic caves is related possibly to the high deformation around the hinge of the mega-fold (point 2B). The second mega-fold anticline 3 is placed in point 1A and it is shown in profile A (Supplementary Figures S1 and S2b). The hinge is parallel to Samaria gorge anticline, NNE-SSW. The area "Kaka kastelia" is situated near point 1A and has a high number of karstic caves within the gray-black dolomites similarly to point 2B (Supplementary Figures S1 and S2a). Thus, both points 2B and 1A are areas where the three lithological and geomorphological conditions are met: the older Plattenkalk autochthonous gray-black dolomites are exposed (anticlines), most of the karstic caves are situated and a change in the aspect is observed. Anticlines 2 and 3 identify the southern and eastern boundaries of the extended karst. In this respect, the delineation of the extended karst southern boundary is redrawn to include the karstic cave system along the anticlines 2 and 3. Supplementary Figure S1 confirms that the mega-structures are consistent with the overall dipping of the autochthonous sedimentary beds (Plattenkalk). It is clear, especially for anticline 2 that the south dipping beds (red arrows in Supplementary Figure S1) are situated southern than anticline 2 with the overall aspect to be directed south (greenish area in Supplementary Figure S1 Figure S2b,c). Therefore, the extended karst area of Koiliaris River is redrawn based on this analysis and it is presented in Figure 1 and Supplementary Figure S1. The karst area was extended in the southern border by including the caves and sinkholes along the anticline 2 and reduced in the northern border by eliminating the Neogene grabens. The total area of the extended karst remains the same (80 km 2 ) as was estimated from hydrologic modeling [16,17]. Table 2 presents the selected periods for the tidal analysis using the stage data of Koiliaris River from Station R1. Supplementary Figure S3 presents Agios Georgios station tidal analysis for the studied periods 2018-2019. The time periods were chosen where the main source of river flow was Stylos springs and not rainfall events, associated flash floods and the Anavreti spring. In this way, the presence of multiple spikes in the sample data, which may partly hide the tidal signals was avoided. Consequently, six time periods in 2007, 2013, 2014, 2017, 2018, and 2019 were selected. Each dataset was subsequently separated into two key time frames that where identified in the hydrologic analysis: the first corresponds to the early recession controlled by both lower and upper karstic reservoirs (period ER for "early recession") and the second corresponds to the late recession where the flow contribution originates exclusively from the lower reservoir (period LR for "late recession"). The presence of persistent diurnal (~24 h) and semi-diurnal (~12 h) oscillations in the water level was confirmed by using the continuous wavelet transform CWT (Figure 3). The CWT analysis clearly shows the existence of the two key time frames identified from the hydrologic analysis. This is seen from the variation in the amplitude of the oscillating signals (especially the diurnals). This amplitude change coincides with the shift from the ER to the LR time frame and is present in all datasets, with the exception of 2013 which doesn't include the ER period. The duration of each period within the dataset is reported in Table 2. The duration of the ER period ranges between 0 and 79 days and the LR period between 25 and 125 days. An example of the tidal constituents for the ER and LR period of the 2019 and 2018 data is presented in Figure 4 and Supplementary Figure S4, respectively. The harmonic analysis on all datasets confirms the presence of persistent diurnal (~24 h-P1 and K1) and semi-diurnal (~12 h-S2, EPS2 and L2) components in the streamflow (Table 3). Table 3 presents the tidal analysis for period 2019 as an illustration. The analysis also revealed several harmonics of the main tidal components, at periods of about eight (MO3 and SK3), six (SN4, MS4, and S4), and four (2SK5, 3MK7 and M8) hours. These are typical secondary constituents in signals affected by tidal perturbations. The amplitude of the diurnal components ranged from 2.9 mm to 8.7 mm and the semi-diurnal from 0.6 to 1.7 mm. The secondary harmonics amplitude ranged from 0.2 to 0.7 mm. The tidal analysis was also performed separately for the ER and LR periods and the results are presented in Supplementary Tables S1 and S2 for 2019 and 2018 periods. The analysis verified the increase in amplitude for most of the tidal signals during LR period which was visible in the CWT analysis. The amplitude of the diurnal tidal components (K1, P1) are highest, closely followed by the semi-diurnal component S2, consistent for both ER and LR. The amplitude of other harmonics also decreases with period. Another interesting finding is that the amplitude ratio of the two major diurnal constituents K1 and P1 (K1amp/P1amp) is ~0.33 in all datasets and in both ER and LR time frames. Lastly, most of the tidal signals in LR have a leading phase compared to those in ER. The harmonic analysis on all datasets confirms the presence of persistent diurnal (~24 h-P1 and K1) and semi-diurnal (~12 h-S2, EPS2 and L2) components in the streamflow (Table 3). Table 3 presents the tidal analysis for period 2019 as an illustration. The analysis also revealed several harmonics of the main tidal components, at periods of about eight (MO3 and SK3), six (SN4, MS4, and S4), and four (2SK5, 3MK7 and M8) hours. These are typical secondary constituents in signals affected by tidal perturbations. The amplitude of the diurnal components ranged from 2.9 mm to 8.7 mm and the semi-diurnal from 0.6 to 1.7 mm. The secondary harmonics amplitude ranged from 0.2 to 0.7 mm. The tidal analysis was also performed separately for the ER and LR periods and the results are presented in Supplementary Tables S1 and S2 for 2019 and 2018 periods. The analysis verified the increase in amplitude for most of the tidal signals during LR period which was visible in the CWT analysis. The amplitude of the diurnal tidal components (K1, P1) are highest, closely followed by the semi-diurnal component S2, consistent for both ER and LR. The amplitude of other harmonics also decreases with period. Another interesting finding is that the amplitude ratio of the two major diurnal constituents K1 and P1 (K1amp/P1amp) is~0.33 in all datasets and in both ER and LR time frames. Lastly, most of the tidal signals in LR have a leading phase compared to those in ER. From these remarks along with the accumulated knowledge on the hydrologic functioning of the Koiliaris CZO, we postulate three important implications: (a) The periodic signal is definitely related to Earth tide stressing of the karstic reservoirs, as it affects similarly the lower and upper karst. (b) The upper and lower karst have different storage properties, which reflects into the differences in tidal signal amplitudes from ER to LR. The hydromechanical effect of tidal stressing produces a significantly higher discharge oscillation during ER, which indicates that the upper reservoir is more compliant (in a mechanical sense) and therefore more capacitive [81]. This observation is based on the diurnal fluctuation of discharge. A small change in the stage during the ER period creates a large change in the discharge of the river because of the non-linear effect of the rating curve. In theory, this prevails for the porous or fissured karstic matrix that is connected to main conduits and ultimately to the springs. (c) The presence of multiple harmonics of the diurnal and semi-diurnal body tides is intriguing, especially for the fact that the 8-h signal is more prominent than it would be expected from Earth tide loading only. It is possible that even a small ocean loading (ocean mass bending the crust by its weight, not to be mixed with tidal effects in estuaries) may stress the karst on top of the preponderant body loading. Indeed, ocean tides typically present overtones, also called "overtides" [82]. Another possible mechanism is the presence of nonlinearities in the hydromechanical coupling. Agnew [82] showed that if a nonlinear relationship including energy dissipation would characterize the stress-strain relationship in a rock, then an anomalous tidal line should appear at 3 times the main frequency. In our case, this coincides with the prominent eight-hour harmonic, which is even more pronounced during LR. This tells us that body tides might cause very large strains on the rock mass (probably more than 10 nstrain, [82]) at Koiliaris, creating hysterical loops in stress-strain relationships and making the overall mechanical behavior depart from the classical linear poroelasticity framework [81].

Flow Frequency Analysis
Summarizing, the tidal analysis confirmed the partitioning of the karst into two reservoirs that have different hydromechanical properties and behavior, which is in agreement with previous studies [15]. The analysis also underlined the fact that both reservoirs have a well-developed secondary porosity able to feed major fault or karstic conduits and act as a quickly accessible storage volume. Finally, the tidal loading by creating an oscillating hydraulic gradient and an oscillating spring flow on a sub daily basis contributes significantly to the energy and mixing of the water in the karst supporting the modeling hypothesis of well mixed groundwater reservoirs.

Hydrologic Processes
Fifteen years of hydro-meteorological data (2004-2019) were used to simulate the hydrology of Koiliaris CZO using the Karst-SWAT model. The goodness of fit of the Karst-SWAT model calibration is presented in Supplementary Figure S5. The statistical indices Nash-Sutcliffe efficiency (NSE), percent bias (PBIAS) and RMSE-observations standard deviation ratio (RSR), proposed by Moriasi et al. [83] were satisfactory on a daily time step (NSE equal to 0.67, PBIAS equal to 0.2% and RSR equal to 0.57) and the coefficient of determination between daily observed and modeled flows was 0.86 and the slope 1.07. Figure 5 depicts the mean annual hydrologic budget of the basin, including the extended karstic region. Based on the hydrologic simulation, the precipitation in the area (1428.6 mm) is partially evapotranspirated (480 mm, 33.6% of rainfall) and partially becomes surface runoff (177 mm, 12.4% of rainfall). The surface runoff is further subdivided into runoff from the Keramianos tributary (51.6 mm) and runoff from the surface streams of Anavreti (122.7 mm). Small amounts of water are added along the stream's way until the exit (2.7 mm). The flow of all springs is estimated at 541.5 mm (37.9% of rainfall). A substantial percentage of rainfall (15.3%) becomes surface runoff in the extended karst area (218.2 mm). Fifteen years of hydro-meteorological data (2004-2019) were used to simulate the hydrology of Koiliaris CZO using the Karst-SWAT model. The goodness of fit of the Karst-SWAT model calibration is presented in Supplementary Figure S5. The statistical indices Nash-Sutcliffe efficiency (NSE), percent bias (PBIAS) and RMSE-observations standard deviation ratio (RSR), proposed by Moriasi et al. [83] were satisfactory on a daily time step (NSE equal to 0.67, PBIAS equal to 0.2% and RSR equal to 0.57) and the coefficient of determination between daily observed and modeled flows was 0.86 and the slope 1.07. Figure 5 depicts the mean annual hydrologic budget of the basin, including the extended karstic region. Based on the hydrologic simulation, the precipitation in the area (1428.6 mm) is partially evapotranspirated (480 mm, 33.6% of rainfall) and partially becomes surface runoff (177 mm, 12.4% of rainfall). The surface runoff is further subdivided into runoff from the Keramianos tributary (51.6 mm) and runoff from the surface streams of Anavreti (122.7 mm). Small amounts of water are added along the stream's way until the exit (2.7 mm). The flow of all springs is estimated at 541.5 mm (37.9% of rainfall). A substantial percentage of rainfall (15.3%) becomes surface runoff in the extended karst area (218.2 mm). In terms of mean annual flow, the contribution of the Anavreti flow is 44.9% to the total spring flow (55.1% for the Stylos springs) and 33.9% to the total flow of the basin (41.5% for the Stylos springs). The mean annual contribution of Keramianos to the total flow of the basin is 7.2% and for the case of the surface streams at Anavreti, their contribution is 16.9%. The contribution of the hydrologic pathways during a dry and a wet year changes are presented in Supplementary Table S3. More specifically, the stream of Keramianos is more sensitive to the meteorological conditions of the year of interest and its contribution to total flow varies from 5.2% (dry year) to 10.6% (wet year). The intermittent stream of Keramianos drains a sub-basin of 33.1 km 2 (25.1% of the main Koiliaris CZO) before entering a karstic gorge. The Keramianos stream exits the gorge only when daily precipitation at the meteorological station M1 (Figure 1) exceeds 120 mm [17]. Therefore, on an annual basis, about 57% of the Keramianos water amount is lost in the karstic gorge. Therefore, during wet years, the contribution of Keramianos to the total flow of the basin almost doubles. In terms of mean annual flow, the contribution of the Anavreti flow is 44.9% to the total spring flow (55.1% for the Stylos springs) and 33.9% to the total flow of the basin (41.5% for the Stylos springs). The mean annual contribution of Keramianos to the total flow of the basin is 7.2% and for the case of the surface streams at Anavreti, their contribution is 16.9%. The contribution of the hydrologic pathways during a dry and a wet year changes are presented in Supplementary Table S3. More specifically, the stream of Keramianos is more sensitive to the meteorological conditions of the year of interest and its contribution to total flow varies from 5.2% (dry year) to 10.6% (wet year). The intermittent stream of Keramianos drains a sub-basin of 33.1 km 2 (25.1% of the main Koiliaris CZO) before entering a karstic gorge. The Keramianos stream exits the gorge only when daily precipitation at the meteorological station M1 (Figure 1) exceeds 120 mm [17]. Therefore, on an annual basis, about 57% of the Keramianos water amount is lost in the karstic gorge. Therefore, during wet years, the contribution of Keramianos to the total flow of the basin almost doubles.
These hydrologic pathways are of great interest during extreme events, as was the case of the recent flood of February 25th 2019 (Supplementary Figure S6). Supplementary Table S4 presents the modeled flows for every stream and spring in the basin during the peak daily flow of the extreme event. As it is obvious from the Supplementary Table S4, the floods are mainly caused by the surface runoff, and the surface streams of Anavreti and Keramianos have the highest contribution to total flow (37.2% and 26.5% respectively), followed by the spring of Anavreti, which has a slightly bigger contribution than the Stylos springs (19.0% and 17.3%).
The hydrologic budget during the extreme hydrologic event is calculated during a time span of 5 days precipitation commenced on 24 February, it produced the highest daily rainfall amount on the 25th and the following three days were also included so that surface and spring flow could manifest to a satisfactory degree. The total precipitation was 400 mm during these five days. The surface runoff was estimated at 97.9 mm and the surface runoff of the extended karst area at 179.6 mm. More than 5 million m 3 of water were exported from the Koiliaris watershed in one day generating a significant flood in the area [42]. This means that 70% of the rainfall becomes surface runoff. 10% of rainfall becomes spring flow during these five days (38.4 mm), increasing the flood magnitude. Finally, only 2.8 mm are evapotranspirated during the five days of the extreme event (0.7%), in comparison to 33.6% (evapotranspiration to rainfall) on a mean annual level.
The hydrologic analysis provided insights on the hydrologic behavior of the karst during flood events and confirmed that, when the infiltration rate into the groundwater is exceeded by precipitation, the potential for significant and catastrophic floods is certain. Table 4 presents the average chemical parameter values for the nine representative sampling locations (Agios Georgios (R1), Stylos spring (S1), Armenoi spring (S3), Anavreti spring (S5), Zourpos spring (S4), Katochori spring (S2), Kampoi well (G3), Dug well (G2), Katochori temporary river (R2)) using data from September 2010 to August 2019.

Geochemical Characterization and Weathering Fluxes
For the categorization of the geochemistry of Koiliaris CZO, a PCA was conducted using the average values of different chemical species (Ca, Mg, K, Si, Ti, B, Ba, Sr, Na, Li, pH, EC, DO, N-NO 3 , N-NH 4 , P-PO 4 , Cl, DN, SO 4 ) for the years between September 2010 and August 2019. The results of the PCA are presented in Figure 6. The score plot of the first two components of the PCA and the samples can be categorized into three distinct groups. The first three components describe 95% of the chemical variability. Stations R1, S1, S3, and S5 belong to Group 1 which has high pH and dissolved oxygen (DO), low nutrient, chloride and anion and cation concentrations. PCA component 1 distinguishes the karstic water from other origin or impacted water. Stations S2 and R2 in alluvial deposits, and G2 belong to Group 2 which is characterized by high values for nutrients like (N-NO Table 4). The water samples from the shallow aquifers have higher load of chemical species and this is related to higher load of nutrients from agricultural activities and the higher reactivity of brittle thin limestone found in alluvial formations. They also present higher values of Cl on average that range between 45 and 65 ppm.  Finally, G3 is situated in alluvial deposits originated from the weathered products of metamorphic schists. S4 has been affected from sea water and belongs to Group 3 which has high concentrations of Mg 2+ , Na + , Cl − and other constituents. G3 and S4 show high concentrations of Mg (42 ppm and 23 ppm, respectively), Na (41 ppm and 177 ppm respectively), Cl (78 and 275 ppm respectively) and SO4 (124 and 47 ppm respectively) ( Table 4).
The biplot analysis that identifies the chemical constituent that separates the three groups is presented in Supplementary Figure S7. In addition, a PCA was conducted using all data (and not only the average values) for the year 2012-2013 (Supplementary Figure S7) and the results in terms of characterization were the same as in Figure 6 discussed above.
In order to assess the geochemical and anthropogenic impact on water quality, the time series of the data were normalized to conservative values and compared with scientific literature. For instance, the molar (Ca + Mg)/Na ratio in Koiliaris River (R1) is 5 ± 1 and Ca/Na molar ratio is 4 ± 1, suggesting calcite dominated waters. When the molar ratio of (Ca + Mg)/Na in the rivers is greater than 10, the water is influenced by dolostone dissolution [84]. In addition, the molar (Ca + Mg)/SO4 ratio was calculated 44 ± 44 suggesting that there is no influence of sulphide oxidation or gypsum dissolution.
A positive correlation between Na + and Cl − concentrations in the water bodies of Koiliaris CZO was observed, with Cl/Na ratios to be higher than the ratio of 1.04 resulting from the atmospheric data of Finokalia station (FKL) (https://deims.org/d16cc035-054d-4b0e-9e19-6e9ec10aecd0) situated on the north coast of Crete (Figure 7a). This suggests that there is an additional source of Cl, through enrichment that is probably derived from fertilizers or manure input. S4 is not depicted in the Figure 7a as it is out of range of the curve between Na + and Cl − , presenting high values of these species due to sea water intrusion. In Figure  7b, where molar Ca/Na ratio is plotted as a function of the molar NO3/Na ratio, Groups 1, 2, and 3 of the different water bodies also appeared. Al, V, Cr, Ni, As, Se, Cd, Hg, and Pb values for all sites were extremely low or below detection limit of the analytical methods and they are not reported in Table 4. Finally, G3 is situated in alluvial deposits originated from the weathered products of metamorphic schists. S4 has been affected from sea water and belongs to Group 3 which has high concentrations of Mg 2+ , Na + , Cl − and other constituents. G3 and S4 show high concentrations of Mg (42 ppm and 23 ppm, respectively), Na (41 ppm and 177 ppm respectively), Cl (78 and 275 ppm respectively) and SO 4 (124 and 47 ppm respectively) ( Table 4).
The biplot analysis that identifies the chemical constituent that separates the three groups is presented in Supplementary Figure S7. In addition, a PCA was conducted using all data (and not only the average values) for the year 2012-2013 (Supplementary Figure S7) and the results in terms of characterization were the same as in Figure 6 discussed above.
In order to assess the geochemical and anthropogenic impact on water quality, the time series of the data were normalized to conservative values and compared with scientific literature. For instance, the molar (Ca + Mg)/Na ratio in Koiliaris River (R1) is 5 ± 1 and Ca/Na molar ratio is 4 ± 1, suggesting calcite dominated waters. When the molar ratio of (Ca + Mg)/Na in the rivers is greater than 10, the water is influenced by dolostone dissolution [84]. In addition, the molar (Ca + Mg)/SO 4 ratio was calculated 44 ± 44 suggesting that there is no influence of sulphide oxidation or gypsum dissolution.
A positive correlation between Na + and Cl − concentrations in the water bodies of Koiliaris CZO was observed, with Cl/Na ratios to be higher than the ratio of 1.04 resulting from the atmospheric data of Finokalia station (FKL) (https://deims.org/d16cc035-054d-4b0e-9e19-6e9ec10aecd0) situated on the north coast of Crete (Figure 7a). This suggests that there is an additional source of Cl, through enrichment that is probably derived from fertilizers or manure input. S4 is not depicted in the Figure 7a as it is out of range of the curve between Na + and Cl − , presenting high values of these species due to sea water intrusion. In Figure 7b, where molar Ca/Na ratio is plotted as a function of the molar NO 3 /Na ratio, Groups 1, 2, and 3 of the different water bodies also appeared. Al, V, Cr, Ni, As, Se, Cd, Hg, and Pb values for all sites were extremely low or below detection limit of the analytical methods and they are not reported in Table 4. Carbonate weathering rates were calculated using observed and modeled monthly flows. Specifically, observed monthly flows were used for Keramianos ephemeral river and modeled monthly flows were used for S1 and S3 springs. The fluxes were normalized to the contributing surface area of the karst for each spring and the surface area of the sub-catchment of Keramianos. The extended karst area was included in the calculations of the S1 spring weathering rates. Annual carbonate weathering rates range from 10 to 130 t/km 2 /year and 50 to 170 t/km 2 /year for Keramianos ephemeral river and the springs of Koiliaris CZO respectively (Figure 8), which corresponds to about 25 mm/1000 years and 39 mm/1000 years of carbonate being removed by chemical weathering for Keramianos ephemeral river and the springs, respectively (assuming a calcite rock density of 2.5 t/m 3 ). These values are lower to similar estimates reported in the Carbonate weathering rates were calculated using observed and modeled monthly flows. Specifically, observed monthly flows were used for Keramianos ephemeral river and modeled monthly flows were used for S1 and S3 springs. The fluxes were normalized to the contributing surface area of the karst for each spring and the surface area of the sub-catchment of Keramianos. The extended karst area was included in the calculations of the S1 spring weathering rates. Annual carbonate weathering rates range from 10 to 130 t/km 2 /year and 50 to 170 t/km 2 /year for Keramianos ephemeral river and the springs of Koiliaris CZO respectively (Figure 8), which corresponds to about 25 mm/1000 years and 39 mm/1000 years of carbonate being removed by chemical weathering for Keramianos ephemeral river and the springs, respectively (assuming a calcite rock density of 2.5 t/m 3 ). These values are lower to similar estimates reported in the literature. For instance, Jura Mountains carbonate chemical weathering was estimated to be 100 mm/1000 years [74]. We postulate that the lower weathering rates are due to reduced impact from vegetation, the periodicity of precipitation, shallow soil development and high karstification rates. However, Jura Mountains are heavily forested (over 54% of the watershed area) all along its altitudinal range while Koiliaris is sparsely vegetated at the mid elevation and has alpine conditions at the high elevations. The periodicity of precipitation is also different between Jura Mountains (rains thought out the year) and Koiliaris (has a 5-6 months dry period). Figure 8 also presents the monthly based weathering rates for Keramianos and the springs. The slope of the curves is maintained between annual and monthly values and the degree of linearity (R2) is very strong for both and of similar values. Finally, we can observe that there is no interannual variation. Both annual and monthly weathering rates were also calculated for two other springs (Agia and Meskla springs) of the wider area of Koiliaris CZO and the results were consistent with the rates calculated in Keramianos and the springs S1 and S3.
literature. For instance, Jura Mountains carbonate chemical weathering was estimated to be 100 mm/1000 years [74]. We postulate that the lower weathering rates are due to reduced impact from vegetation, the periodicity of precipitation, shallow soil development and high karstification rates. However, Jura Mountains are heavily forested (over 54% of the watershed area) all along its altitudinal range while Koiliaris is sparsely vegetated at the mid elevation and has alpine conditions at the high elevations. The periodicity of precipitation is also different between Jura Mountains (rains thought out the year) and Koiliaris (has a 5-6 months dry period). Figure 8 also presents the monthly based weathering rates for Keramianos and the springs. The slope of the curves is maintained between annual and monthly values and the degree of linearity (R2) is very strong for both and of similar values. Finally, we can observe that there is no interannual variation. Both annual and monthly weathering rates were also calculated for two other springs (Agia and Meskla springs) of the wider area of Koiliaris CZO and the results were consistent with the rates calculated in Keramianos and the springs S1 and S3.
Supplementary Figure S8 presents the carbonate weathering intensity as a function of runoff in Keramianos ephemeral river and in S1 and S3 springs. With increasing discharge, concentrations do not follow a pure dilution trend and remain constant over the whole range of discharge (Supplementary Figure  S8). The ratio between Ca/Mg was calculated to be 4 ± 2 suggesting calcite dissolution (Ca/Mg ratio < 1.5 suggests influence of dolomite). Moreover, the air temperature does not affect the dissolution intensity and does not follow the general boomerang curve observed in the GLORICH dataset [84]. Supplementary Figure S8 presents the carbonate weathering intensity as a function of runoff in Keramianos ephemeral river and in S1 and S3 springs. With increasing discharge, concentrations do not follow a pure dilution trend and remain constant over the whole range of discharge (Supplementary Figure S8). The ratio between Ca/Mg was calculated to be 4 ± 2 suggesting calcite dissolution (Ca/Mg ratio < 1.5 suggests influence of dolomite). Moreover, the air temperature does not affect the dissolution intensity and does not follow the general boomerang curve observed in the GLORICH dataset [84].

Conclusions
The Koiliaris CZO is a karstic Mediterranean watershed on the island of Crete and is part of the European LTER Network and the LTER-Greece Network. Significant research has been conducted over the past 15 years in the areas of stream and groundwater chemistry, soil degradation and formation, hydrologic and geochemical modeling, sustainable management of water and soil resources coupled with high frequency data acquisition. In this study, we combine the knowledge generated in previous studies and examine geomorphological and hydrogeochemical monitoring data to provide evidence and test hypotheses to further elucidate the hydrologic and geochemical response of Koiliaris CZO.
The extended karst area of Koiliaris River was redrawn based on geomorphologic and tectonic analyses. It was extended in the southern border by including the caves and sinkholes along an anticline, and was reduced in the northern border by eliminating the Neogene grabens. The total area of the extended karst remained the same as it was previously estimated.
The harmonic analysis of the fluctuation of the spring flow in Koiliaris shows that it is related to Earth tide stressing of the karstic reservoirs, as it affects both the lower and upper reservoirs. The analysis confirmed the partitioning of the karst into two reservoirs that have different hydro-mechanical properties and behavior. The oscillating hydraulic gradient generated by the tidal loading contributes significantly to the energy and mixing of the water in the karst supporting the modeling hypothesis of well-mixed groundwater reservoirs.
The geochemistry of Koiliaris was categorized into three distinct groups: (a) group 1 identifies the karstic water from other origin or impacted water, (b) group 2 is characterized by high values for nutrients, and (c) group 3 has high concentrations of Mg 2+ , Na + , Cl − and other constituents. In general, the water of Koiliaris CZO is influenced by calcite dissolution. Carbonate weathering rates were found to be lower than similar studies, possibly due to a reduced impact from vegetation, the periodicity of precipitation, shallow soil development, and high karstification of the carbonate bedrock.
The conceptualization of a two reservoir, well-mixed karstic system was confirmed by both the geomorphologic and the tidal analysis. The hydrologic analysis provided insights on the hydrologic behavior of the karst during flood events and confirmed that when the infiltration capacity of the karst is exceeded by precipitation, the potential for significant and catastrophic floods is certain. During extreme events, 70% of the precipitation becomes surface runoff, creating tremendous flood events that need to be managed. Modeling such events is of the upmost importance, considering the impacts of climate change [40] and the need for adaptation measures to climate change and risk management [42]. Finally, the geochemical response of the watershed is closely linked to the geomorphology and hydrologic response, justifying in this way the weathering rates determined for the catchment. The studies at Koiliaris CZO illustrate the importance of critical zone science and they show how such transdisciplinary studies can be used to develop tools that can be scaled up for water and soil management at regional scales as well as for the development of strategies for the mitigation of climate change impacts.  Figure S7: Groups resulting from PCA in water bodies of Koiliaris CZO using the average values of chemical species, Figure S8: Relationship between weathering intensity and instantaneous water discharge, Table S1: Tidal signals comparison for period 2019 ER, LR. Tidal constituents' amplitude and phase with 95% CI estimates; phases are at central time, Table S2: Tidal signals comparison for period 2018 ER, LR. Tidal constituents' amplitude and phase with 95% CI estimates; phases are at central time, Table S3: Mean contribution (%) of each stream and spring to the total flow (at the basin exit) and contribution during a wet and dry year, Table S4

Conflicts of Interest:
The authors declare no conflict of interest.