ENSO representation in climate models: from CMIP3 to CMIP5

We analyse the ability of CMIP3 and CMIP5 coupled ocean–atmosphere general circulation models (CGCMs) to simulate the tropical Pacific mean state and El Niño-Southern Oscillation (ENSO). The CMIP5 multi-model ensemble displays an encouraging 30 % reduction of the pervasive cold bias in the western Pacific, but no quantum leap in ENSO performance compared to CMIP3. CMIP3 and CMIP5 can thus be considered as one large ensemble (CMIP3 + CMIP5) for multi-model ENSO analysis. The too large diversity in CMIP3 ENSO amplitude is however reduced by a factor of two in CMIP5 and the ENSO life cycle (location of surface temperature anomalies, seasonal phase locking) is modestly improved. Other fundamental ENSO characteristics such as central Pacific precipitation anomalies however remain poorly represented. The sea surface temperature (SST)-latent heat flux feedback is slightly improved in the CMIP5 ensemble but the wind-SST feedback is still underestimated by 20–50 % and the shortwave-SST feedbacks remain underestimated by a factor of two. The improvement in ENSO amplitudes might therefore result from error compensations. The ability of CMIP models to simulate the SST-shortwave feedback, a major source of erroneous ENSO in CGCMs, is further detailed. In observations, this feedback is strongly nonlinear because the real atmosphere switches from subsident (positive feedback) to convective (negative feedback) regimes under the effect of seasonal and interannual variations. Only one-third of CMIP3 + CMIP5 models reproduce this regime shift, with the other models remaining locked in one of the two regimes. The modelled shortwave feedback nonlinearity increases with ENSO amplitude and the amplitude of this feedback in the spring strongly relates with the models ability to simulate ENSO phase locking. In a final stage, a subset of metrics is proposed in order to synthesize the ability of each CMIP3 and CMIP5 models to simulate ENSO main characteristics and key atmospheric feedbacks.

Mid-holocene and last glacial maximum climate simulations with the IPSL model-part I: comparing IPSLCM5A to IPSLCM4

The climates of the mid-Holocene (MH), 6,000 years ago, and of the Last Glacial Maximum (LGM), 21,000 years ago, have extensively been simulated, in particular in the framework of the Palaeoclimate Modelling Intercomparion Project. These periods are well documented by paleo-records, which can be used for evaluating model results for climates different from the present one. Here, we present new simulations of the MH and the LGM climates obtained with the IPSL_CM5A model and compare them to our previous results obtained with the IPSL_CM4 model. Compared to IPSL_CM4, IPSL_CM5A includes two new features: the interactive representation of the plant phenology and marine biogeochemistry. But one of the most important differences between these models is the latitudinal resolution and vertical domain of their atmospheric component, which have been improved in IPSL_CM5A and results in a better representation of the mid-latitude jet-streams. The Asian monsoon’s representation is also substantially improved. The global average mean annual temperature simulated for the pre-industrial (PI) period is colder in IPSL_CM5A than in IPSL_CM4 but their climate sensitivity to a CO2 doubling is similar. Here we show that these differences in the simulated PI climate have an impact on the simulated MH and LGM climatic anomalies. The larger cooling response to LGM boundary conditions in IPSL_CM5A appears to be mainly due to differences between the PMIP3 and PMIP2 boundary conditions, as shown by a short wave radiative forcing/feedback analysis based on a simplified perturbation method. It is found that the sensitivity computed from the LGM climate is lower than that computed from 2 × CO2 simulations, confirming previous studies based on different models. For the MH, the Asian monsoon, stronger in the IPSL_CM5A PI simulation, is also more sensitive to the insolation changes. The African monsoon is also further amplified in IPSL_CM5A due to the impact of the interactive phenology. Finally the changes in variability for both models and for MH and LGM are presented taking the example of the El-Niño Southern Oscillation (ENSO), which is very different in the PI simulations. ENSO variability is damped in both model versions at the MH, whereas inconsistent responses are found between the two versions for the LGM. Part 2 of this paper examines whether these differences between IPSL_CM4 and IPSL_CM5A can be distinguished when comparing those results to palaeo-climatic reconstructions and investigates new approaches for model-data comparisons made possible by the inclusion of new components in IPSL_CM5A.

The Climate-system Historical Forecast Project: do stratosphere-resolving models make better seasonal climate predictions in boreal winter?

