Soil moisture is a key state variable that controls runoff formation, infiltration and partitioning of radiation into latent and sensible heat. However, the experimental characterisation of near surface soil moisture patterns and their controls on runoff formation remains a challenge. This subject was one aspect of the BMBF-funded OPAQUE project (operational discharge and flooding predictions in head catchments). As part of that project the focus of this dissertation is on: (1) testing the methodology and feasibility of the Spatial TDR technology in producing soil moisture profiles along TDR probes, including an inversion technique of the recorded signal in heterogeneous field soils, (2) the analysis of spatial variability and temporal dynamics of soil moisture at the field scale including field experiments and hydrological modelling, (3) the application of models of different complexity for understanding soil moisture dynamics and its importance for runoff generation as well as for improving the prediction of runoff volumes. To fulfil objective 1, several laboratory experiments were conducted to understand the influence of probe rod geometry and heterogeneities in the sampling volume under different wetness conditions. This includes a detailed analysis on how these error sources affect retrieval of soil moisture profiles in soils. Concerning objective 2 a sampling strategy of two TDR clusters installed in the head water of the Wilde Weißeritz catchment (Eastern Ore Mountains, Germany) was used to investigate how well “the catchment state” can be characterised by means of distributed soil moisture data observed at the field scale. A grassland site and a forested site both located on gentle slopes were instrumented with two Spatial TDR clusters that consist of up to 39 TDR probes. Process understanding was gained by modelling the interaction of evapotranspiration and soil moisture with the hydrological process model CATFLOW. A field scale irrigation experiment was carried out to investigate near subsurface processes at the hillslope scale. The interactions of soil moisture and runoff formation were analysed using discharge data from three nested catchments: the Becherbach with a size of 2 km², the Rehefeld catchment (17 km²) and the superordinate Ammelsdorf catchment (49 km²). Statistical analyses including observations of pre-event runoff, soil moisture and different rainfall characteristics were employed to predict stream flow volume. On the different scales a strong correlation between the average soil moisture and the runoff coefficients of rainfall-runoff events could be found, which almost explains equivalent variability as the pre-event runoff. Furthermore, there was a strong correlation between surface soil moisture and subsurface wetness with a hysteretic behaviour between runoff soil moisture. To fulfil objective 3 these findings were used in a generalised linear model (GLM) analysis which combines state variables describing the catchments antecedent wetness and variables describing the meteorological forcing in order to predict event runoff coefficients. GLM results were compared to simulations with the catchment model WaSiM ETH. Hereby were the model results of the GLMs always better than the simulations with WaSiM ETH. The GLM analysis indicated that the proposed sampling strategy of clustering TDR probes in typical functional units is a promising technique to explore soil moisture controls on runoff generation and can be an important link between the scales. Long term monitoring of such sites could yield valuable information for flood warning and forecasting by identifying critical soil moisture conditions for the former and providing a better representation of the initial moisture conditions for the latter. / Abflussentwicklung, Infiltration und die Umverteilung von Strahlung in latenten und sensiblen Wärmestrom werden maßgeblich durch die Bodenfeuchte der vadosen Zone gesteuert. Trotz allem, gibt s wenig Arbeiten die sich mit der experimentellen Charakterisierung der Bodenfeuchteverteilung und ihre Auswirkung auf die Abflussbildung beschäftigen. Der Fokus dieser Dissertation wurde darauf ausgerichtet: (1) die Methode des Spatial TDR und deren Anwendbarkeit einschließlich der Inversion des TDR Signals in heterogenen Böden zu prüfen, (2) die Analyse der räumlichen und zeitlichen Dynamik der Bodenfeuchte auf der Feldskala einschließlich Feldexperimenten und hydrologischer Modellierung, (3) der Aufbau verschiedener Modellanwendungen unterschiedlicher Komplexität um die Bodenfeuchtedynamiken und die Abflussentwicklung zu verstehen und die Vorhersage des Abflussvolumens zu verbessern. Um die Zielsetzung 1 zu erreichen, wurden verschiedene Laborversuche durchgeführt. Hierbei wurde der Einfluss der Sondenstabgeometrie und verschiedener Heterogenitäten im Messvolumen bei verschiedenen Feuchtegehalten untersucht. Dies beinhaltete eine detaillierte Analyse wie diese Fehlerquellen die Inversion des Bodenfeuchteprofils beeinflussen. Betreffend der Zielsetzung 2, wurden 2 TDR-Cluster in den Quellgebieten der Wilden Weißeritz installiert (Osterzgebirge) und untersucht, wie gut der Gebietszustand mit räumlich hochaufgelösten Bodenfeuchtedaten der Feldskala charakterisiert werden kann. Um die Interaktion zwischen Evapotranspiration und Bodenfeuchte zu untersuchen wurde das hydrologische Prozessmodell CATFLOW angewendet. Ein Beregnungsversuch wurde durchgeführt um die Zwischenabflussprozesse auf der Hangskala zu verstehen. Die Interaktion zwischen Bodenfeuchte und Abflussentwicklung wurde anhand von drei einander zugeordneten Einzugsgebieten analysiert. Statistische Analysen unter Berücksichtigung von Basisabfluss, Bodenvorfeuchte und verschiedenen Niederschlagscharakteristika wurden verwendet, um auf das Abflussvolumen zu schließen. Auf den verschiedenen Skalen konnte eine hohe Korrelation zwischen der mittleren Bodenfeuchte und dem Abflussbeiwert der Einzelereignisse festgestellt werden. Hierbei konnte die Bodenfeuchte genauso viel Variabilität erklären wie der Basisabfluss. Im Hinblick auf Zielsetzung 3 wurden “Generalised liner models” (GLM) genutzt. Dabei wurden Prädiktorvariablen die den Gebietszustand beschreiben und solche die die Meteorologische Randbedingungen beschreiben genutzt um den Abflussbeiwert zu schätzen. Die Ergebnisse der GLMs wurden mit Simulationsergebnissen des hydrologischen Gebietsmodells WaSiM ETH verglichen. Hierbei haben die GLMs eindeutig bessere Ergebnisse geliefert gegenüber den WaSiM Simulationen. Die GLM Analysen haben aufgezeigt, dass die verwendete Messstrategie mehrerer TDR-Cluster in typischen funktionalen Einheiten eine viel versprechende Methode ist, um den Einfluss der Bodenfeuchte auf die Abflussentwicklung zu verstehen und ein Bindeglied zwischen den Skalen darstellen zu können. Langzeitbeobachtungen solcher Standorte sind in der Lage wichtige Zusatzinformationen bei der Hochwasserwarnung und -vorhersage zu liefern durch die Identifizierung kritischer Gebietszustände für erstere und eine bessere Repräsentation der Vorfeuchte für letztere.

