Co-occurring anthropogenic activities influence coastal ecosystems around the world. Notions of ecological exposure are promising indicators to better understand environmental status and enhance ecosystem protection. This study characterized anthropogenic exposure in the context of multiple human activities on coastal benthic ecosystems at a scale of <100 km. Using a particle diffusion model and fishing event data, we developed an exposure index for seven human activities (aquaculture, artificial structures, dredging, fisheries, runoff, sewers and shipping) in a Canadian industrial harbour area. A generally low cumulative exposure was obtained, with the highest values observed directly in front of the city and industrial areas. Derived exposure indices explained a portion of the benthic community structure (R2 = 0.22), suggesting an ecological link between the exposure of species and their vulnerability to human activities. Such tools are relevant in data-poor environments where proxies are required to assess the state of an ecosystem, facilitating the application of ecosystem-based management.

Management of coastal marine ecosystems requires efficient monitoring of ecosystem components, including human activities, in order to accurately guide environmental conservation initiatives. This need is especially true in the face of intensifying and ever diversifying human activities in marine ecosystems, with the omnipresence of their impacts (Halpern et al., 2019). Environmental assessments should thus consider the cumulative effects of multiple co-occurring human activities to best describe the current and anticipated states of marine ecosystems (Crain et al., 2008; Brown et al., 2014; Côté et al., 2016). With 40% of humanity living less than 100 km from coasts (Socioeconomic Data and Applications Center, 2020), coastal habitats and communities are influenced by a wide variety of human activities from terrestrial, freshwater and marine realms (Feist and Levin, 2016; Micheli et al., 2016). There is thus an urgent need to better understand how these ecosystems may by impacted by anthropogenic influences to accurately support their protection.

Integrative approaches, such as ecosystem-based management and marine spatial planning, are important tools for assessing, monitoring and managing human activities in coastal ecosystems (Margules and Pressey, 2000; Link, 2002; Pikitch et al., 2004; Levin et al., 2009; Santos et al., 2019). Ecosystems are complex entities composed of interconnected components, including biological communities, habitats and human activities, each governed by their own dynamics. An attempt to capture the complexity of this network of interactions has been made with the description of socio-ecological systems, often used to address community resilience issues (Berkes et al., 2000; Redman et al., 2004; Young et al., 2006; Díaz et al., 2011; Glaser et al., 2012). However, efficient environmental monitoring, particularly in the context of multiple human activities, is still a challenge and requires new tools to improve management capacity (Crain et al., 2008; Darling and Côté, 2008; Séguin et al., 2014; Piggott et al., 2015; Galic et al., 2018; Hodgson et al., 2019; Carrier-Belleau et al., 2021).

Many studies have focused on the influence of cumulative impacts on ecosystems around the globe (Korpinen and Andersen, 2016), including in-situ field studies (e.g., Ocaña et al., 2019; D’Alessandro et al., 2020), geospatial modelling (e.g., Ban et al., 2010; Stelzenmüller et al., 2010; Parravicini et al., 2012; Okey et al., 2015; Beauchesne et al., 2020) and experimental manipulations (e.g., Beermann et al., 2018; Carrier-Belleau et al., 2022). In particular, Halpern et al. (2008), later updated by Halpern et al. (2019), proposed a comprehensive cumulative impact score for marine ecosystems. This score was calculated by combining spatial data on the exposure (co-occurrence between ecosystems and human pressures) and vulnerability (how ecosystem components react to this pressure) of ecosystems to 17 human activities (Wilson et al., 2005; Halpern et al., 2007, 2008). The various facets of this score represent important addition to the cumulative impacts literature, highlighting a ubiquitous anthropogenic footprint on marine ecosystems (Halpern et al., 2019). However, several limitations have been identified in these studies, such as the inclusion of very diverse pressures in a common metric, the proper assessment of spatial and temporal variability, the description of how ecosystem components respond to impacts, the inclusion of non-additive effects and the establishment of non-impacted reference conditions (Halpern and Fujita, 2013; Korpinen and Andersen, 2016; Hodgson et al., 2019).

Relating exposure with ecological indicators at a fine spatial resolution is a way to overcome some of these challenges and to address some shortcomings of cumulative impacts studies, especially ecosystem responses to pressure and possible emergent effects. Many relationships between biodiversity and local anthropogenic influence have been observed worldwide (Millenniun Environmental Assessment, 2005; Andersen et al., 2015; Solan and Whiteley, 2016; Ellis et al., 2017), such that ecological indicators provide relevant insights to understand the effects of human activities on the environment. As defined by Pinto et al. (2009), indicators describe ecosystems to determine a status with quantitative data, for example by using key habitat variables or the abundance of characteristic species (Borja et al., 2012; Teixeira et al., 2016). Macrobenthic invertebrates, a highly diverse biological component whose links with human activities have been described in a variety of ecosystems, may serve as ecological indicators (Pearson and Rosenberg, 1978; Grall and Glémarec, 1997; Teixeira et al., 2016). Various benthic species are characterized by a sedentary lifestyle and a relatively long lifespan which tends to reflect medium-term environmental conditions, resulting in adaptation or local extinction when disturbed (e.g., Dauer, 1993; Borja et al., 2000; Wei et al., 2020).

Linking cumulative pressure and biodiversity assessments at a spatial extent below 100 km, this study evaluates the influence of anthropogenic activities on local coastal ecosystems. In this context, the specific objectives of this study are to (i) model the exposure of benthic ecosystems to multiple anthropogenic activities at a local scale and (ii) evaluate how well this index reflects benthic conditions. We expect that the structure of biological communities within high exposure areas (“anthropogenic hotspots”) will present evidence of disturbance, such as lower diversity and the presence of opportunistic species, compared to the rest of the study area.

2.1. Study area

As a case study, we focused on the industrial harbour area of Sept-Îles (Québec, Canada; Figure 1). Located in the Gulf of St. Lawrence, a priority management area identified by Fisheries and Oceans Canada and a major economic region for Québec (Department of Fisheries and Oceans, 2009; Beauchesne et al., 2016; Daigle et al., 2017; Schloss et al., 2017), Sept-Îles is the fourth largest Canadian port in terms of total exchanged goods and the second largest in Québec in 2020 (Statistics Canada, 2011; Binkley, 2020; Ferrario et al., 2021; Port de Sept-Îles, 2021). Available ecological data on coastal ecosystems in this region were limited, which supported the need to characterize benthic ecosystems and their relation to coastal human activities (Snelgrove et al., 2012; Carrière, 2018; Dreujou et al., 2020b, 2021).

Figure 1.

Map of stations sampled in Baie des Sept Îles (Canada). Isolines correspond to the regional bathymetry (depth in meters).

Figure 1.

Map of stations sampled in Baie des Sept Îles (Canada). Isolines correspond to the regional bathymetry (depth in meters).

Close modal

The area includes Baie des Sept Îles and the archipelago at its entrance, covering approximately 200 km2 (Figure 1). Bathymetry is shallow within the bay, with a maximum depth of 50 m at its entrance, then becoming deeper (up to 200 m) in the archipelago (Dutil et al., 2012). The general sediment profile is sandy-silty, with a small fraction of gravel. Benthic communities are diverse with a high density of annelids, arthropods and mollusks (Dreujou et al., 2020b). This region has sub-Arctic environmental conditions, with sea ice formation in November–December and substantial freshwater runoff due to snowmelt in April (Demers et al., 2018). The area is characterized by strong tidal currents, resulting in an estuarine circulation within the bay, along with freshwater inputs from multiple streams (Shaw, 2019).

Local industrial operations include aluminium production in plants at the Pointe-Noire sector and the southeastern part of the city of Sept-Îles, international shipping of iron ore through bulk carriers (reaching 33.1 MT in 2020) and coastal fisheries targeting fishes (Atlantic herring Clupea harengus, Atlantic cod Gadus morhua), crustaceans (snow crab Chionoecetes opilio, rock crab Cancer irroratus, northern shrimp Pandalus borealis) and mollusks (whelk Buccinum sp, Arctic surf clam Mactromeris polynyma; Department of Fisheries and Oceans, 2019; Port de Sept-Îles, 2021).

2.2. Sources of anthropogenic activities

Relevant human activities were identified through a literature review and a compilation of data from local organizations (Port de Sept-Îles, Ville de Sept-Îles and Institut Nordique de Recherche en Environnement et en Santé au Travail) and integrative databases (Beauchesne et al., 2020). This data screening process identified seven relevant human activities: mussel aquaculture, dredging of sediment, city and industrial runoff, sewer discharge, commercial vessel movements and operations (shipping), artificial structures and commercial fisheries.

The distribution and intensity of human activities were characterized using R v4.2 and packages raster and sf (Pebesma, 2018; Hijmans, 2020; R Core Team, 2022). Data for anthropogenic sources consisted of spatial objects (multipoints, multilines and multipolygons), where the relative importance of each component was determined by comparing sources metadata, such as water discharge volume or number of fishing events, to grant standardized weighting coefficients. Data were handled according to confidentiality policies of each source provider.

2.3. Exposure of ecosystems to anthropogenic activities

Because characterization of vulnerability requires extensive data on the physiological responses of species and how influence translates to impact, we focused on the exposure of benthic communities to human activities (i.e., component Sj,x in the score by Halpern et al., 2019). We thus developed an index of exposure E for each considered human activity to describe the anthropogenic footprint in the study area. A “static” environment without spatial or temporal dynamics was considered, such that the index represents a “snapshot of exposure.” E was computed differently for land/sea-based activities and for fisheries, as explained below.

2.3.1. Land/sea-based human activities

Indices of exposure for aquaculture, dredging, runoff, sewers, shipping and structures were obtained using a diffusion model, implemented in the absence of a complete circulation model for the Baie des Sept Îles. We developed a unique model that uses theoretical particles set to diffuse within a defined area. These particles result from an activity (such as contaminants or sediment) introduced by point or line sources in the environment. The length of the journey from the source(s) of activity to a location D was used as a proxy of exposure: when D is low, particle density is high (being close to the source), thus indicating a high exposure of the ecosystem to this activity, and vice versa. We identified 11 sources of human activity in the study area, from punctual sources, e.g., sewer drains, to diffuse sources, e.g., coastal runoff from the city (Figure 2), acting as sources of particles in the diffusion model.

Figure 2.

Maps of the considered sources of land- and sea-based human activities in the study area. The human activities (in blue) are (A) aquaculture, (B) dredging, (C) runoff, (D) sewers, (E) structures, and (F) shipping.

Figure 2.

Maps of the considered sources of land- and sea-based human activities in the study area. The human activities (in blue) are (A) aquaculture, (B) dredging, (C) runoff, (D) sewers, (E) structures, and (F) shipping.

Close modal

Distance D was obtained using package gdistance (van Etten, 2017). A 100 × 100 m grid was created for the study area, where we established a connectivity matrix in a chess queen configuration (each cell to its eight direct neighbours using horizontal, vertical and diagonal directions). The cost of moving from one cell to another was computed with two constraints: (i) particles only diffuse in the marine environment and (ii) particles sink according to gravity and settle on the seafloor. To implement these aspects of the model, we used coastlines as boundaries (cost to select land cells is infinite) and bathymetry (movement of particles is primarily downward, while upward movement is secondary and hindered by topography) in the transition function. A least-cost pathfinding algorithm computed distance D from the source(s) of human activity to a specific grid cell (Dijkstra, 1959; van Etten, 2017).

Exposure indices E were calculated for each cell using D and a Gaussian kernel function (exponential quadratic relationship) to account for dispersion in a 2D environment while reducing the contribution of the highest values. The equation for E is thus:

where i is a cell, j is a human activity, A is the decay coefficient and r is the spatial extent of the grid.

