output: word_document always_allow_html: yes

Goals

Pressures

"The pressure score, p, describes the cumulative impact of ecological and social stressors in a given year and region which tend to depress the goal score in future years. Pressure scores range from 0 to 1, and include both ecological (p_E) and social pressures (p_S), such that:

\[p=γp_E+(1-γ)p_S\]

where γ=0.5 is the relative weight for ecological vs. social pressures categories. We default to equal weighting as little evidence was available to justify or quantify unequal weights between ecological and social pressures categories. It may be that future work can inform unequal weighting terms γ for individual goals.

For each goal, subgoal, or goal element (e.g. specific habitat), we calculated pressures as an impact-weighted cumulative impact for each pressure category \(p_{ecol}^{goal}\) and \(p_{soc}^{goal}\). Impact weights are based on a goal’s sensitivity \(w_(s_i)^{goal}\) to specific stressors s_i ranked as low \((w_(s_i)^{goal}=1)\), medium \((w_(s_i)^{goal}=2)\), high \((w_(s_i)^{goal}=3)\), or no impact \((w_(s_i)^{goal}=NA)\), as determined by peer-reviewed literature and expert judgment (Fig. SM2 shows the matrix of stressors, goals, and weights). The denominator represents the maximum stressor impact weight for that category and goal. If cumulative pressure load for a goal/component combination exceeds the maximum possible stressor intensity, we cap it to 1.0, i.e. the equivalent to an individual stressor at maximum stress and intensity. \[p_{cat}^{goal}=min(\frac{1}{(w_{cat,max}^{goal}}\sum_{i=1}^N w_{s_i}^{goal} ×s_i, 1)\]" OHIBC

Ecological Pressure

We included five subcategories of ecological stressors/pressures: water pollution, climate change, habitat destruction, fishing pressure, and species pollution

– more about this here and equations; unsure if the way BC calculated it is the same –

Climate Change

cc_acid

Ocean Acidification

Ocean acidification data were obtained from the East Coast Ocean Acidification Product Suite (ECOAPS)1, which is produced by the NOAA Coral Health and Monitoring Program. Changes in the aragonite saturation state (\(\Omega_{arag}\)) can be attributed to changes in the concentration of CO2 and thus we use this measurement as a proxy for increasing CO2. This layer is rescaled between 0-1 using two accepted thresholds; 1.5\(\Omega_{arag}\) the level when organism show mild dissolution2 and 1\(\Omega_{arag}\) the lethal saturation level3. References

cc_sst

Sea Surface Temperature

Sea surface temperature (SST) data were obtained from the Coral Reef Temperature Anomaly Database (CoRTAD) (Casey et al. 2015), which is produced by the NOAA National Center for Environmental Information (NCEI) using 4.6 km (nominally 21 km2 at the equator) Advanced Very High Resolution Radiometer (AVHRR) Pathfinder Version 6 (v6 data) SST data (http://www.nodc.noaa.gov/sog/cortad/). Weekly SST data are used to compute the standard deviation (SD) of SST’s per pixel across all years. We define an anomaly as exceeding the standard deviation of SSTs from the climatology for that location (i.e., grid cell) and week of the year. The frequency of weekly anomalies was calculated for each year in the dataset. We then quantified the difference between the number of anomalies in the 5 most recent years and the 5 oldest years in the dataset. The 99.99th quantile of raster values from all years was used as the reference point to rescale the layer from 0 to 1, and the mean value of the raster cells within each OHI region was calculated.

cc_slr

Sea Level Rise

The sea level rise pressure layer is derived from satellite altimetry data4. Monthly mean sea level anomalies since 1993 track changes in sea level (mm) compared to a reference period from 1993-2012. Raw monthly data are provided on a 0.25x0.25 degree grid. These data were clipped to within 3 nautical miles of the coast, and monthly data layers were aggregated and averaged across pixels to compute mean sea level anomalies. We are using a reference point 110% larger than the maximum value in the entire Northeast region from 1993-2017 to scale the layer from 0-1. The current maximum value is unlikely to be the future max, thus setting this reference point allows for an increase in sea level rise pressure into the future. All negative values were set to zero (i.e., no negative pressure), such that only positive sea level rise values effected the pressure. The mean value of the raster cells within each OHI region was calculated for regional pressure score.

References

Fishing

prs_fishing

Habitat Destruction

hd_intertidal

Trampling

Beaches and coastlines are very sensitive to disturbances, even a few individuals can cause damage. As population densities along the coastlines increases, so does the mounting pressure on intertidal habitats from the presence of humans utilizing them.

We used 2010 census data downloaded from NASA’s Socioeconomic Data and Applications Center5. The GeoTiff files represent population densities in 1 km2 grids. To focus on populations that would contribute coastal trampling pressure, we selected areas within 25 miles of the coastline. The distribution of the population densities were very skewed, due to large cities like New York and Boston, so we log transformed the population densities, since a small number of people can have a very large effect on habitats. The 99.99th quantile of raster values was used as the reference point to rescale the layer from 0 to 1, and the mean value of the raster cells within each OHI region was calculated.

References

hd_subtidal

prs_hardening

Coastal Hardening

Artificially hardening coastlines is a reinforcement method to protect developments and land areas along the coast from erosion6. However this is often a short sighted solution. Coastal hardening costs money to maintain and is detrimental to the habitat, such as narrowing of the intertidal zone7 or loss of salt marsh habitats.

Data from the NOAA Environmental Response Management Application8 portal was used to calculate the percentage of armored coastline in each of the regions. We included all man made structures such as seawalls, riprap, jetties, bulkheads, etc within a 1 km distance of the coastline. The percentage hardened in each region was used directly as the pressure score for this layer.

References

prs_aqua_hab

Aquaculture Development

The effects of aquaculture facility location and production on habitats differ greatly across regions, methods, and species produced. Generally, shellfish aquaculture is the least harmful, because there is less change to the habitat, not much extra structure needs to build, and they actual clean the water through filter feeding9. On the other hand, shrimp production in some areas require the removal of mangroves and the growing process creates poor water quality conditions10.

Sustainability scores for farmed species were obtained from the Monterey Bay Aquarium Seafood Watch Seafood Recommendations11. This assessment incorporates 10 criteria: availability and robustness of data, effluent discharge, effects on habitats, use of chemicals, escapes during production, disease (pathogen, and parasite interaction), stock source, predator and wildlife moralities, and escapes during transit when calculating the overall sustainability score. We used the effects on habitats criteria to inform aquaculture production’s contribution to the habitat destruction layer. The criteria assesses the impact or risk of impact the aquaculture farm has on the “ecosystem services” provided by habitat within the farm the immediate/neighboring areas.

To calculate the score for this layer we rescaled the criteria score from 0-10 to 0-1, scaled the production by multiplying it with the new criteria score, summed up all the scaled production in each region, and then divided it by the total area of the region to distribute the pressures. Since we know that aquaculture production is in the early stages and will likely grow much larger than the current production levels, we used a reference point of twice the current maxium of the scaled distributed production in each region, and then rescaled the pressure score 0-1.

References

prs_sasi

Species Pollution

prs_disease

Disease Amplification

One of the observed side effects during aquaculture production is the increased prevalence of disease, pathogens, and parasite spreading in the farms due to high densities of individuals12. These can spread out through neighboring water bodies and infect wild stocks.

Sustainability scores for farmed species were obtained from Monterey Bay Aquarium Seafood Watch Seafood Recommendations13. This assessment incorporates 10 criteria: availability and robustness of data, effluent discharge, effects on habitats, use of chemicals, escapes during production, diseases (pathogen, and parasite interaction), stock source, predator and wildlife mortalities, and escapes during transit when calculating the overall sustainability score. We used the disease, pathogen, and parasite interaction criteria to inform the contribution of aquaculture production to the species pollution layer. The criteria assesses the increases in transmission or retransmission of diseases or pathogens to wild fish stocks due to the amplification caused by aquaculture production.

To calculate the score for this layer we rescaled the criteria score from 0-10 to 0-1, scaled the production by multiplying it with the new criteria score, summed up all the scaled production in each region, and then divided it by the total area of the region to distribute the pressures. Since we know that aquaculture production is in the early stages and will likely grow much larger than the current production levels, we used a reference point of twice the current maxium of the scaled distributed production in each region, and then rescaled the pressure score 0-1.

References

sp_alien

Nonindigenous Species

The USGS Nonindigenous Aquatic Species database was used to inform this layer. Locally specific occurrence data was used to calculate the total number of occurrences in each region that have become established and eradicated14. This layer was calculated by dividing the number of species eradicated by the maximum number of occurances from 1970-2017.

References

sp_genetic

Aquaculture Escapes

There is the risk of some farmed individuals escaping into the wild stocks or spawning events introducing farmed genetic recruits into the wild gene pool from aquaculture facilities. This is generally considered detrimental to wild stocks, and that diluting the native gene pool results in lower fitness stocks15.

Sustianability scores for farmed species were obtained from the Monterey Bay Aquarium Seafood Watch Seafood Recommendations16. This assessment incorportates 10 criteria: avaliablity and robustness of data, effluent discharge, effects on habitats, use of chemicals, escapes during production, diseases (pathogen, and parasite interaction), stock source, predator and wildlife mortalities, and escapes during transit when calculating the overall sustainability score. We used the escapes during production criteria to inform contribution of aquaculture to the species pollution layer. The criteria assesses the risk level of impacts on wild stock populations’s fitness and habitat caused by the aquaculture escapes of non-native or genetically distinct fish.

To calculate the score for this layer we rescaled the criteria score from 0-10 to 0-1, scaled the production by multiplying it with the new criteria score, summed up all the scaled production in each region, and then divided it by the total area of the region to distribute the pressures. Since we know that aquculture production is in the early stages and will likely grow much larger than the current prodcution levels, we used a reference point of twice the current maxium of the scaled distributed production in each region, and then rescaled the pressure score 0-1.

References

Social Pressure

prs_social

Resilience

“Resilience for each goal and region, r, is based on three components: ecological integrity, r_ecol, and regulatory efforts that target specific ecological pressures, r_reg, which combine to address resilience to ecological pressures; and social integrity, r_soc, which addresses social pressures. Each resilience category contains one or more layers reflecting the magnitude of resilience within each region for each year; layers are “activated” to address specific pressures acting on specific goals based on a resilience matrix (Fig. matrix), and active layers are summed to determine a score for each resilience category. Each layer is constrained from 0 to 1.

\[r=γ((\frac{r_{ecol}+r_{reg})}{2})+(1-γ)r_{soc}\]

These components are weighted such that resilience to ecological pressures (i.e., ) and resilience to social pressures (i.e., ) reflect the proportional contribution of ecological and social pressures in the pressures model, i.e. \(γ=0.5\) (BC OHI)”

I. Ecological Integrity

“An intact biodiverse ecosystem provides general resilience to ecological pressures by ensuring the system’s ability to maintain functionality in the face of stressors imposed by human activity and climate change. For OHIBC, we consider the area-weighted average conservation status of all species found in the coastal zone (3 nmi offshore) (as resilience to coastal pressures) and found within the entire EEZ (as resilience to pressures not limited to the coast). (BC OHI)”

II. Regulatory Efforts/Resilience

“Regulatory resilience describes the institutional measures (e.g., rules, regulations, and laws) designed to address ecological pressures. The regulatory resilience datasets are grouped into five categories that address the 5 pressure categories: water pollution, climate change, fishing pressure, habitat destruction and species introductions (invasive species and genetic escapes)(OHI Goals website).” The resilience layers for each category are composed of 3 components: presence of regulations, implementation/enforcement of regulations, and effectiveness/compliance of regulations. We searched for data that would serve as direct measurements or proxy indicators of the extent and successfulness of these components. Federal or state set targets were used as references for calculating resilience scores, and when this information was not available we used choose references we felt reasonable. In cases where multiple data layers are used to calculate one component, they were averaged, and then the three regulatory resilience categories are equally weighted in the overall resilience layer score. (reg metrics table)

\[regulatory resilience_{layer x} = \frac{presence + implemenation/enforcemnet + effectiveness/compliance}{3}\]

1. Water Pollution

res_water_pollution

The data used to measure water pollution resilience are from the Environmental Protection Agency (EPA) Enforcement and Compliance History Online (ECHO) database. The EPA’s management of water discharge facilities is used as an indicator for other water pollution prevention regulations, implementation and effectiveness in each region. This database includes historic and current information on facilities that discharge waste and pollution. There are 4 data layers used in the water pollution resilience layer.

Presence

  1. The existence of the EPA and the ECHO database imply the presence of regulations, so each region received a score of 1.

Implementation/Enforcement

  1. We calculated the total number of discharge facilities that were inspected in each region and the total number of discharge facilities. We used 50% of facilities in each region inspected as our target, since that is the EPA established “National Goals”17.

Effectiveness/Compliance

  1. We calculated the total number of facility violations in each region and assigned a binary value to each facility, at least one violation that year or no violation so as to not over penalize regions were a single facility had multiple violations. There are 4 different types of violation, but we did not differentiate between them. We used a reference point of 0%, no violations1.

  2. Each facility has a permit which requires discharge monitoring reports at various intervals depending on the permit. We calculated the percent of expected discharge monitoring reports submitted and rescaled it to the target of 95% submitted, which is based on the EPA established “National Goals”1.

References

2. Climate Change

res_climate_change

Multiple data sources were used to calculate the climate change resilience layer including American Council for Energy Efficiency (ACEE), America’s Pledge Report, US Energy Information Administration (EIA), and the Environmental Protection Agency (EPA). There are 4 data layers used in the climate change resilience layer.

Presence

  1. The state’s scores from the annual ACEE State Energy Efficiency Scorecard, which scores each state on their energy efficient policies, were compiled from 2005-201718. We used the maximum score possible of 50 as the reference point, and then rescaled 0-1.

  2. We incorporated state level climate friendly actions since the United States exit from the Paris Climate Agreement in this layer. This information was taken from America’s Pledge Report which lists the number of climate friendly actions in 33 different categories. For this layer, we removed the energy efficiency categories as to not double count policies in this area, since we used the ACEE scorecard are a more comprehensive assessment of energy polices. The states are scored based on actions taken in 30 different categories, and the target is at least one action in each of the category19.

Implementation/Enforcement

  1. The EIA energy consumption data is used to resume how well each state is doing with regards to implementation of their clean energy consumption targets. Scores are calculated by taking the EIA energy consumption data, calculating the percent of energy consumed that comes from clean sources, and dividing it by the target emissions. Clean energy consumption targets are calculated by performing a linear regression from 2004 to the target year of their consumption goal. While there are a variety of energy plans and portfolios in the US that differ from state to state, every state in our region of interest has a “renewable energy portfolio standards”, so we do not include nuclear energy production as part of the total renewable energy consumption target20 (Table X21).
Table X. State Set Renewable Energy Consumption Targets
State Goal
Maine 40% renewable by 2017
New Hampshire 25.2% renewable by 2025
Connecticut 44% renewable by 2030
Massachusettes 35% renewable 2030
New York 50% renewable by 2030
Rhode Island 38.5% renewable by 2035
  1. The EPA carbon emissions data is used to resume how well each state is doing with regards to implementation of their carbon emissions targets. Scores are calculated by taking the EPA emissions data and dividing it by the target emissions. Total annual carbon emissions were pulled from the EPA website22. Carbon emission targets are calculated by performing a linear regression from 2004 to the target year of state set emissions goal (Table X23).
## Warning: Missing column names filled in: 'X3' [3]
Table X. State Set CO2 Emissions Targets
State Goal
Maine 80% reduction below 1990 levels by 2050
New Hampshire 80% reduction below 1990 levels by 2050
Connecticut 80% reduction below 2001 levels by 2050
Massachusetts 80% reduction below 1990 levels by 2050
New York 100% reduction below 1990 levels by 2050
Rhode Island 80% reduction below 1990 levels by 2050

Effectiveness/Compliance

The same data as Implementation/Enforcement are used as a proxy indicator for this layer.

References

3. Fishing Pressure

res_fishing

Multiple data sources were used to calculate the fishing pressure resilience layer including National Oceanographic and Atmospheric Administration (NOAA) Fisheries, NOAA Office of Law Enforcement, NOAA Fisheries Observer Program, and the United States Geographic Survey (USGS).There are 7 data layers used in the water pollution resilience layer.

Presence

  1. Of the species reported in the NOAA landings data, we calculated how many have had a formal stock assessment24. The target is that every species caught has been assessed at least once. However, since the more heavily fished species are more likely to have been assessed, we weighted these scores by the amount of catch. By doing so we still count species that are rarely caught and have not yet been assessed, but they don’t “overrule” major fishery species that are assessed.

  2. Stock assessments are very valuable, but need to be regularly reassessed to monitor any changes in the stocks1. We used the NOAA stock assessment to calculated the percentage of stocks that are adequately assessed. Using expert knowledge, we set the target of an assessment at least once every 5 years.

  3. Using the USGS Protected Areas Database we identified areas where all fishing is prohibited year-round; areas where fish would have refuge from fishing25. We calculated percentage of area in each region that fit these criteria and then, using a target of 30%, rescaled the pressure 0-1. We used a target of 30%, because this is the observed level of ocean closure needed provide fisheries resilience26.

  4. Most species that are fished do not have a formal stock assessment, but they are managed or regulated in some way, either at the federal, state, or local level. We compiled a list of all species caught and whether or not they are included in management plans or regulatory policies. From this list we calculated what percentage of species landed have either a management plan or some type of regulation. We rescaled from 0-1 using the reference of 100%.

Implementation/Enforcement

  1. Using data obtained from the NOAA Fisheries Observer Program, we calculated the number of sea days a fisheries observer spent on board a fishing vessel27. Then, using past SBRM Annual Discharge Reports with Sea Day Allocation, we pulled the allotted number of sea days for each year and used those as the reference. The allotted sea days was recalculated annually based on the previous years number of sea days and catch the number of allotted sea days. We then rescaled from 0-1, with the goal to reach the allotted number of sea days.

  2. We extracted the total number of staff, patrols, and outreach events OLE partakes from the NOAA OLE Northeast Region Annual reports. These are summed up and rescaled 0-1 using the highest historic total as the reference point28.

Effectiveness/Compliance

  1. Using data from the NOAA OLE Northeast Region Annual report, this layer was calculated by taking the number of investigations and dividing it by the number of enforcement actions29.

References

4. Habitat Destruction

res_hab_destruction

The data sources used to calculate the habitat destruction resilience layer includes Surfrider Foundation, NOAA Greater Atlantic Regional Fishing Office (GARFO), and NOAA Office of Law Enforcement.There are 4 data layers used in the water pollution resilience layer.

Presence

  1. The Surfrider Foundation releases an annual State of the Beach Reportcard (since 2017) which scores each coastal state on their policies to protect their beaches from coastal erosion, haphazard development, and sea level for a total potential score of 10. We rescaled these scores 0-130.

  2. NOAA GARFO shapefiles were used to identify areas where fisheries regulations are in place that prohibit trawling and dredging year-round31. We used a target of 30% protected total area in each region.

Implementation/Enforcement

  1. The total number of staff, patrols, and outreach events OLE partakes in is pulled from the NOAA OLE Northeast Region Annual report. These are summed by year and divided by the highest historic total32.

Effectiveness/Compliance

  1. Using data from the NOAA OLE Northeast Region Annual report, this layer was calculated by taking the number of investigations and dividing it by the number of enforcement actions3.

References

5. Species Pollution

res_alien

We used data/information from Coast Guard, NOAA Office of Law Enforcement (OLE), USGS Nonindigenous Aquatic Species database to calculate the species pollution resilience layer. There are a total of 4 data layers used in the species pollution resilience layer.

Presence

  1. The existence of the Ballast Water Management Plan (BWMP) Requirements implemented and monitored by the US Coast Guard imply the presence of regulations, so each region received a score of 1. It appears that information exists on the number of vessels inspected and how many are in violation of the BWMP, however this is not available to the public. A declassified version or a similar data set would be optimal for quantifying each component of this resilience layer, if it were to become available.

  2. The existence of aquaculture permitting system imply the presence of regulations, so each region received score of 1. For calculation purposes this layer does not change the overall score in this layer, but is used more as a placeholder for future iterations of the OHI as aquaculture becomes more streamlined, regulated, and large-scale in the Northeast and better data becomes available.

Implementation/Enforcement

  1. We extracted the total number of staff, patrols, and outreach events OLE partakes in from the NOAA OLE Northeast Region Annual reports. These are summed up and rescaled 0-1 using the highest historic total as the reference point. This information was given to us in PDF report form for 4 years33. This data in raw form, over the entire period of time would be optimal, if this data becomes available.

Effectiveness/Compliance

  1. The USGS Nonindigenous Aquatic Species database was used to inform this layer34. We used locally specific occurrence data to calculate the total number of occurrences in each region that have become established and eradicated. This layer was calculated by dividing the number of species eradicated by the number of new occurrences in each region each year.

References

III. Social Integrity

res_social

“Social resilience describes the social integrity of coastal communities that allow for adaptive responses to social and ecological pressures (BC)” This layer was calculated using data from the US Opportunity Index, Beacon Hill Competitiveness Reports, League of Conservation Voters Scorecards, and the NOAA Office of Science and Technology Social Indicators scores.

\[socialintegrity = \frac{OpportunityIndex + CompetitivenessReport + LCV + SocialIndicators}{4}\]

US Opportunity Index:

“The Opportunity Index was created by Opportunity Nation35. The Opportunity Index aims to use indicators beyond economics and to measure the access of opportunities in communities across the country. The index measures four dimensions of community well-being: economy, education, health, and community. Index scores from 2012, and 2014-2017 were used. The Opportunity Index scores each state on a scale of 0-100, which we then rescaled 0-1 and applied directly into each state resilience measure.

Beacon Hill Competitiveness Report:

The State Competitiveness Report is produced by Beacon Hill Institute on an annual basis and aggregates key microeconomic variables into a single index to measure the economic competitiveness of each state. A state is considered by Beacon Hill Institute to be competitive “if it has in place the policies and conditions that ensure and sustain a high level of per capita income and its continued growth. Sub-indicators used to assess each state include: Government and fiscal policies, Security, Infrastructure, Human resources, Technology, Business incubation, Openness, and Environmental Policy36.” Index scores from 2005-2017 were used. The State Competitiveness Report scores each state on a scale of 0-10, which we rescaled 0-1 and applied directly into each state resilience measure.

League of Conservation Voters Scorecards:

The League of Conservation Voters is an environmental group that “advocates for sound environmental laws and policies, holds elected officials accountable for their votes and actions, and elects pro-environment candidates.”37. One of their methods for doing so is through and annual scorecard grading Congress members on how they vote on environmentally related issues (for or against). We use these scores to determine how the people in each state value their environments health and protection. How state elected officials vote on these issues is a proxy for the priorities of their constituents who voted them into office.

NOAA Office of Science and Technology Social Indicators:

NOAA aims to research and better understand the impacts fisheries management decisions have on fishing communities, so they have developed a new method measuring social indicators of the fishing communities to quantify vulnerability and resilience in US coastal regions. Each community is scored in 4 indicators: Social Vulnerability, Gentrification Pressure, Sea Level Rise Vulnerability, Fishing Engagement and Reliance38. The factor analysis technique was used to construct the indicators which “consisted of a principal component analysis with a varimax rotation (methods doc).” Since there was no cap on the maximum value for vulnerability and resilience, we averaged each indicator by region and then rescaled 0-1 using the maximum value as a 1. We then averaged all four rescaled indicators to calculate each regions resilience measure.

References

Resilience Data Assessment

insert the matrix for quality of data here. use the highlighting functions in kable for the green, yellow, red


  1. NOAA Coral Health and Monitoring Program, Ocean Acidification Product Suite, https://www.coral.noaa.gov/accrete/oaps.html

  2. Bednaršek, Nina, et al. “Systematic Review and Meta-Analysis Toward Synthesis of Thresholds of Ocean Acidification Impacts on Calcifying Pteropods and Interactions With Warming.” Frontiers in Marine Science, vol. 6, 2019, doi:10.3389/fmars.2019.00227.

  3. Ries, J. B., et al. “Marine Calcifiers Exhibit Mixed Responses to CO2-Induced Ocean Acidification.” Geology, vol. 37, no. 12, 2009, pp. 1131–1134., doi:10.1130/g30210a.1.

  4. AVISO+ Satelite Altimetry Data, MSLA - Monthly Mean and Climatology Maps of Sea Level Anomalies, http://www.aviso.altimetry.fr/en/data/products/sea-surface-height-products/global/msla-mean-climatology.html

  5. “U.S. Census Grids.” Summary File 1, v1: U.S. Census Grids | SEDAC, NASA Socioeconomic Data and Applications Center (SEDAC), 2010, doi.org/10.7927/H4TB14TN.

  6. Gittman, Rachel K, et al. “Engineering Away Our Natural Defenses: an Analysis of Shoreline Hardening in the US.” The Ecological Society of America, John Wiley & Sons, Ltd, 1 Aug. 2015, esajournals.onlinelibrary.wiley.com/doi/full/10.1890/150065.

  7. Dugan, Jenifer E., et al. “Ecological Effects of Coastal Armoring on Sandy Beaches.” Wiley Online Library, John Wiley & Sons, Ltd (10.1111), 3 June 2008, onlinelibrary.wiley.com/doi/full/10.1111/j.1439-0485.2008.00231.x.

  8. NOAA Office of Response and Restoration. Environmental Response Management Application (ERMA) available at https://erma.noaa.gov/atlantic/erma.html#/layers=1+13763+16973+35410+491&x=-78.22196&y=38.23326&z=6&panel=layer

  9. Rose, Julie M., et al. “Comparative Analysis of Modeled Nitrogen Removal by Shellfish Farms.” Marine Pollution Bulletin, vol. 91, no. 1, 2015, pp. 185–190., doi:10.1016/j.marpolbul.2014.12.006.

  10. PaEz-Osuna, Federico. “The Environmental Impact of Shrimp Aquaculture: Causes, Effects, and Mitigating Alternatives.” Environmental Management, vol. 28, no. 1, 2001, pp. 131–140., doi:10.1007/s002670010212.

  11. Monterrey Bay Aquarium Seafood Watch, available at: https://www.seafoodwatch.org/seafood-recommendations

  12. Bondad-Reantaso, Melba G., et al. “Disease and Health Management in Asian Aquaculture.” Veterinary Parasitology, vol. 132, no. 3-4, 2005, pp. 249–272., doi:10.1016/j.vetpar.2005.07.005.

  13. Monterrey Bay Aquarium Seafood Watch, available at: https://www.seafoodwatch.org/seafood-recommendations

  14. Fuller, P., and Benson, A., 2016, Nonindigenous Aquatic Species Database: U.S. Geological Survey, http://nas.er.usgs.gov.

  15. Mcginnity, Philip, et al. “Fitness Reduction and Potential Extinction of Wild Populations of Atlantic Salmon, Salmo Salar, as a Result of Interactions with Escaped Farm Salmon.” Proceedings of the Royal Society of London. Series B: Biological Sciences, vol. 270, no. 1532, 2003, pp. 2443–2450., doi:10.1098/rspb.2003.2520.

  16. Monterey Bay Aquarium Seafood Watch, avaliable at: https://www.seafoodwatch.org/seafood-recommendations

  17. U.S. Environmental Protection Agency 2019, The Integrated Compliance Information System (ICIS) Federal Enforcement and Compliance (FE&C), viewed June 2019, https://echo.epa.gov/tools/data-downloads

  18. American Council for an Energy-Efficient Economy 2019, The State Energy Efficiency Scorecard. Washington D.C.: American Council for an Energy-Efficient Economy.

  19. America’s Pledge, “America’s Pledge Phase 1 Report: States, Cities, and Businesses in the United States Are Stepping Up on Climate Action,” November 2017.

  20. U.S. Energy Information Administration 2019, Total Energy Data, Viewed July 2019, https://www.eia.gov/

  21. Database of State Incentives for Renewables and Efficiency, Renewable Portfolio Standards. (Raleigh, N.C. State University, 2017.

  22. U.S. Environmental Protection Agency 2019, CO2 Emissions from Fossil Fuel Combustion - Million Metric Tons CO2 (MMTCO2), Viewed August 2019, https://www.epa.gov/sites/production/files/2017-09/documents/co2ffc_2015.pdf

  23. Center for Climate and Energy Solutions, “U.S. State Greenhouse Gas Emissions Targets”, 8 July 2019, www.c2es.org/document/greenhouse-gas-emissions-targets/.

  24. Fisheries landings and stock data, how to cite?

  25. U.S. Geological Survey (USGS) Gap Analysis Project (GAP), 2018, Protected Areas Database of the United States (PAD-US): U.S. Geological Survey data release, https://doi.org/10.5066/P955KPLE.

  26. Gaines, S. D., et al. “Designing Marine Reserve Networks for Both Conservation and Fisheries Management.” Proceedings of the National Academy of Sciences, vol. 107, no. 43, 2010, pp. 18286–18293., doi:10.1073/pnas.0906473107.

  27. Fisheries Observer data, how to cite? Debrah says there isn’t a specific way for us to cite it

  28. (National Oceanographic Atmospheric Association Office of Law Enforcement, Northeast Enforcement Division: Annual Report Fiscal Year 2018-2015, 2018-2015)

  29. (National Oceanographic Atmospheric Association Office of Law Enforcement, Northeast Enforcement Division: Annual Report Fiscal Year 2018-2015, 2018-2015)

  30. Surfrider Foundation (2017), State of the Beach Report. San Clemente, CA

  31. NOAA Greater Atlantic Region Fisheries Office, GIS Datasets. https://archive.fisheries.noaa.gov/garfo/educational_resources/gis/data/indextab.html

  32. (National Oceanographic Atmospheric Association Office of Law Enforcement, Northeast Enforcement Divison: Annual Report Fiscal Year 2018-2015, 2018-2015)

  33. (National Oceanographic Atmospheric Association Office of Law Enforcement, Northeast Enforcement Division: Annual Report Fiscal Year 2018-2015, 2018-2015)

  34. Fuller, P., and Benson, A., 2016, Nonindigenous Aquatic Species Database: U.S. Geological Survey, http://nas.er.usgs.gov.

  35. Opportunity Nation. 2012, 2014-2017. Opportunity Index. Available at http://opportunityindex.org/

  36. Beacon Hill Institute. 2005-2017. State Competitiveness Report. Available at http://beaconhill.org/economic-competitiveness/

  37. League of Conservation Voters, 2017 National Environmental LCV Scorecard. Washington D.C.: League of Conservation Voters. (Same citation for 2005-2017, different reports)

  38. NOAA Office of Science and Technology. Social Indicators of Fishing Community Vulnerability and Resilience in the U.S. Southeast and Northeast Regions. Ava liable at https://www.st.nmfs.noaa.gov/humandimensions/social-indicators/map