Prioritizing the protection and creation of natural and nature-based features for coastal resilience using a GIS-based ranking framework – an exportable approach

Increasing the preservation and creation of natural and nature-based features (NNBF), like wetlands, living shorelines, beaches, dunes and other natural features to improve community resilience in the face of increasing coastal flooding may be achieved by highlighting the locally relevant benefits that these features can provide. Here we present a novel application of the least-cost geospatial modeling approach to generate inundation pathways that highlight landscape connections between NNBF and vulnerable infrastructure. Inundation pathways are then used to inform a ranking framework that assesses NNBF based on their provision of benefits and services to vulnerable infrastructure and for the broader community including 1) the flooding mitigation potential of NNBF, 2) the relative impact of those NNBF on local infrastructure, and 3) co-benefits for the broader community linked to incentive programs like nutrient reduction crediting and the Federal Emergency Management Agency’s Community Rating System. Inundation pathways are also used to identify locations lacking in benefits from NNBF as target areas for NNBF restoration or creation. This approach, applied here for coastal Virginia, with project outputs available via an interactive map viewer1, can be customized for application in any community to identify high-priority NNBF that are particularly beneficial for preservation and to identify target areas for new or restored features.


General
Coastal flooding resulting from high tides and storm surge, exacerbated by accelerating rates of sea-level rise (SLR), causes significant impacts on infrastructure of coastal communities in Virginia (Ezer, 2018) and globally (Vitousek et al., 2017). In fact, coastal Virginia SLR rates are some of the highest in the U.S. (5.4mm/ yr) surpassed only by locales on the U.S. Gulf Coast (Boon et al., 2018). Preserving existing and/or creating new natural and naturebased features (NNBF) are increasingly recognized as beneficial and cost-effective solutions for coastal management (Temmerman et al., 2013;Bridges et al., 2015;Bilkovic et al., 2016;Reguero et al., 2018;Sutton-Grier et al., 2018;Powell et al., 2019;Bridges et al., 2021).
There are state and federal policies and incentives designed to support the preservation and restoration of NNBF in the U.S., and specifically living shorelines, as a preferred management strategy for shoreline erosion control in appropriate settings. However, there remains in Virginia the perception that NNBF are not as effective at protecting shorelines as more traditional shoreline armoring approaches (Center for Coastal Resources Management, Bilkovic et al., 2016). NNBF are being lost throughout coastal zones, in part due to the lack of detailed and local actionable information about the many services and benefits that these features can provide (Arkema et al., 2017). It is suggested that providing targeted information can help bolster local decisions that preserve and create new NNBF (Center for Coastal Resources Management, Spalding et al., 2014b;Restore America's Estuaries, 2015;Sutton-Grier et al., 2015;Arkema et al., 2017).
Here we present an adaptable and generalized framework, informed by input from community partners that ranks existing NNBF and enables identification of areas where new or restored NNBF would provide maximum benefits for community resilience to coastal flooding. With the goal of providing actional information for local decision-making via an interactive map viewer 1 , this framework utilizes a geographic information system (GIS) approach that links NNBF with buildings using high resolution elevation data to create inundation pathways. While benefit frameworks can be modified to evaluate the provision of different combinations of services (VanZomeren and Acevedo-Mackey, 2019), our approach focuses on two primary benefits: 1) mitigating coastal flooding impacts on built infrastructure, and 2) allowing communities to take advantage of Federal and state incentives for communities to reduce flood risk and improve water quality, specifically the Federal Emergency Management Agency's (FEMA) National Flood Insurance Program (NFIP) Community Rating System (CRS) (hereafter CRS) Credit Program, and the water quality load reduction credits for approved best management practices of the Chesapeake Bay Program Total Maximum Daily Load (TMDL).

