Variable responses within epiphytic and benthic microalgal communities to nutrient enrichment

We examined the spatial extent of nitrogen (N) and phosphorus (P) limitation of each of the major benthic primary producer groups in Florida Bay (seagrass, epiphytes, macroalgae, and benthic microalgae) and characterized the shifts in primary producer community composition following nutrient enrichment. We established 24 permanent 0.25-m2 study plots at each of six sites across Florida Bay and added N and P to the sediments in a factorial design for 18 mo. Tissue nutrient content of the turtlegrass Thalassia testudinum revealed a spatial pattern in P limitation, from severe limitation in the eastern bay (N:P > 96:1), moderate limitation in two intermediate sites (approximately 63:1), and balanced with N availability in the western bay (approximately 31:1). P addition increased T. testudinum cover by 50-75% and short-shoot productivity by up to 100%, but only at the severely P-limited sites. At sites with an ambient N:P ratio suggesting moderate P limitation, few seagrass responses to nutrients occurred. Where ambient T. testudinum tissue N:P ratios indicated N and P availability was balanced, seagrass was not affected by nutrient addition but was strongly influenced by disturbance (currents, erosion). Macroalgal and epiphytic and benthic microalgal biomass were variable between sites and treatments. In general, there was no algal overgrowth of the seagrass in enriched conditions, possibly due to the strength of seasonal influences on algal biomass or regulation by grazers. N addition had little effect on any benthic primary producers throughout the bay. The Florida Bay benthic primary producer community was P limited, but P-induced alterations of community structure were not uniform among primary producers or across Florida Bay and did not always agree with expected patterns of nutrient limitation based on stoichiometric predictions from field assays of T. testudinum tissue N:P ratios.

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.

Forecasting Responses of Seagrass Distributions to Changing Water Quality Using Monitoring Data

Extensive data sets on water quality and seagrass distributions in Florida Bay have been assembled under complementary, but independent, monitoring programs. This paper presents the landscape-scale results from these monitoring programs and outlines a method for exploring the relationships between two such data sets. Seagrass species occurrence and abundance data were used to define eight benthic habitat classes from 677 sampling locations in Florida Bay. Water quality data from 28 monitoring stations spread across the Bay were used to construct a discriminant function model that assigned a probability of a given benthic habitat class occurring for a given combination of water quality variables. Mean salinity, salinity variability, the amount of light reaching the benthos, sediment depth, and mean nutrient concentrations were important predictor variables in the discriminant function model. Using a cross-validated classification scheme, this discriminant function identified the most likely benthic habitat type as the actual habitat type in most cases. The model predicted that the distribution of benthic habitat types in Florida Bay would likely change if water quality and water delivery were changed by human engineering of freshwater discharge from the Everglades. Specifically, an increase in the seasonal delivery of freshwater to Florida Bay should cause an expansion of seagrass beds dominated by Ruppia maritima and Halodule wrightii at the expense of the Thalassia testudinum-dominated community that now occurs in northeast Florida Bay. These statistical techniques should prove useful for predicting landscape-scale changes in community composition in diverse systems where communities are in quasi-equilibrium with environmental drivers.

Nutrient content of the seagrass Thalassia testudinum reveals regional patterns of relative availability of nitrogen and phosphorus in the Florida Keys USA

Between 1992 and 2000, we sampled 504 randomly chosen locations in theFlorida Keys, Florida, USA, for the elemental content of green leaves of theseagrass Thalassia testudinum. Carbon content ranged from29.4–43.3% (dry weight), nitrogen content from 0.88–3.96%, andphosphorus content from 0.048–0.243%. N and P content of the samples werenot correlated, suggesting that the relative availability of N and P variedacross the sampling region. Spatial pattern in C:N indicated a decrease in Navailability from inshore waters to the reef tract 10 km offshore;in contrast, the pattern in C:P indicated an increase in P availability frominshore waters to the reef tract. The spatial pattern in N:P was used to definea P-limited region of seagrass beds in Florida Bay and near shore, and anN-limited region of seagrass beds offshore. The close juxtaposition ofN–and P-limited regions allows the possibility that N loading from thesuburban Florida Keys could influence the offshore, N-limited seagrass bedswithout impacting the more nearshore, P-limited seagrass beds. Carbonate - Nutrient lim

