Xiaojing Fu

Pattern formation and coarsening dynamics in three-dimensional convective mixing in porous media

Geological carbon dioxide (CO2) sequestration entails capturing and injecting CO2 into deep saline aquifers for long-term storage. The injected CO2 partially dissolves in groundwater to form a mixture that is denser than the initial groundwater. The local increase in density triggers a gravitational instability at the boundary layer that further develops into columnar plumes of CO2-rich brine, a process that greatly accelerates solubility trapping of the CO2. Here, we investigate the pattern-formation aspects of convective mixing during geological CO2 sequestration by means of high-resolution three-dimensional simulation. We find that the CO2 concentration field self-organizes as a cellular network structure in the diffusive boundary layer at the top boundary. By studying the statistics of the cellular network, we identify various regimes of finger coarsening over time, the existence of a non-equilibrium stationary state, and a universal scaling of three-dimensional convective mixing.

Improved characterization of heterogeneous permeability in saline aquifers from transient pressure data during freshwater injection

Managing recharge of freshwater into saline aquifers requires accurate estimation of the heterogeneous permeability field for maximizing injection and recovery efficiency. Here, we present a methodology for subsurface characterization in saline aquifers that takes advantage of the density difference between the injected freshwater and the ambient saline groundwater. We combine high resolution forward modeling of density-driven flow with an efficient Bayesian geostatistical inversion algorithm. In the presence of a density difference between the injected and ambient fluids due to differences in salinity, the pressure field is coupled to the spatial distribution of salinity. This coupling renders the pressure field transient: the time evolution of the salinity distribution controls the density distribution which then leads to a time-evolving pressure distribution.We exploit this coupling between pressure and salinity to obtain an improved characterization of the permeability field without multiple pumping tests or additional salinity measurements. We show that the inversion performance improves with an increase in the mixed convection ratio–the relative importance between viscous forces from injection and buoyancy forces from density difference. Our work shows that measuring transient pressure data at multiple sampling points during freshwater injection into saline aquifers can be an effective strategy for aquifer characterization, key to the successful management of aquifer recharge.

Uncertainty in modeled and observed climate change impacts on American Midwest hydrology

An important potential consequence of climate change is the modification of the water cycle in agricultural areas, such as the American Midwest. Soil moisture is the integrand of the water cycle, reflecting dynamics of precipitation, evapotranspiration, and runoff in space and time, and a key determinant of yield. Here we present projected changes in the hydrologic cycle over a representative area of the American Midwest from regional climate model experiments that sample a range of model configurations. While significant summer soil moisture drying is predicted in some ensemble members, others predict soil moisture wetting, with the sign of soil moisture response strongly influenced by choice of boundary conditions. To resolve the contradictory predictions of soil moisture across ensemble members, we assess an extensive and unique observational data set of the water budget in Illinois. No statistically significant monotonic trends are found in observed soil moisture, precipitation, streamflow, groundwater level, or 2 m air temperature over a recent 26 year period (soil moisture 25 years). Based on this analysis of model simulations and observations, we conclude that the sign of climate change impacts on the regional hydrology of the American Midwest remains uncertain.

Understanding the Impact of Boundary and Initial Condition Errors on the Solution to a Thermal Diffusivity Inverse Problem

In this work, we consider simulation of heat flow in the shallow subsurface. As sunlight heats up the surface of soil, the thermal energy received dissipates dow nward into the ground. This process can be modeled using a partial differential equation known a s the heat equation. The spatial distribution of soil thermal conductivities is a key factor in the mode ling process. Prior to this study, temperature profiles were recorded at different dep ths at various times. This work is motivated by trying to match these temperature profiles using a simulation-b ased approach and analytic approaches in the context of an inverse problem. Spe cifically we determine soil thermal conductivities using derivative-free optimization to minimize the non li ear-least square errors between simulation and data profile. Here, we conduct two se s f studies, assuming homogeneous and heterogeneous soil envrionments respectively. W e also study how errors in the initial and boundary conditions propagate over time using both a numerica l approach and an analytical method.