Study area
The study area encompasses coastal areas that are contiguous to the Chesapeake Bay, its tributaries, and Atlantic Ocean, and generally lie at less than 3 meters (10 feet) in elevation above a tidal shoreline (Center for Coastal Resources Management, 2019) ( Figure 1). These at-risk areas were identified based on local experience of past storms (e.g., Hurricane Isabel in 2003, Lawrence et al., 2005), and land most vulnerable to coastal flooding (Boon, 2012;Sallenger et al., 2012;Ezer and Atkinson, 2014;Dahl et al., 2017;Boon et al., 2018;Ezer, 2019).
In coastal Virginia, these low-lying and vulnerable areas comprise approximately 2,829 square kilometers of land, approximately 5% of the total area of Virginia (110,784 square kilometers). Within these areas are more than 22,100 kilometers of tidal shoreline, more than 170,000 primary buildings, and all or portions of 45 counties and cities, based on our analysis. For reference, the coastal Virginia localities encompassing the study area are home to 6 million people, or 70% of the state's population, and provide 78% of the gross domestic product (Virginia Coastal Resilience Master Plan, 2021). The coastal land ownership includes significant areas of public ownership (22%), including large military facilities, with the remaining lands in private ownership (78%) (Virginia Natural Heritage Program, 2018).
Virginia's coastal lands are split into five regions by the Chesapeake Bay and its tributaries: Northern Neck, Middle Peninsula, Peninsula, Southern Tidewater, and the Eastern Shore ( Figure 1). It consists of a heterogeneous landscape of low-lying and steeply graded topography, extensive tidal marshes and beach and dune systems, and rural and urban development patterns. These regions, commonly known as Tidewater Virginia, include areas with steep grades and high banks, and areas typified by large tidal rivers, low-relief, and low banks (Ludwig, 2016). Development patterns are also heterogeneous, with a mix of very urbanized and rural localities. Urban areas, as defined by the U.S. Census Bureau (Virginia Department of Health, 2006), include those of the Southern Tidewater, Northern Virginia and limited areas in or near the Richmond metropolitan area (portions east of Richmond are in the study area). Rural localities are located across the Middle Peninsula, Eastern Shore, and the Northern Neck.

Datasets
Assessment of coastal resources and their benefits to infrastructure requires spatial datasets for NNBF, as well as the inputs necessary for developing IPs: source points (building footprints) and elevation data. Available datasets were utilized to capitalize on existing data and expertise while also improving the exportability of the approach for localities with different data availability.
No comprehensive mapped NNBF inventory exists for many coastal areas, including for Virginia. NNBF datasets from multiple local, state, and national sources were collected and updated. Table 1 details data sources for each NNBF type and any processing completed.
Multiple data sources were used to create the primary buildings layer with identifiers for important community infrastructure. The statewide building footprint dataset (Virginia Geographic Information Network, 2017a) was processed to reduce the inclusion of secondary structures like boathouses and freestanding garages by extracting only those footprints that exceeded a minimum area of 80 m2, a threshold determined by test runs in urban and rural locales. The statewide building footprint dataset is based on locality-supplied data, and some localities lack footprint polygons. For those localities, Virginia land cover data (Virginia Geographic Information Network, 2017b) was processed to capture building footprints. The USGS National Structures Dataset (U.S. Geological Survey, 2018) was used to add information to buildings identified as critical community facilities.
Elevation data was extracted from a 1-meter resolution digital elevation model (DEM) (Danielson et al., 2016) and used to create the study area boundary, supply elevation information for NNBF, and for the spatial modeling steps detailed in Section 2.2.

