Plants are able to survive a variety of environmental conditions due in part to their ability to detect and respond to their environment. Although ultraviolet (UV) light has the potential to damage DNA and photosynthetic machinery, plants can minimize its negative effects by synthesizing UV-absorbing compounds (Li et al., 1993; Landry et al., 1995) such as flavonoids, hydroxycinnamic acids (HCAs), and mycosporine-like amino acids (Beggs and Wellmann, 1994; Cockell and Knowland, 1999; Agati et al., 2013). These phenolic secondary metabolites can protect plants from UV damage by decreasing the transmittance of UV photons through tissue. Given that UV-absorbing compounds are energetically expensive to produce (Weinig et al., 2004), most plant species plastically upregulate their production in response to UV-exposure (Lois, 1994; Dixon and Paiva, 1995). When these compounds are energetically expensive yet adaptive (i.e., when they increase fitness in the presence of UV, but decrease fitness in the absence of UV), natural selection will favor individuals that optimally allocate these compounds to cells that shield valuable tissues from UV exposure (Weinig et al., 2004). For example, UV-absorbing compounds are concentrated within trichomes and/or epidermal cells on the adaxial surface of leaves, which decreases the transmittance of UV-light into the mesophyll cells beneath (Cen and Bornman, 1993; Reuber et al., 1996; Schnitzler et al., 1996; Burchard et al., 2000; Tattini et al., 2000; Bilger et al., 2001; Agati et al., 2002). The reduction in epidermal transmittance in turn reduces damage to photosynthetic machinery (Tevini et al., 1991; Kolb et al., 2001) and DNA (Mazza et al., 2000). There is also evidence that plants nonrandomly distribute UV-absorbing compounds within their reproductive organs (Day and Demchik, 1996a, 1996b) and seeds (Griffen et al., 2004).
The concentration of UV-absorbing compounds may also vary among leaves within canopies. Previous studies have subdivided leaves into broad categories (i.e., “unshaded vs. shaded” or “upper vs. lower canopy”) and found higher concentrations of UV-absorbing compounds in “unshaded” or “upper” leaves (Lovelock et al., 1992; McKee, 1995). However, using categorical variables to quantify leaf position does not sufficiently capture the differences in microenvironments experienced by leaves throughout the entire three-dimensional space of a canopy. For example, leaves on the exterior of canopies alter both the quantity and quality of light received by leaves on the interior of canopies (Le Roux et al., 2001; Baldocchi et al., 2002; Valladares, 2003). Although multiple methods exist for quantifying leaf position within a canopy, they are often challenging to implement due to both financial and/or logistical constraints in the field. For example, digitizing methods (i.e., sound propagation: Sinoquet et al., 1991; Room et al., 1996; magnetic field current induction: Sinoquet et al., 1998; Le Roux et al., 2001; Everhart et al., 2011; and LiDAR: Greaves et al., 2015; Magney et al., 2016; Swatantran et al., 2016) are expensive, require electrical power supplies, and fail to access the inner canopy (Magney et al., 2016). Moreover, digitizing methods generate enormous data sets that often require discrete variable distinctions such as flush numbering, stem and branch classifications, or primary and lateral leaf groups to clearly observe trends within the canopy (Sinoquet et al., 1998).
We introduce a simple, inexpensive, and portable method based on a spherical polar coordinate system to quantify the exact location of leaves within the canopy of black mangrove (Avicennia germinans (L.) Stearn) trees. By quantifying leaf position with three continuous variables (radial distance, zenith angle, and azimuthal angle) and a set of mathematical combination variables (radial and vertical canopy depths), we are able to use linear regressions to determine how light quantity changes throughout the canopy and how leaf position within the canopy affects the concentration of UV-absorbing compounds. Black mangrove trees were selected for this study for several reasons. First, as inhabitants of the tropics, black mangroves receive higher levels of UV radiation than plants inhabiting temperate regions (Frederick et al., 1989; Madronich, 1993). Moreover, this species inhabits coastlines and is rarely shaded by taller trees. Consequently, black mangroves have had the opportunity to adapt to this high-UV environment by evolving strategies that minimize UV damage (Lovelock et al., 1992; Weinig et al., 2004). Black mangrove trees are also an ideal species because their canopies are large enough to create a multitude of microenvironments, but small enough to sample thoroughly.
MATERIALS AND METHODS
Plant system—The four black mangrove (A. germinans) trees included in our study were of similar sizes and occurred as single isolated trees along the beach on the west side of Hummingbird Cay, a privately owned island located approximately 13 km west of Georgetown in the Great Exuma island chain of the Bahamas. Geospatial coordinates for the trees are as follows: (23°27′35.2800″N, 75°56′39.8040″W), (23°27′36.8640″N, 75°56′39.5520″W), (23°27′37.2168″N, 75°56′39.8148″W), and (23°27′43.1964″N, 75°56′39.8256″W).
Quantifying leaf position using spherical polar coordinates—A spherical polar coordinate system was used to quantify the spatial location of 49–50 leaves within each canopy. This coordinate system uses three continuous variables (ρ, θ, ϕ) to specify leaf location. Approximately 12 leaf samples were collected per quadrant of each canopy, which yielded spatial locations with many combinations of ρ, θ, and ϕ. ρ is the radial distance or the direct distance of a point from a fixed origin, θ is the zenith angle measured from a fixed vertical direction, and ϕ is the azimuthal angle, measured as the angle relative to a fixed direction on a reference plane that passes through the origin and is orthogonal to the vertical direction (Fig. 1). The origin is defined as the trunk of the tree at the base of the canopy. Radial distance (ρ) was measured by looping a piece of twine around the trunk of the tree at the base of the canopy. A flexible tape measure was then secured to the twine with 0 cm positioned at the base of the canopy. The twine was drawn to the center of the leaf location, and the direct distance from the trunk to the leaf was measured to the nearest centimeter.
To measure the angular components (θ, ϕ) of the spatial location, two measurement protractors were prepared from printed 360° protractor images. The protractors were secured to two layers of folder stock, cut to shape, and then laminated. Materials were selected on the basis of availability, ease of construction, and field-portability. The protractor used to measure zenith angle (θ) was prepared by cutting the first 360° protractor in half along the 0° line. A second cut (∼55 mm) was made, originating from the flat end of the corresponding 180° protractor, along the 90° line (Appendix 1). The protractor was positioned such that the flat edge was flush with the tree trunk and 0° was oriented directly vertical. Measurements were recorded to the nearest degree (with a spacing of 1.75 mm per degree) and ranged from 0–90°. The protractor used to measure azimuthal angle (ϕ) was prepared using the second 360° protractor. A ∼100 mm cut was made along the 0° line toward the center of the protractor. Three additional cuts (∼55 mm) originating from the origin of the protractor were made along the 90°, 180°, and 270° lines (Appendix 2). The single long cut allowed the protractor to be pulled around the base of the canopy, and the three shorter cuts provided space to position the base of the canopy at the origin. The protractor was then secured, leveled using a bubble level, and oriented with 0° positioned at due north using a compass. Measurements of ϕ were recorded to the nearest degree and ranged from 0–360°. The 180° zenith protractor was fitted perpendicular to the azimuthal protractor using the partial cut along 90° line of the zenith protractor (Appendix 3). The angular coordinates were measured using the direct path of the twine to the leaf location relative to the two protractors. Using this method, we are able to spatially map the canopy at the resolution scale of individual leaves. Measurement errors in ρ, θ, and ϕ are ±1 cm, ±1°, and ±1°, respectively, and are below the scale of individual leaves.
Quantifying PAR and UV-absorbing compounds—After measuring the spherical polar coordinates of each leaf, the amount of photosynthetically active radiation (PAR) it received was determined by placing the sensor of a handheld light meter (LI-250; LI-COR, Lincoln, Nebraska, USA) at the adaxial leaf surface. To minimize variation in full sun (unobstructed by leaves) PAR, we only collected data between 12:00 p.m. and 3:00 p.m. on sunny days and we sampled all 49–50 leaves within a single canopy on the same day. Full sun PAR during data collection ranged from 1566 to 1719 µmol m−2 s−1. Collecting data from a subset of leaves from each tree each day would have required that we spend more time sampling each day and consequently would have increased the range of full sun PAR levels during sampling.
Following the measurement of PAR, a leaf disk (31.2 mm2) was collected from the middle of the leaf at its widest part using a cork borer. Disks were placed in individual 5-mL Eppendorf tubes containing 4 mL of a 90 : 1 : 1 methanol : HCl : H2O (v : v : v) extraction solution prepared from methanol (Thermo Fisher Scientific, Waltham, Massachusetts, USA), concentrated HCl (37% HCl; Sigma-Aldrich, St. Louis, Missouri, USA), and deionized H2O (Griffen et al., 2004). Tubes were covered in aluminum, stored at −3°C, and then transported to the laboratory. Aliquots (200 µL) of each leaf extraction solution were then diluted to total volumes of 2000 µL using the 90 : 1 :1 methanol : HCl : H2O solution. The concentration of UV-absorbing compounds in each leaf sample was quantified by measuring absorption at 300 nm (A300) using a UV-Vis spectrometer (Shimadzu UVmini-1240; Shimadzu Corporation, Kyoto, Japan).
Analyses—Although the four black mangroves included in our study were similarly sized, the maximum value of ρ differed among trees. Thus, we standardized this variable by dividing the value of ρ for each leaf within a tree data set by the maximum value of ρ within the data set. In this framework, the maximum ρ for each tree was 1.0 or 100%ρ. The new variable (%ρ) was then used in all calculations of canopy depth. Because ϕ and θ were oriented to global directions, 0° at due north and 0° at vertical, respectively, the angular measurements of individual trees could be combined into a single data set. Measurements in ϕ were converted to radians and transformed from circular to linear variables using a sine function.
We quantified leaf position by calculating radial (RD) and vertical (VD) canopy depths, continuous variables that indicate the thickness of the canopy along the vector formed by θ and directly above each sampled leaf, respectively (Fig. 2). Essentially, VD is the quantitative analog of the coarser distinctions of “upper”/“lower” or “shaded”/“unshaded” canopy leaves, and measures the vertical distance from the canopy exterior to a leaf sample. The expressions for RD and VD were derived by approximating the exterior of the tree canopy as a hemisphere with an arc of radius 100%ρ. RD was calculated as the direct distance from the canopy exterior to the leaf along θ as:
For VD a pair of overlapping right triangles with bases of identical lengths and hypotenuses of 100%ρ and %ρ were used to determine the total thickness of the canopy and the thickness below the leaf position, respectively. VD was calculated as the difference between these thicknesses using trigonometric identities and the values of 100%ρ, RD, and θ as follows:
Vertical canopy depths ranged from a minimum value of 0%ρ (for leaves located on the exterior of the canopy) to a maximum of 100%ρ (for leaves at the base of the canopy with θ = 0°).
Multiple linear regressions were used to determine the effects of VD, RD, θ, and ϕ on PAR and A300. Fixed effects of RD, θ, ϕ, VD, as well as interactions between the spatial variables (RD, θ, ϕ), were examined with tree identity treated as a random effect. Analyses of residuals revealed that model errors were normally distributed and variances were homogenous. All variance inflation factors were less than 1.08, indicating that correlations between variables defining leaf location (RD, θ, ϕ) were weak. The relationships between the predictors and the response variables in our data set were linear; other statistical analyses could be adopted if the relationships between variables are nonlinear.
Analyses were performed using JMP version 9.0 (SAS Institute, Cary, North Carolina, USA). Results are considered to be significant when P < 0.05. Threedimensional spatial visualization was performed using MATLAB version 2015a (MathWorks, Natick, Massachusetts, USA). Spherical polar coordinates were converted to a Cartesian coordinate system for ease of visualization using the following formulas:
PAR and A300 were plotted as color gradient scatter plots on a unitless Cartesian coordinate scale.
RESULTS
Spatial distribution of PAR—The amount of PAR received by the leaves ranged from 76.5 to 1718.7 µmol m−2 s−1 and depended on the leaf position (Table 1, Fig. 3A). Overall, leaves on the exterior of the canopy received more PAR than leaves on the interior. Although PAR was negatively correlated with RD (β= −10.6088, P < 0.0001), it was not significantly correlated with the other spatial variables θ and ϕ or any interaction term combinations of the spatial variables (Table 1). The strong dependence of PAR on RD is consistent with an internal shading model within the canopy, where the pigments within outer leaves reduce the amount of radiation that reaches the interior leaves.
Spatial distribution of UV-absorbing compound concentrations (A300)—The concentration of UV-absorbing compounds varied among leaves throughout the canopy (i.e., individual A300 values ranged from 0.301 to 1.089). This variation was explained in part by the position of the leaf within the canopy (Table 1, Fig. 3B). Much like the observed spatial trends in PAR, A300 was negatively correlated with RD; leaves at the exterior of the tree canopy contained higher concentrations of UV-absorbing compounds than leaves on the interior of the canopy (β = −0.0041, P < 0.0001). Also similar to PAR, A300 was not significantly correlated with any of the other spatial variables or the interaction terms.
Vertical canopy depth (VD)—Both PAR and A300 were negatively correlated with VD (P < 0.0001; Table 2). To understand this observation, we examined the independent spatial variables that can contribute to VD (i.e., RD and θ). Although PAR and A300 were both correlated with RD, neither showed a significant correlation with θ or the RD . θ interaction term in the linear regression analysis, indicating that the effect of RD on PAR and A300 was independent of θ. Consequently, the observed correlation of PAR and A300 with VD was driven by RD and was not significantly influenced by θ.
Table 1.
Results of multiple linear regression models showing the effects of independent variables on PAR and A300 in Avicennia germinans.
DISCUSSION
We have developed a simple, economical, and portable method based on a spherical polar coordinate system to quantify the exact location of leaves within the canopy. By defining leaf position with three continuous variables in spherical coordinates, our three-dimensional spatial mapping technique can detect fine-scale heterogeneity in both environmental conditions (i.e., PAR) and biological traits (i.e., A300) throughout the entire canopy of individual trees. Data collected using this method can be analyzed using multiple linear regressions to determine the degree to which each spatial variable (i.e., RD, θ, and ϕ) and combination of spatial variables (i.e., their interaction terms) explain the variation observed.
In our application of this method, we found that the spatial variation in both PAR and A300 within the canopy of black mangrove (A. germinans) trees was due to RD. The shared RDdependence of PAR and A300 suggests that the spatial distribution of light within canopies could determine the spatial distribution of UV-absorbing compounds. Multiple studies have observed increased production of photoprotective flavonoids and HCAs upon increased light exposure for plant species of diverse geographical origins (McKee, 1995; Schnitzler et al., 1996; Tegelberg et al., 2001; Izaguirre et al., 2007; Sullivan et al., 2007; El Morchid et al., 2014; Nascimento et al., 2015). In a manipulative growth chamber experiment, black mangrove seedlings grown in high light (339 ± 4 µmol m−2 s−1) had 52% more phenolics than those grown in low light (37 ± 1 µmol m−2 s−1) (McKee, 1995). The observed significant negative correlations of PAR and A300 with RD in our study suggest that the light microenvironment and the biosynthesis of UV-absorbing compounds are also coupled in adult trees in the field. Moreover, these results provide a plausible framework for explaining how the canopy distributes UV-absorbing compounds. Rather than produce a uniform concentration of UV-absorbing compounds across the entirety of the canopy, black mangroves preferentially allocate more UV-absorbing compounds to the outermost leaves, which receive more incoming solar radiation.
Other methods have been used to quantify the fine-scale environmental and biological heterogeneity within canopies. For example, leaf position has been quantified using leaf area index (Hirose and Werger, 1987; Pierce and Running, 1988; Ellsworth and Reich, 1993; Bréda, 2003; Weiss et al., 2004; Li et al., 2014; Magney et al., 2016) and vertical canopy depth (Ellsworth and Reich, 1993; Sinoquet et al., 1998; Magney et al., 2016). However, leaf area index methods often require either direct harvesting methods (i.e., stratified clipping or leaf collection within traps), which can be destructive, or indirect methods based on radiative transfer theory, which require assumptions about light propagation within the canopy (Bréda, 2003; Jonckheere et al., 2004; Weiss et al., 2004). Vertical canopy depth, in turn, assumes a vertically stratified arrangement of the canopy, thus ignoring the three-dimensional nature of the canopy structure. Digitizing techniques, such as LiDAR, have shown exceptional promise as universal approaches for quantifying both environmental and biological variation within canopies at resolutions below the scale of individual leaves (Greaves et al., 2015; Magney et al., 2016; Swatantran et al., 2016). Recent LiDAR studies have provided new insights into canopy heterogeneity by coupling a raytracing algorithm with spatial coordinates to determine path length (Magney et al., 2016). This method, similar to our own calculation of RD, requires determining the distance from the canopy exterior to a fixed point within the canopy interior. Our three-dimensional mapping approach, however, provides a less expensive and more portable alternative, capable of resolving canopy heterogeneity at a leaf-scale resolution. In addition, our method offers the ability to quantify leaf position within the interior of the canopy, a functionality not possible via LiDAR.
Table 2.
Results of linear regression models showing the effects of vertical depth (VD) on PAR and A300 in Avicennia germinans. VD is the quantitative analog of existing coarse descriptions of “upper” and “lower” canopies.
Although our three-dimensional mapping approach is inexpensive, easy to implement, and can provide detailed quantitative data regarding environmental and biological heterogeneity throughout the entirety of the canopy, this method is inherently built upon several assumptions. In particular, the method treats the trunk of the tree as a volume-less origin within the coordinate system. Although this treatment does not influence the quantification of ϕ, systematic errors in ρ and θ can occur. In addition, both the calculations of RD and VD rely on the assumption that the canopy shape is hemispherical, ignoring potential asymmetries or gaps within the canopy. Although RD and VD can be measured directly in the field, the determination of absolute depths would substantially increase the collection time of individual leaf parameters. The mathematical calculation of RD and VD, based on a hemispherical canopy model, can be executed rapidly across a large data set. In practice, the spherical polar coordinate system is best suited for small to medium hemispherical canopies that are easily accessible. However, for larger canopies or taller trees, the protractors could be repositioned to multiple locations within the canopy. Sampling data sets could then be mathematically merged using measured distances between sampling origin points. In addition, while our approach does not provide information regarding individual leaf angles, it can easily be coupled with similar inexpensive/portable methods capable of quantifying leaf angle (Escribano-Rocafort et al., 2014) or used in combination with other sampling techniques.
Overall, we have demonstrated that three-dimensional mapping can be used to probe the structure-property relationships of both environmental (i.e., PAR) and biological (i.e., A300) measurements. The significant negative correlation of RD with both PAR and A300 in our study provides vital clues regarding how the canopy structure efficiently balances the need for photoprotection with the corresponding energetic cost of biosynthesis. Using the spherical polar coordinate system to quantify leaf position, however, is not limited to studying the relationship between light and pigments. This simple, inexpensive, and portable method can also be used to more thoroughly understand how plants are affected by and respond to any abiotic or biotic factor, including quantifying how intraplant variation in leaf or floral traits is linked to spatial patterns of herbivory, the colonization and proliferation of microbes, and the effectiveness of pollinators.
ACKNOWLEDGMENTS
The authors thank V. Bernhard and his family for supporting undergraduate research and for allowing access to the collection sites, Shane Heschel for providing statistical advice, and Rachel Cochrane for providing an artistic rendering of the canopy. Funding was provided by the Hummingbird Cay Foundation, the Honors Program, the Department of Organismal and Environmental Biology, and the Department of Molecular Biology and Chemistry at Christopher Newport University.