Using an international, multi-model suite of historical forecasts from the World Climate Research Programme (WCRP) Climate-system Historical Forecast Project (CHFP), we compare the seasonal prediction skill in boreal wintertime between models that resolve the stratosphere and its dynamics (“high-top”) and models that do not (“low-top”). We evaluate hindcasts that are initialized in November, and examine the model biases in the stratosphere and how they relate to boreal wintertime (Dec-Mar) seasonal forecast skill. We are unable to detect more skill in the high-top ensemble-mean than the low-top ensemble-mean in forecasting the wintertime North Atlantic Oscillation, but model performance varies widely. Increasing the ensemble size clearly increases the skill for a given model. We then examine two major processes involving stratosphere-troposphere interactions (the El Niño-Southern Oscillation/ENSO and the Quasi-biennial Oscillation/QBO) and how they relate to predictive skill on intra-seasonal to seasonal timescales, particularly over the North Atlantic and Eurasia regions. High-top models tend to have a more realistic stratospheric response to El Niño and the QBO compared to low-top models. Enhanced conditional wintertime skill over high-latitudes and the North Atlantic region during winters with El Niño conditions suggests a possible role for a stratospheric pathway.

Caracterização de processos morfodinâmicos e hidrodinâmicos do cinturão lacustre meridional da Reserva Biológica do Lago Pirituba, Amapá

The Amapá State has an important natural lake system, known as The Amapá Lakes Region . Most of these lakes are on the southern part of Amapá s coastal plain, which has 300 km of extension and it s composed by holocenic sediments deposited at the northern part of Amazon River to the Orange Cape located on the northern part of Amapá state. This region is under influence of the Amazon River discharge which is the largest liquid discharge of about 209.000 m³/s and biggest sediment budget discharged on the ocean in the order 6.108 ton per day. The climate is influenced by the Intertropical Convergence Zone and El Niño Southern Oscillation which act mainly under precipitation, nebulosity, local rivers and tidal hidrology. In this region lake belts are Ocidental, Oriental and Meridional Lake Belts. The last one is formed by the by the lakes Comprido de Cima, Botos, Bacia, Lodão, Ventos, Mutuco and Comprido de Baixo. These lakes are the closest to the Araguari River and are characterized by pelitic sedimentation associated with fluvial and estuarine flood plains under influence of tides. The lakes are interconnected, suffer influence of flood pulses from the Tartarugal, Tartarugalzinho and Araguari rivers and the hydrodynamic and morphodynamic know edge is poor. Volume and area reduction, natural eutrophication, anthophic influence, hidrodynamic alterations, morphological changes and are factors which can contribute to the closing of such lakes on the Meridional Lake Belt. This belt is inside the boundaries of the Biological Reserve of Piratuba Lake, created in 1980 for integral protection. Due to the fragility of the environment together with the poor knowledge of the system and with the study area relevancy it is necessary to know the hydrodynamic and geoenvironmental processes. This work aims the characterization of morphodynamic and hydrodynamic processes in order to understand the geoambiental context of the Meridional Lake Belt, from the Comprido de Baixo Lake to the dos Ventos Lake, including the Tabaco Igarape. Methodology was based on the hydrodynamic data acquisition: liquid discharge (acoustic method), tides, bathymetry and the interpretation of multitemporal remote sensing images, integrated in a Geographic Information System (GIS). By this method charts of the medium liquid discharges of Lake Mutuco and Tabacco Igarape the maximum velocity of flow were estimated in: 1.1 m/s, 1.6 m/s and 1.6 m/s (rainy season) and 0.6 m/s, 0.6 m/s and 0.7 m/s (dry period), the maximum flow in: 289 m³/s, 297 m³/s and 379 m³/s (rainy season) and 41 m³/s , 79 m³/s and 105 m³/s (dry period), respectively. From the interpretation of multitemporal satellite images, maps were developed together with the analysis of the lakes and Tobaco Igarape evolution from 1972 to 2008, and were classified according to the degree of balance in the area: stable areas, eutrophic areas, areas of gain, and eroded areas. Troughout analysis of the balance of areas, it was possible to quantify the volume of lake areas occupied by aquatic macrophytes. The study sought to understand the hydrodynamic and morphodynamic processes occurring in the region, contributing to the elucidation of the processes which cause and/or favor geoenvironmental changes in the region; all such information is fundamental to making the management of the area and further definition of parameters for environmental monitoring and contributing to the development of the management plan of the Biological Reserve of Lake Piratuba. The work activities is a part of the Project "Integration of Geological, geophysical and geochemical data to Paleogeographic rebuilding of Amazon Coast, from the Neogene to the Recent