Spatial modeling
Characterizing landscape connections between vulnerable infrastructure (in this study, primary buildings) and NNBF was accomplished using a least-cost path analysis. Least-cost approaches are used to delineate the most cost-effective routes across a landscape, where cost can be a function of economic cost Map of study area with regions, and urban and rural areas denoted. (Imagery source: Esri, DigitalGlobe, GeoEye, i-cubed, USDA FSA, USGS, AEX, Getmapping, Aerogrid, IGN, IGP, swisstopo, and the GIS User Community).
(money), time, distance, or, as here, elevation. Using elevation as the cost function has been previously applied in hydrologic characterizations to identify concentrated flowpaths or preferential water flows across a landscape (Etherington, 2016;Wallace et al., 2018;Luo et al., 2019), and in Virginia for evaluating stormwater management and flows (Chesapeake Conservancy, n.d.). Using elevation values as the cost function results in a delineation of the lowest overall elevation connections from source points (in this case, building footprints) to the destination at the shoreline (the lowest elevation cells in the DEM), and were termed inundation pathways (IPs) for this study. The least cost path is best used as a qualitative rather than a quantitative measure of connectivity (Etherington, 2016), and IPs are best interpreted as highlighting areas where coastal flooding is likely to flow through on approach towards buildings. IPs are therefore locations where the presence of an NNBF is highly likely to provide flooding mitigation benefits for community infrastructure.
A series of geoprocessing steps were undertaken to create IPs using tools in ESRI's ArcMap 10.6.1 Spatial Analyst toolbox (ESRI, 2017) using ESRI's guidance for generating least-cost pathways (ESRI). Prior to running least-cost tools, the DEM was clipped to the tidal shoreline digitized from aerial imagery (Center for Coastal Resources Management, 2019) so that the lowest elevation cells of the DEM nearest to a source point coincided with the boundaries of digitized fine-scale hydrologic features like tidal channels and small coastal inlets. This addressed the loss in detail including small tidal marsh channels across relatively flat tidal marsh terrain, a known issue resulting from processing steps undertaken in least-cost analyses (Callow et al., 2007;Kenny et al., 2008).
Creation of the cost surface was accomplished by smoothing imperfections in the raster data then running the Flow Direction tool to calculate elevation differences across cells. Finally, the Cost Path tool was run iteratively for each of the more than 170,000 buildings to create the IPs. Computer processing time for these analyses was significant, so geoprocessing models were created to restrict each Cost Path analysis to a 2000-meter buffer around each building. This step decreased processing time without restricting resultant IP size. To further decrease processing times, geoprocessing was run in batches by region (e.g., for the Eastern Shore, for the Middle Peninsula, etc). Even with these adjustments, the geoprocessing still took months on dedicated computers. While this study utilized the best available tools and knowledge available at the time of analysis, for subsequent applications it is likely that using more current software (e.g., ESRI's ArcPro) or utilizing other scripting (e.g., Python) tools could reduce processing time.
The resultant IPs can be complex in shape and often indicate more than one distinct pathway for a single building ( Figure 2). Each building's collection of flow paths is considered as one IP. IPs lie along low elevation connections anywhere on the land surface ( Figure 2A) and not only directly from the closest tidal shoreline. This reflects the possibility of coastal flooding approaching from other places such as a nearby low-lying area like a wetland. Multiple building IPs often overlap in low-relief areas ( Figure 2B). IPs are a simplified  2022). Project location points were mapped and buffered based on the permitted length of the project. Spatial analysis required polygons on shore, so these were projected to the shoreline digitized in Center for Coastal Resources Management (2019). Elevations of these NNBF were defined as breakwater = 0-3 feet, marsh sill = 0-2 feet, and oyster sill = 0-2 feet, based on best professional judgment.

Scrub-Shrub
Extracted 'scrub/shrub' polygons from Virginia land cover classification data (Virginia Geographic Information Network, 2017b). To reduce misclassification errors, polygons in the water and adjacent to riprap were eliminated.
Tidal Marsh Tidal marshes from the Tidal Marsh Inventory (Center for Coastal Resources Management, 2019).

Wooded
Extracted 'forested' and 'tree' polygons from Virginia land cover classification data (Virginia Geographic Information Network, 2017b). Inspection of the land cover data revealed inconsistencies with how forest and tree areas were defined (i.e., forested areas were defined as areas of more than 1 acre but many areas in the map were coded as forest but smaller than one acre). Because data seemed not to support two separate classifications, these were combined into a 'wooded' class. To reduce other misclassification errors, polygons in the water and adjacent to riprap were eliminated.
representation of the areas with highest probability to intercept coastal flooding flow, and do not map flooding extents or represent timing of flooding.
Assessing the relative provision of benefit of an NNBF to community infrastructure is accomplished by using the Spatial Join tool to count the number of IPs that intersect an NNBF. For example, Figure 2B depicts IPs from building 2 and building 3, both of which intersect a forested wetland and a tidal marsh; each of those NNBF would have two IP intersections recorded. NNBF that intersect IPs for many buildings are identified as highly beneficial and used to score criteria 2 and 3 of the Coastal Benefit Index, described in Section 2.3. The IP approach also allows for identifying buildings that lack benefits from NNBF (i.e., IPs with 0 intersecting NNBF) which defines target areas for the creation of new NNBF, described in Section 2.4.

The coastal benefit index
The Coastal Benefit Index (CBI) was developed to highlight those NNBF that provide the most benefits for communities, so that those features can be prioritized for preservation and/or restoration at the local level. The CBI ranking approach consists of four criteria or indicators intended to evaluate the degree that an NNBF can provide flooding mitigation services for community buildings and allow communities to take advantage of incentive programs that provide financial or other benefits for the community. The four criteria evaluated in this study are (1) NNBF Flooding Mitigation Potential, (2) Coastal Building Benefit, (3) Critical Community Facility Benefit, and (4) Co-Benefits Potential, and are ranked from 1 to 3, with rank 1 representing low benefit provision and rank 3 representing high benefit provision.

FIGURE 2
Maps depicting Inundation Pathways (IPs) and coastal building footprints for a location in the Middle Peninsula. (A) shows buildings and IPs with the underlying elevation raster, illustrating that IPs follow the lowest overall elevation connections between a building and the shoreline. (B) depicts those same buildings and IPs showing how each IP/building may be complex in shape and that buildings can have overlapping IPs. This figure also displays the spatial relationship between IPs and NNBF. (Imagery source: Virginia Geographic Information Network, 2017c).

NNBF flooding mitigation potential
The potential for an NNBF to provide coastal flooding mitigation services is a function of both its intrinsic physical characteristics that can mitigate flooding impacts ('capacity'), and its elevation, which controls the frequency that the feature is exposed to flooding ('opportunity').
NNBF capacity was assigned based on three physical characteristics important for attenuating wave energy and allowing for the infiltration of flooding waters: permeability of ground soils, surface roughness, and vegetation roughness (i.e., stem density), and is shown in Table 2. A simple ranking approach to assigning relative scores for these characteristics was based on literature review and best professional judgment. For example, surface roughness values in Table 2 are informed by Manning's roughness coefficients used in hydrodynamic modelling, but do not directly relate to those values. This study is focused on assessing the relative provision of flooding mitigation actions from NNBF across a heterogeneous mix of landscape types of varying slope and soil composition, across seasons, and for a range of flooding depth scenariosall factors which change the appropriate roughness coefficients (see Ye et al., 2018 for one discussion of roughness coefficient dynamics). Thus, relative scores were used.
The opportunity for an NNBF to provide flooding mitigation services is related to its location in the landscape, and specifically its elevation in relation to coastal flooding events. Rather than evaluating flooding scenarios based on historic events (spatial data for which are often lacking, as here in Virginia) or using hydrodynamic modeling (also not commonly available), water level records from a nearby tidal gauge (Virginia Natural Heritage Program, 2018) were used to calculate the opportunity score (Table 3), which is based on the frequency of high-water events reaching a set of elevation zones.
High-water events were defined as distinct episodes of consecutive water level recordings above 0m. Each NNBF was then assigned an opportunity score based on the minimum elevation band that it occupies. Low elevation NNBF, which encounter flooding waters first and most often, are most frequently able to mitigate coastal flooding impacts, and therefore receive a higher score.
A flooding mitigation potential score was calculated for each NNBF by multiplying capacity and opportunity scores. Resultant scores ranged from 0-1 and the values were divided into three groups containing approximately equal numbers of NNBF (i.e., tertiles), to reflect the relative ranking of low, medium, or high for each individual NNBF to provide flooding mitigation services (Table 4).To assess the benefit of an NNBF to buildings in the study area, IP intersections were counted for each NNBF as shown in Figure 2B, except for dunes. Dune features were evaluated using a slightly different approach because IPs do not intersect the relatively high-elevation dune features. However, because dunes provide important flooding mitigation benefits for communities (Fernańdez-Montblanc et al., 2020), they are important to include. Dunes were scored by counting the number of IPs within 100-feet of the dune polygon.
To score this criterion, NNBF were grouped into tertiles to represent the ranking of low, medium, and high, described in Table 4.

Critical community facility benefit
Of particular importance for community resilience are critical facilities like law enforcement facilities, medical buildings, schools, and emergency response centers. These facilities often serve as community shelters and/or bases for emergency response during storm events. An NNBF that intersects an IP from a critical community facility was assigned a high score for this category, as described in Table 4.

Co-benefits potential
The fourth variable considers the potential for the presence of an NNBF to be used by local governments to take advantage of Federal and state incentive programs. This study considered two programs: 1) FEMA CRS crediting program for participating localities' floodplain management activities and 2) TMDL credits for employing Best Management Practices (BMPs) to improve water quality.
Participating in the CRS program allows localities to lower flood insurance premiums for residents by engaging in floodplain management activities. The amount of discount available is determined by the number of credits a locality receives based on particular activities they undertake, including, under Section 420, preserving natural open space and natural floodplain functions (Federal Emergency Management Agency, 2017). In Virginia, the Chesapeake Bay Preservation Act mandates the conservation and protection of a Riparian Protection Area (RPA) which consists of natural features (e.g., wetlands) that are hydrologically connected to the Chesapeake Bay and a 100-foot buffer (as this dimension is specified in law, we are using the imperial system of measurement) landward of those features. Any NNBF within the RPA buffer and the FEMA-mapped special hazard flood area has the potential to gain CRS credits for the locality (Jarbeau and Stiff, 2017). Taking advantage of NNBF to provide CRS credits for the locality means that NNBF can indirectly provide financial support for residents who own flood-vulnerable property (Stiff, 2017).
Also considered is the potential for an NNBF to provide water quality services through the reduction of sediment and nutrient pollution and thereby qualify to receive water quality credits through the Chesapeake Bay TMDL program. The Bay TMDL established nutrient and sediment targets for Bay jurisdictions and accredited more than 200 BMPs including tidal and non-tidal wetlands and riparian buffers (Chesapeake Bay Program, 2022). Thus, all vegetated NNBF (i.e., NNBF except for beaches, dunes, and hybrid living shorelines with breakwaters) were considered to have the potential to receive water quality credits for localities. NNBF were scored based on how many of the two co-benefits could be provided, as described in Table 4.