Soil conditions under vegetation cover and their spatial and temporal variations from point to catchment scale are crucial for understanding hydrological processes within the vadose zone, for managing irrigation and consequently maximizing yield by precision farming. Soil moisture and soil roughness are the key parameters that characterize the soil status. In order to monitor their spatial and temporal variability on large scales, remote sensing techniques are required. Therefore the determination of soil parameters under vegetation cover was approached in this thesis by means of (multi-angular) polarimetric SAR acquisitions at a longer wavelength (L-band, lambda=23cm). In this thesis, the penetration capabilities of L-band are combined with newly developed (multi-angular) polarimetric decomposition techniques to separate the different scattering contributions, which are occurring in vegetation and on ground. Subsequently the ground components are inverted to estimate the soil characteristics. The novel (multi-angular) polarimetric decomposition techniques for soil parameter retrieval are physically-based, computationally inexpensive and can be solved analytically without any a priori knowledge. Therefore they can be applied without test site calibration directly to agricultural areas. The developed algorithms are validated with fully polarimetric SAR data acquired by the airborne E-SAR sensor of the German Aerospace Center (DLR) for three different study areas in Germany. The achieved results reveal inversion rates up to 99% for the soil moisture and soil roughness retrieval in agricultural areas. However, in forested areas the inversion rate drops significantly for most of the algorithms, because the inversion in forests is invalid for the applied scattering models at L-band. The validation against simultaneously acquired field measurements indicates an estimation accuracy (root mean square error) of 5-10vol.% for the soil moisture (range of in situ values: 1-46vol.%) and of 0.37-0.45cm for the soil roughness (range of in situ values: 0.5-4.0cm) within the catchment. Hence, a continuous monitoring of soil parameters with the obtained precision, excluding frozen and snow covered conditions, is possible. Especially future, fully polarimetric, space-borne, long wavelength SAR missions can profit distinctively from the developed polarimetric decomposition techniques for separation of ground and volume contributions as well as for soil parameter retrieval on large spatial scales. / Zur Verbesserung der hydrologischen Abflussmodellierung, der Flutvorhersage, der gezielten Bewässerung von landwirtschaftlichen Nutzflächen und zum Schutz vor Ernteausfällen ist die Bestimmung der Bodenfeuchte und der Bodenrauhigkeit von grosser Bedeutung. Aufgrund der hohen zeitlichen sowie räumlichen Dynamik dieser Bodenparameter ist eine flächenhafte Erfassung mit hoher Auflösung und in kurzen zeitlichen Abständen notwendig. In situ Messtechniken stellen eine sehr zeit- und personalaufwändige Alternative dar, deshalb werden innovative Fernerkundungsverfahren mit aktivem Radar erprobt. Diese Aufnahmetechniken sind von Wetter- und Beleuchtungsverhältnissen unabhängig und besitzen zudem die Möglichkeit, abhängig von der Wellenlänge, in Medien einzudringen. Mit dem in dieser Arbeit verwendeten polarimetrischen Radar mit synthetischer Apertur (PolSAR) werden die Veränderungen der Polarisationen ausgewertet, da diese aufgrund der physikalischen Eigenschaften der reflektierenden Medien objektspezifisch verändert und gestreut werden. Es kann dadurch ein Bezug zwischen der empfangenen Radarwelle und den dielektrischen Eigenschaften (Feuchtegehalt) sowie der Oberflächengeometrie (Rauhigkeit) des Bodens hergestellt werden. Da vor allem in den gemässigten Klimazonen die landwirtschaftlichen Nutzflächen die meiste Zeit des Jahres mit Vegetation bestanden sind, wurden in dieser Dissertation Verfahren entwickelt, um die Bodenfeuchte und die Bodenrauhigkeit unter der Vegetation erfassen zu können. Um die einzelnen Rückstreubeiträge der Vegetation und des Bodens voneinander zu trennen, wurde die Eindringfähigkeit von längeren Wellenlängen (L-band, lambda=23cm) mit neu entwickelten (multi-angularen) polarimetrischen Dekompositionstechniken kombiniert, um die Komponente des Bodens zu extrahieren und auszuwerten. Für die Auswertung wurden polarimetrische Streumodelle benutzt, um die Bodenkomponente zu modellieren und dann mit der extrahierten Bodenkomponente der aufgenommenen Daten zu vergleichen. Die beste Übereinstimmung von Modell und Daten wurde als die gegebene Bodencharakteristik gewertet und dementsprechend invertiert. Die neu entwickelten, polarimetrischen Dekompositionstechniken für langwelliges polarimetrisches SAR basieren auf physikalischen Prinzipien, benötigen wenig Rechenzeit, erfordern keine Kalibrierung und sind ohne Verwendung von a priori Wissen analytisch lösbar. Um die entwickelten Algorithmen zu testen, wurden in drei verschiedenen Untersuchungsgebieten in Deutschland mit dem flugzeuggetragenen E-SAR Sensor des Deutschen Zentrums für Luft- und Raumfahrt (DLR) polarimetrische SAR Daten aufgenommen. Die Auswertungen der PolSAR Daten haben bestätigt, dass die besten Invertierungsergebnisse mit langen Wellenlängen erzielt werden können (L-Band). Des Weiteren konnten bei der Bestimmung der Bodenfeuchte und der Bodenrauhigkeit hohe Inversionsraten erreicht werden (bis zu 99% der Untersuchungsfläche). Es hat sich gezeigt, dass die polarimetrischen Streumodelle bei der gegebenen Wellenlänge nicht für bewaldete Gebiete geeignet sind, was die Anwendbarkeit des Verfahrens auf landwirtschaftliche Nutzflächen einschränkt. Die Validierung mit Bodenmessungen in den Untersuchungsgebieten, die zeitgleich zu den PolSAR Aufnahmen durchgeführt wurden, hat ergeben, dass eine kontinuierliche Beobachtung des Bodenzustandes (ausgenommen in Zeiten mit gefrorenem oder Schnee bedecktem Boden) mit einer Genauigkeit (Wurzel des mittleren quadratischen Fehlers) von 5-10vol.% für die Bodenfeuchte (in situ Messbereich: 1-46vol.%) und von 0.37-0.45cm für die Bodenrauhigkeit (in situ Messbereich: 0.5-4.0cm) möglich ist. Besonders künftige Fernerkundungsmissionen mit langwelligem, voll polarimetrischem SAR können von den entwickelten Dekompositionstechniken profitieren, um die Vegetationskomponente von der Bodenkomponente zu trennen und die Charakteristik des Oberbodens flächenhaft zu bestimmen.