As the considered human activities do not have the same relative influence on ecosystems, we simulated dispersion patterns with different relationships between distance and exposure by tuning a parameter in the function, the decay coefficient A. Five unique behaviour profiles for particles were established, from very localized (Type I) to ubiquitous (Type V; Figure S1). We performed a literature review to identify physical, chemical and biological pressures to assign behaviour profiles to human activities, considering a Sept-Îles context (Table 1). A pressure is defined here as a consequence of a driver (being natural or anthropogenic) affecting the ecosystem, following the Driver-Pressure-State-Impact-Response (DPSIR) framework (European Environmental Agency, 1999; Gari et al., 2015; Judd et al., 2015; Oesterwind et al., 2016). We queried articles and reviews from dedicated scientific studies about each pressure to obtain spatial and temporal ranges, allowing to select a profile type based on the decision table shown in Figure S2. Finally, the most prevalent behaviour profile between pressures was assigned to the corresponding human activity (Tables 1 and S1).

Table 1.

Description of human activities and related pressures considered in this study

Human ActivityPressuresDescriptionMain References
Exploitation of mussel farms Increase in organic matter concentration Introduced from mussel metabolism and related bacterial activity Christensen et al. (2003); Crawford et al. (2003); Richard et al. (2007); Callier et al. (2009); McKindsey et al. (2011); Heery et al. (2017); Lacoste et al. (2019)  
Modification of particulate matter Changes in composition of POM from farm operation and organism degradation Crawford et al. (2003); McKindsey et al. (2011); Wilding and Nickell (2013); Gallardi (2014)  
Increase in nutrient concentrations Related to metabolism of mussels and associated species Christensen et al. (2003); Cranford et al. (2003); McKindsey et al. (2011); Wilding and Nickell (2013); Gallardi (2014)  
Decrease in dissolved oxygen and sediment redox potential Related to metabolism of mussels, bacteria and associated species Wilding and Nickell (2013); Tičina et al. (2020)  
Introduction of shellfish diseases Inherent or emerging diseases from mussels or related organisms Tičina et al. (2020)  
Introduction of alien species Parasites, bacteria, viruses and other organisms linked to mussels Gallardi (2014); Tičina et al. (2020)  
Collection and dumping of sediment material Modification of sediment grain-size Related to currents and hydrodynamics induced by dredging operations Desprez (2000)  
Modification of sediment topography Changes in local sediment slope and bathymetry Desprez (2000); International Council for the Exploration of the Sea (2001)  
Resuspension of sediments Related to movement of dredged material Desprez (2000); International Council for the Exploration of the Sea (2001)  
Modification of chemical concentrations Chemical elements buried in sediment released back to water column International Council for the Exploration of the Sea (2001)  
Decrease in dissolved oxygen Related to metabolism induced by chemical release from sediment International Council for the Exploration of the Sea (2001)  
Transportation and destruction of organisms Related to sediment modification Desprez (2000); International Council for the Exploration of the Sea (2001)  
Runoff from city and industries Increase in nutrient concentrations Introduced by terrigenous inputs and biological activity from coastal settlements and facilities Müller et al. (2020)  
Increase in heavy metal concentrations Introduced by terrigenous inputs and biological activity from coastal settlements and facilities Müller et al. (2020)  
Increase in organic matter concentration Introduced by terrigenous inputs and biological activity from coastal settlements and facilities Müller et al. (2020)  
Modification of salinity gradients Related to freshwater inputs Müller et al. (2020)  
Wastewater and rainwater from sewers Increase in nutrient concentrations Introduced by wastewater management facilities Cotano and Villate (2006); Bertocci et al. (2019); Culhane et al. (2019)  
Increase in heavy metal concentrations Introduced by wastewater management facilities Bertocci et al. (2019); Culhane et al. (2019)  
Increase in organic matter concentration Introduced by wastewater management facilities and related to bacterial activity Cotano and Villate (2006); Bertocci et al. (2019); Culhane et al. (2019)  
Introduction of exogenic compounds (e.g., drugs) Introduced by wastewater management facilities Islam and Tanaka (2004); Culhane et al. (2019)  
Modification of particulate matter Related to solid and dissolved matter in wastewater outputs Oviatt et al. (1987)  
Decrease in dissolved oxygen Related to bacterial activity Culhane et al. (2019)  
Modification of temperature gradients Related to freshwater inputs Bertocci et al. (2019); Culhane et al. (2019)  
Modification of salinity gradients Related to freshwater inputs Bertocci et al. (2019); Culhane et al. (2019)  
Increase in biological activity Bacteria, viruses and other organisms present in wastewaters Islam and Tanaka (2004); Müller et al. (2020)  
Artificial structures (piers, marina) Modification of hydrodynamics Related to addition of solid structures Bulleri and Chapman (2010); Heery et al. (2017)  
Increase in anthropogenic noise Produced from shipping operations and machinery Heery et al. (2017)  
Increase in artificial light Produced by dedicated structures Bulleri and Chapman (2010); Heery et al. (2017)  
Modification of electromagnetic fields Produced by underwater cables Heery et al. (2017)  
Increase in turbidity Related to sediment resuspension from increased hydrodynamics Mineur et al. (2012); Heery et al. (2017)  
Introduction of exogenic compounds Dilution from paints or materials used Heery et al. (2017)  
Increase in organic matter concentration Related to accumulation of sediments Heery et al. (2017)  
Introduction of alien species Settlement rate increased on solid structures Bulleri and Chapman (2010); Mineur et al. (2012); Heery et al. (2017)  
Modification of species communities Linked to new habitats provided by solid structures and to changes in connectivity Bulleri and Chapman (2010); Bishop et al. (2017); Heery et al. (2017); Momota and Hosokawa (2021)  
Commercial vessels anchoring, movement and operation Introduction of exogenic compounds Dilution from hull paints Jägerbrand et al. (2019); Byrnes and Dunn (2020)  
Increase in hydrocarbons Produced from onboard systems, potential spills and cargo hauling operations Jägerbrand et al. (2019); Byrnes and Dunn (2020)  
Increase in heavy metal concentrations Linked to cargo hauling operations Jägerbrand et al. (2019); Byrnes and Dunn (2020)  
Introduction of marine litter Plastics and other wastes Jägerbrand et al. (2019); Byrnes and Dunn (2020)  
Increase in coastal erosion Disturbance of coastal topography by waves Jägerbrand et al. (2019); Byrnes and Dunn (2020)  
Increase in anthropogenic noise Produced from engine and ship machinery operation Jägerbrand et al. (2019)  
Increase in artificial light Produced from onboard systems Jägerbrand et al. (2019); Byrnes and Dunn (2020)  
Introduction of alien species Non-indigenous species present in ballast waters or fouling external surfaces Jägerbrand et al. (2019); Byrnes and Dunn (2020)  
Modification of sediment topography Use of anchor systems Davis et al. (2018); Jägerbrand et al. (2019); Byrnes and Dunn (2020)  
 Resuspension of sediments Use of anchor systems Davis et al. (2018); Byrnes and Dunn (2020); World Wildlife Fund (2020)  
Human ActivityPressuresDescriptionMain References
Exploitation of mussel farms Increase in organic matter concentration Introduced from mussel metabolism and related bacterial activity Christensen et al. (2003); Crawford et al. (2003); Richard et al. (2007); Callier et al. (2009); McKindsey et al. (2011); Heery et al. (2017); Lacoste et al. (2019)  
Modification of particulate matter Changes in composition of POM from farm operation and organism degradation Crawford et al. (2003); McKindsey et al. (2011); Wilding and Nickell (2013); Gallardi (2014)  
Increase in nutrient concentrations Related to metabolism of mussels and associated species Christensen et al. (2003); Cranford et al. (2003); McKindsey et al. (2011); Wilding and Nickell (2013); Gallardi (2014)  
Decrease in dissolved oxygen and sediment redox potential Related to metabolism of mussels, bacteria and associated species Wilding and Nickell (2013); Tičina et al. (2020)  
Introduction of shellfish diseases Inherent or emerging diseases from mussels or related organisms Tičina et al. (2020)  
Introduction of alien species Parasites, bacteria, viruses and other organisms linked to mussels Gallardi (2014); Tičina et al. (2020)  
Collection and dumping of sediment material Modification of sediment grain-size Related to currents and hydrodynamics induced by dredging operations Desprez (2000)  
Modification of sediment topography Changes in local sediment slope and bathymetry Desprez (2000); International Council for the Exploration of the Sea (2001)  
Resuspension of sediments Related to movement of dredged material Desprez (2000); International Council for the Exploration of the Sea (2001)  
Modification of chemical concentrations Chemical elements buried in sediment released back to water column International Council for the Exploration of the Sea (2001)  
Decrease in dissolved oxygen Related to metabolism induced by chemical release from sediment International Council for the Exploration of the Sea (2001)  
Transportation and destruction of organisms Related to sediment modification Desprez (2000); International Council for the Exploration of the Sea (2001)  
Runoff from city and industries Increase in nutrient concentrations Introduced by terrigenous inputs and biological activity from coastal settlements and facilities Müller et al. (2020)  
Increase in heavy metal concentrations Introduced by terrigenous inputs and biological activity from coastal settlements and facilities Müller et al. (2020)  
Increase in organic matter concentration Introduced by terrigenous inputs and biological activity from coastal settlements and facilities Müller et al. (2020)  
Modification of salinity gradients Related to freshwater inputs Müller et al. (2020)  
Wastewater and rainwater from sewers Increase in nutrient concentrations Introduced by wastewater management facilities Cotano and Villate (2006); Bertocci et al. (2019); Culhane et al. (2019)  
Increase in heavy metal concentrations Introduced by wastewater management facilities Bertocci et al. (2019); Culhane et al. (2019)  
Increase in organic matter concentration Introduced by wastewater management facilities and related to bacterial activity Cotano and Villate (2006); Bertocci et al. (2019); Culhane et al. (2019)  
Introduction of exogenic compounds (e.g., drugs) Introduced by wastewater management facilities Islam and Tanaka (2004); Culhane et al. (2019)  
Modification of particulate matter Related to solid and dissolved matter in wastewater outputs Oviatt et al. (1987)  
Decrease in dissolved oxygen Related to bacterial activity Culhane et al. (2019)  
Modification of temperature gradients Related to freshwater inputs Bertocci et al. (2019); Culhane et al. (2019)  
Modification of salinity gradients Related to freshwater inputs Bertocci et al. (2019); Culhane et al. (2019)  
Increase in biological activity Bacteria, viruses and other organisms present in wastewaters Islam and Tanaka (2004); Müller et al. (2020)  
Artificial structures (piers, marina) Modification of hydrodynamics Related to addition of solid structures Bulleri and Chapman (2010); Heery et al. (2017)  
Increase in anthropogenic noise Produced from shipping operations and machinery Heery et al. (2017)  
Increase in artificial light Produced by dedicated structures Bulleri and Chapman (2010); Heery et al. (2017)  
Modification of electromagnetic fields Produced by underwater cables Heery et al. (2017)  
Increase in turbidity Related to sediment resuspension from increased hydrodynamics Mineur et al. (2012); Heery et al. (2017)  
Introduction of exogenic compounds Dilution from paints or materials used Heery et al. (2017)  
Increase in organic matter concentration Related to accumulation of sediments Heery et al. (2017)  
Introduction of alien species Settlement rate increased on solid structures Bulleri and Chapman (2010); Mineur et al. (2012); Heery et al. (2017)  
Modification of species communities Linked to new habitats provided by solid structures and to changes in connectivity Bulleri and Chapman (2010); Bishop et al. (2017); Heery et al. (2017); Momota and Hosokawa (2021)  
Commercial vessels anchoring, movement and operation Introduction of exogenic compounds Dilution from hull paints Jägerbrand et al. (2019); Byrnes and Dunn (2020)  
Increase in hydrocarbons Produced from onboard systems, potential spills and cargo hauling operations Jägerbrand et al. (2019); Byrnes and Dunn (2020)  
Increase in heavy metal concentrations Linked to cargo hauling operations Jägerbrand et al. (2019); Byrnes and Dunn (2020)  
Introduction of marine litter Plastics and other wastes Jägerbrand et al. (2019); Byrnes and Dunn (2020)  
Increase in coastal erosion Disturbance of coastal topography by waves Jägerbrand et al. (2019); Byrnes and Dunn (2020)  
Increase in anthropogenic noise Produced from engine and ship machinery operation Jägerbrand et al. (2019)  
Increase in artificial light Produced from onboard systems Jägerbrand et al. (2019); Byrnes and Dunn (2020)  
Introduction of alien species Non-indigenous species present in ballast waters or fouling external surfaces Jägerbrand et al. (2019); Byrnes and Dunn (2020)  
Modification of sediment topography Use of anchor systems Davis et al. (2018); Jägerbrand et al. (2019); Byrnes and Dunn (2020)  
 Resuspension of sediments Use of anchor systems Davis et al. (2018); Byrnes and Dunn (2020); World Wildlife Fund (2020)  

