Organic matter dynamics in the Florida coastal Everglades: A molecular marker and isotopic approach

Everglades National Park (ENP) is about to undergo the world's largest wetland restoration with the aim of improving the quality, timing and distribution of water flow. The changes in water flow are hypothesized to alter the nutrient fluxes and organic matter (OM) dynamics within ENP, especially in the estuarine areas. This study used a multi-proxy approach of molecular markers and stable δ 13C isotope measurements, to determine the present day OM dynamics in ENP. ^ OM dynamics in wetland soils/sediments have proved to be difficult to understand using traditional geochemical approaches. These are often inadequate to describe the multitude of OM sources (e.g. higher land plant, emergent vegetation, submerged vegetation) to the soils/sediments and the complex diagenetic processes that can alter the OM characteristics. A multi-proxy approach, however, that incorporates both molecular level and bulk parameter information is ideal to comprehend complex OM dynamics in aquatic environments. Therefore, biomass-specific molecular markers or proxies can be useful in tracing the sources and processing of OM. This approach was used to examine the OM dynamics in the two major drainage basins, Shark River Slough and Taylor River Slough, of ENP. Freshwater to marine transects were sampled in both systems for soils/sediments and suspended particulate organic matter (SPOM) to be characterized through bulk OM analyses, lipid biomarker determinations (e.g. sterols, fatty acids, hydrocarbons and triterpenoids) and compound-specific stable carbon isotope (δ 13C) determinations. ^ One key accomplishment of the research was the assessment of a molecular marker proxy (Paq) to distinguish between emergent/higher plant vegetation from submerged vegetation within ENP. This proxy proved to be quite useful at tracing OM inputs to the soils/sediments of ENP. A second key accomplishment was the development of a 3-way model using vegetation specific molecular markers. This novel, descriptive model was successfully applied to the estuarine areas of Taylor and Shark River sloughs, providing clear evidence of mixing of freshwater, estuarine and marine derived OM in these areas. In addition, diagenetic transformations of OM in these estuaries were found to be quite different between Taylor and Shark Rivers, and are likely a result of OM quality and hydrological differences. ^

Phosphorus budgets in Everglades wetland ecosystems: the effects of hydrology and nutrient enrichment

The Florida Everglades is a naturally oligotrophic hydroscape that has experienced large changes in ecosystem structure and function as the result of increased anthropogenic phosphorus (P) loading and hydrologic changes. We present whole-ecosystem models of P cycling for Everglades wetlands with differing hydrology and P enrichment with the goal of synthesizing existing information into ecosystem P budgets. Budgets were developed for deeper water oligotrophic wet prairie/slough (‘Slough’), shallower water oligotrophic Cladium jamaicense (‘Cladium’), partially enriched C. jamaicense/Typha spp. mixture (‘Cladium/Typha’), and enriched Typha spp. (‘Typha’) marshes. The majority of ecosystem P was stored in the soil in all four ecosystem types, with the flocculent detrital organic matter (floc) layer at the bottom of the water column storing the next largest proportion of ecosystem P pools. However, most P cycling involved ecosystem components in the water column (periphyton, floc, and consumers) in deeper water, oligotrophic Slough marsh. Fluxes of P associated with macrophytes were more important in the shallower water, oligotrophic Cladium marsh. The two oligotrophic ecosystem types had similar total ecosystem P stocks and cycling rates, and low rates of P cycling associated with soils. Phosphorus flux rates cannot be estimated for ecosystem components residing in the water column in Cladium/Typha or Typha marshes due to insufficient data. Enrichment caused a large increase in the importance of macrophytes to P cycling in Everglades wetlands. The flux of P from soil to the water column, via roots to live aboveground tissues to macrophyte detritus, increased from 0.03 and 0.2 g P m−2 yr−1 in oligotrophic Slough and Cladium marsh, respectively, to 1.1 g P m−2 yr−1 in partially enriched Cladium/Typha, and 1.6 g P m−2 yr−1 in enriched Typha marsh. This macrophyte translocation P flux represents a large source of internal eutrophication to surface waters in P-enriched areas of the Everglades.

Relationships among dietary intake, dietary patterns, dietary quality and nutrient adequacy and substance use in Latina mother/daughter dyads

