Groundwater Contribution to Sewer Network Baseflow in an Urban Catchment-Case Study of Pin Sec Catchment, Nantes, France

: Sewer systems affect urban soil characteristics and subsoil water flow. The direct connection observed between baseflow in sewer systems under drainage infiltrations and piezometric levels influences the hydrological behavior of urban catchments, and must consequently be considered in the hydrologic modeling of urban areas. This research studies the groundwater contribution to sewer networks by first characterizing the phenomenon using experimental data recorded on a small urban catchment in Nantes (France). Then, the model MODFLOW was used to simulate the infiltration of groundwater into a sewer network and model dry weather flows at an urban catchment scale. This application of MODFLOW requires representing, in a simplified way, the interactions between the soil and the sewer trench, which acts as a drain. Observed average groundwater levels were satisfactorily simulated by the model while the baseflow dynamics is well reproduced. Nonetheless, soil parameters resulted to be very sensitive, and achieving good results for joint groundwater levels and baseflow was not possible.


Introduction
Urbanization modifies land use and affects soil, sub-soil and subsurface processes in different ways. Urban features as well as surface and underground infrastructure can have a strong impact on groundwater levels. Some causes explaining a reduction in groundwater include an infiltration decrease due to additional imperviousness, groundwater pumping for various urban water uses, and groundwater flow into drainage trenches [1]. On the other hand, leakage from water supply and waste water networks becomes a source of recharge for urban groundwater [2][3][4][5]. Furthermore, leakage from waste water systems is also a possible source of groundwater contamination [6]. Despite aquifer levels and stream flow having proved to be related in various rural contexts [7,8], few studies have focused on examining this relationship in urban catchments [5,[9][10][11]. Soil water infiltration in sewers, however, is a phenomenon affecting urban hydrology, sewerage and waste water treatment plant management [12,13]. Initially, the presence of soil water in wastewater sewerage was revealed by [14,15], when flow variations in separate wastewater sewerage during wet weather periods were commonly attributed to inappropriate or irregular stormwater runoff connections. These authors showed that groundwater seep into the sewers through defective cracks once the water table level reaches the depth at which sewers are buried. This phenomenon has been recently reinforced by [13] for several Flemish catchments; some interesting stable isotopes methods have been used to detect and estimate groundwater intrusions in sewers in Nancy, France [16] and in Brussels, Belgium [17] and temperature-based methods in Trondheim, Norway [18]. Similarly, [9] showed that soil water may explain the variations of runoff coefficients observed on a small urban catchment in the city of Nantes, France, as groundwater drainage flow becomes significant once the water table exceeds a given threshold. The flowrate in sewers is the sum of runoff during rain events, wastewaters from housing (for separative wastewater or combined sewers), and a groundwater infiltration component, representing the baseflow in sewers; this component is often part of "extraneous" water in sewers. Interestingly, [19] suggest the appropriate terminology of urban "karst" to denote the multiple water soil and pipe interactions below the surface. Despite the relevance of improving our understanding of the impact of groundwater drained by urban sewer systems, only a few efforts are currently devoted to the problem. Groundwater drainage by sewers is often noted in literature [20], although is rarely quantified when assessing the urban water budget. Only a few papers deal with the determination of the urban water budget, likely due to the lack of reliable data of all the water budget components (i.e., rainfall, evapotranspiration, surface and groundwater flow rates). This lack of knowledge must be addressed as urban stormwater practices increasingly consider innovative infiltration-based technologies and approaches to mitigate the hydrological impacts of urbanization [21]. Thus, a better understanding and quantification of the fate of urban soil water becomes essential [22,23].
The first attempts to represent the interaction between surface water and the aquifer in urban integrated hydrological models were simple and conceptual. Aquacycle [24] and the Storm Water Management Model (SWMM, [25]) are examples of such models. In a comprehensive review of ten stormwater models, [26] identified the models able to simulate the groundwater baseflow component, mostly through a conceptual linear reservoir. Some of the more recent physically-based modeling efforts consider soil drainage by sewers. An integrated sewer-aquifer model was developed [27] to test the effects of future buried drains on both groundwater flow and sewer infiltration phenomena. This integration proved to be relevant to describe the groundwater-sewer interactions. The Network Exfiltration and Infiltration Model (NEIMO) model [28] was developed in close connection with already available groundwater models such as MODFLOW [29] or FEFLOW [30], to simulate infiltration and exfiltration processes. This modeling approach requires developing an integrated modeling framework like Urban Volume and Quality (UVQ) developed by CSIRO, Australia [31,32]. MODFLOW was used [33] to analyze the relevant infiltration parameters at a local scale, and compared a 1D-infiltration approach and a MODFLOW modeling approach to model groundwater infiltration at a larger scale. Overall, coupling groundwater models with the simulation of infiltration/exfiltration processes is a major challenge, especially due to the high uncertainty associated with these models [34]. The coupling of a hydraulic model was realized by [35] with both a groundwater model and a sewer failure estimation approach to identify the sections potentially affected by infiltration in coastal urban areas. Satisfying results were obtained by [11] with regard to the groundwater infiltration by coupling two commercial models MIKE URBAN for simulation of sewer flow and MIKE SHE for simulation of groundwater transport. The Urban Runoff Branching Structure Model (URBS-MO) [36] was developed to represent surface and subsurface water flows while focusing on the impact of sewer networks. Although the model simulates groundwater drainage by sewer, real observations at a local scale to test this component of the model were not available. Finally, [5] coupled the hydrological/water management WEAP model and MODFLOW to verify a strong stream-aquifer interaction in areas with shallow groundwater, as well as quantify the local recharge associated with pipe leaks and inefficient urban irrigation. Overall, the literature review shows that much more experimental data and modeling using comprehensive 3D modeling tools such as MODFLOW, are needed to better understand and quantify the impact of groundwater infiltration and groundwater levels on low flows in urban environments.
This study investigates the interactions between urban groundwater and baseflow in wastewater and stormwater sewer systems, and its first goal is to better understand their role on the hydrological budget of an urban catchment. These interactions are analyzed and quantified using a field study conducted in the "Pin Sec" catchment, located in the city of Nantes (France) [37]. The second objective consists of testing the ability of a groundwater modeling tool to mimic these interactions. The results are used to build a numerical model focused on simulating the hydrological behavior of the urban soil and low flows in the sewer system. The outline of this paper is as follows: the next section presents the case study and the methods adopted in this work, which focus on using MODFLOW to model baseflow and interactions between the sewer systems and groundwater at catchment scale. The results section presents and compares observed and simulated groundwater levels and baseflow discharges. Finally, the last section summarizes our results and highlights the main conclusions.

Case Study
The experimental area is part of the Observatoire Nantais des Environnements Urbains (ONEVU) initiative, devoted to long-term monitoring of water and pollutant fluxes and soil-atmosphere exchanges in urban settings [37]. The Pin Sec catchment is located in the East of Nantes (France) between the Loire and the Erdre rivers. The area is under an oceanic climate, with frequent but not very intense rainfalls. The Pin Sec neighborhood developed between 1930 and 1970 is mainly composed of single and multi-family housing ( Figure 1). This 31 ha catchment (called urban catchment) belongs to a larger hydrogeological catchment of 120 ha, both of them are represented on Figure 1. This catchment has a gentle slope, the highest altitude of the catchment is 28 m asl and the lowest is 13 m. The Pin sec catchment has an imperviousness of 45% and is equipped with a 50 years old separate sewer network. The wastewater and stormwater sewer systems have total lengths of 7.3 and 4 km, and mean depths of 2 and 2.9 m respectively. Although the wastewater sewer network is denser than the stormwater system, both are mainly superimposed, and both present cracks or faulty sealing joints. Both system outlets, while collecting water from the same geographic area, are not located in the same place. The stormwater system is connected to the Gohards river, whereas the wastewater sewer system drains into a downstream combined sewer system. Gohards river stream is an old river which was buried between 1945 and 2012 and which has been re-opened through a renovation urban planning of this area in 2012. It is a perennial stream. Recorded data include (1) 5 min rainfall records from 1999 to 2010, deduced from the average of 3 rain gauges covering the study area, (2) 5-min flow discharge records at the outlet of both the stormwater and the wastewater systems from September 1st, 2006 to August 31st, 2008, (3) 20-min groundwater levels measured with pressure sensors in eleven piezometers located throughout the catchment from September 1st, 2006 to August 31st, 2008 (one piezometer, PZGN being only used to map the groundwater level contours), and (4) other meteorological records as well as evapotranspiration estimated with the Penman method, measured 12 km away from the catchment outlet, from 1999 to 2010. We focus our analysis on the main period between September 1st, 2006 and August 31st, 2008 (i.e., two hydrological years), but we also used year 2002 for some modeling tests as rainfall characteristics were representative of the average conditions in the area.

Geology and Hydrogeology
The city of Nantes has developed on the Armorican massif, and especially in the south Armorican shear zone. The geology is mainly composed of a mica-schists bedrock and a stack of alluvial deposits along the Loire River. On the Pin sec catchment, the mica-schists bedrock is covered with altered mica-schists, silty eolian deposits and alluvial deposits along the Gohards stream ( Figure 2a) [38,39]. The piezometers are mainly located in the top layers and especially within altered mica-schists, with a drilling depth varying from 4 to 8 m ( Figure 2b); for two of them, the depth was limited due to hard non altered micaschists. The hydraulic conductivities estimated from in-situ measurements by water bail tests [40] using the Hvorslev method [41] varies from 8.5 10 −8 to 1.810 −5 m.s −1 [42]. Piezometer PZGN is out of the hydrogeological catchment and was not considered for further results. From now on groundwater level will designate the level of the saturated zone, which is equivalent to the local groundwater level in our case study. The piezometers were used for mapping the behavior of the water table on the catchment. Groundwater level contours were interpolated using inverse distance weighting (IDW) techniques with Vertical Mapper© in Mapinfo platform, and used to determine groundwater main flow directions. The groundwater levels of the high water table periods are presented in Figure 3. Note that the upstream condition of the water table is given by two piezometers (PZCRI and PZGO) located at the upper limit of the catchment, and that groundwater moves in the northwest-southeast direction towards the stream outlet following the surface slope. The low water table condition is characterized by a similar profile.

Baseflow in Sewer Networks
Both waste-and stormwater systems outlets are equipped with continuous flow rate measurement devices that use triangular weirs. Nonetheless, measuring the lowest flow rates within sewers has proved difficult as dry-weather flows do not necessarily reach the minimum water depth (~3 cm) needed for proper velocity measurement. Diurnal patterns of the wastewater flow rates show a typical domestic double-bump flow during dry weather time (Figure 4), and an increase of this flow during rain weather, due to parasitic waters (due to both wrong sewer connexions and groundwater infiltration). The stormwater flow rate is typical for urban catchments, with a fast hydrological response related to the runoff produced by impervious areas. In this study, baseflow is considered on a daily basis. Stormwater baseflow is basically the minimum flowrate of each day. Because of the small size of the catchment, minimum wastewater flows corresponding to water draining from the water table are assumed to occur in the middle of the night [43][44][45]. Thus, we assume that the flow rate at 4:00 A.M. is representative of this drainage flow. Those days in which rain takes place are removed from the wastewater and stormwater discharge records to assess baseflow at the outlet of both systems during the study period (i.e., 65% of the time).

Groundwater Modeling at the Urban Catchment Scale
Our aim in proposing a model to represent soil-sewer interactions is to improve our understanding of water fluxes in urban catchments. We use the groundwater model MODFLOW [29] to simulate these interactions because of (1) its widely-used extensive modeling capabilities [46], (2) its drainage modeling options to account for soil discontinuities using fine discretizations, and (3) the possibility of easily linking MODFLOW with other models. Few studies, however, have used the model in urban groundwater modeling [5,33,[47][48][49], and in doing so, quite large grid resolutions (100-500 m) are adopted. In our study MODFLOW is used at a much finer resolution to represent sewers or buried trenches and their discontinuities in detail. Simple application conditions have been retained, because the use of MODFLOW in an urban context is not frequent and a cautious approach has been adopted. Moreover, the urban soil is both vertical and horizontal heterogeneous, and therefore not well documented.

Main Modeling Principles and Assumptions
MODFLOW uses a finite-difference approximation and an implicit scheme [29] to solve the 3D groundwater flow diffusion equation. In our implementation, the geological features of the catchment are deduced from the geological characteristics found in the case study (i.e., silt, altered and non-altered mica-schist). Two soil layers are defined from base to top: (1) a basement layer consisting of non-altered mica-schists; (2) a topsoil layer consisting of silty eolian deposits and altered mica-schists on the main part of the area, and alluvial deposits on the downstream part; this topsoil layer contains the sewer trenches.
Three MODFLOW modeling options were tested to represent the sewer trench, as this element cannot be explicitly implemented in the model: (1) a basic modification of the hydrodynamic properties of the cells containing the trench; (2) a field trench parallel to the surface slope with a small drainage capacity, using the drain package (DRAIN); and (3) a representation of the sewer network using a river with the river package (RIV).
A preliminary sensitivity analysis was carried out at a street scale for a 100 × 100 m unit element, both in steady state and transient conditions (not presented here). This analysis helped to understand the impact of the modeling configurations (grid resolution and the sewer trench representation option) on both groundwater levels and flows simulation. Moreover, the influence of the main parameters used in the model (i.e., sewer trench conductance, hydraulic conductivities, specific yield and specific storage) was analyzed. The choices about the parameters or configurations described below were not adopted to replicate experimental results, but to retain realistic characteristics of the subsurface processes to be represented.
The grid resolution in our analysis varied between 1 and 10 m and affected mainly the water flows. Finally, a 4 m grid cell resolution was chosen. This size is quite large to represent the sewer trenches encompassing sewer pipes, but is a good compromise given the area of the study catchment (120 ha). Along the three tested options, the option using the DRAIN package to represent the sewer trench proved best to represent the saturated zone drainage by a sewer, because it can better reproduce the decreasing groundwater near the trench. Specific yields and conductivities are the main sensitive parameters, influencing both the groundwater levels and the flows. Due to the low impact of anisotropy in our sensitivity analysis, conductivities are considered homogenous in the three directions. The water flux increases with the hydraulic conductivities, while the groundwater level decreases. The initial hydraulic conductivities are estimated from in-situ measurements described above ( §2.1.2). The influence of the specific yield of the topsoil layer is significant for the sewer trench outflow: the downstream flow reduces by 7% when specific yield decreases from 20% to 5%, which are typical values assessed by [50]. The field drainage system representing the sewer trench is characterized by a drain hydraulic conductance (L 2 /T) (conductance thereafter), which depends a priori on the material and characteristics of the drain. Since the conductance is usually unknown, it must be estimated by calibration [49,51]. We first assessed the range of variation of the conductance from the sewer baseflow observations. As discussed by [51], the discharge rate to drain cells may be calculated as the product of the conductance and the head gradient (h-d), where h is the hydraulic head and d is the elevation of the drain. We integrated this relation at the whole sewer network scale, which allows estimating a global conductance associated with the total groundwater flowrate in sewers. Assuming that the baseflow in sewers is only due to groundwater infiltrations, the application of this relation for the high water level period during year 2007 and for the total length of flooded sewer pipes generates a conductance of ~0.9 m 2 /day. This parameter only changes the water flows and its effect on groundwater levels is less significant. In the end, the specific storage slightly affects the modeling results, as it ranges theoretically from 3.3 × 10 −6 m −1 (rock) to 2 × 10 2 m −1 (plastic clay) [51]. A medium value of 1 × 10 −5 m −1 was adopted initially.

Urban Groundwater Modeling at the Catchment Scale and Application to the Case Study
The MODFLOW model is used on the Pin Sec catchment to assess soil water-sewer interactions at the urban catchment scale. Furthermore, certain simplifications not considered in the sensitivity analysis for the unit element were adopted. The two soil layers are represented in this way: (1) the basement layer is located in the entire catchment at a mean depth of 18 m, and (2) the topsoil layer is located in the valley area around the old stream bed, with a mean depth of 2 m. Because both the stormwater and wastewater systems are often placed within the same trench below the street surface, both networks have been combined to represent the draining trench in the soil. Thus, the field drain represents both stormwater and wastewater networks.
Boundary conditions are stipulated as follows and summarized on Figure 5: (1) a zero-flux condition is defined in the upstream boundary of the hydrogeological catchment by considering no-flow cells and (2) a river boundary condition representing the Gohards stream is defined as the downstream condition by using the River package (RIV) on a 200 m long stream; it is a head-dependant boundary condition. The combined 'field drain', representing the sewer systems, drains into this stream. A one hydrological year simulation in transient conditions is performed on the catchment using a daily time step. The transient simulation is initialized using the steady state case calibrated to the field data observed the first day of simulation. The recharge is calculated on a daily basis using the urban surface hydrological model URBS-MO, already successfully applied to urban catchments in the area to simulate the various components of the water budget [36]. The recharge is assumed to be uniform throughout the catchment, and equal to the infiltration calculated by URBS-MO (i.e., rainfall intensity minus both the actual evapotranspiration and surface runoff). Rainfall and potential evapotranspiration data available for the Pin Sec catchment were used with URBS-MO. The initial set of parameters for the simulation is deduced from the observations and a first steady state simulation conducted on the catchment (Table 1). Specific yield and specific storage values are the same as those used for the sensitivity analysis; they vary in space along with the two geological layers defined in this catchment. Due to the lack of information about the spatial distribution of cracks and faulty sealing joints in sewers and to an homogeneous age of the pipes, the drain conductance is assumed to be uniform and the initial value is deduced from the groundwater baseflow observed in the catchment sewers presented in next section.

Groundwater Dynamics and Experimental Relationship Between Groundwater Level and Baseflow
The temporal dynamics of groundwater levels are presented in Figure 6 for two of the piezometers; more details may be found in [39]. The mean difference between the high and low water Groundwater infiltration in sewer networks depended on several factors such as the location and shape of the water table, the depth of the sewers, and their conditions (e.g., presence of defects). [9] collected information in a small urban catchment in the metropolitan area where the Pin Sec catchment is located, and noticed that groundwater drainage flow became significant once the water table reached a depth of 1.5 m below the surface, and increased as the water table rose. In fact, the mean sewer depth in the area is ~1.2 m below the surface. .
The connection between the average groundwater level of the catchment H and daily baseflow Q in the sewer systems is evident from Figure 7. A specific piezometer (PZPS) located near the main sewer intersection and representative of the average behaviour of the groundwater level dynamics has been chosen; it is considered to be in the catchment center according to groundwater main flow directions. Figure 7a shows a Q~H 2 relationship for the baseflow in the storm sewer, which is typical of ideal field drains used in rural hydrology [52]. Groundwater drainage is observed at the outlet of the stormwater network during the entire observation period, which demonstrates that sewer infiltrations take place all year round. The behavior is similar for the wastewater sewers (Figure 7b). The wastewater baseflow barely varies during summer when the water table is low (between July and October). Baseflow becomes substantial again when the groundwater level for this piezometer PZPS exceeds a threshold value of ~16.75 m (i.e., 3.2 m below the surface). Wasteand storm-baseflows differ because the drainage density of the wastewater system is higher than that of the stormwater system. This is the reason why the Q~H 2 relationships are more consistent for the wastewater baseflow that the stormwater baseflow, as revealed by the determination coefficient variation R 2 (SW) or R 2 (WW) ( Table 2). Additionally, these relationships are generally more pronounced for the downstream piezometers, closer to the baseflow monitoring locations.   Figure 7 for piezometer PZPS).  Figure 8). This map was realized through the superimposition of groundwater level contours ( Figure 3) and the sewer depth levels deduced from GIS data. The total sewer length soaked during winter is 2.6 km for the wastewater sewers and 1.1 km for the stormwater sewers. Finally, the annual subsurface water volume drained into both sewer systems deduced from the sum of daily base-flows between September 2006 and August 2007 is 28% (wastewater sewer) and 14% (stormwater sewer) of the total rainfall respectively. This result proves that the base flow discharges in artificial sewer systems can be a major component of the urban water budget, as 42% of the total annual rainfall is drained to the sewer systems.

