Evolution of the variability of surface temperature and vegetation density in the great plains

This study focuses on how the variability of land surface temperature and vegetation density at the SGP ARM-CART site changes over episodic (day to day) and seasonal time scales using AVHRR satellite data. Four drying periods throughout the year are analyzed. Land surface temperature had an erratic relationship with time exhibiting no deterministic pattern from day-to-day or season-to-season. Furthermore, it did not exhibit spatial pattern persistence. On the other hand, vegetation density had a consistent spatial pattern and temporal decay during average length drying periods (less than 7 days) as well as within a season. However, there were distinct differences in the seasonal pattern of variation between winter and growing seasons. In addition, the paper highlights a methodology to quantify the relationships that exist at the land surface between the primary parameter of interest and the controlling variables.     

Evaluating the dual‐boundary forcing concept in subsurface–land surface interactions of the hydrological cycle

Subsurface and land surface processes (e.g. groundwater flow, evapotranspiration) of the hydrological cycle are connected via complex feedback mechanisms, which are difficult to analyze and quantify. In this study, the dual‐boundary forcing concept that reveals space–time coherence between groundwater dynamics and land surface processes is evaluated. The underlying hypothesis is that a simplified representation of groundwater dynamics may alter the variability of land surface processes, which may eventually affect the prognostic capability of a numerical model. A coupled subsurface–land surface model ParFlow.CLM is applied over the Rur catchment, Germany, and the mass and energy fluxes of the coupled water and energy cycles are simulated over three consecutive years considering three different lower boundary conditions (dynamic, constant, and free‐drainage) based on groundwater dynamics to substantiate the aforementioned hypothesis. Continuous wavelet transform technique is applied to analyze scale‐dependent variability of the simulated mass and energy fluxes. The results show clear differences in temporal variability of latent heat flux simulated by the model configurations with different lower boundary conditions at monthly to multi‐month time scales (~32–91 days) especially under soil moisture limited conditions. The results also suggest that temporal variability of latent heat flux is affected at even smaller time scales (~1–3 days) if a simple gravity drainage lower boundary condition is considered in the coupled model. This study demonstrates the importance of a physically consistent representation of groundwater dynamics in a numerical model, which may be important to consider in local weather prediction models and water resources assessments, e.g. drought prediction. Copyright © 2015 John Wiley & Sons, Ltd.     

Oceanographic and meteorological forcing of the pelagic ecosystem on the Gulf of Cadiz shelf (SW Iberian Peninsula)

The spatial and temporal distribution of physical, chemical and biological variables of the NE continental shelf of the Gulf of Cadiz were analyzed monthly during almost three annual cycles. This analysis was performed with the aim of deriving the main forcing factors controlling variability at inter-annual, seasonal and short-time scales. Meteorological forcing related to heavy episodes of rainfall that affected river discharges and the wind regime, controlled both the currents along the shelf together and the nutrient concentrations of the surface waters. Meteorological forcing in turn determined the subsequent development and maintenance of phytoplankton blooms. Superimposed on the seasonal cycle typical of temperate latitudes, the inputs of continental nutrients mainly from the Guadalquivir River, along with episodes of upwelling favored by the predominance of westerly winds triggered phytoplankton growth on the shelf, highlighting the markedly relevant role of this large estuary in the control of the biological activity on the shelf.     

Spatio‐temporal variability of water and energy fluxes – a case study for a mesoscale catchment in pre‐alpine environment

Water and energy fluxes at and between the land surface, the subsurface and the atmosphere are inextricably linked over all spatio‐temporal scales. Our research focuses on the joint analysis of both water and energy fluxes in a pre‐alpine catchment (55 km²) in southern Germany, which is part of the Terrestrial Environmental Observatories (TERENO). We use a novel three‐dimensional, physically based and distributed modelling approach to reproduce both observed streamflow as an integral measure for water fluxes and heat flux and soil temperature measurements at an observation location over a period of 2 years. While heat fluxes are often used for comparison of the simulations of one‐dimensional land surface models, they are rarely used for additional validation of physically based and distributed hydrological modelling approaches. The spatio‐temporal variability of the water and energy balance components and their partitioning for dominant land use types of the study region are investigated. The model shows good performance for simulating daily streamflow (Nash–Sutcliffe efficiency > 0.75). Albeit only streamflow measurements are used for calibration, the simulations of hourly heat fluxes and soil temperatures at the observation site also show a good performance, particularly during summer. A limitation of the model is the simulation of temperature‐driven heat fluxes during winter, when the soil is covered by snow. An analysis of the simulated spatial fields reveals heat flux patterns that reflect the distribution of the land use and soil types of the catchment. The water and energy partitioning is characterized by a strong seasonal cycle and shows clear differences between the selected land use types. Copyright © 2016 The Authors Hydrological Processes Published by John Wiley & Sons Ltd.     