Spatiotemporal patterns in community structure of macroinvertebrates inhabiting calcareous periphyton mats

Calcareous floating periphyton mats in the southern Everglades provide habitat for a diverse macroinvertebrate community that has not been well characterized. Our study described this community in an oligotrophic marsh, compared it with the macroinvertebrate community associated with adjacent epiphytic algae attached to macrophytes in the water column, and detected spatial patterns in density and community structure. The floating periphyton mat (floating mat) and epiphytic algae in the water column (submerged epiphyton) were sampled at 4 sites (1 km apart) in northern Shark River Slough, Everglades National Park (ENP), in the early (July) and late (November) wet season. Two perpendicular 90-m transects were established at each site and 100 samples were taken in a nested design. Sites were located in wet-prairie spikerush-dominated sloughs with similar water depths and emergent macrophyte communities. Floating mats were sampled by taking cores (6-cm diameter) that were sorted under magnification to enumerate infauna retained on a 250-μm-mesh sieve and with a maximum dimension >1 mm. Our results showed that floating mats provide habitat for a macroinvertebrate community with higher densities (no. animals/g ash-free dry mass) of Hyalella azteca, Dasyhelea spp., and Cladocera, and lower densities of Chironomidae and Planorbella spp. than communities associated with submerged epiphyton. Densities of the most common taxa increased 3× to 15× from early to late wet season, and community differences between the 2 habitat types became more pronounced. Floating-mat coverage and estimated floating-mat biomass increased 20 to 30% and 30 to 110%, respectively, at most sites in the late wet season. Some intersite variation was observed in individual taxa, but no consistent spatial pattern in any taxon was detected at any scale (from 0.2 m to 3 km). Floating mats and their resident macroinvertebrate communities are important components in the Everglades food web. This community should be included in environmental monitoring programs because degradation and eventual loss of the calcareous periphyton mat is associated with P enrichment in this ecosystem.

Drivers of Decadal-Scale Change in Southern Everglades Wetland Macrophyte Communities of the Coastal Ecotone

Understanding the relationships between hydrology and salinity and plant community structure and production is critical to allow predictions of wetland responses to altered water management, changing precipitation patterns and rising sea-level. We addressed how salinity, water depth, hydroperiod, canal inflows, and local precipitation control marsh macrophyte aboveground net primary production (ANPP) and structure in the coastal ecotone of the southern Everglades. We contrasted responses in two watersheds - Taylor Slough (TS) and C-111 - systems that have and will continue to experience changes in water management. Based on long-term trajectories in plant responses, we found continued evidence of increasing water levels and length of inundation in the C-111 watershed south of the C-111 canal. We also found strong differentiation among sites in upper TS that was dependent on hydrology. Finally, salinity, local precipitation and freshwater discharge from upstream explained over 80 % of the variance in Cladium ANPP at a brackish water site in TS. Moreover, our study showed that, while highly managed, the TS and C-111 watersheds maintain legacies in spatial pattern that would facilitate hydrologic restoration. Based on the trajectories in Cladium and Eleocharis, shifts in plant community structure could occur within 5–10 years of sustained water management change.

Conservação genética da orquídea Cattleya granulosa Lindley