The temporal dynamics of hydrological model performance gives insights into errors that cannot be obtained from global performance measures assigning a single number to the fit of a simulated time series to an observed reference series. These errors can include errors in data, model parameters, or model structure. Dealing with a set of performance measures evaluated at a high temporal resolution implies analyzing and interpreting a high dimensional data set. This paper presents a method for such a hydrological model performance assessment with a high temporal resolution and illustrates its application for two very different rainfall-runoff modeling case studies. The first is the Wilde Weisseritz case study, a headwater catchment in the eastern Ore Mountains, simulated with the conceptual model WaSiM-ETH. The second is the Malalcahuello case study, a headwater catchment in the Chilean Andes, simulated with the physicsbased model Catflow. The proposed time-resolved performance assessment starts with the computation of a large set of classically used performance measures for a moving window. The key of the developed approach is a data-reduction method based on self-organizing maps (SOMs) and cluster analysis to classify the high-dimensional performance matrix. Synthetic peak errors are used to interpret the resulting error classes. The final outcome of the proposed method is a time series of the occurrence of dominant error types. For the two case studies analyzed here, 6 such error types have been identified. They show clear temporal patterns, which can lead to the identification of model structural errors.

Graduation date: 1992

The mining of oil sands in the sub-humid region of Northern Alberta, Canada causes large-scale landscape disturbance, which subsequently requires extensive reclamation to re-establish the surface and subsurface hydrology. The reconstructed watersheds examined in this study are located at the Syncrude Canada Limited mine site, 40 km North of Fort McMurray, Alberta, Canada. The three experimental reconstructed watersheds, with nominal soil thicknesses of 1.0 m, 0.50 m and 0.35 m comprised a thin layer of peat (15-20 cm) over varying thicknesses of secondary (till) soil, have been constructed to cover saline sodic overburden and to provide sufficient moisture storage for vegetation while minimizing surface runoff and deep percolation to the underlying shale overburden. In order to replicate the hydrological behavior, assess the sustainability, and trace the evolution over time of the reclaimed watersheds, a suitable modeling tool is needed.</p> <p>In this research, a model is developed using the system dynamics approach to simulate the hydrological processes in the three experimental reconstructed watersheds and to assess their ability to provide the various watershed functions. The model simulates the vertical and lateral water movement, surface runoff and evapotranspiration within each watershed. Actual evapotranspiration, which plays an important role in the hydrology of the Canadian semi-arid regions, is simulated using an indexed soil moisture method. The movement of water within the various soil layers of the cover is based on parametric relationships in conjunction with conceptual infiltration models. The feedback relationships among the various dynamic hydrologic processes in the watershed are captured in the developed System Dynamic Watershed Model (SDWM). </p> <p>Most hydrological models are evaluated using runoff as the determining criterion for model calibration and validation, while accounting for the movement of moisture in the soil as a water loss. Since one of the primary objectives of a reconstructed watershed is to maintain the natural flora and fauna, it is important to recognize that soil moisture plays an important role in assessing the performance of the reconstructed watersheds. In turn, soil moisture becomes an influential factor for quantifying the health of the reconstructed watershed. The developed model has been calibrated and validated with data for two years (2001-2002), upholding the sensitive relationship between soil moisture and runoff. Accurate calibration of the model based on simulations of soil moisture in the various soil layers improves its overall performance. The model was subsequently used to simulate the three sub-watersheds for five years, with changing the calibrated model parameters to use them as indicators of watershed evolution. The simulated results were compared with the observed values. </p> <p>The results of the study illustrate that all three watersheds are still evolving. Failure to identify a unique parameter set for simulating the watershed response supports the hypothesis of watershed evolution. Soil moisture exchange between the till and peat layers changed with time in all of the watersheds. There was also a modest change in the water movement from the till to shale layers in each of the sub-watersheds. Vegetation is increasing in all of watersheds although there is an indication that one of the sub-watersheds may be sustaining deep rooted vegetation. The results demonstrate the successful application of the system dynamics approach and the developed model in simulating the hydrology of reconstructed watersheds and the potential for using this approach in assessing complex hydrologic systems.

