Mountains are natural dams that impede atmospheric moisture transport and water towers that cool, condense, and store precipitation. They are essential in the western United States where precipitation is seasonal, and snowpack is needed to meet water demand. With anthropogenic climate change increasingly threatening mountain snowpack, there is a pressing need to better understand the driving climatological processes. However, the coarse resolution typical of modern global climate models renders them largely insufficient for this task, and signals a need for an advanced strategy. This paper continues the assessment of variable‐resolution in the Community Earth System Model (VR‐CESM) in modeling mountain hydroclimatology to understand the role of grid‐spacing at 55, 28, 14, and 7 km and microphysics, specifically the Morrison and Gettelman (, MG1, https://doi.org/10.1175/2008JCLI2105.1) scheme versus the Gettelman and Morrison (, MG2, https://doi.org/10.1175/JCLI-D-14-00102.1) scheme. Eight VR‐CESM simulations were performed from 1999 to 2015 with the F_AMIP_CAM5 component set, which couples the atmosphere‐land models and prescribes ocean data. Refining horizontal grid‐spacing from 28 to 7 km with the MG1 scheme did not improve the simulated mountain hydroclimatology. Substantial improvements occurred with the use of MG2 at grid‐spacings ≤28 km compared to MG1 as shown with subsequent statistics. Average SWE bias diminished by 9.4X, 4.9X, and 3.5X from 55 to 7 km. The range in minimum (maximum) DJF spatial correlations increased by 0.1–0.2 in both precipitation and SWE. Mountain windward/leeward distributions and elevation profiles improved across hydroclimate variables, however not always with model resolution alone. Disconcertingly, all VR‐CESM simulations exhibited a systemic mountain cold bias that worsened with elevation and will require further examination.

Major components of hydrologic and elemental cycles reside underground, where their complex dynamics and linkages to surface waters are obscure. We delineated seasonal subsurface flow and transport dynamics along a hillslope in the Rocky Mountains (USA), where precipitation occurs primarily as winter snow and drainage discharges into the East River, a tributary of the Gunnison River. Hydraulic and geochemical measurements down to 10 m below ground surface supported application of transmissivity feedback of snowmelt to describe subsurface flow and transport through three zones: soil, weathering shale, and saturated fractured shale. Groundwater flow is predicted to depths of at least 176 m, although a shallower limit exists if hillslope‐scale hydraulic conductivities are higher than our local measurements. Snowmelt during the high snowpack water year 2017 sustained flow along the weathering zone and downslope within the soil, while negligible downslope flow occurred along the soil during the low snowpack water year 2018. We introduce subsurface concentration‐discharge (C‐Q) relations for explaining hillslope contributions to C‐Q observed in rivers and demonstrate their calculations based on transmissivity fluxes and measured pore water specific conductance and dissolved organic carbon. The specific conductance data show that major ions in the hillslope pore waters, primarily from the weathering and fractured shale, are about six times more concentrated than in the river, indicating hillslope solute loads are disproportionately high, while flow from this site and similar regions are relatively smaller. This methodology is applicable in different representative environments within snow‐dominated watersheds for linking their subsurface exports to surface waters.

Abstract. Long-distance seawater intrusion has been widely observed through the subsurface conduit system in coastal karst aquifers as a source of groundwater contaminant. In this study, seawater intrusion in a dual-permeability karst aquifer with conduit networks is studied by the two-dimensional density-dependent flow and transport SEAWAT model. Local and global sensitivity analyses are used to evaluate the impacts of boundary conditions and hydrological characteristics on modeling seawater intrusion in a karst aquifer, including hydraulic conductivity, effective porosity, specific storage, and dispersivity of the conduit network and of the porous medium. The local sensitivity analysis evaluates the parameters' sensitivities for modeling seawater intrusion, specifically in the Woodville Karst Plain (WKP). A more comprehensive interpretation of parameter sensitivities, including the nonlinear relationship between simulations and parameters, and/or parameter interactions, is addressed in the global sensitivity analysis. The conduit parameters and boundary conditions are important to the simulations in the porous medium because of the dynamical exchanges between the two systems. The sensitivity study indicates that salinity and head simulations in the karst features, such as the conduit system and submarine springs, are critical for understanding seawater intrusion in a coastal karst aquifer. The evaluation of hydraulic conductivity sensitivity in the continuum SEA-WAT model may be biased since the conduit flow velocity is not accurately calculated by Darcy's equation as a function of head difference and hydraulic conductivity. In addition, dispersivity is no longer an important parameter in an advectiondominated karst aquifer with a conduit system, compared to the sensitivity results in a porous medium aquifer. In the end, the extents of seawater intrusion are quantitatively evaluated and measured under different scenarios with the variabilities of important parameters identified from sensitivity results, including salinity at the submarine spring with rainfall recharge, sea level rise, and a longer simulation time under an extended low rainfall period.

No abstract

A hybrid discrete‐continuum numerical model, Variable‐Density Flow and Solute Transport—Conduit Flow Process (VDFST‐CFP), is developed to simulate seawater intrusion to a coastal karst aquifer with a conduit network. The Darcy‐Weisbach equation is applied to simulate the nonlaminar groundwater flow in the conduit system that is conceptualized as pipes, while the Darcy equation is used for laminar groundwater flow in the continuum porous medium. Density‐dependent groundwater flow with appropriate additional density terms in the conduit is analytically derived. The flow and transport equations are coupled, and numerically solved by the finite difference method with an implicit iteration procedure. Two synthetic benchmarks are developed to compare the VDFST‐CFP model results with other numerical models, such as the variable‐density SEAWAT, constant‐density continuum MODFLOW/MT3DMS, and constant‐density discrete‐continuum CFPv2/UMT3D models. The VDFST‐CFP model compares reasonably well with the other model results in both conduit and porous medium domains, and well describes water and salt exchange between the two systems. Under turbulent flow conditions within the conduit, the Darcy‐Weisbach equation calculates the flow rate more accurately without overestimation by the Darcy equation. Sensitivity analysis indicates that conduit diameter, friction factor, matrix hydraulic conductivity, and effective medium porosity are important parameters in the VDFST‐CFP model. The pros and cons of the VDFST‐CFP model are discussed, including the model assumptions and simplifications, limitations of the discrete‐continuum modeling method, and the convergence criteria. In general, the newly developed VDFST‐CFP model provides a new numerical modeling method for simulating seawater intrusion in a coastal karst aquifer with conduits.

In this study, a groundwater flow cycling in a karst springshed and an interaction between two springs, Spring Creek Springs and Wakulla Springs, through a subground conduit network are numerically simulated using CFPv2, the latest research version of MODFLOW-CFP (Conduit Flow Process). The Spring Creek Springs and Wakulla Springs, located in a marine estuary and 11 miles inland, respectively, are two major groundwater discharge spots in the Woodville Karst Plain (WKP), North Florida, USA. A three-phase conceptual model of groundwater flow cycling between the two springs and surface water recharge from a major surface creek (Lost Creek) was proposed in various rainfall conditions. A high permeable subground karst conduit network connecting the two springs was found by tracer tests and cave diving. Flow rate of discharge, salinity, sea level and tide height at Spring Creek Springs could significantly affect groundwater discharge and water stage at Wakulla Springs simultaneously. Based on the conceptual model, a numerical hybrid discrete-continuum groundwater flow model is developed using CFPv2 and calibrated by field measurements. Non-laminar flows in conduits and flow exchange between conduits and porous medium are implemented in the hybrid coupling numerical model. Time-variable salinity and equivalent freshwater head boundary conditions at the submarine spring as well as changing recharges have significant impacts on seawater/freshwater interaction and springs' discharges. The developed numerical model is used to simulate the dynamic hydrological process and quantitatively represent the three-phase conceptual model from June 2007 to June 2010. Simulated results of two springs' discharges match reasonably well to measurements with correlation coefficients 0.891 and 0.866 at Spring Creeks Springs and Wakulla Springs, respectively. The impacts of sea level rise on regional groundwater flow field and relationship between the inland springs and submarine springs are evaluated as well in this study.

Dynamical downscaling is a widely used technique to properly capture regional surface heterogeneities that shape the local hydroclimatology. However, in the context of dynamical downscaling, the impacts on simulation fidelity have not been comprehensively evaluated across many user-specified factors, including the refinements of model horizontal resolution, large-scale forcing datasets, and dynamical cores. Two global-to-regional downscaling methods are used to assess these: specifically, the variable-resolution Community Earth System Model (VR-CESM) and the Weather Research and Forecasting (WRF) Model with horizontal resolutions of 28, 14, and 7 km. The modeling strategies are assessed by comparing the VR-CESM and WRF simulations with consistent physical parameterizations and grid domains. Two groups of WRF Models are driven by either the NCEP reanalysis dataset (WRF_NCEP) or VR-CESM7 results (WRF_VRCESM) to evaluate the effects of large-scale forcing datasets. The simulated hydroclimatologies are compared with reference datasets for key properties including total precipitation, snow cover, snow water equivalent (SWE), and surface temperature. The large-scale forcing datasets are critical to the WRF simulations of total precipitation but not surface temperature, controlled by the wind field and atmospheric moisture transport at the ocean boundary. No significant benefit is found in the regional average simulated hydroclimatology by increasing horizontal resolution refinement from 28 to 7 km, probably due to the systematic biases from the diagnostic treatment of rainfall and snowfall in the microphysics scheme. The choice of dynamical core has little impact on total precipitation but significantly determines simulated surface temperature, which is affected by the snow-albedo feedback in winter and soil moisture estimations in summer.

Five periods of increased electrical conductivity have been found in the karst conduits supplying one of the largest first magnitude springs in Florida with water. Numerous well-developed conduit networks are distributed in the Woodville Karst Plain (WKP), Florida and connected to the Gulf of Mexico. A composite analysis of precipitation and electrical conductivity data provides strong evidence that the increases in conductivity are directly tied to seawater intrusion moving inland and traveling 11 miles against the prevailing regional hydraulic gradient from from Spring Creek Spring Complex (SCSC), a group of submarine springs at the Gulf Coast. A geochemical analysis of samples from the spring vent rules out anthropogenic contamination and upwelling regional recharge from the deep aquifer as sources of the rising conductivity. The interpretation is supported by the conceptual model established by prior researchers working to characterize the study area. This paper documents the first and longest case of seawater intrusion in the WKP, and also indicates significant possibility of seawater contamination through subsurface conduit networks in a coastal karst aquifer.