This study examined the relationships among dietary intake, substance use, socioeconomic and acculturation-related factors among Latinas in Miami-Dade County. Substance abuse is rising among Latinas. A fuller understanding of this problem is needed given the rise of the Hispanic population and the role of women in Latin society. A better understanding between substance use and dietary intake can guide nutrition interventions to reduce negative substance-related health consequences. A purposeful sample of 320 Latina mother/daughter dyads were recruited and interviewed face-to-face as part of the Latino Women's Study. Dietary intake was collected via a 24-hour recall and examined by (1) nutrient intake, (2) dietary patterns using cluster analysis, (3) quality of diet using the Healthy Eating Index (HEI) and (4) the Dietary Reference Intakes to determine nutrient adequacy. Substance use was measured with the Drug Use Frequency and the Healthy and Daily Living Form. Acculturation was measured with the Cultural Identity Scale. Three dietary patterns emerged based on the number of servings from the food groups established in MyPyramid. None were associated with substance use. Latinas who reported using cannabis, cocaine, sedatives without prescription and/or more than five alcoholic drinks on an occasion at least once a month during the previous twelve months had significantly lower HEI scores (64 vs. 60; F = 7.8, p = .005) and consumed fewer fruits (F = 16, p < .001) than non-users. Latinas classified as mothers whom reported consuming cannabis at least 1-7 times a week had significantly lower HEI scores (F = 4.23, p = .015, η2 = .027) than daughters with the same frequency of substance use. One dimension of acculturation, greater familiarity with Latin culture, was associated with good dietary quality (β = .142, p = .012) regardless of any type of substance used or income level. There was a high prevalence of inadequacy of folic acid intake (50-75%) regardless of substance use. Substance users consumed significantly more energy (1,798 vs. 1,615; p = .027) than non-users. Although effect sizes were small, associations between dietary intake and substance use among Latinas deserve further exploration while acknowledging the combined association with acculturation. ^

Sediment and Nutrient Deposition Associated with Hurricane Wilma in Mangroves of the Florida Coastal Everglades

The distribution of mangrove biomass and forest structure along Shark River estuary in the Florida Coastal Everglades (FCE) has been correlated with elevated total phosphorus concentration in soils thought to be associated with storm events. The passage of Hurricane Wilma across Shark River estuary in 2005 allowed us to quantify sediment deposition and nutrient inputs in FCE mangrove forests associated with this storm event and to evaluate whether these pulsing events are sufficient to regulate nutrient biogeochemistry in mangrove forests of south Florida. We sampled the spatial pattern of sediment deposits and their chemical properties in mangrove forests along FCE sites in December 2005 and October 2006. The thickness (0.5 to 4.5 cm) of hurricane sediment deposits decreased with distance inland at each site. Bulk density, organic matter content, total nitrogen (N) and phosphorus (P) concentrations, and inorganic and organic P pools of hurricane sediment deposits differed from surface (0–10 cm) mangrove soils at each site. Vertical accretion resulting from this hurricane event was eight to 17 times greater than the annual accretion rate (0.30± 0.03 cm year−1) averaged over the last 50 years. Total P inputs from storm-derived sediments were equivalent to twice the average surface soil nutrient P density (0.19 mg cm−3). In contrast, total N inputs contributed 0.8 times the average soil nutrient N density (2.8 mg cm−3). Allochthonous mineral inputs from Hurricane Wilma represent a significant source of sediment to soil vertical accretion rates and nutrient resources in mangroves of southwestern Everglades. The gradient in total P deposition to mangrove soils from west to east direction across the FCE associated with this storm event is particularly significant to forest development due to the P-limited condition of this carbonate ecosystem. This source of P may be an important adaptation of mangrove forests in the Caribbean region to projected impacts of sea-level rise.

Spatio-temporal patterns and nutrient status of macroalgae in a heavily managed region of Biscayne Bay, Florida, USA