Comparison between Modeling Results and Observed Data at the Catchment Scale
Three metrics were adopted for the evaluation of the model: the root mean square error (RMSE), a Bias error, and the determination coefficient R 2 : where n is the number of time steps involved, 1 V and 2 V are the observed and simulated values, Absolute error (m) was used as well to estimate the difference between simulated and observed groundwater levels in specific piezometers.

Implementation and Calibration of MODFLOW at the Catchment Scale
In the MODFLOW model of the catchment, we first calibrated the hydraulic conductivity, specific yield and specific storage of each soil type by minimizing the global bias error between observed and simulated groundwater levels. Then we calibrated the drain conductance by maximizing the coefficient of determination R 2 between observed and simulated combined baseflows (i.e., stormwater + wastewater sewers baseflows) at the outlet of both sewer systems. The final calibration presented in Table 1 is achieved through a trial and error method. The obtained hydraulic conductivities are two orders of magnitude higher than the observed hydraulic conductivities using the piezometers. This result is consistent with previous statements following the application of other models to ditch-groundwater interactions [53], and confirms the difficulty to simulate groundwater flows and groundwater levels using a few specific point observations of the hydraulic conductivities. High soil heterogeneity, particularly in urban areas, cannot be well characterized using soil properties measurements sampled in a few locations, which leads to a poor simulation of the general groundwater behavior. Despite the differences between observed and calibrated conductivities, the observed differences between the different types of soil are consistent with the theory. Hence, the highest hydraulic conductivity is obtained for alluvial deposits, and the smallest one for mica-schist. Finally, the calibrated drain conductance is higher than initially estimated, which allows a better simulation of high values and the dynamics of the baseflow.
In order to assess the influence of the sewer system on the evolution of the groundwater levels, a complementary simulation of the catchment was done on the same period and with the same set of parameters by removing the sewer system and the corresponding field drain. This simulation exercise shows that the groundwater reaches the soil surface during the winter 2007-2008 on a significant part of the basin; that result clearly confirmed the importance of the drainage ensured by the sewer system [39].

