Shallow reef ecosystems are increasingly becoming degraded from human activities, coastal development, and commercial activities. As human population grows and expands to coastal areas, threats to reef habitats will continue unabated. Additionally, many human activities increase the likelihood of environmental contamination through introduction of invasive species, further threatening the stability of coastal ecosystems  . To curb this trend, reef habitat conservation and recovery efforts need to employ proven techniques toward that goal, of which protected areas are likely the most promising.
Marine protected areas have proven effective in conservation of several ecosystems, including reef habitats    . Its widespread adoption by local, national, and international regulators attests to its effectiveness in conservation and recovery of coastal habitats for commercially important or emblematic species  . The benefits of MPAs in the conservation of non-target species and the maintenance of ecosystem functions are, however, still not well known  . For such techniques to be effective siting is one of the most important factors to consider   . Locations of MPAs need to include a network of high diversity sites to provide optimum viability to the persistence of such networks long term  . Toward objective, we investigate three areas in terms of diversity of non-target benthic biota to assist in marine conservation policies and recovery efforts of coastal ecosystems.
2. Materials and Methods
Research of benthic diversity was conducted in shallow reef habitats of the Rio de Janeiro and São Paulo states. The study regions comprised of an upwelling area to the north, an adjacent high impacted area, and a more distant marine park to the south (Figure 1). At each region, four sites were chosen and two depths strata at each site. Depth strata were less than 8 m and between 10 and 15 m. At each depth strata, two 20 m long transects were randomly selected for investigations. Each depth stratum was visited once every quarter from December 2013 to June 2015, totaling six data collection visits.
Data collections for biodiversity were done by SCUBA using quadrats 25 cm on a side. Twenty quadrats were placed randomly along each transect before a still digital photograph of each was taken. Photographs were processed to extract percent area covered by benthic organism. Organisms were identified to the lowest taxa possible during photograph processing. No live organisms were taken, to minimize habitat disturbance during data collection events. Data collection for substrate relief were conducted using a 5 m chain placed at random locations in each stratum. The chain was draped over rocks and boulders and the loss of length recorded (no relief implies in a stretched-out chain, i.e., no loss of the 5 m chain length). Loss of length was used as the estimate of substrate complexity. Water temperature data was obtained by placing one continuous temperature logger at each study region.
Data analysis comprised of estimates of biodiversity, substrate relief, and water temperature patterns for the study area. Estimates were obtained using independent analysis of variance (AOV) factorial designs. The Shannon index of
Figure 1. Study area and regions I-III (in circles) used in benthic biodiversity studies along the coast of Rio de Janeiro and São Paulo states, Brazil; Region I: upwelling area, Region II: high-impact area, Region III: marine park.
biodiversity was the response variable for the biodiversity analysis. Average monthly temperature was the response for water temperature pattern estimates. Loss of chain length was the response for the substrate complexity analyses. Region and stratum were the factors selected for AOV analyses. Sites were not included as factors because they had no ecological relevance aside from serving as replicates within study regions.
To avoid biases for biodiversity analyses, only organisms identified to genus or species were included in analyses. In case an organism was only possible to be identified to genus and it was the only species observed throughout the study, it was included in analyses and assigned a species denomination (sp.).
3. Results and Discussion
Fifty-one organisms were sampled at all regions combined (Table 1). The vast majority of the observations were of macroalgae, echinoderms, and poriphera, typical of shallow-water marine ecosystems. The region with most organisms was the upwelling region, followed by the high-impact region, with the marine park last. The shallow stratum always yielded the most organisms (Table 2),
Table 1. List of organisms identified and total counts during biodiversity studies at shallow-water reef habitats on the coasts of São Paulo and Rio de Janeiro states, Brazil; Region I: upwelling area, Region II: high-impact area, Region III: marine park.
Table 2. Mean species count by study region and depth stratum observed during biodiversity studies at shallow-water reef habitats on the coasts of São Paulo and Rio de Janeiro states, Brazil.
which is also expected, due to its higher exposure to sunlight.
The highest diversity was found in the upwelling and the adjacent high-impact regions, with no statistical difference between both. The marine park had the lowest diversity among all three regions (Figure 2). This may at first be counter-intuitive, since the high-impacted region is located in a metropolitan area, with detrimental impacts from human recreation, such as diving and spear-fishing, and pollution. As the high-impact area is south from the upwelling area, there may, however, be a strong overspill from the upwelling area into the high-impacted area, supplying the latter with biological resources. The major driver of the observed overspill effects may be the predominantly north-south current flow at the study region  . Overspill effects are often reported one of the many benefits of marine protected areas, as they do not only preserve habitat and resources locally, but may be a source for open, adjacent areas     .
The shallow stratum consistently had a higher diversity compared with the deep stratum (Figure 2). Unlike the results above by region, this is expected. The shallow stratum was generally around 5 m deep, compared with the deep stratum of around 12 m deep. Near coastal waters at the study area are rich in planktonic and non-skeletal suspensates, which impedes light from penetrating
Figure 2. Analysis of variance least square means ± standard error of Shannon biodiversity indices following biodiversity studies on the coasts of São Paulo and Rio de Janeiro states, Brazil; panel (a): factor Region, panel (b): factor Region X Stratum; panel (c): factor Stratum; lines indicate factor grouping significance (continuous lines under factor indicate non-significance).
higher depths  . What could be surprising is the rate of diversity losses with depth, placing a higher importance of shallow reef habitats in the maintenance of biodiversity.
The mechanisms driving the diversity patterns in the study area were not clear from the results obtained herein. Habitat complexity as indicated by chain length losses seems not to have been a factor driving diversity patterns, as no statistically significant differences were found in substrate relief among regions or depth stratum. Habitat complexity is known to affect diversity of aquatic habitats    , indicating that for this study other factors were at play in determining the observed biodiversity patterns.
One of such factors could be associated with water temperature, as there was a statistically significant difference of monthly average water temperature among region and between depth strata (Figure 3). Water temperature per se, however,
Figure 3. Analysis of variance least square means ± standard error of water temperature following biodiversity studies on the coasts of São Paulo and Rio de Janeiro states, Brazil; panel A: factor Region, panel B: factor Region X Stratum; panel C: factor Stratum; lines indicate factor grouping significance (continuous lines under factor indicate non-significance).
was likely not the driving factor, as biodiversity patterns did not parallel that of temperature. Water temperature was lowest for the deep stratum (low diversity) and highest at the marine park (low diversity). We, therefore, conclude that water temperature is an indicator of upwelling phenomena, serving as a surrogate variable to indicate such regions at our study area. Other factors, possibly more related to nutrient availability from upwelling currents, are possibly at play for determining the observed biodiversity patterns.
The implementation of protected areas has great promise for recovering degraded habitats and preservation of entire ecosystems. Siting of such areas, however, is of prime importance in determining the effectiveness of such conservation practices. Setting aside areas of great diversity may be preferable, as they constitute more effective sources for adjacent, open areas. There is, however, a trade-off, as high-diversity areas are also prime locations for recreation or commercial fishing activities   . We demonstrated in this study that high-diversity areas may benefit adjacent areas of high human impact. As the importance of intact ecosystems becomes more apparent, trade-offs between preservation and economic interests may increasingly tip toward preserving ecosystems, especially if the benefits of spillover effects are demonstrated, as done in this study.
Financial support for this project provided in full by the program Petrobras Socioambiental, Petrobras, Brazil. We also recognize and thank the Ecological Station Tamoios for their continued logistical and technical support throughout this project and Dr. Guilherme Muricy, from the Federal University of Rio de Janeiro National Museum, for his assistance in identification of benthic organisms.
 Kaplan, K., Hart, D., Hopkins, K., Gallager, S., York, A., Taylor, R. and Sullivan, P. (2017) Invasive Tunicate Restructures Invertebrate Community on Fishing Grounds and a Large Protected Area on Georges Bank. Biological Invasions, 20, 87-103.
 Ban, N., Davies, T., Aguilera, S., Brooks, C., Cox, M., Epstein, G., Evans, L., Maxwell, S. and Nenadovic, M. (2017) Social and Ecological Effectiveness of Large Marine Protected Areas. Global Environmental Change, 43, 82-91.
 Friedlander, A., Golbuu, Y., Ballesteros, E., Caselle, J., Gouezo, M., Olsudong, D. and Sala, E. (2017) Size, Age, and Habitat Determine Effectivenss of Palau’s Marine Protected Areas. PLOS ONE, 12, e0174787.
 Welch, H., Pressey, R. and Reside, A. (2018) Using Temporally Explicit Habitat Suitability Models to Assess Threats to Mobile Species and Evaluate the Effectiveness of Marine Protected Areas. Journal of Nature Conservation, 41, 106-115.
 Nicoll, R. and Day, J. (2017) Correct Application of the IUCN Protected Area Management Categories to the CCAMLR Convention Area. Marine Policy, 77, 9-12.
 Mellin, C., MacNeil, M., Cheal, A., Emslie, M. and Caley, M. (2016) Marine Protected Areas Increase Resilience among Coral Reef Communities. Ecology Letters, 19, 626-637.
 Cabral, R., Gaines, S., Lim, M., Atrigenio, M., Mamauag, S., Pedemonte, G. and Alino, P. (2017) Siting Marine Protected Areas Based on Habitat Quality and Extent Provides the Greatest Benefit to Spatially Structured Metapopulations. Ecosphere, 7, e01533.
 Cox, C., Valdivia, A., McField, M., Castillo, K. and Bruno, J. (2017) Establishment of Marine Protected Areas Alone Does Not Restore Coral Reef Communities in Belize. Marine Ecology Progress Series, 563, 65-79.
 Visconti, P. and Elkin, C. (2009) Using Connectivity Metrics in Conservation planning—When Does Habitat Quality Matter? Diversity and Distributions, 15, 602–612.
 Lima, I., Garcia, C. and Moller, O. (1996) Ocean Surface Processes on the Southern Brazilian Shelf: Characterization and Seasonal Variability. Continental Shelf Research, 16, 1307-1317.
 Andrello, M., Mouillot, D., Somot, S., Thuiller, W. and Manel, S. (2014) Additive Effects of Climate Change on Connectivity between Marine Protected Areas and Larval Supply to Fishes Areas. Diversity and Distributions, 21, 139-150.
 Colleter, M., Gascuel, D., Albouy, C., Francour, P., Morais, L., Valls, A. and Le Loch, F. (2014) Fishing Inside or Outside? A Case Study Analysis of Potential Spillover Effect from Marine Protected Areas Using Food-Web Models. Journal of Marine Systems, 139, 383-395.
 Di Lorenzo, M., Claudet, J. and Guidetti, P. (2016) Spillover from Marine Protected Areas to Adjacent Fisheries Has an Ecological and Fishery Component. Journal for Nature Conservation, 32, 62-66.
 Cuervo-Sanchez, R., Madonado, J. and Rueda, M. (2018) Spillover from Marine Protected Areas on the Pacific Coast in Colombia: A Bioeconomic Modeling Approach for Shrimp Fisheries. Marine Policy, 88, 182-188.
 Ding, B., Curole, J., Husemann, M. and Danley, P. (2014) Habitat Complexity Predicts the Community Diversity of Rock-Dwelling Cichlid Fish in Lake Malawi, East Africa. Hydrobiologia, 748, 133-143.
 Oyugi, D., Mavuti, K., Aloo, P., Ojuok, J. and Britton, J. (2014) Fish Habitat Suitability and Community Structure in the Equatorial Lake Naivasha, Kenya. Hydrobiologia, 727, 51-63.
 Martelo, J., Grossman, G., Porto, M. and Magalhaes, M. (2014) Habitat Patchiness Affects Distribution and Microhabitat Use of Endangered Mira Chub Squaliustorgalensis (Actimopterygii, Cypriniformes). Hydrobiologia, 732, 93-109.