Basin Scale Water Resources Systems Modeling Under Cascading Uncertainties

Global change in climate and consequent large impacts on regional hydrologic systems have, in recent years, motivated significant research efforts in water resources modeling under climate change. In an integrated future hydrologic scenario, it is likely that water availability and demands will change significantly due to modifications in hydro-climatic variables such as rainfall, reservoir inflows, temperature, net radiation, wind speed and humidity. An integrated regional water resources management model should capture the likely impacts of climate change on water demands and water availability along with uncertainties associated with climate change impacts and with management goals and objectives under non-stationary conditions. Uncertainties in an integrated regional water resources management model, accumulating from various stages of decision making include climate model and scenario uncertainty in the hydro-climatic impact assessment, uncertainty due to conflicting interests of the water users and uncertainty due to inherent variability of the reservoir inflows. This paper presents an integrated regional water resources management modeling approach considering uncertainties at various stages of decision making by an integration of a hydro-climatic variable projection model, a water demand quantification model, a water quantity management model and a water quality control model. Modeling tools of canonical correlation analysis, stochastic dynamic programming and fuzzy optimization are used in an integrated framework, in the approach presented here. The proposed modeling approach is demonstrated with the case study of the Bhadra Reservoir system in Karnataka, India.

Influence of Spatial Variability on Pavement Responses Using Latin Hypercube Sampling on Two-Dimensional Random Fields

Although uncertainties in material properties have been addressed in the design of flexible pavements, most current modeling techniques assume that pavement layers are homogeneous. The paper addresses the influence of the spatial variability of the resilient moduli of pavement layers by evaluating the effect of the variance and correlation length on the pavement responses to loading. The integration of the spatially varying log-normal random field with the finite-difference method has been achieved through an exponential autocorrelation function. The variation in the correlation length was found to have a marginal effect on the mean values of the critical strains and a noticeable effect on the standard deviation which decreases with decreases in correlation length. This reduction in the variance arises because of the spatial averaging phenomenon over the softer and stiffer zones generated because of spatial variability. The increase in the mean value of critical strains with decreasing correlation length, although minor, illustrates that pavement performance is adversely affected by the presence of spatially varying layers. The study also confirmed that the higher the variability in the pavement layer moduli, introduced through a higher value of coefficient of variation (COV), the higher the variability in the pavement response. The study concludes that ignoring spatial variability by modeling the pavement layers as homogeneous that have very short correlation lengths can result in the underestimation of the critical strains and thus an inaccurate assessment of the pavement performance. (C) 2014 American Society of Civil Engineers.

Modeling non-stationarity in intensity, duration and frequency of extreme rainfall over India

Significant changes are reported in extreme rainfall characteristics over India in recent studies though there are disagreements on the spatial uniformity and causes of trends. Based on recent theoretical advancements in the Extreme Value Theory (EVT), we analyze changes in extreme rainfall characteristics over India using a high-resolution daily gridded (1 degrees latitude x 1 degrees longitude) dataset. Intensity, duration and frequency of excess rain over a high threshold in the summer monsoon season are modeled by non-stationary distributions whose parameters vary with physical covariates like the El-Nino Southern Oscillation index (ENSO-index) which is an indicator of large-scale natural variability, global average temperature which is an indicator of human-induced global warming and local mean temperatures which possibly indicate more localized changes. Each non-stationary model considers one physical covariate and the best chosen statistical model at each rainfall grid gives the most significant physical driver for each extreme rainfall characteristic at that grid. Intensity, duration and frequency of extreme rainfall exhibit non-stationarity due to different drivers and no spatially uniform pattern is observed in the changes in them across the country. At most of the locations, duration of extreme rainfall spells is found to be stationary, while non-stationary associations between intensity and frequency and local changes in temperature are detected at a large number of locations. This study presents the first application of nonstationary statistical modeling of intensity, duration and frequency of extreme rainfall over India. The developed models are further used for rainfall frequency analysis to show changes in the 100-year extreme rainfall event. Our findings indicate the varying nature of each extreme rainfall characteristic and their drivers and emphasize the necessity of a comprehensive framework to assess resulting risks of precipitation induced flooding. (C) 2014 Elsevier B.V. All rights reserved.

