BRIEF RESEARCH REPORT article

Front. Remote Sens., 22 November 2022

Sec. Multi- and Hyper-Spectral Imaging

Volume 3 - 2022 | https://doi.org/10.3389/frsen.2022.1020184

A GEE toolkit for water quality monitoring from 2002 to 2022 in support of SDG 14 and coral health in marine protected areas in Belize

  • 1. Department of Civil and Environmental Engineering, University of California, Los Angeles, Los Angeles, CA, United States

  • 2. Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, United States

  • 3. Department of Geography, University of Georgia, Athens, GA, United States

  • 4. Wildlife Conservation Society Belize, Belize City, Belize

  • 5. Faculty of Science and Technology, University of Belize, Belmopan, Belize

  • 6. Earth System Science Center, University of Alabama in Huntsville, Huntsville, AL, United States

  • 7. Department of Atmospheric and Earth Science, University of Alabama in Huntsville, Huntsville, AL, United States

  • 8. Coastal and Marine Data Centre, Coastal Zone Management Authority and Institute, Belize City, Belize

Abstract

Coral reefs are highly diverse ecosystems that provide many goods and ecosystem services globally. Coral reef ecosystems are also threatened by environmental stressors from anthropogenic sources and shifting climates. The United Nations Sustainable Development Goal 14 (“Life Below Water”) addresses the need to conserve and sustainably use the ocean, seas, and marine ecosystems, including reef systems. Belize’s coral reef system is the second largest in the world, providing sources of income to Belizeans through tourism and fisheries as well as coastline protection. In order to conserve their marine ecosystems, Belize has a network of Marine Protected Areas (MPAs) throughout their coastal waters. Using Aqua MODIS satellite imagery from 2002 to 2022, Google Earth Engine, and RStudio, we present a workflow to calculate stress days on MPAs and a coral vulnerability index based on sea surface temperature (SST) and Kd (490), a proxy of water clarity. The Corozal Bay, Swallow Caye, Port Honduras, and South Water Caye MPAs had the highest percentages of stress days and coral vulnerability stress index score based on these two parameters among the 24 MPAs analyzed. Additionally, SST in the warmest month of the year in Belize were seen to increase across all MPAs from 2002 to 2022 (p < 0.01). This GEE toolkit provides a straightforward and accessible tool to help governments monitor both water quality and risks to coral reefs in accordance with SDG 14.

1 Introduction

Coral reefs have substantial cultural, economic, and environmental value (; ; ). These reefs provide a host of ecosystem services through tourism, biodiversity, fisheries, and coastal protection (). Coral reefs are estimated to have an asset value of approximately $1 trillion dollars, goods and services valued at over $375 billion dollars per year, with benefits reaching around 500 million people in at least 90 nations (). International efforts must continue to conserve these unique habitats, in light of climate change. The 2030 United Nations Agenda for Sustainable Development set forth 17 Sustainable Development Goals (SDGs) which include the most recent international goals for the sustainability and protection of oceans (; ), succeeding the 2015 Millennium Development Goals. SDG 14 entitled “Life Below Water,” most closely outlines targets and indicators for the protection of coral reefs (). In particular, SDG indicators 14.1.1a, 14.2.1, 14.3.1, and 14.5.1 relate directly to the health of coral reefs and surrounding water, as well as coverage of marine protected areas (). Descriptions of these indicators can be found in Supplementary Table S1. In 2021, the United Nations Environment Programme (UNEP) published an updated manual on measuring various indicators of SDG 14 (). This manual provides a comprehensive guide to implementing indicators via a subsection of monitoring parameters and methods under each indicator. However, for many parameters, such as water and habitat quality, ecosystem health, microalgal growth, and management effectiveness of protected areas, the UNEP states that these measurements should be taken when “national capacity to do so exists,” but does not provide any available data or methods on the parameter (). Therefore, there is a need for further research addressing feasible tactics for measuring these SDG 14 indicators on a national scale. The health of coral reefs may also affect the other 16 SDGs through direct benefits to the economy, indirect benefits to society, and general governance (). For example, SDG 12 entitled “Responsible Consumption and Production,” includes the consumption of ecosystem services such as tourism and fishing, which are both directly impacted by the health of coral reefs. SDG 8 entitled “Decent Work and Economic Growth” includes job incomes and economic sectors also relying on tourism and fishing. SDG 2, “Zero Hunger,” is impacted by communities that rely on fish catch, which in turn relies on coral reef health. Obura’s paper assesses the interactions between coral reefs and each SDG, determining that at least 11 of the 17 SDGs are potentially or strongly affected by coral reefs ().

