Seabed Morphology and Sedimentary Regimes deﬁning Fishing Grounds along the Eastern Brazilian Shelf

: Shelf morphology and sedimentary regimes are inﬂuenced by processes operating at different temporal and spatial scales and are important records of sea level changes and sediment supply and/or carbonate production. The northern continental shelf of Esp í rito Santo (Brazil) contains evidence of different sedimentary regimes that distribute diverse and complex marine habitats. Herein, seabed morphology, acoustic images of the seaﬂoor (side scan sonar and sub-bottom proﬁler), and sediment samples were used to investigate the inﬂuence of sedimentary regimes on physical marine habitat distribution. Seabed mapping is also integrated with available data on ﬁsheries to understand the inﬂuence of shelf morphology and sedimentology in the usage of distinct ﬁshing gears. The results indicate ﬁve morpho-sedimentary facies: terrigenous mud, terrigenous sand, rhodolith beds, carbonate gravel with rhodoliths, and hardground. Through an integrated analysis of the geomorphology and sedimentary distribution, two morpho-sedimentary domains were identiﬁed: a sediment-fed shelf adjacent to the Doce River associated with a major mud depocenter and a delta front morphology characterized by gentle slopes and low terrain ruggedness, and a sediment-starved shelf dominated by carbonate sedimentation showing an irregular morphology associated with higher slopes and terrain ruggedness. These contrasting morpho-sedimentary domains are a result of sedimentary responses to sea level ﬂuctuation during Late Quaternary, specially, during the deglaciation processes after the Last Glacial Maximum. The morphological and sedimentary contrasts along the area deﬁne the physical habitat distribution. The sediment supply regime area is associated with a terrigenous ﬁne/muddy sedimentation bed, which control the local morphology and favors coastal and delta front progradation. This physical habitat is a well-known shrimp-ﬁshing ground where intense trawling takes place, as well as gillnet ﬁsheries targeting weakﬁsh and croakers. The accommodation regime or low sediment inﬂux area is characterized by carbonate sedimentation associated with hardgrounds and rhodolith beds. In contrast, this physical habitat with scarce sediment supply, facilitates extensive benthic colonization by crustose coralline algae (CCA), which is primarily associated to line ﬁsheries, longlines, and spearﬁshing. Rhodoliths show a high diversity of CCA and the occurrence of an endemic kelp species. Long-term processes such as relative sea level ﬂuctuations and sediment supply are a legacy for the distribution of benthic habitats, and their resulting morphology can be a surrogate for predicting ﬁshing activities or a ﬁrst-base analysis for marine spatial planning. Available low-resolution bathymetric datasets can be a powerful tool, if applied with caution and in a regional scale approach. Here, terrain variables (terrain slope and ruggedness) derived from an extensive available (low-resolution and interpolated) bathymetric dataset distinguished two contrasting morphological domains characterized by rugged and smooth/ﬂat seabeds.


Materials and Methods
The investigation was carried out using a compilation of bathymetric data, and new acoustic and surface sediment data. The fishing ground information was obtained from workshops performed with five local communities and follows the same procedure adopted in [7]. The data were integrated and interpreted using a Geographic Information System (GIS) platform.

Study Area
The Espírito Santo continental shelf is located along the Eastern coast of Brazil. The shelf width varies from 50 to 200 km. Its widest part is known as the Abrolhos Shelf or the Abrolhos Bank [8][9][10]. To the south of the Abrolhos Bank, the shelf width averages 50 km and the shelf break occurs at 60 to 80 m water depth. The Doce river is the main sediment source to the shelf, with an average load of 133 × 10 6 t during the rainy season (November to April) [11]. Shelf morphology varies significantly with four distinct regions: Abrolhos shelf, Doce river shelf, Paleovalley shelf, and the wide-southern inner shelf [10]. Here, we focus on contrasting the Doce river shelf with the Paleovalley shelf (Figures 1 and 2).

Acoustic Imaging
Data collection was conducted using a Side Scan Sonar Edgetech 4100 with a 560p Digital Acquisition System, which included a 272 TD Towfish with Deep-Tow and a Dual-Frequency (100 and 500 Khz) Transducer coupled to a Garmin GPS and the Acquisition Software Discover 4100 ( Figure 1). Side-Scan Range was set in 200 or 300 m full swath.
Approximately 1600 km of side scan data were collected. Acoustic images were processed using the software SonarWiz 5 (Chesapeak Technology, Mountain View, CA, USA). Sonograms were corrected for bottom track and signal attenuation. Processed sonograms were exported as geotiff files (1m pixel resolution) and interpreted in a GIS platform.

