DETERMINING HOW SEA LEVEL RISE IMPACTS SALTWATER INTRUSION WITHIN COASTAL GROUNDWATER FLOW SYSTEMS IN THE SOUTHERN CHESAPEAKE BAY

A number of coastal cities located at the mouth of the Chesapeake Bay, such as Norfolk and Virginia Beach, are currently experiencing the effects of sea level rise (SLR) and are projected to continue to do so in years to come (e.g., Ezer & Corbet, 2013). Many factors affect relative SLR in the Mid-Atlantic coastal region, including eustatic sea level changes, postglacial rebound, groundwater withdrawal, and stratigraphic subsidence. The composite effect of these factors result in SLR rates that have been calculated to be approximately 3.4 mm/yr - 5.8 mm/yr (e.g., Boon et al., 2010), placing the Hampton Roads area on a list of the cities most affected by rising coastal waters (e.g., Sallenger et al., 2012). The focus of this study is how submarine groundwater discharge (SGD) changes in response to SLR. SGD, or groundwater flowing from the seabed, is important to quantify in hydrologic and geochemical budgets because it serves as a significant pathway for contaminants and nutrients. This investigation considers three different settings often found in coastal regions, which are modeled after the relatively high-relief southern Eastern Shore, riverine-dominated Elizabeth River, and wave-dominated Ocean View beach. Each of these three settings is discretized into a two-dimensional model domain as a shore-perpendicular transect approximately 5 km in the horizontal direction and 10 m - 100 m in vertical extent. SGD is calculated using SUTRA, a USGS hybrid finite-element and finite-difference numerical model that solves the coupled equations for flow and solute transport assuming that density varies with space throughout the model. For each of the three locations, a series of models are run with sea level boundary conditions reflecting different SLR scenarios extending to the year 2100. Results simulate the location of the freshwater/saltwater interface in the three coastal environments at various future times. This forward modeling approach provides a framework to better understand how different rates of SLR can impact current groundwater flow systems within a specific area of interest.