Near-surface sediments in the vadose zone play a fundamental role in the hydrologic system. The shallow vadose zone can act as a buffer to delay or attenuate surface contaminants before they reach the water table. It also acts as a temporary soil moisture reservoir for plant and atmospheric uptake, and regulates the seasonal groundwater recharge process. Over the past few decades, geophysical methods have received unprecedented attention as an effective vadose zone characterization tool offering a range of non-invasive to minimally invasive techniques with the capacity to provide detailed soil moisture information at depths typically unattainable using conventional point-measurement sensors. Ground-penetrating radar (GPR) has received much of this attention due to its high sensitivity to the liquid water phase in geologic media. While much has been learned about GPR soil moisture monitoring and characterization techniques, it has not been evaluated across highly dynamic natural soil conditions. Consequently, GPR’s capacity to characterize a complete range of naturally occurring vadose zone conditions including wetting/drying and freeze/thaw cycles, is not yet fully understood. Further, the nature of GPR response during highly dynamic moisture periods has not been thoroughly investigated. The objective of this thesis is to examine the capacity of various surface GPR techniques and methodologies for the characterization of soil moisture dynamics in the upper few meters of vadose zone, and to develop measurement strategies capable of providing quantitative information about the current and future state of the shallow hydrologic system. To achieve this, an exhaustive soil moisture monitoring campaign employing a range of GPR antenna frequencies and survey acquisition geometries was initiated at three different agricultural field sites located in southern Ontario, Canada, between May 2006 and October 2008. This thesis represents the first attempt to evaluate multiple annual cycles of soil conditions and associated hydrological processes using high-frequency GPR measurements. Summaries of the seven major works embodied in this thesis are provided below. Direct ground wave (DGW) measurements obtained with GPR have been used in a number of previous studies to monitor volumetric water content changes in the root zone; however, these studies have involved controlled field experiments or measurements collected across limited ranges in soil moisture. To further investigate the capacity of the DGW method, multi-frequency (i.e., 225 MHz, 450 MHz and 900 MHz) common-midpoint (CMP) measurements were used to monitor a complete annual cycle of soil water content variations at three sites with different soil textures (i.e., sand, sandy loam and silt loam). CMP surveys permitted characterization of the nature and evolution of the near-surface electromagnetic wavefields, and their subsequent impact on DGW velocity measurements. GPR results showed significant temporal variations in both the near-surface wavefield and multi-frequency DGW velocities corresponding to both seasonal and shorter term variations in soil conditions. While all of the measurement sites displayed similar temporal responses, the rate and magnitude of these velocity variations corresponded to varying soil water contents which were primarily controlled by the soil textural properties. Overall, the DGW measurements obtained using higher frequency antennas were less impacted by near-surface wavefield interference due to their shorter signal pulse duration. The estimation of soil water content using GPR velocity requires an appropriate petrophysical relationship between the dielectric permittivity and volumetric water content of the soil. The ability of various empirical relationships, volumetric mixing formulae and effective medium approximations were evaluated to predict near-surface volumetric soil water content using high-frequency DGW velocity measurements obtained from CMP soundings. Measurements were collected using 225, 450 and 900 MHz antennas across sand, sandy loam and silt loam soil textures over a complete annual cycle of soil conditions. A lack of frequency dependence in the results indicated that frequency dispersion had minimal impact on the data set. However, the accuracy of soil water content predictions obtained from the various relationships ranged considerably. The best fitting relationships did exhibit some degree of textural bias that should be considered in the choice of petrophysical relationship for a given data set. Further improvements in water content estimates were obtained using a field calibrated third-order polynomial relationship and three-phase volumetric mixing formula. While DGW measurements provide valuable information within the root zone, the characterization of vertical moisture distribution and dynamics requires a different approach. A common approach utilizes normal-moveout (NMO) velocity analysis of CMP sounding data. To further examine this approach, an extensive field study using multi-frequency (i.e., 225 MHz, 450 MHz, 900 MHz) CMP soundings was conducted to monitor a complete annual cycle of vertical soil moisture conditions at the sand, sandy loam and silt loam sites. The use of NMO velocity analysis was examined for monitoring highly dynamic vertical soil moisture conditions consisting of wetting/drying and freeze/thaw cycles with varying degrees of magnitude and vertical velocity gradient. NMO velocity analysis was used to construct interval-velocity-depth models at a fixed location collected every 1 to 4 weeks. Time-lapse models were combined to construct temporal interval-velocity fields, which were converted into soil moisture content. These moisture fields were used to characterize the vertical distribution, and dynamics of soil moisture in the upper few meters of vadose zone. Although the use of multiple antenna frequencies provided varying investigation depths and vertical resolving capabilities, optimal characterization of soil moisture conditions was obtained with 900 MHz antennas. The integration of DGW and NMO velocity data from a single CMP sounding could be used to assess the nature of shallow soil moisture coupling with underlying vadose zone conditions; however, a more quantitative analyses of the surface moisture dynamics would require definitive knowledge of GPR sampling depth. Although surface techniques have been used by a number of previous researchers to characterize soil moisture content in the vadose zone, limited temporal sampling and low resolution near the surface in these studies impeded the quantitative analysis of vertical soil moisture distribution and its associated dynamics within the shallow subsurface. To further examine the capacity of surface GPR, an extensive 26 month field study was undertaken using concurrent high-frequency (i.e., 900 MHz) reflection profiling and CMP soundings to quantitatively monitor soil moisture distribution and dynamics within a sandy vadose zone environment. An analysis on the concurrent use of reflection and CMP measurements was conducted over two contrasting annual cycles of soil conditions. Reflection profiles provided high resolution traveltime data between four stratigraphic reflection events while cumulative results of the CMP sounding data set produced precise depth estimates for those reflecting interfaces, which were used to convert interval traveltime data into soil water content estimates. The downward propagation of episodic infiltration events associated with seasonal and transient conditions were well resolved by the GPR data. The GPR data also revealed variations in the nature of these infiltration events between contrasting annual cycles. The use of CMP soundings also permitted the determination of DGW velocities, which enabled better characterization of short-duration wetting/drying and freezing/thawing processes. This higher resolution information can be used to examine the nature of the coupling between shallow and deep moisture conditions. High-resolution surface GPR measurements were used to examine vertical soil moisture distribution and its associated dynamics within the shallow subsurface over a 26 month period. While the apparent ability of surface GPR methods to give high quality estimates of soil moisture distribution in the upper 3 meters of the vadose zone was demonstrated, the nature of these GPR-derived moisture data needed to be assessed in the context of other hydrological information. As a result, GPR soil moisture estimates were compared with predictions obtained from a well-accepted hydrological modeling package, HYDRUS-1D (Simunek et al., 2008). The nature of transient infiltration pulses, evapotranspiration episodes, and deep drainage patterns were examined by comparing them with vertical soil moisture flow simulations. Using laboratory derived soil hydraulic property information from soil samples and a number of simplifying assumptions about the system, very good agreement was achieved between measured and simulated soil moisture conditions without model calibration. The overall good agreement observed between forward simulations and field measurements over the vertical profile validated the capacity of surface GPR to provide detailed information about hydraulic state conditions in the upper few meters of vadose zone. A unique DGW propagation phenomenon was observed during early soil frost formation. High-frequency DGW measurements were used to monitor the seasonal development of a thin, high velocity frozen soil layer over a wet low velocity unfrozen substratum. During the freezing process, the progressive attenuation of a low velocity DGW and the subsequent development of a high velocity DGW were observed. Numerical simulations using GPRMAX2D (Giannopoulos, 2005) showed that low velocity DGW occurring after freezing commenced was due to energy leaking across the frozen layer from the spherical body wave in the unfrozen half space. This leaky phase progressively dissipated until the frozen layer reached a thickness equivalent to one quarter of the dominant wavelength in the frozen ground. The appearance of the high velocity DGW was governed by its destructive interference with the reflection events from the base of the frozen layer. This interference obscured the high velocity DGW until the frozen layer thickness reached one half of the dominant wavelength in the frozen ground. While GPR has been extensively used to study frozen soil conditions in alpine environments, its capacity to characterize highly dynamic shallow freeze-thaw processes typically observed in temperate environments is not well understood. High-frequency reflection profiles and CMP soundings were used to monitor the freezing and thawing process during the winter seasonal period at the sand and silt loam sites. Reflection profiles revealed the long-term development of a very shallow (<0.5 m) soil frost zone overlying unfrozen wet substratum. During the course of the winter season, long-term traveltime analysis yielded physical properties of the frozen and unfrozen layers as well as the spatial distribution of the base of the soil frost zone. Short-term shallow thawing events overlying frozen substratum formed a dispersive waveguide for both the CMP and reflection profile surveys. Inversion of the dispersive wavefields for the CMP data yielded physical property estimates for the thawed and frozen soils and thawed layer thickness. It was shown that GPR can be used to monitor very shallow freezing and thawing events by responding to changes in the relative dielectric permittivity of the soil water phase. The works embodied in this thesis demonstrate the effectiveness of high-frequency GPR as a non-invasive soil moisture monitoring tool under a full range of naturally occurring moisture conditions with the temporal and vertical resolution necessary to quantitatively examine shallow vadose zone moisture dynamics. Because this study encompassed an unprecedented range of naturally occurring soil conditions, including numerous short and long duration wetting/drying and freezing/thawing cycles, complex geophysical responses were observed during highly dynamic soil moisture processes. Analysis and interpretation of these geophysical responses yielded both qualitative and quantitative information about the state of the hydrologic system, and hence, provided a non-invasive means of characterizing soil moisture processes in shallow vadose zone environments. In the future, these GPR soil moisture monitoring strategies should be incorporated into advanced land-surface hydrological modeling studies to improve our understanding of shallow hydrologic systems and its impacts on groundwater resources.