The economics of climate change in Central America: summary 2010

Includes Bibliography

Variabilidade da precipitação em tempo e espaço associada à Zona de Convergência Intertropical

Este estudo visa apresentar uma análise atmosférica da variabilidade espacial e temporal da Zona de Convergência Intertropical (ZCIT) nas cidades de Belém, Jakarta e Nairóbi, que estão localizadas sobre os continentes da América do Sul, Ásia e África, respectivamente. Para isso, foram utilizados dados diários de precipitação observada e radiação de onda longa para o período de 1999 a 2008, e aplicadas as técnicas matemáticas e estatísticas, como a média aritmética e a transformada em ondeletas Morlet. Em geral, os resultados indicam que do ponto de vista espacial, a precipitação mensal varia consideravelmente, pois as três cidades estudadas localizam-se em diferentes continentes da faixa tropical. Isto ocorre principalmente, durante os meses de Janeiro a Maio, período de maior atuação da ZCIT no hemisfério sul. As variações atmosféricas observadas, a partir dos escalogramas de fase, - de ondeleta indicam que as escalas interdecadal, anual, interanual e intrassazonal são moduladoras da precipitação. Tais escalas podem ser representadas pelos mecanismos oceano-atmosfera dos fenômenos El Niño Oscilação Sul e da oscilação intrassazonal de Madden e Julian. A contribuição destes fenômenos na distribuição da chuva nessas regiões é evidente durante o período estudado, sendo que Nairóbi, apesar de estar localizada em latitude semelhante à de Belém, apresenta pouca evidência do ciclo anual e forte na escala interdecadal. No caso de Belém e de Jakarta as oscilações de múltiescala de precipitação concentram-se nas escalas dos mecanismos moduladores da chuva associados com o ciclo anual e intrassazonal, durante todo o período.

A oscilação Madden - Julian na Amazônia Oriental: variáveis superfíciais

Tendo como foco as múltiplas escalas de tempo que atuam na Amazônia, este trabalho foi desenvolvido com o objetivo de investigar a possível influencia da Oscilação Madden – Julian (OMJ) em elementos turbulentos da CLP. A OMJ foi identificada a partir de 30 anos de dados de reanálise de radiação de onda longa (ROL) e componente zonal do vento (u). As grandezas turbulentas foram estudadas a partir da variância, covariância e coeficiente de correlação de um conjunto de dados de resposta rápida coletado na torre micrometeorológica de Caxiuanã (PA), e tratados com a Transformada em Ondeletas (TO) para se obter a contribuição de cada escala para estes momentos estatísticos. A análise dos 30 anos de dados de ROL e u mostrou que a ocorrência da OMJ está ligada com o fenômeno do El Niño/Oscilação Sul (ENOS), bem como influência do ENOS no tempo da região amazônica pode estar associado a presença ou não da OMJ. Foi observado que anos de El Niño tendem a desfavorecer a ocorrência da OMJ e anos de La Niña tendem a favorecer o desenvolvimento da oscilação. Caso uma OMJ se desenvolva durante um episodio de El Niño, a oscilação pode influenciar a temperatura, a velocidade do vento e a precipitação de forma diferente ao do El Niño. A análise por fase da OMJ mostrou que, em Belém, há diferença significativa na temperatura máxima e na precipitação entre cada fase, porém, a temperatura mínima e o módulo do vento apresentaram pouca diferença. Os fluxos cinemáticos turbulentos analisados, por escala, em três horários distintos, foram mais diferentes durante o período diurno, principalmente w’T’ e w’q’. A diferença entre fase ativa e fase inativa foi reduzindo com passar do dia, durante o período de transição dia – noite, poucas escalas tiveram diferença significativa, e durante a noite, nenhuma escala teve nível de confiança acima ou igual a 95%. Estes resultados indicam que a convecção diurna é o mecanismo responsável por esta diferença e como a OMJ atua como uma grande célula convectiva, a convecção local é amplificada, explicando a grande diferença observada entre as fases durante o período diurno.