Groundwater Level Distribution Assessment
The temporal and spatial simulation of groundwater levels is assessed using available piezometric records. Figure 9 illustrates the simulation of the isopiezometric contours for the high water table period (7th March, 2007), and can be compared to the observed isopiezometric contours at the same date shown in Figure 3. At this period, the absolute errors between simulated and observed values may be larger than 2 m in some piezometers, especially on the south catchment boundary. They are however better in the middle of the catchment and near the outlet. In addition, absolute errors are higher around the piezometers distant from the sewer trenches. Note that the simulated contours in Figure 9 are plotted on a surface area larger than the experimental one, as the simulated catchment exceeds the piezometric set boundary. Consequently, the simulated flow directions differ significantly from observed ones near the urban catchment boundary, but are rather similar within the urban catchment. Nonetheless, at this time the simulated saturation level is often higher than the observed one. The simulated vs. observed groundwater level dynamics is assessed by computing the RMSE, bias error and R 2 for all the piezometer locations, both for the calibration and validation periods (Table 3). R 2 values are often smaller than 0.5 for the calibration stage, but increase for the validation stage. In general, the model underestimates the groundwater level for half of the piezometers (PZCPS, PZCS, PZGO, PZUV, PZJV) and overestimates it for the others. The boxplot analysis of the bias error ( Figure 10) shows that the underestimation is mainly concentrated on the upstream part of the catchment (PZCRI, PZGO, PZCS, PZJV, PZUV), and might be explained by either the upstream boundary condition or the uniform recharge assumption. Indeed, the spatial analysis of the imperviousness coefficient of the Pin Sec catchment reveals that underestimated piezometric observations are located in more pervious areas within the catchment (i.e., green public parks). Hence, using a spatially distributed recharge could allow a more accurate modeling of the groundwater levels dynamics in these areas. Overall, the simulated groundwater levels compared properly with observed ones both for the calibration and validation periods (Figure 11), despite an underestimation of the model compared to the observed data at the beginning of the validation period, when a particularly dry autumn took place. In the absence of rain, the groundwater model simulates a groundwater level decrease higher than observed, which to some extent questions both the upstream boundary condition assumption and the recharge estimation. Moreover, the quite stable groundwater level observed during this period could be explained by a groundwater influx coming from upstream, which would disprove the assumed upstream boundary condition. Nonetheless, no observations outside the catchment are available to clarify this issue. Finally, drinking water leakage in the soil can occur, which would increase the recharge throughout the year and prevent the high decrease of groundwater levels observed in the simulation; this phenomenon has been discussed in several cities like Bucharest, Romania [4,49], St Louis, Missouri(US) [54] or Santiago, Chile [5]. Figure 9 shows a small curvature of the simulated groundwater level contours near the sewer trenches caused by the local groundwater drawdown. Groundwater decrease is quite significant in places where sewer density is high (i.e., central and north portions of the catchment). This phenomenon is enhanced when the water table rises above the sewer system, and becomes less significant in summer. For simulated data samples, groundwater drawdown is higher during the wet year 2006-2007 with a decrease of about 0.3 m near the sewer trench, and a zone of influence ranging up to 120 m away from the trench axis.

