Impact of Late Holocene climate variability and anthropogenic activities on Biscayne Bay (Florida, U.S.A.): Evidence from diatoms

Shallow marine ecosystems are experiencing significant environmental alterations as a result of changing climate and increasing human activities along coasts. Intensive urbanization of the southeast Florida coast and intensification of climate change over the last few centuries changed the character of coastal ecosystems in the semi-enclosed Biscayne Bay, Florida. In order to develop management policies for the Bay, it is vital to obtain reliable scientific evidence of past ecological conditions. The long-term records of subfossil diatoms obtained from No Name Bank and Featherbed Bank in the Central Biscayne Bay, and from the Card Sound Bank in the neighboring Card Sound, were used to study the magnitude of the environmental change caused by climate variability and water management over the last ~ 600 yr. Analyses of these records revealed that the major shifts in the diatom assemblage structures at No Name Bank occurred in 1956, at Featherbed Bank in 1966, and at Card Sound Bank in 1957. Smaller magnitude shifts were also recorded at Featherbed Bank in 1893, 1942, 1974 and 1983. Most of these changes coincided with severe drought periods that developed during the cold phases of El Niño Southern Oscillation (ENSO), Atlantic Multidecadal Oscillation (AMO) and Pacific Decadal Oscillation (PDO), or when AMO was in warm phase and PDO was in the cold phase. Only the 1983 change coincided with an unusually wet period that developed during the warm phases of ENSO and PDO. Quantitative reconstructions of salinity using the weighted averaging partial least squares (WA-PLS) diatom-based salinity model revealed a gradual increase in salinity at the three coring locations over the last ~ 600 yr, which was primarily caused by continuously rising sea level and in the last several decades also by the reduction of the amount of freshwater inflow from the mainland. Concentration of sediment total nitrogen (TN), total phosphorus (TP) and total organic carbon (TOC) increased in the second half of the 20th century, which coincided with the construction of canals, landfills, marinas and water treatment plants along the western margin of Biscayne Bay. Increased magnitude and rate of the diatom assemblage restructuring in the mid- and late-1900s, suggest that large environmental changes are occurring more rapidly now than in the past.

Impact of Late Holocene Climate Variability and Anthropogenic Activities on Biscayne Bay (Florida, U.S.A.): Evidence From Diatoms

Shallow marine ecosystems are experiencing significant environmental alterations as a result of changing climate and increasing human activities along coasts. Intensive urbanization of the southeast Florida coast and intensification of climate change over the last few centuries changed the character of coastal ecosystems in the semi-enclosed Biscayne Bay, Florida. In order to develop management policies for the Bay, it is vital to obtain reliable scientific evidence of past ecological conditions. The long-term records of subfossil diatoms obtained from No Name Bank and Featherbed Bank in the Central Biscayne Bay, and from the Card Sound Bank in the neighboring Card Sound, were used to study the magnitude of the environmental change caused by climate variability and water management over the last ~ 600 yr. Analyses of these records revealed that the major shifts in the diatom assemblage structures at No Name Bank occurred in 1956, at Featherbed Bank in 1966, and at Card Sound Bank in 1957. Smaller magnitude shifts were also recorded at Featherbed Bank in 1893, 1942, 1974 and 1983. Most of these changes coincided with severe drought periods that developed during the cold phases of El Niño Southern Oscillation (ENSO), Atlantic Multidecadal Oscillation (AMO) and Pacific Decadal Oscillation (PDO), or when AMO was in warm phase and PDO was in the cold phase. Only the 1983 change coincided with an unusually wet period that developed during the warm phases of ENSO and PDO. Quantitative reconstructions of salinity using the weighted averaging partial least squares (WA-PLS) diatom-based salinity model revealed a gradual increase in salinity at the three coring locations over the last ~ 600 yr, which was primarily caused by continuously rising sea level and in the last several decades also by the reduction of the amount of freshwater inflow from the mainland. Concentration of sediment total nitrogen (TN), total phosphorus (TP) and total organic carbon (TOC) increased in the second half of the 20th century, which coincided with the construction of canals, landfills, marinas and water treatment plants along the western margin of Biscayne Bay. Increased magnitude and rate of the diatom assemblage restructuring in the mid- and late-1900s, suggest that large environmental changes are occurring more rapidly now than in the past.