Cattleya granulosa Lind is a large and endemic orchid in Atlantic Forest fragments in Northeast Brazil. The facility of collecting, uniqueness of their flowers, which have varying colors between green and reddish brown, and distribution in coastal areas of economic interest make their populations a constant target of predation, which also suffer from environmental degradation. Due to the impact on their populations, the species is threatened. In this study, we evaluate the levels of spatial aggregation in a preserved population, analyze the phylogenetic relationships of C. granulosa Lindl. with four other Laeliinae species (Brassavola tuberculata, C. bicolor, C. labiata and C. schofieldiana) and also to evaluate the genetic diversity of 12 remaining populations of C. granulosa Lindl. through ISSR. There was specificity of epiphytic C. granula Lindl. with a single host tree, species of Eugenia sp. C. granulosa Lindl. own spatial pattern, with the highest density of neighbors within up to 5 m. Regarding the phylogenetic relationships and genetic patterns with other species of the genus, C. bicolor exhibited the greatest genetic diversity (HE = 0.219), while C. labiata exhibited the lowest level (HE = 0.132). The percentage of genetic variation among species (AMOVA) was 23.26%. The principal component analysis (PCA) of ISSR data showed that unifoliate and bifoliolate species are genetically divergent. PCA indicated a close relationship between C. granulosa Lindl. and C. schofieldiana, a species considered to be a variety of C. granulosa Lindl. by many researchers. Population genetic analysis using ISSR showed all polymorphic loci. The high genetic differentiation between populations (ФST = 0.391, P < 0.0001) determined the structure into nine groups according to log-likelihood of Bayesian analysis, with a similar pattern in the dendrogram (UPGMA) and PCA. A positive and significant correlation between geographic and genetic distances between populations was identified (r = 0.794, P = 0.017), indicating isolation by distance. Patterns of allelic diversity suggest the occurrence of population bottlenecks in most populations of C. granulosa Lindl. (n = 8). Genetic data indicate that enable the maintenance of genetic diversity of the species is complex and is directly related to the conservation of different units or groups that are spatially distant.

(Table 2) Radiocarbon dates from lakes on Fildes Peninsula, King George Island

Precise relative sea level (RSL) data are important for inferring regional ice sheet histories, as well as helping to validate numerical models of ice sheet evolution and glacial isostatic adjustment. Here we develop a new RSL curve for Fildes Peninsula, South Shetland Islands (SSIs), a sub-Antarctic archipelago peripheral to the northern Antarctic Peninsula ice sheet, by integrating sedimentary evidence from isolation basins with geomorphological evidence from raised beaches. This combined approach yields not only a Holocene RSL curve, but also the spatial pattern of how RSL change varied across the archipelago. The curve shows a mid-Holocene RSL highstand on Fildes Peninsula at 15.5 m above mean sea level between 8000 and 7000 cal a BP. Subsequently RSL gradually fell as a consequence of isostatic uplift in response to regional deglaciation. We propose that isostatic uplift occurred at a non-steady rate, with a temporary pause in ice retreat ca. 7200 cal a BP, leading to a short-lived RSL rise of ~1 m and forming a second peak to the mid-Holocene highstand. Two independent approaches were taken to constrain the long-term tectonic uplift rate of the SSIs at 0.22-0.48 m/ka, placing the tectonic contribution to the reconstructed RSL highstand between 1.4 and 2.9 m. Finally, we make comparisons to predictions from three global sea level models.

(Table 1) Core top Globorotalia truncatulinoides d18O ratios of Atlantic Ocean samples

We investigate whether the oxygen isotope ratio in the test of Globorotalia truncatulinoides can serve as a proxy for intermediate depth (200-500 m) density. Since intermediate depth horizontal density gradients are associated with the vertical shear of upper ocean flows, this proxy could provide a tool for reconstructing past ocean circulation. The spatial pattern of core top Gr. truncatulinoides d18O in the Atlantic Ocean mimics the upper ocean density gradients associated with the major ocean currents. To better constrain the controls on the calcification depth(s) of Gr. truncatulinoides, we attempt to simulate the surface sediment data set using water column temperature and salinity conditions above the core sites. We predicted foraminiferal d18O for each core site assuming (1) the calcification occurs at a single depth and (2) the initial calcification is at the surface and the subsequent calcification is at 800 m water depth. The predicted d18O best resembles measured d18O of Gr. truncatulinoides when using (1) a single depth calcification at 350 m or (2) a two-depth approximation with 30% surface and 70% 800-m calcification. This result gives us confidence in the ability of d18O in Gr. truncatulinoides to proxy lateral density gradients at the intermediate depths associated with upper ocean flow.

