Evaluation of flood risk in response to climate variability

Standard procedures for forecasting flood risk (Bulletin 17B) assume annual maximum flood (AMF) series are stationary, meaning the distribution of flood flows is not significantly affected by climatic trends/cycles, or anthropogenic activities within the watershed. Historical flood events are therefore considered representative of future flood occurrences, and the risk associated with a given flood magnitude is modeled as constant over time. However, in light of increasing evidence to the contrary, this assumption should be reconsidered, especially as the existence of nonstationarity in AMF series can have significant impacts on planning and management of water resources and relevant infrastructure. Research presented in this thesis quantifies the degree of nonstationarity evident in AMF series for unimpaired watersheds throughout the contiguous U.S., identifies meteorological, climatic, and anthropogenic causes of this nonstationarity, and proposes an extension of the Bulletin 17B methodology which yields forecasts of flood risk that reflect climatic influences on flood magnitude. To appropriately forecast flood risk, it is necessary to consider the driving causes of nonstationarity in AMF series. Herein, large-scale climate patterns—including El Niño-Southern Oscillation (ENSO), Pacific Decadal Oscillation (PDO), North Atlantic Oscillation (NAO), and Atlantic Multidecadal Oscillation (AMO)—are identified as influencing factors on flood magnitude at numerous stations across the U.S. Strong relationships between flood magnitude and associated precipitation series were also observed for the majority of sites analyzed in the Upper Midwest and Northeastern regions of the U.S. Although relationships between flood magnitude and associated temperature series are not apparent, results do indicate that temperature is highly correlated with the timing of flood peaks. Despite consideration of watersheds classified as unimpaired, analyses also suggest that identified change-points in AMF series are due to dam construction, and other types of regulation and diversion. Although not explored herein, trends in AMF series are also likely to be partially explained by changes in land use and land cover over time. Results obtained herein suggest that improved forecasts of flood risk may be obtained using a simple modification of the Bulletin 17B framework, wherein the mean and standard deviation of the log-transformed flows are modeled as functions of climate indices associated with oceanic-atmospheric patterns (e.g. AMO, ENSO, NAO, and PDO) with lead times between 3 and 9 months. Herein, one-year ahead forecasts of the mean and standard deviation, and subsequently flood risk, are obtained by applying site specific multivariate regression models, which reflect the phase and intensity of a given climate pattern, as well as possible impacts of coupling of the climate cycles. These forecasts of flood risk are compared with forecasts derived using the existing Bulletin 17B model; large differences in the one-year ahead forecasts are observed in some locations. The increased knowledge of the inherent structure of AMF series and an improved understanding of physical and/or climatic causes of nonstationarity gained from this research should serve as insight for the formulation of a physical-casual based statistical model, incorporating both climatic variations and human impacts, for flood risk over longer planning horizons (e.g., 10-, 50, 100-years) necessary for water resources design, planning, and management.

Use of hydroclimatic forecasts for improved water management in central Texas

Accurate seasonal to interannual streamflow forecasts based on climate information are critical for optimal management and operation of water resources systems. Considering most water supply systems are multipurpose, operating these systems to meet increasing demand under the growing stresses of climate variability and climate change, population and economic growth, and environmental concerns could be very challenging. This study was to investigate improvement in water resources systems management through the use of seasonal climate forecasts. Hydrological persistence (streamflow and precipitation) and large-scale recurrent oceanic-atmospheric patterns such as the El Niño/Southern Oscillation (ENSO), Pacific Decadal Oscillation (PDO), North Atlantic Oscillation (NAO), the Atlantic Multidecadal Oscillation (AMO), the Pacific North American (PNA), and customized sea surface temperature (SST) indices were investigated for their potential to improve streamflow forecast accuracy and increase forecast lead-time in a river basin in central Texas. First, an ordinal polytomous logistic regression approach is proposed as a means of incorporating multiple predictor variables into a probabilistic forecast model. Forecast performance is assessed through a cross-validation procedure, using distributions-oriented metrics, and implications for decision making are discussed. Results indicate that, of the predictors evaluated, only hydrologic persistence and Pacific Ocean sea surface temperature patterns associated with ENSO and PDO provide forecasts which are statistically better than climatology. Secondly, a class of data mining techniques, known as tree-structured models, is investigated to address the nonlinear dynamics of climate teleconnections and screen promising probabilistic streamflow forecast models for river-reservoir systems. Results show that the tree-structured models can effectively capture the nonlinear features hidden in the data. Skill scores of probabilistic forecasts generated by both classification trees and logistic regression trees indicate that seasonal inflows throughout the system can be predicted with sufficient accuracy to improve water management, especially in the winter and spring seasons in central Texas. Lastly, a simplified two-stage stochastic economic-optimization model was proposed to investigate improvement in water use efficiency and the potential value of using seasonal forecasts, under the assumption of optimal decision making under uncertainty. Model results demonstrate that incorporating the probabilistic inflow forecasts into the optimization model can provide a significant improvement in seasonal water contract benefits over climatology, with lower average deficits (increased reliability) for a given average contract amount, or improved mean contract benefits for a given level of reliability compared to climatology. The results also illustrate the trade-off between the expected contract amount and reliability, i.e., larger contracts can be signed at greater risk.