Sub-Bottom Profiling
Sub-bottom profiles were collected to observe the morphology variations along the shelf using the surface reflector only and not for sub-surface or stratigraphic interpretations. Profiling data were acquired using a sub-bottom profiler, SyQuest StrataBox model, operating at 3.5 KHz. Data were received and recorded by the software StrataBox 3.0.4.1 (SyQuest Technology, San Jose, CA, USA) and processed using the software SonarWiz 5. A total of 1400 km of sub-bottom profile data were acquired.

Sediment Sampling
Sediment samples were collected using a Van Veen grab (n = 187 surficial sediment samples) ( Figure 1) and processed for grain size analysis and calcium carbonate content. A laser granulometer (Malvern Mastersizer 2000 (Malvern Panalytical, Nottingham, UK)) was used for particle size

Acoustic Imaging
Data collection was conducted using a Side Scan Sonar Edgetech 4100 with a 560p Digital Acquisition System, which included a 272 TD Towfish with Deep-Tow and a Dual-Frequency (100 and 500 Khz) Transducer coupled to a Garmin GPS and the Acquisition Software Discover 4100 ( Figure 1). Side-Scan Range was set in 200 or 300 m full swath.
Approximately 1600 km of side scan data were collected. Acoustic images were processed using the software SonarWiz 5 (Chesapeak Technology, Mountain View, CA, USA). Sonograms were corrected for bottom track and signal attenuation. Processed sonograms were exported as geotiff files (1m pixel resolution) and interpreted in a GIS platform.

Sub-Bottom Profiling
Sub-bottom profiles were collected to observe the morphology variations along the shelf using the surface reflector only and not for sub-surface or stratigraphic interpretations. Profiling data were acquired using a sub-bottom profiler, SyQuest StrataBox model, operating at 3.5 KHz. Data were received and recorded by the software StrataBox 3.0.4.1 (SyQuest Technology, San Jose, CA, USA) and processed using the software SonarWiz 5. A total of 1400 km of sub-bottom profile data were acquired.

Sediment Sampling
Sediment samples were collected using a Van Veen grab (n = 187 surficial sediment samples) ( Figure 1) and processed for grain size analysis and calcium carbonate content. A laser granulometer (Malvern Mastersizer 2000 (Malvern Panalytical, Nottingham, UK)) was used for particle size analysis of the mud fraction and the sieving method was used for the sand fraction. Calcium carbonate content was determined through dissolution with 10% hydrochloric acid.
Seabed images were obtained with a dropcamera consisting of a high-resolution camera (Go Pro Silver 3) attached to a steel frame, and looking vertically down. The base of the frame is a 60 × 60 cm 2 . analysis of the mud fraction and the sieving method was used for the sand fraction. Calcium carbonate content was determined through dissolution with 10% hydrochloric acid. The samples were described using the classification proposed by [12], based on carbonate content and grain texture, including: lithoclastic/terrigenous (CaCO3 < 30%), lithobioclastic (30 to 50% CaCO3), biolithoclastic (50 to 70% CaCO3), and bioclastic (>70% CaCO3). Lithobioclastic and biolithoclastic are described as mixed sediments.
Seabed images were obtained with a dropcamera consisting of a high-resolution camera (Go Pro Silver 3) attached to a steel frame, and looking vertically down. The base of the frame is a 60 × 60 cm 2 .