The coastal bays of South Florida are located downstream of the Florida Everglades, where a comprehensive restoration plan will strongly impact the hydrology of the region. Submerged aquatic vegetation communities are common components of benthic habitats of Biscayne Bay, and will be directly affected by changes in water quality. This study explores community structure, spatio-temporal dynamics, and tissue nutrient content of macroalgae to detect and describe relationships with water quality. The macroalgal community responded to strong variability in salinity; three distinctive macroalgal assemblages were correlated with salinity as follows: (1) low-salinity, dominated by Chara hornemannii and a mix of filamentous algae; (2) brackish, dominated by Penicillus capitatus, Batophora oerstedii, and Acetabularia schenckii; and (3) marine, dominated by Halimeda incrassata and Anadyomene stellata. Tissue-nutrient content was variable in space and time but tissues at all sites had high nitrogen and N:P values, demonstrating high nitrogen availability and phosphorus limitation in this region. This study clearly shows that distinct macroalgal assemblages are related to specific water quality conditions, and that macroalgal assemblages can be used as community-level indicators within an adaptive management framework to evaluate performance and restoration impacts in Biscayne Bay and other regions where both freshwater and nutrient inputs are modified by water management decisions.

Carbon and nutrient storage in subtropical seagrass meadows: examples from Florida Bay and Shark Bay

Seagrass meadows in Florida Bay and Shark Bay, contain substantial stores of both organic carbon and nutrients. Soils from both systems are predominantly calcium carbonate, with an average of 82.1% CaCO3 in Florida Bay compared to 71.3% in Shark Bay. Soils from Shark Bay had, on average, 21% higher organic carbon content and 35% higher phosphorus content than Florida Bay. Further, soils from Shark Bay had lower mean dry bulk density (0.78 ± 0.01 g mL-1) than those from Florida Bay (0.84 ± 0.02 mg mL-1). The most hypersaline regions of both bays had higher organic carbon content in surficial soils. Profiles of organic carbon and phosphorus from Florida Bay indicate that this system has experienced an increase in P delivery and primary productivity over the last century; in contrast, decreasing organic carbon and phosphorus with depth in the soil profiles in Shark Bay point to a decrease in phosphorus delivery and primary productivity over the last 1000 y. The total ecosystem stocks of stored organic C in Florida Bay averages 163.5 MgCorg ha-1, lower than the average of 243.0 MgCorg ha-1 for Shark Bay; but these values place Shark and Florida Bays among the global hotspots for organic C storage in coastal ecosystems.

Herbivore and nutrient impact on primary producer assemblages in a tropical marine environment

Globally, human populations are increasing and coastal ecosystems are becoming increasingly impacted by anthropogenic stressors. As eutrophication and exploitation of coastal resources increases, primary producer response to these drivers becomes a key indicator of ecosystem stability. Despite the importance of monitoring primary producers such as seagrasses and macroalgae, detailed studies on the response of these benthic habitat components to drivers remain relatively sparse. Utilizing a multi-faceted examination of turtle-seagrass and sea urchin-macroalgae consumer and nutrient dynamics, I elucidate the impact of these drivers in Akumal, Quintana Roo, Mexico. In Yal Ku Lagoon, macroalgae bioindicators signified high nutrient availability, which is important for further studies, but did not consistently follow published trends reflecting decreased δ 15N content with distance from suspected source. In Akumal Bay, eutrophication and grazing by turtles and fishes combine to structure patches within the seagrass beds. Grazed seagrass patches had higher structural complexity and productivity than patches continually grazed by turtles and fishes. Results from this study indicate that patch abandonment may follow giving-up density theory, the first to be recorded in the marine environment. As Diadema antillarum populations recover after their massive mortality thirty years ago, the role these echinoids will have in reducing macroalgae cover and altering ecosystem state remains to be clear. Although Diadema antillarum densities within the coral reef ecosystem were comparable to other regions within the Caribbean, the echinoid population in Akumal Bay was an insufficient driver to prevent dominance of a turf-algal-sediment (TAS) state. After a four year study, declining coral cover coupled with increased algal cover suggests that the TAS-dominated state is likely to persist over time despite echinoid recovery. Studies on macroalgal diversity and nutrients within this same region of echinoids indicated diversity and nutrient content of macroalgae increased, which may further increase the persistence of the algal-dominated state. This study provides valuable insight into the variable effects of herbivores and nutrients on primary producers within a tropical coastal ecosystem. Results from this work challenge many of the currently accepted theories on primary producer response to nutrients and herbivory while providing a framework for further studies into these dynamics.

Changes in mass and nutrient content of wood during decomposition in a south Florida mangrove forest

1. Large pools of dead wood in mangrove forests following disturbances such as hurricanes may influence nutrient fluxes. We hypothesized that decomposition of wood of mangroves from Florida, USA (Avicennia germinans, Laguncularia racemosa and Rhizophora mangle), and the consequent nutrient dynamics, would depend on species, location in the forest relative to freshwater and marine influences and whether the wood was standing, lying on the sediment surface or buried. 2. Wood disks (8–10 cm diameter, 1 cm thick) from each species were set to decompose at sites along the Shark River, either buried in the sediment, on the soil surface or in the air (above both the soil surface and high tide elevation). 3. A simple exponential model described the decay of wood in the air, and neither species nor site had any effect on the decay coefficient during the first 13 months of decomposition. 4. Over 28 months of decomposition, buried and surface disks decomposed following a two-component model, with labile and refractory components. Avicennia germinans had the largest labile component (18 ± 2% of dry weight), while Laguncularia racemosa had the lowest (10 ± 2%). Labile components decayed at rates of 0.37–23.71% month−1, while refractory components decayed at rates of 0.001–0.033% month−1. Disks decomposing on the soil surface had higher decay rates than buried disks, but both were higher than disks in the air. All species had similar decay rates of the labile and refractory components, but A. germinans exhibited faster overall decay because of a higher proportion of labile components. 5. Nitrogen content generally increased in buried and surface disks, but there was little change in N content of disks in the air over the 2-year study. Between 17% and 68% of total phosphorus in wood leached out during the first 2 months of decomposition, with buried disks having the greater losses, P remaining constant or increasing slightly thereafter. 6. Newly deposited wood from living trees was a short-term source of N for the ecosystem but, by the end of 2 years, had become a net sink. Wood, however, remained a source of P for the ecosystem. 7. As in other forested ecosystems, coarse woody debris can have a significant impact on carbon and nutrient dynamics in mangrove forests. The prevalence of disturbances, such as hurricanes, that can deposit large amounts of wood on the forest floor accentuates the importance of downed wood in these forests.

Distribution of Diatoms and Development of Diatom-Based Models for Inferring Salinity and Nutrient Concentrations in Florida Bay and Adjacent Coastal Wetlands of South Florida (USA)

The composition and distribution of diatom algae inhabiting estuaries and coasts of the subtropical Americas are poorly documented, especially relative to the central role diatoms play in coastal food webs and to their potential utility as sentinels of environmental change in these threatened ecosystems. Here, we document the distribution of diatoms among the diverse habitat types and long environmental gradients represented by the shallow topographic relief of the South Florida, USA, coastline. A total of 592 species were encountered from 38 freshwater, mangrove, and marine locations in the Everglades wetland and Florida Bay during two seasonal collections, with the highest diversity occurring at sites of high salinity and low water column organic carbon concentration (WTOC). Freshwater, mangrove, and estuarine assemblages were compositionally distinct, but seasonal differences were only detected in mangrove and estuarine sites where solute concentration differed greatly between wet and dry seasons. Epiphytic, planktonic, and sediment assemblages were compositionally similar, implying a high degree of mixing along the shallow, tidal, and storm-prone coast. The relationships between diatom taxa and salinity, water total phosphorus (WTP), water total nitrogen (WTN), and WTOC concentrations were determined and incorporated into weighted averaging partial least squares regression models. Salinity was the most influential variable, resulting in a highly predictive model (r apparent 2  = 0.97, r jackknife 2  = 0.95) that can be used in the future to infer changes in coastal freshwater delivery or sea-level rise in South Florida and compositionally similar environments. Models predicting WTN (r apparent 2  = 0.75, r jackknife 2  = 0.46), WTP (r apparent 2  = 0.75, r jackknife 2  = 0.49), and WTOC (r apparent 2  = 0.79, r jackknife 2  = 0.57) were also strong, suggesting that diatoms can provide reliable inferences of changes in solute delivery to the coastal ecosystem.

A preliminary analysis of the correlation of food-web characteristics with hydrology and nutrient gradients in the southern Everglades

We estimated trophic position and carbon source for three consumers (Florida gar, Lepisosteus platyrhincus; eastern mosquitofish, Gambusia holbrooki; and riverine grass shrimp, Palaemonetes paludosus) from 20 sites representing gradients of productivity and hydrological disturbance in the southern Florida Everglades, U.S.A. We characterized gross primary productivity at each site using light/dark bottle incubation and stem density of emergent vascular plants. We also documented nutrient availability as total phosphorus (TP) in floc and periphyton, and the density of small fishes. Hydrological disturbance was characterized as the time since a site was last dried and the average number of days per year the sites were inundated for the previous 10 years. Food-web attributes were estimated in both the wet and dry seasons by analysis of δ15N (trophic position) and δ13C (food-web carbon source) from 702 samples of aquatic consumers. An index of carbon source was derived from a two-member mixing model with Seminole ramshorn snails (Planorbella duryi) as a basal grazing consumer and scuds (amphipods Hyallela azteca) as a basal detritivore. Snails yielded carbon isotopic values similar to green algae and diatoms, while carbon values of scuds were similar to bulk periphyton and floc; carbon isotopic values of cyanobacteria were enriched in C13compared to all consumers examined. A carbon source similar to scuds dominated at all but one study site, and though the relative contribution of scud-like and snail-like carbon sources was variable, there was no evidence that these contributions were a function of abiotic factors or season. Gar consistently displayed the highest estimated trophic position of the consumers studied, with mosquitofish feeding at a slightly lower level, and grass shrimp feeding at the lowest level. Trophic position was not correlated with any nutrient or productivity parameter, but did increase for grass shrimp and mosquitofish as the time following droughts increased. Trophic position of Florida gar was positively correlated with emergent plant stem density.

Leaf Gas Exchange and Nutrient Use Efficiency Help Explain the Distribution of Two Neotropical Mangroves under Contrasting Flooding and Salinity

Rhizophora mangle and Laguncularia racemosa cooccur along many intertidal floodplains in the Neotropics. Their patterns of dominance shift along various gradients, coincident with salinity, soil fertility, and tidal flooding. We used leaf gas exchange metrics to investigate the strategies of these two species in mixed culture to simulate competition under different salinity concentrations and hydroperiods. Semidiurnal tidal and permanent flooding hydroperiods at two constant salinity regimes (10 g L−1 and 40 g L−1) were simulated over 10 months. Assimilation ( ), stomatal conductance ( ), intercellular CO2 concentration ( ), instantaneous photosynthetic water use efficiency (PWUE), and photosynthetic nitrogen use efficiency (PNUE) were determined at the leaf level for both species over two time periods. Rhizophora mangle had significantly higher PWUE than did L. racemosa seedlings at low salinities; however, L. racemosa had higher PNUE and and, accordingly, had greater intercellular CO2 (calculated) during measurements. Both species maintained similar capacities for A at 10 and 40 g L−1 salinity and during both permanent and tidal hydroperiod treatments. Hydroperiod alone had no detectable effect on leaf gas exchange. However, PWUE increased and PNUE decreased for both species at 40 g L−1 salinity compared to 10 g L−1. At 40 g L−1 salinity, PNUE was higher for L. racemosa than R. mangle with tidal flooding. These treatments indicated that salinity influences gas exchange efficiency, might affect how gases are apportioned intercellularly, and accentuates different strategies for distributing leaf nitrogen to photosynthesis for these two species while growing competitively.

Distribution of Diatoms Along Environmental Gradients in the Charlotte Harbor, Florida (USA), Estuary and Its Watershed: Implications for Bioassessment of Salinity and Nutrient Concentrations

The relative abundance of diatom species in different habitats can be used as a tool to infer prior environmental conditions and evaluate management decisions that influence habitat quality. Diatom distribution patterns were examined to characterize relationships between assemblage composition and environmental gradients in a subtropical estuarine watershed. We identified environmental correlates of diatom distribution patterns across the Charlotte Harbor, Florida, watershed; evaluated differences among three major river drainages; and determined how accurately local environmental conditions can be predicted using inference models based on diatom assemblages. Sampling locations ranged from freshwater to marine (0.1–37.2 ppt salinity) and spanned broad nutrient concentration gradients. Salinity was the predominant driver of difference among diatom assemblages across the watershed, but other environmental variables had stronger correlations with assemblages within the subregions of the three rivers and harbor. Eighteen indicator taxa were significantly affiliated with subregions. Relationships between diatom taxon distributions and salinity, distance from the harbor, total phosphorus (TP), and total nitrogen (TN) were evaluated to determine the utility of diatom assemblages to predict environmental values using a weighted averaging-regression approach. Diatom-based inferences of these variables were strong (salinity R 2 = 0.96; distance R 2 = 0.93; TN R 2 = 0.83; TP R 2 = 0.83). Diatom assemblages provide reliable estimates of environmental parameters on different spatial scales across the watershed. Because many coastal diatom taxa are ubiquitous, the diatom training sets provided here should enable diatom-based environmental reconstructions in subtropical estuaries that are being rapidly altered by land and water use changes and sea level rise.