Season- and depth-dependent variability of a demersal fish assemblage in a large fjord estuary (Puget Sound, Washington)

Fjord estuaries are common along the northeast Pacific coastline, but little information is available on fish assemblage structure and its spatiotemporal variability. Here, we examined changes in diversity metrics, species biomasses, and biomass spectra (the distribution of biomass across body size classes) over three seasons (fall, winter, summer) and at multiple depths (20 to 160 m) in Puget Sound, Washington, a deep and highly urbanized fjord estuary on the U.S. west coast. Our results indicate that this fish assemblage is dominated by cartilaginous species (spotted ratfish [Hydrolagus colliei] and spiny dogfish [Squalus acanthias]) and therefore differs fundamentally from fish assemblages found in shallower estuaries in the northeast Pacific. Diversity was greatest in shallow waters (<40 m), where the assemblage was composed primarily of flatfishes and sculpins, and lowest in deep waters (>80 m) that are more common in Puget Sound and that are dominated by spotted ratf ish and seasonally (fall and summer) by spiny dogfish. Strong depth-dependent variation in the demersal fish assemblage may be a general feature of deep fjord estuaries and indicates pronounced spatial variability in the food web. Future comparisons with less impacted fjords may offer insight into whether cartilaginous species naturally dominate these systems or only do so under conditions related to human-caused ecosystem degradation. Information on species distributions is critical for marine spatial planning and for modeling energy flows in coastal food webs. The data presented here will aid these endeavors and highlight areas for future research in this important yet understudied system.

Age validation, growth, mortality, and demographic modeling of spotted gully shark (Triakis megalopterus) from the southeast coast of South Africa

This study documents validation of vertebral band-pair formation in spotted gully shark (Triakis megalopterus) with the use of fluorochrome injection and tagging of captive and wild sharks over a 21-year period. Growth and mortality rates of T. megalopterus were also estimated and a demographic analysis of the species was conducted. Of the 23 OTC (oxytetracycline) -marked vertebrae examined (12 from captive and 11 from wild sharks), seven vertebrae (three from captive and four from wild sharks) exhibited chelation of the OTC and fluoresced under ultraviolet light. It was concluded that a single opaque and translucent band pair was deposited annually up to at least 25 years of age, the maximum age recorded. Reader precision was assessed by using an index of average percent error calculated at 5%. No significant differences were found between male and female growth patterns (P>0.05), and von Bertalanffy growth model parameters for combined sexes were estimated to be L∞=1711.07 mm TL, k=0.11/yr and t0=–2.43 yr (n=86). Natural mortality was estimated at 0.17/yr. Age at maturity was estimated at 11 years for males and 15 years for females. Results of the demographic analysis showed that the population, in the absence of fishing mortality, was stable and not significantly different from zero and particularly sensitive to overfishing. At the current age at first capture and natural mortality rate, the fishing mortality rate required to result in negative population growth was low at F>0.004/ yr. Elasticity analysis revealed that juvenile survival was the principal factor in explaining variability in population growth rate.

Ontogenetic and temporal variability in the fat content and fatty acid composition of Atlantic herring (Clupea harengus) from the Bay of Fundy, Canada

Atlantic herring (Clupea harengus) is an ecologically and economically valuable species in many food webs, yet surprisingly little is known about the variation in the nutritional quality of these fish. Atlantic herring collected from 2005 through 2008 from the Bay of Fundy, Canada, were examined for variability in their nutritional quality by using total lipid content (n=889) and fatty acid composition (n=551) as proxies for nutritional value. A significant positive relationship was found between fish length and total lipid content. Atlantic herring also had significantly different fatty acid signatures by age. Fish from 2005 had significantly lower total lipid content than fish from 2006 through 2008, and all years had significantly different fatty acid signatures. Summer fish were significantly fatter than winter fish and had significantly different fatty acid signatures. For all comparisons (ontogenetic, annual, and seasonal) percent concentrations of omega-3, -6, and long-chain monounsaturated fatty acids were the most important for distinguishing between the fatty acid signatures of fish. This study underscores the importance of quantifying variation in prey quality synoptically with prey quantity in food webs over ontogenetic and temporal scales when evaluating the effect of prey nutritional quality on predators and on modeling trophic dynamics.