Earth observation data has been vital for assessing coral reef health and coverage (; ; ). Two important water quality parameters to monitor for coral reef health include sea surface temperature (SST) and turbidity (). SST is important to measure since anomalous warm temperatures can lead to coral bleaching, coral disease, coral death, loss of coral cover, and shifts in biodiversity (; ; ; ). Aqua MODIS SST imagery has been used for coral reefs monitoring (; ). Daily SST remote sensing products are used in nowcast prediction methods such as calculating HotSpot, Degree Heating Weeks, and ReefTemp (). In the literature, Degree Heating Weeks (DHW) is often used for assessing heat stress on coral. The NOAA Coral Reef Watch DHW product is widely used (; ; ; ), but some studies have utilized Aqua MODIS SST to calculate DHW as well (; ).

Water clarity or a measure of the diffuse attenuation coefficient is also important to map since corals depend on light penetration to grow () among other water quality parameters such as turbidity, nutrients, and sedimentation (). Kd (490) products are used as a proxy for water clarity (; ; ). The use of water clarity is important to use in addition to thermal stress data, especially in places such as Belize, where thermally induced bleaching events may play a lesser role compared to anthropogenic causes (). The Aqua MODIS Kd (490) product using the NASA operational algorithm has been used in global reefs (; ), Brazil (; ; ), and Colombian Caribbean (). DWH and Kd (490) have been used in combination for coral health (; ).

Although literature on remote sensing applications for water quality monitoring is prevalent, there is a lack of research specifically using SST and Kd (490) to measure progress toward SDG 14 indicators. A systematic literature review of papers mentioning terms related to remote sensing, water quality monitoring, SDGs, Kd (490), and SST showed results for only six published papers (Supplementary Table S2). Of the six papers, only four were relevant to the criteria after closer analysis. In one of these studies, in situ samples validated Sentinel-2 MSI satellite data showing decreased water quality in Vembanad Lake, India, after the demolition of four high rise buildings on the shore of the lake (). This study measures sea surface temperature along with other variables such as salinity, DO, and pH. It also emphasizes that its use of lake water quality monitoring has potential to support SDGs 3, 6, 10, and 14. However, no specific indicators are mentioned and SDGs are only mentioned briefly in the paper. Additionally, this paper focuses on the involvement of the public in water monitoring rather than advancing remote sensing applications. Another study used ArcGIS 10.3.1 to map the spatial distribution of in-situ groundwater quality in Thatta, Sindh, in support of SDG 6. The study collected samples from pumps within the district that were then analyzed for parameters such as pH and turbidity (). Variation for each parameter was then interpolated using the “Kriging” tool to visualize water quality indicators throughout the entire district. Although this study used spatial variation to monitor groundwater, no satellite data were used, and SDGs were only mentioned briefly once. A third study targeting Vembanad Lake compared in situ water quality measurements collected from scientists with samples and observations determined by citizens through a mobile application (). The study suggests that the future use of community derived data could be used to validate satellite data, and emphasizes its potential for contribution to SDG indicator 6.3.2, “proportion of bodies of water with good ambient water quality.” A final study proposed a framework involving geo-spatial maps created via remote sensing data and GIS techniques to compare water quality parameters in 13 districts of the Uttarakhand state in India in support of SDG 6 (). Satellite data were applied to machine learning methods such as the random forest model in support of a water quality index. Turbidity was among the water quality parameters observed in this study, and SDG 6 was mentioned briefly as a potential use for the water quality index proposed. None of the studies in the criteria-based literature search specifically targeted monitoring of both SST and Kd (490) through satellite data in support of particular SDG 14 parameters. Therefore, further research on the use remote sensing water quality indices in support of specific SDG parameters is necessary for the expansion of accessible water quality monitoring.

