
Coastal Flooding Risk Assessment for Prince William County
Occoquan River, Occoquan Bay, Marumsco Creek - Prince William County, Virginia
Flood Hazards analysis conducted by:
George Mason University, Department of Civil, Environmental and Infrastructure Engineering (CEIE) 445/645 Flood Hazards Engineering students
for
Prince William County, Department of Emergency Services
CLIENT REQUEST
Photo of the Town of Occoquan from the Great Flood of 2008 https://historicoccoquan.com/history/2008/GreatFlood.htm
The Prince William County (PWC) Office of Emergency Management is seeking storm surge scenario data to validate and support evacuation planning assumptions for the coastal area. Specific areas to be studied are the Town of Occoquan, Veterans Park, and coastal region that sits between these two locations. Their project request asks that riverine, coastal, and sea level rise flooding hazards be incorporated to create a compound flood model. This project is intended to advise evacuation plan development, our task is to provide models that demonstrate which populated areas are at risk of flooding. The PWC Office of Emergency Management can assess the risk to critical building and transportation infrastructure in the region using the provided models and GIS Files.
PWC Emergency Management Coordinator Brian Meisner and Deputy Emergency Management Coordinator Katie Kitzmiller provided feedback on some specific expectations for this model- the county is interested in models based on storms with moderate return periods and longer lasting storm durations (as opposed to flashy storm durations).
The primary objective of this flood hazard project is to generate scenario data that will validate and support evacuation planning efforts for Prince William County, a region increasingly vulnerable to coastal flooding. The project aims to provide comprehensive models using HEC-RAS and ArcGIS software, addressing the specific areas of concern identified by the client. In addition to the client’s request for models focused on individual flooding hazards, our team proposes to develop compound flooding models that integrate both coastal and pluvial hazards for a more robust analysis. Different scenarios with varying return periods and intensities will offer the client a wider perspective of the risks that the region faces. This data will be invaluable for assessing the need for flood evacuation routes, thus enhancing overall preparedness in the region for inevitable flooding events.
BACKGROUND
The Town of Occoquan currently faces challenges with coastal flooding, and the region anticipates increased future flooding due to rising sea levels and intensified storms. Historically, the Occoquan and the Woodbridge areas have experienced inundation during hurricanes, particularly from coastal flooding:
1972
Combined flooding from Hurricane Agnes destroys businesses, homes, and critical infrastructure in the county, including the Route 123 bridge.
2003
Hurricane Isabel , one of the costliest and most severe hurricane in the recorded history of Virginia.
2008
Heavy rainfall across the region caused significant pluvial and fluvial flooding in the county.
Aside from these extreme weather events, smaller storms frequently result in riverine flooding along the banks of the Occoquan River. It is anticipated that sea level rise will exacerbate flooding in the town, making rain events more likely to flood surrounding areas. There is a concern that the Town of Occoquan, Veteran’s Park, and Bay St. are at risk to coastal storm surge and sea level rise. Our client, the Prince William County Office of Emergency Management, has also observed recent flooding occurring after rainfall events.
The study area encompasses approximately six miles of river shoreline that stretches from the Occoquan Dam to the Neabsco Creek railroad crossing.
The first area that the client has requested to be studied is the Town of Occoquan between River Mill Park and Ox Road Bridge (Virginia Route 123).
Observing the current Federal Emergency Management Agency (FEMA) Flood Hazard Zones available via the agency's National Flood Hazards Layer (NFHL) , which model the existing and future 100 and 500 yr flood plain(s) (1% and 0.2% annual chance respectively) for the National Flood Insurance Program (NFIP) , a significant area of the waterfront town is already in a known regulatory flood plain.
The second area is located further south and studies the effects of compound flooding on Veterans Park from Marumsco Creek.
Observing the existing FEMA NFHL, a major concern aside from structures and residences being directly within the flood plain is the multiple wash outs and wetlands that cross Veteran's Drive, which can lead to people becoming trapped if not evacuated prior to a significant flooding event.
The final area examines the impact of flood hazards on coastal residents on Bay Street, also located near Marumsco Creek, which is subject coastal, and pluvial flooding. A number of residences in Bay Street are completely within the NFHL flood plain in addition to being potentially trapped due to flooding along Veteran's Drive.
Additionally, although not directly relevant to this project, both Veterans Park and Bay Street areas are located adjacent to the Occoquan Bay National Wildlife Refuge (Former Woodbridge Research Facility) (FWRF), an EPA superfund site .
Existing Federal Emergency Management Agency (FEMA) coastal and pluvial Flood Hazard zones (National Flood Hazards Layer) (NFHL).
MODELING
To prepare Virginia for the every-growing risk of destructive storms and sea level rise due to climate change, the Virginia Department of Recreation and Conservation (DCR) adopted the Coastal Resilience Master Plan (CRMP) in 2021. As a part of the CRMP, the DCR produced flood models for 1,830 sub-watersheds adjacent to Virginia’s eastern coast. Two of these sub-watersheds are used to model the region.
As apart of their models, the DCR includes ...
- Terrain (.hdf, .topo, .vrt): 10ft-by-10ft resolution of the elevation profile, one for each of the two sub-watersheds.
- Geometry (.hdf): Includes a mesh that describes the underlying terrain. Breaklines are included to provide additional resolution to critical areas like waterways and major roads. The geometry is pre-set with coastal and pluvial boundaries that allow for time-dependent data like precipitation to be put into the model. One geometry exists for each sub-watershed.
- Curve Number (.hdf, .tif): 1ft-by-1ft resolution infiltration data from the 2022 Chesapeake Land Use and Land Cover (LULC) Database to model what proportion of water becomes runoff. One curve number file exists for each sub-watershed.
- Manning’s N (.hdf., .tif): 30m-by-30m resolution friction data from the 2019 National Land Cover Database (NLCD) to model how quickly runoff occurs. Each sub-watershed region has a Manning's N file that pertains to it.
- Hyetographs (.dss): Includes 63 precipitation plans. 2-hour, 6-hour, and 24-hour storm durations with varying whole number intensities (2in, 3in, 4in, etc). Derived from NOAA Atlas 14 data. Various scenarios were applied to each of the sub-watersheds.
These DCR models served as the basis for this project’s models. The flood modeling was performed with HEC-RAS 6.5, a flood modeling application produced by the US Army Corps of Engineers. Two separate DCR sub-watersheds captured the vast majority of the area specified by the client, leaving only a portion of the Occoquan Bay National Wildlife Refuge unanalyzed.
The meat of the modeling is the boundary conditions; this is what permits the input of riverine, coastal, and pluvial data into the model to be applied to the five aforementioned model components (terrain, geometry, etc.).
Coastal elevation in feet over 24 hours
- Riverine: Due to a lack of access to Occoquan River data, it was decided to represent the riverine boundary condition with a baseflow of 0ft because water is released from the Occoquan Dam only on occasions.
- Coastal: For 24-hour storms, the coastal elevation over time used a sine function such that the maximum was the MHW (0.87ft), the minimum was the MLW (-0.72ft), and two full tide cycles were achieved in the 24-hour period. An elevation was calculated for and applied to the models every five minutes of the 24 hour period using this sine formula. To model higher elevation patterns, a constant value was added to each 5-minute elevation. The MHW and MLW are from the NOAA Tides and Currents Dahlgren station, taken October 2024. For 6-hour storms, a constant whole number coastal elevation was applied.
Pluvial: The precipitation vs time data was taken directly from the DCR hyetographs
To decide which storm scenario to provide to the client, a scenario table (shown below) was established based off the client’s request for storm scenarios with moderate rainfall intensity, a 6-hr or 24-hr storm duration, and a high coastal elevation. The group agreed that due to being restricted to uploading a single sub-watershed per map, more than two models would be too much to succinctly present to the client in this report.
Scenario Table
To understand the models, it is important to understand each part of the scenario table...
Model Properties
- Storm Duration: Simply, it is the duration of the storm in the model. 6-hr and 24-hr were chosen over other storm durations because these are the durations that the DRC precipitation datasets last for.
- Rainfall Intensity: This is the collective depth of rainfall during the entire storm durations. The rainfall intensity at any exact moment within the storm duration varies just as it does in a real storm.
- Coastal Elevation - NAVD88: All coastal water elevations in the model are with reference to the NAVD-88 datum. For 6-hr storms, we chose a static coastal boundary, meaning whatever elevation "y" is equal to in the table is what the coastal elevation was modeled as. For 24-hr storms, the aforementioned non-static, sinusoidal coastal boundary was imployed. Whatever "y" is equal to in the table is the maximum elevation of the coastal boundary throughout the storm duration. Because the maximum elevation is relative to MHW (0.87ft) and the lowest elevation is relative to MLW (-0.72ft), the lowest 24-hr coastal boundary gets is y-0.87ft-0.72ft. For example, a scenario with a coastal boundary of non-static y=3ft, the highest the coastal boundary gets is 3ft, and the lowest the coastal boundary gets is 1.41ft.
Design Scenario
Location of NOAA Atlas 14 intensity, duration, and frequency data used for the project
Design storms are the "once every 100 years" storms. A design storm is defined by its rainfall intensity, storm duration, and frequency (25 years, 50 years, etc.). Durations and frequencies are standardized; it is the intensity for any specific duration and frequency that varies from region to region. An example from the scenario table: a 24-hr, 25-yr storm in the Town of Occoquan has an intensity of 6.04in. In our models, a 24-hr, 6in storm is 0.04 inches short from replicating a 25-yr storm. Therefore, the first highlighted scenario in the scenario table can be referenced as being approximately a 25-yr storm. Which models that were chosen for this presentation was also motivated by how closely they resembled specific design scenarios.
NOAA Atlas 14 intensity, duration, and frequency data for the Town of Occoquan https://hdsc.nws.noaa.gov/pfds/pfds_map_cont.html?bkmrk=va
It is important to note that no considerations for existing stormwater systems in the Town of Occoquan were incorporated into our models. This likely results in an overestimate in the flooding that the Town of Occoquan region experiences. Ballywhack Branch Creek, a stream that approaches the town from the west, has a noticeable impact on the flooding that the town experiences. Beginning with the team's earliest model runs, it was noted that Ballywhack Creek was the Town of Occoquan's main source of excessive flooding. As depicted in the results, Union Street suffers the most from Ballywhack overflowing.
RESULTS
Flood Model Simulation Results: Scenario 1
This visual representation illustrates the impact of Scenario 1 on buildings and roadways. Scenario 1 simulates a storm with a 24-hour duration and a 25-year return period, combined with a 2-foot non-static coastal elevation. This scenario was selected to align with our client's expectations for mid-range climate events.
Key Findings:
- Affected Buildings: 6,833 buildings intersect with the flood extents across the entire model domain.
- AOI 1 has more buildings affected than AOI 2 and 3
An important note about our models:
Due to the nature of the DCR Terrain and Geometry files, some man-made infrastructure projects, like bridges and culverts, are not incorporated. The most noticeable example of this is the Ox Road Bridge spanning the Occoquan River. Although the models show this bridge to be flooded, this actually just represents water flowing under the bridge.
Building Intersections - 24-hour duration, 25-year frequency storm +2 ft non static coastal elevation - AOI 1
Building Intersections - 24-hour duration, 25-year frequency storm +2 ft non static coastal elevation - AOI 2 &3
- Roadways: We have assumed that any flooding depth greater than 0.5 feet can cause roadblocks.
- AOI 1 (Historic Town of Occoquan): Parts of Union Street, Mill Street, Commerce Street, Polar Alley and Center Lane seems to be experiencing roadblocks at portions highlighted in red.
Roadblocks - 24-hour duration, 25-year frequency storm +2 ft non static coastal elevation - AOI 1
- AOI 2 (Veteran's Park): This area experiences fewer roadblocks on the park's frontage. One significant portion is on Veteran's drive south of the park. However more roadblocks are visible on the west side of the AOI.
- AOI 3 (Bay Street): No Significant Roadblock in this area
Roadblocks - 24-hour duration, 25-year frequency storm +2 ft non static coastal elevation - AOI 2 AND 3
Please note, this simulation does not account for the existing storm drain network and represents an exaggerated version of the actual conditions.
We have also included a second Scenario; a 24-hour storm with a 50-year return period in conjunction with 3 ft non-static coastal elevation. We included this scenario because we wanted to see a more severe case, keeping in mind that the area has previously encountered such heavy storms.
Building Intersections - 24-hour duration, 50-year frequency storm +3 ft non static coastal elevation - AOI 1
Building Intersections - 24-hour duration, 50-year frequency storm +3 ft non static coastal elevation - AOI 2,3
Roadblocks - 24-hour duration, 50-year frequency storm +3 ft non static coastal elevation - AOI 1
Roadblocks - 24-hour duration, 50-year frequency storm +3 ft non static coastal elevation - AOI 2,3
The result from this scenario follows the same pattern as the previous one, with slightly bigger impact.
Inter-Scenario Comparison for Roadblocks - AOI 1 (Scenario 2 on the right)
CONCLUSION
Flooding in the Town of Occoquan (photo: https://www.dcr.virginia.gov/dam-safety-and-floodplains/document/cfpf/App-Studytoevaluatestormwaterandfloodresilience.pdf )
Historically, the Town of Occoquan, Veterans Park, and the coastal area in between have been subject to major flooding from coastal and pluvial hazards. With this known risk, the Prince William County Office of Emergency Management has requested flood models for moderate storms to better anticipate infrastructure damage and plan evacuation routes.
Our models use a combination of pluvial and coastal conditions, which model current 25-50 year storms with varying static coastal elevations, show flooding through the County’s identified Areas of Interest. The results obtained from these simulations can now be used to plan emergency evacuation maps as roadblocks and affected buildings have been identified for each storm.
The Featherstone Neighborhood along the banks of Occoquan Bay (photo: https://www.neighborhoods.com/featherstone-woodbridge-va )
According to the results, out of the three study areas, the Town of Occoquan experienced the most infrastructure damage due to flooding since it is the most densely populated/developed area. However, during a 50-year storm, the roads and developed areas surrounding Veterans Park were also heavily impacted.
In the future, the accuracy of our model can be further enhanced by including the existing storm drain infrastructure in our model. Other models may also consider riverine flooding downstream of the Occoquan Dam, or take into account a non-static coastal boundary for a more robust simulation.