Age validation, growth modeling, and mortality estimates for striped trumpeter (Latris lineata) from southeastern Australia: making the most of patchy data

Age estimates for striped trumpeter (Latris lineata) from Tasmanian waters were produced by counting annuli on the transverse section of sagittal otoliths and were validated by comparison of growth with known-age individuals and modal progression of a strong recruitment pulse. Estimated ages ranged from one to 43 years; fast growth rates were observed for the first five years. Minimal sexual dimorphism was shown to exist between length, weight, and growth characteristics of striped trumpeter. Seasonal growth variability was strong in individuals up to at least age four, and growth rates peaked approximately one month after the observed peak in sea surface temperature. A modified two-phase von Bertalanffy growth function was fitted to the length-at-age data, and the transition between growth phases was linked to apparent changes in physiological and life history traits, including offshore movement as fish approach maturity. The two-phase curve was found to represent the mean length at age in the data better than the standard von Bertalanffy growth function. Total mortality was estimated by using catch curve analysis based on the standard and two-phase von Bertalanffy growth functions, and estimates of natural mortality were calculated by using two empirical models, one based on longevity and the other based on the parameters L∞ and k from both growth functions. The interactions between an inshore gillnet fishery targeting predominately juveniles and an offshore hook fishery targeting predominately adults highlight the need to use a precautionary approach when developing harvest strategies.

Modeling North Pacific temperature and pressure changes from coastal tree-ring chronologies

Climate modeling using coastal tree-ring chronologies has yielded the first summer temperature reconstructions for coastal stations along the Gulf of Alaska and the Pacific Northwest. These land temperature reconstructions are strongly correlated with nearby sea surface temperatures, indicating large-scale ocean-atmospheric influences. Significant progress has also been made in modeling winter land temperatures and sea surface temperatures from coastal and shipboard stations. In addition to temperature, the pressure variability center over the central North Pacific Ocean (PAC), which is related to the strength and location of the Aleutian Low pressure system, could be extended using coastal tree rings.

Variability of nanoindentation responses of bone and artificial bone-like composites

Nanoindentation is a popular technique for measuring the intrinsic mechanical response of bone and has been used to measure a single-valued elastic modulus. However, bone is a composite material with 20-80 nm hydroxyapatite plates embedded in a collagen matrix, and modern instrumentation allows for measurements at these small length scales. The present study examines the indentation response of bone and artificial gelatin-apatite nanocomposite materials across three orders of magnitude of lengthscale, from nanometers to micrometers, to isolate the composite phase contributions to the overall response. The load-displacement responses were variable and deviated from the quadratic response of homogeneous materials at small depths. The distribution of apparent elastic modulus values narrowed substantially with increasing indentation load. Indentation of particulate nanocomposites was simulated using finite element analysis. Modeling results replicated the convergence in effective modulus seen in the experiments. It appears that the apatite particles are acting as the continuous ("matrix") phase in bone and nanocomposites. Copyright © 2004 by ASME.

Sediment accumulation rates in subarctic lakes: insights into age-depth modeling from 22 dated lake records from the Northwest Territories, Canada

Age-depth modeling using Bayesian statistics requires well-informed prior information about the behavior of sediment accumulation. Here we present average sediment accumulation rates (represented as deposition times, DT, in yr/cm) for lakes in an Arctic setting, and we examine the variability across space (intra- and inter-lake) and time (late Holocene). The dataset includes over 100 radiocarbon dates, primarily on bulk sediment, from 22 sediment cores obtained from 18 lakes spanning the boreal to tundra ecotone gradients in subarctic Canada. There are four to twenty-five radiocarbon dates per core, depending on the length and character of the sediment records. Deposition times were calculated at 100-year intervals from age-depth models constructed using the ‘classical’ age-depth modeling software Clam. Lakes in boreal settings have the most rapid accumulation (mean DT 20 ± 10 years), whereas lakes in tundra settings accumulate at moderate (mean DT 70 ± 10 years) to very slow rates, (>100 yr/cm). Many of the age-depth models demonstrate fluctuations in accumulation that coincide with lake evolution and post-glacial climate change. Ten of our sediment cores yielded sediments as old as c. 9,000 cal BP (BP = years before AD 1950). From between c. 9,000 cal BP and c. 6,000 cal BP, sediment accumulation was relatively rapid (DT of 20 to 60 yr/cm). Accumulation slowed between c. 5,500 and c. 4,000 cal BP as vegetation expanded northward in response to warming. A short period of rapid accumulation occurred near 1,200 cal BP at three lakes. Our research will help inform priors in Bayesian age modeling.