The Soil Moisture and Ocean Salinity (SMOS) mission is an Earth Explorer Opportunity mission from the European Space Agency (ESA). It was a direct response to the global observations of soil moisture and ocean salinity. Its goal is to produce global of these parameters using a dual-polarization L-band interferometric radiometer the Microwave Imaging Radiometer by Aperture Synthesis (MIRAS). This instrument is a new polarimetric two-dimensional (2-D) Y-shaped synthetic aperture interferometric radiometer based on the techniques used in radio-astronomy to obtain high resolution avoiding large antenna structures. MIRAS measures remotely the brightness temperature (TB) emitted by the Earth's surface, which is not isotropic, since it depends on the incidence angle and polarization, the Soil Moisture (SM) or the Sea Surface (SSS), the surface roughness etc. among others. The scope of this doctoral thesis is the study of some potential improvements could eventually be implemented in future interferometric radiometers. To validate improvements a ground-based instrument concept demonstrator the Passive Advanced Unit Synthetic Aperture or (PAU-SA) has being designed and implemented. Both MIRAS and PAU-SA are Y-shaped array, but the receiver topology and the processing unit are different. This Ph.D. thesis has been developed in the frame of The European Investigator Awards (EURYI) 2004 project entitled "Passive Advanced Unit (PAU): Hybrid L-band Radiometer, GNSS Refectometer and IR-Radiometer for Passive Sensing of the Ocean", and supported by the European Science Foundation (ESF).