Present-day Pattern Variations in the Central and Eastern Pacific El Nino

El Niño and the Southern Oscillation (ENSO) is a cycle that is initiated in the equatorial Pacific Ocean and is recognized on interannual timescales by oscillating patterns in tropical Pacific sea surface temperatures (SST) and atmospheric circulations. Using correlation and regression analysis of datasets that include SST’s and other interdependent variables including precipitation, surface winds, sea level pressure, this research seeks to quantify recent changes in ENSO behavior. Specifically, the amplitude, frequency of occurrence, and spatial characteristics (i.e. events with maximum amplitude in the Central Pacific versus the Eastern Pacific) are investigated. The research is based on the question; “Are the statistics of ENSO changing due to increasing greenhouse gas concentrations?” Our hypothesis is that the present-day changes in amplitude, frequency, and spatial characteristics of ENSO are determined by the natural variability of the ocean-atmosphere climate system, not the observed changes in the radiative forcing due to change in the concentrations of greenhouse gases. Statistical analysis, including correlation and regression analysis, is performed on observational ocean and atmospheric datasets available from the National Oceanographic and Atmospheric Administration (NOAA), National Center for Atmospheric Research (NCAR) and coupled model simulations from the Coupled Model Inter-comparison Project (phase 5, CMIP5). Datasets are analyzed with a particular focus on ENSO over the last thirty years. Understanding the observed changes in the ENSO phenomenon over recent decades has a worldwide significance. ENSO is the largest climate signal on timescales of 2 - 7 years and affects billions of people via atmospheric teleconnections that originate in the tropical Pacific. These teleconnections explain why changes in ENSO can lead to climate variations in areas including North and South America, Asia, and Australia. For the United States, El Niño events are linked to decreased number of hurricanes in the Atlantic basin, reduction in precipitation in the Pacific Northwest, and increased precipitation throughout the southern United Stated during winter months. Understanding variability in the amplitude, frequency, and spatial characteristics of ENSO is crucial for decision makers who must adapt where regional ecology and agriculture are affected by ENSO.

Present-day Pattern Variations in the Central and Eastern Pacific El Nino

El Niño and the Southern Oscillation (ENSO) is a cycle that is initiated in the equatorial Pacific Ocean and is recognized on interannual timescales by oscillating patterns in tropical Pacific sea surface temperatures (SST) and atmospheric circulations. Using correlation and regression analysis of datasets that include SST’s and other interdependent variables including precipitation, surface winds, sea level pressure, this research seeks to quantify recent changes in ENSO behavior. Specifically, the amplitude, frequency of occurrence, and spatial characteristics (i.e. events with maximum amplitude in the Central Pacific versus the Eastern Pacific) are investigated. The research is based on the question; “Are the statistics of ENSO changing due to increasing greenhouse gas concentrations?” Our hypothesis is that the present-day changes in amplitude, frequency, and spatial characteristics of ENSO are determined by the natural variability of the ocean-atmosphere climate system, not the observed changes in the radiative forcing due to change in the concentrations of greenhouse gases. Statistical analysis, including correlation and regression analysis, is performed on observational ocean and atmospheric datasets available from the National Oceanographic and Atmospheric Administration (NOAA), National Center for Atmospheric Research (NCAR) and coupled model simulations from the Coupled Model Inter-comparison Project (phase 5, CMIP5). Datasets are analyzed with a particular focus on ENSO over the last thirty years. Understanding the observed changes in the ENSO phenomenon over recent decades has a worldwide significance. ENSO is the largest climate signal on timescales of 2 - 7 years and affects billions of people via atmospheric teleconnections that originate in the tropical Pacific. These teleconnections explain why changes in ENSO can lead to climate variations in areas including North and South America, Asia, and Australia. For the United States, El Niño events are linked to decreased number of hurricanes in the Atlantic basin, reduction in precipitation in the Pacific Northwest, and increased precipitation throughout the southern United Stated during winter months. Understanding variability in the amplitude, frequency, and spatial characteristics of ENSO is crucial for decision makers who must adapt where regional ecology and agriculture are affected by ENSO.

Individual foraminifera variability from the western and eastern equatorial Pacific during the Late Holocene and Last Glacial Maximum

