We developed a high-resolution MIKE SHE/MIKE 11 model of a 231.3 ha headwater catchment in the granitic uplands of the French Massif Central to estimate the contribution of groundwater upwelling to the water balance of the Dauges mire, an acidic valley mire of international importance for nature conservation. We estimated that groundwater upwelling from the underlying weathered granite formations – mostly an approximately 55 m deep fissured zone – provides 27.1 % of total long-term inflows to the mire. This contribution increases to 37.2 % in September when total inflows are small. Overland boundary inflow accounts for an average of 40.2 % of total inflows. However, most of this originates from groundwater seepage through mineral soils along the mire margins or in small non-channelised valleys upslope of the mire. A sensitivity analysis showed that model performance in terms of the simulation of mire groundwater levels was most sensitive to parameters describing the mineral soils and weathered granite formations rather than the overlying peat layer. Variation partitioning demonstrated that groundwater upwelling was the most important factor driving simulated monthly groundwater table depth within the mire. Sustained groundwater upwelling maintains the mire water table close to or at ground level for most of the year. As a result, precipitation and overland boundary inflows quickly leave the wetland as saturation-excess runoff. There was close agreement between the observed distribution of mire habitats and areas where the simulated long-term groundwater seepage rate was larger than zero in September. Our results demonstrate that, contrary to the assumed small contribution of groundwater to the hydrology of hard-rock regions, groundwater upwelling from underlying weathered formations can be a quantitatively important and functionally critical element of the water balance of valley mires in granitic headwater catchments. These results have important legal and management implications.
Wetlands, and in particular mires (i.e. wetlands actively accumulating peat, Rydin and Jeglum, 2006), are widely recognised as providing multiple ecosystem services (Okruszko et al., 2011). For example, mires are more efficient than any other terrestrial ecosystems at sequestering carbon. Despite occupying only 3 % of the world's land area, they contain twice as much carbon as all forests (Frolking et al., 2011; Lindsay, 2010; Parish et al., 2008; Worrall et al., 2011; Yu et al., 2011). Mires are characterised by very distinctive environmental conditions and thereby support unique ecosystems and many specialised species (Parish et al., 2008). The long-term provision of the ecosystem services provided by mires depends on their peat remaining waterlogged and therefore on the long-term stability of the wetlands' water balance. Accurate characterisation and evaluation of the water balance of a mire, and indeed any wetland, is a vital prerequisite for understanding ecosystem functioning, developing management strategies and assessing potential impacts of future perturbations that may include climate change or other anthropogenic activities within the wetland or its catchment (e.g. Al-Khudhairy et al., 1999; Baker et al., 2009; House et al., 2016b; Thompson et al., 2017).
Groundwater inflow has long been recognised as having a critical role in the water balance of many mires in sedimentary contexts (Boeye and Verheyen, 1992; Gilvear et al., 1993; House et al., 2016a; Koerselman, 1989; Rossi et al., 2012; Siegel and Glaser, 1987; Wassen et al., 1990). In contrast, the importance of groundwater for the water balance of wetlands in hard-rock regions, especially in upland and mountainous areas, has most often been considered, if not negligible, then at least less significant. This is based on the assumption, common to both hydrogeologists and hillslope hydrologists, that beneath the soil and saprolite (i.e. the loose in situ products of weathering), the bedrock is mostly impermeable (Banks et al., 2009; Gabrielli et al., 2012; Haria and Shand, 2004; Lachassagne, 2008; Tromp-van Meerveld et al., 2007). However, in the last 2 decades, hydrogeologists have progressively realised that classical conceptual models of groundwater flow in hard-rock regions, in particular within granitic areas, may need to be revised. Granite weathering over sufficiently long periods and under suitable environmental conditions can produce a stratiform lateritic profile, which, in addition to the saprolite, includes a densely fissured layer typically 40–50 m deep (Dewandel et al., 2006; Lachassagne, 2008). This layer is characterised by numerous sub-horizontal fractures of weathering origin (Lachassagne et al., 2011; Maréchal et al., 2004), the density and hydraulic connectivity of which decrease with depth until the fresh bedrock (Guihéneuf et al., 2014). In parallel, hillslope hydrologists have questioned the “impermeable bedrock dogma” (Tromp-van Meerveld et al., 2007) and re-evaluated the role of “deep” groundwater in streamflow generation in granitic areas (Gabrielli et al., 2012; Haria and Shand, 2004; Kosugi et al., 2006, 2011; Uchida et al., 2003).
The contribution of “deep” groundwater to the water balance of many headwater mires in hard-rock regions may, therefore, have been overlooked, in particular where these mires are located in topographic lows. A small number of observational studies have suggested that such mires correspond to areas of sustained groundwater upwelling and seepage (Branfireun and Roulet, 1998; Morley et al., 2011; Šanda et al., 2014). However, direct long-term quantification of groundwater upwelling at a high spatial resolution through field measurements is currently impractical, and as a result these studies provide limited insights into the degree of groundwater dependence of these systems. Physically based, integrated, spatially distributed modelling approaches provide an alternative way to constrain estimates based on limited data and are often associated with smaller errors than estimates based solely on field observations (Gilvear and Bradley, 2009). This approach also provides a means of investigating potential impacts on future ecohydrological conditions using alternatively specified scenarios (e.g. Clilverd et al., 2016; House et al., 2016b; Thompson et al., 2009, 2017; Thompson, 2012). However, the use of physically based spatially distributed models at high resolution in hard-rock regions is complicated by the large contrast in hydrodynamic properties between the largely impermeable matrix and the water-bearing fractures and by the small-scale variability in the density and connectivity of fractures in the fissured zone (Levison et al., 2014; Singhal and Gupta, 2010). Integrated hydrological models able to represent groundwater flow in fractured media are still scarce (Brunner and Simmons, 2012), and hydrodynamic properties of the fissured zone are extremely difficult to characterise at high resolution (Singhal and Gupta, 2010). As a result the vast majority of modelling studies have used the equivalent porous medium approach, assuming that the fractured medium behaves as a continuous porous medium with more or less homogeneous hydrodynamic properties at larger spatial scales (Long et al., 1982). There is a relative consensus that this approach provides acceptable results at the regional scale but poorly reproduces local flow systems (Singhal and Gupta, 2010). The applicability of the equivalent porous medium approach to provide high-resolution estimates of groundwater inflow to small headwater wetlands in hard-rock regions is therefore unclear.
Here we apply the MIKE SHE/MIKE 11 modelling system to an acidic valley mire in the granitic uplands of the French Massif Central, as a case study that is representative of valley mires common in the granitic uplands of the Variscan Belt (Cubizolle, 2019; Etlicher, 2005; Tanneberger et al., 2017). Our objectives are (i) to test the ability of an equivalent porous medium approach with limited data on the hydrodynamic properties of the weathered granite formations to reproduce high-resolution spatial and temporal patterns in groundwater seepage and groundwater table depth within the mire, (ii) to quantify the mire water balance including its dependence on groundwater inflows from weathered granite formations, and (iii) to investigate the hydrological processes driving groundwater table depth in the mire.
The 231.3 ha catchment that is modelled broadly matches the boundaries of
the Dauges National Nature Reserve (NNR), located near
Saint-Léger-la-Montagne, in the administrative department of
Haute-Vienne, Nouvelle-Aquitaine, France (latitude: 46
The Dauges mire, its catchment and the hydrological monitoring network.
Groundwater level was monitored from December 2010 to October 2013 in a network of 16 shallow dipwells installed within the mire and on the lower mineral slopes (Fig. 1). These were equipped with automatic loggers (Mini-Diver, Schlumberger) during all or part of the monitoring period. The 15 min records were corrected for atmospheric pressure recorded by a barometric logger (Baro-Diver, Schlumberger) installed at the site and for logger drift based on regular manual checks and aggregated to daily mean values (Duranel, 2015).
Discharge was measured from January 2011 to December 2013 at the outlet of the main mire area (Pont de Pierre, Fig. 1) using an Orpheus-Mini (OTT) level logger. The 15 min stage records were corrected for logger drift based on regular manual checks and converted to discharge on the basis of a rating curve established from spot discharge measurements undertaken using dilution gauging or an electromagnetic flowmeter and the velocity–area method (Duranel, 2015). Discharge was also measured using sharp-crested V-notch weirs and float loggers (Thalimèdes, OTT) from January 2011 to June 2013 at three locations in the upper reaches of the mire (Girolles, Marzet and Rocher, Fig. 1). Stage records were corrected for logger drift, converted to discharge using the equation provided by Dingman (1994, p. 544) and aggregated to daily mean flow.
Radiation, air temperature, relative humidity, wind speed and rainfall were measured 2 m above ground at the centre of the mire using an Enerco 404 (Cimex) meteorological station from June 2010 to March 2013. Generalised least-square regression models were used to reconstruct long-term (1 August 1998–31 December 2013) daily meteorological time series based on records from nearby Météo-France permanent weather stations, including a station located 4.2 km from the centre of the research site. The reconstructed dataset was used to extend the observed time series and infill missing records (Duranel, 2015). Reference evapotranspiration (ETo) was calculated using the FAO Penman–Monteith method (Allen et al., 1998).
A 3D geological model was built based on a range of field investigations including topographic surveying, geological drilling, electrical resistivity tomography, manual augering and probing, and analysis of existing outcrops. This model is detailed in Duranel (2015) and is summarised here.
Three different datasets were aggregated to produce a digital elevation
model (DEM). Within the mire, surface elevation was measured using a
differential geo-positioning system (DGPS) at approximately 5 m resolution
and with a vertical accuracy of around 10 cm. On mineral soils in the southernmost part of the catchment, topographic data were extracted from the IGN BD Topo, a nationwide DEM with a resolution of 25 m and a root mean square error (RMSE) of 4.8 m against DGPS ground-truthing points within the site.
Finally, on mineral soils in the rest of the catchment, the DEM was based on
mass points and contour lines at 1 m intervals digitised from existing
Electrical resistivity tomography (ERT) was used to estimate the depth of weathered granite formations along three 380 to 868 m long transects across the mire and the lower hillslopes and one 315 m long transect located on the southern hilltop. A detailed description of these profiles is given in Duranel (2015). Figure 2 shows the two most informative profiles, the locations of which are shown in Fig. 1. A Schlumberger configuration was used in both cases, with an electrode spacing of 5 m and a maximum penetration depth of 60 m below ground level. The inversion was performed using the Res2Dinv software (Loke, 2013; Loke and Barker, 1996) and a L1 norm optimisation, leading to RMSEs of 5.3 % for Profile A and 2.2 % for Profile B. Both profiles are characterised by strong contrasts in resistivity values. The most striking feature is the presence on each side of the profiles of highly resistive material overlaying more conductive material. Highly resistive material is only recorded beneath mineral soils on hilltops and slopes and does not occur beneath the mire. The transition between these layers is very sharp, its altitude is largest beneath hilltops, and its depth greatest beneath steep slopes. Out of the four ERT profiles, the only location where conductive material was recorded above resistive material is on one hilltop (visible on the left-hand side of Profile A). These observations were interpreted as demonstrating that most of the material investigated corresponds to the fissured layer of a truncated weathered granite profile. Highly resistive and conductive materials correspond to the unsaturated and saturated fissured layers, respectively. Taking into account the ERT positional accuracy, the margins of the mire coincided on all transects with the location where the inferred groundwater table reaches the ground surface, suggesting a determining role of groundwater seepage in the mire water balance and a high degree of hydrological connectivity between the mire and the fissured layer. The increase in resistivity with depth at the very bottom of Profile B at a depth of around 55 m below ground level was interpreted as a decrease in fissure density (and therefore bulk porosity) and a transition towards unweathered granite. The shallow (less than 2 m deep) superficial layer of conductive material on the hilltop on the left-hand side of Profile A is consistent with (probably unsaturated) saprolite, which is more conductive than the fissured granite layer due to its larger clay and water contents (Baltassat et al., 2005). This configuration was not seen anywhere else within the study area, leading to the conclusion that on most hilltops and hillslopes the majority of saprolite has been eroded away and the combined thickness of the soil, periglacial deposits and remaining in situ saprolite is too small (less than a metre) to have been detected during the ERT survey. There is no indication of the presence of a substantial saprolite layer in the valley bottom beneath the mire; however, the complete saturation of the profile may make the distinction between fissured granite and saprolite impossible in this area.
Electrical resistivity tomography images obtained after the inversion of apparent resistivity data acquired using a Schlumberger configuration. The survey was conducted using 80 electrodes for Profile A and 144 electrodes for Profile B, with an electrode spacing of 5 m. The inversion was performed using the Res2Dinv software and a L1 norm optimisation, leading to RMSEs of 5.3 % for Profile A and 2.2 % for Profile B. White dots show the centres of model blocks, the black line the surface topography and the black vertical marks the electrode positions. The profile locations are shown in Fig. 1.
Large-scale underground uranium mining was undertaken in the area during the
second part of the 20th century. Logs of seven geological boreholes
drilled in 1973 by the Commissariat à l'Energie Atomique for the purpose
of uranium exploration were provided by the mining company Areva (now
Orano). These were all drilled along the eastern margin of the mire, across
presumed mineralised faults located using ground surface radiological
surveys. They were drilled at an angle of 49.5
Geological boreholes drilled by the Commissariat à l'Energie Atomique, showing the granite weathering grade and the core recovery percentage. The borehole locations are shown in Fig. 1. The “highly weathered” and “weathered” layers were interpreted as corresponding to the fissured layer, while the “slightly weathered” to “unweathered” layers were interpreted as corresponding to the bedrock.
Our interpretation of the ERT and geological data agrees with the current understanding of granite weathering processes, granite landscape geomorphology and granite hydrogeology in the Massif Central in general and in the Monts d'Ambazac specifically (Desire-Marchand and Klein, 1986; Dewandel et al., 2006; Godard et al., 2001; Klein, 1978; Klein et al., 1990; Mauroux et al., 2009). The absence of substantial saprolite formations from hilltops and slopes, and presumably from the downstream part of the valley, is in line with Mauroux et al. (2009), who used airborne spectral radiometry to demonstrate that saprolite has been largely eroded away in the Monts d'Ambazac.
Existing sections cutting through periglacial deposits and weathered formations and resulting from small-scale grus extraction or from road construction were also described (Duranel, 2015). These investigations
showed that periglacial deposits are very patchy, thin where they occur, and
can be considered to be hydrodynamically equivalent to in situ saprolite. Pedological pits dug by Verger (1998) and Gratia (2012 in Duranel, 2015)
showed that soils outside the mire are generally relatively deep
(40–70 cm), even on hilltops and relatively steep slopes such as those
found on Puy Rond. According to the French pedological classification (Baize
and Girard, 2009), they mostly belong to acidic
Peat deposits were mapped using manual augering and probing. The average
peat depth within the mire was 0.80 m (standard deviation 0.49 m), with a
maximum of 3.45 m at its centre (Duranel, 2015). Stratigraphical surveys and
slug tests showed a rapid increase in peat humification and a decrease in
hydraulic conductivity with depth, which is consistent with observations
from many mires worldwide. Ingram (1978), based on concepts developed by
Ivanov (1953) and Romanov (1968), proposed the definition of “acrotelm” as the upper layer characterised by a low humification and high hydraulic
conductivity and porosity, and “catotelm” the lower layer with the opposite
characteristics. At the Dauges site, peat hydraulic conductivities measured
between 0.1 and 0.2 m below ground level had a median value of
The MIKE SHE/MIKE 11 modelling system (Graham and Butts, 2005; Refsgaard et al., 2010) was used to model the mire and its catchment. MIKE SHE/MIKE
11 is a deterministic, fully distributed modelling system able to model
actual evapotranspiration, unsaturated zone flow and storage, saturated zone
flow and levels, overland flow and channel flow processes in an integrated
manner. It has previously been successfully applied in a wide range of
wetland environments (e.g. Al-Khudhairy et al., 1999; Bourgault et al.,
2014; Clilverd et al., 2016; Hammersmark et al., 2010; House et al., 2016a;
Thompson, 2004; Thompson et al., 2004, 2017). Originally fully
physics-based, it has evolved to include a suite of modules that also allow
for the use of simplified physical equations or even conceptual models based
on distributed, semi-distributed or lumped approaches (Refsgaard et al.,
2010). In the present study, a combination of both physics-based and
conceptual modelling approaches was used. A 10 m
Saturated flow was modelled using an iterative implicit finite-difference technique to solve the 3D Darcy equation including within the fissured granite for which an equivalent porous medium approach (Long et al., 1982) was used. In the development of the saturated flow model, we used the most plausible geological conceptualisation based on the data detailed in Sect. 2.3 while making a number of reasonable simplifying assumptions to account for constraints associated with practical computing times, available data and the risk of over-parameterisation. These assumptions, further justified in Table 1, were that (i) saturated flow occurs predominantly in the fissured granite layer and in the peat; (ii) mineral soils, saprolite, periglacial deposits and alluvium can be neglected as far as saturated flow is concerned; (iii) the hydrologically active fissured granite layer is 55 m deep, follows the surface topography, and has homogeneous properties throughout; and (iv) peat deposits have homogeneous properties throughout.
Justification of simplifying assumptions of the saturated flow model component.
Only two geological formations (peat and fissured granite) were therefore represented in the saturated flow model, with peat only occurring within the mire extent. The saturated flow model was nevertheless divided into two computational layers throughout the model extent, because MIKE SHE does not allow for discontinuous computational layers. Throughout the model, the lower layer represented the fissured granite, with a constant depth of 55 m below the topographic surface. Outside of the mire, the upper computational layer had identical properties to the lower layer and also represented the (upper) fissured granite. In MIKE SHE the depth of the upper computational layer of the saturated flow component cannot be smaller than the depth of the root zone in the unsaturated flow component. Therefore, outside of the mire, the upper layer was defined as 2 m deep to match the estimated root depth of woodlands. Within the mire, where shallow-rooted open vegetation dominates, the upper layer represented peat deposits. Its depth followed that measured in the field, with a minimum value of 0.5 m to avoid numerical instabilities. A no-flow boundary was assumed around the boundaries of the topographic catchment. At the downstream end of the catchment, a constant-gradient boundary was used, assuming that the slope of the phreatic surface broadly follows that of the topographic surface. Overland flow was modelled using a 2D, finite-difference, diffusive wave approximation of the Saint-Venant equations, solved using the successive over-relaxation method.
The empirical two-layer unsaturated zone model (Yan and Smith, 1994) was
selected. This was principally to reduce computation time but was also justified by the lack of detailed data on unsaturated zone properties since
the lower number of parameters required by the two-layer model facilitates
their estimation through calibration. Due to the two-layer model's
limitations in representing interception in woodland, evaporation from
interception was modelled outside MIKE SHE using the same equations as those
implemented in the HYLUC model (Calder, 2003), with a distinction being made
for tall vegetation (woodland, heath and shrubs) between wet-time
evaporation (corresponding to evaporation from interception, modelled
outside MIKE SHE) and dry-time evapotranspiration (which includes
transpiration and any other evaporation from the ground surface, simulated
using MIKE SHE's two-layer model). Calder's
Parameterisation of the evapotranspiration model.
The MIKE SHE model was coupled to a MIKE 11 1D river model of the main channels within the catchment. In order to achieve logistically feasible run times during calibration and validation, channel flow was modelled using a kinematic routing method and water levels estimated a posteriori based on Manning's equation. This method is fast and generally relatively accurate for fast-flowing streams where no backwater occurs (DHI, 2009b). It is not expected to give an accurate representation of stream stage, and as a consequence it was assumed that exchanges between the overland flow component and watercourses occur towards the latter only (i.e. there is no over-bank flooding). This simplification is acceptable since within the research site over-bank flooding only occurs within a very narrow band along the main stream and in a limited area in the downstream part of the mire.
The Dauges MIKE SHE/MIKE 11 model used an adaptive time step, whilst results were aggregated to a daily resolution. Model outputs were extracted for the period 1 January 2001–31 December 2013, after an approximately 2.5-year warm-up period beginning on 1 August 1998 (corresponding to the installation of the nearest permanent automatic meteorological station).
The model was manually calibrated and validated against observed discharges and groundwater table depths. To accommodate slightly different data availabilities, the calibration period was set as 1 January 2011–30 June 2012 for discharge and 1 January 2011–31 December 2012 for groundwater table depth. The validation period was set as 1 July 2012–31 December 2013 for discharge and 1 January–31 December 2013 for groundwater table depth. Model performance was evaluated using the Nash–Sutcliffe efficiency (NSE) and percent bias (PBIAS) statistics recommended by Moriasi et al. (2007) and a visual assessment of the ability of the model to reproduce seasonal groundwater table patterns. Calibration parameters together with their calibrated values are shown in Table 3.
Calibration parameters of the MIKE SHE/MIKE 11 model of the Dauges catchment, including their final values and the range used in sensitivity analyses.
The ability of the calibrated model to simulate realistic groundwater tables across the catchment was further investigated by validating a raster map of the mean groundwater table depth simulated across the site for the 2001–2013 period against the observed mire boundaries, mapped using a combination of botanical (Durepaire and Guerbaa, 2008) and pedological (Duranel, 2015) criteria. The groundwater table depth threshold best discriminating between mire and other habitats was found by optimising a Cohen's kappa agreement function (Congalton, 1991) between simulated and observed mire boundaries.
The local sensitivity of the model around the calibrated parameter values
was evaluated by assessing the rate of change in the root mean square error
for simulated discharge and groundwater table depth when a single parameter
was perturbed by a defined proportion (Rochester, 2010). The sensitivity of
the model to a parameter
A water balance was computed using the MIKE SHE water balance tool for two areas upstream of the main mire outlet corresponding to mineral and peat soils, respectively. These water balance domains cover 36.9 and 125.0 ha, respectively (Fig. 4). The water balance was not calculated over the entire model domain for two reasons: firstly to reduce potential errors caused by uncertainty in boundary conditions at the north-western end of the model domain as well as the relatively sparse hydrogeological and hydrological data downstream of Pont de Pierre; and secondly to constrain the focus of the analysis on the main designated mire area upstream of Pont de Pierre. The stream flowing through the mire was not included in the water balance of the mire (i.e. only fluxes between the stream and peat soils or the surface of the mire were considered, with flow in and out of the stream where it enters and exits the mire excluded).
Water balance domains used within the MIKE SHE/MIKE 11 model of the Dauges catchment.
Simulated groundwater table depth, upward saturated flow from the lower computational layer (hereafter referred to as groundwater upwelling) and exchange from the saturated zone to overland flow (hereafter referred to as groundwater seepage) were mapped across the entire catchment.
Processes driving changes in groundwater table depth in peat soils were investigated using variation partitioning (Legendre and Legendre, 1998), based on simulated spatially averaged monthly means to focus on timescales relevant to vegetation (Wheeler et al., 2004). Explanatory variables included precipitation, groundwater upwelling, overland boundary inflow, actual evapotranspiration and groundwater table depth in the preceding month. An exploratory analysis showed that the relations between groundwater table depth and explanatory variables were mostly non-linear; therefore, variation partitioning was undertaken in a generalised additive modelling (GAM) framework (Wood, 2016).
Figures 5 and S1 in the Supplement show
observed and simulated daily discharge at the four gauging stations during
the calibration and validation periods. According to the guidelines of
Moriasi et al. (2007), model performance with regard to discharge at the
mire outlet (Pont de Pierre) was good to very good (NSE close to 0.75 during
both calibration and validation and overall, percent bias lower than
Observed and simulated stream discharge and model performance
statistics for two selected locations within the Dauges catchment
(1 January 2011–31 December 2013). Note different
Model performance with regard to groundwater levels was variable
(Figs. 6 and S2). It
should be assessed relative to the site's topography: the altitude range
from the bottom of the mire to the top of the hills surrounding it and
delineating the model boundaries is 110 m. The performance was therefore
clearly good to very good for a large number of dipwells (such as D7, D13,
D15, D18, D20 and D21) for which RMSEs of the order of
Observed and simulated groundwater table depth and model
performance statistics for selected dipwells within the Dauges catchment
(1 January 2011–31 December 2013). Note different
Continued.
The model performance was less good, but still satisfactory, for some dipwells (e.g. D3, D9 and D10) which are close to the mire boundary or within the mineral soils just outside it. In D3 for instance, the much shallower groundwater table levels observed in 2012 and 2013 compared to 2011 were not very well reproduced. Conversely in D10, simulated levels were too shallow during the summers of 2012 and 2013. Dipwells D3, D9 and D10 are located where discrepancies between the actual ground elevation of the observation point and that interpolated at the corresponding DEM grid cell are the largest. This is a result of the relatively steep slopes, and it is at these locations that the approximations caused by the model discretisation are the largest. During calibration, it became evident that simulated water table depths at these dipwells were more sensitive to the parameterisation of the fissured zone than those located at the centre of the mire. While it was possible to improve the fit for dipwells on one side of the mire, this improvement was balanced by deterioration in performance for dipwells located on the other sides. This may be explained by the fact that the fissured zone was modelled as a homogeneous layer, which facilitates calibration and reduces the risk of over-parameterisation at the model scale but does not allow for fine-scale adjustment to local variations in hydrodynamic properties (see Sect. 4.1).
Relatively poor performance was achieved for dipwells located close to the stream in the lower part of the mire (D8, D16, D17). Dipwells D16 and D17 are inserted into highly permeable alluvial gravel beneath or within the peat along the stream, and as a consequence groundwater table in these dipwells is largely driven by water levels in the nearby stream. The poor performance of the model in these locations was therefore a direct consequence of the absence of alluvial deposits in the model and the likely poor simulation of water levels in watercourses resulting from the kinematic routing method (Duranel, 2015). Finally, performance was poor for Dipwell 22 as a consequence of the presence of a 30 cm deep superficial layer of oxidised peat in this part of the mire resulting from historical agricultural drainage. This layer was not represented in the model due to its small extent and depth, and, as a result, the deeper and more variable observed groundwater level was not adequately simulated (Fig. 6). However, both the oxidised peat and the gravel deposits occupy only a small proportion of the mire extent, and their impact on the overall model performance is therefore considered to be small.
The model predicts the observed extent of mire habitats very well. The long-term mean annual groundwater level threshold that best discriminates between mire and non-mire habitats was found to be 0.166 m below ground level. This gives a value for Cohen's kappa agreement function of 0.841, corresponding to 95.0 % of grid cells being correctly classified (Table 4, Fig. 7). The good match between observed and simulated mire boundaries was not only observed within the main extent of the mire, but also along the narrow valley downstream and in the small sub-basins upstream (locations labelled (b), (d), (g) and (h) in Fig. 7). These small basins are located 10–30 m above the main mire, suggesting that the model simulates groundwater levels well even in the upper part of the catchment where no groundwater level time series were available for calibration. Small, localised discrepancies between the observed and simulated mire boundaries may be explained by different factors. Some false positives (grid cells wrongly predicted to support mire vegetation) may in part be due to errors in the mire vegetation map itself, for instance along the narrow valley in the north-east of the catchment (location (a) in Fig. 7) where difficult access and tree cover may have hindered mapping efforts. Location (c) coincides with relatively large errors produced by the MIKE SHE overland flow component. DEM inaccuracies may explain discrepancies at locations (e) and (f) along the south-eastern margin of the mire, since this area corresponds to the boundary between the low-resolution, low-accuracy IGN BD Topo DEM on the south-eastern hillslopes and the high-resolution, high-accuracy DGPS DEM within the mire.
Confusion matrix between the observed mire distribution based on
botanical (Durepaire and Guerbaa, 2008) and pedological (Duranel, 2015)
criteria and the predicted mire distribution based on model results
(simulated mean annual groundwater level
Observed mire boundary based on botanical (Durepaire and Guerbaa,
2008) and pedological (Duranel, 2015) criteria and the predicted mire
distribution based on model results (simulated 2001–2013 mean groundwater
table depth higher than 0.166 m below ground level).
Nevertheless, despite the issues discussed above, the satisfactory to very good performance at a large number of monitoring locations for both discharge and groundwater table depth and the excellent fit between observed and simulated mire boundaries suggest that the model reproduces the dominant hydrological characteristics of the mire, including its water balance which is discussed below.
Figure 8 provides the mean absolute scaled sensitivity coefficient for each model parameter that was tested with respect to groundwater level, discharge at the mire outlet (Pont de Pierre in Fig. 1) and discharges of the small upstream reaches (Girolles, Marzet and Rocher in Fig. 1). The parameters to which the calibrated model is most sensitive with regard to groundwater level are first and foremost those related to the hydrological characteristics of mineral formations: the horizontal saturated hydraulic conductivity, the depth and specific yield of the fissured zone, and the available water capacity of the unsaturated mineral soils. This is the case whether dipwells are located on mineral soils or within the mire. The peat parameter to which the model is most sensitive is the vertical saturated hydraulic conductivity, which controls upwelling/downwelling fluxes between the peat layer and the underlying mineral formations.
Summary of model sensitivity analysis results (see Table 3 for parameter codes). Parameters are arranged in order of increasing mean sensitivity from bottom to top.
The horizontal saturated hydraulic conductivity of the fissured zone is also the parameter to which simulated discharges of upstream reaches are by far most sensitive. Other influential parameters determine unsaturated flow in mineral soils: specific yield, available water capacity and the bypass fraction. Discharge at the mire outlet is also strongly influenced by parameters describing fluxes in mineral formations: available water capacity and specific yield in mineral soils and horizontal saturated hydraulic conductivity and depth of the fissured zone. Stream Manning's n is also very influential. Parameters describing the peat layer are less important.
Figure 9 shows the simulated long-term (2001–2013)
mean annual water balance of the mire (corresponding to the peat soil domain
in Fig. 4). Detailed water balances of both the
peat and mineral soil domains are provided in Table S1.
The water balance computational error was 9.9 % of total annual inflows
and was almost entirely due to the overland flow component. Extensive
attempts were made to reduce this computational error by tightening the
convergence criteria; however, this resulted in relatively small gains but a substantial increase in computing time. This is considered to be a minor
issue because, as a result of the quasi-constant saturation, very little
infiltration is simulated within the mire. Consequently, the overland flow
component error mainly affects simulated overland outflow to the river,
which is not the main focus of this study and has little impact on other water balance terms. The water balance error within the mineral soils domain
was only 1.6 % of total inflows. This is, again, mostly caused by the
overland flow component. The three main sources of water to the mire are
precipitation (32.1 % of total inflows), groundwater upwelling from the
underlying mineral formations (27.1 %), and overland boundary inflow
(i.e. non-channelised runoff from the mire catchment). The latter accounts for 40.2 % of total inflows; however, a large proportion of this originates from edge-focused groundwater seepage on mineral soils just
upstream of the mire (see Sect. 3.4). Lateral
saturated boundary inflows from the mineral catchment and from the river are
negligible (0.5 % and
Simulated mean annual water balance of the Dauges mire
(mm yr
Figure 10 breaks down the mire's long-term water balance by calendar month. The largest total inflows to the mire occur in winter and early spring. They then gradually decline throughout the spring and summer, reaching a minimum in September. Total outflows, dominated by overland runoff to the stream and to a lesser extent by evapotranspiration, follow the same pattern. Changes in storage are relatively small, reflecting the very shallow groundwater levels and the rapid removal of any water surplus through saturation-excess runoff, particularly in late autumn, winter and early spring. Due to the site's position at the transition between a degraded oceanic climate and a mountainous climate (Joly et al., 2010), precipitation is, on average, relatively well distributed throughout the year, although there is relatively large inter-annual variability. Overland boundary inflow is highly variable seasonally. It reaches its maximum in winter and early spring (accounting for 49.6 % of total inflows in January and February) before gradually declining to a minimum in September (22.9 %). It also exhibits large inter-annual variability, reflecting the similar variability in precipitation. Upward saturated flow from the lower computational layer (i.e. groundwater upwelling) is relatively constant from December to May before gradually declining to reach a minimum in September. It is characterised by much smaller seasonal and inter-annual variabilities than precipitation and overland boundary inflow and, therefore, generally provides a larger proportion of total inflows in summer and early autumn (37.2 % in September, Fig. S3). As noted above, groundwater upwelling accounts for, on average, 92.4 % of total inflows to the mire saturated zone. Seasonally this figure peaks at 95.5 % in January and declines to 88.1 % in August. At this time receding groundwater levels allow for the formation of an unsaturated zone across a larger proportion of the mire and for precipitation and overland boundary inflow to infiltrate and percolate, respectively, instead of producing saturation-excess runoff.
Simulated mean monthly water balance of the Dauges mire (2001–2013). Figures follow the MIKE SHE water balance convention: inputs are negative, outputs positive, change in storage positive when storage increases, and the water balance error is the sum of all inputs and outputs. Whiskers show the standard deviation. Top plot: see titles of bottom plots for variable codes.
Figure 11 shows the spatial distribution of
simulated long-term (2001–2013) mean annual groundwater table depth and
groundwater seepage rates in January and September when long-term mean
seepage rates are highest and lowest, respectively (data for other calendar
months are shown in Fig. S4). Spatial and temporal
patterns in groundwater upwelling rates are very similar to those in seepage
rates and are therefore not shown. As discussed above, the mire boundaries
are closely associated with shallow groundwater table depths. Groundwater
upwelling and seepage rates are highest along or immediately upstream of the
mire boundary. This is a clear example of edge-focused discharge (Richardson
et al., 2001), whereby groundwater discharge at the ground surface is
enhanced where there is a break in the slope of the surface topography and,
consequently, of the water table gradient. The lower topographic slope
within the mire itself leads to a rapid decrease in the vertical hydraulic
gradient towards the centre of the mire; however, the model suggests that groundwater seepage occurs throughout the mire even in the driest month. The
extent of the area where groundwater seepage occurs peaks in January, when
it extends largely to mineral soils surrounding the mire and to the small
non-channelised valleys and slope breaks in the upper part of the catchment. This area progressively shrinks, reaching its minimum extent in September. It is at this time of the year that the agreement between spatial patterns
in simulated groundwater seepage and observed mire boundaries is highest
(kappa
Simulated mean annual groundwater table depth and groundwater seepage rates in January and September (2001–2013). SZ: saturated zone; OL: overland.
Figure 12 shows time series of spatially averaged monthly mean groundwater table depth and selected water balance items in the peat soils. Groundwater
table depth shows low inter-annual variability in winter, being consistently
within a few centimetres of ground level. Inter-annual variability is larger
in summer. The grey bands in Fig. 12 highlight summers with the deepest
(2005) and shallowest (2007) mean groundwater tables. These are concomitant
with unusually low and high groundwater upwelling rates and overland
boundary inflows, which themselves followed relatively dry and wet winters
and springs, respectively. Evapotranspiration shows relatively little
inter-annual variability, which can be explained by the relatively shallow
groundwater table depth within the mire and the absence of water
availability limitation. Groundwater upwelling and overland boundary inflow
are highly correlated (Spearman's rank correlation
Approximate significance of smooth terms in the full generalised additive model.
GWT
Monthly means of simulated groundwater table depth and selected water balance items in the peat soil domain (2001–2013). Grey bands highlight two contrasting summers further described in Sect. 3.5.
Variation partitioning shows that, out of the four explanatory variables in the reduced model, groundwater upwelling explains the largest proportion of variation in groundwater table depth overall (70.3 %; see Fig. 14). However, due to collinearity, a large proportion of this variation is jointly explained by other variables, in particular groundwater table depth in the preceding month. Evapotranspiration explains 35.0 % of the variation in groundwater table depth, with 14.0 % of this being attributable to this variable independently of others.
Smooth terms of the final generalised additive model of simulated
monthly mean groundwater table depth on the response scale. The curve shows
the fitted value of the response, conditional upon the other explanatory
variables being held at their sample mean. The shaded area is the
approximate 95 % confidence interval. The time series used are simulated spatially averaged monthly means in the peat soil domain for the period
2001–2013: groundwater table depth (GWT), groundwater table depth in the
preceding month (GWT
Venn diagram of the variation partitioning results. The diagram
shows the proportion of deviance in groundwater table depth (GWT) explained independently or jointly by groundwater table depth in the preceding month (GWT
Groundwater modelling in hard-rock regions is often described as difficult due to the heterogeneity of the substratum (e.g. Levison et al., 2014;
Lubczynski and Gurwin, 2005; Neuman, 2005; Singhal and Gupta, 2010). In the
Dauges catchment, poorer model performance in some locations (e.g. D4 and
D5), an apparent trade-off between performance at dipwells located on three
opposite margins of the mire (D3, D9, D10, D24 and D25), and detailed
analysis of ERT results and geological drilling logs (Duranel, 2015) suggest
that, at the very local scale, representation of the fissured zone as an
homogenous entity may be an over-simplification. This is particularly the
case for its hydraulic conductivity. The model was moderately sensitive to
vertical saturated hydraulic conductivity of the fissured zone. Although the
calibrated value was
In comparison to some long-term hydrological modelling studies of larger catchments, the hydrological time series assembled to calibrate and validate our model may appear relatively short. However, calibration and validation periods covering a total duration of 3 years (and often much less) are the norm rather than the exception in physically based hydrological modelling studies of wetlands (e.g. Ala-aho et al., 2017; Armandine Les Landes et al., 2014; Haahti et al., 2016; House et al., 2016a; Levison et al., 2014; Li et al., 2019; Quillet et al., 2017; Thompson et al., 2004). In many cases this is the result of the unfortunate exclusion of wetland environments from formal hydrometric networks (e.g. Hollis and Thompson, 1998) despite their ecological and socio-economic significance.
During calibration, the bed resistance of channels represented in MIKE 11 had to be increased substantially and beyond tabulated values given in the literature (Dingman, 1994) to improve performance with regard to discharge at the Pont de Pierre mire outlet. This reflects issues with the kinematic routing representation of discharge and stage within the mire. The low slope gradient, the shallow and ill-defined channel, the presence of shrubs and dense vegetation within the channel, and, in places, the subterranean course of the stream make the flow sluggish and create constant backwater effects. As a consequence, flooding and exchanges through highly permeable alluvial gravel deposits below the peat could not be accurately modelled, as evidenced by the poor performance for dipwells D8, D16 and D17. However, as shown by the good performance obtained for dipwells D18 and D19, the impact of the poor simulation of stream stage does not extend beyond the area where substantial alluvial deposits occur.
Mires are frequently characterised by a hummock and hollow micro-topography
caused by peat accumulation and plant growth rates that vary at a very small
spatial scale (Belyea and Clymo, 2001; Malmer et al., 1994). Purple moor-grass (
During model calibration, the peat specific yield and available water capacity had to be reduced substantially to 0.05 to reproduce observed groundwater table dynamics. These values are lower than those generally cited in the literature even for sapric peat. For example, based on a literature review, Letts et al. (2000) provided a mean value of 0.13 for the specific yield and 0.49 for the available water capacity of sapric peat. Specific yields similar to the value we obtained by calibration have, however, been reported elsewhere in degraded peatlands. Boelter (1968), Malloy and Price (2014) and Mustamo et al. (2016) measured specific yields of 0.08 in decomposed peats in northern Minnesota, circa 0.05 in a cut-over bog in Canada and 0.067 in a cultivated peatland in northern Finland, respectively. It is well known that peat physical properties are highly variable (Boelter, 1964, 1968; Brandyk et al., 2002; Gobat, 1990; Gobat et al., 1986; Grosvernier et al., 1999), and the lower specific yield produced by model calibration at the Dauges site may result from compaction by cattle or from a higher humification rate caused by drier climatic conditions encountered in central France compared to more northern latitudes where mires are more common and to which most of the available literature refers. The low calibrated value for available water capacity is more puzzling as it is an order of magnitude lower than those usually reported for peat soils (Letts et al., 2000). Very low values have sometimes been reported from degraded peat soils (Mustamo et al., 2016). It may also be an artefact of the two-layer unsaturated zone/evapotranspiration model, whereby the allowable range for soil moisture in the unsaturated zone has to be fully depleted before evapotranspiration can proceed from the saturated zone (DHI, 2009c) and impact groundwater levels. In wetlands where the groundwater level is usually within the root zone, representation of transpiration as a process drawing water from both unsaturated and saturated zones simultaneously might be more realistic.
Groundwater table depths within peat soils and discharge at the Pont de
Pierre mire outlet were moderately sensitive to the saturated hydraulic
conductivity of the peat layer. The calibrated value was almost identical to
the median value obtained from slug tests in the catotelm (
This study demonstrates that groundwater upwelling from underlying mineral formations can be a quantitatively important and functionally critical element of the water balance of valley mires in granitic headwater catchments. Runoff inputs to the Dauges mire (a substantial part of which actually originates from edge-focused groundwater seepage just upstream of the mire boundaries) are quantitatively slightly larger than precipitation and groundwater upwelling. However, both runoff and precipitation do not contribute substantially to the water balance of the peat layer since they quickly leave the mire as saturation-excess runoff due to the shallow water table. Water within the saturated peat overwhelmingly originates from groundwater upwelling from the underlying weathered granite formations. At the monthly timescale, groundwater upwelling explains the largest proportion of the variation in groundwater table depth within the mire. The influence of hydrological processes occurring in mineral formations upslope of the mire on the mire's water balance is further demonstrated by the much higher sensitivity of simulated groundwater levels in the peat to parameters describing hydrodynamic properties of weathered granite formations and mineral soils compared to those describing peat soils. It is very clear that there is a high degree of hydrological connectivity between weathered granite formations within the mire catchment and the mire itself.
There is a close match between spatial patterns of simulated groundwater
upwelling and seepage and the observed distribution of mire habitats. Mire
habitats seem to be particularly associated with areas where the long-term
mean seepage rate in September is greater than zero. Good agreement between
spatial patterns in wetland vegetation and groundwater upwelling and seepage
has been reported for a number of wetlands that depend on large aquifers in
sedimentary contexts (Faulkner et al., 2016; Gerla, 1999; Reeve and Gracz,
2008), but much more rarely in hard-rock regions without substantial sedimentary cover. As far as we are aware, the only study demonstrating this
inter-relationship in a granitic context is that of Ala-aho et al. (2017).
Using an integrated physics-based hydrological model of the
3.2 km
The findings from this study have important management and legal implications. Indeed, the potential for anthropogenic pressures on groundwater bodies to indirectly impact groundwater-dependent terrestrial ecosystems (GWDTE) is now recognised in water resource legislation in many countries. In the European Union for instance, the Water Framework Directive (2000/60/EC) requires member states to ensure that water bodies, including groundwater bodies, achieve `good ecological status' which comprises a combination of good chemical and quantitative condition. Annex V of the Directive states that a groundwater body cannot achieve good status if the water table level is subject to anthropogenic alteration that would result in significant damage to terrestrial ecosystems that depend directly on the groundwater body. In practice the status of potential GWDTEs, in particular in hard-rock upland regions, is rarely or only superficially considered due to lack of understanding of the (i) impact of anthropogenic pressures on upland hard-rock aquifers; (ii) hydrological connectivity between hard-rock aquifers and terrestrial ecosystems; and (iii) hydrological controls on ecological conditions (Vernoux et al., 2010; Whiteman et al., 2010). Our study demonstrates that the Dauges mire depends directly on groundwater upwelling from the underlying mineral formations. Even though further work is required to upscale our conclusions, there is no reason to believe that the Dauges mire and its catchment are a unique case among valley mires in granitic headwaters. On the contrary, the Dauges site is representative of valley mires across the granitic Massif Central, and shares many similarities with other such systems found across the world, in particular throughout the Variscan Belt (Etlicher, 2005; Godard et al., 2001; Valadas, 1984, 1998). The results of this study show that the ecohydrological conditions of valley mires must be considered when assessing the status of hard-rock groundwater bodies under the Water Framework Directive.
In common with the Dauges mire, many headwater valley mires in hard-rock regions are designated for nature conservation under national or international legislation (Muller, 2018). Anthropogenic pressures on groundwater recharge within the catchments of these mires must be considered when defining conservation measures for these habitats or undertaking environmental impact assessments of policies and planning decisions. Article 6 of the EU 92/43/EEC Habitats Directive states that `any plan or project […] likely to have a significant effect thereon, either individually or in combination with other plans or projects, shall be subject to appropriate assessment of its implications for the site in view of the site's conservation objectives'. The present study demonstrates that any activities likely to result in substantial changes in groundwater flow within weathered granite formations upstream of the mire, even beyond SAC boundaries, do fall within the remit of this article. Examples of such activities include changes in land use such as large-scale afforestation and groundwater abstraction. Afforestation has, for example, been widespread within the Massif Central (Derrière et al., 2013; Dodane, 2009). The resulting changes in interception and evapotranspiration, with consequent impacts on surface and subsurface hydrological flow paths (Andréassian, 2004; Farley et al., 2005; Katzensteiner et al., 2011), may have had impacts on the region's mire ecosystems that should be further investigated.
Assuming that the relative importance of precipitation, runoff from mineral soil and groundwater seepage in the mire's overland component water balance reflects their relative contributions to overland inputs to water courses, 23 % of stream discharge at Pont de Pierre is water originating from the underlying mineral formations that has seeped through the peat layer. Excluding evapotranspiration, 46 % of water that flowed out of the saturated zone through seepage to overland flow, direct seepage to watercourses and groundwater flow out of the area upstream of Pont de Pierre did so through the peat layer. The vast majority (96 %) of the stream discharge at Pont de Pierre flowed at some point through the mire, either as overland flow or saturated flow. These figures highlight the importance of the mire as an interface between weathered granite formations and watercourses and hint at its potential role in mitigating the transfer of potential harmful compounds such as arsenic or uranium derived from the local geology (Mauroux et al., 2009) to water courses. Due to their organic soils and reductive conditions, pristine mires have been shown to constitute major sinks for almost all metallic and metalloid contaminants, and to efficiently reduce loads in downstream watercourses (Brown et al., 2000; González Acevedo et al., 2006; Lidman et al., 2012; Sobolewski, 1999). However, most of these studies have focussed on load reduction in watercourses as they flow through riparian mires, and further research is needed on load reduction in groundwater as it seeps through groundwater-fed mires.
Using an integrated MIKE SHE/MIKE 11 hydrological/hydraulic model, this study has shown that groundwater upwelling from weathered granite formations is a quantitatively important and functionally critical element of the water balance of the Dauges mire, a headwater valley mire in the granitic uplands of the French Massif Central. This is contrary to the “impermeable bedrock dogma” that has long been the dominant view in hard-rock hydrogeology and hillslope hydrology (Tromp-van Meerveld et al., 2007). Model performance in terms of simulated water table levels and stream discharge was at least satisfactory and in many locations good to very good. It included the replication of the seasonal rise and fall in groundwater levels. Based on model simulations, groundwater upwelling provided 27.1 % of total long-term inflows to the mire, rising to 37.2 % in September when total inflows are small. Furthermore, groundwater upwelling provided 92.4 % of total long-term inflows to the mire saturated zone. Variation partitioning showed that groundwater upwelling explained 70.3 % of the variation in monthly mean groundwater levels within the mire, although a large proportion of this is shared variation jointly explained by other variables, in particular groundwater table depth in the preceding month. The distribution of mire habitats within the catchment appears to be strongly controlled by groundwater upwelling. It was shown to closely coincide with areas where the simulated long-term mean groundwater seepage rate in September, the driest month associated with the lowest water table levels, is greater than zero. Model results demonstrate that the Dauges mire is a groundwater-dependent system closely connected to the granitic aquifer. This has important implications for the management and conservation of this and other similar ecosystems. However, whilst the Dauges mire is representative of acidic valley mires commonly found in many Variscan uplands, further work is required to upscale the conclusions of this study and to investigate how groundwater upwelling rates influence the ecology and biogeochemistry of such systems.
MIKE SHE/MIKE 11 is a commercial modelling software produced by DHI. Demo versions are available from
With the exception of secondary data derived under licence from third-party datasets, data used in this study are available upon reasonable request from the first author (arnaud.duranel@univ-st-etienne.fr). The meteorological data can be obtained from Météo-France.
The supplement related to this article is available online at:
AD designed and managed the study, installed the hydrometeorological network, collected all the data with the exception of manual checks, developed the geological conceptual model and MIKE SHE/MIKE 11 models, interpreted the results and wrote the first draft of the manuscript. JRT, HB and HC were AD's PhD supervisors and provided advice on the research design and on data analysis. JRT provided substantial expert advice on MIKE SHE/MIKE 11 and edited the initial manuscript. HC was principal investigator on the FEDER Loire and AELB grants that funded the hydrometeorological equipment and most field expenses. PD manages the Dauges National Nature Reserve. He contributed to the maintenance of the hydrometeorological network, completed regular manual checks, and provided substantial technical and logistical support during fieldwork. SG provided the ERT equipment, processed the ERT data and contributed to their interpretation. RW made a major contribution to the interpretation of the geological and geomorphological data and to the development of the geological model.
The authors declare that they have no conflict of interest.
The Conservatoire des Espaces Naturels de Nouvelle-Aquitaine and the Réserve Naturelle Nationale de la Tourbière des Dauges dedicated staff time to the collection of data used to develop the hydrological model. Météo-France provided part of the meteorological data. The Institut National de l'Information Géographique et Forestière, Areva (now Orano), the Conservatoire des Espaces Naturels de Nouvelle-Aquitaine and the Réserve Naturelle Nationale de la Tourbière des Dauges provided geological, topographic and cartographic data.
This research has been supported by a UK Natural Environmental Research Council doctoral studentship to Arnaud Duranel (grant no. NE/H525203/1) and a grant to UMR 5600 CNRS EVS from the European Fund for Regional Development (FEDER Loire) and the Agence de l'Eau Loire-Bretagne, obtained as part of the Plan Loire Grandeur Nature managed by the Établissement Public Loire (
This paper was edited by Philippe Ackerer and reviewed by two anonymous referees.