A Suite of Prey Traits Determine Predator and Nutrient Enrichment Effects in a Tri-Trophic Food Chain

Predation, predation risk, and resource quality affect suites of prey traits that collectively impact individual fitness, population dynamics, and community structure. However, studies of multi-trophic level effects generally focus on a single prey trait, failing to capture trade-offs among suites of covarying traits that govern population responses and emergent community patterns. We used structural equation models (SEM) to summarize the non-lethal and lethal effects of crayfish, Procambarus fallax, and phosphorus (P) addition, which affected prey food quality (periphyton), on the interactive effects of behavioral, morphological, developmental, and reproductive traits of snails, Planorbella duryi. Univariate and multivariate analyses suggested trade-offs between production (growth, reproduction) and defense (foraging behavior, shell shape) traits of snails in response to non-lethal crayfish and P addition, but few lethal effects. SEM revealed that non-lethal crayfish effects indirectly limited per capita offspring standing stock by increasing refuge use, slowing individual growth, and inducing snails to produce thicker, compressed shells. The negative effects of non-lethal crayfish on snails were strongest with P addition; snails increased allocation to shell defense rather than growth or reproduction. However, compared to ambient conditions, P addition with non-lethal crayfish still yielded greater per capita offspring standing stock by speeding individual snail growth enabling them to produce more offspring that also grew faster. Increased refuge use in response to non-lethal crayfish led to a non-lethal trophic cascade that altered the spatial distribution of periphyton. Independent of crayfish effects, snails stimulated periphyton growth through nutrient regeneration. These findings illustrate the importance of studying suites of traits that reveal costs associated with inducing different traits and how expressing those traits impacts population and community level processes.

Movement Effects on the Flow Physics and Nutrient Delivery in Engineered Valvular Tissues

Mechanical conditioning has been shown to promote tissue formation in a wide variety of tissue engineering efforts. However the underlying mechanisms by which external mechanical stimuli regulate cells and tissues are not known. This is particularly relevant in the area of heart valve tissue engineering (HVTE) owing to the intense hemodynamic environments that surround native valves. Some studies suggest that oscillatory shear stress (OSS) caused by steady flow and scaffold flexure play a critical role in engineered tissue formation derived from bone marrow derived stem cells (BMSCs). In addition, scaffold flexure may enhance nutrient (e.g. oxygen, glucose) transport. In this study, we computationally quantified the i) magnitude of fluid-induced shear stresses; ii) the extent of temporal fluid oscillations in the flow field using the oscillatory shear index (OSI) parameter, and iii) glucose and oxygen mass transport profiles. Noting that sample cyclic flexure induces a high degree of oscillatory shear stress (OSS), we incorporated moving boundary computational fluid dynamic simulations of samples housed within a bioreactor to consider the effects of: 1) no flow, no flexure (control group), 2) steady flow-alone, 3) cyclic flexure-alone and 4) combined steady flow and cyclic flexure environments. We also coupled a diffusion and convention mass transport equation to the simulated system. We found that the coexistence of both OSS and appreciable shear stress magnitudes, described by the newly introduced parameter OSI-t , explained the high levels of engineered collagen previously observed from combining cyclic flexure and steady flow states. On the other hand, each of these metrics on its own showed no association. This finding suggests that cyclic flexure and steady flow synergistically promote engineered heart valve tissue production via OSS, so long as the oscillations are accompanied by a critical magnitude of shear stress. In addition, our simulations showed that mass transport of glucose and oxygen is enhanced by sample movement at low sample porosities, but did not play a role in highly porous scaffolds. Preliminary in-house in vitro experiments showed that cell proliferation and phenotype is enhanced in OSI-t environments.

Edward Castaneda retrieving a soil core to assess sedimentation and nutrient accumulation rates in a dwarf mangrove forest, Taylor Slough

This material is based upon work supported by the National Science Foundation through the Florida Coastal Everglades Long-Term Ecological Research program under Cooperative Agreements #DBI-0620409 and #DEB-9910514. This image is made available for non-commercial or educational use only.