Application of an improved global-scale groundwater model for water table estimation across New Zealand

Many studies underline the importance of groundwater assessment at the larger, i.e., global, scale. The large-scale models used for these assessments are often simplified and typically not used for smaller-scale, i.e., catchment-scale, studies, because hydrology and water policy are traditionally best constrained at the catchment scale, and because large-scale models are too uncertain for that scale. However, smaller-scale groundwater models can still have considerable uncertainty, especially in data-sparse areas. There is a potential for larger-scale models to constrain the uncertainty for small-scale models. That 5 is because they can provide an extra source of information in data-sparse areas, such as the initial estimate of hydraulic head. Large-scale models, often quick and simple, can thus take away some of the computational burden of local and more sophisticated applications. The problem of this approach is that model uncertainty of large-scale models is often too large, because the quality of their, coarse and global-scale, input data is large, and often inconsistent with the input data of local models. What is needed is an approach where large-scale and local models can meet ‘in the middle’. 10 This study uses an existing, global-scale, groundwater flow model. It feeds that model with national input data of New Zealand terrain, geology, and recharge. It then builds the first New-Zealand national-scale groundwater model. The resulting nationwide maps of hydraulic head and water table depths show that the model points out the main alluvial aquifers with fine spatial detail (200m grid resolution). The national input data and finer spatial detail result in better and more realistic variations of water table depth than the original, global-scale, model outputs. In two regional case studies in New Zealand, the hydraulic 15 head matches the available groundwater level data well. The nationwide water tables show that the model is mostly driven by the elevation (gravity) and impeded by the geology (permeability). The use of this first New Zealand-wide model can aid in provision of water table estimates in data-sparse regions. The national model can also be used to solve inconsistency of models in areas of trans-boundary aquifers, i.e., aquifers that cover more than one region in New Zealand. Shortcomings of the model are caused by the simplified model properties, but also by 20 the accuracy of input data. Future research should therefore not only focus on further improvements of model equations, but also improved estimation of hydraulic conductivity and the digital elevation model, especially in areas of shallow groundwater level. We further surmise that the findings of this study, i.e., application of a global-scale models at smaller-scales, will lead to subsequent improvement of the global-scale model equations. 2 Hydrol. Earth Syst. Sci. Discuss., https://doi.org/10.5194/hess-2018-235 Manuscript under review for journal Hydrol. Earth Syst. Sci. Discussion started: 18 May 2018 c © Author(s) 2018. CC BY 4.0 License.


Introduction
Groundwater is a key water resource to many countries in the world, providing water for irrigation, domestic consumption, and industry.Groundwater also provides baseflow to streams and rivers, which, in drier times, helps sustain ecology when rainfall runoff is low.Existing studies underline the importance of global-scale assessment of groundwater (Wada et al., 2012(Wada et al., , 2014)), its role in a changing climate (Taylor et al., 2012;Green et al., 2011), and its current depletion (Gleeson et al., 2012;Richey et al., 2015).
Groundwater models used for global-scale assessments are simplified.Because these models have to cover the entire globe, they generally use coarse input data or embed a simplified model algorithm.For example, Fan et al. (2013a) apply a global model of water table estimates with grid cells of approximately 1 km, using simplified subsurface parametrisation and flow equations that calculate an equilibrium water table (explained in more detail in Sect.2).De Graaf et al. (2015) apply a model with approximately 10 km resolution and global-scale input data and also look for the equilibrium water table, albeit using a MODFLOW environment, with subsequent improvement to a two-layer MODFLOW application to better include confining layers (de Graaf et al., 2017).
Advanced numerical advanced groundwater flow studies are mostly at the catchment scale (e.g.Gusyev et al., 2013;Wang et al., 2008;Oude Essink et al., 2010), where hydrology is traditionally best constrained.Water policy also typically fits the catchment scale, and, in most countries, decision-making is based on hydrological studies of this scale.In New Zealand, that is not much different: water resources and allocation are managed by regional councils, and regions mostly have the same boundaries as the catchments.Advanced numerical groundwater flow models at the catchment scale in New Zealand have been developed by, e.g., White (2018), Toews et al. (2016), and Wöhling et al. (2018).
National-scale groundwater information in New Zealand is also important.Guidelines on groundwater allocation are defined by the national government (Ministry for the Environment, 2013Environment, , 2008)), thus requiring a nationwide overview.Furthermore, national-scale hydrology models in New Zealand (TOPNET: Bandaragoda et al., 2004) require groundwater information to calculate the interaction of groundwater with surface water, such as baseflow or flow loss from rivers.However, many challenges still need to be overcome in order to obtain a national overview of New Zealand's groundwater resources.Substantial parts of New Zealand are not covered by advanced numerical groundwater flow models, and the regions that do have this coverage do not have agreed-upon and consistent modelling approaches, data formats, and data availability (Westerhoff et al., 2018), causing "trans-boundary" issues, such as drinking water source protection and catchment boundary definition.Where some regions are relatively data-rich, others do not have enough information to obtain a detailed groundwater flow model.Development of individual advanced groundwater flow models for each region, compatible across regions, is recommended but would take many years to complete in the current policy framework, with significant efforts required in aligning regional and national policymakers.Also, the combined model runtime of those advanced, and therefore computationally demanding, models would be a large computational burden.
This study proposes a solution to obtain a nationwide overview of groundwater that bridges the gap between (too-)expensive advanced local models and (too-)simple globalscale models.It describes the development of a groundwater model at the national scale of New Zealand.Our model is inspired by the global-scale Equilibrium Water Table (EWT) model (e.g.Fan et al., 2013a), with improved input data (i.e.national input datasets of elevation, geology, and recharge) and adjusted resolution and computational efficiency.It is therefore called the National Water Table (NWT) model.A review of the EWT model is first given, after which the improved NWT model method is explained.Results simulated at the national scale are presented, where sensitivity to the input parameters of recharge, and hydraulic conductivity is also assessed.The NWT is furthermore evaluated at the catchment scale by comparing simulated groundwater depths to observations in two regional case studies, where the role of a simplified calibration is assessed.In addition, this paper discusses the relevance of the NWT model to solving inconsistency between regional models, and the strength and weakness of the simplifying assumptions in the model.We also discuss the importance of model input components (e.g.terrain model and hydraulic conductivity).