2.3.2. Fisheries

The exposure index E for fisheries was calculated by considering the number of fishing events by gear type: areas with a high number of events indicate a high exposure, and vice versa. Data were extracted from the eDrivers platform in the industrial harbour area of Sept-Îles for events recorded between 2010 and 2015 (Beauchesne et al., 2020). Fishing events were compiled in a raster file for four types of fishing gear: traps, bottom-trawls, nets and dredges. We averaged the number of events to obtain a proxy of fishing intensity per gear G. We obtained the exposure index E by combining G from the four gear types using the following equation:

where i is a cell and k is a gear type.

2.3.3. Cumulative exposure

The seven indices of exposure described above were standardized between 0 (lowest exposure) and 1 (highest exposure) then summed to provide a cumulative exposure score C. Because the relative importance between human activities is unknown for our study area, we considered each activity in the cumulative score to be of equivalent importance (i.e., no weighting parameters). The equation for C is thus:

where i is a cell and j is a human activity.

2.4. Habitat and biological samples

To understand how the calculated exposure indices relate to benthic ecosystems, we used several ecological datasets in the study area from Dreujou et al. (2021). Samples were collected in July 2017 within the bay and archipelago of Sept-Îles, where a total of 108 stations were selected using a semi-randomization algorithm. Stations were constrained between depths of 0 m and 80 m with an increased sampling effort in areas where sources of human activities were present (Figure 1).

Station depth was obtained from a navigation sonar, then corrected with respect to tide height at the time of sampling. A Ponar grab (0.05 m2) was deployed at each station from a boat with two independent casts. This first cast collected sediment samples for chemical and physical analyses while sediment from the second cast was sieved on a 0.5 mm mesh size for macrofauna identification. This approach allowed the collection of data on organic matter content, sediment grain-size percentages (gravel, sand, silt, clay), heavy-metal concentrations (arsenic, cadmium, chromium, copper, iron, manganese, mercury, lead, zinc), and the density and wet biomass for each taxon identified.

2.5. Statistical analyses

Correlation between exposure indices and habitat parameters was assessed using Spearman’s rank coefficient (Quinn and Keough, 2002). Building on the concept of the Ecological Quality Ratio by van de Bund and Solimini (2007), we computed an exposure ratio ER at each sampled station based on their cumulative exposure score. This ratio compares the value of the score to its extrema to provide integrative information on the severity of the cumulative exposure:

where Ci is the cumulative exposure score at station i, Rhigh is the reference value for a high exposure (here, 7) and Rlow is the reference value for a low exposure (here, 0). Using ER, we assigned an exposure status to each station: a low status corresponds to a high exposure ratio and a high status to a low exposure ratio. Five categories were defined for this status, similar to those used for the Ecological Quality Status: “bad,” “low,” “moderate,” “good” and “high” status (van de Bund and Solimini, 2007).

Characteristic taxa were assessed for each exposure status by computing the indicator value score on benthic assemblages (IndVal, 1000 randomization iterations; Dufrêne and Legendre, 1997). We used benthic community descriptors and environmental indicators calculated by Dreujou et al. (2021): total density, total biomass, taxa richness, Shannon index (base e logarithm), Pielou evenness, Multivariate AZTI Marine Biotic Index (M-AMBI), BENTIX score and Benthic Opportunistic Polychaete/Amphipod ratio (BOPA; Legendre and Legendre, 1998; Simboura and Zenetos, 2002; Dauvin and Ruellet, 2007; Muxika et al., 2007; Magurran and McGill, 2011; Dauvin et al., 2016). M-AMBI, BENTIX and BOPA are calculated using relative abundance of species groups, established based on tolerance to perturbation. Phylum mean density and mean biomass were calculated for each exposure status to evaluate their relative variation.

Relationships between exposure indices (predictors) and benthic descriptors (independent variables) were evaluated using multiple regression models. Variables were transformed (log(x + 1) or square root) if the assumptions of normality and homoscedasticity were not respected (Quinn and Keough, 2002). We also explored relationships between the taxa assemblage and exposure indices using non-parametric multivariate regression with distance-based linear modelling (DistLM, 9,999 permutations; McArdle and Anderson, 2001). In both regression analyses, we added depth as a covariate to account for bathymetric variation between stations. Statistical analyses were done using R v4.2 with package vegan and PRIMER-E v6 software (Clarke and Gorley, 2006; Oksanen et al., 2022; R Core Team, 2022).

3.1. Exposure indices

Our literature review highlighted 25 unique pressures linked to human activities in the Baie des Sept Îles (Table 1), including biological, physical and chemical effects on ecosystems. Most pressures were associated with a localized profile (Type II, 17), followed by diffused (Type III, 12) and very diffused (Type IV, 9) profiles.

Overall, bay-wide average exposure indices were low to moderate, varying between 0.05 (fisheries) and 0.6 (sewers). Only stations close to sources of activity presented high index values, with a limited area of influence for each human activity (Figure 3). The fisheries exposure index was highly localized, similar to those for aquaculture and dredging (Figure 3). Sewers had the most extensive footprint in the bay, with a highest exposure index at the intersection of the zone of influence of multiple sewer drains. Exposure due to runoff and shipping both decreased with increasing distance from their sources (Figure 3). Correlation tests on exposure indices showed positive relationships between runoff and structures (correlation coefficient ρ = 0.87), dredging and runoff (ρ = 0.76), sewers and shipping (ρ = 0.74); all other pairwise tests showed low to moderate relationships (0.28 < |ρ| < 0.7).

Figure 3.

Exposure indices calculated in the study area for each human activity. Histograms represent the number of stations along the value of the index, and colours correspond to exposure classes, from low (0.0) to high (1.0). The human activities considered are (A) aquaculture, (B) dredging, (C) runoff, (D) sewers, (E) structures, (F) shipping, and (G) fisheries.

Figure 3.

Exposure indices calculated in the study area for each human activity. Histograms represent the number of stations along the value of the index, and colours correspond to exposure classes, from low (0.0) to high (1.0). The human activities considered are (A) aquaculture, (B) dredging, (C) runoff, (D) sewers, (E) structures, (F) shipping, and (G) fisheries.

Close modal

The cumulative exposure score calculated at each station varied between 0.383 and 3.698 across the bay, with an average of 1.865 (standard error of 0.08, theoretical maximum value of 7; Figure 4A). The highest values were detected close to the main industrial and urban activity sources, especially in front of the city of Sept-Îles and the Pointe-Noire sector (Figure 4A). The vast majority of the stations sampled for macrofauna were assigned an exposure status of “high” (n = 40) or “good” (n = 50), while 18 presented a “moderate” status (Figure 4B).

Figure 4.

Values of the cumulative exposure score in the study area. (A) Bay-scale variability of the score and (B) exposure status obtained at each station based on this score.

Figure 4.

Values of the cumulative exposure score in the study area. (A) Bay-scale variability of the score and (B) exposure status obtained at each station based on this score.

Close modal

3.2. Relationships with benthic ecosystems

Correlations between exposure indices and habitat parameters are summarized in Table 2. Overall, sediment parameters had lower correlations with exposure indices than did heavy metals. Organic matter showed positive relationships with sewers (ρ = 0.42) and shipping (ρ = 0.37) and a negative one with fisheries (ρ = −0.49). Sand and silt presented relationships for fisheries, runoff, sewers and shipping, with low to moderate coefficients (0.19 < |ρ| < 0.38). Notably, gravel was poorly related to fisheries and shipping, while all correlations for clay were non-significant. Concerning heavy metals, three patterns were observed: aquaculture and fisheries had moderate to high negative coefficients for all heavy metals (–0.24 < ρ < –0.61); dredging and structures had moderate positive coefficients for some metals (0.23 < ρ < 0.57); and sewers and shipping had moderate to high positive coefficients for all metals (0.27 < ρ < 0.75). For environmental indicators, most coefficients were not significant but some low relationships (|ρ| < 0.25) were evident. Finally, the cumulative exposure score showed low to moderate positive relationships with organic matter and most heavy metals (0.19 < ρ < 0.58).

Table 2.

Spearman rank correlation coefficients between human activity exposure scores and ecological variablesa

Human Activitiesb
Ecological VariablesAqDrFiRuSeShStCE
Sediment parameters 
 Organic matter — 0.24 –0.49 — 0.42 0.37 — 0.28 
 Gravel — — 0.2 — — –0.22 — — 
 Sand — — 0.38 0.3 –0.39 –0.23 — — 
 Silt — — –0.38 –0.19 0.34 0.24 — — 
 Clay — — — — — — — — 
Heavy metal concentrations 
 Arsenic –0.28 — –0.57 — 0.64 0.42 — 0.19 
 Cadmium –0.24 — –0.54 — 0.59 0.27 — — 
 Chromium –0.34 0.23 –0.55 — 0.69 0.46 0.28 0.36 
 Copper –0.48 0.37 –0.61 0.3 0.75 0.57 0.47 0.57 
 Iron –0.48 0.57 –0.58 0.28 0.67 0.57 0.51 0.58 
 Manganese –0.48 0.37 –0.59 — 0.76 0.55 0.4 0.47 
 Mercury –0.27 — –0.54 — 0.64 0.41 — 0.2 
 Lead –0.27 — –0.56 — 0.71 0.42 — 0.3 
 Zinc –0.41 0.28 –0.61 — 0.74 0.54 0.33 0.45 
Environmental indicators 
 M-AMBI — — 0.2 — — — — — 
 BENTIX — — –0.23 — 0.21 0.23 — — 
 BOPA –0.19 — — — 0.25 — — — 
Human Activitiesb
Ecological VariablesAqDrFiRuSeShStCE
Sediment parameters 
 Organic matter — 0.24 –0.49 — 0.42 0.37 — 0.28 
 Gravel — — 0.2 — — –0.22 — — 
 Sand — — 0.38 0.3 –0.39 –0.23 — — 
 Silt — — –0.38 –0.19 0.34 0.24 — — 
 Clay — — — — — — — — 
Heavy metal concentrations 
 Arsenic –0.28 — –0.57 — 0.64 0.42 — 0.19 
 Cadmium –0.24 — –0.54 — 0.59 0.27 — — 
 Chromium –0.34 0.23 –0.55 — 0.69 0.46 0.28 0.36 
 Copper –0.48 0.37 –0.61 0.3 0.75 0.57 0.47 0.57 
 Iron –0.48 0.57 –0.58 0.28 0.67 0.57 0.51 0.58 
 Manganese –0.48 0.37 –0.59 — 0.76 0.55 0.4 0.47 
 Mercury –0.27 — –0.54 — 0.64 0.41 — 0.2 
 Lead –0.27 — –0.56 — 0.71 0.42 — 0.3 
 Zinc –0.41 0.28 –0.61 — 0.74 0.54 0.33 0.45 
Environmental indicators 
 M-AMBI — — 0.2 — — — — — 
 BENTIX — — –0.23 — 0.21 0.23 — — 
 BOPA –0.19 — — — 0.25 — — — 

aOnly significant relationships are presented.

bAquaculture (Aq), dredging (Dr), fisheries (Fi), runoff (Ru), sewers (Se), shipping (Sh), structures (St), cumulative exposure (CE), Multivariate AZTI Marine Biotic Index (M-AMBI), Benthic Opportunistic Polychaete/Amphipod ratio (BOPA).