Aspectos regionais da variabilidade de precipitação no estado do Pará: estudo observacional e modelagem climática em alta resolução

O objetivo do presente trabalho foi agregar diferentes redes de estações meteorológicas de superfície para a criação de um novo banco de dados integrado, a partir do qual foi gerada uma climatologia recente (1978-2007) para a precipitação do estado do Pará em alta resolução espacial – 30 km, permitindo melhor identificar a variabilidade climática regional, sobretudo influenciada pelos aspectos da fisiografia e em função de mecanismos climáticos de grande escala dos oceanos Pacífico e Atlântico. Buscou-se, ainda, estabelecer uma configuração otimizada do modelo climático RegCM3 utilizando duas diferentes parametrizações de cumulus: RegCM3/Grell e RegCM3/MIT. Foram realizadas 26 simulações (1982/83 a 2007/08) durante a estação chuvosa na Amazônia oriental (dezembro a maio) para cada esquema de parametrização convectiva, utilizando 30 km de resolução espacial. O modelo mostrou-se capaz de capturar os sinais de anomalia na presença de forçantes climáticas extremas, como o El Niño-Oscilação Sul e o dipolo do Atlântico. O RegCM3/MIT obteve ótimo desempenho na região de Altamira/PA e performance razoável nos setores Nordeste (região de Belém), Leste ( região de Marabá), Sudeste (região de Conceição do Araguaia), e Noroeste (região de Tiriós). O RegCM3/Grell destacou-se nas regiões Nordeste, Leste, Sudeste e Noroeste, com desempenho razoável. O setor Norte (região de Macapá) foi o mais problemático, com pouca ou nenhuma sensibilidade apresentada pelo modelo. Embora o RegCM3 tenha obtido resultados razoáveis na maior parte do domínio, foram detectados erros sistemáticos nas simulações, com viés seco para o RegCM3/Grell e viés úmido para o RegCM3/MIT na porção Sul e viés seco na porção Norte. Estas características denotam a necessidade de ajustes às condições regionais dos esquemas de convecção.

Study of the inter-annual variability of particle vertical fluxes in two moorings in the Ross Sea (Antarctica)

The most ocean - atmosphere exchanges take place in polar environments due to the low temperatures which favor the absorption processes of atmospheric gases, in particular CO2. For this reason, the alterations of biogeochemical cycles in these areas can have a strong impact on the global climate. With the aim of contributing to the definition of the mechanisms regulating the biogeochemical fluxes we have analyzed the particles collected in the Ross Sea in different years (ROSSMIZE, BIOSESO 1 and 2, ROAVERRS and ABIOCLEAR projects) in two sites (mooring A and B). So it has been developed a more efficient method to prepare sediment trap samples for the analyses. We have also processed satellite data of sea ice, chlorophyll a and diatoms concentration. At both sites, in each year considered, there was a high seasonal and inter-annual variability of biogeochemical fluxes closely correlated with sea ice cover and primary productivity. The comparison between the samples collected at mooring A and B in 2008 highlighted the main differences between these two sites. Particle fluxes at Mooring A, located in a polynia area, are higher than mooring B ones and they happen about a month before. In the mooring B area it has been possible to correlate the particles fluxes to the ice concentration anomalies and with the atmospheric changes in response to El Niño Southern Oscillations. In 1996 and 1999, years subjected to La Niña, the concentrations of sea ice in this area have been less than in 1998, year subjected to El Niño. Inverse correlation was found for 2005 and 2008. In the mooring A area significant differences in mass and biogenic fluxes during 2005 and 2008 has been recorded. This allowed to underline the high variability of lateral advection processes and to connect them to the physical forcing.

Response of nutrient cycles of an old-growth montane forest in Ecuador to experimental low-level nutrient amendments