Analysis of atmospheric transport and the effects of seasonal and interannual transport variability on measurements made at the Pico Mountain observatory

A non-hierarchical K-means algorithm is used to cluster 47 years (1960–2006) of 10-day HYSPLIT backward trajectories to the Pico Mountain (PM) observatory on a seasonal basis. The resulting cluster centers identify the major transport pathways and collectively comprise a long-term climatology of transport to the observatory. The transport climatology improves our ability to interpret the observations made there and our understanding of pollution source regions to the station and the central North Atlantic region. I determine which pathways dominate transport to the observatory and examine the impacts of these transport patterns on the O3, NOy, NOx, and CO measurements made there during 2001–2006. Transport from the U.S., Canada, and the Atlantic most frequently reaches the station, but Europe, east Africa, and the Pacific can also contribute significantly depending on the season. Transport from Canada was correlated with the North Atlantic Oscillation (NAO) in spring and winter, and transport from the Pacific was uncorrelated with the NAO. The highest CO and O3 are observed during spring. Summer is also characterized by high CO and O3 and the highest NOy and NOx of any season. Previous studies at the station attributed the summer time high CO and O3 to transport of boreal wildfire emissions (for 2002–2004), and boreal fires continued to affect the station during 2005 and 2006. The particle dispersion model FLEXPART was used to calculate anthropogenic and biomass-burning CO tracer values at the station in an attempt to identify the regions responsible for the high CO and O3 observations during spring and biomass-burning impacts in summer.

North-south palaeohydrological contrasts in the central Mediterranean during the Holocene: tentative synthesis and working hypotheses