The analysis of phylum mean density and mean biomass varied by exposure status (Figure 5). The biomass of annelids was greatest for lower status classes, increasing from 3.7% in “high” status to 38.9% in “moderate” status, while density stayed between 28% and 40% of the total community. Arthropod density peaked at “good” status (53.3%) and was also quite high at “high” and “moderate” status (30.3% and 40.2%, respectively), whereas biomass proportion did not exceed 3%. Mollusc biomass also peaked at “good” status (40.8%) then dropped to 16.9% at “moderate” status and 8.6% at “high” status. Mollusc density stayed around 15% for each status. Echinoderm biomass and nematode density were highest at “high” status stations (86% and 20.3%, respectively) then dropped significantly at “moderate” status stations (as low as 32.7% and 0%, respectively).

Figure 5.

Proportions of each phylum for the five exposure statuses. Proportions are based on (A) mean density and (B) mean biomass of each phylum (colour-coded, inset legend). Numbers on top of each bar correspond to the number of stations in each status.

Figure 5.

Proportions of each phylum for the five exposure statuses. Proportions are based on (A) mean density and (B) mean biomass of each phylum (colour-coded, inset legend). Numbers on top of each bar correspond to the number of stations in each status.

Close modal

The calculation of IndVal yielded few significant characteristic taxa for each exposure status. Stations with a “high” status presented four characteristic taxa: Nematoda (p < 0.001), the gastropod Ameritella agilis (p = 0.011), Nephtyidae polychaetes (p = 0.021) and the amphipod Byblis gaimardii (p = 0.033). Only harpacticoid copepods were identified as significant for stations with a “good” status (p = 0.029), while no characteristic taxa were related to stations with a “moderate” status.

Multiple regression analyses between exposure indices and benthic descriptors showed that predictive power was highest for Shannon index (adjusted R2 = 0.29) followed by taxa richness (adjusted R2 = 0.16) and Pielou evenness (adjusted R2 = 0.14). Predictive power for total density and biomass was quite low (adjusted R2 < 0.04) (Table 3). Marginal tests for depth indicated statistically significant relationships for Shannon index (standardized coefficient = 0.54), Pielou evenness (0.46) and taxa richness (0.25). Exposure indices mostly showed positive effects on total density, taxa richness, Shannon index and Pielou evenness, with coefficients seldom going over 0.15 (Table 3). Different patterns were detected for total biomass, where most activities had negative effects, especially sewers (–0.58) and runoff (–0.48), while structures presented a strong positive effect (0.54). Post-analysis diagnostics for homoscedasticity, normality of residuals and independence were quite robust, especially for taxa richness, Shannon index and Pielou evenness.

Table 3.

Standardized predictor coefficients (and standard error) from multiple linear regression models between depth and human activity exposure indices and benthic community descriptors

Human Activitiesb
Community DescriptorsaDepthAqDrFiRuSeShStR2adjc
N –0.19 (0.11) 0.05 (0.14) –0.12 (0.13) 0.12 (0.12) 0.19 (0.22) 0.15 (0.19) –0.09 (0.13) –0.18 (0.25) 0.02 
p-value 0.1024 0.6995 0.3483 0.2918 0.3952 0.4198 0.4884 0.4879  
B –0.21 (0.11) –0.25 (0.13) –0.01 (0.13) –0.1 (0.11) –0.48 (0.22) –0.58 (0.18) 0.09 (0.13) 0.54 (0.25) 0.04 
p-value 0.0659 0.0598 0.9393 0.3936 0.0301d 0.0022 0.4761 0.034  
S 0.25 (0.1) 0.14 (0.12) –0.18 (0.12) 0.2 (0.1) 0.26 (0.2) –0.16 (0.17) 0.24 (0.12) –0.15 (0.23) 0.20 
p-value 0.0172 0.2602 0.1297 0.0531 0.1994 0.3362 0.0459 0.5159  
H 0.54 (0.1) 0.17 (0.12) 0.01 (0.11) 0.03 (0.1) 0.41 (0.19) 0.01 (0.16) 0.11 (0.11) –0.34 (0.22) 0.29 
p-value <0.0001 0.1499 0.9584 0.7924 0.0287 0.9877 0.3289 0.1150  
J 0.46 (0.11) 0.07 (0.13) 0.14 (0.12) –0.13 (0.11) 0.32 (0.21) 0.07 (0.17) –0.05 (0.12) –0.33 (0.23) 0.14 
p-value <0.0001 0.6073 0.2627 0.2378 0.1236 0.7013 0.7063 0.1715  
Human Activitiesb
Community DescriptorsaDepthAqDrFiRuSeShStR2adjc
N –0.19 (0.11) 0.05 (0.14) –0.12 (0.13) 0.12 (0.12) 0.19 (0.22) 0.15 (0.19) –0.09 (0.13) –0.18 (0.25) 0.02 
p-value 0.1024 0.6995 0.3483 0.2918 0.3952 0.4198 0.4884 0.4879  
B –0.21 (0.11) –0.25 (0.13) –0.01 (0.13) –0.1 (0.11) –0.48 (0.22) –0.58 (0.18) 0.09 (0.13) 0.54 (0.25) 0.04 
p-value 0.0659 0.0598 0.9393 0.3936 0.0301d 0.0022 0.4761 0.034  
S 0.25 (0.1) 0.14 (0.12) –0.18 (0.12) 0.2 (0.1) 0.26 (0.2) –0.16 (0.17) 0.24 (0.12) –0.15 (0.23) 0.20 
p-value 0.0172 0.2602 0.1297 0.0531 0.1994 0.3362 0.0459 0.5159  
H 0.54 (0.1) 0.17 (0.12) 0.01 (0.11) 0.03 (0.1) 0.41 (0.19) 0.01 (0.16) 0.11 (0.11) –0.34 (0.22) 0.29 
p-value <0.0001 0.1499 0.9584 0.7924 0.0287 0.9877 0.3289 0.1150  
J 0.46 (0.11) 0.07 (0.13) 0.14 (0.12) –0.13 (0.11) 0.32 (0.21) 0.07 (0.17) –0.05 (0.12) –0.33 (0.23) 0.14 
p-value <0.0001 0.6073 0.2627 0.2378 0.1236 0.7013 0.7063 0.1715  

aTotal density (N), total biomass (B), taxa richness (S), Shannon index (H), Pielou evenness (J).

bAquaculture (Aq), dredging (Dr), fisheries (Fi), runoff (Ru), sewers (Se), shipping (Sh), structures (St).

cAdjusted R squared (R2adj) of the model.

dSignificant p-values of marginal tests on predictors are highlighted in bold.

DistLM regression on the taxa assemblage had an R2 of 0.22, and the ancillary constrained ordination is shown on Figure 6. The first two axes explained 14.9% of the variance, where two negatively correlated groups of exposure indices were obtained: sewers/shipping to one side, aquaculture/fisheries/runoff to the other (dredging and structures had lower influence). Interestingly, depth did not correlate with exposure indices, although it correlated strongly with benthic community structure, as expected. No structure may be described based on the similarity between stations, except for a group of stations with a higher cumulative exposure score at the centre of the biplot (Figure 6).

Figure 6.

Constrained ordination with a distance-based Redundancy Analysis on the taxa assemblage. Axes percentages are the proportion of variance explained, and colours correspond to the cumulative exposure score. Grey lines indicate depth and human activities: aquaculture (Aq), dredging (Dr), fisheries (Fi), runoff (Ru), sewers (Se), shipping (Sh), and structures (St).

Figure 6.

Constrained ordination with a distance-based Redundancy Analysis on the taxa assemblage. Axes percentages are the proportion of variance explained, and colours correspond to the cumulative exposure score. Grey lines indicate depth and human activities: aquaculture (Aq), dredging (Dr), fisheries (Fi), runoff (Ru), sewers (Se), shipping (Sh), and structures (St).

Close modal

This study presents a modelling framework to study anthropogenic influences on coastal benthic ecosystems at a local scale of <100 km using proxies such as distance from sources of activity and fishing events. This tool represents a relevant addition for environmental managers by providing a way to estimate exposure to human activities that may be related to existing ecological data and by covering a spatial resolution that is usually less investigated in the cumulative impact literature.

We detected patterns between benthic communities and exposure indices for seven human activities: mussel aquaculture, dredging of sediment, runoff from city and industries, sewer discharge, commercial vessels movement and operation, artificial structures and commercial fisheries. When studying benthic species assemblages, characteristic species and phylum composition changed along a cumulative exposure gradient. The most striking result is an increase of annelid mean biomass in stations with a higher cumulative exposure score, with the dominance of Nephtys insica, Praxillella praetermissa and Maldanidae spp. Many studies have highlighted the use of certain species to indicate ecological status, providing useful indicators on the state of disturbance (Pearson and Rosenberg, 1978; Grall and Glémarec, 1997; Borja et al., 2000). When classifying species according to Borja et al. (2000), abundance of species sensitive to disturbance (Type I) and tolerant to disturbance (Type III) varied between “high” status and “moderate” status stations. However, no significant trend could be detected because of very different sample sizes between conditions. These results suggest a certain link between the station classification based on their cumulative exposure and the detection of possible disturbance effects (highlighted as a proxy of organic matter enrichment; Pearson and Rosenberg, 1978; Borja et al., 2000). An increased sampling effort and the consideration of other pressures would be needed to strengthen this conclusion.

Models explaining community characteristics (i.e., species richness) by environmental variables and exposure indices provided valuable information on anthropogenic influences on benthic communities. While the predictive power of these regressions was moderate with relatively low coefficients, depth was a significant predictor of community composition. This finding is coherent with patterns of biodiversity in marine ecosystems that distribute along a gradient from richest to poorest with increasing depth (Gray and Elliott, 2009; Levinton, 2013; Piacenza et al., 2015). These results advocate for a complementarity of covariates from abiotic and anthropogenic sources in the study of anthropogenic influence, where the explained variance increases when both types of predictors are considered. That being said, community characteristics may not be the best descriptors to properly account for disturbance, in particular because of their univariate nature (Drouin et al., 2011). Previous works in Baie des Sept Îles by Carrière (2018), Dreujou et al. (2021) and Ferrario et al. (2022) used various methods, including ecological indicators to assess the ecological status of the region. These studies showed that the overall ecological status of the bay is high, which may represent another hypothesis to explain the lack of a strong disturbance gradient in the area.

All sampled stations were influenced by at least two different sources of exposure, which reinforces the importance of studying the cumulative effects of human activities in an integrative way (Dreujou et al., 2020a; Carrier-Belleau et al., 2021). Classification based on the cumulative exposure score showed that 83% of the sampled stations were assigned a “high” to “good” status (i.e., C < 2.8). This result suggests a relatively low anthropogenic influence in the majority of the Baie des Sept Îles, which is coherent with the state of the benthic communities described above. Stations reaching higher cumulative exposure scores (“moderate,” “poor” and “bad” status) can be considered anthropogenic “hotspots” (Crain et al., 2008; Darling and Côté, 2008; Côté et al., 2016; Galic et al., 2018), as co-occurring human activities have an increased probability to produce possible emergent effects. Hotspot identification may guide environmental protection and sustainable development as a way to target conservation areas or to prioritize management resources where they will be the most impactful.

Overall, exposure indices showed a moderate relationship with community composition in multivariate regression models, except concerning a group of stations with a high cumulative exposure score where communities were similar. Interestingly, the same predictive power was obtained when using only environmental parameters (R2 = 0.24), which reinforces the pertinence of the exposure scores to explain benthic communities as complements of the environmental variables. Dominant taxa at stations with high cumulative exposure scores included Bipalponephtys neotena and Macoma calcarea, species tolerant to disturbance and found within disturbed areas by Dreujou et al. (2020b). These assemblages may indicate a disturbed profile, where sensitive taxa are rare relative to the dominance of disturbance-tolerant taxa without opportunistic species (Pearson and Rosenberg, 1978; Grall and Glémarec, 1997), reinforcing previous conclusions. Most importantly, hotspots of cumulative exposure are located in areas where communities have a moderately disturbed profile, as detected by Dreujou et al. (2020b; stations from cluster A). These results describe a “snapshot” of ecosystem exposure as they do not include seasonal or temporal variation, which may drastically modify the structure of benthic communities (Dreujou et al., 2018). Establishing long-term monitoring protocols will thus be necessary to increase the robustness of the analysis. Furthermore, relative weight of anthropogenic influence and the integration of the functioning of the ecosystem need to be carefully addressed when calculating cumulative indices, using expert opinion and dedicated research.