El Niño-Southern Oscillation (ENSO) is a major source of global interannual variability, but its response to climate change is uncertain. Paleoclimate records from the Last Glacial Maximum (LGM) provide insight into ENSO behavior when global boundary conditions (ice sheet extent, atmospheric partial pressure of CO2) were different from those today. In this work, we reconstruct LGM temperature variability at equatorial Pacific sites using measurements of individual planktonic foraminifera shells. A deep equatorial thermocline altered the dynamics in the eastern equatorial cold tongue, resulting in reduced ENSO variability during the LGM compared to the Late Holocene. These results suggest that ENSO was not tied directly to the east-west temperature gradient, as previously suggested. Rather, the thermocline of the eastern equatorial Pacific played a decisive role in the ENSO response to LGM climate.

Age determination of North Atlantic sediment cores

Available overwash records from coastal barrier systems document significant variability in North Atlantic hurricane activity during the late Holocene. The same climate forcings that may have controlled cyclone activity over this interval (e.g., the West African Monsoon, El Niño-Southern Oscillation (ENSO)) show abrupt changes around 6000 yrs B.P., but most coastal sedimentary records do not span this time period. Establishing longer records is essential for understanding mid-Holocene patterns of storminess and their climatic drivers, which will lead to better forecasting of how climate change over the next century may affect tropical cyclone frequency and intensity. Storms are thought to be an important mechanism for transporting coarse sediment from shallow carbonate platforms to the deep-sea, and bank-edge sediments may offer an unexplored archive of long-term hurricane activity. Here, we develop this new approach, reconstructing more than 7000 years of North Atlantic hurricane variability using coarse-grained deposits in sediment cores from the leeward margin of the Great Bahama Bank. High energy event layers within the resulting archive are (1) broadly correlated throughout an offbank transect of multi-cores, (2) closely matched with historic hurricane events, and (3) synchronous with previous intervals of heightened North Atlantic hurricane activity in overwash reconstructions from Puerto Rico and elsewhere in the Bahamas. Lower storm frequency prior to 4400 yrs B.P. in our records suggests that precession and increased NH summer insolation may have greatly limited hurricane potential intensity, outweighing weakened ENSO and a stronger West African Monsoon-factors thought to be favorable for hurricane development.

Recent Changes in Central and Eastern Pacific El Niño

Recent research indicates that characteristics of El Niño and the Southern Oscillation (ENSO) have changed over the past several decades. Here, I examined different flavors of El Niño in the observational record and the recent changes in the character of El Niño events. The fundamental physical processes that drive ENSO were described and the Eastern Pacific (EP) and Central Pacific (CP) types or flavors of El Niño were defined. Using metrics from the peer-reviewed literature, I examined several historical data sets to interpret El Niño behavior from 1950-2010. A Monte Carlo Simulation was then applied to output from coupled model simulations to test the statistical significance of recent observations surrounding EP and CP El Niño. Results suggested that EP and CP El Niño had been occurring in a similar fashion over the past 60 years with natural variability, but no significant increase in CP El Niño behavior.

(Table 1) Age control points and of sediment core MD01-2444

The marine isotopic stage 3 (MIS3) at Ocean Drilling Program (ODP) Site 1060 (Gulf Stream) shows both sharp onset and end of interstadials, the existence of very short lived warm events during stadials, and points to differences in detail between the sea surface temperature (SST) record from the western North Atlantic and the atmospheric temperature record inferred from d18O in Greenland ice. Investigating MIS3 and obtaining comparable data from other locations appears crucial. The eastern Atlantic provides well-documented records of climate changes. We have selected a core from off Portugal and use it to examine Dansgaard/Oeschger events (D/O) at centennial-scale resolution (139 years on average between two data points). We have obtained a faunal data set for core MD01-2444, 37°N, 10°W, 2600 m water depth and use a group of species (Globigerina bulloides + Globigerinita glutinata) as a proxy of upwelling intensity driven by trade winds intensity changes. We tentatively relate the variation of this group to a North Atlantic Oscillation-like phenomenon (NAO) off Portugal. We observe that it resembles the rainfall index in the Caribbean as recorded at ODP Site 1002 (Cariaco Basin) which traces the Intertropical Convergence Zone (ITCZ) location through changes of terrigenous inputs. The driest intervals (ITZC to the south) at Site 1002 correspond to intervals of increased upwelling in MD01-2444 as well as the driest periods identified during stadials on similar cores in the area. Because the ITZC to the south is consistent with an El Niño-Southern Oscillation (ENSO+) situation, our study suggests a positive correlation between ENSO-like conditions and NAO-like conditions at a millennial timescale. During interstadial intervals when increased wetness over Cariaco is recorded (ITCZ to the north) and the upwelling in MD01-2444 is decreased, we see from both SSTs and faunal tropical indicators that MD01-2444 site is warm. In addition, interstadials are equally warm through each so-called Bond cycle. This contrasts with the Greenland Ice Core Project (GRIP) record where interstadial peaks are successively cooler through each Bond cycle. This record confirms a link between tropical climate linked to ITCZ position and the climate of southern Europe at millennial timescales, in spite of showing a very good correlation with polar latitudes (GRIP) through d18O on Globigerina bulloides. In addition, because the warmest SSTs and the d18O on G. bulloides are so remarkably different, our work points to changes in seasonality as a strong control over the climatic pattern of the North Atlantic area and the marked influence of winter conditions.