On the basis of a multi-proxy approach and a strategy combining lacustrine and marine records along a north–south transect, data collected in the central Mediterranean within the framework of a collaborative project have led to reconstruction of high-resolution and well-dated palaeohydrological records and to assessment of their spatial and temporal coherency. Contrasting patterns of palaeohydrological changes have been evidenced in the central Mediterranean: south (north) of around 40° N of latitude, the middle part of the Holocene was characterised by lake-level maxima (minima), during an interval dated to ca. 10 300–4500 cal BP to the south and 9000–4500 cal BP to the north. Available data suggest that these contrasting palaeohydrological patterns operated throughout the Holocene, both on millennial and centennial scales. Regarding precipitation seasonality, maximum humidity in the central Mediterranean during the middle part of the Holocene was characterised by humid winters and dry summers north of ca. 40° N, and humid winters and summers south of ca. 40° N. This may explain an apparent conflict between palaeoclimatic records depending on the proxies used for reconstruction as well as the synchronous expansion of tree species taxa with contrasting climatic requirements. In addition, south of ca. 40° N, the first millennium of the Holocene was characterised by very dry climatic conditions not only in the eastern, but also in the central- and the western Mediterranean zones as reflected by low lake levels and delayed reforestation. These results suggest that, in addition to the influence of the Nile discharge reinforced by the African monsoon, the deposition of Sapropel 1 has been favoured (1) by an increase in winter precipitation in the northern Mediterranean borderlands, and (2) by an increase in winter and summer precipitation in the southern Mediterranean area. The climate reversal following the Holocene climate optimum appears to have been punctuated by two major climate changes around 7500 and 4500 cal BP. In the central Mediterranean, the Holocene palaeohydrological changes developed in response to a combination of orbital, ice-sheet and solar forcing factors. The maximum humidity interval in the south-central Mediterranean started ca. 10 300 cal BP, in correlation with the decline (1) of the possible blocking effects of the North Atlantic anticyclone linked to maximum insolation, and/or (2) of the influence of the remnant ice sheets and fresh water forcing in the North Atlantic Ocean. In the north-central Mediterranean, the lake-level minimum interval began only around 9000 cal BP when the Fennoscandian ice sheet disappeared and a prevailing positive NAO-(North Atlantic Oscillation) type circulation developed in the North Atlantic area. The major palaeohydrological oscillation around 4500–4000 cal BP may be a non-linear response to the gradual decrease in insolation, with additional key seasonal and interhemispheric changes. On a centennial scale, the successive climatic events which punctuated the entire Holocene in the central Mediterranean coincided with cooling events associated with deglacial outbursts in the North Atlantic area and decreases in solar activity during the interval 11 700–7000 cal BP, and to a possible combination of NAO-type circulation and solar forcing since ca. 7000 cal BP onwards. Thus, regarding the centennial-scale climatic oscillations, the Mediterranean Basin appears to have been strongly linked to the North Atlantic area and affected by solar activity over the entire Holocene. In addition to model experiments, a better understanding of forcing factors and past atmospheric circulation patterns behind the Holocene palaeohydrological changes in the Mediterranean area will require further investigation to establish additional high-resolution and well-dated records in selected locations around the Mediterranean Basin and in adjacent regions. Special attention should be paid to greater precision in the reconstruction, on millennial and centennial timescales, of changes in the latitudinal location of the limit between the northern and southern palaeohydrological Mediterranean sectors, depending on (1) the intensity and/or characteristics of climatic periods/oscillations (e.g. Holocene thermal maximum versus Neoglacial, as well as, for instance, the 8.2 ka event versus the 4 ka event or the Little Ice Age); and (2) on varying geographical conditions from the western to the eastern Mediterranean areas (longitudinal gradients). Finally, on the basis of projects using strategically located study sites, there is a need to explore possible influences of other general atmospheric circulation patterns than NAO, such as the East Atlantic–West Russian or North Sea–Caspian patterns, in explaining the apparent complexity of palaeoclimatic (palaeohydrological) Holocene records from the Mediterranean area.

Changing correlation structures of the Northern Hemisphere atmospheric circulation from 1000 to 2100 AD

Atmospheric circulation modes are important concepts in understanding the variability of atmospheric dynamics. Assuming their spatial patterns to be fixed, such modes are often described by simple indices from rather short observational data sets. The increasing length of reanalysis products allows these concepts and assumptions to be scrutinised. Here we investigate the stability of spatial patterns of Northern Hemisphere teleconnections by using the Twentieth Century Reanalysis as well as several control and transient millennium-scale simulations with coupled models. The observed and simulated centre of action of the two major teleconnection patterns, the North Atlantic Oscillation (NAO) and to some extent the Pacific North American (PNA), are not stable in time. The currently observed dipole pattern of the NAO, its centre of action over Iceland and the Azores, split into a north–south dipole pattern in the western Atlantic with a wave train pattern in the eastern part, connecting the British Isles with West Greenland and the eastern Mediterranean during the period 1940–1969 AD. The PNA centres of action over Canada are shifted southwards and over Florida into the Gulf of Mexico during the period 1915–1944 AD. The analysis further shows that shifts in the centres of action of either teleconnection pattern are not related to changes in the external forcing applied in transient simulations of the last millennium. Such shifts in their centres of action are accompanied by changes in the relation of local precipitation and temperature with the overlying atmospheric mode. These findings further undermine the assumption of stationarity between local climate/proxy variability and large-scale dynamics inherent when using proxy-based reconstructions of atmospheric modes, and call for a more robust understanding of atmospheric variability on decadal timescales.

Decadal Variability in the Northeast Pacific in a Physical-Ecosystem Model: Role of Mixed Layer Depth and Trophic Interactions

