Point-scale multi-objective calibration of the Community Land Model (version 5.0) using in situ observations of water and energy fluxes and variables

This study evaluates water and energy fluxes and variables in combination with parameter optimization of version 5 of the state-of-the-art Community Land Model (CLM5) land surface model, using 6 years of hourly observations of latent heat flux, sensible heat flux, groundwater recharge, soil moisture and soil temperature from an agricultural observatory in Denmark. The results show that multi-objective calibration in combination with truncated singular value decomposition and Tikhonov regularization is a powerful method to improve the current practice of using lookup tables to define parameter values in land surface models. Using measurements of turbulent fluxes as the target variable, parameter optimization is capable of matching simulations and observations of latent heat, especially during the summer period, whereas simulated sensible heat is clearly biased. Of the 30 parameters considered, the soil texture, monthly leaf area index (LAI) in summer, stomatal conductance and root distribution have the highest influence on the local-scale simulation results. The results from this study contribute to improvements of the model characterization of water and energy fluxes. This work highlights the importance of performing parameter calibration using observations of hydrologic and energy fluxes and variables to obtain the optimal parameter values for a land surface model.</p

Transboundary geophysical mapping of geological elements and salinity distribution critical for the assessment of future sea water intrusion in response to sea level rise

Geophysical techniques are increasingly being used as tools for characterising the subsurface, and they are generally required to develop subsurface models that properly delineate the distribution of aquifers and aquitards, salt/freshwater interfaces, and geological structures that affect groundwater flow. In a study area covering 730 km2 across the border between Germany and Denmark, a combination of an airborne electromagnetic survey (performed with the SkyTEM system), a high-resolution seismic survey and borehole logging has been used in an integrated mapping of important geological, physical and chemical features of the subsurface. The spacing between flight lines is 200–250 m which gives a total of about 3200 line km. About 38 km of seismic lines have been collected. Faults bordering a graben structure, buried tunnel valleys, glaciotectonic thrust complexes, marine clay units, and sand aquifers are all examples of geological structures mapped by the geophysical data that control groundwater flow and to some extent hydrochemistry. Additionally, the data provide an excellent picture of the salinity distribution in the area and thus provide important information on the salt/freshwater boundary and the chemical status of groundwater. Although the westernmost part of the study area along the North Sea coast is saturated with saline water and the TEM data therefore are strongly influenced by the increased electrical conductivity there, buried valleys and other geological elements are still revealed. The mapped salinity distribution indicates preferential flow paths through and along specific geological structures within the area. The effects of a future sea level rise on the groundwater system and groundwater chemistry are discussed with special emphasis on the importance of knowing the existence, distribution and geometry of the mapped geological elements, and their control on the groundwater salinity distribution is assessed

Impact of urban geology on model simulations of shallow groundwater levels and flow paths

This study examines the impact of urban geology and spatial discretization on the simulation of shallow groundwater levels and flow paths at the city scale. The study uses an integrated hydrological model based on the MIKE SHE code that couples surface water and 3D groundwater simulations with a leaky sewer system. The effect of the geological configuration was analyzed by applying three geological models to an otherwise identical hydrological model. The effect of spatial discretization was examined by using two different horizontal discretizations for the hydrological models of 50 and 10 m, respectively. The impact of the geological configuration and spatial discretization was analyzed based on model calibration, simulations of high water levels, and particle tracking. The results show that a representation of the subsurface infrastructure, and near-terrain soil types, in the geological model impacts the simulation of the high water levels when the hydrological model is simulated in a 10 m discretization. This was detectable even though the difference between the geological models only occurs in 7 % of the volume of the geological models. When the hydrological model was run in a 50 m horizontal discretization, the impact of the urban geology on the high water levels was smoothed out. Results from particle tracking show that representing the subsurface infrastructure in the hydrological model changed the particles' flow paths and travel time to sinks in both the 50 and 10 m horizontal discretization of the hydrological model. It caused less recharge to deeper aquifers and increased the percentage of particles flowing to saturated-zone drains and leaky sewer pipes. In conclusion, the results indicate that even though the subsurface infrastructure and fill material only occupy a small fraction of the shallow geology, it affects the simulation of local water levels and substantially alters the flow paths. The comparison of the spatial discretization demonstrates that, to simulate this effect, the spatial discretization needs to be of a scale that represents the local variability in the shallow urban geology.</p

Transboundary geophysical mapping of geological elements and salinity distribution critical for the assessment of future sea water intrusion in response to sea level rise

Geophysical techniques are increasingly being used as tools for characterising the subsurface, and they are generally required to develop subsurface models that properly delineate the distribution of aquifers and aquitards, salt/freshwater interfaces, and geological structures that affect groundwater flow. In a study area covering 730 km&lt;sup&gt;2&lt;/sup&gt; across the border between Germany and Denmark, a combination of an airborne electromagnetic survey (performed with the SkyTEM system), a high-resolution seismic survey and borehole logging has been used in an integrated mapping of important geological, physical and chemical features of the subsurface. The spacing between flight lines is 200–250 m which gives a total of about 3200 line km. About 38 km of seismic lines have been collected. Faults bordering a graben structure, buried tunnel valleys, glaciotectonic thrust complexes, marine clay units, and sand aquifers are all examples of geological structures mapped by the geophysical data that control groundwater flow and to some extent hydrochemistry. Additionally, the data provide an excellent picture of the salinity distribution in the area and thus provide important information on the salt/freshwater boundary and the chemical status of groundwater. Although the westernmost part of the study area along the North Sea coast is saturated with saline water and the TEM data therefore are strongly influenced by the increased electrical conductivity there, buried valleys and other geological elements are still revealed. The mapped salinity distribution indicates preferential flow paths through and along specific geological structures within the area. The effects of a future sea level rise on the groundwater system and groundwater chemistry are discussed with special emphasis on the importance of knowing the existence, distribution and geometry of the mapped geological elements, and their control on the groundwater salinity distribution is assessed

Climate model uncertainty versus conceptual geological uncertainty in hydrological modeling

Projections of climate change impact are associated with a cascade of uncertainties including in CO₂ emission scenarios, climate models, downscaling and impact models. The relative importance of the individual uncertainty sources is expected to depend on several factors including the quantity that is projected. In the present study the impacts of climate model uncertainty and geological model uncertainty on hydraulic head, stream flow, travel time and capture zones are evaluated. Six versions of a physically based and distributed hydrological model, each containing a unique interpretation of the geological structure of the model area, are forced by 11 climate model projections. Each projection of future climate is a result of a GCM–RCM model combination (from the ENSEMBLES project) forced by the same CO₂ scenario (A1B). The changes from the reference period (1991–2010) to the future period (2081–2100) in projected hydrological variables are evaluated and the effects of geological model and climate model uncertainties are quantified. The results show that uncertainty propagation is context-dependent. While the geological conceptualization is the dominating uncertainty source for projection of travel time and capture zones, the uncertainty due to the climate models is more important for groundwater hydraulic heads and stream flow