Target areas for new or restored NNBF
Using the IP analysis, it is possible to identify high-priority areas where buildings currently lack benefit from NNBF and the creation or restoration of NNBF could be particularly beneficial. Target areas were defined as a 100-foot diameter circle centered on the shoreline in locations where IPs with zero intersecting NNBF cross. The 100-foot dimension was chosen to encompass the nearshore area protected by the RPA buffer, as discussed in Section 2.3.4. In some cases, multiple zero IPs hit the shoreline within 100-feet from each other, and in these cases, target area circles were merged into one polygon (Figure 3).
While creating an NNBF along any portion of the IP would confer benefits to the associated building (or buildings), placing an NNBF on or near the shoreline maximizes the potential for the assessed benefits for the community. NNBF on or near the shoreline 1) act as a first line of defense for coastal flooding, 2) exist within the regulated and protected RPA and can therefore be used to capitalize on existing state  and Federal incentive programs as discussed above, and 3) can enhance the connectivity of valuable shoreline habitats and riparian buffers. Another practical reason for centering target areas on the shoreline was the ability to capitalize on existing guidance on suitable shoreline practices for tidal shoreline erosion control via CCRM's Shoreline Management Model (SMM). The SMM is a decisionsupport tool initially developed for and applied to Virginia coastlines, with expanded use for Maryland shorelines and for some Gulf of Mexico localities (Nunez et al., 2022). SMM management recommendations reflect Virginia policy requiring the use of NNBF approaches such as living shorelines, beaches, and vegetated buffers, where suitable, and can guide decision-makers on the potential NNBF practices that are suitable along a shoreline. Relevant information to support decision making regarding potential NNBF placement in the target area was collected and provided. This includes four primary pieces of information: 1) the number of buildings with IPs intersecting the target area to provide context for how impactful a new NNBF may be for the local community, 2) land cover within and adjacent to the target area to identify existing cover types with potential for BMP credits if converted to NNBF (e.g., turf grass), 3) existing adjacent NNBF that provide an indicator of the likelihood of success for NNBF restoration or creation, as well as identifying opportunities for filling gaps along or across shoreline habitats, and 4) SMM recommendations for shoreline within the target area that can be used to guide decisions on potential suitable NNBF practices.