Stable oxygen isotope record of coral core RUS-95 from the Red Sea

A 245-year coral oxygen isotope record from the northern Red Sea (Ras Umm Sidd/Egypt, ~28°N) in bimonthly resolution is presented. The mean annual coral delta18O signal apparently reflects varying proportions of both sea surface temperature and delta18Oseawater variability. In conjunction with instrumental observations of climate the coral record suggests for interannual and longer timescales that colder periods are accompanied by more arid conditions in the northern Red Sea but increased rainfall in the southeastern Mediterranean, whereas warmer periods are accompanied by decreased rainfall in the latter and less arid conditions in the northern Red Sea. A ~70-year oscillation of probably North Atlantic origin dominates the coral time series. Interannual to interdecadal variability is correlated with instrumental indices of the North Atlantic Oscillation (NAO), the El Niño-Southern Oscillation (ENSO), and North Pacific climate variability. The results suggest that these modes contributed consistently to Middle East climate variability since at least 1750, preferentially at a period of ~5.7 years.

Oxygen isotopic measurements from Holcene of sediment core MD02-2529, tropical pacific

Changes in El Niño-Southern Oscillation (ENSO) variability are difficult to extract from paleoceanographic reconstructions because they are superimposed on changes in seasonal variability that modulate the first-order climate signal. Here we address this problem by reconstructing thermocline structure from a marine sediment core retrieved from the eastern equatorial Pacific. At the core location, changes in hydrologic parameters within the thermocline are linked to ENSO activity, with a reduced influence of seasonal variability compared to surface waters. We performed repeated isotopic analyses (d18O) on single specimens of the thermocline-dwelling planktonic foraminifera Neogloboquadrina dutertrei at several targeted time periods over the last 50 ka to extract the total thermocline variance, a parameter supposed to reveal changes in ENSO. No fundamental changes in amplitude and frequency of the events were detected despite differences in climatic background. However, our data suggest that long-term variations in the thermocline variability occurred over the last 50 ka, with the highest and lowest ENSO activities occurring during the last glacial period and the Last Glacial Maximum, respectively.

Planktic foraminifera abundances and winter SST reconstruction for sediment core SO90-39KG

The Indian winter monsoon (IWM) is a key component of the seasonally changing monsoon system that affects the densely populated regions of South Asia. Cold winds originating in high northern latitudes provide a link of continental-scale Northern Hemisphere climate to the tropics. Western Disturbances (WD) associated with the IWM play a critical role for the climate and hydrology in northern India and the western Himalaya region. It is vital to understand the mechanisms and teleconnections that influence IWM variability to better predict changes in future climate. Here we present a study of regionally calibrated winter (January) temperatures and according IWM intensities, based on a planktic foraminiferal record with biennial (2.55 years) resolution. Over the last ~250 years, IWM intensities gradually weakened, based on the long-term trend of reconstructed January temperatures. Furthermore, the results indicate that IWM is connected on interannual- to decadal time scales to climate variability of the tropical and extratropical Pacific, via El Niño Southern Oscillation (ENSO) and Pacific Decadal Oscillation (PDO). However, our findings suggest that this relationship appeared to begin to decouple since the beginning of the 20th century. Cross-spectral analysis revealed that several distinct decadal-scale phases of colder climate and accordingly more intense winter monsoon centered at the years ~1800, ~1890 and ~1930 can be linked to changes of the North Atlantic Oscillation (NAO).