Review of the EWT model
The Equilibrium Water Table (EWT) model is a steady-state groundwater model that calculates water table depth (in metres below ground level: m b.g.l.) and water table elevation (in metres above sea level: m a.s.l.) at the global scale using a variety of ground-based, satellite-observed, and modelled data (Fan et al., 2013a(Fan et al., , 2007;;Fan and Miguez-Macho, 2010b;Miguez-Macho et al., 2007, 2008).The water table ("the level in the saturated zone at which the pressure is equal to the atmospheric pressure"; Heath, 1995) represents a longterm average at a broad scale without human-induced effects, e.g.pumping, draining, and irrigation.The model calculates a single water table, and therefore it generally cannot represent smaller-scale water tables and piezometric surfaces in perched or multi-layered aquifers.To calculate the water table, the long-term balance between the groundwater recharge and horizontal groundwater flow is calculated using a simple groundwater flow equation, based on a mass balance and Darcy's law (e.g.Hendriks, 2010;Dingman, 2002;Freeze and Cherry, 1979).The long-term balance is sought by a flux-based approach, i.e. by adding recharge and calculating groundwater flow, hydraulic head, and groundwater discharge to the surface for each time step unit (e.g.daily) until the model converges, i.e. until the water table no longer changes.The groundwater flow is constrained by the sea level, assuming that the hydrology is in equilibrium with the climate and sea level (Marshall and Clarke, 1999).Model input data include a global topography model, a long-term time series of global groundwater recharge, and a global soil model.Ground-observed groundwater level data have been used for calibration of hydraulic conductivity on a continental scale (Fan et al., 2013b).However, no New Zealand data were used in that calibration.As the EWT method is well described in other studies, we will not provide more detail in the main text of our study.However, the methods are described in Fan et al. (2013b) and summarised in Appendix A.
Westerhoff and White (2013) evaluated water table depths from the EWT model in New Zealand.Generally, they found that the model correctly estimated shallow water tables in New Zealand's alluvial aquifer systems, such as the Hauraki Plains, Heretaunga Plains, and the Canterbury Plains (Fig. 1).Overall, water tables have a tendency to be too shallow, and Westerhoff and White (2013) concluded that was mostly because the model does not incorporate pumping or other human draining features, e.g.tile drainage.However, they also found several issues in the Canterbury Plains, which included a bias towards calculation of shallow water table depth -i.e.EWT depths were too shallow compared to ground observations at many locations -that was possibly driven by wrong model input data (e.g.hydraulic conductivity or recharge).They also described the importance of better terrain models, as the uncertainty of the global-scale datasets is generally 10 m or more (Gesch et al., 1999;Harding et al., 1999;Rodríguez et al., 2006).Therefore, Westerhoff and White (2013) recommended improvement in the method and its input data, which included better model convergence criteria; representation of terrain; estimation of hydraulic conductivity, i.e. a better representation of the underlying geology; and rainfall recharge data that are relevant to the New Zealand climate.

Methodology
The NWT model is based on the EWT model but uses adjustments as proposed by other global-scale groundwater studies (Gleeson et al., 2011;de Graaf et al., 2015) and uses national-scale New Zealand input data.Our model improvements include a national terrain model; national recharge estimates; hydraulic conductivity estimates based on the geological map of New Zealand; and model resolution, initial conditions, and convergence criteria.

National terrain model
The Geographx New Zealand DEM 2.1 (Geographx, 2012) is a national digital elevation model with an 8 m grid resolution which combines New Zealand-based topographic data and satellite data from the Shuttle Radar Topography Mission (SRTM; USGS, 2006).More accurate elevation data are available in New Zealand, but only piecemeal at the local scale of catchments.

National rainfall recharge
Rainfall recharge to groundwater was taken from a recently developed national New Zealand rainfall recharge dataset (Westerhoff et al., 2018).This dataset contains monthly estimates of recharge from January 2000 to December 2014, including an uncertainty estimate.For the purpose of feeding recharge into the NWT model with daily time steps, mean daily recharge was calculated for each NWT model cell (Fig. 2).

Hydraulic conductivity
The NWT model uses data from the national geological map of New Zealand.The steps to derive saturated hydraulic conductivity K values suitable for NWT model input are estimation of near-surface K (Fig. 3), (optional) local calibration of near-surface K, and estimation of K over depth.These steps are described below.

Estimation of near-surface K
The NWT method applies geology data, i.e. the deeper subsurface underlying the soil, from the national 1 : 250 000 Geological Map of New Zealand (QMAP; GNS Science, 2012) to calculate near-surface K.The QMAP is a GIS-based digital map that shows "surface geology", i.e. the geology of the subsurface up to approximately 10 m depth.The polygon attributes of main rock type, secondary rock type, and age were used to estimate K (Westerhoff et al., 2018) and were partly based on a classification method described by Gleeson et al. (2011), who estimated hydraulic permeability κ (m 2 ) and its standard deviation for a range of hydrolithological classes.This method was put in the New Zealand context by Tschritter et al. (2016), who defined 10 hydrolithological classes and their associated κ (Table 1).All QMAP rock type attributes were summarised to 183 "dictionary" values, which were then assigned to permeabilities with a look-up table.Permeabilities of main rock type and secondary rock type were averaged, with main rock type weighted twice as much as secondary rock type.

Local calibration of near-surface K
For local testing purposes, a simple calibration procedure was built in.Values of the original K (Eq. 1) were adjusted (to K adj ) in iterative calibration steps in areas where water table depth and elevation were deemed to be sufficiently known from local measurements: The calibration damping factor ζ was chosen by trial and error to be 0.66 day −1 .The maximum adjustment allowed per calibration time step was 100 m day −1 .This calibration was run nine times, at intervals of 10 % of the standard model runtime (see Sect. 3.4).This procedure was only used in the case where ample ground-observed water level data were available for this study (Canterbury and Waipa; see Sects. 4.2 and 4.3,respectively).Before calibration, the ground-observed water level point data were gridded to an interpolated surface within the area containing ground observations.

Estimation of K over depth
Near-surface K was assumed to be represented by the QMAP-derived near-surface K (Fig. 3).Deeper than 10 m, an exponential decrease of hydraulic conductivity over depth was assumed, similar to Eq. (A1).As cell resolution of the NWT model is 200 m, the calibration constants used in this exponential decrease (a, b, and f min ; see Appendix, Eq.A2) were set equal to an earlier-developed 200 m resolution and locally calibrated EWT model in the Amazon Basin (Fan and Miguez-Macho, 2010a).

Model resolution, initial conditions, and convergence criteria
The model grid cell size was chosen as 200 m, in the New Zealand Transverse Mercator (NZTM) coordinate system.
Hydrol.Earth Syst.Sci., 22, 6449-6472, 2018 www.hydrol-earth-syst-sci.net/22/6449/2018/All input data were gridded at this cell size, either by averaging (if the cell size of the original input data was smaller than the model grid cell size) or by its nearest value (if cell size of the original input data was larger than the model grid cell size).Locally, the model was run at different spatial resolutions for the sake of more in-depth research on model behaviour (see Sect. 4.3).Initial estimates of water table depth depth and elevation for the NWT model were set equal to the EWT water table depth and elevation, i.e. the results of the EWT model.
The model was run with a standard model runtime of 36 525 iterations, i.e. the equivalent of 100 years in daily time steps.The EWT method had a convergence criterion, where the iterative calculation stops when the water table reaches an equilibrium, i.e. when lateral groundwater flow equals recharge for every model cell.This criterion has been changed in the NWT method; the model currently runs for no longer than 36 525 iterations.Appendix B explains that convergence tests showed that running for longer improves model convergence but does not necessarily lead to better results, because most of the changes with increased runtime occur in shallow-water features, which in reality are mostly drained through artificial drainage, i.e. by humans.Additionally, feeding in recharge that contains uncertainty causes the equilibrium to probably never be reached.
All model runs in this study were performed on one core of an Intel Core i7-6700 CPU @ 3.40 GHz with memory (RAM) of 32.0 GB.

National water table maps
National maps of NWT water table depth (Fig. 4) and water table elevation (Fig. 5) show that water table depth is relatively deep in higher mountainous regions, but water table elevation clearly follows the terrain elevation.The NWT water table depths show the locations of the main alluvial aquifers, similar to the EWT model.The NWT water table depth clearly follows the aquifer delineations of White (2001), e.g. in the Canterbury Plains and Heretaunga Plains (Fig. 6), and is detailed enough to possibly improve this delineation in some areas.The finer detail in the NWT model results in a more realistic local variation of water table, compared to the EWT model results (Fig. 1).Many very small shallowwater-table features are found in between areas where water tables are deeper.These features resemble stream valleys, as the more detailed terrain model allows the model to discharge more water to the surface in these valleys.The finer spatial detail is best appreciated when zooming in to smaller regions.Therefore, evaluation analyses are described for two regional case studies (Canterbury region, Sect.4.2; Waipa River catchment, Sect.4.3).

Evaluation of NWT water table in the Canterbury region
The Canterbury Plains is New Zealand's largest alluvial aquifer system.Nationwide, regional groundwater allocation is largest in the Canterbury region (Rajanayaka et al., 2010), with the majority of groundwater use in the plains.Surface geology in the Canterbury region is dominated by metamorphic rock and floodplains (Fig. 7).Metamorphic rock includes a large area of greywacke that forms much of the Southern Alps.Aquifers in the higher plains are mostly unconfined; in the lower plains they are typically confined and supply groundwater to extensive networks of spring-fed streams around Lake Ellesmere and the city of Christchurch.These confining conditions are due, in part, to the occurrence of marine sediments deposited near the current coast during interglacial periods.However, the pattern of groundwater flow is not influenced by any extensive, low-permeability layers of sediment within the upper 100 to 150 m below the water table (Hanson and Abraham, 2009).Therefore, application of the NWT (and EWT) methods, which assume unconfined aquifers, is justified in the Canterbury Plains.
Groundwater level observations from 3286 wells in Canterbury, acquired from 1894 to 2013 (Kelly Palmer, personal communication, 2013), have been used in this research.Most of these wells are located in the Canterbury Plains.It is assumed that these observations represent the water table and can thus be compared to NWT water table depth and water table elevation.All groundwater level observations were carefully quality-checked (as described in Westerhoff and White, 2013).NWT water table depth and elevation covering the locations of these wells were sampled from both the NWT and EWT datasets.Median observed groundwater level for all wells was calculated.From here on, these will be referred to as "ground-observed water level".The national Geographx elevation model was used for conversion of groundwater depth to groundwater elevation.
NWT water table depths are similar to ground-observed water depth: for example, 24 % and 53 % of NWT water table depths are within 1 and 3 m from ground-observed water depth, respectively (Table 2).This is an improvement to the EWT water table depths, where this was 13 % and 42 %, respectively.Also, the NWT model shows significantly fewer large discrepancies than the EWT model.For example, 3.4 %, 0.18 %, and 0.03 % of the NWT water table depths differed by more than 50, 100, and 150 m, respectively, from the ground-observed water depth.For the EWT water table depths, these statistics were 6 %, 1.2 %, and 0.05 %.
Correlation of NWT water table elevation with groundobserved water level (m a.s.l.) is high (R = 0.99, Fig. 8, right) but does show local differences, which are shown more clearly in the crossplot of water table depth (Fig. 8, left).Westerhoff and White (2013) showed that evaluation of the EWT model revealed higher discrepancies in the same area (i.e. the area between the Ashburton and Rakaia rivers, Figs. 9 and 10).Although the discrepancies are smaller for the NWT model than for the EWT model in these areas, they occur at the same locations.Possible explanations for these discrepancies are considered in Sect. 5.

Evaluation of NWT water table in the Waipa River catchment
The Waipa River is a tributary of the Waikato River, in the Waikato region of New Zealand (Fig. 11), with a catchment size of approximately 300 000 ha. Increasing agricultural land use and deforestation, mainly in the low-lying Hamilton Basin, could potentially result in the deterioration of water quality in the Waipa River catchment, according to Rawlinson (2014), who performed a review on existing information in the catchment.They calculated a map of water table elevation using observed groundwater elevation from 758 wells located in the catchment.However, these observations are sparse because, temporally, data mostly consist of few measurements and, spatially, data are concentrated more in some areas than others; i.e. most wells are located in the Hamilton Basin.The NWT model can thus provide a better spatial insight into the water table in the catchment, provided that they correlate well enough with the ground observations.NWT water table depths clearly demonstrate the location of shallow water tables in the low-lying basin area (Fig. 13, left).The spatial pattern of NWT water table elevation is a good visual match to the interpolated ground-observed water table elevation (Fig. 12).In addition, the NWT map shows more detail than the interpolated ground-observed surface and calculates water table depth at many places where there are no ground observations.Because of its higher resolution, the NWT also shows more detail than the much coarser EWT water table (Fig. 13).
Because of the availability of a lidar terrain model for this study, statistics were collected for a multitude of model runs using different terrain models with different resolutions (100 m, 200 m, 500 m, 1 km).The effect of calibration towards the ground observations was also taken into account.Table 3 shows that correlation with ground observations improves with higher model resolution.Calibration towards hydraulic conductivity, albeit a simple implementation, improves the correlation of the modelled water table with ground-observed water level model output.Inclusion of the lidar DEM does not significantly improve the results of the uncalibrated run, but it does improve calibrated model runs.For the best model run (lidar DEM, 100 m resolution, calibrated), correlation of model results to ground-observed data was R = 0.41 and R = 0.97 for water table depth and water table elevation, respectively (Fig. 14).
Correlation of both runs is much higher than that of the EWT model (R = 0.1 and R = 0.91, Fig. 15).Water table depth shows significantly lower correlation than water table elevation, as was also the case in the evaluation in the Canterbury region.The Discussion section describes possible causes for these differences.

Sensitivity of the water table depth to recharge and hydraulic conductivity
A total of 35 (uncalibrated) model runs for the NWT model at the national scale were performed, where the value of recharge and K model input were weighted.Weighting values for recharge (W r ) were in between 0.8 and 1.2 (be-    cause the nominal uncertainty of recharge of New Zealand's aquifers is less than 20 %; Westerhoff et al., 2018).Weighting values for K (W K ) were chosen to follow a logarithmic distribution, i.e. in between 0.1 (10 %) and 10 (1000 %).Results are shown in Fig. 16.The average water table is sensitive both to recharge (an increase in recharge results in an increase in average water table depth) and to K (a decrease in K results in a increase in average water table depth).Spatial plots of two runs (Fig. 17, all 35 model results shown in the Supplement) show that most changes are taking place in the the foothills of aquifer areas.

Discussion
The NWT model estimates water The NWT model includes all catchments of the mainland of New Zealand.Because water is primarily regulated at the regional level, regional models can show different results (e.g.groundwater catchment delineation) at regional boundaries.These inconsistencies can lead to trans-boundary issues, such as catchment boundary definition and source protection.A fundamental role of groundwater science is to identify and characterise groundwater catchments of water supplies.Better delineation of groundwater catchments is therefore essential for source protection.Source protection is clearly an important issue associated with the prevention of waterborne diseases, and a national environmental standard for source protection, including groundwater, has been proposed for New Zealand (Ministry for the Environment, 2009Environment, , 2016)).Where two regional groundwater models show different results, the NWT model can be used as a "third" model to help solve those inconsistencies.This makes the NWT model relevant to solving trans-boundary issues and to assisting in provision of nationally consistent groundwater data.
The advantage of the NWT model, compared to catchment-scale groundwater flow models, is that its simple algorithm makes computation of the water table across all catchments relatively fast.The NWT model can therefore provide useful preliminary water table information for catchment-scale numerical groundwater flow models, in both data-sparse and data-rich regions.Furthermore, the NWT model inputs can provide other initial estimates that are useful to other models.For example, data of hydraulic conductivity and recharge are also nationwide.The disadvantage of simplifying assumptions in the NWT model is that more complex groundwater features are not handled well.Currently, the model does not include confining layers reliably, nor does it incorporate fractures and groundwater age.This is mostly due to the simplifying assumption of K over depth.The assumption of a decrease of K with depth is based on unconfined water movement.Water table depth therefore does not necessarily equal groundwater level in confined aquifers, as the confining layer holds the water under higher pressure at greater depths than the modelled water table depth.The model calculates one water level in this case, i.e. the unconfined groundwater level.In the case of artesian water (i.e. a water table elevation higher than the ground surface) the NWT discharges all water above the ground surface.NWT water table elevation can thus not show artesian water table elevation as the water table is set at the ground surface.Finally, the calibration module of K (see Sect. 3.3.2) was shown to yield better groundwater table elevations in the case studies of Canterbury and Waipa.It was chosen not to further refine this coarse calibration (e.g. by testing K calibration in the log domain or by calibrating more often than only once per 10 years).That is mainly because, in such cases, more advanced local groundwater flow models, including better calibration options, might be preferable.These models could then still be constrained with initial estimates of NWT water tables, hydraulic conductivities, and recharge in data-sparse areas.
As with all groundwater models, the NWT and EWT models are sensitive to input parameters of terrain, recharge and hydraulic conductivity.We have shown that model correlation with ground-observed information improves with model resolution (because rivers are better resolved; see earlier in the Discussion section).Additionally, we have shown that a better-quality terrain model improves results, but only if the model is calibrated.Finally, sensitivity analyses in this study  have quantified the sensitivity of the water table to recharge and K (i.e. the water table increases with increasing recharge and with decreasing K), where sensitivity is highest in the foothills and more elevated aquifer areas.Further study is recommended on using the results of our model tests to further develop a NWT model calibration which includes these findings and is better suited to multiple model input parameters.However, in doing that, we should be careful to avoid too many time-consuming (e.g.Monte Carlo and calibration) runs.After all, we are after a model that remains relatively simple, i.e. bridges the gap between global-scale and localscale models.Furthermore, calibration and validation studies with a model that is sensitive in the foothills might not neces-sarily show improved validation, because validation data in the foothills are scarce.More in-depth calibration and validation studies are therefore recommended that use the findings of this study; these might then also best be done by more advanced local numerical flow models.
Both the NWT and the EWT have a simplified representation of gaining and losing reaches.First, we address the simplifying assumptions about gaining reaches.The EWT and NWT models have a hydraulic-gradient-dependent interaction: if the water reaches the surface, it flows out of the system.That pattern follows the elevation, and thus it follows the rivers.Hence, through the model one will easily spot the rivers when looking at this outflow, showing the groundwater  table equal to the surface elevation.The underlying assumption in the model is that rivers are resolved in the grid used for integration.This also explains why calculations done at high resolution lead to better results and have better correlation with ground-observations.The NWT model thus is an improvement of the EWT model, since its resolution is higher than the original global 30 arcsec grid.However, at about 200 m resolution rivers are still not resolved well enough, and we recommend running the model with even higher resolution.However, before doing that, we need an optimised set of model conditions, such as the balance between required (smaller) model time steps and computational effort.
Bearing in mind the strength of fast and simplistic models, we should then also consider that this might better be done by more advanced local numerical models.Second, we address the simplifying assumptions about losing reaches.Losing rivers are only simulated by the model assuming that this is recharge; i.e. additional recharge from river runoff is not taken into account.If data were available on where rivers are losing, this could be implemented.However, New Zealand has no comprehensive nationwide datasets of losing (or gaining) reaches.Many advanced numerical groundwater models also still have difficulty in implementing gaining and losing river reaches, because of the numerical complexity and com-Hydrol.Earth Syst.Sci., 22, 6449-6472, 2018 www.hydrol-earth-syst-sci.net/22/6449/2018/ putational effort required to do so and the lack of information on losing and gaining reaches.
Water table depths showed lower correlation with ground observations than water table elevation in the Canterbury and Waipa examples.That is mainly because small inaccuracies of water table depth are much more significant than water table elevations at shallow water tables.For example, in an area that lies 100 m a.s.l., the water table depth error could be highly significant (e.g. 2 ± 1 m b.g.l.), where the difference in water table elevation is much less significant (e.g.98 ± 1 m a.s.l.).The most likely candidate for the cause of such inaccuracies is the uncertainty of the terrain model.For example, Rodríguez et al. (2006) concluded that the SRTM data, on which the national terrain model is partly based, can have average absolute height errors of more than 10 m.Comparison of the lidar terrain model in the Waipa evaluation with the national terrain model showed a median absolute difference of 9 m (for areas below 20 m a.s.l., the difference was 4 m).Correlation of water table elevation data is less dependent on the inaccuracy of the terrain models used.However, correlation plots of water table depth proved useful in this study, as they monitored the improvements of model runs by fine-tuning input parametrisation.For example, in the Waipa River catchment correlation of water table depth with ground observations significantly improved when the lidar terrain model was used (Table 3).In both the Canterbury region and Waipa River catchment we used water table depth correlation with test improvement of the calibration for K.
The abovementioned lower correlation with water table depth is amplified by the fact that the model still does not incorporate human pumping or draining; i.e. it still has a bias towards shallow water tables (although less than the EWT due to improvement in model resolution and K).With very high resolution input data of terrain, which would include these small but important drainage features, this issue would be better resolved, as groundwater discharge would be better resolved.However, this bias of the shallow water table is also a correct indicator of the fact that many of these shallowwater-table areas used to be wetlands: approximately 90 % of wetlands have been lost since European settlement in New Zealand (StatsNZ and Ministry for the Environment, 2018), mostly due to drainage for agriculture.
This research confirms that high-quality model input data, such as hydraulic properties and accurate terrain models, are shown to be important to the improvement of groundwater models.More accurate input data than used in the NWT model exist in New Zealand.For example, high-resolution mapping of near-surface geology has been achieved at the regional scale (Westerhoff et al., 2014) through Airborne Electromagnetic (AEM) geophysics, and some regional terrain models are based on lidar data (e.g.Waikato Regional Council, 2016).However, data acquisition and management of those geophysical and terrain models are at the regional scale.There is no comprehensive national-scale AEM-and lidar-based information as of yet.This study therefore rec- ommends continuing efforts of characterisation of these model input data at the national scale.
Improvements of the NWT model can lead to improved insights of the global-scale model approach.For example, the EWT model crucially requires geology data to infer better water table estimates.The approach in this study, i.e. using hydraulic permeability based on a method of Gleeson et al. (2011), was also used by de Graaf et al. (2015) for the global scale.We recommend that the EWT model approach also embed these improved input data and merge other nationalscale geological data from other countries.Furthermore, we recommend that the findings of our sensitivity analyses and our calibration of hydraulic conductivity at locations where water tables are well known be used to further improve the global-scale EWT model.

Conclusions
The NWT model is a revised version of the global-scale EWT model that has been improved for application in New Zealand.The NWT model gives an estimate of water table depth and water table elevation with a 200 m grid resolution.The NWT model uses slightly adjusted model parametrisation, but mostly improved input datasets, amongst which are a national terrain model, a national digital geological dataset, and a nationwide recharge dataset.
Because of the improvements, NWT water table depths show the areas of alluvial aquifers, including their spatial variation of water table depths, better and with higher spatial resolution than the EWT model.The NWT water table elevation shows excellent correlation with ground-observed water level data in the Canterbury region and Waipa River catchment.The NWT model estimates the water table at places where there are no ground observations and therefore shows more detail than other existing interpolated water table surfaces.In terms of correlation and spatial resolution, the NWT model outperforms the EWT model.
Because of its simplified character, the NWT model has the advantage of fast calculation at the national scale.In fact, it is currently the only nationwide groundwater model dedicated to New Zealand application.The NWT water table, as well as its nationwide data components (e.g.recharge and hydraulic conductivity), can be used as an initial estimate for more advanced catchment-scale numerical groundwater flow models where data are sparse.In addition, the NWT model might also be used to solve inconsistency of different regional models at regional boundaries.
Sensitivity analyses at the national scale show that the NWT and EWT models are most sensitive to recharge and hydraulic conductivity in the foothills of aquifers.Calibration tests show that improved model resolution leads to better results (i.e. higher correlation with ground observations) and that a better-quality terrain model (i.e.lidar used in the Waipa local case study) improves model results, but only if the model is calibrated.
Use of the NWT model parametrisation improvements could lead to the improvement of the global-scale EWT model, for example in a better estimation method of hy-draulic conductivity.Also, the findings of our model tests and sensitivity analyses might be extrapolated to the global EWT application.We therefore recommend the findings of this study to be used for improvement of the global-scale EWT model application.
The NWT model does not handle complex groundwater features well -i.e.confining layers, fractures, and groundwater age -because the model contains simplifying assumptions.Because we want the NWT model to remain a simplified model, i.e. the model that bridges the gap between (too-)expensive advanced local models and (too-)simple globalscale models (see Sect. 1), we recommend that state-ofthe-art numerical catchment-scale groundwater flow models should be used in those circumstances if available.However, the NWT model is still useful to provide initial model estimates (i.e.water table, hydraulic conductivity, and recharge) to those more advanced models.
The NWT model still has a bias towards shallow water tables, although less than the EWT model because of the finer model resolution.This study shows that improved input data of hydraulic conductivity (or calibration towards it) further reduces this bias.However, this bias of shallow water table is also a valuable indicator that correctly shows that many of the indicated shallow-water-table areas used to be wetlands: approximately 90 % of wetlands have been drained since European settlement in New Zealand, mostly to develop agriculture.
Possible improvements of the NWT model are the use of better model input components, such as a better terrain model, and improved model calibration.Therefore, this study recommends further efforts in making available high-quality nationwide geophysical, terrain, and water level data at the national scale of New Zealand.Hydrol.Earth Syst.Sci., 22, 6449-6472, 2018 www.hydrol-earth-syst-sci.net/22/6449/2018/ Supplement.The supplement related to this article is available online at: https://doi.org/10.5194/hess-22-6449-2018-supplement.
Author contributions.RW conceived, designed, and conducted the experiments and analysed all data.RW and PW wrote the paper.PW advised in the validation experiments in Canterbrury and Waipa.GM-M wrote parts of the discussion texts where the NWT was compared to the EWT.PW and GM-M furthermore advised in comprehensively addressing all reviewer comments.
Competing interests.The authors declare that they have no conflict of interest.
Special issue statement.This article is part of the special issue "Integration of Earth observations and models for global water resource assessment".It is not associated with a conference.

Figure 1 .
Figure 1.EWT water table depth in New Zealand, after Fan et al. (2013a).Main cities (black text boxes) as well as main alluvial aquifer are pointed out (grey text boxes).

R
. Westerhoff et al.: Application of an improved global-scale groundwater model for water table estimation

Figure 3 .
Figure 3. Hydraulic conductivity estimates for near-surface geology across New Zealand.

Figure 4 .
Figure 4. NWT water table depth in New Zealand.

Figure 5 .
Figure 5. NWT water table elevation in New Zealand.

Figure 8 .
Figure 8. Correlation of NWT water table depth and water table elevation with ground observations in the Canterbury region.The dashed line is the 1 : 1 relation.

Figure 9 .
Figure 9. Water table depth of the (a) NWT and (b) EWT models in the Canterbury region.

Figure 10 .
Figure 10.Correlation of EWT water table depth and water table elevation with ground observations in the Canterbury region.The dashed line is the 1 : 1 relation.

Figure 11 .
Figure 11.Waipa River catchment, Waikato region, New Zealand.The study area is shown in red.The low-lying Hamilton Basin is shown in black (Rawlinson, 2014, their Fig.1.2).

Figure 12 .
Figure 12.(a) NWT water table elevation and (b) the interpolated surface of water table elevation from ground observation in the Waipa catchment (Rawlinson, 2014).All elevations are in m a.s.l.Ground observations used for the interpolated surface are shown in black dots.

Figure 13 .
Figure 13.Water table depth of the (a) NWT and (b) EWT models in the Waipa River catchment.

Figure 14 .
Figure 14.Correlation of NWT water table depth (100 m lidar DEM, calibrated to ground observations) and water table elevation with ground observations in the Waipa River catchment.The dashed line is the 1 : 1 relation.

Figure 15 .
Figure 15.Correlation of EWT water table depth and water table elevation with ground observations in the Waipa River catchment.The dashed line is the 1 : 1 relation.

Figure 16 .
Figure 16.Mean water table difference plot for 35 model runs, with varying weights for recharge (W r ) and K (W K ).Difference is compared to the reference run of W r = 1 and W K = 1.

Figure 17 .
Figure 17.Output of two models runs with different weights for recharge (W r ) and K (W K ).

Figure A1 .
Figure A1.Schematic of the EWT model to simulate the water table at the continental scale, using recharge (R), horizontal groundwater flow (Q), groundwater discharge in rivers (Q r ), and sea level (boundary condition).The blue fading colours indicate the decrease of hydraulic conductivity with depth (Fan and Miguez-Macho, 2010b, their Fig. 4a).

Figure B1 .
Figure B1.Model area of the Mataura catchment, Southland, New Zealand.

Figure B2 .
Figure B2.Convergence test in the Mataura catchment, Southland, New Zealand, with inclusion of uncertainty and a Gaussian distribution over the year to include seasonality of recharge.

Figure B3 .
Figure B3.Convergence test in the Mataura catchment, Southland, New Zealand, with inclusion of uncertainty and a Gaussian distribution over the year to include seasonality of recharge.

Table 2 .
Absolute differences between ground-observed and modelled water table depths ( ) and their percentage of occurrence in the Canterbury region for the NWT and EWT model.

Table 3 .
Statistics on different model runs in the Waipa River catchment.RMSE stands for root mean square error.