Monitoring coral health through the synthesis of SST and Kd (490) data can directly address indicator 14.1.1 (a), Index of coastal eutrophication, since coral decline is often a direct result of eutrophication (). This index can also address indicator 14.2.1, Proportion of countries using ecosystem-based approaches to managing marine areas, as it provides an accessible way of monitoring both water quality and ecosystem health that can be applied to other water systems. Additionally, because acidification reduces the skeletal density of corals and makes them more prone to deterioration, this index could also be used to validate ocean acidification, in relation with SDG indicator 14.3.1, Average marine acidity (pH) measured at agreed suite of representative sampling station (). Specifically, increasing SST increases the rate of benthic respiration while lowering the ratio of productivity to respiration, therefore increasing the rate of CaCO3 sediment dissolution and leading to a lower ocean pH (). Indicator 14.5.1, coverage of protected areas in relation to marine areas, can also be benefitted by a SST and Kd (490) index, as the remote assessment of coral health would allow less accessible and larger areas to be monitored. Therefore, this study provides an important, accessible workflow model for supplementing monitoring of water quality and coral reef habitat in support of at least four indicators of SDG 14 in data scarce areas.

In this study, we created a straightforward workflow based on Google Earth Engine (GEE), RStudio, and Aqua MODIS-derived SST, DHW, and Kd (490) from 2002 to early 2022 and applied it to Marine Protected Areas (MPAs) in Belize. Region-specific thresholds were used for both SST and Kd (490) data to indicate stress days for all MPAs. A coral vulnerability index was created based on SST and Kd (490) to indicate MPAs that have experienced both turbidity and temperature stress in the last 20 years. Statistical analysis was also performed on the warmest months in Belize and visualized as maps.

2 Data and methods

The overall workflow and accompanying code for each step is outlined in Figure 1.

FIGURE 1

2.1 Study sites

Belize is a country located in Central America and home to the Belize Barrier Reef Reserve System (BBRRS). The BBRRS is the second largest barrier reef in the world and was inscribed as a UNESCO World Heritage Site in 1996 (; ). The reef system is around 250 km in length between Mexico and Guatemala. The distance between the reefs and the mainland ranges between 0.5 and 80 km (). The reef system contains reef patches, faros, fringing reefs, cays, and atolls (; ). Marine species within the reef include fish, manatees, invertebrates, and multiple species of sea turtles (; ; ). The coral reef needs clear water and consistent temperature regimes to thrive (; ). Global climate change (), pollution (; ), mining and dredging (), marine transportation (), algal blooms (), and overfishing stand to threaten the reef ecosystems of the BBRRS (). Recreational tourism is also seen as a threat to the reef system () which accounts for over 40% of the nation’s GDP ().

MPAs are managed marine environments with the purpose of conserving biodiversity (). MPAs in Belize were established in the early 1980s beginning with the Half Moon Caye National Monument. Later through community lobbying, more areas were added to the MPA network, such as the Hol Chan Marine Reserve in 1987 (). Eventually an integrated approach was necessary to account for land-based pollution from outside the bounds of the MPAs that threatened the health of the reef (). showed that the MPA network alone had not been enough to promote the restoration of reefs in Belize. In this study, we analyzed 24 MPAs spanning the Belizean coastal lagoon as of 2020, comprising a total area of around 5,000 km2, as shown in Figure 2. Each MPA contains a variety of habitat designations, including fishing and spawning, general use, mangrove, marine and coral reefs, preservation, and special management, with some study areas becoming MPAs as early as 1982 and others as recently as 2020.

FIGURE 2

2.2 Satellite imagery

Images from the Moderate Resolution Imaging Spectroradiometer (MODIS) on the Aqua satellite were accessed through GEE. Images included in the study are Level-3 daily data with a 4 km resolution spanning from 4 July 2002 to 28 February 2022. These images were used to calculate water clarity and sea surface temperature (SST), a combination of water quality parameters used to assess coral health ().

The average vertical diffuse attenuation coefficient for downwelling irradiance at 490 nm, Kd (490), was used as a proxy for water clarity. The NASA operational algorithm was used to calculate Kd (490) (), as we used in our previous study of water clarity in Belizean coastal waters (). This Kd (490) product has been used for coral health in Brazil (; ; ), for an analysis of global corals (), and corals in the Colombian Caribbean ().