Atmospheric nitrogen (N) and phosphorus (P) depositions are expected to increase in the tropicsrnas a consequence of increasing human activities in the next decades. Furthermore, a possiblernshortened El Niño Southern Oscillation cycle might come along with more frequent calcium (Ca)rndepositions on the eastern slope of the Ecuadorian Andes originating from Saharan dust. It isrncrucial to understand the response of the old-growth montane forest in Ecuador to increasedrnnutrient deposition to predict the further development of this megadiverse ecosystem.rnI studied experimental additions of N, P, N+P and Ca to the forest and an untreatedrncontrol, all in a fourfold replicated randomized block design. These experiments were conductedrnin the framework of a collaborative research effort, the NUtrient Manipulation EXperimentrn(NUMEX). I collected litter leachate, mineral soil solution (0.15 and 0.30 m depths), throughfallrnand fine litterfall samples and determined N, P and Ca concentrations and fluxes. This approachrnalso allowed me to assess whether N, P and/or Ca are limiting nutrients for forest growth.rnFurthermore, I evaluated the response of fine root biomass, leaf area index, leaf area and specificrnleaf area, tree diameter growth and basal area increment contributed from a cooperating group inrnthe Ca applied and control treatments.rnDuring the observation period of 16 months after the first fertilizer application, less thanrn10, 1 and 5% of the applied N, P and Ca, respectively, leached below the organic layer whichrncontained almost all roots but no significant leaching losses occurred to the deeper mineral soil.rnDeposited N, P and Ca from the atmosphere in dry and wet form were, on balance, retained in therncanopy in the control treatment. Retention of N, P and Ca in the canopy in their respectiverntreatments was reduced resulting in higher concentrations and fluxes of N, P and Ca inrnthroughfall and litterfall. Up to 2.5% of the applied N and 2% of the applied P and Ca werernrecycled to the soil with throughfall. Fluxes of N, P and Ca in throughfall+litterfall were higher inrnthe fertilized treatments than in the control; up to 20, 5 and 25% of the applied N, P and Ca,rnrespectively, were recycled to the soil with throughfall+litterfall.rnIn the Ca-applied plots, fine root biomass decreased significantly. Also the leaf area of thernfour most common tree species tended to decrease and the specific leaf area increasedrnsignificantly in Graffenrieda emarginata Triana, the most common tree species in the study area.rnThese changes are known plant responses to reduced nutrient stress. Reduced aluminium (Al)rntoxicity as an explanation of the Ca effect was unlikely, because of almost complete organocomplexationrnof Al and molar Ca:Al concentration ratios in solution above the toxicity threshold.rnThe results suggest that N, P and Ca co-limit the forest ecosystem functioning in thernnorthern Andean montane forests in line with recent assumptions in which different ecosystemrncompartments and even different phenological stages may show different nutrient limitationsrn(Kaspari et al. 2008). I conclude that (1) the expected elevated N and P deposition will bernretained in the ecosystem, at least in the short term and hence, quality of river water will not bernendangered and (2) increased Ca input will reduce nutrient stress of the forest.

High-resolution reconstruction of Holocene climate variability and environmental changes in the North Pacific using bivalve shells

Bivalve mollusk shells are useful tools for multi-species and multi-proxy paleoenvironmental reconstructions with a high temporal and spatial resolution. Past environmental conditions can be reconstructed from shell growth and stable oxygen and carbon isotope ratios, which present an archive for temperature, freshwater fluxes and primary productivity. The purpose of this thesis is the reconstruction of Holocene climate and environmental variations in the North Pacific with a high spatial and temporal resolution using marine bivalve shells. This thesis focuses on several different Holocene time periods and multiple regions in the North Pacific, including: Japan, Alaska (AK), British Columbia (BC) and Washington State, which are affected by the monsoon, Pacific Decadal Oscillation (PDO) and El Niño/Southern Oscillation (ENSO). Such high-resolution proxy data from the marine realm of mid- and high-latitudes are still rare. Therefore, this study contributes to the optimization and verification of climate models. However, before using bivalves for environmental reconstructions and seasonality studies, life history traits must be well studied to temporally align and interpret the geochemical record. These calibration studies are essential to ascertain the usefulness of selected bivalve species as paleoclimate proxy archives. This work focuses on two bivalve species, the short-lived Saxidomus gigantea and the long-lived Panopea abrupta. Sclerochronology and oxygen isotope ratios of different shell layers of P. abrupta were studied in order to test the reliability of this species as a climate archive. The annual increments are clearly discernable in umbonal shell portions and the increments widths should be measured in these shell portions. A reliable reconstruction of paleotemperatures may only be achieved by exclusively sampling the outer shell layer of multiple contemporaneous specimens. Life history traits (e.g., timing of growth line formation, duration of the growing season and growth rates) and stable isotope ratios of recent S. gigantea from AK and BC were analyzed in detail. Furthermore, a growth-temperature model based on S. gigantea shells from Alaska was established, which provides a better understanding of the hydrological changes related to the Alaska Coastal Current (ACC). This approach allows the independent measurement of water temperature and salinity from variations in the width of lunar daily growth increments of S. gigantea. Temperature explains 70% of the variability in shell growth. The model was calibrated and tested with modern shells and then applied to archaeological specimens. The time period between 988 and 1447 cal yrs BP was characterized by colder (~1-2°C) and much drier (2-5 PSU) summers, and a likely much slower flowing ACC than at present. In contrast, the summers during the time interval of 599-1014 cal yrs BP were colder (up to 3°C) and fresher (1-2 PSU) than today. The Aleutian Low may have been stronger and the ACC was probably flowing faster during this time.

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.