Nong Thale Pron - a key site from southern Thailand for studying monsoon variability during the past 15000 years

Studies of marine sediments, cave speleothemes, annually laminated corals, and tree rings from Asian monsoon regions have added knowledge to our understanding of the factors that control inter-annual to millennial monsoon variability in the past and have provided important constraints for climate modeling scenarios. In contrast, the spatial and temporal pattern of sub-millennial scale monsoon variability and its impact on land cover in SE Asia are still unresolved. This shortcoming stems from the fact that temporally well-resolved paleo-environmental studies are missing from large parts of SE Asia, especially from Thailand. Given that global and regional climate models are increasingly using terrestrial paleo- data to test their performance, past changes in land cover are therefore important variables to better understand feedbacks between different Earth systems. We obtained sediments from Lake Nong Thale Pron, in southern Thailand (8º 10`N, 99 º23`E; 380 m.asl). The aim of our study is to reconstruct lake status changes and to evaluate whether the extent of these changes are linked to known shifts in monsoon intensity and variability. Preliminary results show that lake infilling started more than 15,000 years ago and that the sediments cover the last deglaciation and the Holocene. Current analyses include Itrax XRF core scanning, loss-on-ignition (LOI at 950 and 550ºC), CN elemental and isotopic composition. We expect that our results will be able to give a picture of how the lake's status has changed over time and whether the extent of these changes is linked to known shifts in monsoon intensity and variability.

Modeling the Western Iberian margin circulation: present and future

The purpose of this work is to carry out a comprehensive study on the Western Iberian Margin (WIM) circulation my means of numerical modeling, and to postulate what this circulation will be in the future. The adopted approach was the development of a regional ocean model configuration with high resolution, capable of reproducing the largeand small-scale dynamics of the coastal transition zone. Four numerical experiences were carried out according to these objectives: (1) a climatological run, in order to study the system’s seasonal behavior and its mean state; (2) a run forced with real winds and fluxes for period 2001-2011 in order to study the interannual variability of the system; (3) a run forced with mean fields from Global Climate Models (GCMs) for the present, in order to validate GCMs as adequate forcing for regional ocean modeling; (4) a similar run (3) for period 2071-2100, in order to assess possible consequences of a future climate scenario on the hydrography and dynamics of the WIM. Furthermore, two Lagrangian particle studies were carried out: one in order to trace the origin of the upwelled waters along the WIM; the other in order to portrait the patterns of larval dispersal, accumulation and connectivity. The numerical configuration proved to be adequate in the reproduction of the system’s mean state, seasonal characterization and an interannual variability study. There is prevalence of poleward flow at the slope, which coexists with the upwelling jet during summer, although there is evidence of its shifting offshore, and which is associated with the Mediterranean Water flow at deeper levels, suggesting a barotropic character. From the future climate scenario essay, the following conclusions were drawn: there is general warming and freshening of upper level waters; there is still poleward tendency, and despite the upwellingfavorable winds strengthening in summer the respective coastal band becomes more restricted in width and depth. In what concerns larval connectivity and dispersion along the WIM, diel vertical migration was observed to increase recruitment throughout the domain, and while smooth coastlines are better suppliers, there is higher accumulation where the topography is rougher.

DPSIR framework application supported by integrated watershed modeling for irrigation from Alqueva dam (Portugal) on Guadiana river estuary

Dissertação de mest., Gestão da Água e da Costa, Faculdade de Ciências e Tecnologia, Universidade do Algarve, 2009

Dynamics and variability of the alongshore flows on the northwestern iberian margin

Tese de doutoramento, Ciências Geofísicas e da Geoinformação (Oceanografia), Universidade de Lisboa, Faculdade de Ciências, 2014