This study proposes a means to model exposure using relatively few environmental data, mainly spatial information and surveys with local stakeholders (types of human activity present, location and intensity of sources), which is well suited to describe anthropogenic influence in areas where historical ecological data are scarce. While the score from Halpern et al. (2008) is useful for characterizing cumulative effects globally, the global scale of this assessment prevented the use of fine-scale local data from which an environmental assessment would benefit. Focusing on gradients of exposure is promising as it allows the quantification of anthropogenic influences without needing to define reference conditions, which are often biased due to a lack of historical data or pristine ecosystems (Borja et al., 2012; Korpinen and Andersen, 2016). The addition of other relevant activities, such as tourism or recreational boating, or environmental drivers (e.g., freshwater or terrigenous inputs) would greatly increase its general applicability to other regions. Future works should also consider emerging trends, such as antagonistic or synergistic effects (Korpinen and Andersen, 2016; Galic et al., 2018; Carrier-Belleau et al., 2021), as they influence the description of vulnerability and ecosystem responses to perturbation, thus drastically impacting environmental assessment outcomes.

Our results contribute to a better understanding of sub-Arctic ecosystems and how multiple human activities may influence them, which is of tremendous importance in the context of global climate change. By including ecological data at a high spatial resolution, we propose a tool that allows direct evaluation of anthropogenic exposure and its possible effects on benthic communities. Such initiatives could increase the efficiency of adaptive environmental management, especially where data are still scarce, by being incorporated within standardized frameworks that are applicable to a wide range of coastal regions.

Data and code used for this study are available online at https://doi.org/10.5683/SP3/5XZKP3.

The supplemental files for this article can be found as follows:

Figure S1, Figure S2, Table S1.DOCX

The authors would like to acknowledge and thank all the people that helped during the field campaigns, lab work and data analysis: Claudy Dechêsne, Serge Gallienne, Cindy Grant, Lisa Tréau de Coeli, Laure de Montety, Philippe-Olivier Dumais, Raphaël Bouchard, Simon Bélanger, Raphaël Mabit, and Carlos Araujo. They also thank Andrew MacDonald for his highly constructive comments on the early development of the pathfinding algorithm. This study is a contribution to the research program of FRQNT Strategic Network Québec-Océan, which also provided scientific support.

This research was sponsored by NSERC Canadian Healthy Oceans Network (CHONe) and its partners: Fisheries and Oceans Canada and INREST (representing the Port of Sept-Îles and City of Sept-Îles), by the FRQNT Strategic Network Québec-Océan and by the Joint CNRS-UL Laboratory Takuvik to PA.

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Contributed to conception and design: ED, DB, RMD, CWM, PA.

Contributed to acquisition of data: ED.

Contributed to analysis and interpretation of data: ED, DB, RMD, FN.

Drafted and/or revised the article: ED, DB, RMD, JC, FN, CWM, PA.

Approved the submitted version for publication: ED, DB, RMD, JC, FN, CWM, PA.

