- 1Marine Biology Section, Department of Biology, University of Copenhagen, Helsingør, Denmark
- 2Biological and Environmental Sciences and Engineering Division, King Abdullah University of Science and Technology, Thuwal, Saudi Arabia
Reef building corals are efficient biological collectors of solar radiation and consist of a thin stratified tissue layer spread over a light scattering calcium carbonate skeleton surface that together construct complex three dimensional (3D) colony structures forming the foundation of coral reefs. They exhibit a vast diversity of structural forms to maximize photosynthesis of their dinoflagellate endosymbionts (Symbiodiniaceae), while simultaneously minimizing photodamage, offer resistance to hydrodynamic stress, reduce attack by predators and increase prey capture and heterotrophic feeding. The symbiosis takes place in the presence of dynamic gradients of light, temperature and chemical species that are affected by the interaction of incident irradiance and water flow with the coral colony. We developed a multiphysics modelling approach to simulate the microscale spatial distribution of light, temperature and O2 in a coral fragment with its morphology determined by 3D scanning techniques. Model results compared well with spatial measurements of light, O2 and temperature under similar flow and light conditions. The model enabled us to infer the effect of coral morphology and light scattering in tissue and skeleton on the internal light environment experienced by the endosymbionts, as well as the combined contribution of light, water flow and ciliary movement on O2 and temperature distributions in the coral.
1 Introduction
Reef building, scleractinic corals construct a complex three dimensional (3D) calcium carbonate skeleton, which is the framework for coral reef ecosystems of considerable biological and socio-economic importance (Costanza et al., 2014). Tropical, scleractinic corals rely on endosymbiotic dinoflagellate algae from the family Symbiodiniaceae, which are hosted in the coral endoderm tissue and excrete photosynthates to the host, providing the majority (up to >90%) of carbon needed for coral animal respiration along with O2 (Muscatine, 1973). Symbiont photosynthesis can also enhance coral calcification (Chalker, 1981; Goreau et al., 1996). In exchange, the host provides the endosymbionts with a conducive and protected environment and a steady, albeit limited supply of nutrients.
However, the coral-algal symbiosis can rapidly deteriorate under environmental stress (Moberg and Folke, 1999; LaJeunesse et al., 2018). Environmental factors related to ongoing climate change such as ocean acidification (Jackson et al., 2001; van der Zande et al., 2020), deoxygenation (Altieri et al., 2017; Hughes et al., 2020) and warming (Hughes et al., 2017), increasingly result in severe mass bleaching and mortality of corals (Knowlton et al., 2021). Coral bleaching and mortality varies as a consequence of solar irradiation levels (Dunne and Brown, 2001; Mumby et al., 2001), local seawater temperature variation (Teneva et al., 2012), coral colony morphology (Loya et al., 2001; Van Woesik et al., 2012), water flow around the coral (Nakamura et al., 2003; Jimenez et al., 2008), and O2 levels in the surrounding water (Altieri et al., 2017; Johnson et al., 2018; Hughes et al., 2020). Hence, it is important to improve the mechanistic understanding of how the photobiology and physiological activity of the coral holobiont are affected by macroscale colony morphology and microscale variations in tissue/skeleton architecture, and how interactions with the incident solar irradiance and flowing seawater modulate the physico-chemical microenvironment and metabolic activity of corals. This will enable us to better estimate hotspots for coral bleaching, the coral reef’s response to environmental stressors, and the relationship between the host microenvironment and microbiome/endosymbiont populations. Understanding light propagation, heat and mass transfer of O2 in relation to the 3D morphology of a coral colony can also provide insight into the mechanisms governing spatial and temporal variability in coral bleaching and recovery within and between coral species.
Corals are efficient biological collectors of solar radiation (Wangpraseurt et al., 2014a) and consist of a thin (from ~100 μm to several mm) stratified tissue layer spread over a light-scattering skeleton matrix (Davy et al., 2012). They are exposed to a wide range of light environments ranging from super-saturating, broadband (UV to near infrared) solar irradiance in shallow reefs (Veal et al., 2010; Jimenez et al., 2012; Wangpraseurt et al., 2014b) to very low irradiance in the blue-green spectral range at mesophotic depths (Eyal et al., 2015; Tamir et al., 2019). The morphology of the colony and the surrounding reef, along with the optical properties of the host tissue and skeleton, can lead to variable light conditions experienced by different parts of a coral colony and hence the endosymbionts (Kaniewska et al., 2011; Wangpraseurt et al., 2014b).
Corals have developed a vast diversity of structural forms and mechanisms at cellular, polyp and colony level, to maximize photosynthesis of their endosymbionts under spatial and temporal variation in light availability (Anthony and Hoegh-Guldberg, 2003; Todd, 2008; Kaniewska et al., 2011), while simultaneously minimizing photo-damage (Terán et al., 2010; Brodersen et al., 2014) and oxidative stress (Pacherres et al., 2022), and maximizing solute and heat exchange (Jimenez et al., 2011). Consequently, corals can reach high photosynthetic quantum efficiencies close to 0.1 O2 photon-1, approaching theoretical limits, under moderate irradiance levels (Dubinsky et al., 1984; Brodersen et al., 2014). Structural diversity also helps coral improve their growth rates via alleviating hydrodynamic stress (Veron, 2000), enhance autotrophy by improving mass transfer of gasses, nutrients and metabolites across coral tissue (Falter et al., 2004), improved prey capture for heterotrophic feeding, as well as reducing attack by predators(Smee et al., 2010; Lenihan et al., 2015).
Light absorption by corals not only drives symbiont photosynthesis in hospite, but can also affect the thermal microenvironment via local heating (Jimenez et al., 2008), as most of the absorbed light energy is dissipated as heat (Brodersen et al., 2014). Light is mainly absorbed in the coral tissue and, to some extent, in the skeleton by endolithic algae and photosynthetic bacteria (Kühl et al., 2008; Ricci et al., 2023). However, both coral morphology (Kaniewska et al., 2011; Kaniewska and Sampayo, 2022; Kramer et al., 2022) and the inherent scattering properties of coral skeleton and tissue modulate light propagation and thus the light exposure of the microalgal symbionts (Enríquez et al., 2005; Lyndby et al., 2016; Wangpraseurt et al., 2016; Bollati et al., 2022).
The morphology and the degree of branching in coral colonies affect not only the light field but also the water flow around the colonies. The flow determines the thickness and shape of boundary layers for momentum transfer (momentum boundary layer, MBL, i.e., water velocity change near the coral surface), mass transfer (diffusive boundary layer, DBL, i.e., concentration change in the diffusion-dominated region adjacent to the coral tissue) (Chan et al., 2016), as well as heat transfer (thermal boundary layer, TBL, referring to temperature changes) (Jimenez et al., 2008; Ong et al., 2019). The coral consumes part of the O2 produced by the symbionts under light (Al-Horani et al., 2003), while an eventual surplus of oxygen and heat generated in the tissue is transported into the surrounding water column and the coral skeleton. By introducing a resistance to mass and heat transfer, boundary layers decrease the mass and heat exchange rates between the bulk water and the coral, to the extent in which solute diffusion (Shashar et al., 1993) and heat conduction (Ong et al., 2012a; Ong et al., 2012b; Ong et al., 2017; Ong et al., 2019) can become bottlenecks. These transfer resistances can be alleviated to some extent by actively decreasing boundary layer thickness through an intensified flow, as achieved by tissue and polyp level changes like contraction-expansion and ciliary movement that can enhance mass (Shapiro et al., 2014; Pacherres et al., 2020; Pacherres et al., 2022) and heat (Ong et al., 2017) transfer, along with affecting light penetration (Wangpraseurt et al., 2014a; Wangpraseurt et al., 2017). The resulting spatial and temporal variations in temperature distribution and concentration of chemical species (e.g., dissolved oxygen, inorganic carbon and pH) affect the coral microenvironment and can form distinct microhabitats within a given coral colony. Such fine scale ecological niche heterogeneity may affect the composition of both endosymbionts and microbiomes across the different compartments of the coral holobiont, but the microscale fitness landscape of corals remains largely unexplored, in part due to technical challenges (Hughes et al., 2022).
Numerical modelling is a powerful way to integrate the physical, chemical and biological complexity (at several spatial and temporal scales) into a systematic framework, with the aim to describe, understand and ultimately predict coral responses to a changing environment. Modelling the growth of coral colonies in response to environmental parameters has been pioneered by the work of Kaandorp and coworkers (Chindapol et al., 2013; Kaandorp, 2013). Several studies have also modeled the interaction of light and fluid flow with colony morphology and the resulting surface temperature (Ong et al., 2012a; Ong et al., 2012b; Ong et al., 2017; Ong et al., 2018; Ong et al., 2019) and mass transfer at the coral surface (Chang et al., 2014; Shapiro et al., 2014). However, these models have not included radiative transfer in the tissue-skeleton matrix, or simulations of internal gradient of light, temperature and chemical parameters across the tissue and skeleton compartments.
Empirical and theoretical models have shown how skeleton-dependent scattering can enhance the local light field and absorption of Symbiodiniaceae in hospite (Enríquez et al., 2005; Swain et al., 2016). Ray tracing (Ong et al., 2018) and probabilistic Monte Carlo (MC) modeling techniques for light propagation (Terán et al., 2010; Wangpraseurt et al., 2016) assuming a simple two-layer tissue-skeleton geometry have also been employed. Recent advances in estimating inherent optical properties of different species of living corals (Jacques et al., 2019; Spicer et al., 2019; Wangpraseurt et al., 2019) now enable more accurate simulation of internal light fields in specific corals.
We recently developed a numerical model of a stratified coral tissue on top of skeleton to link internal light distribution to light absorption, radiative heat dissipation, heat transfer and mass transfer of photosynthetically produced O2 (Taylor Parkins et al., 2021). The model couples a Monte Carlo (MC) simulation of light propagation with numerical modelling of heat production and metabolism inside the coral to simulate irradiance, temperature and O2 microprofiles at small scale in simplified, schematic geometries. In the present study, we have expanded this multiphysics modelling approach to simulate microscale distribution of light, temperature and O2 in and around branched coral fragments (cm scale) with complex, natural morphology, as determined by 3D scanning techniques. We validate the model by comparing simulated microenvironmental parameters over the coral topography with the corresponding microscale measurements on the same fragment under similar irradiance and laminar flow conditions. This enabled us to describe the effect of coral morphology on the internal light environment and the combined effects of light, water flow and ciliary movement on the O2 and temperature distribution in the coral.
2 Materials and methods
The methodological approach involved an experimental and a theoretical part. First, microsensor measurements of light, O2 and temperature were performed on a fragment of the branched coral Stylophora pistillata (Figures 1A, B) in a flow chamber under defined irradiance and flow conditions. The fragment was subsequently 3D scanned and meshed (Figures 1C, D), followed by multiphysics 3D modeling of the radiative, heat and mass transfer of the meshed structure (Figure 1F) and subsequent comparison of simulated distributions of light, O2 and temperature with the experimental measurements (Figure 1E).
Figure 1 Overview of experimental and modeling approach. Experimental flow chamber set-up for microsensor measurements (A) on a coral fragment of Stylophora pistilla (B, C). Tetrahedral mesh based on the 3D scan of the coral (D), used for multi-physics simulation. Comparison of microsensor measurements and simulations (here illustrated with measured and simulated O2 concentration profiles) in particular areas of the coral (E). Simulated surface values of normalized scalar irradiance, O2 concentration and temperature at the coral tissue-water interface (F).
2.1 Coral specimens and husbandry
Colonies of the branched coral Stylophora pistillata were obtained from a sustainable, commercial provider (Dejong Marine Life, Netherlands). Colonies were fragmented into smaller nubbins that were mounted in small plastic caps using AF Gel Fix (https://aquaforest.eu) and kept in dedicated aquaria supplied with continuous flow of artificial reef water at 25-26°C, salinity of 36 g kg-1, and moderate levels of downwelling photon irradiance (150 - 200 µmol photons m-2 s-1; 400–700 nm) provided over a 12:12 hour dark-light cycle by a programmable aquarium lamp (MAXSPECT R5X 150W; Decocean Aquarium).
2.2 Experimental setup
A small fragment of S. pistillata (Figure 1B) was placed in a custom-made black, acrylic flow chamber, and measurements were done with two different orientations of the fragment relative to the laminar flow. The samples were placed for at least 1 hour in the chamber prior to microsensor measurements to ensure steady-state conditions. The coral fragment was continuously flushed with aerated seawater at 26°C and salinity of 35 g kg-1. An average flow velocity of ~1 mm s-1 was maintained by a water pump connected to the flow chamber and submerged in the thermostated aquarium reservoir.
The sample was illuminated with a constant downwelling photon irradiance of 1300 µmol photons m-2 s-1 from a fibre-optic tungsten halogen lamp equipped with a heat filter and a collimating lens (KL-2500LCD, Schott GmbH, Germany), positioned vertically above the flow chamber (Figure 1A). The downwelling photon irradiance (Ed in units of µmol photons m-2 s-1) of photosynthetically active radiation (PAR; 400–700 nm) was measured with a calibrated photon irradiance meter (ULM-500, Walz GmbH, Germany) equipped with a planar cosine collector (LI-192S, LiCor, USA) positioned in the light path approximately at the same distance as the coral. We note that the incident light conditions in the reef would most likely shift between collimated and diffuse light fields, and would only reach such high irradiance e.g. on reef flats during mid-day, with clear skies and calm waters (Jimenez et al., 2012; Tagliafico et al., 2022). Such high incident light conditions were chosen for simplicity and to be able to measure temperature increments in the coral fragment above the detection limit of the sensor. Effects of diffuse light field on coral light harvesting and energy budgets (e.g., Lichtenberg et al., 2017) can, however, be incorporated both in the simulations and experimental set-up for future studies.
2.3 Microsensor measurements
Measurements were conducted on the connective tissue, i.e., coenosarc, between individual polyps, which had a more even topography and less contractile tissue as compared to polyp tissue. For measurements, the microsensors were mounted on a motorized micromanipulator (MU-1, PyroScience, GmbH), which was interfaced to a PC and controlled by dedicated data acquisition software (ProFix, Pyros-Science GmbH). All measurements were made at an angle of 45° relative to the vertically incident light beam to avoid shading. Scalar irradiance measurements were done by placing the sensor tip at the coral tissue surface. Positioning of microsensor tips relative to the coral tissue surface was monitored visually by observation under a dissection scope. The O2 and temperature profiles were measured from the water column into the coenosarc in vertical steps of 100 µm, as described previously (Jimenez et al., 2008). The profiles (two profiles each) were made at the coenosarc tissue in three different regions of the coral fragment as shown in Figure 1D.
Light measurements: Spectral scalar irradiance was measured using a fibre-optic scalar irradiance microprobe with a spherical tip diameter of ~100 µm (Rickelt et al., 2016) connected to a fibre-optic spectrometer (USB 2000+, Ocean Optics, USA). The incident downwelling irradiance was measured at the same height as the coral surface by placing the light sensor in a black (non-reflective) light well under the vertically incident light in the flow chamber. All spectral scalar irradiance measurements at the coral tissue surface were normalized to the incident spectral irradiance. The results are represented as mean ± standard deviation, averaged over 12 measurement points taken close to each other in a given region of interest (as indicated by squares in Figure 2C) and for the 2 different orientations of the fragment with respect to the flow.
Figure 2 Comparison between measured and simulated scalar irradiance (636 nm; ratio of scalar irradiance to incident irradiance). (A) Simulated scalar irradiance distribution over the coral tissue-water interface; (B) Cross-section of simulated scalar irradiance depth distribution along the middle section of the fragment; (C) Measured scalar irradiance in similar areas, from which simulated values were extracted. Optical properties (OP1) used for simulation are listed in Table 1. The values shown here represent means ± standard deviation, averaged over 12 different points in each area of interest as indicated by the squares.
Oxygen measurements: Oxygen concentrations were measured with Clark-type microelectrodes (tip diameter approx. 25 µm; OX-25, Unisense A/S, Aarhus, Denmark)(Revsbech, 1989). The microsensor was connected to a pA-meter (Oxymeter, Unisense A/S) interfaced to a PC. The O2 microsensors were linearly calibrated from sensor signal measurements in aerated and anoxic seawater at experimental temperature and salinity.
Temperature measurements: Thermocouple microsensors (tip diameter approx. 50 mm; T50, Unisense A/S) were connected to a thermocouple meter (Unisense A/S). A high precision thermometer (Testo 110, Testo AG, Germany) was used to linearly calibrate the temperature microsensor signal from readings in seawater at different temperatures. The temperature profiles are represented as the temperature increase relative to the water temperature (in the free flowing region) directly above the measurement point.
2.4 3D scanning and optical coherence tomography imaging
We used a PC-controlled, structured light 3D scanner (Einscan-SP, Shining 3D) to generate non-texture 3D scans of coral fragments in air. The resolution of the 3D scanner is ~50 µm. A sample was placed on the turn table of the calibrated scanner. Data acquisition and processing was controlled by the system software (EXScan S_V3; Shining 3D). The brightness setting of the scanner was adjusted to avoid saturation of the camera signal. A 360° scan of the fragment was made under different orientations, to capture the complete morphology of the fragment avoiding missing faces. The different scans were aligned using the “align by feature” mode in the scanner software. A watertight mesh of the model was exported as an STL file. The fragment was scanned with and without the tissue layer to estimate an average tissue thickness, by dividing the volume and surface area difference between the two scans. The tissue was removed by immersing the fragment in 30% hydrogen peroxide solution overnight. This procedure led to a tissue thickness estimation of ~900 µm. However, this could be an overestimate as the reconstruction software fills in holes in the scan by extrapolation, in places where it could not scan the object correctly (e.g., due to complex shape and highly absorbing or reflective surfaces), especially for the scan with tissue. Another reason for overestimation of the tissue thickness, could be due to the immersion of coral fragment in 30% H2O2 overnight, which might have led to some skeleton dissolution (Holcomb et al., 2015).
For an alternative determination of the coral tissue thickness, we used a 930 nm spectral domain optical coherence tomography (OCT) system (Ganymed II, Thorlabs, Germany) equipped with an objective lens with an effective focal length of 18 mm and a working distance of 7.5 mm (LSM02-BB; Thorlabs GmbH, Dachau, Germany) for imaging of corals immersed in seawater with a maximal axial and lateral resolution in water of 5.8 μm and 8 μm, respectively (Wangpraseurt et al., 2017). Two-dimensional OCT B-scans were acquired at a fixed pixel size of 581 x 1024. The scans were used to estimate the tissue thickness (Figure S1B), however, due to shadowing effects and uncertainties in the exact refractive index of the coral tissue (here assumed to be 1.33), the OCT measurements underestimated tissue thickness (~400 µm). For the simulations, we therefore selected a tissue thickness of 650 μm, which was an average of the values from the OCT and 3D scans.
2.5 Numerical modeling
A three-dimensional mathematical model was constructed with the aim of simulating the spatial distribution of O2 and temperature around and within the coral, as influenced by the light transport and by the water flow. Model predictions were compared to the measured profiles of light, O2 and temperature on the coral fragment.
For the simulations, the scanned 3D coral geometry was assumed to consist of a tissue layer with an uniform thickness created on top of the skeleton (Figures S1A, B), which was placed in a rectangular flow chamber (Figure S1C). Different sub-layers within the coral tissue (Taylor Parkins et al., 2021), with varying material properties and supporting different chemical reactions, were represented indirectly as a function of the distance from the tissue-water surface. The model first computes the light field within and around the coral (Figure S1D), then the laminar flow of water around the coral in the flow chamber. Subsequently, mass and heat balances allow calculation of the radiation-driven O2 and temperature distribution in the stratified coral domain and the water column (Figures S1E, F). All the model parameters are listed in the supplementary information (Table 1).
2.5.1 Model geometry
The model geometry consisted of a solid domain, the 3D scan of the coral fragment, enclosed by the liquid in a box matching the dimensions of the experimental flow chamber (Figure S1C). This geometry was built in COMSOL Multiphysics (v6, COMSOL Inc., Burlington, MA) based on the STL file provided by the coral 3D scanning. Two different orientations of the fragment with respect to the flow (Figure S1) were simulated.
An essential part of the model includes the effect of coral tissue sub-layers with various properties, thus assigning different material properties and reactions as a function of depth within the tissue, as done in (Taylor Parkins et al., 2021) albeit in a much simpler geometry. However, an explicit partitioning of the large 3D coral domain in several very thin subdomains on the skeleton surface proved to be very difficult computationally, especially regarding the accurate meshing of these thin layers. We therefore adopted the solution of representing these layers implicitly, by assigning different material properties, based on the distance (d) from the coral/water surface, as shown in Figure S1A (i.e. mucus layer for d from 0 to 80 µm, epidermis for d from 80 to 180 µm and so on). The perpendicular distance (d) within the coral from the coral/water surface was computed in COMSOL (the wall-distance interface) by a modified Eikonal equation (Fares and Schröder, 2002). Consequently, the optical, mass and heat transfer properties characteristic for different tissue sub-layers were defined as a function of this distance, d, as explained in further detail in the following sections.
2.5.2 Radiative transfer
The scalar irradiance at different positions in the coral and water was determined by calculating the radiative transfer (Figure S1D) using the 3D Monte Carlo (MC) approach implemented in the free software ValoMC (Leino et al., 2019), as described in our previous work (Taylor Parkins et al., 2021). The simulation of photon transport was performed by launching photons with a wavelength of 636 nm. This wavelength is within the absorption band for Chl c and Chl a of the coral microalgal symbionts. The present model ignores fluorescence and takes only scattering and absorption into account. Various tissue layers and skeleton optical properties (OP1, Table 1) were assigned as a function of the calculated wall distance d, i.e. tissue optical properties for d< 650 µm and skeleton optical properties for d > 650 µm
The simulated point cloud of scalar irradiance values from ValoMC were imported into COMSOL and mapped over the 3D geometry by solving Poisson’s diffusion equation in weak form, in order to smoothen the variations in the scalar irradiance data inherent due to the stochastic nature of the MC simulation:
where Is is the smoothed normalised scalar irradiance, I0 is the initial normalised scalar irradiance, h is the local mesh element size, and fs is a smoothing factor (0.1; 0< fs< 1). A zero light flux condition was applied to all boundaries.
The simulated scalar irradiance at different positions in the geometry was normalized to the incident irradiance. The results are represented as mean ± standard deviation, averaged over 12 different points, sampled close to each other in each investigated area, as indicated by the squares in Figure 2A.
2.5.3 Fluid flow
Stationary incompressible Navier-Stokes equations for laminar flow (Reynolds number of ~9) were used to simulate the water flow around the coral fragment using COMSOL:
where u is the velocity vector, p the pressure, the dynamic viscosity and the density of water. The water inflow had an average velocity of 1 mm s-1, while a zero gauge pressure was set in the outflow. Top, bottom and lateral walls of the flow cell were no-slip (zero-velocity). Two cases were assumed for flow at the coral surface: zero-velocity (without ciliary movement) and a set velocity (with ciliary movement). Ciliary movement as proposed in other studies (Shapiro et al., 2014; Pacherres et al., 2020) was included as cilia-induced currents by assuming an oscillating horizontal velocity component at the coral surface as (Pacherres et al., 2022):
with and y being the distances from surface, the maximum ciliate beating velocity here set to 150 µm s-1 as measured experimentally (Shapiro et al., 2014), and the characteristic length scale of the vortices (assumed 1200 µm, close to the average calyx size of the coral). We note that our representation of ciliary movement is a rough estimate. The length of the vortices might vary depending on the location on the fragment, flow conditions, ambient water temperature etc. To model this in more detail, we need more fundamental insight to e.g. cilia coverage, dimensions and beating frequency and motility over coral tissues and across different species. It’s in scope of future work to model this in more detail, as more experimental data become available.
The simulated flow profile around the coral was similar with and without the base onto which the coral was glued (Figure S2). Hence, in order to simplify the geometry and the computational burden, all the simulations presented here were made without the base.
2.5.4 Oxygen transport and reactions
The dissolved oxygen concentration, cO2,in the coral and surrounding water was computed from a stationary material balance, written in a general form as eq. (6):
Diffusive transport was applicable in both domains with diffusion coefficient DO2 taking the value DO2,w for dissolved O2 in water, but reduced values in the coral (0.5·DO2,w in the tissue and 0.01·DO2,w in the skeleton). In the water domain there were no reactions involving oxygen (RO2=0), while in the coral tissue/skeleton there was no convective transport (u=0). A constant O2 concentration c0,O2 was imposed in the water inflow and the classical convection-only condition was assumed in the water outlet, with the rest of flow cell walls insulated (no flux of O2). The flux and concentration continuity were assumed on the coral/water interface.
The net photosynthetic O2 production rates in oral and aboral gastrodermis (gas) were calculated as a function of scalar irradiance Ip using the average light absorption coefficient of the coral tissue :
Here, the quantum efficiency of symbiont photosynthesis (mol O2 produced per mol photon) is represented as a function decreasing with increasing scalar irradiance (Figure S3; Brodersen et al., 2014). The rate is furthermore limited by the ratio of photosynthesis to respiration rate, PR=3.5 (Cooper et al., 2011). fgas is a switch function depending on the distance d from the coral/water interface, which allows defining RO2 only in the gastrodermis, i.e., = 1 if the distance is 260 – 360 µm for oral, 460 – 560 µm for aboral gastrodermis layer, and = 0 elsewhere.
The respiration rate in the epidermis layer was calculated as function of O2 concentration cO2:
with the coral host maximum respiration rate and KO2 the half-saturation coefficient for O2 limitation for respiration. Again, the switch function activates this rate ( = 1) only at a distance d between 80 and 160 µm from the coral surface.
The O2 consumption in the coral skeleton by endolithic bacteria follows a similar rate:
with Rcp,skel the maximum O2 consumption rate and the same half-saturation coefficient KO2. switches on this rate for d > 650 µm and off elsewhere.
2.5.5 Heat generation and transport
The spatial temperature T distribution in coral and surrounding water was computed from the general heat balance equation including conduction, convection and source terms:
where k is the thermal conductivity of the specific layer, Cp and ρ are the specific heat capacity and density of water, Q is the heat source and u is the vector of water velocity. In particular, heat conduction and generation terms were considered in the coral fragment, while conduction and convection were included in the water domain (neglecting heat generated from light absorption in water). The heat source, Q, originates from absorbed light only, which is proportional with the scalar irradiance, Ip, and the light absorption coefficient, μa, of the specific layer:
where switch functions correspond to specific layers at different distanced from the coral/water surface. fheat is the fraction of light dissipated as heat in the oral and aboral gastrodermis (Figure S4; Brodersen et al., 2014), where fheat<1 because part of the light energy is used in photosynthesis, while fheat was set to 1 in the other tissue layers and skeleton.
The inflow of water had a temperature of T0, while a zero-temperature-gradient condition was set for the water outflow. Thermal insulation was set at all other flow-cell walls, while heat flux continuity was assumed at the coral/water interface. The thermal properties of the skeleton were taken from Jimenez et al. (Jimenez et al., 2008), while the tissue thermal properties were assigned values similar to human tissue (Hasgall et al., 2022) (due to lack of other measurements).
2.6 Comparison between simulated and measured data
For comparison of measured and simulated scalar irradiance values, scalar irradiance microsensor measurements at 636 nm (normalized to incident irradiance) at the coral tissue/water interface were averaged in different areas of interest over the coral fragment, and were then compared to simulated scalar irradiance values in the same areas of interest (Figure 2). For comparison of measured and simulated O2 and temperature measurements, we compared measured profiles in particular positions over the coral structure with extracted simulated vertical profiles of O2 concentration and temperature over the same regions (Figure S5).
3 Results and discussion
We used a combination of microsensor measurements, 3D scanning and numerical modeling to investigate the influence of coral structure and morphology on the light, O2 and temperature microenvironment in and around a coral fragment (Figure 1). Microsensor measurements are local and it is tedious to map the entire sample and impossible to fully account for spatial heterogeneity and hot spots with such one-dimensional measurement approach, and intr-tissue/skeleton measurements are very complicated, if not impossible. Supplementing high resolution microsensor measurements with simulations of the light, O2 concentration and temperature distribution over the full 3D morphology of actual coral samples can help identifying areas of interest for detailed microsensor measurements and lead to better data interpretation. On the other hand, supporting simulations with measurements is imperative to build a realistic model. Such combination of simulations and microsensor measurements provides a powerful toolset for exploring the coral microenvironment under different flow and light scenarios.
3.1 Light microenvironment
The 3D light simulation generates the distribution of the scalar irradiance within the whole computational domain including water, coral tissue and skeleton. Computed scalar irradiance is displayed at the surface (Figure 2A) and within the coral fragment (Figure 2B). Both the simulations and the light measurements demonstrated the presence of a heterogeneous light field over and within the coral fragment. The computed scalar irradiance shows a fair match with the measurements (Figure 2C) at the same wavelength (636 nm).
Both simulations and the measurements showed that the surface scalar irradiance values (quantified as the ratio of the scalar irradiance over the incident irradiance; mean ± standard deviation, n=12 in each area of interest) were highest at the center of the fragment (measured: 1.40 ± 0.13; simulated: 1.50 ± 0.04; Figures 2A, C), where the tissue surface was relatively flat. At the tips of the fragment, the relative scalar irradiance values were lower because the surface was more inclined with respect to the incoming light.
The rear (shaded) sides of the fragment exhibited much lower scalar irradiance values around half of the incident irradiance (measured and simulated). The simulations indicate that the highest scalar irradiance was close to the tissue-skeleton interface in the middle of the fragment, reaching up to ~1.7 times the incident irradiance (Figure 2B). Overall, the simulated scalar irradiance at the coral surface compared well to actual measured values in the same regions of the investigated coral. We note that the light simulations relied on literature values of the inherent optical parameters of coral tissue and skeleton (Jacques et al., 2019) which are complex to estimate as it involves e.g. various tissue sublayers with different refractive indices (e.g. Wangpraseurt et al., 2014a), scattering and absorption properties (Wangpraseurt et al., 2019).
To illustrate how the model responds to differences in optical parameters, we also simulated the light field in another fragment of S. pistillata fragment using optical parameters determined from OCT measurements (Wangpraseurt et al., 2019) (Figures S6, S7). While the absolute values of scalar irradiance enhancement changed between simulations with different optical parameters, the overall light distribution remained more or less identical between the simulations and reflected the pattern in the experimental light measurements. This highlights the fact that we are still lacking a thorough fine scale characterization of inherent optical properties of coral tissue and skeleton.
The light simulation and measurements revealed an interplay between skeleton and tissue optics that may be important in enhancing coral light harvesting. The light scattering from the skeleton leads to light enhancement in the tissue and a distribution of incident light to shaded areas. For a given direction of the incident light, the distribution of light can vary significantly, even across a small region of a coral colony, leading to hotspots and shaded regions. Light enhancement in a coral fragment is due to a combination of surface morphology and scattering in the tissue and skeleton, which enhances and redistributes light within the fragment (Wangpraseurt et al., 2016). Regions of a coral with relatively flat surfaces (angle of incident light relative to the tissue surface close to 0°) will experience less loss of light due to surface reflections arising from refractive index mismatch, as compared to more inclined surface areas in the colony experiencing lower incident irradiance (due to the cosine dependence of Fresnel reflection on the angle of incidence). However, even shaded regions in a coral colony can receive some light due to lateral distribution of light via the tissue and skeleton (Wangpraseurt et al., 2014a; Wangpraseurt et al., 2016; Enríquez et al., 2017).
The light field in corals is affected by the optical properties of tissue and skeleton, as well as tissue thickness and composition (e.g. distribution of symbionts and coral host pigments), and overall colony morphology. Hence the irradiance distribution between two morphologically similar coral fragments can vary, and in a reef environment a coral colony could experience varying external and internal light fields throughout the day, depending on the sun angle. Such heterogeneity in the light microenvironment across tissue and over the coral colony surface might present different optical niches that can drive phenotypic or genotypic diversification of symbionts (Lichtenberg et al., 2016) with different light adaptation and bleaching resistance across the same colony or between colonies with different morphology. It is, however, difficult to account for such spatial heterogeneity with microscale light measurements, especially for intra-tissue measurements that often rely on making a small incision in the coral tissue for probe insertion (Wangpraseurt et al., 2012). We show here that simulation of the spatial light distribution in corals with a known morphology and tissue structure is a powerful supplement to fine scale measurements of coral light fields. Such simulations can predict internal light fields even in the coral skeleton compartment, which is very difficult to measure experimentally and involves invasive micro drilling (Magnusson et al., 2007).
3.2 Oxygen microenvironment
The dissolved O2 distribution in and around a coral colony for a given light field or in darkness is affected by the orientation of the coral with respect to the water flow field. Thus, different parts of the colony will experience changes in the thickness of the MBL and DBL depending on their flow exposure, which will translate into differences in O2 concentration within and at the surface of tissue. The 3D simulations of light-driven O2 production in the coral fragment within the flow-cell were executed under two fragment orientations with respect to the water flow: with the groove between the two branches shaded against the flow by one branch (Figures 3A–D), and with the groove and both branches directly exposed to the flow (Figures 3E–H). The simulations showed higher O2 concentration in parts of the fragment that were shaded from the flow and exhibited a thicker DBL, as compared to more exposed areas (Figures 4A, D). The O2 concentration was high at the skeleton/tissue interface, due to the low diffusivity of the skeleton matrix (Figures 4B, C, E, F).
Figure 3 Computed flow field around the coral fragment in two orientations relative to the flow. (A–D) aligned with the flow; (E–H) across the flow. (A, E) 3D streamlines of water flowing around the coral fragment. b,f: flow velocity magnitude (color map) and streamlines in the planar sections indicated in (A, E) with red dash-dot lines. (C, D, G, H) details of velocity magnitude and streamlines in the neighborhood of the coral surface, with and without ciliary movement (same section planes as in b and f, respectively).
Figure 4 Computed dissolved O2 concentration in and around the coral fragment. (A–C) coral fragment aligned with the flow; (D–F) coral fragment across the flow. (A, D) 3D O2 distribution on the coral (tissue/water) surface. (B, E) 2D cross-sections of O2 concentration without ciliary movement. (C, F) 2D cross-sections of oxygen concentration with ciliary movement. The flow direction and 2D sectioning planes are indicated in Figures 3A, E.
When the coral fragment was oriented with one branch facing the flow (named aligned with the flow in Figure 3A), the groove and the second branch exhibited a thicker DBL (Figure 3B), as compared to the orientation across the flow (Figures 3E, F) where the groove of the fragment was more flow-exposed. By constructing iso-velocity lines (IVL) at, for example, 1 mm s-1, the hydrodynamic boundary layer, in the more closed-groove orientation appears thicker (Figure S8A). This indicates a retarded water flow, as compared to a more open-groove orientation (Figure S8B) that showed IVL more conformal and closer to the coral surface. Vertical O2 concentration profiles extracted from the simulated O2 distributions (Figure S5) were compared with O2 microsensor measurements done within similar areas of the fragment, and showed a fairly good match (Figure 5) for the two different orientations of the fragment with respect to the water flow.
Figure 5 Comparison between microsensor-measured O2 profiles and the corresponding simulated profiles through the water column and tissue, with and without considering the surface ciliary motion. Coral fragment, (A–C) aligned with the flow; (D–F) across the flow. “center”, “tip 1” and “tip 2” measurement positions are indicated in Figure 4. The tissue/water boundary indicated here was determined from the simulated DO profile and compared with the measurements.
Under low or no flow conditions, mass transfer and therefore O2 concentration in corals can also be affected by the movement of cilia covering the coral ectoderm, which can create vortices and some advective transport at the coral tissue surface (Pacherres et al., 2020; Pacherres et al., 2022). We investigated the role of the ciliated coral tissue surface by comparing O2 microenvironment simulations with and without ciliary beating (Figures 4B, C). When no cilia beating was included, the simulations indicate a pronounced O2 accumulation at the center of the fragment (Figure 4B). Whereas, when cilia movement was included, the simulated O2 distribution was more homogenous due to enhanced advective transport by the ciliary vortex formation (Figure 4C). This in line with experimental data (Pacherres et al., 2020; Ahmerkamp et al., 2022; Pacherres et al., 2022) showing that ciliary movements can to some extend homogenize the O2 concentration field over the coral surface. However, the role of cilia for coral mass transfer has only been investigated in two coral species (Shapiro et al., 2014; Pacherres et al., 2022) using 2D images of flow and O2 concentration, and there is a need for more detailed experimental studies of the 3D flow and O2 concentration fields over coral surfaces to test our model predictions. The comparison between simulated O2 profiles and the corresponding measurements (Figures 5A–C) also showed a closer match when cilia beating was included in the model. However, when the coral fragment was oriented across the flow, the central area was more exposed and the MBL was relatively thin. Consequently, the measured and simulated O2 concentration profiles indicated a smaller impact of the ciliary movements on mass transfer (Figures 4D–F, 5D–F). This suggests that cilia-induced advective transport has a more significant effect on the O2 transport between the coral tissue and the surrounding seawater in colony regions with slow flow. We note that movements in the coral such as tissue-contraction and expansion and tentacle movement, which are presently not included in the model, could also potentially influence the O2 mass transport (Patterson, 1992; Malul et al., 2020).
3.3 Temperature microenvironment
Heat transfer simulations and temperature measurements were executed for both orientations of the coral sample with respect to the flow (Figure 6 - aligned; Figure S9 - across). The simulated temperature profiles with and without ciliary movement were very similar, showing an increase of less than 0.2°C, as visible in the 2D cross-sectional images of the computed temperature distribution (Figures 6B, C and Figures S9B, C) and extracted profiles (Figures 6D–F and Figures S9D–F). Temperature profiles measured at three different locations (Figures 6A, S9A) in the fragment were within the range of simulated temperatures, but with considerable noise. The water flow affected local heat transfer, where the temperature simulations show higher tissue heating from the absorbed light in fragment regions exposed to slower flow, albeit the absolute temperature increase was very moderate for the investigated branched coral.
Figure 6 Comparison between measured and simulated temperature distribution for the coral fragment aligned with the flow. The temperature differences ΔT are between the local values and the inflow temperature. (A) computed 3D ΔT on the coral (tissue/water) surface; (B, C) ΔT in 2D cross-sections through the coral without and with ciliary movement. Sectioning planes are indicated in Figure 3A. (D–F) measured and simulated temperature difference profiles through the water column and tissue, with and without considering the surface ciliary motion.
Simulations indicate that the slower heat transfer (thicker MBL and TBL) over tissue regions less exposed to flow was not significantly enhanced by ciliary movement. The simulated temperature increase was highest near the skeleton/tissue interface, due to heat-insulating properties of the skeleton (Jimenez et al., 2008). Clearly, it is difficult to make a systematic comparison between simulated and measured temperature profiles in this case (Figures 6D–F, S9D–F) due to the very small temperature gradient between coral tissue and water close to the noise level of the microsensor measurements. The small temperature differences in the branched coral fragment are due to a high surface area to volume ratio (Jimenez et al., 2008), allowing a quick heat dissipation into the environment. There is thus a need for more accurate temperature microsensors, but it could also be relevant to compare simulated and measured temperatures on massive corals that have been shown to exhibit a stronger heating than branched corals (Jimenez et al., 2008; Jimenez et al., 2011; Jimenez et al., 2012). There is a lack of thermal property measurements in the different layers in coral tissue, and our simulation of the temperature microenvironment thus largely relied on values for tissue thermal properties from the biomedical literature.
3.4 Outlook
We developed a multiphysics modelling approach that enables simulation of the physico-chemical microenvironment of corals with a known 3D structure under defined irradiance and laminar flow conditions similar to commonly used flow chamber experimental setups. The model results were evaluated on 3D scanned coral samples, previously characterized with microsensor measurements of light, temperature and O2 in a laminar flow chamber. Such comparison generally showed a good agreement between measured and simulated data, but several model improvements can be considered for future work. The current model only considers light of individual wavelengths at a time and does not include inelastic scattering (e.g. conversion of light energy due to fluorescence). Broadband spectral simulations would e.g. enable simulations of how the spectral red shifts caused by fluorescent host pigments (Salih et al., 2000; Alieva et al., 2008), wavelength-specific reflectivity by chromoproteins (Alieva et al., 2008; D’Angelo et al., 2008) affect the coral light field and photosynthesis (Ben-Zvi et al., 2021), internal heat generation (Quick et al., 2018), as well as coral bleaching and recovery (Grinblat et al., 2018). However, an important prerequisite for more detailed simulations would be a better quantification of the inherent optical and thermal properties of coral tissue and skeleton, which are still very scarce in the literature. The present models therefore partially rely on assumed optical and thermal property values taken from biomedical tissue studies.
The simple laminar flow scenario currently implemented in the present study does certainly not cover all flow scenarios of corals in situ, and is more representative of calm sea conditions on reef flats and inside coral patches (Jimenez et al., 2011; Jimenez et al., 2012), as well as in many experimental flow chambers used for ecophysiological studies of coral metabolism. It is possible to include more complex flow scenarios in our modeling approach such as turbulent and oscillating flows (Ong et al., 2012a; Ong et al., 2012b), albeit at higher computational costs. Furthermore, it would be interesting to compare simulated and measured flow fields around corals, e.g. using particle imaging velocimetry (Pacherres et al., 2020) or combined measurements of flow and O2 fields around corals (Ahmerkamp et al., 2022; Pacherres et al., 2022). The latter could also enable implementation of a more detailed account for the role of ciliary movement for coral mass and heat transfer.
The present model only accounts for photosynthetic O2 production of symbionts using a simple approach based on reported quantum efficiencies, while O2 consumption by respiration in different tissue layers is partially based on published values and P/R relationships. The model could be expanded to better account for photosynthetic light saturation and photoinhibition. Considering other chemical species (inorganic and organic carbon, various acid/base couples) and reactions in our model would enable the representation of calcification (Willard et al., 2023) and carbon transfer between symbionts and host, and such work is in progress. Another desired expansion could be a more detailed account of metabolic processes in the coral skeleton, performed by the endolithic algae and microbes (Ricci et al., 2019; Tandon et al., 2022).
Furthermore, the present model assumes homogenous tissue thickness over the coral colony and does not include an accurate representation of fine scale topographic and anatomic features of coral tissue and skeleton, such as the polyps with tentacles and a complex internal gastrovascular system, and the more homogenous connective tissue between polyps (Hughes et al., 2022. Our model also represents coral tissue as a static structure. This is a simplification, as corals exhibit pronounced tissue plasticity via contraction and expansion, which can strongly affect surface area to volume ratios (Patterson, 1992) and optical properties (Wangpraseurt et al., 2017), which in turn modulate mass, heat and radiative transfer processes. More precise tomographic mapping of coral tissue and skeleton morphology and thickness in colonies e.g. with µCT scanning or OCT in combination with including fluid-structure interaction models with moving interfaces (Taherzadeh et al., 2012) could allow for simulating effects of coral topography and mechanics on the coral microenvironment. Likewise, such simulations can also be expanded to encompass larger, more complex coral colonies, but this will require substantial computational resources.
Finally, the presented 3D modelling approach can also be used to simulate different time-dependent environmental conditions (e.g., variable flow, day-night hypoxia, different solar irradiation regimes), which can help evaluating mechanisms driving coral stress responses as well as basic niche shaping factors for symbionts and microbiomes in the coral holobiont. The model could also be useful in more applied research such as in the ongoing attempts to create bionic corals (Wangpraseurt et al., 2020; Wangpraseurt et al., 2022) or for optimization of other 3D bioprinted constructs (Krujatz et al., 2022), where different designs can be evaluated and optimized. Our approach is also relevant for simulating structure-function relationships in other benthic systems such as stratified photosynthetic biofilms and aquatic plant tissue, and can also be adapted to other sessile organisms such as symbiont-bearing giant clams, ascidians, jellyfish or foraminifera.
4 Conclusion
We developed a 3D multiphysics model to simulate the spatial distribution of light, O2 and temperature (and the corresponding radiative, mass and heat transfer) across stratified coral tissue, skeleton and around natural coral morphologies exposed to a defined flow field and incident irradiance. Model results compared well with spatial measurements of light, O2 and temperature under similar flow and light conditions. The model reveals how the interaction between incident irradiance, water flow and complex coral morphology leads to pronounced spatial heterogeneity and microenvironments, both across tissue layers and between different areas and compartments in coral colonies. Such simulation of the physico-chemical microenvironmental landscape using real 3D scanned coral structures as an input, can i) give insights to coral tissue and skeleton compartments that are difficult to reach with existing sensor technology, and ii) identify hotspots of activity that can inform detailed measurements e.g. with microsensors and/or various bioimaging techniques for mapping structure and function. The model can simulate effects of different environmental (e.g. light, temperature, hypoxia or flow) and structural (e.g. symbiont and host pigment density and distribution, and the distribution of endolithic microbes) factors on the coral microenvironment and metabolic activity, which can be tested experimentally. We argue that such combination of modeling and experimental investigation is a strong tool set for unravelling structure-function relations and basic regulatory mechanisms in coral biology and stress responses, including effects of climate change and other anthropogenic threats.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
SM, CP and MK designed the research. SM developed the coral model, performed computer simulations and measurements. SM, CP and MK analyzed the data. SM, CP and MK wrote the manuscript. CP and MK provided research infrastructure. All authors contributed to the article and approved the submitted version.
Funding
This study was funded by the Gordon and Betty Moore Foundation through grant no. GBMF9206 to MK (https://doi.org/10.37807/GBMF9206). CP acknowledges support and access to computational resources offered by the KAUST Computing Center.
Acknowledgments
We acknowledge technical assistance with coral husbandry by Sofie Jakobsen, Victoria Thuesen and Caroline Vigsbo Christensen.
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/fmars.2023.1160208/full#supplementary-material
References
Ahmerkamp S., Jalaluddin F. M., Cui Y., Brumley D. R., Pacherres C. O., Berg J. S., et al. (2022). Simultaneous visualization of flow fields and oxygen concentrations to unravel transport and metabolic processes in biological systems. Cell Rep. Methods 2 (5), 100216. doi: 10.1016/j.crmeth.2022.100216
Al-Horani F. A., Al-Moghrabi S. M., de Beer D. (2003). The mechanism of calcification and its relation to photosynthesis and respiration in the scleractinian coral Galaxea fascicularis. Mar. Biol. 142 (3), 419–426. doi: 10.1007/s00227-002-0981-8
Alieva N. O., Konzen K. A., Field S. F., Meleshkevitch E. A., Hunt M. E., Beltran-Ramirez V., et al. (2008). Diversity and evolution of coral fluorescent proteins. PloS One 3 (7), e2680. doi: 10.1371/journal.pone.0002680
Altieri A. H., Harrison S. B., Seemann J., Collin R., Diaz R. J., Knowlton N. (2017). Tropical dead zones and mass mortalities on coral reefs. Proc. Natl. Acad. Sci. 114 (14), 3660–3665. doi: 10.1073/pnas.1621517114
Anthony K. R. N., Hoegh-Guldberg O. (2003). Variation in coral photosynthesis, respiration and growth characteristics in contrasting light microhabitats: an analogue to plants in forest gaps and understoreys? Funct. Ecol. 17 (2), 246–259. doi: 10.1046/j.1365-2435.2003.00731.x
Ben-Zvi O., Wangpraseurt D., Bronstein O., Eyal G., Loya Y. (2021). Photosynthesis and bio-optical properties of fluorescent mesophotic corals. Front. Mar. Sci. 8, 651601. doi: 10.3389/fmars.2021.651601
Bollati E., Lyndby N. H., D’Angelo C., Kühl M., Wiedenmann J., Wangpraseurt D. (2022). Green fluorescent protein-like pigments optimise the internal light environment in symbiotic reef-building corals. Elife 11, e73521. doi: 10.7554/eLife.73521
Bowie G., Mills W., Porcella D., Campbell C., Chamberlin C. (1985). Rates, constants, and kinetics formulations in surface water quality modeling. Washington, D.C.: SECOND EDITION. U.S. Environmental Protection Agency.
Bragg W. L.(1924). The structure of aragonite. Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character. 105 (729), 16–39. doi: 10.1098/rspa.1924.0002
Brodersen K. E., Lichtenberg M., Ralph P. J., Kühl M., Wangpraseurt D. (2014). Radiative energy budget reveals high photosynthetic efficiency in symbiont-bearing corals. J. R. Soc. Interface 11 (93), 20130997. doi: 10.1098/rsif.2013.0997
Chalker B. E. (1981). Simulating light-saturation curves for photosynthesis and calcification by reef-building corals. Mar. Biol. 63 (2), 135–141. doi: 10.1007/BF00406821
Chan N. C. S., Wangpraseurt D., Kühl M., Connolly S. R. (2016). Flow and coral morphology control coral surface pH: implications for the effects of ocean acidification. Front. Mar. Sci. 3, 10. doi: 10.3389/fmars.2016.00010
Chang S., Iaccarino G., Ham F., Elkins C., Monismith S. (2014). Local shear and mass transfer on individual coral colonies: computations in unidirectional and wave-driven flows. J. Geophysical Research: Oceans 119 (4), 2599–2619. doi: 10.1002/2013JC009751
Chindapol N., Kaandorp J. A., Cronemberger C., Mass T., Genin A. (2013). Modelling growth and form of the scleractinian coral Pocillopora verrucosa and the influence of hydrodynamics. PloS Comput. Biol. 9 (1), e1002849. doi: 10.1371/journal.pcbi.1002849
Cooper T. F., Ulstrup K. E., Dandan S. S., Heyward A. J., Kühl M., Muirhead A., et al. (2011). Niche specialization of reef-building corals in the mesophotic zone: metabolic trade-offs between divergent Symbiodinium types. Proc. R. Soc. B: Biol. Sci. 278 (1713), 1840–1850. doi: 10.1098/rspb.2010.2321
Costanza R., de Groot R., Sutton P., van der Ploeg S., Anderson S. J., Kubiszewski I., et al. (2014). Changes in the global value of ecosystem services. Global Environ. Change 26, 152–158. doi: 10.1016/j.gloenvcha.2014.04.002
D’Angelo C., Denzel A., Vogt A., Matz M., Oswald F., Salih A., et al. (2008). Blue light regulation of host pigment in reef-building corals. Mar. Ecol. Prog. Ser. 364, 97–106. doi: 10.3354/meps07588
Davy S. K., Allemand D., Weis V. M. (2012). Cell biology of cnidarian-dinoflagellate symbiosis. Microbiol. Mol. Biol. Rev. 76 (2), 229–261. doi: 10.1128/MMBR.05014-11
DeWitt D. P., Incropera F. P. (1990). Fundamentals of Heat and Mass Transfer (3 ed.). John Wiley & Sons.
Dubinsky Z., Falkowski P. G., Porter J. W., Muscatine L., Smith D. C. (1984). Absorption and utilization of radiant energy by light- and shade-adapted colonies of the hermatypic coral Stylophora pistillata. Proc. R. Soc. London. Ser. B. Biol. Sci. 222 (1227), 203–214. doi: 10.1098/rspb.1984.0059
Dunne R. P., Brown B. E. (2001). The influence of solar radiation on bleaching of shallow water reef corals in the Andaman Sea 1993-1998. Coral Reefs 20 (3), 201–210. doi: 10.1007/s003380100160
Enríquez S., Méndez E. R., Hoegh-Guldberg O., Iglesias-Prieto R. (2017). Key functional role of the optical properties of coral skeletons in coral ecology and evolution. Proc. R. Soc. B. 284, 20161667. doi: 10.1098/rspb.2016.1667
Enríquez S., Méndez E. R., Prieto R. I. (2005). Multiple scattering on coral skeletons enhances light absorption by symbiotic algae. Limnol. Oceanogr. 50 (4), 1025–1032. doi: 10.4319/lo.2005.50.4.1025
Eyal G., Wiedenmann J., Grinblat M., D’Angelo C., Kramarsky-Winter E., Treibitz T., et al. (2015). Spectral diversity and regulation of coral fluorescence in a mesophotic reef habitat in the red Sea. PloS One 10 (6), e0128697. doi: 10.1371/journal.pone.0128697
Falter J. L., Atkinson M. J., Merrifield M. A. (2004). Mass-transfer limitation of nutrient uptake by a wave- dominated reef flat community. Limnology Oceanography 49), 1820–1831. doi: 10.4319/lo.2004.49.5.1820
Fares E., Schröder W. (2002). A differential equation for approximate wall distance. Int. J. Numerical Methods Fluids 39 (8), 743–762. doi: 10.1002/fld.348
Ghosh G.(1999). Dispersion-equation coefficients for the refractive index and birefringence of calcite and quartz crystals. Optics Communications 163 (1), 95–102, doi: 10.1016/S0030-4018(99)00091-7
Goreau T. J., Goreau N. I., Trench R. K., Hayes R. L. (1996). Calcification rates in corals. Science 274 (5284), 117–117. doi: 10.1126/science.274.5284.117-b
Grinblat M., Fine M., Tikochinski Y., Loya Y. (2018). Stylophora pistillata in the red Sea demonstrate higher GFP fluorescence under ocean acidification conditions. Coral Reefs 37 (1), 309–320. doi: 10.1007/s00338-018-1659-0
Hasgall P. A., Baumgartner C., Neufeld E., Lloyd B., Gosselin M. C., Payne D., et al (2018). IT’IS Database for thermal and electromagnetic parameters of biological tissues. doi: 10.13099/VIP21000-04-0
Hasgall P. A., Di Gennaro F., Baumgartner C., Neufeld E., Lloyd B., Gosselin M. C., Payne D., Klingenböck A., Kuster N. (2022). “IT’IS database for thermal and electromagnetic parameters of biological tissues,” Version 4.1, Feb 22, 2022. doi: 10.13099/VIP21000-04-1. itis.swiss/database
Haynes W. M., Lide D. R., Bruno T. J. (2016). CRC handbook of chemistry and physics: a ready-reference book of chemical and physical data. 2016-2017, 97th Edition / Boca Raton, Florida, CRC Press.
Holcomb M., DeCarlo T. M., Schoepf V., Dissard D., Tanaka K., McCulloch M. (2015). Cleaning and pre-treatment procedures for biogenic and synthetic calcium carbonate powders for determination of elemental and boron isotopic compositions. Chem. Geology 398, 11–21. doi: 10.1016/j.chemgeo.2015.01.019
Hughes D. J., Alderdice R., Cooney C., Kühl M., Pernice M., Voolstra C. R., et al. (2020). Coral reef survival under accelerating ocean deoxygenation. Nat. Climate Change 10 (4), 296–307. doi: 10.1038/s41558-020-0737-9
Hughes T. P., Kerry J. T., Álvarez-Noriega M., Álvarez-Romero J. G., Anderson K. D., Baird A. H., et al. (2017). Global warming and recurrent mass bleaching of corals. Nature 543 (7645), 373–377. doi: 10.1038/nature21707
Hughes D. J., Raina J.-B., Nielsen D. A., Suggett D. J., Kühl M. (2022). Disentangling compartment functions in sessile marine invertebrates. Trends Ecol. Evol. 37 (9), 740–748. doi: 10.1016/j.tree.2022.04.008
Jackson J. B. C., Kirby M. X., Berger W. H., Bjorndal K. A., Botsford L. W., Bourque B. J., et al. (2001). Historical overfishing and the recent collapse of coastal ecosystems. Science 293 (5530), 629–637. doi: 10.1126/science.1059199
Jacques S. L., Wangpraseurt D., Kühl M. (2019). Optical properties of living corals determined with diffuse reflectance spectroscopy. Front. Mar. Sci. 6, 472. doi: 10.3389/fmars.2019.00472
Jimenez I. M., Kühl M., Larkum A. W. D., Ralph P. J. (2008). Heat budget and thermal microenvironment of shallow-water corals: do massive corals get warmer than branching corals? Limnol. Oceanogr. 53 (4), 1548–1561. doi: 10.4319/lo.2008.53.4.1548
Jimenez I. M., Kühl M., Larkum A. W. D., Ralph P. J. (2011). Effects of flow and colony morphology on the thermal boundary layer of corals. J. R. Soc. Interface 8 (65), 1785–1795. doi: 10.1098/rsif.2011.0144
Jimenez I. M., Larkum A. W. D., Ralph P. J., Kühl M. (2012). In situ thermal dynamics of shallow water corals is affected by tidal patterns and irradiance. Mar. Biol. 159 (8), 1773–1782. doi: 10.1007/s00227-012-1968-8
Johnson M. D., Rodriguez L. M., Altieri A. H. (2018). Shallow-water hypoxia and mass mortality on a Caribbean coral reef. Bull. Mar. Sci. 94 (1), 143–144. doi: 10.5343/bms.2017.1163
Kaandorp J. A. (2013). Macroscopic modelling of environmental influence on growth and form of sponges and corals using the accretive growth model. Int. Scholarly Res. Notices 2013, 1–14. doi: 10.1155/2013/159170
Kaniewska P., Magnusson S. H., Anthony K. R., Reef R., Kühl M., Hoegh-Guldberg O. (2011). Importance of macro- versus microstructure in modulating light levels inside coral colonies. J. Phycol 47 (4), 846–860. doi: 10.1111/j.1529-8817.2011.01021.x
Kaniewska P., Sampayo E. M. (2022). Macro- and micro-scale adaptations allow distinct Stylophora pistillata-symbiodiniaceae holobionts to optimize performance across a broad light habitat. J. Phycol 58 (1), 55–70. doi: 10.1111/jpy.13215
Kirk J. T. O. (1994). Light and Photosynthesis in Aquatic Ecosystems (2 ed.). Cambridge University Press. doi: 10.1017/CBO9780511623370
Knowlton N., Corcoran E., Felis T., de Goeij J., Grottoli A., Harding S., et al. (2021). Rebuilding coral reefs: a decadal grand challenge (International Coral Reef Society and Future Earth Coasts), 56. doi: 10.53642/NRKY9386
Kramer N., Guan J., Chen S., Wangpraseurt D., Loya Y. (2022). Morpho-functional traits of the coral Stylophora pistillata enhance light capture for photosynthesis at mesophotic depths. Commun. Biol. 5 (1), 861. doi: 10.1038/s42003-022-03829-4
Krujatz F., Dani S., Windisch J., Emmermacher J., Hahn F., Mosshammer M., et al. (2022). Think outside the box: 3D bioprinting concepts for biotechnological applications - recent developments and future perspectives. Biotechnol. Adv. 58, 107930. doi: 10.1016/j.biotechadv.2022.107930
Kühl M., Holst G., Larkum A. W., Ralph P. J. (2008). Imaging of oxygen dynamics within the endolithic algal community of the massive coral Porites lobata. J. Phycol 44 (3), 541–550. doi: 10.1111/j.1529-8817.2008.00506.x
LaJeunesse T. C., Parkinson J. E., Gabrielson P. W., Jeong H. J., Reimer J. D., Voolstra C. R., et al. (2018). Systematic revision of symbiodiniaceae highlights the antiquity and diversity of coral endosymbionts. Curr. Biol. 28 (16), 2570–2580.e2576. doi: 10.1016/j.cub.2018.07.008
Leino A. A., Pulkkinen A., Tarvainen T. (2019). ValoMC: a Monte Carlo software and MATLAB toolbox for simulating light transport in biological tissue. Osa Continuum 2 (3), 957–972. doi: 10.1364/Osac.2.000957
Lenihan H. S., Hench J. L., Holbrook S. J., Schmitt R. J., Potoski M. (2015). Hydrodynamics influence coral performance through simultaneous direct and indirect effects. Ecology 96 (6), 1540–1549. doi: 10.1890/14-1115.1
Lichtenberg M., Brodersen K. E., Kühl M. (2017). Radiative energy budgets of phototrophic surface-associated microbial communities and their photosynthetic efficiency under diffuse and collimated light. Front. Microbiol. 8, 452. doi: 10.3389/fmicb.2017.00452
Lichtenberg M., Larkum A. W. D., Kühl M. (2016). Photosynthetic acclimation of Symbiodinium in hospite depends on vertical position in the tissue of the scleractinian coral Montastrea curta. Front. Microbiol. 7, 230. doi: 10.3389/fmicb.2016.00230
Loya Y., Sakai K., Yamazato K., Nakano Y., Sambali H., Van Woesik R. (2001). Coral bleaching: the winners and the losers. Ecol. Lett. 4 (2), 122–131. doi: 10.1046/j.1461-0248.2001.00203.x
Lyndby N. H., Kühl M., Wangpraseurt D. (2016). Heat generation and light scattering of green fluorescent protein-like pigments in coral tissue. Sci. Rep. 6 (1), 26599. doi: 10.1038/srep26599
MacDougall J. D., McCabe M. (1967). Diffusion coefficient of oxygen through tissues. Nature 215 (5106), 1173–1174. doi: 10.1038/2151173a0
Magnusson S. H., Fine M., Kühl M. (2007). Light microclimate of endolithic phototrophs in the scleractinian corals Montipora monasteriata and Porites cylindrica. Mar. Ecol. Prog. Ser. 332, 119–112. doi: 10.3354/meps332119
Malul D., Holzman R., Shavit U. (2020). Coral tentacle elasticity promotes an out-of-phase motion that improves mass transfer. Proc. R. Soc. B: Biol. Sci. 287 (1929), 20200180. doi: 10.1098/rspb.2020.0180
Moberg F., Folke C. (1999). Ecological goods and services of coral reef ecosystems. Ecol. Economics 29 (2), 215–233. doi: 10.1016/S0921-8009(99)00009-9
Mumby P. J., Chisholm J. R. M., Edwards A. J., Andrefouet S., Jaubert J. (2001). Cloudy weather may have saved society island reef corals during the 1998 ENSO event. Mar. Ecol. Prog. Ser. 222, 209–216. doi: 10.3354/meps222209
Muscatine L. (1973). “4. - NUTRITION OF CORALS,” in Biology and geology of coral reefs. Eds. Jones O. A., Endean R. (Academic Press), 77–115. doi: 10.1016/B978-0-12-395526-5.50012-2
Nakamura T., Yamasaki H., Van Woesik R. (2003). Water flow facilitates recovery from bleaching in the coral Stylophora pistillata. Mar. Ecol. Prog. Ser. 256, 287–291. doi: 10.3354/meps256287
Ong R. H., King A. J. C., Caley M. J., Mullins B. J. (2018). Prediction of solar irradiance using ray-tracing techniques for coral macro- and micro-habitats. Mar. Environ. Res. 141, 75–87. doi: 10.1016/j.marenvres.2018.08.004
Ong R. H., King A. J. C., Kaandorp J. A., Mullins B. J., Caley M. J. (2017). The effect of allometric scaling in coral thermal microenvironments. PloS One 12 (10), e0184214. doi: 10.1371/journal.pone.0184214
Ong R. H., King A. J. C., Mullins B. J., Caley M. J. (2019). The effect of small-scale morphology on thermal dynamics in coral microenvironments. J. Thermal Biol. 86, 102433. doi: 10.1016/j.jtherbio.2019.102433
Ong R. H., King A. J. C., Mullins B. J., Caley M. J., Cooper T. F. (2012b) CFD Simulation of low Reynolds-number turbulence models in coral thermal microenvironments. 18th Australasian Fluid Mechanics Conference, Launceston, Australia, Dec 3-7 2012 (Australasian Fluid Mechanics Society). Available at: http://hdl.handle.net/20.500.11937/13242.
Ong R. H., King A. J., Mullins B. J., Cooper T. F., Caley M. J. (2012a). Development and validation of computational fluid dynamics models for prediction of heat transfer and thermal microenvironments of corals. PloS One 7 (6), e37842. doi: 10.1371/journal.pone.0037842
Pacherres C. O., Ahmerkamp S., Koren K., Richter C., Holtappels M. (2022). Ciliary flows in corals ventilate target areas of high photosynthetic oxygen production. Curr. Biol. 32 (19), 4150–4158. doi: 10.1016/j.cub.2022.07.071
Pacherres C. O., Ahmerkamp S., Schmidt-Grieb G. M., Holtappels M., Richter C. (2020). Ciliary vortex flows and oxygen dynamics in the coral boundary layer. Sci. Rep. 10 (1), 7541. doi: 10.1038/s41598-020-64420-7
Patterson M. R. (1992). A chemical engineering view of cnidarian symbioses. Am. Zoologist 32 (4), 566–582. doi: 10.1093/icb/32.4.566
Quick C., D’Angelo C., Wiedenmann J. (2018). Trade-offs associated with photoprotective green fluorescent protein expression as potential drivers of balancing selection for color polymorphism in reef corals. Front. Mar. Sci. 5, 11. doi: 10.3389/fmars.2018.00011
Revsbech N. P. (1989). An oxygen microsensor with a guard cathode. Limnology Oceanography 34 (2), 474–478. doi: 10.4319/lo.1989.34.2.0474
Ricci F., Rossetto Marcelino V., Blackall L. L., Kühl M., Medina M., Verbruggen H. (2019). Beneath the surface: community assembly and functions of the coral skeleton microbiome. Microbiome 7 (1), 159. doi: 10.1186/s40168-019-0762-y
Ricci F., Tandon K., Mosshammer M., Cho E., Blackall L., Kühl M., et al. (2023). Fine-scale mapping of physicochemical and microbial landscapes of the coral skeleton. Environ. Microbiol. doi: 10.1111/1462-2920.16369
Rickelt L. F., Lichtenberg M., Trampe E. C. L., Kühl M. (2016). Fiber-optic probes for small-scale measurements of scalar irradiance. Photochem. Photobiol. 92 (2), 331–342. doi: 10.1111/php.12560
Salih A., Larkum A., Cox G., Kühl M., Hoegh-Guldberg O. (2000). Fluorescent pigments in corals are photoprotective. Nature 408 (6814), 850–853. doi: 10.1038/35048564
Shapiro O. H., Fernandez V. I., Garren M., Guasto J. S., Debaillon-Vesque F. P., Kramarsky-Winter E., et al. (2014). Vortical ciliary flows actively enhance mass transport in reef corals. Proc. Natl. Acad. Sci. 111 (37), 13391–13396. doi: 10.1073/pnas.1323094111
Shashar N., Cohen Y., Loya Y. (1993). Extreme diel fluctuations of oxygen in diffusive boundary layers surrounding stony corals. Biol. Bull. 185 (3), 455–461. doi: 10.2307/1542485
Smee D. L., Ferner M. C., Weissburg. M. J. (2010). Hydrodynamic sensory stressors produce nonlinear predation patterns. Ecology 91), 1391–1314. doi: 10.1890/09-0017.1
Spicer G. L. C., Eid A., Wangpraseurt D., Swain T. D., Winkelmann J. A., Yi J., et al. (2019). Measuring light scattering and absorption in corals with inverse spectroscopic optical coherence tomography (ISOCT): a new tool for non-invasive monitoring. Sci. Rep. 9 (1), 14148. doi: 10.1038/s41598-019-50658-3
Swain T. D., DuBois E., Gomes A., Stoyneva V. P., Radosevich A. J., Henss J., et al. (2016). Skeletal light-scattering accelerates bleaching response in reef-building corals. BMC Ecol. 16, 10. doi: 10.1186/s12898-016-0061-4
Tagliafico A., Baker P., Kelaher B., Ellis S., Harrison D. (2022). The effects of shade and light on corals in the context of coral bleaching and shading technologies. Front. Mar. Sci. 9, 919382. doi: 10.3389/fmars.2022.919382
Taherzadeh D., Picioreanu C., Horn H. (2012). Mass transfer enhancement in moving biofilm structures. Biophys. J. 102 (7), 1483–1492. doi: 10.1016/j.bpj.2012.02.033
Tamir R., Eyal G., Kramer N., Laverick J. H., Loya Y. (2019). Light environment drives the shallow-to-mesophotic coral community transition. Ecosphere 10 (9), e02839. doi: 10.1002/ecs2.2839
Tandon K., Pasella M. M., Iha C., Ricci F., Hu J., O’Kelly C. J., et al. (2022). Every refuge has its price: Ostreobium as a model for understanding how algae can live in rock and stay in business. Semin. Cell Dev. Biol. 134, 27–36. doi: 10.1016/j.semcdb.2022.03.010
Taylor Parkins S. K., Murthy S., Picioreanu C., Kühl M. (2021). Multiphysics modelling of photon, mass and heat transfer in coral microenvironments. J. R. Soc. Interface 18 (182), 20210532. doi: 10.1098/rsif.2021.0532
Teneva L., Karnauskas M., Logan C. A., Bianucci L., Currie J. C., Kleypas J. A. (2012). Predicting coral bleaching hotspots: the role of regional variability in thermal stress and potential adaptation rates. Coral Reefs 31 (1), 1–12. doi: 10.1007/s00338-011-0812-9
Terán E., Méndez E. R., Enríquez S., Iglesias-Prieto R. (2010). Multiple light scattering and absorption in reef-building corals. Appl. Optics 49 (27), 5032–5042. doi: 10.1364/AO.49.005032
Todd P. A. (2008). Morphological plasticity in scleractinian corals. Biol. Rev. Camb Philos. Soc. 83 (3), 315–337. doi: 10.1111/j.1469-185x.2008.00045.x
van der Zande R. M., Achlatis M., Bender-Champ D., Kubicek A., Dove S., Hoegh-Guldberg O. (2020). Paradise lost: end-of-century warming and acidification under business-as-usual emissions have severe consequences for symbiotic corals. Global Change Biol. 26 (4), 2203–2219. doi: 10.1111/gcb.14998
Van Woesik R., Irikawa A., Anzai R., Nakamura T. (2012). Effects of coral colony morphologies on mass transfer and susceptibility to thermal stress. Coral Reefs 31 (3), 633–639. doi: 10.1007/s00338-012-0911-2
Veal C. J., Carmi M., Dishon G., Sharon Y., Michael K., Tchernov D., et al. (2010). Shallow-water wave lensing in coral reefs: a physical and biological case study. J. Exp. Biol. 213 (24), 4304–4312. doi: 10.1242/jeb.044941
Wangpraseurt D., Jacques S., Lyndby N., Holm J. B., Pages C. F., Kühl M. (2019). Microscale light management and inherent optical properties of intact corals studied with optical coherence tomography. J. R. Soc. Interface 16, 20180567. doi: 10.1101/376723
Wangpraseurt D., Jacques S. L., Petrie T., Kühl M. (2016). Monte Carlo Modeling of photon propagation reveals highly scattering coral tissue. Front. Plant Sci. 7, 1404. doi: 10.3389/fpls.2016.01404
Wangpraseurt D., Larkum A. W., Franklin J., Szabó M., Ralph P. J., Kühl M. (2014a). Lateral light transfer ensures efficient resource distribution in symbiont-bearing corals. J. Exp. Biol. 217 (Pt 4), 489–498. doi: 10.1242/jeb.091116
Wangpraseurt D., Larkum A., Ralph P., Kühl M. (2012). Light gradients and optical microniches in coral tissues. Front. Microbiol. 3, 316. doi: 10.3389/fmicb.2012.00316
Wangpraseurt D., Polerecky L., Larkum A. W. D., Ralph P. J., Nielsen D. A., Pernice M., et al. (2014b). The in situ light microenvironment of corals. Limnol Oceanogr 59 (3), 917–926. doi: 10.4319/lo.2014.59.3.0917
Wangpraseurt D., Sun Y., You S., Chua S.-T., Noel S. K., Willard H. F., et al. (2022). Bioprinted living coral microenvironments mimicking coral-algal symbiosis. Adv. Funct. Materials 32 (35), 2202273. doi: 10.1002/adfm.202202273
Wangpraseurt D., Wentzel C., Jacques S. L., Wagner M., Kühl M. (2017). In vivo imaging of coral tissue and skeleton with optical coherence tomography. J. R. Soc. Interface 14, 1420161003. doi: 10.1098/rsif.2016.1003
Wangpraseurt D., You S., Azam F., Jacucci G., Gaidarenko O., Hildebrand M., et al. (2020). Bionic 3D printed corals. Nat. Commun. 11 (1), 1748. doi: 10.1038/s41467-020-15486-4
Keywords: heat transfer, mass transfer, microenvironment, numerical simulation, radiative transfer
Citation: Murthy S, Picioreanu C and Kühl M (2023) Modeling the radiative, thermal and chemical microenvironment of 3D scanned corals. Front. Mar. Sci. 10:1160208. doi: 10.3389/fmars.2023.1160208
Received: 06 February 2023; Accepted: 02 May 2023;
Published: 22 May 2023.
Edited by:
Aldo Cróquer, The Nature Conservancy, Dominican RepublicReviewed by:
Ian C. Enochs, Atlantic Oceanographic and Meteorological Laboratory (NOAA), United StatesDirk de Beer, Max Planck Society, Germany
Copyright © 2023 Murthy, Picioreanu and Kühl. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Swathi Murthy, c3dhdGhpLm11cnRoeUBiaW8ua3UuZGs=