A basin-wide interdecadal change in both the physical state and the ecology of the North Pacific occurred near the end of 1976. Here we use a physical-ecosystem model to examine whether changes in the physical environment associated with the 1976-1977 transition influenced the lower trophic levels of the food web and if so by what means. The physical component is an ocean general circulation model, while the biological component contains 10 compartments: two phytoplankton, two zooplankton, two detritus pools, nitrate, ammonium, silicate, and carbon dioxide. The model is forced with observed atmospheric fields during 1960-1999. During spring, there is a similar to 40% reduction in plankton biomass in all four plankton groups during 1977-1988 relative to 1970-1976 in the central Gulf of Alaska (GOA). The epoch difference in plankton appears to be controlled by the mixed layer depth. Enhanced Ekman pumping after 1976 caused the halocline to shoal, and thus the mixed layer depth, which extends to the top of the halocline in late winter, did not penetrate as deep in the central GOA. As a result, more phytoplankton remained in the euphotic zone, and phytoplankton biomass began to increase earlier in the year after the 1976 transition. Zooplankton biomass also increased, but then grazing pressure led to a strong decrease in phytoplankton by April followed by a drop in zooplankton by May: Essentially, the mean seasonal cycle of plankton biomass was shifted earlier in the year. As the seasonal cycle progressed, the difference in plankton concentrations between epochs reversed sign again, leading to slightly greater zooplankton biomass during summer in the later epoch.

A 13,600-year diatom oxygen isotope record from the South Carpathians (Romania): Reflection of winter conditions and possible links with North Atlantic circulation changes

This study provides a continuous lateglacial and Holocene record of diatom silica oxygen isotope changes (delta O-18(DIAT)) in a subalpine lake sediment sequence obtained from the Retezat Mts (Taul dintre Brazi, 1740 m a.s.l.). This through-flow, shallow, high-altitude lake with a surface area of only 0.4 ha has short water residence time and is predominantly fed by snowmelt and rainwater. Its delta O-18(DIAT) record principally reflects the oxygen isotope composition of the winter and spring precipitation, as diatom blooms occur mainly in the spring and early summer. Hence, changes in delta O-18(DIAT) are interpreted as seasonal scale changes: in the amount of winter precipitation. Low oxygen isotope values (27-28.5 parts per thousand) occurred during the lateglacial until 12,300 cal BP, followed by a sharp increase thereafter. In the Holocene delta O-18(DIAT) values ranged from 29 to 31 parts per thousand until 3200 cal BP, followed by generally lower values during the late Holocene (27-30 parts per thousand). Short-term decreases in the isotopic values were found between 10,140-9570, 9000-8500, 7800-7300, 6300-5800, 5500-5000 and at 8015, 4400, 4000 cal BP. After 3200 cal BP a decreasing trend was visible with the lowest values between 3100-2500 and after 2100 cal BP The general trend in the record suggests that contribution of winter precipitation was generally lower between 11,680 and 3200 cal BP, followed by increased contribution during the last millennia. The late Holocene decrease in delta O-18(DIAT) shows good agreement with the speleothem delta O-18, lake level and testate amoebae records from the Carpathian Mountains that also display gradual delta O-18 decrease and lake level/mire water table level rise after 3200 cal BR Strong positive correlation with North Atlantic circulation and solar activity proxies, such as the Austrian and Hungarian speleothem records, furthermore suggested that short-term increases in the isotopic ratios in the early and mid Holocene are likely connectable to high solar activity phases and high frequency of positive North Atlantic Oscillation indexes that may have resulted in decreased winter precipitation in this region.

Freshening of the Alaska Coastal Current Recorded By Coralline Algal Ba/Ca Ratios