Figure 2.
Digital Terrain model of the study region. Continental data was obtained from Topodata database (http://www.dsr.inpe.br/topodata/), and bathymetrica data was digitized from sounding sheets provided by the Brazilian Hydrographic Office. The seabed morphology has a spatial resolution of 200m along the continental shelf. PVS: Paleovalley Shelf; DRS: Doce River Shelf. Datum is WGS 84.

Shelf Morphology Analysis
The bathymetric map is the result of a compilation of available data (printed sounding sheets and digital files) from the Brazilian Hydrographic Office. Ten sounding sheets were scanned, digitized, and gridded using Surfer 10 (Golden Software, Golden, CO, USA) and ArcGIS. A total of more than 150,000 points were digitized. Sounding sheets represent data collected by single beam echosounders from the 60's to the 80's. Positioning accuracy may vary due to changes in coordinate acquisition systems and methodologies (for more details see [10]). Original datum in the sounding sheets was Corrego Alegre that was transformed and corrected to WGS84, following the available digital data. Original sounding sheets scales varied from 1:10,000 to 1:160,000. All the bathymetric data obtained by the Brazilian Navy follows a protocol, so by not compiling data from other sources (such as satellite derived data), no mismatch was observed. To avoid mismatch, we limited our map to the 100 m isobaths (available sounding data). A data quality analysis was performed qualitatively considering the regional scale of the study. The only observed suspicious features were broad-scale

Shelf Morphology Analysis
The bathymetric map is the result of a compilation of available data (printed sounding sheets and digital files) from the Brazilian Hydrographic Office. Ten sounding sheets were scanned, digitized, and gridded using Surfer 10 (Golden Software, Golden, CO, USA) and ArcGIS. A total of more than 150,000 points were digitized. Sounding sheets represent data collected by single beam echosounders from the 60's to the 80's. Positioning accuracy may vary due to changes in coordinate acquisition systems and methodologies (for more details see [10]). Original datum in the sounding sheets was Corrego Alegre that was transformed and corrected to WGS84, following the available digital data. Original sounding sheets scales varied from 1:10,000 to 1:160,000. All the bathymetric data obtained by the Brazilian Navy follows a protocol, so by not compiling data from other sources (such as satellite derived data), no mismatch was observed. To avoid mismatch, we limited our map to the 100 m isobaths (available sounding data). A data quality analysis was performed qualitatively considering the regional scale of the study. The only observed suspicious features were broad-scale circular lines to the north of the study area. These features seem to be an artifact from the original distribution of the sounding data.
Terrain analysis was carried out using an ArcGIS toolbox, the Benthic Terrain Modeler-BTM [13]. The bathymetric map used for the morphometric calculations has a resolution of 200 m, following the dataset minimum resolution. Slope and Terrain Ruggedness were the input variables used in BTM. Slope represents the inclination of the seabed expressed in degrees and Terrain Ruggedness was calculated using the Vector Ruggedness measure-VRM. This variable represents terrain ruggedness as the variation in three-dimensional orientation of grid cells within a neighborhood, effectively capturing slope and aspect variability into a single measure [14,15]. Slope and VRM calculations is 3 × 3 pixels, following the toolbox default. The morphometric analysis was restricted to seafloor morphology, slope, and terrain ruggedness.

Shelf Geomorphology
The continental shelf presents a contrasting morphology characterized by a flat deltaic lobe with gentle slopes and an irregular relief marked by paleovalleys with greater slopes (Figures 2 and 3). A deltaic lobe is observed associated to the Doce River mouth. In this area, the shelf shows a flatbed (gradients < 0.2 • ) with a minor step at 30 m depth. Offshore, below 50 m water depth, an increase in seabed gradient (>0.8 • ) is observed where the shelf breaks. To the south, the inner shelf still presents a gently sloping seabed to Barra do Riacho, but the 30 m isobath gets closer to the coastline. In contrast, adjacent to the Piraque Açu River mouth, the seabed morphology is characterized by an irregular bottom with a wide variation in slope (up to 5 • associated with paleovalleys' wall). The shelf morphology in this latter region is characterized by paleovalleys and marine abrasion terraces. circular lines to the north of the study area. These features seem to be an artifact from the original distribution of the sounding data. Terrain analysis was carried out using an ArcGIS toolbox, the Benthic Terrain Modeler-BTM [13]. The bathymetric map used for the morphometric calculations has a resolution of 200 m, following the dataset minimum resolution. Slope and Terrain Ruggedness were the input variables used in BTM. Slope represents the inclination of the seabed expressed in degrees and Terrain Ruggedness was calculated using the Vector Ruggedness measure-VRM. This variable represents terrain ruggedness as the variation in three-dimensional orientation of grid cells within a neighborhood, effectively capturing slope and aspect variability into a single measure [14,15]. Slope and VRM calculations is 3 × 3 pixels, following the toolbox default. The morphometric analysis was restricted to seafloor morphology, slope, and terrain ruggedness.

Shelf Geomorphology
The continental shelf presents a contrasting morphology characterized by a flat deltaic lobe with gentle slopes and an irregular relief marked by paleovalleys with greater slopes (Figures 2 and 3). A deltaic lobe is observed associated to the Doce River mouth. In this area, the shelf shows a flatbed (gradients < 0.2°) with a minor step at 30 m depth. Offshore, below 50 m water depth, an increase in seabed gradient (>0.8°) is observed where the shelf breaks. To the south, the inner shelf still presents a gently sloping seabed to Barra do Riacho, but the 30 m isobath gets closer to the coastline. In contrast, adjacent to the Piraque Açu River mouth, the seabed morphology is characterized by an irregular bottom with a wide variation in slope (up to 5° associated with paleovalleys' wall). The shelf morphology in this latter region is characterized by paleovalleys and marine abrasion terraces.  Analysis of the terrain slope and ruggedness indicate two contrasting areas. The deltaic lobe shelf has a lower variability in slopes with low values of terrain ruggedness. Higher slopes are related to deltaic lobe fronts, dipping seaward. In contrast, the shelf with paleovalleys shows a higher value of terrain ruggedness, with higher slopes.

Morpho-Sedimentary Facies-Physical Habitats
The shelf presents a heterogeneous substrate in morphologically distinct seabeds. Five morphosedimentary facies (MSF) were identified ( Figure 4) based on shelf morphology and sedimentary facies: Terrigenous Mud, Terrigenous Sand, Bioclastic Gravel with Rhodoliths, Rhodoliths, and Coastal Hardground. These MSF are treated herein as physical marine habitats and are described as follows. Analysis of the terrain slope and ruggedness indicate two contrasting areas. The deltaic lobe shelf has a lower variability in slopes with low values of terrain ruggedness. Higher slopes are related to deltaic lobe fronts, dipping seaward. In contrast, the shelf with paleovalleys shows a higher value of terrain ruggedness, with higher slopes.

Morpho-Sedimentary Facies-Physical Habitats
The shelf presents a heterogeneous substrate in morphologically distinct seabeds. Five morphosedimentary facies (MSF) were identified (Figure 4) based on shelf morphology and sedimentary facies: Terrigenous Mud, Terrigenous Sand, Bioclastic Gravel with Rhodoliths, Rhodoliths, and Coastal Hardground. These MSF are treated herein as physical marine habitats and are described as follows.

Terrigenous Mud
This domain is primarily observed adjacent to the Doce River mouth, extending throughout the inner shelf; it was also observed in a small stretch of the middle and outer shelf ( Figure 4). This habitat comprised 29% of the study region. The Terrigenous Mud MSF represents the deltaic lobe morphology, with a gentle slope and extending up to 25-30 m water depth (Figure 3). It is composed basically of terrigenous mud and sandy mud. The lithoclastic composition along the coastal region near the Doce River mouth reflects the continental inflow, representing riverine sediments.

Terrigenous Sand
The terrigenous sandy bottom MSF comprises the largest habitat of the studied region (41%). This physical habitat is observed in the inner and mid shelf, but shows a quite irregular distribution ( Figure 5). Along the shelf adjacent to the Doce River and to its north, the sandy domain only occurs below 30 m water depth. Southwards, away from riverine mud influence, sandy bottoms occur along the inner shelf, in shallower waters. In morpho-sedimentary terms, this habitat is again characterized by lithoclastic sands, with grain sizes ranging from coarse to fine, but predominantly medium and fine sands.

Terrigenous Mud
This domain is primarily observed adjacent to the Doce River mouth, extending throughout the inner shelf; it was also observed in a small stretch of the middle and outer shelf ( Figure 4). This habitat comprised 29% of the study region. The Terrigenous Mud MSF represents the deltaic lobe morphology, with a gentle slope and extending up to 25-30 m water depth (Figure 3). It is composed basically of terrigenous mud and sandy mud. The lithoclastic composition along the coastal region near the Doce River mouth reflects the continental inflow, representing riverine sediments.

Terrigenous Sand
The terrigenous sandy bottom MSF comprises the largest habitat of the studied region (41%). This physical habitat is observed in the inner and mid shelf, but shows a quite irregular distribution ( Figure 5). Along the shelf adjacent to the Doce River and to its north, the sandy domain only occurs below 30 m water depth. Southwards, away from riverine mud influence, sandy bottoms occur along the inner shelf, in shallower waters. In morpho-sedimentary terms, this habitat is again characterized by lithoclastic sands, with grain sizes ranging from coarse to fine, but predominantly medium and fine sands. The occurrence of these sandy bottoms is also associated with a gentle slope and smooth flatbed along the shelf adjacent to the Doce River, but with a more irregular morphology in the inner shelf adjacent to the Piraque Açu river mouth. Remarkably, this domain is related to a second lobe-like morphology feature along the mid-outer shelf adjacent to the Doce River mouth. In this area, this MSF occurs offshore from the mud deposit (Figure 4).

Rhodolith Beds
The Rhodolith Beds MSF were consistently observed in the outer shelf up to the shelf break (Figure 4), composing 14% of the study area. This facies is represented by a gravelly carbonate fraction, composed from granules to rhodoliths. The extensive rhodolith area indicates that terrigenous sedimentation influence in this region is negligible, since the CCA that form rhodoliths are more likely to grow in areas with low sedimentary inflow [16]. The occurrence of these sandy bottoms is also associated with a gentle slope and smooth flatbed along the shelf adjacent to the Doce River, but with a more irregular morphology in the inner shelf adjacent to the Piraque Açu river mouth. Remarkably, this domain is related to a second lobe-like morphology feature along the mid-outer shelf adjacent to the Doce River mouth. In this area, this MSF occurs offshore from the mud deposit (Figure 4).

Rhodolith Beds
The Rhodolith Beds MSF were consistently observed in the outer shelf up to the shelf break (Figure 4), composing 14% of the study area. This facies is represented by a gravelly carbonate fraction, composed from granules to rhodoliths. The extensive rhodolith area indicates that terrigenous sedimentation influence in this region is negligible, since the CCA that form rhodoliths are more likely to grow in areas with low sedimentary inflow [16].
The occurrence of a rhodolith habitat is related to water depth and is adjacent to an irregular bottom morphology. Rhodolith beds are associated to a strong acoustic signature in side scan sonar data. This habitat predominates in depths >40-45 m, being associated with paleochannels and terraces, mainly along the southern part of the shelf. Locally, the formation of hardground or crusts by coalescence of CCA is also observed. Individual rhodoliths exhibited a typical nodular shape, with surface irregularities and diameters ranging between 1 and 10 cm.

Bioclastic Gravel with Rhodoliths
The bioclastic gravel with rhodoliths MSF composes 16% of the study area and is located along much of the shelf adjacent to the Piraque Açu river mouth, within the Costa das Algas Marine Protected Area (Figure 1). This domain is characterized by coarse fragments of biogenic material, including CCA and carbonate skeletons (including shell fragments). Rhodoliths occur along this domain, intermingled with bioclastic gravel. The seafloor is dominated by gravelly sediments and rhodoliths, forming a very irregular morphology, and is associated to paleochannels and hardground terraces.

Coastal Hardground
Coastal hardground is composed mainly by ferrigenous crusts (lateritic crusts) incrusted by organisms. It was only observed at the mouth of the Piraquê-Açu river, although it is a common feature in the intertidal zone along the southern coast of the area. These reefs are formed initially by sediment laterization, followed by marine abrasion during different cycles of sea level changes. Once a hardground, it becomes a substrate for incrustation.

Fishing Activities
Distribution of fishing grounds was obtained from workshops performed with five local communities. Five major fishery types were described (Table 1). The relationship between fisheries with distinct MSF and the percentage of each fishery activity over each habitat domain is shown in Figure 6. Most trawling (96%) and all bottom-associated gillnet fisheries are carried out in the terrigenous sediment flatbed domains. Over the carbonate domain with higher terrain ruggedness and slope, fisheries activities are related to line and longline fishing, as well as to spearfishing. Across the Doce river inner shelf, for example, the areas identified as muddy bottom habitats are associated with shrimp trawling artisanal fisheries [17], as well as bottom-associated gillnet fisheries targeting carnivorous fishes such as sciaenids, carangids, and elasmobranchs [18]. According to [19], from Barra do Riacho to the mouth of the Doce River, trawling is performed by both artisanal and industrial boats, mostly targeting shrimp. In a more recent survey from workshops with local fishermen, [20] also highlighted that in this southern region adjacent to the mouth of the Doce River, shrimp fishing is widespread due to the concentration of sandy-muddy bottoms, whereas on seafloors with a muddy composition, whiting fishing is more frequent. Along this stretch, the fishing area is continuous and is used from the coastline to the 40 m isobaths.

Discussion
The five MSF indicate the occurrence of two main sedimentary regimes: terrigenous supply regime and accommodation/carbonate sedimentation regime. The sedimentation pattern along the inner shelf adjacent to the Doce River is characterized by modern terrigenous sediment input, forming a submarine muddy delta lobe, representing a supply regime area. To the south and offshore from the terrigenous mud deposits, very low or no significant sediment input takes place. Shelf morphology is rather irregular due to erosion or resembles a relict morphological feature, characterizing an accommodation regime [5]. This area is characterized by relict/palimsestic sediments (terrigenous sands) and mixed sediments. Along the outer shelf and parts of the mid shelf, the sedimentary regime is predominantly carbonate, with development of bioclastic gravels and rhodolith beds. It is noteworthy that the carbonate sedimentation develops over an accommodation area characterized by low sediment input and relict morphological features, such as unfilled

Discussion
The five MSF indicate the occurrence of two main sedimentary regimes: terrigenous supply regime and accommodation/carbonate sedimentation regime. The sedimentation pattern along the inner shelf adjacent to the Doce River is characterized by modern terrigenous sediment input, forming a submarine muddy delta lobe, representing a supply regime area. To the south and offshore from the terrigenous mud deposits, very low or no significant sediment input takes place. Shelf morphology is rather irregular due to erosion or resembles a relict morphological feature, characterizing an accommodation regime [5]. This area is characterized by relict/palimsestic sediments (terrigenous sands) and mixed sediments. Along the outer shelf and parts of the mid shelf, the sedimentary regime is predominantly carbonate, with development of bioclastic gravels and rhodolith beds. It is noteworthy that the carbonate sedimentation develops over an accommodation area characterized by low sediment input and relict morphological features, such as unfilled paleochannels, formed during the last lowstand [10]. Thus, carbonate sedimentation is defined by biologically-controlled precipitation (e.g., crustose coralline algae-rhodoliths), but also by bioclastic particles transport and deposition.
The distribution of these regimes is influenced by processes that operate at different time and spatial scales. Typically, sedimentation is controlled by a combination of autogenic and allogenic processes, which determine the distribution of elements and geological variables in a depositional system [1,5]. Shelf morphology, at this scale, is a product of long-term processes, such as sea level changes.
The alteration of such controls introduces an imbalance between sediment supply and modifications within the basin, characterizing local sedimentary regimes and leading to depositional processes of aggradation or erosion, as well as coastal displacement [1,5,21]. Typically, if the rates of sediment supply are higher than sea level fall, the sedimentary regime is dominated by supply and accompanied by progradation or shoreline regression [1]. Conversely, under no, or minimal sediment supply, the trend is erosion, modification of the shoreline by wave action and, therefore, potential shoreline transgression. This process establishes a sedimentary regime dominated by accommodation accompanied by retrogradation of the coast and marine facies overlying the terrigenous facies [1,5,22].
The sedimentary facies distribution is characterized by lithoclastic (terrigenous), mixed, and bioclastic (carbonate) sediments, and is closely associated to seabed morphology ( Figure 6). Terrigenous mud and sand (primarily sandy mud and medium to fine sands) are observed along the deltaic gentle slope flatbeds, with low terrain ruggedness. A muddy patch is found adjacent to the river mouth. Areas with mixed, biolithoclastic, or lithobioclastic sedimentation are also identified; deposits were less common and might represent a transition from terrigenous to carbonate sedimentation. These sediments show a wide range of grain sizes, from sandy mud to gravelly sands. Bioclastic sediments comprise a higher concentration of gravelly and sandy fractions, which dominate the mid-outer shelf in the southern portion of the study area and the outer shelf in its northern part. Rhodolith beds occur along the outer shelf. These carbonate-dominated sediments are associated with higher terrain ruggedness, which is expressed by an irregular shelf morphology with marine abrasion terraces and paleovalleys.
Thus, combining spatial changes in bed slope, which represents the contrasting shelf geomorphology, with sedimentation regimes (terrigenous and carbonate), we can define two major Morpho-Sedimentary Domains (MSD): Terrigenous (terrigenous mud + terrigenous sand facies) and Carbonate (rhodoliths + bioclastic gravel with rhodoliths facies). The Carbonate MSD is related to higher slope, terrain ruggedness, and deeper depths, with carbonate-dominated facies, whereas the Terrigenous MSD is associated with gently inclined bottom, low ruggedness, and shallower areas with siliciclastic sediments. The two MSD are clearly observed in Figure 7, where long and crossshelf bathymetric profiles show changes in morpho-sedimentary characteristics. MSD distribution in relation to a terrain variable (slope) is shown in Figure 8.
In the Doce River area, the sedimentary input is dominated by supply, characterized by the progradation of the coastline. Typically, the deposits in the region are composed of fine to very fine terrigenous sediments. Such deposits are described as regressive deposits associated to high sediment input, promoting rapid accumulation. Other parts of the Brazilian shelf also exhibit a similar sedimentary regime and morphological characteristics, being dominated by deltas in between starved areas [23]. Four main deltas occur (from south to north): the Paraíba do Sul river delta, Doce river delta, Jequitinhonha river delta, and São Francisco river delta. Sedimentation patterns in the adjacent shelf of these deltas are very similar. Along the narrow eastern Brazilian shelf, [24] identified fine terrigenous sediment facies associated to river mouths or shelf incised valleys and carbonate dominated facies. For the São Francisco River shelf (AL), [25] described a progradation zone in the shelf, characterized by mud and muddy sands.
The presence of muddy bottom and terrigenous sand habitats in the study region suggests that the geological characteristics of the substrate drive demersal and benthic communities structured by scavengers. Substrates with unconsolidated marine deposits are typically associated with flat bottoms and a thick sedimentary cover [26], which was identified in both domains, particularly in the submarine delta lobes of the Doce River. Such morpho-sedimentary properties attract scavenging communities (e.g., shrimps, crabs, catfishes), since the thickness of the terrigenous sediments is directly associated with the levels of organic carbon supplied. In the coastal zone, this carbon is partially recycled and partially deposited with the terrigenous sediments [26]. Such sediments with high concentration of organic carbon provide an important food supply for high-biomass demersal and benthic communities.
bottoms and a thick sedimentary cover [26], which was identified in both domains, particularly in the submarine delta lobes of the Doce River. Such morpho-sedimentary properties attract scavenging communities (e.g., shrimps, crabs, catfishes), since the thickness of the terrigenous sediments is directly associated with the levels of organic carbon supplied. In the coastal zone, this carbon is partially recycled and partially deposited with the terrigenous sediments [26]. Such sediments with high concentration of organic carbon provide an important food supply for high-biomass demersal and benthic communities. Unlike the Doce River region, the southern section is characterized by a contrasting sedimentary regime (Figure 4). The continental shelf is deficient in sediment supply, with very low or negligible inflow, and is characterized by a sedimentary accommodation-dominated regime. The seabed habitats comprise coarse sediments, predominantly coarse sand, bioclastic gravel fractions, and extensive occurrence of rhodolith beds. These aspects support the transgressive characteristics in the area [10].
With the rapid rise in sea level, waves and tides eroded the coastline and shoreface substrates, which generated coarse particle deposits that were later modified and produced low sedimentation rates [5,22]. The low volume of accumulated and deposited sediments characterizes the inner and middle shelf of the region as a "starved" shelf with thinner deposits, coarse relict sediments, and carbonate sedimentation related to the observed habitats (rhodolith beds) [10].
Thus, the irregular morphology and the presence of paleovalleys are a result of marine transgressions under negligible sediment input. These unfilled paleochannels control the habitat Unlike the Doce River region, the southern section is characterized by a contrasting sedimentary regime (Figure 4). The continental shelf is deficient in sediment supply, with very low or negligible inflow, and is characterized by a sedimentary accommodation-dominated regime. The seabed habitats comprise coarse sediments, predominantly coarse sand, bioclastic gravel fractions, and extensive occurrence of rhodolith beds. These aspects support the transgressive characteristics in the area [10].
With the rapid rise in sea level, waves and tides eroded the coastline and shoreface substrates, which generated coarse particle deposits that were later modified and produced low sedimentation rates [5,22]. The low volume of accumulated and deposited sediments characterizes the inner and middle shelf of the region as a "starved" shelf with thinner deposits, coarse relict sediments, and carbonate sedimentation related to the observed habitats (rhodolith beds) [10].
Thus, the irregular morphology and the presence of paleovalleys are a result of marine transgressions under negligible sediment input. These unfilled paleochannels control the habitat morphology and likely act as sediment pathways from the shelf to the slope. A similar situation was recorded in the central continental shelf in Bahia, Northeastern Brazil [27], which was classified as a starving shelf. This classification was especially evident for the middle and outer portions, where several valleys dissecting the shelf were observed. In south Espírito Santo, based on surveys by [7], the accommodation system also dominates this portion of the shelf, which comprises low terrigenous mud and sand contents, while carbonate gravel and rhodolith dominate the shelf.  [27], which was classified as a starving shelf. This classification was especially evident for the middle and outer portions, where several valleys dissecting the shelf were observed. In south Espírito Santo, based on surveys by [7], the accommodation system also dominates this portion of the shelf, which comprises low terrigenous mud and sand contents, while carbonate gravel and rhodolith dominate the shelf.  The outer shelf of the studied area was dominated by an extensive rhodolith bed. The development of rhodolith beds along the eastern Brazilian shelf results from the rising sea level over the last 20 thousand years [28]. The development of the extensive accumulations of CCA on the tropical Brazilian shelf is attributed to ecological conditions that include low terrigenous input, good light penetration, efficient water circulation, and a relatively stable substrate [28,29]. Furthermore, the structural pattern of these rhodolith beds along depth gradients may also be related to the combined extension and inclination of the continental shelf [30,31]. In Abrolhos, for example, [32,33] the geological characteristics of the shelf, combined with meteoceanographic constraints, promote favorable conditions for rhodolith growth.
Rhodoliths highly influence the associated organisms and increase biodiversity, primarily by increasing structural complexity [30,34], being therefore considered ecosystem engineers [16]. The rhodolith structure provides a hard surface and three-dimensional substrate that serves as a microhabitat for a wide variety of associated invertebrates, algae and fish, many of which are economically and ecologically important [35][36][37]. Moreover, rhodoliths represent one of the world's largest deposits of calcium carbonate (CaCO 3 ), with an estimated 2 × 10 11 t [28,38].
The Terrigenous Mud and Terrigenous Sand shelter less diverse benthic assemblages. In the region near the coastal plain of the Doce River, the benthic habitat is not amenable to benthic marine macroalgae fixation, [37] due to the large continental sediment contribution that prevents settlement and development. In contrast, the region further south, adjacent to the Piraque-Açu River, is considered the richest Brazilian region in algal species, with CCA species endemism [39,40] in southern Espírito Santo. Bioclastic Gravel with Rhodoliths and Rhodoliths beds represent an important habitat for several commercially important and threatened fish species [41]. These ecosystems also provide an important habitat for benthic communities and maintain 25% of the known macroalgae species on the Brazilian coast [30]. Thus, it appears that the low sedimentary inflow in the Aracruz and adjacent shelves facilitates the colonization of vast rhodolith beds, conversely to the inner shelf adjacent to the Doce River. Recent studies highlight the presence of 74 fish species in mesophotic rhodoliths beds, including an endemic fish genus and many endangered species [42]. At least 18 commercially important species are associated with this habitat, most of them captured by lines and longlines, such as porgies (Pagrus pagrus and Calamus sp.), snappers (Lutjanus buccanella, L. cyanopterus and Ocyurus chrysurus), groupers (six species), and amberjacks (Seriola spp.) [42].
The shelf morphology and sedimentary distribution over the study area is strongly controlled by long-term processes, such as sea level changes and sediment input. Thus, fishing grounds are also controlled or related to this long-term response, i.e., a legacy of the Late Pleistocene-Holocene sedimentary evolution. Short-term and seasonal oceanographic processes are important, especially in terms of water column variability (temperature, salinity, turbidity, dissolved oxygen, primary production, etc.) and bed disturbance. However, the physical setting controlling the benthic habitat is the seabed morphology and sediment distribution, which is related to the geological background.
Seabed characterization is a major source of environmental information [7], such as environmentally relevant areas and fishing grounds. For instance, trawling activities can be less impacting in the northern areas where muddy bottoms are more widespread than in the southern area. Angling and gillnetting in the south area can help shift the currently high fishing effort over the shallow lateritic cuirasses towards less fished rhodolith bed and sand bottom. In this context, understanding the sedimentation regime is important because it is a key factor that directly influences resources' distribution. Moreover, shelf morphology and terrain analysis reflect sedimentation regimes in the study area and can be used as a primary surrogate by managers and decision-makers in a first step for marine spatial planning.
Finally, it is also relevant to highlight the importance of using available bathymetric datasets, even if they are old or interpolated, but reach minimum data quality. High-resolution multibeam bathymetric models are not commonly available for regional seabed analysis. In this case, available sounding sheets or regional digital single-beam data are the only sources for a morphological or geomorphological analysis of the continental margin. Here, we have used an extensive dataset with a minimum quality and a reasonable accuracy to conduct a terrain analysis from an interpolated digital terrain or bathymetric model. The available data and the derived map are powerful tools for a first or basic analysis for marine management, so the use of these sounding sheets became a must in Brazil and elsewhere. Results can be severely impacted by eventual inaccuracies of the data, such as, errors in data positioning due to changes in technology, datum transformation, or distortion from the printed sounding sheets, if caution is not applied. The results, either terrain variables or the digital terrain model, must be used with caution as limitations apply for seabed morphometric analysis. Usually, interpolated digital terrain models will be able to separate distinct morpho-sedimentary domains in a regional scale [43], but can be very inaccurate if one decides to undertake a high-resolution analysis. Herein, we applied the available dataset with caution and carried out a regional analysis distinguishing end-terms morpho-sedimentary domains; rugged from flat seabeds. The regional contrasting regimes and seabed morphology allowed this interpretation, showing that, in a regional scale, available bathymetric dataset can be very useful.

Conclusions
Morphological and sedimentological analysis of a sedimentary regime-contrasting continental shelf revealed the role played by long-term processes in marine physical habitat distribution. Areas with sediment supply regime exhibited physical habitats associated with terrigenous fine sedimentary facies with local prograding morphology. Terrigenous sedimentation was associated with gently inclined beds with low terrain ruggedness. Areas with very low sediment input represent the accommodation regime, in this case, a carbonate dominated shelf. Physical habitats are characterized by bioclastic gravels and living rhodoliths, occurring in deeper areas with irregular morphology, higher slopes, and terrain ruggedness. These distinct habitats explain fishery systems and exploited resources. Understanding the prevailing oceanographic and biological conditions is essential for the management of marine resources, but the study of the geological processes that drive sedimentation regimes allows for a broad understanding of the origin and distribution of marine habitats and fishery resources. In the studied region, the contrasting geomorphology is a reliable surrogate for fishing activities. Terrain analysis involving spatial variability in slope and ruggedness may be more widely used by managers and planners, not only as indicatives of sedimentary regimes, but also for mapping fishing grounds and for designing more effective biodiversity assessments. Moreover, if used with caution, available low-resolution bathymetric datasets can be applied for geomorphometric analysis in a regional scale approach. Here, terrain variables (terrain slope and ruggedness) derived from an extensive available (low-resolution and interpolated) bathymetric dataset distinguished two contrasting morphological domains characterized by rugged and smooth/flat seabeds.