Analytical and Numerical Modeling of Solute Intrusion,Recovery, and Rebound in Fractured Bedrock

Contaminated groundwater in fractured bedrock can expose ecosystems to undesired levels of risk for extended periods due to prolonged back-diffusion from rock matrix to permeable fractures. Therefore, it is key to characterize the diffusive mass loading (intrusion) of contaminants into the rock matrix for successful management of contaminated bedrock sites. Even the most detailed site characterization techniques often fail to delineate contamination in rock matrix. This study presents a set of analytical solutions to estimate diffusive mass intrusion into matrix blocks, it is recovered by pumping and concentration rebound when pumping ceases. The analytical models were validated by comparing the results with (1) numerical model results using the same model parameters and (2) observed chloride mass recovery, rebound concentration, and concentration in pumped groundwater at a highly fractured bedrock site in Alberta, Canada. It is also demonstrated that the analytical solutions can be used to estimate the total mass stored in the fractured bedrock prior to any remediation thereby providing insights into site contamination history. The predictive results of the analytical models clearly show that successful remediation by pumping depends largely on diffusive intrusion period. The results of initial mass from the analytical model was used to successfully calibrate a three-dimensional discrete fracture network numerical model further highlighting the utility of the simple analytical solutions in supplementing the more detailed site numerical modeling. Overall, the study shows the utility of simple analytical methods to support long-term management of a contaminated fractured bedrock site including site investigations and complex numerical modeling.     

Transient and Transition Factors in Modeling Permafrost Thaw and Groundwater Flow

Permafrost covers approximately 24% of the Northern Hemisphere, and much of it is degrading, which causes infrastructure failures and ecosystem transitions. Understanding groundwater and heat flow processes in permafrost environments is challenging due to spatially and temporarily varying hydraulic connections between water above and below the near-surface discontinuous frozen zone. To characterize the transitional period of permafrost degradation, a three-dimensional model of a permafrost plateau that includes the supra-permafrost zone and surrounding wetlands was developed. The model is based on the Scotty Creek basin in the Northwest Territories, Canada. FEFLOW groundwater flow and heat transport modeling software is used in conjunction with the piFreeze plug-in, to account for phase changes between ice and water. The Simultaneous Heat and Water (SHAW) flow model is used to calculate ground temperatures and surface water balance, which are then used as FEFLOW boundary conditions. As simulating actual permafrost evolution would require hundreds of years of climate variations over an evolving landscape, whose geomorphic features are unknown, methodologies for developing permafrost initial conditions for transient simulations were investigated. It was found that a model initialized with a transient spin-up methodology, that includes an unfrozen layer between the permafrost table and ground surface, yields better results than with steady-state permafrost initial conditions. This study also demonstrates the critical role that variations in land surface and permafrost table microtopography, along with talik development, play in permafrost degradation. Modeling permafrost dynamics will allow for the testing of remedial measures to stabilize permafrost in high value infrastructure environments.     

Simulation of the Collection of Catch Crops for the Recovery of Agricultural Resources using Geographic and Statistical Data

The cultivation of catch crops has great potential in agricultural areas for the recovery of nutrients, mitigation of groundwater contamination, and secondary utilization of harvestry. To estimate the potential volume of catch crops that can be collected at prospective locations, we simulated the cultivation and collection process at a regional scale using geographic and statistical data. Three types of geographic data, namely the locations of greenhouses, collecting stations and road networks, were used in a geographic information system to compute the volume of collected catch crops and carrying distance between greenhouses and collecting stations. Carbon emission from transportation of the catch crop was calculated using an improved ton‐km method, and the results were compared with the carbon content of the catch crop to evaluate the carbon balance. We found that the total fresh weight of the collected catch crops was 67900 t, of which 70% was collected at the top 15 of 73 stations. Carbon emissions from transportation ranged from >1 to 12% of the carbon content of the catch crop. The analytical method used in this study can readily be applied to other environmental studies concerning the collection of agricultural products and other biomass resources.     

Integrated Surface and Subsurface Hydrological Modeling with Snowmelt and Pore Water Freeze–Thaw

For the simulation of winter hydrological processes a gap in the availability of flow models existed: one either had the choice between (1) physically-based and fully-integrated, but computationally very intensive, or (2) simplified and compartamentalized, but computationally less expensive, simulators. To bridge this gap, we here present the integration of a computationally efficient representation of winter hydrological processes (snowfall, snow accumulation, snowmelt, pore water freeze–thaw) in a fully-integrated surface water-groundwater flow model. This allows the efficient simulation of catchment-scale hydrological processes in locations significantly influenced by winter processes. Snow accumulation and snowmelt are based on the degree-day method and pore water freeze–thaw is calculated with a vertical heat conduction approach. This representation of winter hydrological processes is integrated into the fully-coupled surface water-groundwater flow model HydroGeoSphere. A benchmark for pore water freeze–thaw as well as two illustrative examples are provided.     

A Review of Coastal Fog Microphysics During C-FOG

Boundary-Layer Meteorology - Our goal is to provide an overview of the microphysical measurements made during the C-FOG (Toward Improving Coastal Fog Prediction) field project. In addition, we...     

Measuring saturated hydraulic conductivity and anisotropy of peat by a modified split-container method

Long-term changes in the physical and hydraulic properties of peat result from the decomposition and consolidation processes. The saturated hydraulic conductivity (K _s) of peat is depth-dependent and could increase or decrease with depth; therefore, K _s determination on a large number of samples is required to more accurately assess field variability. The cube method is a popular laboratory procedure to determine K _s because it allows use of smaller sample dimensions, while minimizing edge effects. This article describes the design and use of an alternate split container to enclose undisturbed peat cubes during measurement of K _s and the application of this method to a field site in Northwest Territories, Canada. The method allows for simpler and more controlled K _s measurements while also permitting anisotropy measurements. Matched tests on identical samples (K _s range 21–314 m/d), using the split-container and wax methods, showed a good agreement for intermediate K _s values; however, significant deviations occurred for low and high K _s values. Hydraulic-conductivity measurements taken on field samples showed a decrease with depth and exhibited anisotropic hydraulic conductivity as expected based on previous studies.     

Characterizing Atmospheric Aerosols off the Atlantic Canadian Coast During C-FOG

Boundary-Layer Meteorology - Marine aerosols play an important role in the Earth’s climate, but their effects remain highly uncertain due to a poor understanding of their sources, properties,...