Arctic Ocean freshening can exert a controlling influence on global climate, triggering strong feedbacks on ocean-atmospheric processes and affecting the global cycling of the world's oceans. Glacier-fed ocean currents such as the Alaska Coastal Current are important sources of freshwater for the Bering Sea shelf, and may also influence the Arctic Ocean freshwater budget. Instrumental data indicate a multiyear freshening episode of the Alaska Coastal Current in the early 21st century. It is uncertain whether this freshening is part of natural multidecadal climate variability or a unique feature of anthropogenically induced warming. In order to answer this, a better understanding of past variations in the Alaska Coastal Current is needed. However, continuous long-term high-resolution observations of the Alaska Coastal Current have only been available for the last 2 decades. In this study, specimens of the long-lived crustose coralline alga Clathromorphum nereostratum were collected within the pathway of the Alaska Coastal Current and utilized as archives of past temperature and salinity. Results indicate that coralline algal Mg/Ca ratios provide a 60 year record of sea surface temperatures and track changes of the Pacific Decadal Oscillation, a pattern of decadal-to-multidecadal ocean-atmosphere climate variability centered over the North Pacific. Algal Ba/Ca ratios (used as indicators of coastal freshwater runoff) are inversely correlated to instrumentally measured Alaska Coastal Current salinity and record the period of freshening from 2001 to 2006. Similar multiyear freshening events are not evident in the earlier portion of the 60 year Ba/Ca record. This suggests that the 21st century freshening of the Alaska Coastal Current is a unique feature related to increasing glacial melt and precipitation on mainland Alaska.

Benthic foraminiferal stable isotope record from ODP Site 138-849 on the East Pacific Rise

Benthic foraminifer and delta13C data from Site 849, on the west flank of the East Pacific Rise (0°11 'N, 110°31'W; 3851 m), give relatively continuous records of deep Pacific Ocean stable isotope variations between 0 and 5 Ma. The mean sample spacing is 4 k.y. Most analyses are from Cibicides wuellerstorfi, but isotopic offsets relative to Uvigerina peregrina appear roughly constant. Because of its location west of the East Pacific Rise, Site 849 yields a suitable record of mean Pacific Ocean delta13C, which approximates a global oceanic signal. The ~100-k.y.-period climate cycle, which is prevalent in delta18O does not dominate the long-term delta13C record. For delta13C, variations in the ~400- and 41-k.y. periods are more important. Phase lags of delta13C relative to ice volume in the 41- and 23-k.y. bands are consistent with delta13C as a measure of organic biomass. A model-calculated exponential response time of 1-2 k.y. is appropriate for carbon stored in soils and shallow sediments responding to glacial-interglacial climate change. Oceanic delta13C leads ice volume slightly in the 100-k.y. band, and this suggests another process such as changes in continental weathering to modulate mean river delta13C at long periods. The delta13C record from Site 849 diverges from that of Site 677 in the Panama Basin mostly because of decay of 13C-depleted organic carbon in the relatively isolated Panama Basin. North Atlantic to Pacific delta13C differences calculated using published data from Sites 607 and 849 reveal variations in Pliocene deep water within the range of those of the late Quaternary. Maximum delta13C contrast between these sites, which presumably reflects maximum influx of high-delta13C northern source water into the deep North Atlantic Ocean, occurred between 1.3 and 2.1 Ma, well after the initiation of Northern Hemisphere glaciation. Export of high-delta13C North Atlantic Deep Water from the Atlantic to the circumpolar Antarctic, as recorded by published delta13C data from Subantarctic Site 704, appears unrelated to the North Atlantic-Pacific delta13C contrast. To account for this observation, we suggest that deep-water formation in the North Atlantic reflects northern source characteristics, whereas export of this water into the circumpolar Antarctic reflects Southern Hemisphere wind forcing. Neither process appears directly linked to ice-volume variations.

Neogene calcareous nannofossils in the Caribbean and the eastern equatorial Pacific

Principal components analysis of nannofossil assemblages in five ODP sites in the Caribbean and the eastern equatorial Pacific was used to monitor Neogene surface water circulation changes. In 15.83-10.71 Ma, floral distributions clearly show the existence of the Circum-tropical Current between the Caribbean and the eastern equatorial Pacific and no surface water communication between the northern and southern Caribbean. In 10.71-9.36 Ma, the Circum-tropical Current weakened and the northward intra-Caribbean current had been initiated. Northern and southern Caribbean surface waters again became separated in 8.35-3.65 Ma and the Circum-tropical Current was regenerated. After 2.76 Ma, the northward intra-Caribbean current has completely been established and the Circum-tropical Current disappeared. This suggests that a barrier to the Circum-tropical Current had formed, indicating the final closure of the Isthmus of Panama.

Isotope measurements from late Paleocene - early Eocene of 5 cores from Atlantic and Pacific Ocean