(Table 2) Seasonal flux data and percentages of major bulk components of total flux at the mesotrophic sediment trap site CB from 1988 - 2012

A more than two-decadal sediment trap record from the Eastern Boundary Upwelling Ecosystem (EBUE) off Cape Blanc, Mauritania, is analysed with respect to deep ocean mass fluxes, flux components and their variability on seasonal to decadal timescales. The total mass flux revealed interannual fluctuations which were superimposed by fluctuations on decadal timescales. High winter fluxes of biogenic silica (BSi), used as a measure of marine production (mostly by diatoms) largely correspond to a positive North Atlantic Oscillation (NAO) index (December-March). However, this relationship is weak. The highest positive BSi anomaly was in winter 2004-2005 when the NAO was in a neutral state. More episodic BSi sedimentation events occurred in several summer seasons between 2001 and 2005, when the previous winter NAO was neutral or even negative. We suggest that distinct dust outbreaks and deposition in the surface ocean in winter and occasionally in summer/autumn enhanced particle sedimentation and carbon export on short timescales via the ballasting effect. Episodic perturbations of the marine carbon cycle by dust outbreaks (e.g. in 2005) might have weakened the relationships between fluxes and large-scale climatic oscillations. As phytoplankton biomass is high throughout the year, any dry (in winter) or wet (in summer) deposition of fine-grained dust particles is assumed to enhance the efficiency of the biological pump by incorporating dust into dense and fast settling organic-rich aggregates. A good correspondence between BSi and dust fluxes was observed for the dusty year 2005, following a period of rather dry conditions in the Sahara/Sahel region. Large changes of all bulk fluxes occurred during the strongest El Niño-Southern Oscillation (ENSO) in 1997-1999 where low fluxes were obtained for almost 1 year during the warm El Niño and high fluxes in the following cold La Niña phase. For decadal timescales, Bakun (1990) suggested an intensification of coastal upwelling due to increased winds (''Bakun upwelling intensification hypothesis''; Cropper et al., 2014) and global climate change. We did not observe an increase of any flux component off Cape Blanc during the past 2 and a half decades which might support this. Furthermore, fluxes of mineral dust did not show any positive or negative trends over time which might suggest enhanced desertification or ''Saharan greening'' during the last few decades.

Observed Teleconnections in Northern Winter: Subseasonal Evolution and Tropical Linkages

Teleconnections refer to the climate variability links between non-contiguous geographic regions, and tend to be associated with variability in both space and time of the climate’s semi-permanent circulation features. Teleconnections are well-developed in Northern winter, when they influence subseasonal-to-seasonal climate variability, notably, in surface temperature and precipitation. This work is comprised of four independent studies that improve understanding of tropical-extratropical teleconnections and their surface climate responses, subseasonal teleconnection evolution, and the utility of teleconnections in attribution of extreme climate events. After an introduction to teleconnection analysis as well as the major teleconnection patterns and associated climatic footprints manifest during Northern winter, the lagged impact of the Madden-Julian Oscillation (MJO) on subseasonal climate variability is presented. It is found that monitoring of MJO-related velocity potential anomalies is sufficient to predict MJO impacts. These impacts include, for example, the development of significant positive temperature anomalies over the eastern United States one to three weeks following an anomalous convective dipole with enhanced (suppressed) convection centered over the Indian Ocean (western Pacific). Subseasonal teleconnection evolution is assessed with respect to the Pacific-North America (PNA) pattern and the North Atlantic Oscillation (NAO). This evolution is analyzed both in the presence and absence of MJO-related circulation anomalies. It is found that removal of the MJO results only in small shifts in the centers of action of the NAO and PNA, and that in either case there is a small but significant lag in which the NAO leads a PNA pattern of opposite phase. Barotropic vorticity analysis suggests that this relationship may result in part from excitation of Rossby waves by the NAO in the Asian waveguide. An attempt is made to elegantly differentiate between the MJO extratropical response and patterns of variability more internal to the extratropics. Analysis of upper-level streamfunction anomalies is successful in this regard, and it is suggested that this is the preferred method for the real time monitoring of tropical-extratropical teleconnections. The extreme 2013-2014 North American winter is reconstructed using teleconnection analysis, and it is found that the North Pacific Oscillation-West Pacific (NPO/WP) pattern was the leading contributor to climate anomalies over much of North America. Such attribution is cautionary given the propensity to implicate the tropics for all midlatitude climate anomalies based on the El Niño-Southern Oscillation (ENSO) paradigm. A recent hypothesis of such tropical influence is presented and challenged.