Baseflow Rate Assessment
The simulated soil water flux combines two components simulated in MODFLOW: groundwater flux within the soil and groundwater drainage inside the sewer trench. Both components flow in the same general direction forced by the catchment topography and the downstream condition imposed by the Gohards stream. Groundwater drainage may be compared with the baseflow discharge observed in both waste-and storm-sewer systems during the 2006-2007 simulation (note that the period 2007-2008 is not presented due to the lack of data). Figure 12 shows that the daily groundwater drainage component varies significantly, ranging from 50 m 3 /day during very dry periods (i.e., the end of summer 2007) to 2600 m 3 /day during the particularly humid 2006-2007 winter. Indeed, ~630 mm of rain (i.e., ~80% of the catchment mean annual rainfall) fell between November 2006 and April 2007 in the area. Unfortunately, a direct comparison between the simulated groundwater drainage flow and observed baseflow is difficult because of the lack of flow rate data. Considering the data actually available, the comparison focuses on the dynamic evolution of the baseflow. Simulated groundwater drainage and in-situ observations over the simulated period vary in a quite similar manner. Nonetheless the overestimation of the observed flows could not be reduced by further calibration, as drain conductance variations do not affect groundwater drainage significantly. As seen in the sensitivity analysis, soil hydraulic conductivities have more impact on the simulation than drain conductance. Because soil hydraulic conductivities affect not only groundwater drainage but also groundwater level, a simultaneous enhancement of both is not possible.  Special attention must be paid to the relationship between baseflow rates and groundwater levels. As discussed above and shown in Figure 7, both are strongly linked through a power-law function. Figure 13 overlaps simulated and observed values when considering the combined baseflow, and shows that the model can satisfactorily reproduce this groundwater level-baseflow relationship despite of the baseflow overestimation during high water table periods. Although only 32% of the experimental results during the 2006-2007 period are available simultaneously for both sewers, and apart from the low groundwater levels occurring in summer and discussed above, the validity of the relationship holds quite well.