High-resolution stable carbon isotope records for upper Paleocene - lower Eocene sections at Ocean Drilling Program Sites 1051 and 690 and Deep Sea Drilling Project Sites 550 and 577 show numerous rapid (40 - 60 kyr duration) negative excursions of up to 1 per mill. We demonstrate that these transient decreases are the expected result of nonlinear insolation forcing of the carbon cycle in the context of a long carbon residence time. The transients occur at maxima in Earth's orbital eccentricity, which result in high-amplitude variations in insolation due to forcing by climatic precession. The construction of accurate orbital chronologies for geologic sections older than ~ 35 Ma relies on identifying a high-fidelity recorder of variations in Earth's orbital eccentricity. We use the carbon isotope records as such a recorder, establishing a robust orbitally tuned chronology for latest Paleocene-earliest Eocene events. Moreover, the transient decreases provide a means of precise correlation among the four sites that is independent of magnetostratigraphic and biostratigraphic data at the <10^5-year scale. While the eccentricity-controlled transient decreases bear some resemblance to the much larger-amplitude carbon isotope excursion (CIE) that marks the Paleocene/Eocene boundary, the latter event is found to occur near a minimum in the ~400-kyr eccentricity cycle. Thus the CIE occurred during a time of minimal variability in insolation, the dominant mechanism for forcing climate change on 104-year scales. We argue that this is inconsistent with mechanisms that rely on a threshold climate event to trigger the Paleocene/Eocene thermal maximum since any threshold would more likely be crossed during a period of high-amplitude climate variations.

(Table 2) Compressional and shear wave velocites from basalts and sediments of the Atlantic and Pacific Ocean

Compressional (Vp) and shear (Vs) wave velocities have been measured to 1.0 kbar for 14 cores of well-consolidated sedimentary rock from Atlantic and Pacific sites of the Deep Sea Drilling Project. The range of VP (2.05-5.38 km/sec at 0.5 kbar) shows significant overlap with the range of oceanic layer-2 seismic velocities determined by marine refraction surveys, suggesting that sedimentary rocks may, in some regions, constitute the upper portion of layer 2. Differing linear relationships between VP and Vs for basalts and sedimentary rocks, however, may provide a method of resolving layer-2 composition. This is illustra ted for a refraction survey site on the flank of the Mid-Atlantic Ridge where layer-2 velocities agree with basalt, and two sites on the Saya de Malha Bank in the Indian Ocean where layer-2 velocities appear to represent sedimentary rock.

Stable carbon isotope composition of benthic foraminifera from sediments of the Skagerrak, North Sea

The sediment cores 225514 and 225510 were recovered from 420 and 285 m water depth, respectively. They were investigated for their benthic foraminiferal delta13C during the last 500 years. Both cores were recovered from the southern flank of the Skagerrak. The delta13C values of Uvigerina mediterranea and other shallow infaunal species in both cores indicate that organic matter rain rates to the seafloor varied around a mean value until approximately AD 1950 after which they increased. This increase might result from changes in the North Atlantic Current System and a co-occurring persistently high North Atlantic Oscillation index state in the 1980s to 1990s, rather than from anthropogenic eutrophication. Using delta13C mean values of multiple species, we reconstruct delta13C gradients of dissolved inorganic carbon (DIC) within pore waters for the time periods AD 1500 to 1950 and AD 1950 to 2000. The calculated delta13CDIC ranges, interpreted as indicating total organic matter remineralization due to respiration, are generally bigger in Core 225514 than in Core 225510. Since mean delta13C values of U. mediterranea suggest that organic matter rain rates were similar at both locations, differences in total organic matter remineralization are attributed to differing oxygen availability. However, oxygen concentrations in the overlying bottom water masses are not likely to have differed significantly. Thus, we suggest that organic matter remineralization was controlled by oxygen availability within the sediments, reflecting strong differences in sedimentation rates at the two investigated core sites. Based on the assumptions that tests of benthic foraminiferal species inhabiting the same microhabitat depth should show equal delta13C values unless they are affected by vital effects and that Globobulimina turgida records pore water delta13CDIC, we estimate microhabitat-corrected vital effects for several species with respect to G. turgida: >0.7 per mil for Cassidulina laevigata, >1.3 per mil for Hyalinea balthica, and >0.7 per mil for Melonis barleeanus. Melonis zaandami seems to closely record pore water delta13CDIC.