The role of the Indian Ocean sector and sea surface salinity for prediction of the coupled Indo-Pacific system

The purpose of this dissertation is to evaluate the potential downstream influence of the Indian Ocean (IO) on El Niño/Southern Oscillation (ENSO) forecasts through the oceanic pathway of the Indonesian Throughflow (ITF), atmospheric teleconnections between the IO and Pacific, and assimilation of IO observations. Also the impact of sea surface salinity (SSS) in the Indo-Pacific region is assessed to try to address known problems with operational coupled model precipitation forecasts. The ITF normally drains warm fresh water from the Pacific reducing the mixed layer depths (MLD). A shallower MLD amplifies large-scale oceanic Kelvin/Rossby waves thus giving ~10% larger response and more realistic ENSO sea surface temperature (SST) variability compared to observed when the ITF is open. In order to isolate the impact of the IO sector atmospheric teleconnections to ENSO, experiments are contrasted that selectively couple/decouple the interannual forcing in the IO. The interannual variability of IO SST forcing is responsible for 3 month lagged widespread downwelling in the Pacific, assisted by off-equatorial curl, leading to warmer NINO3 SST anomaly and improved ENSO validation (significant from 3-9 months). Isolating the impact of observations in the IO sector using regional assimilation identifies large-scale warming in the IO that acts to intensify the easterlies of the Walker circulation and increases pervasive upwelling across the Pacific, cooling the eastern Pacific, and improving ENSO validation (r ~ 0.05, RMS~0.08C). Lastly, the positive impact of more accurate fresh water forcing is demonstrated to address inadequate precipitation forecasts in operational coupled models. Aquarius SSS assimilation improves the mixed layer density and enhances mixing, setting off upwelling that eventually cools the eastern Pacific after 6 months, counteracting the pervasive warming of most coupled models and significantly improving ENSO validation from 5-11 months. In summary, the ITF oceanic pathway, the atmospheric teleconnection, the impact of observations in the IO, and improved Indo-Pacific SSS are all responsible for ENSO forecast improvements, and so each aspect of this study contributes to a better overall understanding of ENSO. Therefore, the upstream influence of the IO should be thought of as integral to the functioning of ENSO phenomenon.

Probabilistic methods for seasonal forecasting in a changing climate: Cox-type regression models.

For climate risk management, cumulative distribution functions (CDFs) are an important source of information. They are ideally suited to compare probabilistic forecasts of primary (e.g. rainfall) or secondary data (e.g. crop yields). Summarised as CDFs, such forecasts allow an easy quantitative assessment of possible, alternative actions. Although the degree of uncertainty associated with CDF estimation could influence decisions, such information is rarely provided. Hence, we propose Cox-type regression models (CRMs) as a statistical framework for making inferences on CDFs in climate science. CRMs were designed for modelling probability distributions rather than just mean or median values. This makes the approach appealing for risk assessments where probabilities of extremes are often more informative than central tendency measures. CRMs are semi-parametric approaches originally designed for modelling risks arising from time-to-event data. Here we extend this original concept beyond time-dependent measures to other variables of interest. We also provide tools for estimating CDFs and surrounding uncertainty envelopes from empirical data. These statistical techniques intrinsically account for non-stationarities in time series that might be the result of climate change. This feature makes CRMs attractive candidates to investigate the feasibility of developing rigorous global circulation model (GCM)-CRM interfaces for provision of user-relevant forecasts. To demonstrate the applicability of CRMs, we present two examples for El Ni ? no/Southern Oscillation (ENSO)-based forecasts: the onset date of the wet season (Cairns, Australia) and total wet season rainfall (Quixeramobim, Brazil). This study emphasises the methodological aspects of CRMs rather than discussing merits or limitations of the ENSO-based predictors.