The use of delayed drainage in land subsidence modeling greatly complicates model calibration, particularly when the thickness of the fine-grained interbeds varies throughout the modeled region. This thesis documents two separate projects (chapters) related to the use of delayed drainage in groundwater flow and subsidence modeling with parameter estimation. The overall goal of these projects was to better understand how delayed drainage affects accurate parameter estimation and how it is currently affecting the subsidence processes occurring in Las Vegas Valley. Chapter 1 describes an investigation on the value of subsidence data for groundwater model calibration considering delayed drainage. The calibration results of 13 hydraulic parameters of a synthetic conceptual model evaluated for 24 test cases indicate that (1) the inverse of the square of the observation values is a reasonable method to weight the observations, (2) spatially abundant subsidence data typically produce superior parameter estimates even with observation error under constant and cyclical pumping, (3) when subsidence data are limited and combined with drawdown data, outstanding results are obtained for constant pumping conditions. However, for cyclical pumping with observation errors, the best parameter estimates are achieved when multiple years of seasonal subsidence data are provided. The results provide useful suggestions for real-world calibration problems. Chapter 2 outlines the development of an updated flow and subsidence model for Las Vegas Valley covering the entire period of development of the basin. The new model includes a subsidence package that takes into account delayed drainage of fine-grained interbeds. Previous models used subsidence packages that assumed instantaneous equilibration of heads across all hydrogeologic units. The new model resulted in an agreement with measured water-level and improved the simulation of land subsidence. The analysis shows that the typical residual subsidence in Las Vegas Valley can be accurately simulated by incorporating delayed drainage in a long-term model. The study also indicates the need for more sophisticated modeling practices that use delayed drainage with parameter estimation processes to accurately calibrate flow and subsidence models. / Master of Science

The need to consider groundwater and surface water as a single resource has fostered the interest of the scientific community on the interactions between surface water and groundwater. The region below and alongside rivers where surface hydrology and subsurface hydrology concur is the hyporheic zone. This is the region where water exchange determines many biogeochemical and ecological processes of great impact on the functioning of rivers. However, the complex processes taking place in the hyporheic zone require a multidisciplinary approach. The combination of innovative point and distributed techniques originally developed in separated disciplines is of great advantage for the indirect identification of water exchange in the hyporheic zone. Distributed techniques using temperature as a tracer such as fiber-optic distributed temperature sensing can identify the different components of groundwater-surface water interactions based on their spatial and temporal thermal patterns at the sediment-water interface. In particular, groundwater, interflow discharge and local hyporheic exchange flows can be differentiated based on the distinct size, duration and sign of the temperature anomalies. The scale range and resolution of fiber-optic distributed temperature sensing are well complemented by geophysics providing subsurface structures with a similar resolution and scale. Thus, the use of fiber-optic distributed temperature sensing to trace flux patterns supported by the exploration of subsurface structures with geophysics enables spatial and temporal investigation of groundwater-surface water interactions with an unprecedented level of accuracy and resolution. In contrast to the aforementioned methods that can be used for pattern identification at the interface, other methods such as point techniques are required to quantify hyporheic exchange fluxes. In the present PhD thesis, point methods based on hydraulic gradients and thermal profiles are used to quantify hyporheic exchange flows. However, both methods are one-dimensional methods and assume that only vertical flow occurs while the reality is much more complex. The study evaluates the accuracy of the available methods and the factors that impact their reliability. The applied methods allow not only to quantify hyporheic exchange flows but they are also the basis for an interpretation of the sediment layering in the hyporheic zone. For upscaling of the previous results three-dimensional modelling of flow and heat transport in the hyporheic zone combines pattern identification and quantification of fluxes into a single framework. Modelling can evaluate the influence of factors governing groundwater-surface water interactions as well as assess the impact of multiple aspects of model design and calibration of high impact on the reliability of the simulations. But more importantly, this modelling approach enables accurate estimation of water exchange at any location of the domain with unparalleled resolution. Despite the challenges in 3D modelling of the hyporheic zone and in the integration of point and distributed data in models, the benefits should encourage the hyporheic community to adopt an integrative approach comprising from the measurement to the upscaling of hyporheic processes.