The mining of oil sands in the sub-humid region of Northern Alberta, Canada causes large-scale landscape disturbance, which subsequently requires extensive reclamation to re-establish the surface and subsurface hydrology. The reconstructed watersheds examined in this study are located at the Syncrude Canada Limited mine site, 40 km North of Fort McMurray, Alberta, Canada. The three experimental reconstructed watersheds, with nominal soil thicknesses of 1.0 m, 0.50 m and 0.35 m comprised a thin layer of peat (15-20 cm) over varying thicknesses of secondary (till) soil, have been constructed to cover saline sodic overburden and to provide sufficient moisture storage for vegetation while minimizing surface runoff and deep percolation to the underlying shale overburden. In order to replicate the hydrological behavior, assess the sustainability, and trace the evolution over time of the reclaimed watersheds, a suitable modeling tool is needed.</p> <p>In this research, a model is developed using the system dynamics approach to simulate the hydrological processes in the three experimental reconstructed watersheds and to assess their ability to provide the various watershed functions. The model simulates the vertical and lateral water movement, surface runoff and evapotranspiration within each watershed. Actual evapotranspiration, which plays an important role in the hydrology of the Canadian semi-arid regions, is simulated using an indexed soil moisture method. The movement of water within the various soil layers of the cover is based on parametric relationships in conjunction with conceptual infiltration models. The feedback relationships among the various dynamic hydrologic processes in the watershed are captured in the developed System Dynamic Watershed Model (SDWM). </p> <p>Most hydrological models are evaluated using runoff as the determining criterion for model calibration and validation, while accounting for the movement of moisture in the soil as a water loss. Since one of the primary objectives of a reconstructed watershed is to maintain the natural flora and fauna, it is important to recognize that soil moisture plays an important role in assessing the performance of the reconstructed watersheds. In turn, soil moisture becomes an influential factor for quantifying the health of the reconstructed watershed. The developed model has been calibrated and validated with data for two years (2001-2002), upholding the sensitive relationship between soil moisture and runoff. Accurate calibration of the model based on simulations of soil moisture in the various soil layers improves its overall performance. The model was subsequently used to simulate the three sub-watersheds for five years, with changing the calibrated model parameters to use them as indicators of watershed evolution. The simulated results were compared with the observed values. </p> <p>The results of the study illustrate that all three watersheds are still evolving. Failure to identify a unique parameter set for simulating the watershed response supports the hypothesis of watershed evolution. Soil moisture exchange between the till and peat layers changed with time in all of the watersheds. There was also a modest change in the water movement from the till to shale layers in each of the sub-watersheds. Vegetation is increasing in all of watersheds although there is an indication that one of the sub-watersheds may be sustaining deep rooted vegetation. The results demonstrate the successful application of the system dynamics approach and the developed model in simulating the hydrology of reconstructed watersheds and the potential for using this approach in assessing complex hydrologic systems.