VHF radar observations of surface currents off the northern Opal coast in the eastern English Channel

Two very high-frequency radars (VHFR) operating on the Opal coast of eastern English Channel provided a nearly continuous 35-day long dataset of surface currents over a 500 km² area at 0.6–1.8 km resolution. Argo drifter tracking and CTD soundings complemented the VHFR observations, which extended approximately 25 km offshore. The radar data resolve three basic modes of the surface velocity variation in the area, that are driven by tides, winds and freshwater fluxes associated with seasonal river discharge. The first mode, accounting for 90% of variability, is characterized by an along-shore flow pattern, whereas the second and third modes exhibit cross-shore, and eddy-like structures in the current velocity field. All the three modes show the dominant semi-diurnal variability and low-frequency modulation by the neap-spring tidal cycle. Although tidal forcing provides the major contribution to variability of local currents, baroclinicity plays an important role in shaping the 3D velocity field averaged over the tidal cycle and may strongly affect tracer dynamics on larger time scales. An empirical orthogonal function (EOF) decomposition and a spectral rotary analysis of the VHFR data reveal a discontinuity in the velocity field occurring approximately 10 km offshore which was caused by the reversal in the sign of rotation of the current vector. This feature of local circulation is responsible for surface current convergence on ebb, divergence on flood and strong oscillatory vertical motion. Spectral analysis of the observed currents and the results of the Agro drifter tracking indicate that the line of convergence approximately follows the 30-m isobath. The most pronounced feature of the radar-derived residual circulation is the along-coast intensification of surface currents with velocity magnitude of 0.25 m/s typical for the Regions of Freshwater Influence (ROFI). The analysis has provided a useful, exploratory examination of surface currents, suggesting that the circulation off the Opal coast is governed by ROFI dynamics on the hypertidal background.     

Hydrographic variability on a coastal shelf directly influenced by estuarine outflow

Hydrographic variability on the Alabama shelf just outside of Mobile Bay, a major source of river discharge into the Gulf of Mexico, is examined using time series of water column temperature and surface and bottom salinity from a mooring site with a depth of 20 m in conjunction with a series of across-shelf CTD surveys. The time series data show variability in a range of time scales. The density variation is affected by both salinity and temperature, with its relatively strong annual signal mostly determined by temperature and its year to year variability mostly determined by salinity. Seasonal mean structures of temperature, salinity, and density show a transition from estuarine to shelf conditions in which three regions with distinct seasonal characteristics in their horizontal and vertical gradient structures are identified. Correlation analysis with the available forcing functions demonstrates the influence of Mobile Bay on the variability at the mooring site. At low frequencies, river discharge from Mobile Bay has a varying influence on salinity, which is absent during the periods with unusually low discharge. At shorter synoptic time scales, both the estuarine response to the across-shelf wind stress and the shelf response to the along-shelf wind stress are significantly correlated with temperature/salinity variability: the former becoming important for the surface layer during winter whereas the latter for the bottom layer during both winter and summer. These forcing functions are important players in determining the estuarine-shelf exchange, which in turn is found to contribute to the shelf hydrographic structure.     

Evaluation of sensitivity of land surface hydrology representations with and without land–atmosphere feedback

This paper compares various ways of quantifying the importance of land–atmosphere feedback. A widely used land surface hydrology model is used in coupled (to a planetary boundary layer model) and uncoupled modes to compare the adequacy of different feedback indices. It is found that existing feedback indices are primarily based on ‘one factor at a time’ sensitivity analysis and cannot adequately capture the interaction between land and atmosphere. A new index is used which combines factorial design concepts and traditional sensitivity analysis. This index is shown to capture and quantify the strength of interaction between land surface parameters and atmosphere. To assess the effects of forcing characteristics on the stand alone model sensitivity, several ways to specify near-surface atmospheric conditions are evaluated. It is found that commonly used forcing conditions (e.g. model generated or observed time-series of near-surface atmospheric variables) may not be adequate to mimic the coupled model environment for evaluating the land surface representations. The partially coupled model sensitivity is shown to capture a major feedback loop related to water holding capacity, surface fluxes and near-surface atmospheric processes. These results suggest that sensitivity from the stand alone model should be interpreted with caution and future evaluations should strive to incorporate land–atmosphere feedback, at least within a partially coupled model. © 1997 John Wiley & Sons, Ltd.     