Evaluating the usefulness of soil moisture (forecasts) in streamflow prediction for central Texas

The past decade has brought significant advancements in seasonal climate forecasting. However, water resources decision support and management continues to be based almost entirely on historical observations and does not take advantage of climate forecasts. This study builds on previous work that conditioned streamflow ensemble forecasts on observable climate indicators, such as the El Niño-Southern Oscillation (ENSO) and the Pacific Decadal Oscillation (PDO) for use in a decision support model for the Highland Lakes multi-reservoir system in central Texas operated by the Lower Colorado River Authority (LCRA). In the current study, seasonal soil moisture is explored as a climate indicator and predictor of annual streamflow for the LCRA region. The main purpose of this study is to evaluate the correlation of fractional soil moisture with streamflow using the 1950-2000 Variable Infiltration Capacity (VIC) Retrospective Land Surface Data Set over the LCRA region. Correlations were determined by examining different annual and seasonal combinations of VIC modeled fractional soil moisture and observed streamflow. The applicability of the VIC Retrospective Land Surface Data Set as a data source for this study is tested along with establishing and analyzing patterns of climatology for the watershed study area using the selected data source (VIC model) and historical data. Correlation results showed potential for the use of soil moisture as a predictor of streamflow over the LCRA region. This was evident by the good correlations found between seasonal soil moisture and seasonal streamflow during coincident seasons as well as between seasonal and annual soil moisture with annual streamflow during coincident years. With the findings of good correlation between seasonal soil moisture from the VIC Retrospective Land Surface Data Set with observed annual streamflow presented in this study, future research would evaluate the application of NOAA Climate Prediction Center (CPC) forecasts of soil moisture in predicting annual streamflow for use in the decision support model for the LCRA.

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.

Plurality of tree species responses to drought perturbation in Bornean tropical rain forest

Drought perturbation driven by the El Niño Southern Oscillation (ENSO) is a principal stochastic variable determining the dynamics of lowland rain forest in S.E. Asia. Mortality, recruitment and stem growth rates at Danum in Sabah (Malaysian Borneo) were recorded in two 4-ha plots (trees ≥ 10 cm gbh) for two periods, 1986–1996 and 1996–2001. Mortality and growth were also recorded in a sample of subplots for small trees (10 to <50 cm gbh) in two sub-periods, 1996–1999 and 1999–2001. Dynamics variables were employed to build indices of drought response for each of the 34 most abundant plot-level species (22 at the subplot level), these being interval-weighted percentage changes between periods and sub-periods. A significant yet complex effect of the strong 1997/1998 drought at the forest community level was shown by randomization procedures followed by multiple hypothesis testing. Despite a general resistance of the forest to drought, large and significant differences in short-term responses were apparent for several species. Using a diagrammatic form of stability analysis, different species showed immediate or lagged effects, high or low degrees of resilience or even oscillatory dynamics. In the context of the local topographic gradient, species’ responses define the newly termed perturbation response niche. The largest responses, particularly for recruitment and growth, were among the small trees, many of which are members of understorey taxa. The results bring with them a novel approach to understanding community dynamics: the kaleidoscopic complexity of idiosyncratic responses to stochastic perturbations suggests that plurality, rather than neutrality, of responses may be essential to understanding these tropical forests. The basis to the various responses lies with the mechanisms of tree-soil water relations which are physiologically predictable: the timing and intensity of the next drought, however, is not. To date, environmental stochasticity has been insufficiently incorporated into models of tropical forest dynamics, a step that might considerably improve the reality of theories about these globally important ecosystems.