NNBF and inundation pathways
Inundation pathways were delineated for more than 174,636 buildings and the NNBF they intersected. NNBF cover approximately 70% of the study area, with tidal marshes accounting for the largest extent (28%) of that ( Figure 4). Yet, more than half (63%) of NNBF did not have any intersecting IPs, and the range in number of IPs per NNBF was extreme: from 0 IP intersections to a maximum of 2,468. The NNBF with the maximum number of IP intersections occurred in a forested wetland in a heavily developed area of Virginia Beach, which illustrates the important influence of both development patterns and topography on the IP analysis. Southern Tidewater and the Peninsula, both with areas of highly urbanized development, also tend to be lower relief, and therefore wide swaths of these communities with hundreds of buildings and NNBF were included within the study area. Figure 5 summarizes NNBF by region. As expected, these mostly urban and low relief of Southern Tidewater localities had a higher frequency of NNBF with IPs (52% of the NNBF in the region) than the very rural and higher relief Northern Neck region (12% of NNBF in the region). The more rural areas of the Eastern Shore and Middle Peninsula exhibit low relief but, in contrast with urban locales, contain sparser development which results in a lower frequency of NNBF with IPs (25% and 26%, respectively). The Northern Neck is the low outlier with 12% of Target area in the Peninsula region that is centered on a building IP with zero intersecting NNBF, with existing adjacent tidal marsh and wooded NNBF. This target area has existing impervious surfaces and turf grass, is recommended for a non-structural living shoreline by the SMM, and structure enhancement suggestions include adding natural features to the existing bulkhead. (Imagery source: Virginia Geographic Information Network, 2017c). NNBF summarized by the number of IPs intersecting, per study area region. Summary of areal extent within the less than 10-ft study area, by NNBF Type. The 'Other' category includes NNBF types that constitute less than 1% of the study area, including: scrub-shrub (1,609ha), dune (52 ha), marsh sill (12 ha), breakwater (5 ha), and oyster sill (2 ha).
NNBFs with 1 or more intersecting IPs, given that it is both rural and exhibits relatively higher topographic relief, resulting in a limited study area and fewer buildings and NNBF to be considered overall.
Rankings were calculated for 360,153 existing NNBF. NNBF of all types were represented in the highest overall ranking of 'most benefits provided' (Figure 6).