<p>Topography, vegetation, nutrient dynamics, soil features and hydroclimatic forcing are inherently coupled, with feedbacks occurring over a wide range of temporal and spatial scales. Vegetation growth may be limited by soil moisture, nutrient or solar radiation availability, and in turn influences both soil moisture and nutrient balances at a point. These dynamics are further complicated in a complex terrain, through a series of spatial interactions. A number of experiments has characterized the feedbacks between soil moisture and vegetation dynamics, but a theoretical framework linking short-term leaf-level to interannual plot-scale dynamics has not been fully developed yet. Such theory is needed for optimal management of water resources in natural ecosystems and for agricultural, municipal and industrial uses. Also, it complements the current knowledge on ecosystem response to the predicted climate change.</p><p>In this dissertation, the response of vegetation dynamics to unpredictable environmental fluctuations at multiple space-time scales is explored in a modeling framework from sub-daily to interannual time scales. At the hourly time scale, a simultaneous analysis of photosynthesis, transpiration and soil moisture dynamics is carried out to explore the impact of water stress on different photosynthesis processes at the leaf level, and the overall plant activity. Daily soil moisture and vegetation dynamics are then scaled up to the growing season using a stochastic model accounting for daily to interannual hydroclimatic variability. Such stochastic framework is employed to explore the impact of rainfall patterns and different irrigation schemes on crop productivity, along with their implications in terms of sustainability and profitability. To scale up from point to landscape, a probabilistic representation of local landscape features (i.e., slope and aspect) is developed, and applied to assess the effects of topography on solar radiation. Finally, a minimalistic ecosystem model, including soil moisture, vegetation and nutrient dynamics at the year time scale, is outlined; when coupled to the proposed probabilistic topographic description, the latter model can serve to assess the relevance of spatial interactions and to single out the main biophysical controls responsible for ecohydrological variability at the landscape scale.</p> / Dissertation