Biogenic opal, biogenic barium, calcium carbonate and U/Th isotope data from core-top sediments of SONNE cruise SO202 (INOPEX)

Sedimentary proxies used to reconstruct marine productivity suffer from variable preservation and are sensitive to factors other than productivity. Therefore, proxy calibration is warranted. Here we map the spatial patterns of two paleoproductivity proxies, biogenic opal and barium fluxes, from a set of core-top sediments recovered in the Subarctic North Pacific. Comparisons of the proxy data with independent estimates of primary and export production, surface water macronutrient concentrations and biological pCO2 drawdown indicate that neither proxy shows a significant correlation with primary or export productivity for the entire region. Biogenic opal fluxes, when corrected for preservation using 230Th-normalized accumulation rates, show a good correlation with primary productivity along the volcanic arcs (tau = 0.71, p = 0.0024) and with export productivity throughout the western Subarctic North Pacific (tau = 0.71, p = 0.0107). Moderate and good correlations of biogenic barium flux with export production (tau = 0.57, p = 0.0022) and with surface water silicate concentrations (tau = 0.70, p = 0.0002) are observed for the central and eastern Subarctic North Pacific. For reasons unknown, however, no correlation is found in the western Subarctic North Pacific between biogenic barium flux and the reference data. Nonetheless, we show that barite saturation, uncertainty in the lithogenic barium corrections and problems with the reference datasets are not responsible for the lack of a significant correlation between biogenic barium flux and the reference data. Further studies evaluating the factors controlling the variability of the biogenic constituents in the sediments are desirable in this region.

Juniperus navicularis - a study of Portuguese population structure and fitness

Juniperus navicularis Gand. is a dioecious endemic conifer that constitutes the understory of seaside pine forests in Portugal, areas currently threatened by increasing urban expansion. The aim of this study is to assess the conservation status of previously known populations of this species located on its core area of distribution. The study was performed in south-west coast of Portugal. Three populations varying in size and pine density were analyzed. Number of individuals, population density, spatial distribution and individual characteristics of junipers were estimated. Female cone, seed characteristics and seed viability were also evaluated. Results suggest that J. navicularis populations are vulnerable because seminal recruitment is scarce, what may lead to a reduction of genetic variability due solely to vegetative propagation. This vulnerability seems to be strongly determined by climatic constraints toward increasing aridity. Ratio between male and female shrubs did not differ from 1:1 in any population. Deviations from 1:1 between mature and non-mature plants were found in all populations, denoting population ageing. Very low seed viability was observed. A major part of described Juniperus navicularis populations have disappeared through direct habitat loss to urban development, loss of fitness in drier and warmer locations and low seed viability. This study is the first to address J. navicularis conservation, and represents a valuable first step toward this species preservation. 

Lacustrine Arcellinina (Testate Amoebae) as Bioindicators of Arsenic Contamination