Daily SST averages for each MPA were calculated using Aqua MODIS imagery. The MODIS SST product is commonly used to assess temperature stress on corals and impacts of SST anomalies. In addition, degree heating weeks (DHW) were calculated for all MPAs using the Aqua MODIS SST images. DHW is a common unit of measure to understand long-term thermal stress on corals (; ; ). DHW calculates accumulated thermal stress over a 12-week (84 days) period. Units are in °C—weeks, where 1°C—weeks is a week of SST over the maximum monthly mean. Instead of using a maximum monthly mean, a region specific threshold for coral bleaching in Belize (29.7°C) was used instead (; ). The region specific threshold was originally calculated for Channel Cay based on a 15 years record of NOAA/NASA AVHRR (Advanced Very High Resolution Radiometer) Oceans Pathfinder data (). Using this method, the DHW for a day is calculated over the 12-week running window including that day. The difference between the daily SST and Belizean threshold was calculated and retained for summation whenever the difference was greater or equal to 1°C, where it was then multiplied by 1/7 for the units to be in weeks. The factor of 1/7 is necessary since the development of coral bleaching occurs in the order of weeks (). DHW was calculated using the dbcaDHW package (https://github.com/dbca-wa/dbcaDHW/) in RStudio through the following formula:

2.3 Water clarity and sea surface temperature stress classifications

A literature review on Kd (490) thresholds and SST was conducted to find appropriate stress limits for corals in Belize and the Caribbean. A global coral study found that a Kd (490) value of 0.127 m−1 is too turbid for coral growth from performing a nonlinear least squares regression on the turbidity gradient within the inner Great Barrier Reef (). For Kd (490), “turbid” months were classified based on a coral reef diversity study in the Mesoamerican Reef at 0.30 m−1 (). Very turbid water was defined as Kd (490) values above 0.50 m−1 based on a Caribbean study that used this value to indicate Kd (490) anomalies (). The following was used for Kd (490) value classifications: 0–0.127 m−1 “Little/No Stress,” 0.127–0.3 m−1 “Turbidity Stress,” 0.3–0.5 m−1 “Turbid Month,” 0.5 m−1+ “Very Turbid.” The following are the DHW classifications: 0–1°C—weeks “Little/No Stress,” 1–4°C—weeks “Temperature stress,” 4–8°C—weeks “Bleaching Risk,” 8°C—weeks + “Mortality Risk.” Stress days are calculated as days where SST is greater than 29.7°C and Kd (490) is greater than 0.127 m−1.

2.4 Coral vulnerability index

Indices are often created in order to assess coral health (; ) as well as overall habitat protection (). Z-scores have also been used previously to normalize data for corals (; ). Here the SST and Kd (490) data were normalized and combined to create a coral vulnerability index which was calculated for each MPA (; ). First, the z-score for each parameter was calculated for all MPAs through the following formula:Here x is the mean per MPA per year, μ is the mean across all MPAs per year, and σ is the standard deviation across all MPAs per year. The median z-score was calculated for each MPA from all years. The median z-scores were then assigned an ordinal value from 1 to 6 based on the general distribution of parameters (). The two variables were assumed to have equal weight and were added for each MPA in order to compute the coral vulnerability index from 2 to 12, where higher values indicate a higher vulnerability due to high SST and Kd (490). The following are the assigned values for the corresponding median z-score range: z > 0.4 to “6,” z = 0.2–0.4 to “5,” z = 0–0.2 to “4,” z = −0.2–0 to “3,” z = −0.4 to −0.2 to “2,” and z < −0.4 to “1.”

2.5 Statistical analysis

All subsequent visualization and calculations were performed in RStudio (). To account for dependence between SST temperature values, or how dependent each day’s SST was on the SST of the prior day, the lme4 package was used in R to fit the September SST dataset into a linear mixed effect model (). Linear fixed effect models can be favorable over linear regression in downscaling climate variables (), and have been successfully used by other studies analyzing climate trends such as warming sea temperatures (), decreasing sea ice cover (), and declining carbon sinks in the Amazon ().

3 Results

DHW was categorized based on the potential for accumulated temperature stress to impact corals. Figure 3A depicts SST and DHW where they were both plotted from 2002 to 2022 for select MPAs. Ten MPAs had most of their DHW measurements in the “Little/No Stress” category and occasionally have incidence with more severe risk. Around one-third of the MPAs had incidences of DHW values in the “Bleaching Risk” category and the remaining one-fourth had DHW values in the “Mortality Risk” classification. Port Honduras had the highest incidences of DHW values in the “Mortality Risk” category at 10% of all calculable DHW. In almost all MPAs, an increase in severity in DHW is found following 2016. For SST, seasonality is visible with September being the warmest month across all MPAs (Supplementary Figure S1). Higher SSTs are particularly shown in 2008, 2017, and 2020 (Figure 4). Box plots and linear regression of the September SSTs show increases with time (Supplementary Figure S2).

FIGURE 3

FIGURE 4

Similarly to DHW, Kd (490) values were categorized based on thresholds and plotted as shown in Supplementary Figure S3 for all MPAs. Select plots of MPAs are in Figure 3B. In general, most MPA did not have water clarity stress, and few had consistently high values. Corozal Bay displays significantly higher Kd (490) values compared to all MPAs with most of their classifications as “Very Turbid” for the Caribbean region. Swallow Caye also had most of the Kd (490) values within the “Turbidity Stress” classification or higher and only a handful of days where average Kd (490) would have little stress on corals. Port Honduras and South Water Caye both display seasonal trends with their Kd (490) values where higher classifications of stress occur during Belize’s wet season. Other MPAs such as Laughing Bird Caye and Sapodilla Cayes Marine Reserve display seasonal trends in Kd (490) values, but all within the “Little/No Stress” category.

Using a z-score approach, sites were characterized for heat, turbidity, and a coral vulnerability index combining the two. Supplementary Table S3 lists the median z-scores and indices for SST and Kd (490) separately and the combined coral vulnerability index. Port Honduras and Sapodilla Cayes Marine Reserve have the highest SST index values while Swallow Caye and Corozal Bay have the highest Kd (490) index values. When the variables are combined for the overall index, Port Honduras is the most vulnerable, followed by Swallow Caye and Sapodilla Cayes Marine Reserve. Swallow Caye, Corozal Bay, Port Honduras, and South Water Caye had higher Kd (490) indices. Port Honduras, Sapodilla Cayes Marine Reserve, Glover’s Reef, Laughing Bird Caye, Hol Chan, Gladden Spit, and Gladden Spit and Silk Cayes had higher SST indices.

In accord with a linear fixed effect model (see Section 2.5), all 24 MPAs were assigned as random effects while year was assigned as a fixed effect in order to assess the temporal and regional effects on SST changes. Results depicted in Supplementary Tables S4, S5 show a highly significant SST increase over time (p < 0.0001) for the month of September and all data.

4 Discussion

There are international efforts being made toward the protection and conservation of coral reefs, including the “Life Below Water” goal outlined in UN SDG 14. As a result, scientists have created indices and use remote sensing to monitor environmental stressors impacting coral health. The accessible GEE toolkit presented in this work was applied to identify locations at risk for coral decline in Belize as a proof of concept for its use more broadly. Several locations in the region are shown to have multiple stressors for corals, as has been observed previously at other locations. In 2008, a coral susceptibility model was created using environmental stressors including SST, chlorophyll-a, solar radiation, and other parameters for corals in the western Indian Ocean (). In this study, the north-western regions of the study area had high vulnerability which may be due to high SST and photosynthetically active radiation, where high solar irradiance at the surface indicates potential for heating and photochemical damage. A remote sensing coral stress index was created for corals in Saudi Arabia using the Quickbird satellite, environmental stressors, and water depth (). The areas with higher coral stress index values were areas that were near large towns and cities with fishing pressure accounting for the majority of variation. Most recently, a Google Earth Engine based tool used a combination of stressors and reducers to find strong correlations between their stress exposure score and El Niño bleaching events for corals in the Red Sea, Chagos Archipelago, and Gilbert Islands (). While this tool extensively uses SST in their score along with other factors like wind, they do not include water clarity as the toolkit presented here.

Using our toolkit, we saw that nearshore MPAs such as Swallow Caye, Corozal Bay, Port Honduras, and South Water Caye had the highest z-scores based on Kd (490) alone. This is supported by previous studies that show physical connectivity as evidence through chlorophyll-a between land and the BBRRS (), sewage plumes off Caye Caulker, Belize (), and excess nutrients near Belize City (). Deviating from the typical use of thermal data for indices, Canto et al. created a light-based index for the Great Barrier Reef based on MODIS Kd data. Using a color-coded scoring system, similar to the methodology in our workflow, the study found strong correlations between “wet” years and discharge data with water clarity for inshore locations (). If runoff and discharge indeed play a major role in water clarity patterns nearshore, then we may expect these patterns to change in the future as the southern region of Belize is projected to experience less precipitation and runoff by 2090 ().

employed remotely sensed data in the calculation of coral stress indices including DHW for five major Indian coral reef regions from 1982 to 2018 (). successfully integrated remotely sensed SST data into a coral health monitoring program in the Eastern Arabian Sea. Large scale bleaching events resulting from the marine heatwave associated with the 2014–2016 El Niño-Southern Oscillation were indicated by bleaching indices relying on remotely sensed data and corroborated with underwater coral health surveys. DHWs for this study site were 4.80, 5.09, and 6.92 for the years 2014, 2015, and 2016, respectively. These values are comparable with max DHW values for vulnerable Belizean MPAs such as Corozal Bay (3.17 in 2014 and 5.25 in 2015 and 2016) and Port Honduras (6.08 in 2014, 6.65 in 2015, and 8.34 in 2016). found DHW ranged between 3.80 and 6.55 in Belize in 2015. The highest instances of Mortality and Bleaching Risk instances calculated with this workflow were greatest in 2002, 2003, 2006, 2011, 2012, and the years following 2017, which align well with years reported for having mass bleaching events in the Caribbean (). This shows that using Aqua MODIS-derived SST to calculate DHW is also comparable to NOAA Coral Reef Watch DHW. Future studies and uses of this workflow should verify these types of findings with coral bleaching data or El Niño-Southern Oscillation years.

The upward trend in September SST is in line with what is being observed globally due to climate change. SST has increased on average by 0.14°F each decade for the past 120 years (). Additional studies have shown especially high temperature increases in September, such as a paper comparing coral reef bleaching and SST over a 30 years time period in La Parguera, Puerto Rico (). This study demonstrated recurring severe coral bleaching in coincidence with sharp temperature rises in September, with an overall increasing trend throughout the assessed time period. SST has increased by at least 1.1°C since the pre-industrial area, and an estimated 60% of marine ecosystems have been degraded (). A SST increase of 1.5°C could threaten 70%–90% of coral reefs, while a 2°C increase is predicted to destroy nearly 100% of all coral reefs permanently (United Nations).

In 1998, the government of Belize passed the Coastal Zone Management Act to mitigate issues surrounding rapid development, overfishing, and rapid population growth. The Belize Coastal Zone Management Authority and Institute (CZMAI) is the leading authority in the nation’s management of coastal resources and is responsible for the development of the National Integrated Coastal Zone Management Plan (ICZMP). The goal of the ICZMP involves a set of recommended actions to ensure sustainable use of coastal resources with conservation and social and economic needs in mind (). The most recent plan released in 2016 developed projected scenarios for 2025 based on three approaches: conservation, informed management, and development. Since 2016, additional development activities such as dredging have occurred which may have effects on water clarity along with general climate change effects. CZMAI is currently looking to revise the 2016 ICZMP by looking at impacts of development on their coastline with specific interest on effects following the plan itself. This document is paving the way for the recently launched Marine Spatial Planning (MSP) process, renamed the Belize Sustainable Ocean Plan, which is a partnership between the Government of Belize and The Nature Conservancy. The MSP is part of a set of agreements that will enable the country’s debt conversion for marine conservation and will be legally enforceable. Under this plan, biodiversity protection zones will be increased from 15.9% to 30% of Belize’s open ocean ().

Through this work, we identified Port Honduras, Swallow Caye, Sapodilla Cayes Marine Reserve, and Corozal Bay as recommended MPAs for closer monitoring. In particular, Corozal Bay showed a high level of turbidity and heat stress consistently. Swallow Caye had turbidity stress most of the time, while Port Honduras and South Water Caye showed strong seasonal trends. This GEE workflow has the potential to be used in conjunction with in situ data to continue to monitor climatic and anthropogenic changes in water quality. The flexibility and dynamic nature of the toolkit will enable governmental entities to monitor new areas, especially with the expansion of protection zones under the new MSP.

For future research, it would be beneficial to add other potential variables as indicators of coral vulnerability that could be remotely detected such as chlorophyll-a measurements and toxic heavy metals concentrations. Chlorophyll-a is a key indicator of phytoplankton biomass, which can be used to assess the eutrophic status of water bodies, or the level of nutrients enriched within the water. Some papers have had success integrating chlorophyll-a as a variable among SST and turbidity to assess coral health (; ). Coastal areas such as the MPAs of Belize assessed in this study are especially important for chlorophyll-a monitoring because of their vulnerable ecosystems. Toxic heavy metals are also a major threat to coral health and their existence in marine ecosystems indicates a major concern for marine life. Even at low concentrations, heavy metals released into the oceans via anthropogenic activity such as mercury, lead, and arsenic can kill corals (). Existent studies sample heavy metal concentrations in marine water with respect to coral health (; ). Therefore, with further research, it should be possible to use remote sensing techniques such as those described in these studies to also map chlorophyll-a and heavy metal concentrations as additional indicators of coral health.

Prior to the creation of this workflow, coral indices and remote sensing products used to monitor water quality are mainly on thermal variables and geared toward well-studied reef systems such as the Great Barrier Reef. The workflow outlined here draws attention to certain Belizean MPAs in terms of conservation of coral health using freely-available data on SST and water clarity in Google Earth Engine and downstream analysis in RStudio. The accessibility of the data used and ability to replicate the analysis with ready to use tools and workflows is necessary to have robust management of coastal resources and understanding of climate impacts. Users are able to adapt thresholds to their own needs, add other stressors, and change the weights of the variable in the index. This toolkit is critical in areas that lack long-term measures of water quality parameters like Belize. Aside from the integration of additional coral stress indicators, future improvements of this toolkit could include comparisons of data with in situ samples for further validation. Because current marine sampling data is limited in Belize, it is necessary to collect more in situ data to assess water quality indicators, especially in the significant MPAs in this paper such as Port Honduras, Swallow Caye, Sapodilla Cayes Marine Reserve, and Corozal Bay.

5 Conclusion

The monitoring of water quality is extremely important for assessing the health of coral reefs and protecting coastal and marine resources. Remote sensing is a powerful tool for monitoring thermal stress and water clarity, two very important environmental stressors for coral. The work presented here used both Aqua MODIS-derived SST and Kd (490) with imagery from 2002 to 2022 to identify seasonal and long-term patterns in 24 of Belize’s MPAs. Coral stress was analyzed using two approaches, one based on previously identified thresholds for SST and Kd (490) for suitable conditions for corals, and one based on the distribution of these parameters among all sites in the study. With the use of this workflow, certain MPAs are revealed to suffer from either high SST, turbidity, or both, which was the case with Port Honduras. The hottest month in Belize is September, and an analysis of mean SST for that month over the time period of the study revealed statistically significant upward trends at all sites (p < 0.01). The easy accessibility and reproducibility of results for MPAs and corals in general can help other data scarce nations manage their coastal resources while striving to meet the UN SDGs. Future work should aim to involve other known environmental stressors such as chlorophyll-a as well as in situ data.

Statements

Data availability statement

The full workflow, codes, and data can be found on GitHub (https://github.com/iacallejas/belize-sdg).

Author contributions

IC, KO, JJ, DM, and CL wrote the manuscript. IC and KO developed the scripts and performed the data processing and analysis. CL, DM, RG, and EC conceived the study and CL, JJ, and DM co-advised the research. CL, RG, and JJ acquired funding. EC, AC, RG, NA, AR, and SR contributed to development of the project and manuscript editing.

Funding

This work was supported by the NASA RRNES (Grant #80NSSC20K1746) and NASA ROSES A.8 (cooperative agreement number #80NSSC19K0200), UCLA’s Center for Diverse Leadership in Science, the Joan Doren Family Foundation, NSF NRT: Graduate Traineeship in Integrated Urban Solutions for Food, Energy, and Water Management-DGE-1735325, and the California NanoSystems Institute. This work was performed in part at the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.

Acknowledgments

We would like to thank Wildlife Conservation Society Belize and Coastal Zone Management Authority and Institute for their support during the project.

Conflict of interest

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.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/frsen.2022.1020184/full#supplementary-material

References

Summary

Keywords

remote sensing, turbidity, Kd (490), Google Earth Engine, heat stress, coral vulnerability

Citation

Callejas IA, Osborn K, Lee C, Mishra DR, Auil Gomez N, Carrias A, Cherrington EA, Griffin R, Rosado A, Rosado S and Jay J (2022) A GEE toolkit for water quality monitoring from 2002 to 2022 in support of SDG 14 and coral health in marine protected areas in Belize. Front. Remote Sens. 3:1020184. doi: 10.3389/frsen.2022.1020184

Received

15 August 2022

Accepted

07 November 2022

Published

22 November 2022

Volume

3 - 2022

Edited by

Abel Ramoelo, University of Pretoria, South Africa

Reviewed by

Kgabo Humphrey Thamaga, University of the Western Cape, South Africa

Mahlatse Kganyago, University of Johannesburg, South Africa

Updates

Copyright

*Correspondence: Christine Lee,

This article was submitted to Multi- and Hyper-Spectral Imaging, a section of the journal Frontiers in Remote Sensing

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics