First insights into an exceptionally deep blue hole in the Western Caribbean: The Taam ja’ Blue Hole

Coastal karst formations may remain unexplored by the global marine scientific research community. Knowledge provided by fishermen coupled with scientific research on karstic formations prompted us to explore the geomorphometric and physicochemical features of a submerged blue hole in the large tropical estuary of Chetumal Bay in the southeastern coast of the Yucatán Peninsula, Mexico. Differently from the well-known literature on cenotes in the Yucatán Peninsula, no references are made to blue holes in the region. Sampling and surveying of the blue hole were conducted in September 2021 by scuba, echosounders, CTD profilers, and the collection of water samples. This study revealed a previously unknown maximum water depth for the Taam ja’ Blue Hole, as it is now named, of 274.4 meters below sea level (mbsl), making it the second deepest known blue hole in the world. The Taam ja’ Blue Hole has a nearly circular shape at its surface, encompassing an area of ~13,690 m2. The sides of the blue hole are steep, with slopes >80° that form a large conic structure covered by biofilms, sediments, limestone, and gypsum ledges. The blue hole mouth locates at ~4.5-5.0 mbsl, at which water properties change significantly with temperature and salinity gradients of 4.9°C/m and 9.8 PSU/m, respectively. Hydrographic profiles show a stratified water column inside the blue hole consisting of a hypoxic layer (5-20 mbsl), a chemocline (50-80 mbsl), and an anoxic layer (>110 mbsl). Evaluation of anions and cations indicated that the water inside corresponds mainly to that of diluted seawater.


Introduction
Blue holes, submarine caves, and sinkholes are common features on coastal karst platforms around the world, including the Caribbean Sea and the Yucatań Peninsula, that they contain unique environments (Mylroie, 2008;Mylroie, 2013;van Hengstum et al., 2011;Little et al., 2021). Karstic systems are formed from the dissolution of limestone. Repeated flooding and draining of the karst lithology by sea-level oscillations during glacial and interglacial periods (van Hengstum et al., 2011;van Hengstum et al., 2019), combined with the higher salinity of coastal waters (Kambesis and Coke, 2013) enhanced the development of impressive karst formations, such as the extensive underwater cave system of 'Sac Actun' (Supper et al., 2009;Kambesis and Coke, 2013) or the numerous 'cenotes' (van Hengstum et al., 2008) along the Mexican Caribbean coast (Bauer-Gottwein et al., 2011;van Hengstum et al., 2020). Similarly, large and deep blue holes have developed worldwide, such as the Sansha Yongle Blue Hole in China (~300 m deep), the Dahab Blue Hole in Egypt (~130 m deep), and the Gozo Blue Hole in Malta (~60 m deep) (Li et al., 2018).
Blue holes in The Bahamas (Mylroie and Mylroie, 2007;Tamalavage et al., 2018;Little et al., 2021) and in the western Caribbean are also notable due to their remarkable structures, such as Dean's Blue Hole in The Bahamas (~202 m deep) (Vimpere, 2017) or the Great Blue Hole off the coasts of Belize (~125 m deep) (Jones and Dill, 2002;Gischler et al., 2013;Schmitt et al., 2021). These karst features are relevant to marine scientific research as they represent repositories of paleoenvironmental and paleoclimate information that enhance the understanding of environmental (Husson et al., 2018) and geologic processes , provide support for modeling the mixing zone (Bottrell et al., 1991), the connectivity of aquifers (Martin et al., 2012) and the water origin of these hydrogeologic systems (Smith et al., 2021). They offer windows into regional hydroclimate variability Wallace et al., 2021) and are galleries to observe subfossil benthic meiofauna and modern biodiversity (van Hengstum et al., 2011;Gonzalez et al., 2020;Baremore et al., 2021;Little et al., 2021).
Blue holes in the Western Caribbean have been less investigated (e.g. Gischler et al., 2013;Baremore et al., 2021;Schmitt et al., 2021), in comparison to the well-known literature on karstic structures from cenotes (inland sinkholes) in the Yucatań Peninsula (e.g. Perry et al., 2009;Bauer-Gottwein et al., 2011). However, a system of cenotes have developed near the border between Mexico and Belize with water depths in the range of 40-90 m (Perry et al., 2009;Bauer-Gottwein et al., 2011), together with submerged blue holes, reported by Carrillo et al. (2009b), with depths up to 42 m inside the shallow tropical estuary of Chetumal Bay. In September of 2021 an expedition to the location of a blue hole of unknown depth in the bay was carried out in collaboration with Jesus Artemio Poot Villa (personal communication, March 15, 2021), a local fisherman.
The information gathered revealed a maximum recorded depth of 274.4 meters below sea level (mbsl). This makes the now named Taam ja' Blue Hole (i.e., deep water in Mayan) (TJBH) the world's second known deepest blue hole. Submarine exploration was conducted through echo-sounding, thermohaline profiles, chemical analysis of collected water samples and scuba diving, leading to the first documentation of this geomorphologic feature. This paper aims to boost the research and to introduce to the scientific community a particularly deep blue hole in the Western Caribbean and likely the deepest known blue hole in the region.

Study area
Chetumal Bay is a semi-closed mesohaline tropical estuary of the Western Caribbean coast. It extends over a shared area between Mexico and Belize, with water depths of 2-5 m (Carrillo et al., 2009a) (Figure 1). It has several freshwater inlets and connections to adjacent lagoon systems and underground springs. Salinity and water temperature range between 25-32°C and 13-22 PSU, respectively (Carrillo et al., 2009b;Ruıź-Pineda et al., 2016). Microtidal conditions (<0.2 m) (Kjerfve, 1981) and restricted water exchange with the sea are present in Chetumal Bay, but a highly dynamic system develops due to wind action (Carrillo et al., 2009a;Avalos-Cueva et al., 2017).
The TJBH is located at 378823 m E, 2059390 m N (UTM 16Q) in the central portion of Chetumal Bay, where submerged coastal karstic sinkholes locally named 'pozas' have been reported (Carrillo et al., 2009b) ( Figure 1A). Following Mylroie (2013) and Yao et al. (2020), these geomorphologic features could be recognized as blue holes since they are roughly circular and steep-walled karstic formations, open to the surface but with their depth fully developed below sea-level, and filled with water of fresh, marine and brackish chemistry.

Fieldwork and sampling techniques
Fieldwork was conducted in September 2021 by scuba, the use of echosounders, CTD profiling, and chemical water sample extraction. Diving operations were undertaken to explore the blue hole along its south−southeast side, down along the wall and up to 20 mbsl. The route followed on scuba (~200 m) was parallel to the blue hole perimeter. Images and video recording were performed using an Action Cam EV-07005 (Evorok USA).
Measurement of the TJBH bathymetry and surrounding areas was conducted by means of a Humminbird 899cxi HD SI sonar unit with an integrated GPS (Humminbird USA), a 200/83 kHz DualBeam PLUS ™ sonar system and a side-scan sonar frequency of 455 kHz. This allowed for the detection of georeferenced water depths at the blue hole and for the seabed to be imaged and recorded. A fastCTD Profiler (Valeport UK) with an integrated Hyperion Valeport fluorometer (excitation: 470 nm, detection: 545 nm) was used to measure pressure, temperature, conductivity, and fluorescence along vertical profiles in the water column at a sampling frequency of 16 Hz. The profile's location was determined by GPS with a GPSMap ® 78 (Garmin USA). The CTD profiles were measured up to a depth of 200 mbsl to prevent instrument damage, but no bottom was reached.
Water samples were collected in triplicate at 2, 20, 50, 80 and 110 mbsl inside the TJBH with an Alpha 1120-G42 (Wildco USA) Van Dorn sampling bottle for the analysis of chemical constituents. The depth limit for water sampling was based on the capacity of the field equipment used. Further surface water samples were collected 2 mbsl at locations 500 m to the west and east of the TJBH center. At the same depths and locations, temperature, pH, conductivity, salinity, dissolved oxygen (DO), oxidation-reduction potential (ORP), and total dissolved solids (TDS) were also measured with a multiparametric probe (Horiba U-52) previously calibrated with standard buffer solutions (Hydrion Buffer Kit), obtaining a precision of ± 0.1°C for temperature, ± 0.01 mg/L for dissolved oxygen, and ± 0.01 for pH. The samples were collected from 0.5-L polypropylene bottles previously rinsed with 15% hydrochloric acid. The samples were transported and preserved at a temperature of< 4°C until analysis (APHA, 1998).

Morphometric analyses
The orthomosaic of the TJBH resulted from merging georeferenced raster images created with SonarTRX software and side-scan sonar recordings. A digital elevation model (DEM) of 0.5 m resolution was generated through QGIS software and the interpolation of the measured georeferenced water depths. A spatial moving average low-pass filter was also applied to the interpolated DEM considering a 5 gridding node window to reduce noisy data from possible acoustic scatter of soundwaves bouncing inside the blue hole structure. Bottom profile of the TJBH was extracted from the resulting DEM with QGIS tools. Similarly, lateral projection of the blue hole was obtained with ParaView software, based on the clipping of the three-dimensional model generated with DEMto3D plugin in QGIS software. Both the surface area and volume of the TJBH were also calculated with QGIS software. The surface area (A=∑p i ) was based on the summation of pixel areas (p i ) from the DEM raster file within the blue hole perimeter (P), which was visually estimated from the orthomosaic generated. The total volume of the blue hole (V=∑h i p i ) was determined by the summation of pixel areas (0.5 x 0.5 m) multiplied by the water depth information (h i ).
Morphometric indexes were described by the Gravelius coefficient (K C =0.28·P/A 0.5 ) and the ratio (R) between the major and minor axes for assessing the circularity and elongation proportion of the blue hole (Fragoso-Servoń et al., 2014). The major axis was identified by considering the longest distance between two points of the blue hole perimeter. The minor axis resulted from the longest distance perpendicular to the major axis.
Triplicates of certified reference material for nutrients (ERA 505, 770, 693, and 698 brand WasteWatR ™ ) were processed, obtaining recovery percentages between 95% and 105%. Based on Appelo and Postma (2005), ionic balance (IB) (eq. 1) was evaluated for the ions Ca 2+ , Mg 2+ , Na + , K + , HCO 3-, SO 4 2-, Cland NO 3considering a margin of error ± 5% to verify the reliability of analytical results. All reported data were within the acceptable range. A Piper diagram (Piper, 1944;Appelo and Postma, 2005) was generated from the anion and cation values to describe the physicochemical conditions in the blue hole water column.
The degree of chemical balance between water and minerals in the aquifer could be assessed through saturation indices (SI) (Yuan et al., 2017). Thus, the SI [log(Q/K)] for gypsum and anhydrite minerals were calculated using the PHREEQC software (Parkhurst, 1995), in which Q is the activity product, and K is the equilibrium constant for the reaction. An equilibrium zone of SI=0 ± 0.5 was considered due to the uncertainty in the precision of the data, where negative SI-values correspond to undersaturated water and positive SI-values to supersaturated one (Deutsch, 1997). The SO 4 -2 /Cland Na + /Clratios were also considered to determine the salinity and sulfate dissolution patterns with depth.

Morphometric features
The edge of TJBH lies at a depth of~5.0 mbsl, which is the depth of Chetumal Bay in that location ( Figure 2A). Georeferenced raster images collected from sidescan sonar recordings ( Figure 1B) and drone surveying ( Figure 2B) identified a nearly circular shape for the mouth of the blue hole. A series of discontinuous terraces were detected at water depths<10 mbsl, near the eastern and northwestern walls of the blue hole and before developing steep and almost vertical slopes. The southwestern wall of the blue hole was the steepest, with a 10:1 slope ratio ( Figure 2C).
The major axis of the TJBH at~10 mbsl was estimated to bẽ 151.8 m and was oriented~10.76°clockwise from the north, similar to the main direction of geological faults in the region. The minor axis, perpendicular to the major axis, was~124.4 m across. The surface area covered by the blue hole is~13,690 m 2 (1.4 ha) with an external perimeter of~492.7 m. The structure of the blue hole can be described as a conic shape that is north-skewed with a volume of~1 million m 3 . The maximum water depth of the blue hole exceeds 270 mbsl, with the deeper areas located mostly scattered on the northern side ( Figures 1B, 2B). However, these areas are situated close to a predominantly north-south elongated area inside the blue hole that was of~100 m 2 across and where consistent depths >240 mbsl were found.
Based on the morphometric features of perimeter and area, the Gravelius coefficient (K C ) or circularity index was equal to 1.179. The elongation ratio (R) between the length of the axes was 1.22. These values reflect the round-ellipsoidal surface shape of the blue hole. A scoured seabed area surrounding the blue hole was identified through aerial drone surveying ( Figure 2B). This scour pattern at~4.5 mbsl describes a 'comet tail' shape directed toward the NE, with its origin close to the blue hole and extending up to a bathymetric limit given by a stepped seafloor ( Figure 2B).

Morphological and biological insights
Biofilms, sediments, and irregularly shaped rocky ledges of 2-3 m were identified covering the internal walls of the TJBH (Figure 3). In the water column, between 15 and 20 mbsl, a layer formed by detrital particles was noted with a general diffuse light and decreased visibility. The detritus was not uniformly distributed inside the blue hole and visibility across sections at the same water depth differed. In addition to the diffusion caused by the detrital layer, light penetration was minimal and artificial light was used.
The walls of the blue hole were covered by a fine-grained layer (presumably gypsum) and coated by brown algal biofilms, which easily disaggregated under agitation or pressure. Over the southwestern wall, exposed white limestone areas several meters across and free of biofilms were found. These were most likely due to material slides or wall collapses as angular rock faces were noticed ( Figure 3C). No caves, large fissures or crevices were identified.
Mucoid filaments of~1 m length were observed floating inside the blue hole or attached to its walls ( Figures 3D-F). These filaments had similar features to those reported in blue holes with anoxic water masses (Iwanowicz et al., 2021). Benthic fauna such as worms of~0.01-0.02 m and exposed calcareous structures of barnacles were also detected directly lying over the algal biofilms ( Figure 3G). These organisms' structures coated areas along the first 5-10 m of the blue hole walls. Pelagic organisms were not observed during the dives.
The bottom sediment surrounding the TJBH, i.e., the Chetumal Bay seabed ( Figure 3C), consists of unconsolidated sand and fine sediment normally resuspended by wind mixing processes. Turbid waters in the bay reduced visibility to between 1 to 2 m. Nevertheless, 0.5 m before reaching the seabed, a clear differentiation between the waters of Chetumal Bay and that of the blue hole was noticed, improving the visibility (5-6 m) and coinciding with the pycnocline location.

Hydrogeochemical characteristics of the water inside the blue hole
Significant variations in salinity and temperature were identified inside the blue hole through vertical CTD profiling ( Figure 4). Profiles measured 500 m to the east, west, north, and south of the blue hole are also compared in Figure 4C. Salinity increased with water depth, ranging from 13.99 PSU at the water surface to 22.65 PSU at 192 mbsl. Simultaneously, the water temperature decreased from 30.0 to 24.7°C. The identification of boundaries (clines) between water layers with different physical properties was clearly revealed by the vertical profiles and their gradients ( Figure 4D). The maximum temperature and salinity gradients (i.e., the thermocline and halocline) depicted a variation of 4.9°C/m and 9.8 PSU/m, respectively, and both were found at 4.5 mbsl and close to the depth of the Chetumal Bay seabed. These depths also coincided with the change in visibility observed during dive operations. These gradients described a boundary that separates the water properties inside the blue hole from those from Chetumal Bay estuary.
Thermoclines and haloclines were also noticed at 99.3 and 136.9 mbsl, where the gradients were 0.03-0.04°C/m and 0.11-0.19 PSU/ m. Despite not being of the same magnitude as those observed close to the entrance of the blue hole, these gradients are related to an almost continuous rate of change of the physical properties as a function of the depth, with the temperature decreasing while salinity increased ( Figure 4C). In contrast, water layers at depths of 4.5-80.9 and 139.5-185.7 mbsl had average gradients of<0.003°C/ m and<0.005 PSU/m, reflecting homogeneous conditions with almost no changes in temperature and salinity. These gradients are below the accuracy error provided by the CTD instrument (i.e. ± 0.01°C and ±0.014 PSU).
Surrounding the TJBH, vertical temperature profiles depicted similar features with depth ( Figure 4C). The mean temperature was 29.8°C ± 0.53, including measurements at the blue hole center. However, changes in the temperature profile below 3.7 mbsl are similar to the profile collected east of the blue hole, depicting a continuous temperature decrease greater than 1.5°C. In terms of salinity, the northern profile had similar values to those at the blue hole center, with 13.83 PSU ± 0.16. The western, southern and eastern profiles recorded salinity values of 12.6-13.25 PSU. However, changes in the eastern profile below 3.7 mbsl, where the salinity increases (> 1.2 PSU), are analogous to those at the blue hole center ( Figure 4C). Therefore, the physical properties of  Frontiers in Marine Science frontiersin.org temperature and salinity in the northeastern surrounding area seem to be more strongly influenced by the blue hole conditions and possibly related to the orientation of the scour pattern shown in Figure 2B. Fluorescence measurements, herein considered as a proxy for fluorescent organic matter inside the water column, reached their maximum in the surface layer at 1.0 m with 0.163 µg/L. These values decreased continuously with depth ( Figure 4C), with values >0.12 µg/L observed in the upper 5 m of the water column above the limit provided by the thermocline and halocline. Below that depth, fluorescence showed homogeneous conditions down to the bottom of the measured profile ( Figure 4D).
The vertical distribution and concentration of additional hydrogeochemical parameters at water depths of 2, 20, 50, 80, and 110 m are described in Figure 5A and Table 1. Total dissolved solids (TDS) increased from 15.77 g/L at the surface to 20.30 g/L at 110 m. In contrast, the concentration of dissolved oxygen (DO) was 7.22 mg/L at the surface. The DO concentration inside the blue hole rapidly decreased with depth, reaching 3.54 mg/L at 20 m depth and 1.93 mg/L at 110 m depth. The same vertical distribution was observed with pH, where the highest value of 8.13 occurred in surface water, decreasing with depth to 6.94 within the anoxic layer at 110 m. In addition, oxidation-reduction potential (ORP), which describes the redox state of the water, was positive in the surface water above the blue hole mouth (216.4 mV) and negative inside the blue hole, with values from -49 mV at 20 m depth to -137 mV at 110 m depth.
Regarding nutrients, the highest concentrations of nitrates NO 3-(0.014 mg/L) and ammonium NH 4+ (0.050 mg/L) were recorded in the surface water above the blue hole, while almost no nitrite NO 2was detected. Inside the TJBH, NH 4+ was the dominant form of dissolved inorganic nitrogen (>97%), where the highest value of 0.037 mg/L was recorded at 110 m depth. NO 3presented a concentration of 0.001 mg/L throughout the water column, and NO 2was undetectable. In contrast, the concentration of orthophosphate, PO 4 3was low in the surface water at 0.003 mg/ L but increased with depth within the blue hole up to 0.010 mg/L at 80 m depth and then decreased further to 0.008 mg/L at 110 m depth. The opposite occurred with silicates, SiO 3 -2 , where the highest concentration of 6.03 mg/L was recorded in the surface water, and then it decreased with depth to 4.57 mg/L at 110 m (Table 1). Chlorophyll-a presented a multimodal behavior, with a high concentration at 20 m depth (2.62 µg/L ± 0.22), which decreased at 50 m (1.90 µg/L ± 0.11) and reached its maximum value at 110 m depth (2.73 µg/L ± 0.29). All anions and cations evaluated presented lower concentrations in the surface water but increased inside the blue hole with respect to depth, and the same was observed with alkalinity and hardness ( Figure 5A, Table 1). The Piper diagram ( Figure 5B) categorizes these parameters and indicates that the water inside the blue hole corresponds to the sodium-chlorinated type, characterized by the primary influence of seawater, with SO 4 2and Mg +2 enrichment corresponding to depth.
Cluster analysis of the evaluated physicochemical variables revealed three layers of stratification within the TJBH ( Figure 6A, Table 1). The first corresponds to a hypoxic layer that lies primarily between water depths of 5 m and 20 m, in which an increase in salinity and a decrease in temperature, dissolved oxygen and pH were recorded. The second is the chemocline that occurs between depths of 50 m and 80 m and constitutes a mixing zone due to the difference between the upper hypoxic layer and the lower nearly anoxic layer, as particularly seen in the measured concentrations of ammonium, orthophosphates, sulfates, sodium, and chlorophyll-a. The third is an anoxic layer that extends to a depth of 110 m,  characterized by low concentrations of DO and phosphates and an increase in ammonium and chlorophyll-a (Table 1). In addition to the stratification given by the physicochemical variables, stiff diagrams of ionic composition (Figures 6B, C) were developed to represent the equivalent concentrations of cations (Ca 2+ , Mg 2+ , Na + , K + ) and anions (HCO 3-, SO 4 2-, Cl -, NO 3-, CO 3 2-) of water samples. An increase in the concentrations of sulfate and calcium ions with respect to depth in TJBH was noticed, which could be due to the dissolution of gypsum that results in high levels of both ions ( Figure 6C). This result coincides with reported dissolved and mineralized gypsum contents in groundwater (almost saturated in sulfates) in the southeast portion of the Yucatań Peninsula and northeastern Belize (Perry et al., 2002;Perry et al., 2019;Perry et al., 2021). Increasing ionic concentrations, for Na + , Cl -, as a function of the water depth were also noticed, coinciding with the increase on salinity. The dissolution capacity is confirmed in this study with the values given by the saturation index (SI), which indicate that the water is undersaturated with respect to gypsum/anhydrite (Table 1).
The TJBH is also characterized by high sulfate concentration (52.65-56.03 meq/L) and high ratios of sulphate and chlorine, 100 x SO 4 -2 /Cl -(16.18-17.90). The increase in salinity with depth in the blue hole would indicate that the primary water source is seawater, but vertical and horizontal water fluxes should be measured. This water source is confirmed by the vertical salinity profile ( Figure 4C) but also by the values obtained from the ratios 100 x SO 4 -2 /Cland Na + /Cl - (Table 1), close to those for seawater of 10.3 and 1.17, respectively (Perry et al., 2002).