Arcellininids (testate amoebae) were examined from 61 surface sediment samples collected from 59 lakes in the vicinity of former gold mines, notably Giant Mine, near Yellowknife, Northwest Territories, Canada to determine their utility as bioindicators of arsenic (As), which occurs both as a byproduct of gold extraction at mines in the area and ore-bearing outcrops. Cluster analysis (Q-R-mode) and detrended correspondence analysis (DCA) reveal five arcellininid assemblages, three of which are related to varying As concentrations in the sediment samples. Redundancy analysis (RDA) showed that 14 statistically significant environmental parameters explained 57 % of the variation in faunal distribution, while partial RDA indicated that As had the greatest influence on assemblage variance (10.7 %; p < 0.10). Stress-indicating species (primarily centropyxids) characterized the faunas of samples with high As concentrations (median = 121.7 ppm, max > 10000 ppm, min = 16.1 ppm, n = 32), while difflugiid dominated assemblages were prevalent in substrates with relatively low As concentrations (median = 30.2 ppm, max = 905.2 ppm, min = 6.3 ppm, n = 20). Most of the lakes with very high As levels are located downwind (N and W) of the former Giant Mine roaster stack where refractory ore was roasted and substantial quantities of As were released (as As2O3) to the atmosphere in the first decade of mining. This spatial pattern suggests that a significant proportion of the observed As, in at least these lakes, are industrially derived. The results of this study highlight the sensitivity of Arcellinina to As and confirm that the group has considerable potential for assessing the impact of As contamination on lakes.

DYNAMICS OF THE SHORT TERM AND LONG TERM AVIAN DISTRIBUTIONS IN NORTH AMERICA - THE ROLES OF VEGETATION HEIGHT HETEROGENEITY AND CLIMATIC FACTORS

Understanding how biodiversity spatially distribute over both the short term and long term, and what factors are affecting the distribution, are critical for modeling the spatial pattern of biodiversity as well as for promoting effective conservation planning and practices. This dissertation aims to examine factors that influence short-term and long-term avian distribution from the geographical sciences perspective. The research develops landscape level habitat metrics to characterize forest height heterogeneity and examines their efficacies in modelling avian richness at the continental scale. Two types of novel vegetation-height-structured habitat metrics are created based on second order texture algorithms and the concepts of patch-based habitat metrics. I correlate the height-structured metrics with the richness of different forest guilds, and also examine their efficacies in multivariate richness models. The results suggest that height heterogeneity, beyond canopy height alone, supplements habitat characterization and richness models of two forest bird guilds. The metrics and models derived in this study demonstrate practical examples of utilizing three-dimensional vegetation data for improved characterization of spatial patterns in species richness. The second and the third projects focus on analyzing centroids of avian distributions, and testing hypotheses regarding the direction and speed of these shifts. I first showcase the usefulness of centroids analysis for characterizing the distribution changes of a few case study species. Applying the centroid method on 57 permanent resident bird species, I show that multi-directional distribution shifts occurred in large number of studied species. I also demonstrate, plain birds are not shifting their distribution faster than mountain birds, contrary to the prediction based on climate change velocity hypothesis. By modelling the abundance change rate at regional level, I show that extreme climate events and precipitation measures associate closely with some of the long-term distribution shifts. This dissertation improves our understanding on bird habitat characterization for species richness modelling, and expands our knowledge on how avian populations shifted their ranges in North America responding to changing environments in the past four decades. The results provide an important scientific foundation for more accurate predictive species distribution modeling in future.

Assessing the influence of abiotic factors and leaf-level properties on the stability of growing-season canopy greenness in a deciduous forest