Target areas for new or restored NNBF
Target areas are locations where a building (or buildings) lack NNBF benefits (i.e., building IPs have no intersecting NNBF); see Figure 3. In the study area, only 3% of buildings had zero associated NNBF (5,637 buildings). Most buildings had one or more NNBF intersecting their IP, and buildings most commonly had either 2 or 3 associated NNBF (90% of buildings) (Figure 7). From the buildings that lack NNBF benefits, 867 unique target areas were identified. These target areas encompass 104 km (341,059 feet) of shoreline and 2.9 square kilometers (725 acres) of land area. This area is likely a conservative estimate of potential new NNBF, if implemented projects extend beyond the 100-foot diameter, which is possible, especially in less developed areas of suburban and rural settings. New or restored NNBF placed in those locations have the potential to provide benefits for 18,745 buildings in coastal Virginia. Importantly, 28 of the NNBF target areas would provide additional benefits to critical community facilities.

Discussion
Application of the GIS least-cost path analysis to inform the ranking of NNBF for provision of ecosystem services is a novel application of this modeling tool. Generating a visualization of the pathways for coastal water flooding allows for observations not readily made without the IPs, such as the inundation pathways shared by hundreds of buildings in urban settings. This approach also highlights instances where IPs indicate potential flooding pathways not only from the closest or readily observable tidal waters, but from other directions or some distance away.
The fine scale used in this analysis is uniquely applicable for local planning, with NNBF impacts measured on an individual building scale, because these decisions are commonly made on a parcel-byparcel basis. NNBF targets are also at the building or parcel scale, identifying potential small-scale restoration opportunities. For example, the majority of target areas in Virginia (92%) currently include some amount of turf grass land cover, which if converted to natural cover (i.e., NNBF) would provide additional opportunity for the locality to capitalize on CRS credits and water quality incentives. When coupled with the CBI ranking for existing NNBF, the approach allows consideration of the potential cumulative benefits across the entire system that can be provided by preserving individual existing or creating new NNBF projects within a community (de Vries et al, 2021). Thus, the value of this GIS approach affords analysis and understanding at parcel, local, and regional perspectives.
The benefits ranking of NNBF also enables informed decisions to protect and preserve these features via local and state decisionmaking processes including for tidal and nontidal wetlands permit reviews, RPA buffer impact applications, site plan reviews, land conservation planning and acquisition, and comprehensive planning. NNBF with the highest CBI ranking are those that have the most potential to provide the most benefits assessed in this analysis and are therefore of highest priority for those preservation actions. In this study, NNBF types with the highest average CBI ranking were tidal marshes and hybrid living shorelines with marsh (i.e., projects that plant vegetation along with an oyster sill or marsh sill). More than 90% of tidal marshes and 100% of living shorelines with marsh receive the highest ranking. In general, the highest CBI rankings were achieved by vegetated NNBF at low elevation, because these features have higher flooding mitigation potential scores and the greatest potential to qualify for programmatic incentive co-benefits. Other non-tidal wetland types (emergent, scrub-shrub, and forested) were also ranked highly, with slightly more than half of those NNBF (66%, 61%, and 67%, respectively) receiving the highest rank. Upland NNBF like wooded and scrub-shrub areas, and sandy features with little to no vegetation (e.g., beaches, dunes, and hybrid living shorelines with breakwaters), tended to receive a moderate ranking. Overall, more than 125,000, or 35%, of NNBF covering 1,523 km2 receive the highest ranking of 'most benefits provided'.
In addition, identifying gaps in NNBF benefits via the target areas can be used to direct limited resources towards new or restored NNBF in locations that increase the resilience of communities. Because of this study's focus on benefits for buildings, target areas are skewed towards relatively developed areas. This is reflected in the high prevalence of target areas (65%) where the SMM indicates that the shoreline is within a highly modified area and recommends seeking expert advice for shoreline management to address the complex landscape setting. Despite this, most target areas (71%) are adjacent to existing NNBF like tidal marshes and upland wooded areas, and 34% of target areas had SMM recommendations to create nonstructural or hybrid living shorelines as a preferred approach to addressing shoreline erosion. Creating same type NNBF within these areas would benefit habitat connectivity, and the presence of existing features could be an indicator of suitable conditions for an NNBF project.
An advantage of the indicator ranking framework and geospatial analysis is that it is flexible and can be adapted and customized to the community of focus by including additional ecosystem, economic, and/or cultural services. The approach can be modified to incorporate additional physical parameters of flooding mitigation action that NNBF can provide, additional vulnerable infrastructure such as roads, and an array of additional socio-economic criteria such as social vulnerability, agricultural commodities, consideration of the housing stock within the study area or economic valuation of coastal buildings, among others.
While this study was completed based on available data, future datasets and research could improve the methodology and geospatial modeling. Improvements include the addition of built flood water conveyance structures such as culverts and storm drains when delineating IPs, particularly in urban areas, as these types of structures can be conveyances for coastal flooding (Loftis et al., 2019;Habel et al., 2020). In coastal Virginia, adequately detailed datasets for these structures are incomplete and/or outdated. Another improvement to the overall methodology would be developing a model to directly link important offshore NNBF to communities they benefit. NNBF such as sand spits, submerged aquatic vegetation (SAV, e.g., seagrass), and barrier islands, are demonstrated to have flooding mitigation benefits for backshore communities (Pinsky et al., 2013;Hanley et al., 2014;Altman et al, 2021), and are therefore important as part of the overall picture of NNBF and community resilience.
To maximize the utility of outputs for communities, project outputs were developed in consultation with and for use primarily by local decision-makers. Locality representatives were engaged during data development and modeling, as well as post modeling, Buildings summarized by the number of NNBFs intersecting their IPs. Project outputs are served on the interactive mapping tool located at AdaptVA 1 . This screenshot shows the ranked NNBF layer, located under the Lands for Protection tab, with informational popup containing information useful to managers. Project outputs are served on the interactive mapping tool located at AdaptVA 1 . This screenshot shows target areas, the layer located under Restoration Opportunities, with its informational popup designed to guide decisions on potential NNBF suitable for that location.
to provide direction for data service and project documentation. Project information are served via an interactive map tool 1 housed on a collaborative website, AdaptVA. Figures 8, 9 show screenshots of the ranked NNBF output and informational popup located under the Lands for Protection section, and target areas and associated popup located under the Restoration Opportunities group. Project information and NNBF fact sheets are served on the CCRM website 2 and accompanied by outreach and training to potential data users from local and state agencies.