This thesis research consists of two separate studies. The first study presents the assessment and representation of the effects of soil macropores on the soil hydraulic properties in land surface models for more accurate simulations of soil moisture and surface hydrology. Hydraulic properties determine the soil water content and its transport in the soil. They are provided in most current climate models as empirical formulas by functions of the soil texture. Such is not realistic if the soil contains a substantial amount of macropores. A two-mode soil pore size distribution is incorporated into a land surface model and tested using an observational dataset at a tropical forest site with aggregated soils. The result showed that the existence of macropores greatly affects the estimation of hydraulic properties. Their influence can be included in land models by adding a second function to the pore-size distribution. A practical hydraulic scheme with macropore considerations was proposed given that the existing schemes are not applicable for large-scale simulations. The developed scheme was based on the physical attributes of the water in soil capillary pores and the statistics of several global soil databases. The preliminary test showed that it captures part of soil macropore hydraulic features without sacrificing the estimation accuracy of hydraulic properties of water in soil matrix. The second study presents the development of an integrated land/ecosystem model by combining the advanced features of a biophysically based land model, the Community Land Model, and an ecosystem biochemical model. The results from tests of the integrated model at four forest sites showed that the model reasonably captures the seasonal and interannual dynamics of leaf area index and leaf nitrogen control on carbon assimilation across different environments. With being coupled to an atmospheric general circulation model (AGCM), the integrated model showed a strong ability to simulate terrestrial ecosystem carbon fluxes together with heat and water fluxes. Its simulated land surface physical variables are reasonable in both geographic distribution and temporal variation with considering the interactive vegetation parameters.