Conclusions
This study presents an experimental analysis and modeling of the interactions between groundwater and urban sewer systems in a 31 ha urban catchment located near Nantes, France. Groundwater contributions to flow in stormwater and wastewater systems due to sewers' defects were deduced from observed flow rates, whereas piezometric records were used to characterize the dynamics of groundwater levels in the catchment.
Overall, the experimental analysis reveals a strong co-fluctuation of groundwater levels and sewer baseflow, which is more pronounced in downstream piezometers. This study highlights that two or three piezometers located near the catchment outlet are likely adequate to estimate the groundwater-baseflow relationship in an urban catchment. The total baseflow volume in wastewater sewers is larger than in stormwater sewers because of its larger coverage. On an annual basis, the total volume of soil water drained by both sewers is 42% of the total rainfall.
A MODFLOW model was implemented in the study catchment to simulate the groundwater dynamics and contribution to the sewer systems on a daily basis. This application proves the model ability to represent satisfactorily the spatial and temporal evolutions of groundwater levels, although low groundwater levels were underestimated. Spatial differences between modeled and observed groundwater levels in some locations could be explained by the spatial variability of the land uses and the urban soil configuration around each piezometer. The simulation of combined waste and storm water sewer baseflows is not entirely relevant because the model fails to reproduce the high values of baseflow. The calibration of the model lead to hydraulic conductivities values higher than those observed in the piezometers. In fact, the calibrated values are similar to those reported in the literature, and they are attributed to possible secondary porosities and hydraulic conductivity distributions comparable to those of a karstic system [55]. Hydraulic conductivity varies significantly in urban soils due to differences in their characteristics, the existence of different land-uses, and the presence of previous buried constructions.
More work is needed to improve our knowledge of urban soil water flow paths. Baseflow must be better measured, as flow meter sensors are not always reliable to be used with low flows. In addition, more observations at a local scale are needed to better understand groundwater-sewer interactions and improve large scale modeling efforts. Such data will come from either (1) physical models relying on better estimations of soil parameters, which can represent the influence area of a sewer trench, or (2) sensors located near the sewers. Groundwater models applied to urban catchments will require using a graded grid size with mesh refinements near sewers and more urbanized areas. Furthermore, better spatially distributed recharge estimations can also improve the performance of these models. This could be done by coupling the groundwater and surface water models, adopting better modeling frameworks like the modeling chain proposed by [56], or implementing other specific coupling modeling approaches. Finally, a better understanding of the occurrence and spatial distribution of leaks from drinking water pipes is also needed to improve the characterization of the recharge.
Author Contributions: This study was designed by F.R., A.-L.L.D., and H.A. The manuscript was prepared by F.R. and revised by H.A. and J.G. All authors have read and agreed to the published version of the manuscript.
Funding: This study received financial support from the French research programs ANR and INSU-EC2CO within the project Hy 2 ville, and from the Région Pays de la Loire. Funding from CONICYT/FONDAP grants 15110017 and 15110020 were also available.