Maps depicting spatial pattern in the stability of summer greenness could advance understanding of how forest ecosystems will respond to global changes such as a longer growing season. Declining summer greenness, or “greendown”, is spectrally related to declining near-infrared reflectance and is observed in most remote sensing time series to begin shortly after peak greenness at the end of spring and extend until the beginning of leaf coloration in autumn,. Understanding spatial patterns in the strength of greendown has recently become possible with the advancement of Landsat phenology products, which show that greendown patterns vary at scales appropriate for linking these patterns to proposed environmental forcing factors. This study tested two non-mutually exclusive hypotheses for how leaf measurements and environmental factors correlate with greendown and decreasing NIR reflectance across sites. At the landscape scale, we used linear regression to test the effects of maximum greenness, elevation, slope, aspect, solar irradiance and canopy rugosity on greendown. Secondly, we used leaf chemical traits and reflectance observations to test the effect of nitrogen availability and intrinsic water use efficiency on leaf-level greendown, and landscape-level greendown measured from Landsat. The study was conducted using Quercus alba canopies across 21 sites of an eastern deciduous forest in North America between June and August 2014. Our linear model explained greendown variance with an R2=0.47 with maximum greenness as the greatest model effect. Subsequent models excluding one model effect revealed elevation and aspect were the two topographic factors that explained the greatest amount of greendown variance. Regression results also demonstrated important interactions between all three variables, with the greatest interaction showing that aspect had greater influence on greendown at sites with steeper slopes. Leaf-level reflectance was correlated with foliar δ13C (proxy for intrinsic water use efficiency), but foliar δ13C did not translate into correlations with landscape-level variation in greendown from Landsat. Therefore, we conclude that Landsat greendown is primarily indicative of landscape position, with a small effect of canopy structure, and no measureable effect of leaf reflectance. With this understanding of Landsat greendown we can better explain the effects of landscape factors on vegetation reflectance and perhaps on phenology, which would be very useful for studying phenology in the context of global climate change

Simulating the effects of vegetation and landscape structure on fire behavior in northeastern Portugal: the case of holm oak (Quercus rotundifolia)

O fogo é um processo frequente nas paisagens do norte de Portugal. Estudos anteriores mostraram que os bosques de azinheira (Quercus rotundifolia) persistem após a passagem do fogo e ajudam a diminuir a sua intensidade e taxa de propagação. Os principais objetivos deste estudo foram compreender e modelar o efeito dos bosques de azinheira no comportamento do fogo ao nível da paisagem da bacia superior do rio Sabor, localizado no nordeste de Portugal. O impacto dos bosques de azinheira no comportamento do fogo foi testado em termos de área e configuração de acordo com cenários que simulam a possível distribuição destas unidades de vegetação na paisagem, considerando uma percentagem de ocupação da azinheira de 2.2% (Low), 18.1% (Moderate), 26.0% (High), e 39.8% (Rivers). Estes cenários tiveram como principal objetivo testar 1) o papel dos bosques de azinheira no comportamento do fogo e 2) de que forma a configuração das manchas de azinheira podem ajudar a diminuir a intensidade da linha de fogo e área ardida. Na modelação do comportamento do fogo foi usado o modelo FlamMap para simular a intensidade de linha do fogo e taxa de propagação do fogo com base em modelos de combustível associados a cada ocupação e uso do solo presente na área de estudo, e também com base em fatores topográficos (altitude, declive e orientação da encosta) e climáticos (humidade e velocidade do vento). Foram ainda usados dois modelos de combustível para a ocupação de azinheira (áreas interiores e de bordadura), desenvolvidos com base em dados reais obtidos na região. Usou-se o software FRAGSATS para a análise dos padrões espaciais das classes de intensidade de linha do fogo, usando-se as métricas Class Area (CA), Number of Patches (NP) e Large Patches Index (LPI). Os resultados obtidos indicaram que a intensidade da linha de fogo e a taxa de propagação do fogo variou entre cenários e entre modelos de combustível para o azinhal. A intensidade média da linha de fogo e a taxa média de propagação do fogo decresceu à medida que a percentagem de área de bosques de azinheira aumentou na paisagem. Também foi observado que as métricas CA, NP e LPI variaram entre cenários e modelos de combustível para o azinhal, decrescendo quando a percentagem de área de bosques de azinheira aumentou. Este estudo permitiu concluir que a variação da percentagem de ocupação e configuração espacial dos bosques de azinheira influenciam o comportamento do fogo, reduzindo, em termos médios, a intensidade da linha de fogo e a taxa de propagação, sugerindo que os bosques de azinhal podem ser usados como medidas silvícolas preventivas para diminuir o risco de incêndio nesta região.