Conclusion
In summary, the inundation pathways presented here represent a novel application of a geospatial approach to assessing NNBF at a community level. IPs are used to inform an NNBF ranking framework, the Coastal Benefits Index, that considers multiple benefits of NNBF for resilience, including flooding mitigation services to flood-vulnerable and important critical community buildings, water quality improvement, and community financial benefits in terms of potential CRS credits. This approach is also used to identify priority target areas to restore or create new NNBF that would maximize benefits for communities.
While the approach detailed in this study has been tailored for coastal Virginia communities, it can be readily adapted to any community of focus through the inclusion of additional or alternative ecosystem, economic, social, and cultural services that NNBF can provide. In addition, using a relative assessment of these benefits rather than a quantitative measure of risk reduction or mitigation means that this approach can be adjusted based on newly available science and/or changes in the understanding of the risk mitigation that NNBF provide. Therefore, the best available science can be used to support local decision making to preserve, protect, restore, and create NNBF for community resilience.

Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: William & Mary ScholarWorks: https://doi.org/10.25773/d9pv-gb12.

Author contributions
CH and PM contributed to conception and design of the study. JHen and JHer performed the geospatial analyses. JHen wrote the first draft of the manuscript. JHen, PM, and JHer wrote sections of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.

Funding
This work was funded by the National Oceanic and Atmospheric Administration (NOAA) through grant # NA17NOS4730142. These data and related items have not been formally disseminated by NOAA, and do not represent any agency determination, view, or policy.