Discussion and concluding remarks
The research resulted in the discovery of the second known deepest blue hole formation worldwide, with maximum recorded water depths of 274.4 mbsl. The development of an exceptionally deep blue hole inside the shallow tropical Chetumal Bay is noteworthy and likely associated with glacio-eustatic changes (Siddall et al., 2003;Raymo et al., 2011;Rohling et al., 2021) in the  (Smart et al., 2006). The obtained high concentrations and ratios of sulphate and chlorine are representative of the Evaporite Region, which is geomorphologically distinctive due to restricted permeability (Perry et al., 2002). This region extends from the southeastern Yucatań Peninsula to the border with Belize, where the TJBH is located. The results represent the first research on the TJBH that, without being exhaustive, allowed for the exploration of its geomorphometric and physicochemical features. Stratification based on temperature and chemical characteristics over the water column of the TJBH was found to be similar to that described for the Sansha Yongle Blue Hole (SYBH) in China (Bi et al., 2018;Xie et al., 2019) or the Amberjack Blue Hole in the Gulf of Mexico (Patin et al., 2021); i.e. decreasing temperature and increasing salinity with depth, as well as the development of an anoxic deep water layer (at depths >80 mbsl). It is expected that each submerged blue hole might develop particular features along the water column according to surface surrounding waters and their existing or absent groundwater connections. For instance, the TJBH surface salinity values are representative of the brackish conditions of the Chetumal Bay estuary, while salinity values at deep layers inside the TJBH showed that the main water source is seawater. Nevertheless, the anoxic conditions in deep layers of the TJBH reflect the lack of 1) wind-driven mixing, 2) photosynthetic processes and 3) groundwater circulation. This may led to speculate that saline water in deep layers of the TJBH could represent marine inundations during the rising Holocene sea, as described for the Great Blue Hole of Belize (Gischler et al., 2013). However, the origin and geological evolution of the TJBH deserves further investigation.
Bi et al. (2018) described that the distribution of biogeochemical elements could be controlled not only by physical properties such as temperature and salinity but also by biological and geochemical processes (e.g., photosynthesis, respiration, redox reactions). This could be the case for the distribution of DO in the TJBH, whose high concentration in the surface water layer might occur due to the dissolution of atmospheric oxygen and photosynthetic activity favored by wind-mixing processes, whereas the decrease in DO in the lower water layers would probably occur by predominant organic matter decomposition and respiration processes. Simultaneously, organic carbon oxidation processes and organic matter mineralization by heterotrophic bacteria favoring water acidification (Xie et al., 2019) could be responsible for the pH decrease with depth in the TJBH (Table 1).
The maximum value of chlorophyll-a and TDS was recorded inside the TJBH below 110 m depth, where the lowest value of DO was obtained. According to Bi et al. (2018), this could be caused by an enrichment of particles containing chlorophyll-a or the production of chlorophyll-a by cyanobacteria under anaerobic and aphotic conditions. In the TJBH, detrital particles containing chlorophyll-a could be retained at the thermo-and haloclines developed around the 110 m depth (Figure 4), thus increasing the chlorophyll-a values or the bacterial production processes (microbial layer), as similarly described in Björnerås et al. (2020) for the blue holes in The Bahamas. Other authors also describe that Dendrogram describing the clustering of water samples with similar chemical features as a function of the water depth (A). Stiff diagrams of ionic composition at 500 m west and east from the center (B), as well as 2, 20, 50, 80 and 110 m depth at the center of the blue hole (C).
Alcé rreca-Huerta et al. 10.3389/fmars.2023.1141160 Frontiers in Marine Science frontiersin.org bacterial and primary production significantly contributes to the high variability in biogeochemical characteristics and food webs inside blue holes (Macalady et al., 2008;Gonzalez et al., 2011;Xie et al., 2019;Björnerås et al., 2020). Low or undetectable concentrations of NO 3and NO 2inside the TJBH could be due to denitrification and anammox processes (reduction of NO 2to N 2 ), similar to the nitrogen cycle described by Xie et al. (2019) for the SYBH. Additionally, the levels of NH 4+ , representing more than 97% of the DIN (dissolved inorganic nitrogen), result from ammonification processes in the blue hole. Further studies should be conducted, including the analysis of the microbial diversity and hydrogen sulfide measurements, since the reduction of sulfates is among the main metabolic pathways in karstic aquifers. Suspended solids, turbidity, photosynthetically active radiation (PAR), biological components, and sediment analysis are also suggested to allow the identification of the biogeochemical processes concerning the stratification of blue holes (Xie et al., 2019;Björnerås et al., 2020;He et al., 2020;Patin et al., 2021). Isotopes could be explored to identify the age and further features of the water (Bottrell et al., 1991) and the sediment deposits  inside the TJBH. XRD data could also provide a more detailed analysis of the mineralogical composition of the karst structure. Measurements from remotely operated vehicles (ROVs) or multibeam systems are also required to obtain more accurate 3D topographic data and to identify possible interactions between the morphometry, the hydrobiochemical and geological features of the region (Li et al., 2018;Gonzalez et al., 2020). Moreover, monitoring of the temporal changes in the physicochemical and biological features inside the TJBH is needed to understand the possible variability in biogeochemical processes and to explore the water exchange with the surrounding estuary of Chetumal Bay.
The close location of the TJBH to the coasts of Mexico and Belize would inevitably imply an increase of local and foreign visitors, the development of economic activities, and very probably of environmental pressure. Regulatory measures might be required although the blue hole forms part of the state reserve Sanctuary of Manatee. Currently, the local population is not aware of the existence of the TJBH. Therefore, the authors encourage the scientific community to explore, monitor, and broaden the research of the TJBH to settle an adequate basis for responsible social appropriation in the near future.

Data availability statement
The datasets presented in this article are not readily available because the data belong to a project funded by the authors. Once published, Data eventually will be shared on the institutional data reservoir. Requests to access the datasets should be directed to Dr. LC, lcarrillo@ecosur.mx.

Funding
The first author personally funded fieldwork expenses during the field survey. The APC fees were personally funded by the authors. The fourth author acknowledges the CONACYT Mexican Scholarship Program (Grant to CVU 1111960) for the funding of postgraduate studies.