An assessment of the mixed layer salinity budget in the tropical Pacific Ocean. Observations and modelling (1990–2009)

This paper investigates mechanisms controlling the mixed-layer salinity (MLS) in the tropical Pacific during 1990–2009. We use monthly 1°?×?1° gridded observations of salinity, horizontal current and fresh water flux, and a validated ocean general circulation model with no direct MLS relaxation in both its full resolution (0.25° and 5 days) and re-sampled as the observation time/space grid resolution. The present study shows that the mean spatial distribution of MLS results from a subtle balance between surface forcing (E???P, evaporation minus precipitation), horizontal advection (at low and high frequencies) and subsurface forcing (entrainment and mixing), all terms being of analogous importance. Large-scale seasonal MLS variability is found mainly in the Intertropical and South Pacific Convergence Zones due to changes in their meridional location (and related heavy P), in the North Equatorial Counter Currents, and partly in the subsurface forcing. Maximum interannual variability is found in the western Pacific warm pool and in both convergence zones, in relation to El Niño Southern Oscillation (ENSO) events. In the equatorial band, this later variability is due chiefly to the horizontal advection of low salinity waters from the western to the central-eastern basin during El Niño (and vice versa during La Niña), with contrasted evolution for the Eastern and Central Pacific ENSO types. Our findings reveal that all terms of the MLS equation, including high-frequency (<1 month) salinity advection, have to be considered to close the salinity budget, ruling out the use of MLS (or sea surface salinity) only to directly infer the mean, seasonal and/or interannual fresh water fluxes.     

Seasonal variation of land–atmosphere coupling strength over the West African monsoon region in an atmospheric general circulation model

Abstract

The seasonal variation of land–atmosphere coupling strength has been examined using an extended series of atmospheric general circulation model (AGCM) simulations. In the Western Sahel of Africa, strong coupling strength for precipitation is found in April and May, just prior to and at the beginning of the monsoon season. At this time, heat and water fluxes from the surface are strongly controlled by land conditions, and the unstable conditions in the lower level of the troposphere, as induced by local land state, allow the surface fluxes to influence the variability of convective precipitation—and thus the timing of monsoon onset.

Editor Z. W. Kundzewicz

Citation Yamada, T.J., Kanae, S., Oki, T., and Koster, R.D., 2013. Seasonal variation of land–atmosphere coupling strength over the West African monsoon region in an atmospheric general circulation model. Hydrological Sciences Journal, 58 (6), 1276–1286.     

First‐year sea ice spring melt transitions in the Canadian Arctic Archipelago from time‐series synthetic aperture radar data, 1992–2002

This paper synthesizes 10‐years' worth of interannual time‐series space‐borne ERS‐1 and RADARSAT‐1 synthetic aperture radar (SAR) data collected coincident with daily measurement of snow‐covered, land‐fast first‐year sea ice (FYI) geophysical and surface radiation data collected from the Seasonal Sea Ice Monitoring and Modeling Site, Collaborative‐Interdisciplinary Cryospheric Experiment and 1998 North Water Polynya study over the period 1992 to 2002. The objectives are to investigate the seasonal co‐relationship of the SAR time‐series dataset with selected surface mass (bulk snow thickness) and climate state variables (surface temperature and albedo) measured in situ for the purpose of measuring the interannual variability of sea ice spring melt transitions and validating a time‐series SAR methodology for sea ice surface mass and climate state parameter estimation. We begin with a review of the salient processes required for our interpretation of time‐series microwave backscatter from land‐fast FYI. Our results suggest that time‐series SAR data can reliably measure the timing and duration of surface albedo transitions at daily to weekly time‐scales and at a spatial scales that are on the order of hundreds of metres. Snow thickness on FYI immediately prior to melt onset explains a statistically significant portion of the variability in timing of SAR‐detected melt onset to pond onset for SAR time‐series that are made up of more than 25 images. Our results also show that the funicular regime of snowmelt, resolved in time‐series SAR data at a temporal resolution of approximately 2·5 images per week, is not detectable for snow covers less than 25 cm in thickness. Copyright © 2007 John Wiley & Sons, Ltd.     

Soil moisture variability and scale-dependency of nonlinear parameterizations in coupled land–atmosphere models

Many of the relationships used in coupled land–atmosphere models to describe interactions between the land surface and the atmosphere have been empirically parameterized and thus are inherently dependent on the observational scale for which they were derived and tested. However, they are often applied at scales quite different than the ones they were intended for due to practical necessity. In this paper, a study is presented on the scale-dependency of parameterizations which are nonlinear functions of variables exhibiting considerable spatial variability across a wide range of scales. For illustration purposes, we focus on parameterizations which are explicit nonlinear functions of soil moisture. We use data from the 1997 Southern Great Plains Hydrology Experiment (SGP97) to quantify the spatial variability of soil moisture as a function of scale. By assuming that a parameterization keeps its general form the same over a range of scales, we quantify how the values of its parameters should change with scale in order to preserve the spatially averaged predicted fluxes at any scale of interest. The findings of this study illustrate that if modifications are not made to nonlinear parameterizations to account for the mismatch of scales between optimization and application, then significant systematic biases may result in model-predicted water and energy fluxes.     

Causes of the Regional Variability in Observed Sea Level,Sea Surface Temperature and Ocean Colour Over the Period 1993–2011

We analyse the regional variability in observed sea surface height (SSH), sea surface temperature (SST) and ocean colour (OC) from the ESA Climate Change Initiative datasets over the period 1993–2011. The analysis focuses on the signature of the ocean large-scale climate fluctuations driven by the atmospheric forcing and do not address the mesoscale variability. We use the ECCO version 4 ocean reanalysis to unravel the role of ocean transport and surface buoyancy fluxes in the observed SSH, SST and OC variability. We show that the SSH regional variability is dominated by the steric effect (except at high latitude) and is mainly shaped by ocean heat transport divergences with some contributions from the surface heat fluxes forcing that can be significant regionally (confirming earlier results). This is in contrast with the SST regional variability, which is the result of the compensation of surface heat fluxes by ocean heat transport in the mixed layer and arises from small departures around this background balance. Bringing together the results of SSH and SST analyses, we show that SSH and SST bear some common variability. This is because both SSH and SST variability show significant contributions from the surface heat fluxes forcing. It is evidenced by the high correlation between SST and buoyancy-forced SSH almost everywhere in the ocean except at high latitude. OC, which is determined by phytoplankton biomass, is governed by the availability of light and nutrients that essentially depend on climate fluctuations. For this reason, OC shows significant correlation with SST and SSH. We show that the correlation with SST displays the same pattern as the correlation with SSH with a negative correlation in the tropics and subtropics and a positive correlation at high latitude. We discuss the reasons for this pattern.     

Simulation of crop growth and energy and carbon dioxide fluxes at different time steps from hourly to daily

Understanding the exchange processes of energy and carbon dioxide (CO₂) in the soil–vegetation–atmosphere system is important for assessing the role of the terrestrial ecosystem in the global water and carbon cycle and in climate change. We present a soil–vegetation–atmosphere integrated model (ChinaAgrosys) for simulating energy, water and CO₂ fluxes, crop growth and development, with ample supply of nutrients and in the absence of pests, diseases and weed damage. Furthermore, we test the hypotheses of whether there is any significant difference between simulations over different time steps. CO₂, water and heat fluxes were estimated by the improving parameterization method of the coupled photosynthesis–stomatal conductance–transpiration model. Soil water evaporation and plant transpiration were calculated using a multilayer water and heat‐transfer model. Field experiments were conducted in the Yucheng Integrated Agricultural Experimental Station on the North China Plain. Daily weather and crop growth variables were observed during 1998–2001, and hourly weather variables and water and heat fluxes were measured using the eddy covariance method during 2002–2003. The results showed that the model could effectively simulate diurnal and seasonal changes of net radiation, sensible and latent heat flux, soil heat flux and CO₂ fluxes. The processes of evapotranspiration, soil temperature and leaf area index agree well with the measured values. Midday depression of canopy photosynthesis could be simulated by assessing the diurnal change in canopy water potential. Moreover, the comparisons of simulated daily evapotranspiration and net ecosystem exchange (NEE) under different time steps indicated that time steps used by a model affect the simulated results. There is no significant difference between simulated evapotranspiration using the model under different time steps. However, simulated NEE produces large differences in the response to different time steps. Therefore, the accurate calculation of average absorbed photosynthetic active radiation is important for the scaling of the model from hourly steps to daily steps in simulating energy and CO₂ flux exchanges between winter wheat and the atmosphere. Copyright © 2007 John Wiley & Sons, Ltd.     

Transport over a narrow shelf: Exuma Cays, Bahamas

Current meter and temperature data were collected over a 402-day period from an outer shelf and a tidal channel study site in the Exuma Cays, Bahamas. The shelf width is less than 2 km, and floods and ebbs through a nearby tidal channel extend across the entire shelf and reduce coherence of wind forcing and along-shelf flow. The data are used in perturbation analyses to investigate the across-shelf turbulent transport of heat and momentum over seasonal time scales. Data show a net landward transport of both heat and momentum over the course of the study, but the perturbation products contain distinct seasonal cycles. In fall and winter months, across-shelf heat and momentum fluxes are landward, while during spring and summer months fluxes are seaward. Comparison of shelf-water temperature with the temperature of bank water leaving on the ebb suggests that seasonal cycles of across-shelf heat and momentum in shelf waters are influenced by the seasonal export of relatively warm and cool water from Great Bahama Bank.Responsible Editor: Iris Grabemann     

Nonstationary modeling of a long record of rainfall and temperature over Rome

A long record (1862–2004) of seasonal rainfall and temperature from the Rome observatory of Collegio Romano are modeled in a nonstationary framework by means of the Generalized Additive Models in Location, Scale and Shape (GAMLSS). Modeling analyses are used to characterize nonstationarities in rainfall and related climate variables. It is shown that the GAMLSS models are able to represent the magnitude and spread in the seasonal time series with parameters which are a smooth function of time. Covariate analyses highlight the role of seasonal and interannual variability of large-scale climate forcing, as reflected in three teleconnection indexes (Atlantic Multidecadal Oscillation, North Atlantic Oscillation, and Mediterranean Index), for modeling seasonal rainfall and temperature over Rome. In particular, the North Atlantic Oscillation is a significant predictor during the winter, while the Mediterranean Index is a significant predictor for almost all seasons.     

土壤含水量空间尺度分布非均匀对中尺度脉动量和通量的影响研究

本文采用RAMS（Regional Atmospheric Modeling System）模式模拟研究了没有背景风的情况下,土壤湿度非均匀分布的长度尺度分别为40km逐渐减小到2km时,地表通量的分布和大气边界层的响应.运用二维傅里叶变换,分析了地表通量、中尺度脉动量和中尺度通量的二维幅度谱分布,初步探讨大尺度模式中非均匀地表条件下的边界层参数化问题.分析结果显示各试验的地表水、热通量和中尺度脉动量的幅度谱的极大值都出现在与各自非均匀尺度相对应的波数处,当有不同尺度的非均匀斑块共存时,最大的非均匀尺度占主导.但是中尺度水、热通量的结果有所不同,除了在与各自非均匀尺度相对应的波数处有峰值之外,在其他波数还有多个峰值.这些结果表明地表水、热通量的空间分布尺度与非均匀尺度之间存在较好的对应关系,而中尺度水、热通量与非均匀尺度的关系并不明显,说明地表水、热通量的网格平均值的代表性较好,但是不能反映次网格脉动的影响,而中尺度通量的网格平均值的代表性较差.     

Decadal variability and scales of the sea surface structure in the northern Ionian

The variability and scales of the sea surface structure of the northern Ionian Sea from January 1993 to December 2007 were studied by means of altimeter remotely-sensed weekly Sea Level Anomaly (SLA) objective maps. Variability in the sea surface structure was addressed by means of empirical orthogonal function (EOF) analysis and, assuming an exponential correlation model, scales of the SLA field were quantified as e-folding distances of the SLA autocorrelation function. The variability in the sea surface structure, described by the first three EOFs, which cumulatively explain 60.3% of the data set variance, is characterized by a large-scale structure with variability on a time scale of ∼10-13 years and, on shorter scales, an eddy system with variability on an annual scale. The variability in the large-scale structure describes an overturning of the SLA field, which took place in 1997, and determines a reversal of the geostrophic upper-layer circulation. As the large-scale circulation transition takes place, time-dependent spectral analysis of EOF coefficients shows a redistribution of the spectral energy from inter-annual to semi-annual and monthly components. Spatial scales display variability on an annual and inter-annual time scale. On the annual time scale, variability in spatial scales is characterized by longer values in summer-fall and shorter in winter-spring. Inter-annual variability in spatial scales is demonstrated by a remarkable drop in the values during fall in the period 1998-2000. We propose an explanation of the variability in horizontal scales in terms of the redistribution of water masses and related modifications of the vertical structure of the water column associated with different regimes of the basin-scale circulation.     

Some aspects of the HAPEX-MOBILHY programme: The data base and the modelling strategy

The main thrust of the HAPEX-MOBILHY experiment was towards investigating techniques involved in integrating the turbulent surface fluxes measured at local sites to a larger scale approaching that used in general circulation models.Some aspects of the field data collected at various times and spatial scales are presented. Annual cycle of the soil moisture at many sites is discussed in relation with outputs of a large scale hydrological model. At shorter time scales, the spatial variability of surface energy partition is examined with regard to spatial contrasts in albedo, surface roughness and plant properties related to the two main vegetation classes found in the HAPEX square: A pine forest and the nearby agricultural area.Finally, examples of daily spatial integration with an atmospheric mesoscale model including a comprehensive treatment of land surface processes are presented.     

Simulating hydrologic and hydraulic processes throughout the Amazon River Basin

Presented here is a model framework based on a land surface topography that can be represented with various degrees of resolution and capable of providing representative channel/floodplain hydraulic characteristics on a daily to hourly scale. The framework integrates two models: (1) a water balance model (WBM) for the vertical fluxes and stores of water in and through the canopy and soil layers based on the conservation of mass and energy, and (2) a routing model for the horizontal routing of surface and subsurface runoff and channel and floodplain waters based on kinematic and diffusion wave methodologies. The WBM is driven by satellite‐derived precipitation (TRMM_3B42) and air temperature (MOD08_M3). The model's use of an irregular computational grid is intended to facilitate parallel processing for applications to continental and global scales. Results are presented for the Amazon Basin over the period Jan 2001 through Dec 2005. The model is shown to capture annual runoff totals, annual peaks, seasonal patterns, and daily fluctuations over a range of spatial scales (>1, 000 to < 4·7M km²). For the period of study, results suggest basin‐wide total water storage changes in the Amazon vary by approximately + /? 5 to 10 cm, and the fractional components accounting for these changes are: root zone soil moisture (20%), subsurface water being routed laterally to channels (40%) and channel/floodplain discharge (40%). Annual variability in monthly water storage changes by + /? 2·5 cm is likely due to 0·5 to 1 month variability in the arrival of significant rainfall periods throughout the basin. Copyright © 2009 John Wiley & Sons, Ltd.     

Effects of multi-annual climate variability on the hydrodynamic evolution (1833 to present) in a shallow aquifer system in northern Belgium

Abstract

Using a groundwater flow model and long historical meteorological time series data, the evolution of the groundwater flow regime in a multi-layered groundwater flow basin in northern Belgium during the last one and a half centuries (since 1833) is reconstructed. Model output parameters such as piezometric levels, depth to water table, seepage fluxes in the valleys and calculated baseflow to the river system are presented and inter-annual and decadal variations are evaluated against seasonal fluctuations. The main time-varying boundary condition in the model is the aquifer recharge which was estimated using the method of Thornthwaite and Mather based on precipitation and temperature data. The model does not take into account changes in boundary conditions due to changes in land use (deforestation, drainage of cultivated land) or groundwater exploitation. Variations in model output parameters are therefore only due to climatological forcing. Only the natural non-exploited state of the aquifer is considered. Although few historical piezometric measurements are available to verify model output, the results give an indication of the natural hydrodynamic variations on a time scale of decades.

Citation Van Camp, M., Coetsiers, M., Martens, K. & Walraevens, K. (2010) Effects of multi-annual climate variability on the hydrodynamic evolution (1833 to present) in a shallow aquifer system in northern Belgium. Hydrol. Sci. J. 55(5), 763–779.