Imaging lateral groundwater flow in the shallow subsurface using stochastic temperature fields

Although temperature has often been used as an indication of vertical groundwater movement, its usefulness for identifying horizontal fluid flow has been limited by the difficulty of obtaining sufficient data to draw defensible conclusions. Here we use stochastic simulation to develop a high-resolution image of fluid temperatures in the shallow subsurface at Borax Lake, Oregon. The temperature field inferred from the geostatistical simulations clearly shows geothermal fluids discharging from a group of fault-controlled hydrothermal springs, moving laterally through the subsurface, and mixing with shallow subsurface flow originating from nearby Borax Lake. This interpretation of the data is supported by independent geochemical and isotopic evidence, which show a simple mixing trend between Borax Lake water and discharge from the thermal springs. It is generally agreed that stochastic simulation can be a useful tool for extracting information from complex and/or noisy data and, although not appropriate in all situations, geostatistical analysis may provide good definition of flow paths in the shallow subsurface. Although stochastic imaging techniques are well known in problems involving transport of species, e.g. delineation of contaminant plumes from soil gas survey data, we are unaware of previous applications to the transport of thermal energy for the purpose of inferring shallow groundwater flow.