Rainfall is spatially and temporally varying, reflecting spatial differences in geology and temporal variation in rainfall. The derivation of a mean annual recharge surface for the using a chloride mass balance approach is described in Section 2.1.3 of companion product 2.1-2.2 for the Hunter subregion (). The long-term temporal variation of rainfall is described in Section 1.1.2 of the Hunter context statement (companion product 1.1 ()). This is normalised so its average throughout the period 1983 to 2012 is one, and the resultant time series is multiplied by the spatial variation to yield the final recharge applied to the model. Recharge is applied as a source of water of prescribed rate to the land surface of the model. To account for in both the temporal and spatial variation of recharge, its magnitude is varied in the uncertainty analysis (see companion submethodology M07 (as listed in Table 1) for modelling ()).
Evapotranspiration is represented by a sink of applied across the entire land surface of the model, excluding the parts under the ocean. Generally in groundwater models, evapotranspiration from groundwater is assumed to be maximum when the phreatic surface is close to the ground surface, or above it (e.g. in the case of ponding). This gives rise to the so-called ‘potential evapotranspiration’, PET, which is the maximum possible rate of evapotranspiration. PET is temporally and spatially varying, PET = PET(x,y,t), which reflects seasonal changes, and different vegetation and surface expression of geology. Conversely, evapotranspiration from groundwater is generally assumed to be zero when the phreatic surface is deep below the ground surface since plant roots cannot draw water from the deep groundwater reserves. This gives rise to the so-called ‘extinction depth’, d, which parameterises the depth of plant roots. Clearly d depends on the type of vegetation, as discussed below. When the phreatic surface is somewhere between the extinction depth and the ground surface, evapotranspiration is a fraction of PET.
The Hunter groundwater model uses the daily PET time series produced by the Australian Water Resource Assessment landscape (AWRA-L) model (Bioregional Assessment Programme, ), which was based on one particular future climate scenario (See Section 22.214.171.124 in companion product 2.6.1 for the Hunter subregion (Zhang et al., 2018)). The groundwater model is run using monthly time steps. Appropriate averaging of the daily PET time series provides data that can be used in the groundwater model, and this is conveniently represented by:
Generally in groundwater models the extinction depth, d, is set between 0 m (zero vegetation) and 20 m (deep rooted trees) with a value of 5 m being typical (). In the Hunter groundwater model d varies between 0 m and 10 m. It is assumed that d is proportional to vegetation height, V, with the rationale being that taller trees typically have deeper roots than smaller shrubs or grasses (). Vegetation height is quantified in Figure 18 of companion product 1.1 for the () and details of the vegetation height (Caltech/JPL, ) can be found in Section 2.1.1 of companion product 2.1-2.2 for the Hunter subregion (). Since it varies between 0 m and 39 m throughout the subregion, the groundwater model uses:
At any point on the surface of the Hunter groundwater model, evapotranspiration varies as a cubic function of groundwater head, h, at that point, as shown in Figure 15. When h is high (the groundwater level is above the topography) ET = PET, but when h is low (the groundwater level is well below the topography) ET = 0. The transition from the cubic to PET is chosen to be at h = 2 m. If this transition was chosen to be at h = 0, as is common, and a point happened to have d = 0 (which is quite common when using the vegetation height to define d), the numerical solver may only slowly converge to the solution for groundwater levels due to the sudden transition from PET to zero.
ET = PET for h>2 m as indicated by the dashed lines.
Product Finalisation date
- 126.96.36.199 Methods
- 188.8.131.52 Review of existing models
- 184.108.40.206 Model development
- 220.127.116.11 Boundary and initial conditions
- 18.104.22.168 Implementation of the coal resource development pathway
- 22.214.171.124 Parameterisation
- 126.96.36.199 Observations and predictions
- 188.8.131.52 Uncertainty analysis
- 184.108.40.206 Limitations and conclusions
- Currency of scientific results
- Contributors to the Technical Programme
- About this technical product