Andersen
,
JH
,
Halpern
,
BS
,
Korpinen
,
S
,
Murray
,
C
,
Reker
,
J.
2015
.
Baltic Sea biodiversity status vs. cumulative human pressures
.
Estuarine, Coastal and Shelf Science
161
:
88
92
. DOI: http://dx.doi.org/10.1016/j.ecss.2015.05.002.
Ban
,
NC
,
Alidina
,
HM
,
Ardron
,
JA.
2010
.
Cumulative impact mapping: Advances, relevance and limitations to marine management and conservation, using Canada’s Pacific waters as a case study
.
Marine Policy
34
(
5
):
876
886
. DOI: http://dx.doi.org/10.1016/j.marpol.2010.01.010.
Beauchesne
,
D
,
Daigle
,
RM
,
Vissault
,
S
,
Gravel
,
D
,
Bastien
,
A
,
Bélanger
,
S
,
Bernatchez
,
P
,
Blais
,
M
,
Bourdages
,
H
,
Chion
,
C
,
Galbraith
,
PS
,
Halpern
,
BS
,
Lavoie
,
C
,
McKindsey
,
CW
,
Mucci
,
A
,
Pineault
,
S
,
Starr
,
M
,
Ste-Marie
,
AS
,
Archambault
,
P.
2020
.
Characterizing exposure to and sharing knowledge of drivers of environmental change in the St. Lawrence system in Canada
.
Frontiers in Marine Science
7
:
1
14
. DOI: http://dx.doi.org/10.3389/fmars.2020.00383.
Beauchesne
,
D
,
Grant
,
C
,
Gravel
,
D
,
Archambault
,
P.
2016
.
L’évaluation des impacts cumulés dans l’estuaire et le golfe du Saint-Laurent: Vers une planification systémique de l’exploitation des ressources
.
Le Naturaliste canadien
140
(
2
):
45
55
. DOI: http://dx.doi.org/10.7202/1036503ar.
Beermann
,
AJ
,
Elbrecht
,
V
,
Karnatz
,
S
,
Ma
,
L
,
Matthaei
,
CD
,
Piggott
,
JJ
,
Leese
,
F.
2018
.
Multiple-stressor effects on stream macroinvertebrate communities: A mesocosm experiment manipulating salinity, fine sediment and flow velocity
.
Science of The Total Environment
610–611
:
961
971
. DOI: http://dx.doi.org/10.1016/j.scitotenv.2017.08.084.
Berkes
,
F
,
Folke
,
C
,
Colding
,
J.
2000
.
Linking social and ecological systems: Management practices and social mechanisms for building resilience
.
Cambridge, UK
:
Cambridge University Press
.
Bertocci
,
I
,
Dell’Anno
,
A
,
Musco
,
L
,
Gambi
,
C
,
Saggiomo
,
V
,
Cannavacciuolo
,
M
,
Lo Martire
,
M
,
Passarelli
,
A
,
Zazo
,
G
,
Danovaro
,
R.
2019
.
Multiple human pressures in coastal habitats: Variation of meiofaunal assemblages associated with sewage discharge in a post-industrial area
.
Science of The Total Environment
655
:
1218
1231
. DOI: http://dx.doi.org/10.1016/j.scitotenv.2018.11.121.
Binkley
,
A.
2020
.
Where is the federal port modernization review headed?
Canadian Sailings
:
12
16
.
Available at
https://www.canadiansailings.ca/documents/weekly_issue/feb_2020/feb_24_2020.pdf.
Accessed January 23, 2023
.
Bishop
,
MJ
,
Mayer-Pinto
,
M
,
Airoldi
,
L
,
Firth
,
LB
,
Morris
,
RL
,
Loke
,
LHL
,
Hawkins
,
SJ
,
Naylor
,
LA
,
Coleman
,
RA
,
Chee
,
SY
,
Dafforn
,
KA.
2017
.
Effects of ocean sprawl on ecological connectivity: Impacts and solutions
.
Journal of Experimental Marine Biology and Ecology
492
:
7
30
. DOI: http://dx.doi.org/10.1016/j.jembe.2017.01.021.
Borja
,
Á
,
Dauer
,
DM
,
Grémare
,
A.
2012
.
The importance of setting targets and reference conditions in assessing marine ecosystem quality
.
Ecological Indicators
12
(
1
):
1
7
. DOI: http://dx.doi.org/10.1016/j.ecolind.2011.06.018.
Borja
,
Á
,
Franco
,
J
,
Pérez
,
V.
2000
.
A Marine Biotic Index to establish the ecological quality of soft-bottom benthos within European estuarine and coastal environments
.
Marine Pollution Bulletin
40
(
12
):
1100
1114
. DOI: http://dx.doi.org/10.1016/S0025-326X(00)00061-8.
Brown
,
CJ
,
Saunders
,
MI
,
Possingham
,
HP
,
Richardson
,
AJ.
2014
.
Interactions between global and local stressors of ecosystems determine management effectiveness in cumulative impact mapping
.
Diversity and Distributions
20
(
5
):
538
546
. DOI: http://dx.doi.org/10.1111/ddi.12159.
Bulleri
,
F
,
Chapman
,
MG.
2010
.
The introduction of coastal infrastructure as a driver of change in marine environments
.
Journal of Applied Ecology
47
(
1
):
26
35
. DOI: http://dx.doi.org/10.1111/j.1365-2664.2009.01751.x.
Byrnes
,
TA
,
Dunn
,
RJK.
2020
.
Boating- and shipping-related environmental impacts and example management measures: A review
.
Journal of Marine Science and Engineering
8
(
11
):
1
49
. DOI: http://dx.doi.org/10.3390/jmse8110908.
Callier
,
MD
,
Richard
,
M
,
McKindsey
,
CW
,
Archambault
,
P
,
Desrosiers
,
G.
2009
.
Responses of benthic macrofauna and biogeochemical fluxes to various levels of mussel biodeposition: An in situ “benthocosm” experiment
.
Marine Pollution Bulletin
58
(
10
):
1544
1553
. DOI: http://dx.doi.org/10.1016/j.marpolbul.2009.05.010.
Carrier-Belleau
,
C
,
Drolet
,
D
,
McKindsey
,
CW
,
Archambault
,
P.
2021
.
Environmental stressors, complex interactions and marine benthic communities’ responses
.
Scientific Reports
11
(
1
):
1
15
. DOI: http://dx.doi.org/10.1038/s41598-021-83533-1.
Carrier-Belleau
,
C
,
Pascal
,
L
,
Tiegs
,
SD
,
Nozais
,
C
,
Archambault
,
P.
2022
.
From organismal physiology to ecological processes: Effects of nutrient enrichment and salinity variation in a freshwater ecosystem
.
Limnology and Oceanography
68
:
S115
S130
. DOI: http://dx.doi.org/10.1002/lno.12269.
Carrière
,
J.
2018
.
Observatoire environnemental de la Baie de Sept Îles.
Sept-Îles, Quebec
:
Institut Nordique de Recherche en Environnement et en Santé au Travail
.
Christensen
,
PB
,
Glud
,
RN
,
Dalsgaard
,
T
,
Gillespie
,
P.
2003
.
Impacts of longline mussel farming on oxygen and nitrogen dynamics and biological communities of coastal sediments
.
Aquaculture
218
(
1–4
):
567
588
. DOI: http://dx.doi.org/10.1016/S0044-8486(02)00587-2.
Clarke
,
KR
,
Gorley
,
RN.
2006
. PRIMER v6: User manual/tutorial.
Plymouth, UK
:
PRIMER-E Ltd
.
Cotano
,
U
,
Villate
,
F.
2006
.
Anthropogenic influence on the organic fraction of sediments in two contrasting estuaries: A biochemical approach
.
Marine Pollution Bulletin
52
(
4
):
404
414
. DOI: http://dx.doi.org/10.1016/j.marpolbul.2005.09.027.
Côté
,
IM
,
Darling
,
ES
,
Brown
,
CJ.
2016
.
Interactions among ecosystem stressors and their importance in conservation
.
Proceedings of the Royal Society B: Biological Sciences
283
(
1824
):
1
9
. DOI: http://dx.doi.org/10.1098/rspb.2015.2592.
Crain
,
CM
,
Kroeker
,
K
,
Halpern
,
BS.
2008
.
Interactive and cumulative effects of multiple human stressors in marine systems
.
Ecology Letters
11
:
1304
1315
. DOI: http://dx.doi.org/10.1111/j.1461-0248.2008.01253.x.
Cranford
,
P
,
Dowd
,
M
,
Grant
,
J
,
Hargrave
,
B
,
McGladdery
,
S.
2003
. Ecosystem level effects of marine bivalve aquaculture, in
Hargrave
,
BT
,
Cranford
,
P
,
Dowd
,
M
,
Grant
,
B
,
McGladdery
,
S
,
Burridge
,
LE
eds.,
A scientific review of the potential environmental effects of aquaculture in aquatic ecosystems. Canadian Technical Report of Fisheries and Aquatic Sciences 2450
.
Ottawa, Canada
:
Fisheries and Oceans Canada
:
51
95
.
Crawford
,
CM
,
Macleod
,
CKA
,
Mitchell
,
IM.
2003
.
Effects of shellfish farming on the benthic environment
.
Aquaculture
224
(
1–4
):
117
140
. DOI: http://dx.doi.org/10.1016/S0044-8486(03)00210-2.
Culhane
,
FE
,
Briers
,
RA
,
Tett
,
P
,
Fernandes
,
TF.
2019
.
Response of a marine benthic invertebrate community and biotic indices to organic enrichment from sewage disposal
.
Journal of the Marine Biological Association of the United Kingdom
99
(
8
):
1721
1734
. DOI: http://dx.doi.org/10.1017/S0025315419000857.
D’Alessandro
,
M
,
Porporato
,
EMD
,
Esposito
,
V
,
Giacobbe
,
S
,
Deidun
,
A
,
Nasi
,
F
,
Ferrante
,
L
,
Auriemma
,
R
,
Berto
,
D
,
Renzi
,
M
,
Scotti
,
G
,
Consoli
,
P
,
Del Negro
,
P
,
Andaloro
,
F
,
Romeo
,
T.
2020
.
Common patterns of functional and biotic indices in response to multiple stressors in marine harbours ecosystems
.
Environmental Pollution
259
:
113959
. DOI: http://dx.doi.org/10.1016/j.envpol.2020.113959.
Daigle
,
RM
,
Archambault
,
P
,
Halpern
,
BS
,
Lowndes
,
JSS
,
Côté
,
IM.
2017
.
Incorporating public priorities in the Ocean Health Index: Canada as a case study
.
PLoS One
12
(
5
):
e0178044
. DOI: http://dx.doi.org/10.1371/journal.pone.0178044.
Darling
,
ES
,
Côté
,
IM.
2008
.
Quantifying the evidence for ecological synergies
.
Ecology Letters
11
(
12
):
1278
1286
. DOI: http://dx.doi.org/10.1111/j.1461-0248.2008.01243.x.
Dauer
,
DM.
1993
.
Biological criteria, environmental health and estuarine macrobenthic community structure
.
Marine Pollution Bulletin
26
(
5
):
249
257
. DOI: http://dx.doi.org/10.1016/0025-326X(93)90063-P.
Dauvin
,
JC
,
Andrade
,
H
,
De-La-Ossa-Carretero
,
JA
,
Del-Pilar-Ruso
,
Y
,
Riera
,
R.
2016
.
Polychaete/amphipod ratios: An approach to validating simple benthic indicators
.
Ecological Indicators
63
:
89
99
. DOI: http://dx.doi.org/10.1016/j.ecolind.2015.11.055.
Dauvin
,
JC
,
Ruellet
,
T.
2007
.
Polychaete/amphipod ratio revisited
.
Marine Pollution Bulletin
55
(
1–6
):
215
224
. DOI: http://dx.doi.org/10.1016/j.marpolbul.2006.08.045.
Davis
,
SJ
,
Ó hUallacháin
,
D
,
Mellander
,
PE
,
Kelly
,
AM
,
Matthaei
,
CD
,
Piggott
,
JJ
,
Kelly-Quinn
,
M.
2018
.
Multiple-stressor effects of sediment, phosphorus and nitrogen on stream macroinvertebrate communities
.
Science of The Total Environment
637–638
:
577
587
. DOI: https://doi.org/10.1016/j.scitotenv.2018.05.052.
Demers
,
KA
,
Le Hénaff
,
A
,
Carrière
,
J.
2018
. État des glaces, in
Carrière
,
J
ed.,
Observatoire environnemental de la Baie de Sept-Îles
.
Sept-Îles, Quebec
:
INREST
:
593
612
.
Department of Fisheries and Oceans
.
2009
.
Development of a framework and principles for the biogeographic classification of Canadian marine areas
.
Canadian Science Advisory Secretariat Advisory Report 2009/056:17
.
Ottawa, Canada
:
Fisheries and Oceans Canada
.
Available at
https://waves-vagues.dfompo.gc.ca/library-bibliotheque/338958.pdf.
Department of Fisheries and Oceans
.
2019
. Update of stock status indicators for Northern Shrimp in the Estuary and Gulf of St. Lawrence.
Canadian Science Advisory Secretariat Science Response 2019/005:7
.
Ottawa, Canada
:
Fisheries and Oceans Canada
.
Desprez
,
M.
2000
.
Physical and biological impact of marine aggregate extraction along the French coast of the Eastern English Channel: Short-and long-term post-dredging restoration
.
ICES Journal of Marine Science
57
(
5
):
1428
1438
. DOI: http://dx.doi.org/10.1006/jmsc.2000.0926.
Díaz
,
S
,
Quétier
,
F
,
Cáceres
,
DM
,
Trainor
,
SF
,
Pérez-Harguindeguy
,
N
,
Bret-Harte
,
MS
,
Finegan
,
B
,
Peña-Claros
,
M
,
Poorter
,
L.
2011
.
Linking functional diversity and social actor strategies in a framework for interdisciplinary analysis of nature’s benefits to society
.
Proceedings of the National Academy of Sciences
108
(
3
):
895
902
. DOI: http://dx.doi.org/10.1073/pnas.1017993108.
Dijkstra
,
EW.
1959
.
A note on two problems in connexion with graphs
.
Numerische Mathematik
1
(
1
):
269
271
. DOI: http://dx.doi.org/10.1007/BF01386390.
Dreujou
,
E
,
Carrier-Belleau
,
C
,
Goldsmit
,
J
,
Fiorentino
,
D
,
Ben-Hamadou
,
R
,
Muelbert
,
JH
,
Godbold
,
JA
,
Daigle
,
RM
,
Beauchesne
,
D.
2020
a.
Holistic environmental approaches and Aichi biodiversity targets: Accomplishments and perspectives for marine ecosystems
.
PeerJ
8
:
e8171
. DOI: http://dx.doi.org/10.7717/peerj.8171.
Dreujou
,
E
,
Desroy
,
N
,
Carrière
,
J
,
Tréau de Coeli
,
L
,
McKindsey
,
CW
,
Archambault
,
P.
2021
.
Determining the ecological status of benthic coastal communities: A case in an anthropized sub-Arctic area
.
Frontiers in Marine Science
8
:
637546
. DOI: http://dx.doi.org/10.3389/fmars.2021.637546.
Dreujou
,
E
,
McKindsey
,
CW
,
Grant
,
C
,
Tréau de Coeli
,
L
,
St-Louis
,
R
,
Archambault
,
P.
2020
b.
Biodiversity and habitat assessment of coastal benthic communities in a sub-Arctic industrial harbor area
.
Water
12
(
9
):
2424
. DOI: http://dx.doi.org/10.3390/w12092424.
Dreujou
,
E
,
Paquette
,
L
,
Grant
,
C
,
Archambault
,
P
,
Carrière
,
J.
2018
. Caractérisation de la faune benthique, in
Carrière
,
J
ed.,
Observatoire environnemental de la Baie de Sept-Îles
.
Sept-Îles, Quebec
:
INREST
:
377
452
.
Drouin
,
A
,
Archambault
,
P
,
Sirois
,
P.
2011
.
Distinction of nektonic and benthic communities between fish-present (Salvelinus fontinalis) and natural fishless lakes
.
Boreal Environmental Research
16
:
101
114
.
Dufrêne
,
M
,
Legendre
,
P.
1997
.
Species assemblages and indicator species: The need for a flexible asymmetrical approach
.
Ecological Monographs
67
(
3
):
345
366
. DOI: http://dx.doi.org/10.1890/0012-9615(1997)067[0345:SAAIST]2.0.CO;2.
Dutil
,
J-D
,
Proulx
,
S
,
Galbraith
,
PS
,
Chasse
,
J
,
Lambert
,
N
,
Laurian
,
C.
2012
.
Coastal and epipelagic habitats of the estuary and Gulf of St. Lawrence
.
Canadian Technical Report of Fisheries and Aquatic Sciences 3009
.
Mont-Joli, Canada
:
Fisheries and Oceans Canada
:
1
87
.
Ellis
,
JI
,
Clark
,
D
,
Atalah
,
J
,
Jiang
,
W
,
Taiapa
,
C
,
Patterson
,
M
,
Sinner
,
J
,
Hewitt
,
J.
2017
.
Multiple stressor effects on marine infauna: Responses of estuarine taxa and functional traits to sedimentation, nutrient and metal loading
.
Scientific Reports
7
(
1
):
12013
. DOI: http://dx.doi.org/10.1038/s41598-017-12323-5.
European Environmental Agency
.
1999
.
Environmental indicators: Typology and overview
.
Copenhagen, Denmark
:
European Union
.
Report No. 25
.
Available at
https://www.eea.europa.eu/publications/TEC25.
Accessed January 23, 2023
.
Feist
,
BE
,
Levin
,
PS.
2016
.
Novel indicators of anthropogenic influence on marine and coastal ecosystems
.
Frontiers in Marine Science
3
:
1
13
. DOI: http://dx.doi.org/10.3389/fmars.2016.00113.
Ferrario
,
F
,
Araújo
,
CAS
,
Bélanger
,
S
,
Bourgault
,
D
,
Carrière
,
J
,
Carrier-Belleau
,
C
,
Dreujou
,
E
,
Johnson
,
LE
,
Juniper
,
SK
,
Mabit
,
R
,
McKindsey
,
CW
,
Ogston
,
L
,
Picard
,
MMM
,
Saint-Louis
,
R
,
Saulnier-Talbot
,
E
,
Shaw
,
JL
,
Templeman
,
N
,
Therriault
,
TW
,
Tremblay
,
JE
,
Archambault
,
P.
2022
.
Holistic environmental monitoring in ports as an opportunity to advance sustainable development, marine science, and social inclusiveness
.
Elementa: Science of the Anthropocene
10
(
1
):
1
21
. DOI: http://dx.doi.org/10.1525/elementa.2021.00061.
Ferrario
,
F
,
Archambault
,
P
,
Templeman
,
N.
2021
. A scan of environmental monitoring in top ports around the globe.
Canadian Technical Report of Fisheries and Aquatic Sciences 3428
.
Ottawa, Canada
:
Fisheries and Oceans Canada
:
vii + 36
.
Galic
,
N
,
Sullivan
,
LL
,
Grimm
,
V
,
Forbes
,
VE.
2018
.
When things don’t add up: Quantifying impacts of multiple stressors from individual metabolism to ecosystem processing
.
Ecology Letters
21
(
4
):
568
577
. DOI: http://dx.doi.org/10.1111/ele.12923.
Gallardi
,
D.
2014
.
Effects of bivalve aquaculture on the environment and their possible mitigation: A review
.
Fisheries and Aquaculture Journal
05
(
03
). DOI: http://dx.doi.org/10.4172/2150-3508.1000105.
Gari
,
SR
,
Newton
,
A
,
Icely
,
JD.
2015
.
A review of the application and evolution of the DPSIR framework with an emphasis on coastal social-ecological systems
.
Ocean & Coastal Management
103
:
63
77
. DOI: http://dx.doi.org/10.1016/j.ocecoaman.2014.11.013.
Glaser
,
M
,
Krause
,
G
,
Ratter
,
BMW
,
Welp
,
M.
2012
.
Human-nature interactions in the Anthropocene
. 1st ed.
New York, NY
:
Routledge
.
Grall
,
J
,
Glémarec
,
M.
1997
.
Using biotic indices to estimate macrobenthic community perturbations in the Bay of Brest
.
Estuarine, Coastal and Shelf Science
44
:
43
53
. DOI: http://dx.doi.org/10.1016/S0272-7714(97)80006-6.
Gray
,
JS
,
Elliott
,
M.
2009
.
Ecology of marine sediments: From science to management
.
Oxford, UK
:
Oxford University Press
.
Available at
https://academic.oup.com/book/41808.
Accessed January 23, 2023
.
Halpern
,
BS
,
Frazier
,
M
,
Afflerbach
,
J
,
Lowndes
,
JS
,
Micheli
,
F
,
O’Hara
,
C
,
Scarborough
,
C
,
Selkoe
,
KA.
2019
.
Recent pace of change in human impact on the world’s ocean
.
Scientific Reports
9
(
1
):
11609
. DOI: http://dx.doi.org/10.1038/s41598-019-47201-9.
Halpern
,
BS
,
Fujita
,
R.
2013
.
Assumptions, challenges, and future directions in cumulative impact analysis
.
Ecosphere
4
(
10
):
1
11
. DOI: http://dx.doi.org/10.1890/ES13-00181.1.
Halpern
,
BS
,
Selkoe
,
KA
,
Micheli
,
F
,
Kappel
,
CV.
2007
.
Evaluating and ranking the vulnerability of global marine ecosystems to anthropogenic threats
.
Conservation Biology
21
(
5
):
1301
1315
. DOI: http://dx.doi.org/10.1111/j.1523-1739.2007.00752.x.
Halpern
,
BS
,
Walbridge
,
S
,
Selkoe
,
KA
,
Kappel
,
CV
,
Micheli
,
F
,
D’Agrosa
,
C
,
Bruno
,
JF
,
Casey
,
KS
,
Ebert
,
C
,
Fox
,
HE
,
Fujita
,
R
,
Heinemann
,
D
,
Lenihan
,
HS
,
Madin
,
EMP
,
Perry
,
MT
,
Selig
,
ER
,
Spalding
,
M
,
Steneck
,
R
,
Watson
,
R.
2008
.
A global map of human impact on marine ecosystems
.
Science
319
(
5865
):
948
952
. DOI: http://dx.doi.org/10.1126/science.1149345.
Heery
,
EC
,
Bishop
,
MJ
,
Critchley
,
LP
,
Bugnot
,
AB
,
Airoldi
,
L
,
Mayer-Pinto
,
M
,
Sheehan
,
EV
,
Coleman
,
RA
,
Loke
,
LHL
,
Johnston
,
EL
,
Komyakova
,
V
,
Morris
,
RL
,
Strain
,
EMA
,
Naylor
,
LA
,
Dafforn
,
KA.
2017
.
Identifying the consequences of ocean sprawl for sedimentary habitats
.
Journal of Experimental Marine Biology and Ecology
492
:
31
48
. DOI: http://dx.doi.org/10.1016/j.jembe.2017.01.020.
Hijmans
,
RJ.
2020
.
Raster: Geographic data analysis and modeling
.
Available at
https://CRAN.R-project.org/package=raster.
Accessed January 23, 2023
.
Hodgson
,
EE
,
Halpern
,
BS
,
Essington
,
TE.
2019
.
Moving beyond silos in cumulative effects assessment
.
Frontiers in Ecology and Evolution
7
:
211
. DOI: http://dx.doi.org/10.3389/fevo.2019.00211.
International Council for the Exploration of the Sea
.
2001
.
ICES Annual Report
.
International Council for the Exploration of the Sea
. DOI: http://dx.doi.org/10.17895/ices.pub.7436.
Islam
,
MS
,
Tanaka
,
M.
2004
.
Impacts of pollution on coastal and marine ecosystems including coastal and marine fisheries and approach for management: A review and synthesis
.
Marine Pollution Bulletin
48
(
7–8
):
624
649
. DOI: http://dx.doi.org/10.1016/j.marpolbul.2003.12.004.
Jägerbrand
,
AK
,
Brutemark
,
A
,
Barthel Svedén
,
J
,
Gren
,
IM.
2019
.
A review on the environmental impacts of shipping on aquatic and nearshore ecosystems
.
Science of The Total Environment
695
:
133637
. DOI: http://dx.doi.org/10.1016/j.scitotenv.2019.133637.
Judd
,
AD
,
Backhaus
,
T
,
Goodsir
,
F.
2015
.
An effective set of principles for practical implementation of marine cumulative effects assessment
.
Environmental Science and Policy
54
:
254
262
. DOI: http://dx.doi.org/10.1016/j.envsci.2015.07.008.
Korpinen
,
S
,
Andersen
,
JH.
2016
.
A global review of cumulative pressure and impact assessments in marine environments
.
Frontiers in Marine Science
3
:
153
. DOI: http://dx.doi.org/10.3389/fmars.2016.00153.
Lacoste
,
É
,
Weise
,
AM
,
Lavoie
,
M-F
,
Archambault
,
P
,
McKindsey
,
CW.
2019
.
Changes in infaunal assemblage structure influence nutrient fluxes in sediment enriched by mussel biodeposition
.
Science of The Total Environment
692
:
39
48
. DOI: http://dx.doi.org/10.1016/j.scitotenv.2019.07.235.
Legendre
,
P
,
Legendre
,
LFJ.
1998
.
Numerical ecology
. 2nd ed.
Amsterdam, the Netherlands
:
Elsevier
.
Levin
,
PS
,
Fogarty
,
MJ
,
Murawski
,
SA
,
Fluharty
,
D.
2009
.
Integrated ecosystem assessments: Developing the scientific basis for ecosystem-based management of the ocean
.
PLoS Biology
7
(
1
):
e1000014
. DOI: http://dx.doi.org/10.1371/journal.pbio.1000014.
Levinton
,
JS.
2013
.
Marine biology: Function, biodiversity, ecology
. 4th ed.
Oxford, UK
:
Oxford University Press
.
Link
,
JS.
2002
.
What does ecosystem-based fisheries management mean?
Fisheries
27
(
4
):
18
21
.
Magurran
,
AE
,
McGill
,
BJ.
2011
.
Biological diversity: Frontiers in measurement and assessment
. 1st ed.
Oxford, UK
:
Oxford University Press
.
Margules
,
CR
,
Pressey
,
RL.
2000
.
Systematic conservation planning
.
Nature
405
(
6783
):
243
253
. DOI: http://dx.doi.org/10.1038/35012251.
McArdle
,
BH
,
Anderson
,
MJ.
2001
.
Fitting multivariate models to community data: A comment on distance-based redundancy analysis
.
Ecology
82
(
1
):
290
297
. DOI: http://dx.doi.org/10.1890/0012-9658(2001)082[0290:FMMTCD]2.0.CO;2.
McKindsey
,
CW
,
Archambault
,
P
,
Callier
,
MD
,
Olivier
,
F.
2011
.
Influence of suspended and off-bottom mussel culture on the sea bottom and benthic habitats: A review
.
Canadian Journal of Zoology
89
(
7
):
622
646
. DOI: http://dx.doi.org/10.1139/z11-037.
Micheli
,
F
,
Heiman
,
KW
,
Kappel
,
CV
,
Martone
,
RG
,
Sethi
,
SA
,
Osio
,
GC
,
Fraschetti
,
S
,
Shelton
,
AO
,
Tanner
,
JM.
2016
.
Combined impacts of natural and human disturbances on rocky shore communities
.
Ocean and Coastal Management
126
:
42
50
. DOI: http://dx.doi.org/10.1016/j.ocecoaman.2016.03.014.
Millenniun Environmental Assessment
.
2005
.
Ecosystems and human well-being: Biodiversity synthesis
.
Washington, DC
:
Island Press
.
Mineur
,
F
,
Cook
,
E
,
Minchin
,
D
,
Bohn
,
K
,
Macleod
,
A
,
Maggs
,
C.
2012
.
Changing coasts: Marine aliens and artificial structures
.
Oceanography and Marine Biology: An Annual Review
50
:
198
243
. DOI: http://dx.doi.org/10.1201/b12157-6.
Momota
,
K
,
Hosokawa
,
S.
2021
.
Potential impacts of marine urbanization on benthic macrofaunal diversity
.
Scientific Reports
11
(
1
):
1
12
. DOI: http://dx.doi.org/10.1038/s41598-021-83597-z.
Müller
,
A
,
Österlund
,
H
,
Marsalek
,
J
,
Viklander
,
M.
2020
.
The pollution conveyed by urban runoff: A review of sources
.
Science of The Total Environment
709
:
136125
. DOI: http://dx.doi.org/10.1016/j.scitotenv.2019.136125.
Muxika
,
I
,
Borja
,
Á
,
Bald
,
J.
2007
.
Using historical data, expert judgement and multivariate analysis in assessing reference conditions and benthic ecological status, according to the European Water Framework Directive
.
Marine Pollution Bulletin
55
(
1–6
):
16
29
. DOI: http://dx.doi.org/10.1016/j.marpolbul.2006.05.025.
Ocaña
,
FA
,
Pech
,
D
,
Simões
,
N
,
Hernández-Ávila
,
I.
2019
.
Spatial assessment of the vulnerability of benthic communities to multiple stressors in the Yucatan Continental Shelf, Gulf of Mexico
.
Ocean and Coastal Management
181
:
104900
. DOI: http://dx.doi.org/10.1016/j.ocecoaman.2019.104900.
Oesterwind
,
D
,
Rau
,
A
,
Zaiko
,
A.
2016
.
Drivers and pressures—Untangling the terms commonly used in marine science and policy
.
Journal of Environmental Management
181
:
8
15
. DOI: http://dx.doi.org/10.1016/j.jenvman.2016.05.058.
Okey
,
TA
,
Alidina
,
HM
,
Agbayani
,
S.
2015
.
Mapping ecological vulnerability to recent climate change in Canada’s Pacific marine ecosystems
.
Ocean and Coastal Management
106
:
35
48
. DOI: http://dx.doi.org/10.1016/j.ocecoaman.2015.01.009.
Oksanen
,
J
,
Simpson
,
GL
,
Blanchet
,
FG
,
Kindt
,
R
,
Legendre
,
P
,
Minchin
,
PR
,
O’Hara
,
RB
,
Solymos
,
P
,
Stevens
,
MHH
,
Szoecs
,
E
,
Wagner
,
H
,
Barbour
,
M
,
Bedward
,
M
,
Bolker
,
B
,
Bocard
,
D
,
Carvahlo
,
G
,
Chirico
,
M
,
De Caceres
,
M
,
Durand
,
S
,
Evangelista
,
HBA
,
FitzJohn
,
R
,
Friendly
,
M
,
Furnaeaux
,
B
,
Hannigan
,
G
,
Hill
,
MO
,
Lahti
,
L
,
McGlinn
,
D
,
Ouelette
,
MH
,
Cunha
,
ER
,
Smith
,
T
,
Stier
,
A
,
Ter Braak
,
CJF
,
Weedon
,
J.
2022
.
Vegan: Community ecology package
.
Available at
https://CRAN.R-project.org/package=vegan.
Accessed January 23, 2023
.
Oviatt
,
C
,
Quinn
,
J
,
Maughan
,
J
,
Ellis
,
J
,
Sullivan
,
B
,
Gearing
,
J
,
Gearing
,
P
,
Hunt
,
C
,
Sampou
,
P
,
Latimer
,
J.
1987
.
Fate and effects of sewage sludge in the coastal marine environment: A mesocosm experiment
.
Marine Ecology Progress Series
41
:
187
203
. DOI: http://dx.doi.org/10.3354/meps041187.
Parravicini
,
V
,
Rovere
,
A
,
Vassallo
,
P
,
Micheli
,
F
,
Montefalcone
,
M
,
Morri
,
C
,
Paoli
,
C
,
Albertelli
,
G
,
Fabiano
,
M
,
Bianchi
,
CN.
2012
.
Understanding relationships between conflicting human uses and coastal ecosystems status: A geospatial modeling approach
.
Ecological Indicators
19
:
253
263
. DOI: http://dx.doi.org/10.1016/j.ecolind.2011.07.027.
Pearson
,
TH
,
Rosenberg
,
R.
1978
.
Macrobenthic succession in relation to organic enrichment and pollution of the marine environment
.
Oceanography and Marine Biology: An Annual Review
16
:
229
311
.
Pebesma
,
E.
2018
.
Simple features for R: Standardized support for spatial cector data
.
The R Journal
10
(
1
):
439
. DOI: http://dx.doi.org/10.32614/RJ-2018-009.
Piacenza
,
SE
,
Barner
,
AK
,
Benkwitt
,
CE
,
Boersma
,
KS
,
Cerny-Chipman
,
EB
,
Ingeman
,
KE
,
Kindinger
,
TL
,
Lee
,
JD
,
Lindsley
,
AJ
,
Reimer
,
JN
,
Rowe
,
JC
,
Shen
,
C
,
Thompson
,
KA
,
Thurman
,
LL
,
Heppell
,
SS.
2015
.
Patterns and variation in benthic biodiversity in a large marine ecosystem
.
PLoS One
10
(
8
):
1
23
. DOI: http://dx.doi.org/10.1371/journal.pone.0135135.
Piggott
,
JJ
,
Townsend
,
CR
,
Matthaei
,
CD.
2015
.
Reconceptualizing synergism and antagonism among multiple stressors
.
Ecology and Evolution
5
(
7
):
1538
1547
. DOI: http://dx.doi.org/10.1002/ece3.1465.
Pikitch
,
EK
,
Santora
,
C
,
Babcock
,
EA
,
Bakun
,
A
,
Bonfil
,
R
,
Conover
,
DO
,
Dayton
,
P
,
Doukakis
,
P
,
Fluharty
,
D
,
Heneman
,
B
,
Houde
,
ED
,
Link
,
J
,
Livingston
,
PA
,
Mangel
,
M
,
McAllister
,
MK
,
Pope
,
J
,
Sainsbury
,
KJ.
2004
.
Ecosystem-based fishery management
.
Science
305
(
5682
):
346
347
. DOI: http://dx.doi.org/10.1126/science.1098222.
Pinto
,
R
,
Patrício
,
J
,
Baeta
,
A
,
Fath
,
BD
,
Neto
,
JM
,
Marques
,
JC.
2009
.
Review and evaluation of estuarine biotic indices to assess benthic condition
.
Ecological Indicators
9
(
1
):
1
25
. DOI: http://dx.doi.org/10.1016/j.ecolind.2008.01.005.
Port de Sept-Îles
.
2021
.
Statistiques relatives aux activités portuaires
.
Available at
http://www.portsi.com/wp-content/uploads/2021/01/2020-12-PSI-Stats.pdf.
Accessed January 23, 2023
.
Quinn
,
GP
,
Keough
,
MJ.
2002
.
Experimental design and data analysis for biologists
.
Cambridge, UK
:
Cambridge University Press
.
R Core Team
.
2022
.
R: A language and environment for statistical computing
.
Vianna, Austria
:
R Foundation for Statistical Computing
.
Redman
,
CL
,
Grove
,
JM
,
Kuby
,
LH.
2004
.
Integrating social science into the Long-Term Ecological Research (LTER) network: Social dimensions of ecological change and ecological dimensions of social change
.
Ecosystems
7
(
2
):
161
171
. DOI: http://dx.doi.org/10.1007/s10021-003-0215-z.
Richard
,
M
,
Archambault
,
P
,
Thouzeau
,
G
,
Desrosiers
,
G.
2007
.
Summer influence of 1 and 2 yr old mussel cultures on benthic fluxes in Grande-Entrée lagoon, Îles-de-la-Madeleine (Québec, Canada)
.
Marine Ecology Progress Series
338
:
131
143
. DOI: http://dx.doi.org/10.3354/meps338131.
Santos
,
C
,
Ehler
,
CN
,
Agardy
,
T
,
Andrade
,
F
,
Orbach
,
MK
,
Crowder
,
LB.
2019
.
Marine spatial planning
, in
Sheppard
,
C
ed.,
World seas: An environmental evaluation
. 2nd ed.
Cambridge, MA
:
Academic Press
:
571
592
. DOI: http://dx.doi.org/10.1016/B978-0-12-805052-1.00033-4.
Schloss
,
I
,
Archambault
,
P
,
Beauchesne
,
D
,
Bourgault
,
D
,
Cusson
,
M
,
Dumont
,
D
,
Ferreyra
,
G
,
Levasseur
,
M
,
Pelletier
,
É
,
St-Louis
,
R
,
Tremblay
,
R.
2017
. Impacts potentiels cumulés des facteurs de stress liés aux activités humaines sur l’écosystème marin du Saint-Laurent dans Les hydrocarbures dans le golfe du Saint-Laurent—Enjeux sociaux, économiques et environnementaux, in
Archambault
,
P
,
Schloss
,
I
,
Grant
,
C
,
Plante
,
S
eds.,
Les hydrocarbures dans le Golfe du Saint-Laurent—Enjeux sociaux, économiques et environnementaux
.
Rimouski, Quebec
:
Notre Golfe
:
132
165
.
Séguin
,
A
,
Gravel
,
D
,
Archambault
,
P.
2014
.
Effect of disturbance regime on alpha and beta diversity of rock pools
.
Diversity
6
(
1
):
1
17
. DOI: http://dx.doi.org/10.3390/d6010001.
Shaw
,
J-L.
2019
.
Hydrodynamique de la Baie de Sept-Îles [M.Sc. thesis]
.
Rimouski, Quebec
:
Université du Québec à Rimouski
.
Simboura
,
N
,
Zenetos
,
A.
2002
.
Benthic indicators to use in ecological quality classification of Mediterranean soft bottom marine ecosystems, including a new biotic index
.
Mediterranean Marine Science
3
(
2
):
77
111
. DOI: https://dx.doi.org/10.12681/mms.249.
Snelgrove
,
PVR
,
Archambault
,
P
,
Kim Juniper
,
S
,
Lawton
,
P
,
Metaxas
,
A
,
Pepin
,
P
,
Rice
,
JC
,
Tunnicliffe
,
V.
2012
.
Canadian Healthy Oceans Network (CHONe): An academic–government partnership to develop scientific guidelines for conservation and sustainable usage of marine biodiversity
.
Fisheries
37
(
7
):
296
304
. DOI: http://dx.doi.org/10.1080/03632415.2012.696002.
Socioeconomic Data and Applications Center
.
2020
.
Percentage of total population living in coastal areas
.
Available at
https://www.un.org/esa/sustdev/natlinfo/indicators/methodology_sheets/oceans_seas_coasts/pop_coastal_areas.pdf.
Accessed January 23, 2023
.
Solan
,
M
,
Whiteley
,
N.
2016
.
Stressors in the marine environment
.
Oxford, UK
:
Oxford University Press
.
Statistics Canada
.
2011
.
Shipping in Canada: 2011
.
Ottawa, Canada
:
Statistics Canada
.
Available at
https://www150.statcan.gc.ca/n1/en/pub/54-205-x/54-205-x2011000-eng.pdf?st=rUehQMzc.
Accessed January 23, 2023
.
Stelzenmüller
,
V
,
Lee
,
J
,
South
,
A
,
Rogers
,
SI.
2010
.
Quantifying cumulative impacts of human pressures on the marine environment: A geospatial modelling framework
.
Marine Ecology Progress Series
398
:
19
32
. DOI: http://dx.doi.org/10.3354/meps08345.
Teixeira
,
H
,
Berg
,
T
,
Uusitalo
,
L
,
Fürhaupter
,
K
,
Heiskanen
,
AS
,
Mazik
,
K
,
Lynam
,
CP
,
Neville
,
S
,
Rodriguez
,
JG
,
Papadopoulou
,
N
,
Moncheva
,
S
,
Churilova
,
T
,
Kryvenko
,
O
,
Krause-Jensen
,
D
,
Zaiko
,
A
,
Verissimo
,
H
,
Pantazi
,
M
,
Carvalho
,
S
,
Patrício
,
J
,
Uyarra
,
MC
,
Borja
À.
2016
.
A catalogue of marine biodiversity indicators
.
Frontiers in Marine Science
3
:
207
. DOI: http://dx.doi.org/10.3389/fmars.2016.00207.
Tičina
,
V
,
Katavić
,
I
,
Grubišić
,
L.
2020
.
Marine aquaculture impacts on marine biota in oligotrophic environments of the Mediterranean Sea—A review
.
Frontiers in Marine Science
7
:
1
11
. DOI: http://dx.doi.org/10.3389/fmars.2020.00217.
van de Bund
,
W
,
Solimini
,
AG.
2007
.
Ecological quality ratios for ecological quality assessment in inland and marine waters
.
Luxembourg
:
European Commission Joint Research Centre, Institute for Environment; Sustainability
.
Available at
http://scholar.google.com/scholar?hl=en&btnG=Search&q=intitle:Ecological+Quality+Ratios+for+ecological+quality+assessment+in+inland+and+marine+waters#8.
Accessed January 23, 2023
.
van Etten
,
J.
2017
.
R package gdistance: Distances and routes on geographical grids
.
Journal of Statistical Software
76
(
13
). DOI: http://dx.doi.org/10.18637/jss.v076.i13.
Wei
,
CL
,
Cusson
,
M
,
Archambault
,
P
,
Belley
,
R
,
Brown
,
T
,
Burd
,
BJ
,
Edinger
,
E
,
Kenchington
,
E
,
Gilkinson
,
K
,
Lawton
,
P
,
Link
,
H
,
Ramey-Balci
,
PA
,
Scrosati
,
RA
,
Snelgrove
,
PVR.
2020
.
Seafloor biodiversity of Canada’s three oceans: Patterns, hotspots and potential drivers
.
Diversity and Distributions
26
(
2
):
226
241
. DOI: http://dx.doi.org/10.1111/ddi.13013.
Wilding
,
TA
,
Nickell
,
TD.
2013
.
Changes in benthos associated with mussel (Mytilus edulis L.) farms on the west-coast of Scotland
.
PLoS One
8
(
7
):
e68313
. DOI: http://dx.doi.org/10.1371/journal.pone.0068313.
Wilson
,
K
,
Pressey
,
RL
,
Newton
,
A
,
Burgman
,
M
,
Possingham
,
H
,
Weston
,
C.
2005
.
Measuring and incorporating vulnerability into conservation planning
.
Environmental Management
35
(
5
):
527
543
. DOI: http://dx.doi.org/10.1007/s00267-004-0095-9.
World Wildlife Fund
.
2020
.
The impacts of shipping on benthic habitats
.
World Wildlife Fund
.
Available at
https://wwf.ca/wp-content/uploads/2021/02/WWF-MPA-3-Impacts-Benthic-Habitat-v4.pdf.
Accessed January 23, 2023
.
Young
,
OR
,
Berkhout
,
F
,
Gallopin
,
GC
,
Janssen
,
MA
,
Ostrom
,
E
,
Van Der Leeuw
,
S.
2006
.
The globalization of socio-ecological systems: An agenda for scientific research
.
Global Environmental Change
16
(
3
):
304
316
. DOI: http://dx.doi.org/10.1016/j.gloenvcha.2006.03.004.

How to cite this article: Dreujou, E, Beauchesne, D, Daigle, RM, Carrière, J, Noisette, F, McKindsey, CW, Archambault, P. 2023. Multiple human activities in coastal benthic ecosystems: Introducing a metric of cumulative exposure. Elementa: Science of the Anthropocene 11(1). DOI: https://doi.org/10.1525/elementa.2023.00024

Domain Editor-in-Chief: Jody W. Deming, University of Washington, Seattle, WA, USA

Associate Editor: Laurenz Thomsen, Department of Marine Sciences, University of Gothenburg, Gothenburg, Sweden

Knowledge Domain: Ocean Science

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International License (CC-BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. See http://creativecommons.org/licenses/by/4.0/.

Supplementary Material