ALBERT

Description: Although catchment storage is an intrinsic control on the rainfall-runoff response of streams, direct measurement remains a major challenge. Coupled models that integrate long-term hydrometric and isotope tracer data are useful tools that can provide insights into the dynamics of catchment storage and the volumes of water involved. In this study, we use a tracer-aided hydrological model to characterize catchment storage as a dynamic control on system function related to streamflow generation, which also allows direct estimation of the non-stationarity of water ages. We show that in a wet Scottish upland catchment dominated by runoff generation from riparian peats (histosols) with high water storage, non-stationarity in water age distributions are only clearly detectable during more extreme wet and dry periods. This is explained by the frequency and longevity of hydrological connectivity and the associated relative importance of flow paths contributing younger or older waters to the stream. Generally, these saturated riparian soils represent large mixing zones that buffer the time variance of water age and integrate catchment-scale partial mixing processes. Although storage simulations depend on model performance, which is influenced by input variability and the degree of isotopic damping in the stream, a longer-term storage analysis of this model indicates a system which is only sensitive to more extreme hydroclimatic variability. This article is protected by copyright. All rights reserved.

Description: ABSTRACT A reactive transport modeling framework is presented that allows simultaneous assessment of groundwater flow, water quality evolution including δ 13 C, and 14 C activity or “age”. Through application of this framework, simulated 14 C activities can be directly compared with measured 14 C activities. This bypasses the need for interpretation of a 14 C age prior to flow simulation through factoring out processes other than radioactive decay, which typically involves simplifying assumptions regarding spatial and temporal variability in reactions, flow, and mixing. The utility of the approach is demonstrated for an aquifer system with spatially variable carbonate mineral distribution, multiple organic carbon sources, and transient boundary conditions for 14 C activity in the recharge water. In this case the simulated 14 C age was shown to be relatively insensitive to isotopic fractionation during DOC oxidation and variations in assumed DOC degradation behaviour. We demonstrate that the model allows quantitative testing of hypotheses regarding controls on groundwater age and water quality evolution for all three carbon isotopes. The approach also facilitates incorporation of multiple environmental tracers and combination with parameter optimization techniques. This article is protected by copyright. All rights reserved.

Description: Water scarcity is likely to increase in the coming years, making improvements in irrigation efficiency increasingly important. An emerging technology that promises to increase irrigation efficiency substantially is a wireless irrigation sensor network that uploads sensor data into irrigation management software, creating an integrated system that allows real-time monitoring and control of moisture status that has been shown in experimental settings to reduce irrigation costs, lower plant loss rates, shorten production times, decrease pesticide application, and increase yield, quality, and profit. We use an original survey to investigate likely initial acceptance, ceiling adoption rates, and profitability of this new sensor network technology in the nursery and greenhouse industry. We find that adoption rates for a base system and demand for expansion components are decreasing in price, as expected. The price elasticity of the probability of adoption suggests that sensor networks are likely to diffuse at a rate somewhat greater than that of drip irrigation. Adoption rates for a base system and demand for expansion components are increasing in specialization in ornamental production: Growers earning greater shares of revenue from greenhouse and nursery operations are willing to pay more for a base system and are willing to purchase larger numbers of expansion components at any given price. We estimate that growers who are willing to purchase a sensor network expect investment in this technology to generate significant profit, consistent with findings from experimental studies. This article is protected by copyright. All rights reserved.

Description: Glass–ceramic tubes were developed from Bayan Obo tailings and fly ash by the centrifugal method. This paper investigates the corrosion rate of the glass–ceramics exposed to hydrochloric acid solution. The leaching concentration of cations into the solution is analyzed by ICP. XRD analysis and Raman spectra were used for identifying the phase transformation of corrosion reactions. The changes in morphological and chemical composition were investigated using SEM combined with EDS. The results revealed that the formed corrosion layer was amorphous phase. It was also indicated that Fe 2 O 3 and CaF 2 exerted some negative effect on the corrosion resistance in HCl solution.

Description: Herein, we report on the crystal structures of Nb 2 AlC and TiNbAlC—actual composition (Ti 0.45 ,Nb 0.55 ) 2 AlC—compounds determined from Rietveld analysis of neutron diffraction patterns in the 300–1173 K temperature range. The average linear thermal expansion coefficients of a Nb 2 AlC sample in the a and c directions are, respectively, 7.9(5) × 10 −6 and 7.7(5) × 10 −6  K −1 on one neutron diffractometer and 7.3(3) × 10 −6 and 7.0(2) × 10 −6  K −1 on a second diffractometer. The respective values for the (Ti 0.45 ,Nb 0.55 ) 2 AlC composition—only tested on one diffractometer—are 8.5(3) × 10 −6 and 7.5(5) × 10 −6  K −1 . These values are relatively low compared to other MAX phases. Like other MAX phases, however, the atomic displacement parameters (APDs) show that the Al atoms vibrate with higher amplitudes than the Ti and C atoms, and more along the basal planes than normal to them. When the predictions of the APDs obtained from density functional theory are compared to the experimental results, good quantitative agreement is found for the Al atoms. In case of the Nb and C atoms, the agreement was more qualitative.

Description: In this study, lithium disilicate (LS 2 ) glass samples with different particle sizes ranging from less than 105 to 850 μm were prepared. These specimens were inserted in a Pt-Rh DSC crucible and heated to 850°C at different rates (ϕ = 0.5–30 K/min) to identify their crystallization peaks. The activation energies for the overall crystallization ( E ) and the Avrami coefficient ( n ) were evaluated using different nonisothermal models. Specifically, n was evaluated using the Augis–Benett model and the Ozawa method, and E was evaluated using the Kissinger and Ligero methods. As expected, the coarse particles mainly crystallized in the volume, while surface crystallization was predominant in the samples with particle sizes of less than 350 μm. This result was confirmed through SEM analysis of the double stage heat-treated samples. In contrast with previous studies, our results demonstrated that the activation energy decreased as the particle size increased. In addition, no clear correlation between the peak intensity (δ T p ) and the particle size was observed.

Description: Selective Laser Melting (SLM) is a candidate for on-site manufacturing as its characteristics of energy source and powder-based fabrication process are suitable for use with in situ material. The feasibility of using lunar regolith simulant to create objects with SLM process is investigated in this study. The process parameters are optimized and multiple objects are fabricated. A qualitative chemical analysis is carried out with scanning electron microscopy using energy-dispersive X-ray emission. Lastly, properties such as particle size distribution, particle shape, and crystal structure of the lunar simulant powder as well as the crystallinity and hardness of the fabricated objects are investigated.

Description: This technical note presents a useful methodology for studying how the variance of hydraulic and/or reactive attributes of an aquifer are linked to the multi-scaled and hierarchical sedimentary architecture of the aquifer. A new recursive equation is derived which quantitatively describes how the variance is related to sedimentary facies defined at all scales across an entire stratal hierarchy. As compared to prior published equations that emphasize differences in means among facies populations within a hierarchical level, it emphasizes differences across levels. Because of the hierarchical relationships among the terms of the equation, we find it to be useful for conducting a holistic analysis of the relative contributions to the variance arising from all facies types defined across all scales. The methodology is demonstrated using appropriate field data, and is shown to be useful in defining parsimonious classification systems.

Description: ABSTRACT A primary concern for geologic carbon storage is the potential for leakage of stored carbon dioxide (CO 2 ) into the shallow subsurface where it could degrade the quality of groundwater and surface water. In order to predict and mitigate the potentially negative impacts of CO 2 leakage, it is important to understand the physical processes that CO 2 will undergo as it moves through naturally heterogeneous porous media formations. Previous studies have shown that heterogeneity can enhance the evolution of gas phase CO 2 in some cases, but the conditions under which this occurs have not yet been quantitatively defined, nor tested through laboratory experiments. This study quantitatively investigates the effects of geologic heterogeneity on the process of gas phase CO 2 evolution in shallow aquifers through an extensive set of experiments conducted in a column that was packed with layers of various test sands. Soil moisture sensors were utilized to observe the formation of gas phase near the porous media interfaces. Results indicate that the conditions under which heterogeneity controls gas phase evolution can be successfully predicted through analysis of simple parameters, including the dissolved CO 2 concentration in the flowing water, the distance between the heterogeneity and the leakage location, and some fundamental properties of the porous media. Results also show that interfaces where a less permeable material overlies a more permeable material affect gas phase evolution more significantly than interfaces with the opposite layering.

Description: Fluvial sediment loads are frequently calculated with rating curves fit to measured sediment transport rates. Rating curves are often treated as statistical representations in which the fitted parameters have little or no physical meaning. Such models, however, may produce large errors when extrapolation is needed, and they provide no insight into the sediment transport process. It is shown that log-linear least squares, the usual method for fitting rating curves, does not generally produce physically meaningful parameter values. In addition, it cannot accommodate data that include zero-transport samples. Alternative fitting methods based non-linear least squares and on maximum likelihood parameter estimation are described and evaluated. The maximum likelihood approach is shown to fit synthetic data better than linear or non-linear least squares, and to perform well with data that include zero-transport samples. In contrast, non-linear least squares methods produce large errors in the parameter estimates when zero-transport samples are present or when the variance structure of the data is incorrectly specified. Analyses with fractional bedload data from a mountain stream suggest that bedload transport rates are gamma distributed, that the arrivals of bedload particles in a sampler conform to a Poisson distribution, and that the variance of non-zero samples can be expressed as a power function of the mean. Preliminary physical interpretations of variations in the rating curve parameters fit to fractional bedload data with the maximum likelihood method are proposed, and their relation to some previous interpretations of rating curve parameters are briefly discussed. This article is protected by copyright. All rights reserved.

Description: SrLn 2 Al 2 O 7 (Ln = La, Nd, Sm) ceramics belonging to Ruddlesden–Popper phase with n  =   2 in space group I 4 /mmm have been prepared by a standard solid-state reaction process. Their crystal structure refinement and infrared reflection spectra analysis have been carried out to investigate the correlations among the crystal structure, polar phonon modes, and microwave dielectric characteristics. Heavy deviations of cation–oxygen bonds in (Sr,Ln)O 12 , (Sr,Ln)O 9 , and AlO 6 polyhedra are observed, and the deviations significantly affect the microwave dielectric characteristics of the present ceramics. The polar phonon modes with lower frequencies due to translation vibration of Sr/La ions against the AlO 6 octahedron give the primary contribution to dielectric constant and intrinsic dielectric loss at microwave frequencies. The calculated Q  ×  f value is significantly larger than the measured one, and this result is very obvious for SrLa 2 Al 2 O 7 ceramics where the calculated one (130,300 GHz) is about 2 times of the measured one (68,800 GHz).

Description: In this letter, 2-μm Pb 0.97 La 0.02 (Zr 0.75 Sn 0.18 Ti 0.07 )O 3 antiferroelectric thick film with tetragonal structure was prepared. The effects of operating electric field, temperature, and frequency on the thermal–electrical energy harvesting capacity of the film were studied by using the Olsen cycle. The results demonstrated that giant energy harvesting effect could be realized in the antiferroelectric thick film. The maximum harvestable energy density per cycle of the film was about 7.8 J/cm 3 at 1 kHz, which was the largest reported value to date. The corresponding energy harvesting efficiency was 0.53%. Moreover, the film had a low leakage current density (about 7.3 × 10 −7 and 3.9 × 10 −5  A/cm 2 at 25 and 200°C, respectively), which was favorable for its application in the devices of the thermal–electrical energy harvesting.

Description: Discharged energy properties of PbO–SrO–Na 2 O–Nb 2 O 5 –SiO 2 glass-ceramics with crystallization time from 1 to 1000 min were investigated by measuring their hysteresis loops (described as quasi-static measuring method) and pulse-discharge current-time curves (described as dynamic measuring method). The results show the same trend for both measuring methods: With the increment of crystallization time, the discharged energy density increases gradually, while the energy efficiency decreases. The highest energy efficiencies were obtained in the sample with crystallization time of 1 min, which are 96.3% and 82.4%, corresponding quasi-static and dynamic measurement, respectively. The reduction of energy efficiency with crystallization time is attributed to combined effect of ferroelectric polarization and interfacial polarization, and part of the corresponding energy could not release in the pulse-discharge process.

Description: A Ba 0.6 Sr 0.4 TiO 3 –ZnO–B 2 O 3 composite ink was prepared and used for the manufacturing of fully inkjet-printed metal-insulator-metal varactors. The dielectric thick films were co-fired with printed silver electrodes at 850°C and show a fine grained microstructure. The films have a relative permittivity of ε r  = 129 and a dielectric loss of tan δ = 0.043 at f  =   3 GHz. Printed varactors with different dielectric film thickness were prepared. The characterization of the printed structures and 3D electromagnetic simulations of the layout reveal the strong influence of electrode inductance and fringing effects on the properties of the components. The printed varactors reach tunabilities between 14.4% and 16.4% by applying a tuning field of 5 V/μm. To demonstrate the capability of the inkjet printing process for the preparation of tunable microwave devices, a fully inkjet-printed phase shifter was fabricated. It is based on seven pairs of the printed varactors and reaches a phase shift of 180° and a figure of merit of 19°/dB at 3.4 GHz.

Description: New lead-free perovskite solid solution ceramics of (1  − x )( Bi 1/2 Na 1/2 ) TiO 3 – x Ba ( Ni 1/2 Nb 1/2 ) O 3 [(1− x )BNT– x BNN, x  =   0.02–0.06) were prepared and their dielectric, ferroelectric, piezoelectric, and electromechanical properties were investigated as a function of the BNN content. The X-ray diffraction results indicated that the addition of BNN has induced a morphotropic phase transformation from rhombohedral to pseudocubic symmetry approximately at x  =   0.045, accompanying an evolution of dielectric relaxor behavior as characterized by enhanced dielectric diffuseness and frequency dispersion. In the proximity of the ferroelectric rhombohedral and pseudocubic phase coexistence zone, the x  =   0.045 ceramics exhibited optimal piezoelectric and electromechanical coupling properties of d 33 ~121 pC/N and k p ~0.27 owing to decreased energy barriers for polarization switching. However, further addition of BNN could cause a decrease in freezing temperatures of polar nanoregions till the coexistence of nonergodic and ergodic relaxor phases occurred near room temperature, especially for the x  =   0.05 sample which has negligible negative strains and thus show the maximum electrostrain of 0.3% under an external electric field of 7 kV/mm, but almost vanished piezoelectric properties. This was attributed to the fact that the induced long-range ferroelectric order could reversibly switch back to its original ergodic state upon removal of external electric fields.

Description: Light transmission in polycrystalline magnesium fluoride was studied as a function of the mean grain size at different wavelengths. The mean grain size was varied by annealing hot-pressed billets in argon atmosphere at temperatures ranging from 600°C to 800°C for 1 h. The grain-size and grain-orientation distributions were characterized by electron back scatter diffraction. The scattering coefficients were calculated from the in-line transmittance measured at various wavelengths. The scattering coefficient of polycrystalline magnesium fluoride increased linearly with the mean grain size and inversely with the square of the wavelength of light. It is shown that these trends are consistent with theoretical models based on both a limiting form of the Raleigh–Gans–Debye (RGD) theory of particle scattering and light retardation theories that take refractive index variations along the light path. Quantitative predictions of the theories are, however, subject to uncertainly due to the restrictive assumptions made in the theories and difficulties in representing the microstructure in the theoretical models. In particular, grain-size distribution has a significant influence on the scattering coefficient calculated using particle scattering models.

Description: Determining the absolute chronology of ceramic artifacts has significant implications for archeological and historical research. Wilson, Hall et al . recently suggested a new technique for direct absolute dating of archeological ceramics based on a moisture-induced chemical reaction, called rehydroxylation (RHX) dating. RHX dating proceeds by measuring the mass of chemically combined water in the ceramics in the form of OH hydroxyls, and the mass gain rate at the Effective Lifetime Temperature (ELT) that the ceramics experienced over its lifetime. To date, ELT determinations have been based on estimates of the ceramic's lifetime temperature history; taking into account weather and climate data and the depth at which the artifact was found. The uncertainty in determining the ELT can be a major component of the overall dating uncertainty. Here, we propose an alternative method which relies minimally on weather and climate data, and provides more precise determinations of the ELT and the ceramic age . The proposed method (SAS: Same Age Samples) involves a minimum of four measurements of the RHX mass gain rate constant for two ceramic samples of the same age at two temperatures. We show via simulations that the proposed SAS method can determine the ELT with a precision of 0.2 K which is comparable to the best ELT determination based on lifetime temperature history, and also comparable to available microbalance temperature resolutions of around 0.1 K. The corresponding percent age error is then 1.4%, or 43 yr for a 3000-yr-old ceramic. The proposed SAS method should be tested with ceramic samples of different ages, whose ELT are well-known.

Description: Glasses in the Na 2 O–CaO–SrO–ZnO–SiO 2 system have previously been investigated for suitability as a reagent in Al-free glass polyalkenoate cements (GPCs). These materials have many properties that offer potential in orthopedics. However, their applicability has been limited, to date, because of their poor strength. This study was undertaken with the aim of increasing the mechanical properties of a series of these Zn-based GPC glasses by doping with nitrogen to give overall compositions of: 10Na 2 O–10CaO–20SrO–20ZnO–(40−3 x )SiO 2 – x Si 3 N 4 ( x is the no. of moles of Si 3 N 4 ). The density, glass-transition temperature, hardness, and elastic modulus of each glass were found to increase fairly linearly with nitrogen content. Indentation fracture resistance also increases with nitrogen content according to a power law relationship. These increases are consistent with the incorporation of N into the glass structure in threefold coordination with silicon resulting in extra cross-linking of the glass network. This was confirmed using 29 Si MAS-NMR which showed that an increasing number of Q 2 units and some Q 3 units with extra bridging anions are formed as nitrogen content increases at the expense of Q 1 units. A small proportion of Zn ions are found to be in tetrahedral coordination in the base oxide glass and the proportion of these increases with the presence of nitrogen.

Description: Hercynite nanoparticles were obtained by pulsed laser ablation of hercynite in liquid media technique. Colloidal solutions of hercynite nanoparticles were synthesized in three different liquids by ablation using 532-nm output of a pulsed Nd:YAG laser with three laser fluencies. Transmission electron microscopy, selected area electron diffraction, energy-dispersive X-ray analysis, and X-ray photoelectron spectroscopy were used for the characterization of these nanoparticles. The nanoparticles obtained were of spherical shapes with cubic crystalline structure. It was observed that the laser fluence and liquid media could affect the morphology and size of hercynite nanoparticles.

Description: This paper presents the results of a comprehensive model-based analysis of a uranyl [U(VI)] tracer test conducted at the U.S. DOE Hanford 300 Area (300A) IFRC. Despite the highly complex field conditions the numerical three-dimensional multi-component reactive transport model was able to capture most of the spatiotemporal variations of the observed U(VI) concentrations. A multi-model analysis was performed to interrogate the relative importance of various processes and factors for controlling field-scale reactive transport during the uranyl tracer test. The results indicate that multi-rate sorption/desorption, surface complexation reactions, and initial concentrations were the most important processes and factors controlling U(VI) migration. On the other hand, cation exchange reactions, the choice of the surface complexation model, and dual-domain mass transfer processes played less important roles under the prevailing field-test condition. Further analysis of the modeling results demonstrates that these findings are conditioned to the relatively stable groundwater chemistry and the selected length of the field experimental duration (16 days). The model analysis also revealed the crucial role of the intraborehole flow that occurred within the long-screened monitoring wells and thus affected both field measurements and simulated U(VI) concentrations as a combined effect of aquifer heterogeneity and dynamic flow conditions. This study provides the first highly data-constrained uranium transport simulations under highly dynamic flow conditions. It illustrates the value of reactive transport modeling for elucidating the relative importance of individual processes in controlling uranium transport under specific field-scale conditions.

Description: The dynamics of drying processes from porous media are critically influenced by the intensity of an adjacent free flow and by processes at the interface between free flow and the porous medium. In this paper, the influence of hydraulic properties of a porous medium and of the interaction between fluids and porous medium on the drying dynamics during the capillary-flow dominated stage-1 and transition to the diffusion-dominated stage-2 are studied using a coupled free-flow - porous-medium flow model on the REV scale. We present a detailed model concept that considers mass balance equations, an energy balance equation and the coupling to the adjacent free flow. Key microscale processes are identified and incorporated in the macroscale description of the evaporation process. Own experimental results are used to illustrate main features of the modeling framework. We demonstrate that the use of a homogeneous distribution of soil parameters without consideration of pore-scale induced nonlinearities in the numerical simulations results in a rather constant drying rate in stage-1, which was not observed for the high evaporative demand in the experiments. To account for the dependency of the drying rate on the surface moisture content, special conditions based on the work of Haghighi et al. [2013] and Schlünder [1988] are analyzed for their applicability on the REV scale. Typical features of a drying process, such as different stages of the drying rate, could be reproduced.

Description: Making useful predictions in ungauged basins is an incredibly difficult task given the limitations of hydrologic models to represent physical processes appropriately across the heterogeneity within and among different catchments. Here, we introduce a new method for this challenge, Bayes empirical Bayes, that allows for the statistical pooling of information from multiple donor catchments and provides the ability to transfer parametric distributions rather than single parameter sets to the ungauged catchment. Further, the methodology provides an efficient framework with which to formally assess predictive uncertainty at the ungauged catchment. We investigated the utility of the methodology under both synthetic and real data conditions, and with respect to its sensitivity to the number and quality of the donor catchments used. This study highlighted the ability of the hierarchical Bayes empirical Bayes approach to produce expected outcomes in both the synthetic and real data applications. The method was found to be sensitive to the quality (hydrologic similarity) of the donor catchments used. Results were less sensitive to the number of donor catchments, but indicated that predictive uncertainty was best constrained with larger numbers of donor catchments (but still adequate with fewer donors)

Description: Autumn is a season of dynamic change in forest streams of the northeastern USA due to effects of leaf fall on both hydrology and biogeochemistry. Few studies have explored how interactions of biogeochemical transformations, various nitrogen sources, and catchment flowpaths affect stream nitrogen variation during autumn. To provide more information on this critical period, we studied 1) the timing, duration, and magnitude of changes to stream nitrate, dissolved organic nitrogen (DON), and ammonium concentrations; 2) changes in nitrate sources and cycling; and 3) source areas of the landscape that most influence stream nitrogen. We collected samples at higher temporal resolution for a longer duration than typical studies of stream nitrogen during autumn. This sampling scheme encompassed the patterns and extremes that occurred during baseflow and stormflow events of autumn. Baseflow nitrate concentrations decreased by an order of magnitude from 5.4 to 0.7 μmol L -1 during the week when most leaves fell from deciduous trees. Changes to rates of biogeochemical transformations during autumn baseflow explained the low nitrate concentrations; in-stream transformations retained up to 72% of the nitrate that entered a stream reach. A decrease of in-stream nitrification coupled with assimilatory nitrate uptake was a primary factor in the seasonal nitrate decline. The period of low nitrate concentrations ended with a storm event in which stream nitrate concentrations increased by 25 fold. In the ensuing weeks, stormflow nitrate concentrations progressively decreased over closely-spaced, yet similarly sized events. Most stormflow nitrate originated from nitrification in near-stream areas with occasional, large inputs of unprocessed atmospheric nitrate, which has rarely been reported for non-snowmelt events. A maximum input of 33% unprocessed atmospheric nitrate to the stream occurred during one event. The large inputs of unprocessed atmospheric nitrate show direct and rapid effects on forest streams that may be widespread, although undocumented, throughout nitrogen-polluted temperate forests. In contrast to a week-long nitrate decline during peak autumn litterfall, baseflow DON concentrations increased after leaf fall and remained high for two months. Dissolved organic nitrogen was hydrologically flushed to the stream from riparian soils during stormflow. In contrast to distinct seasonal changes in baseflow nitrate and DON concentrations, ammonium concentrations were typically at or below detection limit, similar to the rest of the year. Our findings reveal couplings among catchment flow paths, nutrient sources and transformations that control seasonal extremes of stream nitrogen in forested landscapes.

Description: During the recent years there has been an increasing interest in multivariate frequency analysis of hydrological variables, e.g. those describing extreme events like rainfall, floods or droughts. The multivariate analysis provides a better understanding of the phenomena under investigation and an additional insight about the interrelationships between the different variables (e.g. peak, volume and duration of the flood), exploiting the complete structure of the problem and making a full use of the available data. However, while the developments on multivariate analysis of hydrological data has produced a large body of literature, a clear assessment of the use of these methods in the design and risk assessment of hydraulic structures is still a matter of debate. In the present work we illustrate a general, structure-based framework for the design and/or risk assessment of hydraulic structures in a bivariate environment; we also compare it to recently proposed methods which are based on the assumption of hydrological design events (as is customary in the univariate context). For illustration purposes, both the structure-based and the design event-based approaches are applied to the design of an idealized structure, thus exploring the differences among the methods as function of the parameters involved. Our work highlights that the return period of structure failure in a multivariate environment strictly depends on the particular structure under design, and in most cases the design of an hydraulic structure cannot be based on a single, hydrological multivariate design event. This acts as a warning for practitioners against the use of design methods based on single hydrological events, as usually done in the context of univariate hydrology, thus neglecting the interplay between the structure and the hydrological loads acting on it.

Description: Phosphorus (P) is one of the major limiting nutrient in many freshwater ecosystems. During the last decade, attention has been focused on the fluxes of suspended sediment and particulate P through freshwater drainage systems because of severe eutrophication effects in aquatic ecosystems. Hence, the analysis and prediction of phosphorus and sediment dynamics constitutes an important element for ecological conservation and restoration of freshwater ecosystems. In that sense, the development ofa suitable prediction model is justified and the present work is devoted to the validation and application of a predictive Soluble Reactive Phosphorus (SRP) uptake and sedimentation models, to a real riparian system of the middle Ebro river floodplain. Both models are coupled to a fully distributed 2D shallow water flow numerical model. The SRP uptake model is validated using data from three field experiments. The model predictions show a good accuracy for SRP concentration, where the linear regressions between measured and calculated of the three experiments were significant ( r 2  ≥ 0.62; p  ≤ 0.05), and a Nash-Sutcliffe coefficient (E) that ranged from 0.54 to 0.62. The sedimentation model is validated using field data collected during 2 real flooding events within the same river reach. The comparison between calculated and measured sediment deposition showed a significant linear regression ( p  ≤ 0.05; r 2  = 0.97) and an E that ranged from 0.63 to 0.78. Subsequently, the complete model that includes flow dynamics, solute transport, SRP uptake and sedimentation is used to simulate and analyze floodplain sediment deposition, river nutrient contribution and SRP uptake. According to this analysis, the main SRP uptake process appears to be the sediment sorption. The analysis also reveals the presence of a lateral gradient of hydrological connectivity that decreases with distance from the river, and controls the river matter contribution to the floodplain. This article is protected by copyright. All rights reserved.

Description: The inherent effects of global sea surface temperature (SST) anomalies on hydrological cycle and vegetation cover complicate the structure of tropical climate at the regional scale. Assessing hydrological processes related to climate forcing is important in Central America because it is surrounded by both the Pacific and Atlantic oceans and two continental landmasses. In this study, the use of high-resolution remote sensing imagery in wavelet analysis helps identify nonstationary characteristics of hydrological and ecological responses. The wavelet-based empirical orthogonal function (WEOF) further reflects the nonlinear relationship between the Atlantic and Pacific SST and the greenness of a pristine forested site in Panama, La Amistad International Park. Integrated WEOF and descriptive statistics for data analysis reveal a higher temporal variability in terrestrial precipitation relative to in-situ land surface temperature and its probable effects on the presence of dry periods. Such teleconnection signals of SST were identified as a driving force of decline in tropical forest greenness during dry periods. The results of our remote sensing-based wavelet analysis showed intrannual high frequency and biennial to triennial low frequency signals between enhanced vegetation index (EVI)/precipitation datasets and SST indices in both Atlantic and Pacific oceans. A spatiotemporal priority search further confirmed the importance of the effects of the El Niño-Southern Oscillation (ENSO) over terrestrial responses in the selected study site. Coincidence of the effect of ENSO teleconnection patterns on precipitation and vegetation suggest possible impacts of El Niño-associated droughts in Central America, accompanied by reduced rainfall, especially during the first months of rainy season (June, July and August), and decline in vegetation cover during the dry season (March and April). This article is protected by copyright. All rights reserved.

Description: Proliferation of evapotranspiration (ET) products warrants comparison of these products. The study objective was to assess uncertainty in ET output from four land surface models (LSMs), Noah, Mosaic, VIC, and SAC in NLDAS-2, two remote sensing-based products, MODIS and AVHRR, and GRACE-inferred ET from a water budget with precipitation from PRISM, monitored runoff, and total water storage change (TWSC) from GRACE satellites. The three cornered hat method, which does not require a priori knowledge of the true ET value, was used to estimate ET uncertainties. In addition, TWSC or total water storage anomaly (TWSA) from GRACE was compared with water budget estimates of TWSC from a flux-based approach or TWSA from a storage-based approach. The analyses were conducted using data from three regions (humid – arid) in the South Central US as case studies. Uncertainties in ET are lowest in LSM ET (~5 mm/month), moderate in MODIS- or AVHRR-based ET (10 – 15 mm/month), and highest in GRACE-inferred ET (20 – 30 mm/month). There is a tradeoff between spatial resolution and uncertainty, with lower uncertainty in the coarser-resolution LSM ET (~14 km) relative to higher uncertainty in the finer-resolution (~ 1 ‒ 8 km) RS ET. Root-mean-square (RMS) of uncertainties in water budget estimates of TWSC is about half of RMS of uncertainties in GRACE-derived TWSC for each of the regions. Future ET estimation should consider a hybrid approach that integrates strengths of LSMs and satellite-based products to constrain uncertainties.

Description: It is well documented that deforestation results in an increase in landslide frequency due to the control that forest roots have on slope stability. The loss of forest vegetation leads to a reduction in soil cohesion and a decrease in the shear strength of the soil profile. As a result, the slope becomes more susceptible to landsliding and the return time of landslides decreases. When a landslide removes the soil profile, there may not be adequate time for seedlings to grow and enhance soil stability. In this study, we investigate whether bistable dynamics emerge from the interaction of forest vegetation with the formation and accumulation of colluvial deposits in soil-mantled landscapes. To that end, we develop deterministic and stochastic models of landslide occurrence with a dynamic vegetation component. Results show that bistability exists for the deterministic case for both steep and shallow hollows under event and supply limited conditions. However, for the stochastic case, the randomness of landslide occurrence largely changed the states of the system such that the system only exhibited one stable state, which was the fully vegetated condition. Examining different management practices under stochastic conditions showed that the system eventually recovered; however, management practices influenced the recovery time of the forest. Thus, different management practices could render the land in a state of low vegetation over economically significant time periods.

Description: Seepage flux from ephemeral streams can be an important component of the water balance in arid and semi-arid regions. An emerging technique for quantifying this flux involves the measurement and simulation of a flood wave as it moves along an initially dry channel. This study investigates the usefulness of including surface water and groundwater data to improve model calibration when using this technique. We trialed this approach using a controlled flow event along a 1387 m reach of artificial stream channel. Observations were then simulated using a numerical model that combines the diffusion wave approximation of the Saint-Vénant equations for streamflow routing, with Philips’ infiltration equation and the groundwater flow equation. Model estimates of seepage flux for the upstream segments of the study reach, where streambed hydraulic conductivities were approximately 10 1 m d -1 , were on the order of 10 -4 m 3 d -1 m -2 . In the downstream segments, streambed hydraulic conductivities were generally much lower but highly variable (~10 -3 – 10 -7 m d -1 ). The Latin Hypercube Monte Carlo sensitivity analysis showed that the flood front timing, surface water stage, groundwater heads and the predicted stream flow seepage were most influenced by specific yield. Furthermore, inclusion of groundwater data resulted in a higher estimate of total seepage estimates than if the flood front timing were used alone.

Description: This paper evaluates the IBIS land surface model using daily soil moisture data over a three-year period (2005–2007) at a semi-arid site in southeastern Australia, the Stanley catchment, using the Monte-Carlo GLUE approach. The model was satisfactorily calibrated for both the surface 30 cm and full profile 90 cm. However, full-profile calibration was not as good as that for the surface, which results from some deficiencies in the evapotranspiration component in IBIS. Relatively small differences in simulated soil moisture were associated with large discrepancies in the predictions of surface runoff, drainage and evapotranspiration. We conclude that while land surface schemes may be effective at simulating heat fluxes they may be ineffective for prediction of hydrology unless the soil moisture is accurately estimated. Sensitivity analyses indicated that the soil moisture simulations were most sensitive to soil parameters, and the wilting point was the most identifiable parameter. Significant interactions existed between three soils parameters: porosity, saturated hydraulic conductivity and Campbell “b” exponent so they could not be identified independent of each other. There were no significant differences in parameter sensitivity and interaction for different hydroclimatic years. Even though the data record contained a very dry year and another year with a very large rainfall event, this indicated that the soil model could be calibrated without the data needing to explore the extreme range of dry and wet conditions. IBIS was much less sensitive to vegetation parameters. The leaf area index (LAI) could affect the mean of daily soil moisture time series when LAI 〈 1, while the variance of the soil moisture time series was sensitive to LAI 〉 1. IBIS was insensitive to the Jackson rooting parameter, suggesting that the effect of the rooting depth distribution on predictions of hydrology was insignificant. This article is protected by copyright. All rights reserved.

Description: We present a probabilistic sediment cascade model to simulate sediment transfer in a mountain basin (Illgraben, Switzerland) where sediment is produced by hillslope landslides and rockfalls and exported out of the basin by debris flows and floods. The model conceptualizes the fluvial system as a spatially lumped cascade of connected reservoirs representing hillslope and channel storages where sediment goes through cycles of storage and remobilization by surface runoff. The model includes all relevant hydrological processes that lead to runoff formation in an Alpine basin, such as precipitation, snow accumulation, snow melt, evapotranspiration, and soil water storage. Although the processes of sediment transfer and debris flow generation are described in a simplified manner, the model produces complex sediment discharge behavior which is driven by the availability of sediment and antecedent wetness conditions (system memory) as well as the triggering potential (climatic forcing). The observed probability distribution of debris flow volumes and their seasonality in 2000-2009 are reproduced. The stochasticity of hillslope sediment input is important for reproducing realistic sediment storage variability, although many details of the hillslope landslide triggering procedures are filtered out by the sediment transfer system. The model allows us to explicitly quantify the division into transport and supply-limited sediment discharge events. We show that debris flows may be generated for a wide range of rainfall intensities because of variable antecedent basin wetness and snowmelt contribution to runoff, which helps to understand the limitations of methods based on a single rainfall threshold for debris flow initiation in Alpine basins.

Description: This study presents a thermodynamic analysis to predict the type of initial, amorphous oxide overgrowth (i.e., am - Al 2 O 3 or am - SiO 2 ) on bare Al – Si alloy substrates. This analysis have taken into account the energies associated with both its interfaces (interface between the Al – Si alloy substrate and the thin oxide film and interface between the thin oxide film and vacuum) along with the bulk Gibbs free energy of oxide formation. This developed analysis is then applied for various parameters, such as, Si alloying element content at the substrate/oxide interface, the growth temperature, the oxide film thickness (up to 1 nm), and various low-index crystallographic surfaces of the substrate. It is found that am - SiO 2 overgrowth is thermodynamically preferred for a combination of lower oxide film thickness, lower growth temperature, and lower Si alloying content at the alloy/oxide interface. This is because of the overcompensation of the lower energies of both the interfaces over the bulk Gibbs free energy. Furthermore, it is found that for all cases, am - Al 2 O 3 forms a more stable interface with Al – Si alloy than am - SiO 2.

Description: In this study, translucent CeF 3 ceramics was firstly fabricated by hot-pressing, aiming to develop an economical, radiation hard scintillation material. After sintering at 950°C and pressure exertion of 300MPa for 2 h, a transmittance of 10–25% was achieved in the visible region, through a thickness of 0.75 mm. Microstructure investigations indicated the average grain size was around 42.7 μm, with no obvious pores observed on the fracture surface. The X-ray luminescence spectra, decay time, and thermoluminescence spectra were measured, and the results of the ceramics were similar to its corresponding single crystals, which proved it to be a promising scintillation material.

Description: Characterizing the complex geometries and the heterogeneity of the deposits in meandering river systems is a long-standing issue for the 3D modeling of alluvial formations. Such deposits are important sources of accessible groundwater in alluvial aquifers throughout the world and also play a major role as hydrocarbons reservoirs. In this paper we present a method to generate meandering river centerlines that are stochastic, geologically realistic, connected and conditioned to local observations or global geomorphological characteristics. The method is based on fast 1D multiple-point statistics in a transformed curvilinear domain: the succession in directions observed in a real world meandering river (the analog) is considered as statistical model for multiple-point statistics simulation. The integration of local data is accomplished by an inverse procedure ensuring that the channels pass through a given set of locations while conserving the high-order spatial characteristics of an analog. The methodology is applied on seven real world case studies. This work demonstrates the flexibility and the applicability of multiple-point statistics outside the standard paradigm that considers the simulation of a 2D or 3D variable with spatial coordinates.

Description: Large asymmetric bedforms commonly develop in rivers. The turbulence associated with flow separation that develops over their steep lee side is responsible for the form shear stress which can represent a substantial part of total shear stress in rivers. This paper uses the Delft3D modeling system to investigate the effects of bedform geometry and forcing conditions on flow separation length and associated turbulence, and bedform shear stress over angle-of-repose (30° lee side angle) bedforms. The model was validated with laboratory measurements that showed sufficient agreement to be used for a systematic analysis. The influence of flow velocity, bed roughness, relative height (bedform height / water depth) and aspect ratio (bedform height / length) on the variations of the normalized length of the flow separation zone, the extent of the wake region (where the turbulent kinetic energy (TKE) was more than 70% of the maximum TKE), the average TKE within the wake region and the form shear stress were investigated. Form shear stress was found not to scale with the size of the flow separation zone but to be related to the product of the normalized extent of the wake region (extent of the wake region / extent of water body above the bedform) and the average TKE within the wake region. The results add to understanding of the hydrodynamics of bedforms and may be used for the development of better parameterizations of small-scale processes for application in large-scale studies.

Description: We explore the bankfull width ( W bf ) vs. drainage area ( A da ) relationship across a range of climatic and geologic environments, and ask (1) is the relationship between ln ( W bf ) and ln ( A da ) best described by a linear function and (2) can a reliable relationship be developed for predicting W bf with A da as the only independent variable. The principal dataset for this study was compiled from regional curve studies and other reports that represent 1,018 sites (1 m ≤ W bf ≤ 110 m and 0.50 km 2 ≤ A da ≤ 22,000 km 2 ) in the continental U.S. Two additional datasets were used for validation. After dividing the data into small-, medium-, and large-size basins which, respectfully, correspond to A da 〈 4.95 km 2 , 4.95 km 2 ≤ A da 〈 337 km 2 , and A da ≥ 337 km 2 , regression lines from each dataset were compared using one-way analysis of covariance (ANCOVA). A second ANCOVA was performed to determine if mean annual precipitation ( P ) is an extraneous factor in the W bf vs. A da relationship. The ANCOVA results reveal that using A da alone does not yield a reliable W bf vs. A da relationship that is applicable across a wide range of environments and that P is a significant extraneous factor in the relationship. Considering data for very small basins ( A da ≤ 0.49 km 2 ) and very large basins ( A da ≥ 1.0×10 5 km 2 ) we conclude that a two-segment linear model is the most probable form of the ln ( W bf ) vs. ln ( A da ) relationship. This study provides useful information for building complex multivariate models for predicting W bf .

Description: The effect of direct integration of PZT into LTCC by sintering on its polarization behavior is studied. Due to the typically small thermal expansion of PZT, compressive stresses are expected after cooling, which should promote poling. However, measurements show a strongly reduced remanent polarization compared to free PZT. Analysis of thermal expansion reveals that thermal shrinkage of the PZT is larger than that of LTCC. The impact of this finding on the poling behavior is simulated for the LTCC/PZT module, confirming the experimental observations and allowing for conclusions concerning the suitability of the PZT for integration into LTCC structures.

Description: Over the past few decades ground water has become an essential commodity due to increased demand as a result of growing population, industrialization, urbanization etc. The water supply situation is expected to become more severe in the future because of continued unsustainable water use and projected change in hydro-meteorological parameters due to climate change. This study is based on the integrated approach of Remote Sensing (RS), Geographical Information System (GIS) and Multi-Criteria Decision Making (MCDM) techniques to determine the most important contributing factors that affect the ground water resources and to delineate the ground water potential zones. Ten thematic layers viz. geomorphology, geology, soil, topographic elevation (DEM), land use/land cover, drainage density, lineament density, proximity of surface water bodies, surface temperature and post-monsoon ground water depth were considered for the present study. These thematic layers were selected for ground water prospecting based on literature, discussion with the experts of Central Ground Water Board (CGWB), Government of India, field observations, geophysical investigation and multivariate techniques. The thematic layers and their features were assigned suitable weights on the Saaty's scale according to their relative significance for ground water occurrence. The assigned weights of the layers and their features were normalized by using Analytic Hierarchy Process (AHP) and eigenvector method. Finally, the selected thematic maps were integrated using weighted linear combination method to create the final ground water potential zone map. The final output map shows different zones of ground water potential, viz., very good (16%), good (35%), moderate (28%) low (17%) and very low (2.1%). The ground water potential zone map was finally validated using the discharge and ground water depth data from 28 and 98 pumping wells respectively which showed good correlation. This article is protected by copyright. All rights reserved.

Description: Bank erosion is the main source of suspended sediment (SS) and diffuse total phosphorus (TP) in many lowland catchments. This study compared a physically based sediment routing method (Physical method), which distinguishes between stream bed and bank erosion, with the original sediment routing method (Original method) within the Soil and Water Assessment Tool (SWAT) version 2009, for simulating SS and TP losses from a lowland catchment. A SWAT model was set up for the lowland River Odense catchment in Denmark and calibrated against observed stream flow and phosphate (PO 4 ) loads. Based on an initial calibration of hydrological and PO 4 parameters, the SWAT model with the Original method (Original model) and the SWAT model with the Physical method (Physical model) were calibrated separately against observed SS and TP loads. The SWAT model simulated daily stream flow well, but underestimated PO 4 loads. The Physical model simulated daily SS and TP better than the Original model. The simulated contribution of bank erosion to SS in the Physical model (99%) was close to the estimated contribution from in-situ erosion measurements (90-94%). Compared with the Original method, the Physical method is not only more conceptually correct, but also improves model performance. This article is protected by copyright. All rights reserved.

Description: This paper focuses on surface-subsurface water exchange in a steep coarse-bedded stream with step-pool morphology. We use both flume experiments and numerical modelling to investigate the influence of stream discharge, channel slope and sediment hydraulic conductivity on hyporheic exchange. The model step-pool reach, whose topography is scaled from a natural river, consists of 3 step-pool units with 0.1 m step heights, discharges ranging between base and over bankfull flows (scaled values of 0.3-4.5 L/s) and slopes of 4 and 8%. Results indicate that the deepest hyporheic flow occurs with the steeper slope and at moderate discharges and downwelling fluxes at the base of steps are highest at the largest stream discharges. In contrast to pool-riffle morphology, these findings show that steep slopes cause deeper surface-subsurface exchanges than gentle slopes. Numerical simulation results show that the portion of the hyporheic zone influenced by surface water temperature increases with sediment hydraulic conductivity. These experiments and numerical simulations emphasize the importance of topography, sediment permeability and roughness elements along the channel surface in governing the locations and magnitude of downwelling fluxes and hyporheic exchange. Our results show that hyporheic zones in these steep streams are thicker than previously expected by extending the results from streams with pool-riffle bed forms. This article is protected by copyright. All rights reserved.

Description: Recent studies show that multimodel combinations improve hydroclimatic predictions by reducing model uncertainty. Given that climate forecasts are available from multiple climate models, which could be ingested with multiple watershed models, what is the best strategy to reduce the uncertainty in streamflow forecasts? To address this question, we consider three possible strategies: (1) reduce the input uncertainty first by combining climate models and then use the multimodel climate forecasts with multiple watershed models (MM-P) (2) ingest the individual climate forecasts (without multimodel combination) with various watershed models and then combine the streamflow predictions that arise from all possible combinations of climate and watershed models (MM-Q)(3) combine the streamflow forecasts obtained from multiple watershed models based on strategy (1) to develop a single streamflow prediction that reduces uncertainty in both climate forecasts and watershed models (MM-PQ) . For this purpose, we consider synthetic schemes that generate streamflow and climate forecasts, for comparing the performance of three strategies with the true streamflow generated by a given hydrologic model. Results from the synthetic study show that reducing input uncertainty first ( MM-P ) by combining climate forecasts results in reduced error in predicting the true streamflow compared to the error of multimodel streamflow forecasts obtained by combining streamflow forecasts from all-possible combination of individual climate model with various hydrologic models ( MM-Q ). Since the true hydrologic model structure is unknown, it is desirable to consider MM-PQ as an alternate choice that reduces both input uncertainty and hydrologic model uncertainty. Application on two watersheds in NC also indicates that reducing the input uncertainty first is critical before reducing the hydrologic model uncertainty.

Description: The understanding of reasons leading to non-uniqueness of soil erosion susceptibility is still inadequate, yet indispensable for establishing general relations between runoff volume and sediment yield. To obtain relevant insights, we performed a series of numerical simulations with a detailed hydrodynamic model using synthetic storms of varying intensity, duration, and lag time between events as representations of different hydrologic response conditions in a zero-order catchment. The design targeted to generate surface flow and ‘perturb’ soil substrate by a first rainfall event, creating a set of initial conditions in terms of flow and deposited sediment prior to the onset of a subsequent rainfall event. Due to the differential effect of (re)detachment and (re)entrainment processes on soil particles of varying sizes, the deposited sediment mass formed shielding layer. One of the essential results is that unless the initial condition of flow and sediment is identical, the same volume of runoff can generate different total sediment yields and their variation can reach up to ~200%. The effect is attributed to two major conflicting effects exerted by the deposited ‘initialization’ (soil antecedent condition) sediment mass: erosion enhancement, because of supply of highly erodible sediment, and erosion impediment, because of constrain on the availability of lighter particles by heavier sediment. Consistently with this inference, long-term simulations with continuous rainfall show that a peculiar feature of sediment yield series is the existence of maximum before the steady-state is reached. The two characteristic time scales, the time to peak and the time to steady-state, separate three characteristic periods that correspond to flow-limited, source-limited, and steady-state regimes. These time scales are log-linearly and negatively related to the spatially averaged Shields parameter: the smaller the rainfall input and the heavier a given particle is, the larger the two scales are. The results provide insights on how the existence of shield operates on erosion processes, possibly implying that accurate short-term predictions of geomorphic events from headwater areas may never become a tractable problem: the latter would require a detailed spatial characterization of particle size distribution prior to precipitation events.

Description: Aquifer hydraulic properties such as hydraulic conductivity ( K ) are ubiquitously heterogeneous and typically only a statistical characterization can be sought. Additionally statistical anisotropy at typical characterization scales is the rule. Thus, regardless of the processes governing solute transport at the local (pore) scale, transport becomes non-Fickian. Mass-transfer models provide an efficient tool that reproduces observed anomalous transport; in some cases though, these models lack predictability as model parameters cannot readily be connected to the physical properties of aquifers. In this study we focus on a multi-rate mass-transfer model (MRMT), and in particular the apparent capacity coefficient (β), which is a strong indicator of the potential of immobile zones to capture moving solute. We aim to find if the choice of an apparent β can be phenomenologically related to measures of statistical anisotropy. We analyzed an ensemble of random simulations of three-dimensional log-transformed multi-Gaussian permeability fields with stationary anisotropic correlation under convergent flow conditions. It was found that apparent β also displays an anisotropic behavior, physically controlled by the aquifer directional connectivity, which in turn is controlled by the anisotropic correlation model. A high hydraulic connectivity results in large β values. These results provide new insights into the practical use of mass-transfer models for predictive purposes.

Description: Sintered nanoceramics of Pr -doped lanthanum hafnate, La 2 Hf 2 O 7 : Pr , were prepared by means of a high-pressure sintering technique using nanopowders made by Pechini method. Structure, morphology, and spectroscopic properties of the ceramics compared to the starting powder are presented and discussed. Emission and excitation spectra recorded at room temperature as well as at 7 K using synchrotron radiation are presented together with results of luminescence kinetics measurements. In ceramics, at 7 K, the Pr 3+ luminescence from 3 P 0 (blue-green, green, and red region) and 1 D 2 (red) levels is accompanied by a broad-band emission located in the 380–530 nm range of wavelengths, whereas powders gives only the Pr 3+ -related luminescence. Depending on the excitation wavelength, the broad-band emission maximum moves between 430 and 470 nm indicating superposition of at least two components. In sintered nanoceramics, the lifetimes of Pr 3+ emissions from 3 P 0 and 1 D 2 levels were by 10%–20% shorter compared to the powder. The existence of different luminescence centers was proved by the selective emission decays examination. The fast 5 d → 4 f luminescence of Pr 3+ was not observed from either of the two types of La 2 Hf 2 O 7 :Pr materials.

Description: In this paper, an operational algorithm is proposed for the mapping of surface moisture over the northern and central parts of Tunisia, in North Africa. A change detection approach is applied, using 160 multi-incidence Envisat ASAR Wide Swath images acquired in the horizontal polarization over a 7-year period. Parameterization of this algorithm is considered for three classes of vegetation cover density (NDVI〈0.25, 0.25〈NDVI〈0.5 and NDVI〉0.5), retrieved from SPOT-VGT decadal images. A relative soil moisture index, ranging between 0 (for the driest surfaces) and 1 (for saturated soils), is proposed for each date, with a resolution of 1 km. The retrieved soil moistures are validated by means of ground measurements based on continuous thetaprobe measurements, as well as low resolution (25 km) ERS and ASCAT soil moisture products from the Vienna University of Technology (TU Wien). A qualitative relationship between spatio-temporal variations of moisture and precipitation is also discussed.

Description: The objective of this analysis is to demonstrate the feasibility of using a composite L2 SMOS soil moisture product for determining drought conditions by taking advantage of its spatial and temporal resolution. The work investigates the potential relationships between soil moisture anomalies and two drought indices, the Standardized Precipitation Index (SPI) and the Standardized Precipitation Evapotranspiration Index (SPEI), both calculated on a ten-day basis. As the two drought indices can be applied to different time scales for precipitation series, the influence of time scale on the drought definition is also studied. The anomalies were calculated both for the in situ soil moisture by REMEDHUS (Soil Moisture Measurement Stations Network, Spain) and from the SMOS L2 soil moisture product. In general, in situ anomalies exhibit higher correlation coefficients for the drought indices than those of SMOS, except for the shortest time scale. As expected, the short-term remotely sensed anomalies have a high response to precipitation events. This effect may be due to the greater sensitivity of SMOS data to rainfall, as well as to the spatial averaged nature of its observations. The optimal time scale was one month for the SMOS values and ranged between 30 and 50 days for the in situ values. The use of evapotranspiration in the calculation of the indices did not improve the description of the anomalies. The relationship between indices and soil moisture conditions provides encouraging results. Indeed, this method generates preliminary but valuable insights for future satellite products. This article is protected by copyright. All rights reserved.

Description: TiO 2 -doped NiFe 2 O 4 samples were prepared via ball-milling and two-step sintering processes. Besides NiFe 2 O 4 phase, two new phases, NiTiO 3 and Fe 2 TiO 5 , formed in TiO 2 -doped samples. The temperature of sintering onset for 1.0 wt% TiO 2 -doped samples is 230°C lower than that of undoped samples. Early-stage synthesis process of TiO 2 -doped NiFe 2 O 4 ceramics is controlled by grain boundary diffusion mechanism. Increasing TiO 2 content from 0 to 1.0 wt%, the apparent activation energy decreased from 813.919 KJ/mol to 639.361 KJ/mol. The values of relative density and bending strength reached their maximum value with 1.0 wt% TiO 2 . Saturation magnetization, residual magnetization ratio and coercivity decrease with increasing TiO 2 content.

Description: This article demonstrates the energy harvesting capabilities of bulk lead-free ferroelectric materials using Olsen cycle. Lead-free compositions K[(Nb 0.90 Ta 0.10 ) 0.99 Mn 0.01 ]O 3 (KNTM) and 0.5Ba(Zr 0.2 Ti 0 . 8 )O 3 -0.5(Ba 0.7 Ca 0.3 )TiO 3 (BZT-50BCT) have been explored for their energy harvesting capabilities. Further, we studied the variation in energy density over wide range of applied temperature and electric field. The maximum harness-able energy density for KNTM and BZT-50BCT compositions are found to be 629 J/L (629 kJ/m 3 ) and 87 J/L (87 kJ/m 3 ), respectively. This estimated energy density, obtained for the bulk samples, is comparable with the highest energy density reported to-date (888 J/L for lanthanum-doped lead zirconate titanate (8/65/35 PLZT) thick films).

Description: Remote estimation of river discharge from river width variations is an intriguing method for gauging rivers without conventional measurements. Entirely cloud-free imagery of an entire river reach is often rare, but partial coverage is more frequent. Discharge is estimated from spatially discontinuous imagery via construction of multiple width-discharge rating curves within a 62 km reach of the Tanana River, Alaska. The resulting discharge error is as low as 6.7% RMSE. Imagery covering 〈20% of the study reach can be used. This article is protected by copyright. All rights reserved.

Description: Zirconium diboride–based ceramics were studied to explore the effects of high secondary phase content and the relative merits of silicon carbide and zirconium carbide. Densification and grain growth were not influenced by heating rate or composition. Hardness increased with silicon carbide or zirconium carbide content from 18 to 26GPa and 18 to 21.5GPa respectively. Flexural strength was controlled by grain size and mode of crack propagation, and maximized at 380MPa at 20wt% silicon carbide and 590 MPa at 50wt% zirconium carbide. Zirconium carbide enhanced flexural strength better, due to its smaller coefficient of thermal expansion mismatch with zirconium diboride.

Description: Using the association of low-cost screen-printing technology with the sacrificial layer method, the feasibility of totally released piezoelectric thick-film microceramics of gold-electroded PZT type is studied. After the deposition of the sacrificial layer on an alumina substrate and subsequent printing and drying of gold/PZT/gold layers, the final firing is performed at low temperature. This is followed by the releasing step of the Au/PZT/Au in a diluted acidic solution. Impedance analysis shows that the electrical properties and electromechanical coefficients of poled PZT thick films are still lower than those of PZT ceramics. This result is correlated with the high porosity rate of the PZT layer. However, these piezoelectric microceramics present desirable electromechanical behavior and can therefore be used as dynamic-mode sensors or as actuators to generate vibrations in a structure on which they are bonded. Moreover, the successful fabrication and the electromechanical impedance (EMI) levels obtained on a metallic structure suggest possible structural health monitoring applications of these screen-printed PZT microceramics.

Description: We have investigated the electromechanical response of potassium sodium niobate ( K 0.5 Na 0.5 NbO 3 or KNN) thick films. The high-field strain hysteresis loops and weak-field converse piezoelectric d 33 coefficient of the films were measured and compared with those of KNN bulk ceramics under the same electric field conditions. The converse d 33 values of the thick films and bulk ceramics were equal to 82.5 and 138 pm/V, respectively, at 0.4 kV/mm. The fundamental difference between the piezoelectric response of the KNN films and the ceramics was studied in terms of the effective (“clamped”) piezoelectric d 33 coefficient. The reduction in the piezoelectric d 33 coefficient of the KNN films, resulting from the clamping by the substrate, was compared to lead-based ferroelectric thick films, including Pb ( Zr , Ti ) O 3 (PZT) and (1 −  x ) Pb ( Mg 1/3 Nb 2/3 ) O 3 − x PbTiO 3 (PMN-PT). We propose a possible explanation, based on the particular elastic properties of KNN, for the small relative difference observed between the “clamped” and “unclamped” (“bulk”) d 33 of KNN, in comparison with lead-based systems.

Description: Lead zirconate titanate (PbZr 1 −  x Ti x O 3 , PZT)/epoxy composites with one- dimensional epoxy in PZT matrix (called 3-1 type piezocomposites) have been fabricated by tert-butyl alcohol (TBA)-based directional freeze casting of PZT matrix and afterward infiltration of epoxy. The composites with PZT volume fraction ranging from 0.36 to 0.69 were obtained by adjusting initial solid loading in freeze-casting slurry. The effect of poling voltage on piezoelectric properties of the composites was studied for various volume fraction of PZT phase. With the increasing of PZT volume fraction, relative permittivity (ε r ) increased linearly and piezoelectric coefficient ( d 33 and d 31 ) increased step by step. The resultant composites with 0.57 PZT volume fraction possessed the highest hydrostatic piezoelectric strain coefficient ( d h ) value (184 pC/N), voltage coefficient ( g h ) value (13.6 × 10 −3  V/m Pa), and hydrostatic figure of merit (HFOM) value (2168 × 10 −15  Pa −1 ).

Description: Multiphase borosilicate glass-ceramics represent one candidate to contain radioactive nuclear waste separated from used nuclear fuel. In this work, the thermophysical properties from room temperature to 1273 K were investigated for four different borosilicate glass-ceramic compositions containing waste loadings from 42 to 60 wt% to determine the sensitivity of these properties to waste loading, as-fabricated microstructure, and potential evolutions in microstructure brought about by temperature transients. The thermal expansion, specific heat capacity, thermal diffusivity, and thermal conductivity are presented. The impact of increasing waste loading is shown to have a small but measurable effect on the thermophysical properties between the four compositions, contrasted to a much greater impact observed when transitioning from predominantly crystalline to amorphous systems. Thermal cycling below 1273 K was not found to measurably impact the thermophysical properties of the compositions investigated here.

Description: This paper presents the results of a data based comparative study of several hundred catchments across continental United States belonging to the MOPEX dataset to systematically explore the connection between the flood frequency curve and mean annual water balance. Mean annual water balance is expressed in terms of two similarity measures: (i) the climatic aridity index, AI , which is a measure of the competition between energy and water availability; and (ii) the baseflow index, BFI , which is a measure of total runoff partitioning into surface and subsurface components at the annual time scale. The data analyses showed that the aridity index, AI , has a first order control on the shape of the flood frequency curve (also known as the growth curve), as expressed in terms of both the mean and coefficient of variation ( C v ) of the annual maximum floods, once normalized by catchment size (i.e., specific flood discharge) While the mean annual (specific) flood discharge decreases with increasing aridity, C v increases with increasing aridity. On the other hand, the BFI was found to be a second order control on the flood frequency curve. Higher BFI , meaning higher contributions of subsurface flow to total streamflow, leads to a decrease of the mean annual (specific) flood discharge, and vice versa. The statistically significant relationship between AI and the flood frequency curve and the consistent shift of the growth curves with AI support the use of AI as a similarity measure for regionalization of flood frequency.

Description: The performance of glacio-hydrological models which simulate catchment response to climate variability depends to a large degree on the data used to force the models. The forcing data become increasingly important in high elevation, glacierised catchments where the interplay between extreme topography, climate and the cryosphere is complex. It is challenging to generate a reliable forcing dataset that captures this spatial heterogeneity. In this paper, we analyze the results of a one year field campaign focusing on air temperature and precipitation observations in the Langtang Valley in the Nepalese Himalayas. We use the observed time series to characterize both temperature lapse rates (LRs) and precipitation gradients (PGs). We study their spatial and temporal variability, and we attempt to identify possible controlling factors. We show that very clear LRs exist in the valley and that there are strong seasonal differences related to the water vapor content in the atmosphere. Results also show that the LRs are generally shallower than the commonly used environmental lapse rates. The analysis of the precipitation observations reveal that there is great variability in precipitation over short horizontal distances. A uniform valley wide PG cannot be established and several scale-dependent mechanisms may explain our observations. We complete our analysis by showing the impact of the observed LRs and PGs on the outputs of the TOPKAPI-ETH glacio-hydrological model. We conclude that LRs and PGs have a very large impact on the water balance composition and that short-term monitoring campaigns have the potential to improve model quality considerably.

Description: Many physically-based hydrological/hydrogeological models used for predicting groundwater seepage areas, including topography-based index models such as TOPMODEL, rely on the Dupuit assumption. To ensure the sound use of these simplified models, knowledge of the conditions under which they provide a reasonable approximation is critical. In this study, a Dupuit solution for the seepage length in hillslope cross-sections is tested against a full-depth solution of saturated groundwater flow. In homogeneous hillslopes with horizontal impervious base and constant-slope topography, the comparison reveals that the validity of the Dupuit solution depends not only on the ratio of depth to hillslope length d/L (as might be expected), but also on the ratio of hydraulic conductivity to recharge K/R and on the topographic slope s . The validity of the Dupuit solution is shown to be in fact a unique function of another ratio, the ratio of depth to seepage length d/L S . For d/L S 〈 0.2, the relative difference between the two solutions is quite small (〈 14% for the wide range of parameter values tested), whereas for d/L S 〉 0.2, it increases dramatically. In practice, this criterion can be used to test the validity of Dupuit solutions. When d/L S increases beyond that cut-off, the ratio of seepage length to hillslope length L S /L given by the full-depth solution tends towards a non-zero asymptotic value. This asymptotic value is shown to be controlled by (and in many cases equal to) the parameter R/sK . Generalization of the findings to cases featuring heterogeneity, non-horizontal impervious base and variable-slope topography is discussed.

Description: Isothermal compositional flow models require coupling transient compressible flows and advective transport systems of various chemical species in subsurface porous media. Building such numerical models is quite challenging and may be subject to many sources of uncertainties because of possible incomplete representation of some geological parameters that characterize the system's processes. Advanced data assimilation methods, such as the ensemble Kalman filter (EnKF), can be used to calibrate these models by incorporating available data. In this work, we consider the problem of estimating reservoir permeability using information about phase pressure as well as the chemical properties of fluid components. We carry out state-parameter estimation experiments using joint and dual updating schemes in the context of the EnKF with a two-dimensional single-phase compositional flow model (CFM). Quantitative and statistical analyses are performed to evaluate and compare the performance of the assimilation schemes. Our results indicate that including chemical composition data significantly enhances the accuracy of the permeability estimates. In addition, composition data provide more information to estimate system states and parameters than do standard pressure data. The dual state-parameter estimation scheme provides about 10% more accurate permeability estimates on average than the joint scheme when implemented with the same ensemble members, at the cost of twice more forward model integrations. At similar computational cost, the dual approach becomes only beneficial after using large enough ensembles.

Description: Groundwater-fed irrigation has been shown to deplete groundwater storage, decrease surface water runoff and increase evapotranspiration. Here we simulate soil moisture dependent groundwater-fed irrigation with an integrated hydrologic model. This allows for direct consideration of feedbacks between irrigation demand and groundwater depth. Special attention is paid to system dynamics in order to characterized spatial variability in irrigation demand and response to increased irrigation stress. A total of 80 years of simulation are completed for the Little Washita Basin in Southwestern Oklahoma, USA spanning a range of agricultural development scenarios and management practices. Results show regionally aggregated irrigation impacts consistent with other studies. However, here a spectral analysis reveals that groundwater-fed irrigation is also shown to amplify the annual streamflow cycle while dampening longer-term cyclical behavior with increased irrigation during climatological dry periods. Feedbacks between the managed and natural system are clearly observed with respect to both irrigation demand and utilization when water table depths are within a critical range. Although the model domain is heterogeneous with respect to both surface and subsurface parameters, relationships between irrigation demand, water table depth and irrigation utilization are consistent across space and between scenarios. Still, significant local heterogeneities are observed both with respect to transient behavior and response to stress. Spatial analysis of transient behavior shows that farms with groundwater depths within a critical depth range are most sensitive to management changes. Differences in behavior highlight the importance of groundwater's role in system dynamics in addition to water availability.

Description: Including positive feedbacks in hydrological models has recently been shown to result in complex behavior with multiple steady states. When a large disturbance, say a major drought, is simulated within such models the hydrology changes. Once the disturbance ends the hydrology does not return to that prior to the disturbance, but rather, persists within an alternate state. These multiple steady states (henceforth attractors ) exist for a single model parameterization and cause the system to have a finite resilience to such transient disturbances. A limitation of past hydrological resilience studies is that multiple attractors have been identified using mean annual or mean monthly forcing. Considering that most hydrological fluxes are subject to significant forcing stochasticity and do not operate at such large time scales, it remains an open question whether multiple hydrological attractors can exist when a catchment is subject to stochastic daily forcing. This question is the focus of this paper and it needs to be addressed prior to searching for multiple hydrological attractors in the field. To investigate this, a previously developed semi-distributed hill-slope ecohydrological model was adopted which exhibited multiple steady states under average monthly climate forcing. In this paper, the ecohydrological model was used to explore if feedbacks between the vegetation and a saline water table result in two attractors existing under daily stochastic forcing. The attractors and the threshold between them (henceforth repellor ) were quantified using a new limit cycle continuation technique that up-scaled climate forcing from daily to monthly (model and limit cycle code is freely available). The method was used to determine the values of saturated lateral hydraulic conductivity at which multiple attractors exist. These estimates were then assessed against time-integration estimates, which they agreed with. Overall, multiple attractors where found to exist under stochastic daily forcing. However, changing the climate forcing from monthly to daily did significantly reduce the parameter range over which two attractors existed. This suggests fewer catchments may have multiple attractors than previously considered.

Description: The companion paper showed that multiple steady state groundwater levels can exist within a hill-slope Boussinesq-vegetation model under daily stochastic forcing. Using a numerical limit-cycle continuation algorithm, the steady states (henceforth attractors ) and the threshold between them (henceforth repellor ) were quantified at a range of saturated lateral conductivity values, . This paper investigates if stochastic daily forcing can switch the catchment between both of the attractors. That is, an attractor may exist under average forcing conditions but can stochastic forcing switch the catchment into and out of each of the attractor basins?; i.e. making the attractor emerge . This was undertaken using the model of the companion paper and by completing daily time-integration simulations at six values of the saturated lateral hydraulic conductivity, ; three having two attractors and three having only a deep water table attractor. By graphically analyzing the simulations, and comparing against simulations from a model modified to have only one attractor, multiple attractors were found to emerge under stochastic daily forcing. However, the emergence of attractors was significantly more subtle and complex than that suggested by the companion paper. That is, an attractor may exist but never emerge; both attractors may exist and both may emerge but identifying the switching between attractors was often ambiguous; and only one attractor may exist and but a second temporary attractor may exist and emerge during periods of high precipitation. This subtle and complex emergence of attractors was explained using continuation analysis of the climate forcing rate, and not a model parameter such as . It showed that the temporary attractor existed over a large range of values and this suggests that more catchments may have multiple attractors than suggested by the companion paper. By combining this continuation analysis with the time-integration simulations, hydrological signatures indicative of a switch of multiple attractors were proposed. These signatures may provide a means for identifying actual catchments that have switched between multiple attractors.

Description: ABSTRACT A large-time analytical solution is proposed for the spatial variance and coefficient of variation of the depth-averaged concentration due to instantaneous, cross-sectionally uniform solute sources in pseudo-rectangular open channel flows. The mathematical approach is based on the use of the Green functions and on the Fourier decomposition of the depth-averaged velocities, coupled with the method of the images. The variance spatial trend is characterized by a minimum at the center of the mass and two mobile, decaying symmetrical peaks which, at very large times, are located at the inflexion points of the average Gaussian distribution. The coefficient of variation, which provides an estimate of the expected percentage deviation of the depth-averaged point concentrations about the section-average, exhibits a minimum at the center which decays like t -1 and only depends on the river diffusive time-scale. The defect of cross-sectional mixing quickly increases with the distance from the center, and almost linearly at large times. Accurate numerical Lagrangian simulations were performed to validate the analytical results in pre-asymptotic and asymptotic conditions, referring to a particularly representative sample case for which cross-sectional depth and velocity measurements were known from a field survey. In addition, in order to discuss the practical usefulness of computing large-time concentration spatial moments in river flows, and resorting to directly measured input data, the order of magnitude of section-averaged concentrations and corresponding coefficients of variation was estimated in field conditions and for hypothetical contamination scenarios, considering a unit normalized mass impulsively injected across the transverse section of 81 U.S. rivers.

Description: Reports of abruptly declining flows of Canada's Athabasca River have prompted concern because this large, free-flowing river could be representative for northern North America, provides water for the massive Athabasca oil-sands projects, and flows to the extensive and biodiverse Peace-Athabasca, Slave and Mackenzie River deltas. To investigate historic hydrology along the river and its major tributaries we expanded the time series with interpolations for short data gaps; calculations of annual discharges from early, summer-only records; and by splicing records across sequential hydrometric gauges. These produced composite, century-long records (1913 to 2011) and trend detection with linear Pearson correlation provided similar outcomes to non-parametric Kendall τ−b tests. These revealed that the mountain and foothills reaches displayed slight increases in winter discharges versus larger declines in summer discharges, and consequently declining annual flows (~0.16%/yr at Hinton; p 〈 0.01). Conversely, with contrasting boreal contributions, the Athabasca River at Athabasca displayed no overall trend in monthly or annual flows, but there was correspondence with the Pacific Decadal Oscillation (PDO) that contributed to a temporary flow decline from 1970 to 2000. These findings from century-long records contrast with interpretations from numerous shorter-term studies, and emphasize the need for sufficient time series for hydrologic trend analyses. For Northern Hemisphere rivers, the study interval should be at least 80 years to span two PDO cycles and dampen the influence from phase transitions. Most prior trend analyses considered only a few decades and this weakens interpretations of the hydrologic consequences of climate change. This article is protected by copyright. All rights reserved.

Description: The thermal cycling reliability of ceramic/metal laminates is critical for their applications in microelectronic modules. A eutectic bonding process is used to prepare copper/sapphire bilayers in this study. Due to elastic and thermal mismatches between alumina and copper, the Cu/sapphire bilayer cannot pass a thermal cycling test. The thermal cycling reliability can be improved through the use of a metallic nickel interlayer. During the bonding process, the nickel interlayer was oxidized first, reacted with alumina to form a NiAl 2 O 4 spinel phase. The thermal diffusivity of the bilayer with spinel interphase remains high after the temperature cycling test.

Description: Nano-sized (i) N-doped sodium trititanate and (ii) N and Cu 2+ (Ag + ) co-doped sodium trititanates CuTi 3 NO 6−x (Ag 2 Ti 3 NO 6−x ) were prepared by a solid-state and ion-exchange methods, respectively. The materials were characterized by EDS, PXRD, XPS, FESEM, TEM, UV–visible DRS, and Raman spectroscopy. All the materials were crystallized in monoclinic lattice with P21/m space group. The bandgap energy of all the samples was deduced from their UV–visible DRS profiles. Visible-light-induced photocatalytic oxidation of the methylene blue (MB) and methyl orange (MO), cyclohexene and phenol, was examined. The Ag + co-doped trititanate exhibited the highest photocatalytic activity among the materials investigated.

Description: Recent studies have suggested that the hydrologic connectivity of northern headwater catchments is likely controlled by antecedent moisture conditions and land cover patterns. A water storage model ( EWS ), based on water levels ( WLs ), specific yield ( Sy ) and surface elevation ( SE ) changes was compared to a basic water budget of a small, boreal, patterned fen (13 ha) during the ice-free period. Results showed that the EWS model reproduced well storage variations derived from the water budget. These results suggest that storage variations can be properly represented by the fluctuations of WLs when we consider the heterogeneous soil properties. However, storage deviations occurred at the daily scale and could be explained by a lack of information on water retention in unsaturated layers, canopy interceptions and preferential flows. Despite the significant impact of SE changes on the different peatland cover storage budgets (strings and lawns), using Sy mean values had low impact on storage estimations. This can be explained by the large proportion of pools and high WLs throughout the fen. At the fen scale, high storage in the pools seemed to reduce the Sy difference between strings and lawns. The results of this study provide new insights about the complex hydrological behaviour of northern catchments and allow for conceiving new hydrological modeling perspectives. This article is protected by copyright. All rights reserved.

Description: We report a novel oscillatory pressure-assisted hot-pressing process for preparing high-quality ceramics. Compared with the samples prepared by conventional pressureless sintering (PS) and hot-pressing (HP), the zirconia ceramic prepared by oscillatory pressure-assisted hot-pressing (OPAHP) exhibited a higher density, smaller grain size, and more homogeneous structure. More remarkably, the strength of the OPAHP sample reached 1556 MPa, which is much higher than the samples prepared by other two techniques. The results suggest that OPAHP is a more effective technique for preparing high-quality zirconia, which is likely applicable to other material systems.

Description: Two major causes of salt marsh loss are vertical drowning, when sediment accumulation on the platform cannot keep vertical pace with sea level rise, and horizontal retreat, associated with wave-induced marsh boundary erosion. Despite these processes having been extensively documented and modeled, is unclear which loss modality dominates given a set of environmental parameters. A three-point dynamic model was developed to predict marsh loss as a function of sea level rise, allochthonous sediment supply, wind regime, tidal range, and marsh bank and mudflat erodability. Marsh horizontal and vertical evolution was found to respond in opposing ways to wave induced erosion processes. Marsh horizontal retreat was triggered by large mudflats, strong winds, high erodability of marsh bank and mudflat, whereas the opposite conditions acted to reduce the sediment supply to the marsh platform, promoting marsh loss to drowning. With low and moderate rates of sea level rise (~ 5 mm/yr), retreat was found to be a more likely marsh loss modality than drowning. However, conditions associated with marsh retreat also increase the system resilience by transferring sediment on the marsh platform and preventing drowning. Our results suggest the use of a modular strategy for short-term marsh management: selectively protect extensive salt marsh regions by maintaining healthy vegetation on the platform, while allowing other areas to retreat, leveraging the natural resilience embedded in the lateral loss of marsh extent.

Description: Observational data and the Princeton Urban Canopy Model, with its detailed representation of urban heterogeneity and hydrological processes, are combined to study evaporation and turbulent water vapor transport over urban areas. The analyses focus on periods before and after precipitation events, at two sites in the Northeastern United States. Our results indicate that while evaporation from concrete pavements, building rooftops and asphalt surfaces is discontinuous and intermittent, overall these surfaces accounted for nearly 18% of total latent heat fluxes (LE) during a relatively wet 10-day period. More importantly, these evaporative fluxes have a significant impact on the urban surface energy balance, particularly during the 48 hours following a rain event when impervious evaporation is the highest. Thus, their accurate representation in urban models is critical. Impervious evaporation after rainfall is also shown to correlate the sources of heat and water at the earth surface, resulting in a conditional scalar transport similarity over urban terrain following rain events.

Description: The spatial and temporal dynamics of seasonal snow covers play a critical role for many hydrological, ecological, and climatic processes. This paper presents a new, innovative approach to continuously monitor these dynamics using numerous low-cost, standalone snow monitoring stations (SnoMoS). These stations provide snow and related meteorological data with a high temporal and spatial resolution. Data collected by SnoMoS include: snow depth, surface temperature, air temperature and humidity, total precipitation, global radiation, wind speed, and barometric pressure. A total of 99 sensors were placed over the winters 2010/11 and 2011/12 at multiple locations within three 40 - 180 km² basins in the Black Forest region of Southern Germany. The locations were chosen to cover a wide range of slopes, elevations, and expositions in a stratified sampling design. Furthermore, “paired stations” located in close proximity to each other, one in the open and one underneath various forest canopies, were set up to investigate the influence of vegetation on snow dynamics. The results showed that considerable differences in snow depth and, therefore, snow water equivalent (SWE) are present within the study area despite its moderate temperatures and medium elevation range (400 - 1500 m). The relative impact of topographical factors like elevation, aspect, and of different types of forest vegetation were quantified continuously and were found to change considerably over the winter period. The recorded differences in SWE and snow cover duration were large enough that they should be considered in hydrologic and climate models.

Description: Integrated land surface-groundwater models are valuable tools in simulating the terrestrial hydrologic cycle as a continuous system and exploring the extent of land surface–subsurface interactions from catchment to regional scales. However, the fidelity of model simulations is impacted not only by the vegetation and subsurface parameterizations, but also by the antecedent condition of model state variables, such as the initial soil moisture, depth to groundwater and ground temperature. In land surface modeling, a given model is often run repeatedly over a single year of forcing data until it reaches an equilibrium state: the point at which there is minimal artificial drift in the model state or prognostic variables (most often the soil moisture). For more complex coupled and integrated systems, where there is an increased computational cost of simulation and the number of variables sensitive to initialization is greater than in traditional uncoupled land surface modeling schemes, the challenge is to minimize the impact of initialization while using the smallest spin-up time possible. In this study, multi-criteria analysis was performed to assess the spin-up behavior of the ParFlow.CLM integrated groundwater-surface water-land surface model over a 208 km 2 sub-catchment of the Ringkobing Fjord catchment in Denmark. Various measures of spin-up performance were computed for model state variables such as the soil moisture and groundwater storage, as well as for diagnostic variables such as the latent and sensible heat fluxes. The impacts of initial conditions on surface water–groundwater interactions were then explored. Our analysis illustrates that the determination of an equilibrium state depends strongly on the variable and performance measure used. Choosing an improper initialization of the model can generate simulations that lead to a misinterpretation of land surface-subsurface feedback processes and result in large biases in simulated discharge. Estimated spin-up time from a series of spin-up functions revealed that 20 (or 21) years of simulation were sufficient for the catchment to equilibrate according to at least one criterion at the 0.1% (0.01%) threshold level. Amongst a range of convergence metrics examined, percentage changes in monthly values of groundwater and unsaturated zone storages produced a slow system convergence to equilibrium, whereas criteria based on ground temperature allowed a more rapid spin-up. Slow convergence of unsaturated and saturated zone storages is a result of the dynamic adjustment of the water table in response to a physically arbitrary or inconsistent initialization of a spatially uniform water table. Achieving equilibrium in subsurface storage ensured equilibrium across a spectrum of other variables, hence providing a good measure of system-wide equilibrium. Overall, results highlight the importance of correctly identifying the key variable affecting model equilibrium and also the need to use a multi-criteria approach to achieve a rapid and stable model spin-up.

Description: In this paper, we perform a detailed regional analysis of the link between meteorological drought indices and streamflow for a comprehensive Austrian dataset of 47 small to medium-size catchments in humid-temperate climate. Four drought indices considering different components of the catchment water balance are tested. We assess the quality of the link using rank correlation analysis, and the probability of detecting low flow events by hit-scores. Overall, correlations range between 0.4 and 0.8 and differ significantly between regions. A stratified analysis shows that the link is much stronger (i) for summer low flows and droughts than for anomalies in general, and (ii) for more humid than more arid conditions. Under more humid conditions streamflow droughts of small to medium-size catchments are to a large extent generated by climate forcing and therefore well represented by a simple meteorological index. Under increasingly dry conditions, the climate signal gets less predictive and storage properties of the underground become more important. A simple soil moisture accounting scheme (such as those of the Palmer index) can considerably improve the correlations. Overall, we conclude there is a significant link between meteorological drought and streamflow drought, except for catchments where groundwater storage and snow processes are important. The results are encouraging and provide a wealth of information which can profitably be used to set up statistical prediction models to estimate low flows from meteorological time series.

Description: Water resources in the western United States are contingent on interannual variations in snowpack. Interannual snowpack variability has been attributed to large-scale climate patterns including the El Niño-Southern Oscillation (ENSO), however the contribution of snowfall frequency and extreme snowfall events to this variability are less well quantified. Long term records from Snowpack Telemetry and Cooperative Observer Program stations in the eleven western states were used to investigate these relationships by considering the number of snowfall days and snowfall water equivalent (SFE) of extreme snowfall events. The top decile of snowfall events contributed 20-38% of annual SFE, depending on the region. An average of 65% and 69% of the interannual variability in annual SFE was explained by snowfall days and SFE of top decile snowfall events, respectively, with extreme events being a more significant predictor at most stations. The latitudinal dipole in SFE during ENSO phases results from changes in snowfall frequency and extreme events. In the Pacific Northwest, above normal SFE during La Niña winters was a product of both larger contributions from extremes and more snowfall days, while below normal SFE during El Niño winters was primarily associated with a substantial reduction in extremes. Conversely, annual SFE during ENSO phases in the mountains of Arizona was more closely linked to fluctuations in snowfall days than extremes. Results indicate the importance of extreme snowfall events in shaping interannual variability in water resources and suggest that improved predictive ability may inform better water resource management now and in the coming decades.

Description: The formation of a homogeneous Bi 8 TiO 14 phase was successfully achieved in a specimen calcined at 600°C. However, a Bi 4 Ti 3 O 12 secondary phase also developed in specimens calcined at temperatures higher than 600°C, probably because of Bi 2 O 3 evaporation. For specimens sintered above 800°C, a small amount of the Bi 8 TiO 14 phase melted during sintering, with the liquid phase contributing to the densification of the specimens; however, Bi 4 Ti 3 O 12 and Bi 12 TiO 20 secondary phases were still formed in these specimens. The microwave dielectric properties of the Bi 8 TiO 14 phase were considerably affected by variations in the microstructure of the specimens. When the sintering temperature exceeded 825°C, the amount of secondary phases increased, and this decreased the density and Q×f values of the specimens. Bi 8 TiO 14 ceramics sintered at 825°C exhibited promising microwave dielectric properties, with ε r = 47.4, Q×f  =   5370 GHz, and τ f = −16.01 ppm/°C.

Description: Cell-centered finite volume methods are prevailing in numerical simulation of flow in porous media. However, due to the lack of cell-centered finite volume methods for mechanics, coupled flow and deformation is usually treated either by coupled finite-volume-finite element discretizations, or within a finite element setting. The former approach is unfavorable as it introduces two separate grid structures, while the latter approach loses the advantages of finite volume methods for the flow equation. Recently we proposed a cell-centered finite volume method for elasticity. Herein we explore the applicability of this novel method to provide a compatible finite volume discretization for coupled hydro-mechanic flows in porous media. We detail in particular the issue of coupling terms, and show how this is naturally handled. Furthermore, we observe how the cell-centered finite volume framework naturally allows for modeling fractured and fracturing porous media through internal boundary conditions. We support the discussion with a set of numerical examples: The convergence properties of the coupled scheme are first investigated; Secondly, we illustrate the practical applicability of the method both for fractured and heterogeneous media.

Description: The spatiotemporal distribution of Land Surface Temperature (LST) is linked to the partitioning of the coupled surface water and energy budgets. In watersheds with a strong seasonality in precipitation and vegetation cover, the temporal evolution of LST patterns are a signature of the interactions between the land surface and atmosphere. Nevertheless, few studies have sought to understand the topographical and ecohydrological controls on LST in regions of complex terrain. Numerical watershed models, tested against spatially-distributed field and remote sensing data, can aid in linking the seasonal evolution of LST to meteorology, terrain, soil and vegetation. In this study, we use a distributed hydrologic model to explore LST patterns in a semiarid mountain basin during the transition from a dry spring to the wetter North American monsoon (NAM). By accounting for vegetation greening through remotely-sensed parameters, the model reproduces LST and surface soil moisture observations derived from ground, aircraft and satellite platforms with good accuracy at individual sites and as spatial basin patterns. Distributed simulations reveal how LST varies with elevation, slope and aspect and the role played by the seasonal vegetation canopy in cooling the land surface and increasing the spatial variability in LST. As a result, LST is shown to track well with ecosystem-specific changes in vegetation cover, evapotranspiration and soil moisture during the NAM. Furthermore, vegetation greening is shown to modulate the spatial heterogeneity of LST during the NAM that should be considered in subsequent atmospheric studies in regions of complex terrain.

Description: Models of landscape evolution or hydrological processes typically depend on the accurate determination of upslope drainage area from digital elevation data, but such calculations can be very computationally demanding when applied to high-resolution topographic data. To overcome this limitation, we propose calculating drainage area in an implicit, iterative manner using linear solvers. The basis of this method is a recasting of the flow routing problem as a sparse system of linear equations, which can be solved using established computational techniques. This approach is highly parallelizable, enabling data to be spread over multiple computer processors. Good scalability is exhibited, rendering it suitable for contemporary high-performance computing architectures with many processors, such as graphics processing units (GPUs). In addition, the iterative nature of the computational algorithms we use to solve the linear system creates the possibility of accelerating the solution by providing an initial guess, making the method well suited to iterative calculations such as numerical landscape evolution models. We compare this method with a previously proposed parallel drainage area algorithm and present several examples illustrating its advantages, including a continent-scale flow routing calculation at 3 arcsecond resolution, improvements to models of fluvial sediment yield, and acceleration of drainage area calculations in a landscape evolution model. We additionally describe a modification that allows the method to be used for parallel basin delineation.

Description: Sensing of low concentrations of two nitroaromatic compounds, 1,2-dinitrotoluene and 2-nitrophenol, is presented. The sensing mechanism is based on surface-enhanced Raman scattering (SERS) using nanostructured tin oxide as the SERS-active substrate. The SnO x nanostructures are synthesized by a simple solgel method and doped with Ag and Au. The Raman signal of a low concentration of the analyte, otherwise extremely weak, becomes significant when the analytes are attached to these substrates. Doping of SnO x nanopowders with Ag and Au leads to a further increase in the Raman intensities. This study demonstrates the scope of ceramic–metal nanocomposites as convenient solid-state SERS sensors for low-level detection.

Description: Isostatically pressed dense reaction-sintered zirconia–mullite composites were prepared from zircon flour and reactive alumina with different proportions of MgO and CaO additives. The pressed compacts were sintered at elevated temperatures with variable soaking times. Development of different phases and microstructures of the sintered compacts were analyzed to assess the relative influence of the additives on the formation of the zirconia–mullite.

Description: Minor changes to seasonal air temperature and precipitation can have a substantial impact on the availability of water resources within large watersheds. Two such watersheds, the north-flowing Mackenzie and east-flowing Saskatchewan Basins have been identified as highly vulnerable to such changes, and, therefore, selected for study as part of the Climatic Redistribution of western Canadian Water Resources (CROCWR) project. CROCWR aims to evaluate spatial and temporal changes to water resource distribution through the analysis of a suite of hydroclimatic and streamflow variables. As part of this analysis, dominant summer (May-Oct) circulation patterns at 500-hPa for 1950-2011 are identified using the method of Self-Organizing Maps (SOM). Surface climate variables associated with these patterns are then identified, including both daily air temperature and precipitation, and seasonal Standardized Precipitation-Evapotranspiration Index (SPEI) values. Statistical methods are applied to assess the relationships between dominant circulation patterns and the Southern Oscillation Index (SOI) and Pacific Decadal Oscillation (PDO). Results indicate that mid-summer (Jul-Aug) is dominated by a split-flow blocking pattern, resulting in cool (warm), wet (dry) conditions in the southern (northern) portion of the study area. By contrast, the shoulder season (May and Oct) is dominated by a trough of low-pressure over the North Pacific Ocean. The frequency of weak split-flow blocking is higher during positive SOI and negative PDO, while ridging over the western continent is more frequent during negative SOI and positive PDO. Results from this analysis increase our knowledge of processes controlling the distribution of summer water resources in western Canada. This article is protected by copyright. All rights reserved.

Description: A new solid solution of (1− x ) Pb ( Mg 1/2 W 1/2 ) O 3 – x Pb ( Zn 1/2 W 1/2 )O 3 has been prepared in the form of ceramics by solid-state reaction with composition x up to 30%. It is found that with the substitution of Zn 2+ for Mg 2+ on the B site of the of complex perovskite structure the antiferroelectric (AFE) Curie temperature T C of PMW increases from 40°C ( x  = 0) to 67°C ( x  = 30%), indicating an enhancement of antiferroelectric order, whereas, at the same time, the phase transition becomes more diffuse due to a higher degree of chemical inhomogeneity. X-ray diffraction analysis indicates that the crystal structure adopts an orthorhombic space group ( Pmcn ) with a decrease in lattice parameter a, but an increase in b and c as the Zn 2+ concentration increases. The low dielectric constant (~ 10 2 ), low dielectric loss (tanδ ≈ 10 −3 ), linear-field-induced polarization, and significantly high breakdown field (~ 125 kV/cm) at room temperature make this family of dielectric materials a promising candidate for ceramic insulators.

Description: ABSTRACT Rising development costs and growing concerns over environmental impacts have led many communities to explore more diversified water management strategies. These 'portfolio'-style approaches integrate existing supply infrastructure with other options such as conservation measures or water transfers. Diversified water supply portfolios have been shown to reduce the capacity and costs required to meet demand, while also providing greater adaptability to changing hydrologic conditions. However, this additional flexibility can also cause unexpected reductions in revenue (from conservation) or increased costs (from transfers). The resulting financial instability can act as a substantial disincentive to utilities seeking to implement more innovative water management techniques. This study seeks to design portfolios that employ financial tools (e.g. contingency funds, index insurance) to reduce fluctuations in revenues and costs, allowing these strategies to achieve improved performance without sacrificing financial stability. This analysis is applied to the development of coordinated regional supply portfolios in the 'Research Triangle' region of North Carolina, an area comprising four rapidly growing municipalities. The actions of each independent utility become interconnected when shared infrastructure is utilized to enable inter-utility transfers, requiring the evaluation of regional tradeoffs in up to five performance and financial objectives. Diversified strategies introduce significant tradeoffs between achieving reliability goals and introducing burdensome variability in annual revenues and/or costs. Financial mitigation tools can mitigate the impacts of this variability, allowing for an alternative suite of improved solutions. This analysis provides a general template for utilities seeking to navigate the tradeoffs associated with more flexible, portfolio-style management approaches.

Description: Despite significant advances during the last decades, there are still many processes related to non-equilibrium flow and transport in macroporous soil that are far from completely understood. The use of X-ray for imaging time-lapse 3-D solute transport has a large potential to help advance the knowledge in this field. We visualized the transport of potassium iodide (20 g iodide l -1 H 2 O) front through a small undisturbed soil column (height 3.8 cm, diameter 6.8 cm) under steady-state hydraulic conditions using an industrial X-ray scanner. In addition, the electrical conductivity was measured in the effluent solution during the experiment. We attained a series of seventeen 3-D difference images which we related to iodide concentrations using a linear calibration relationship. The solute transport through the soil mainly took place in two cylindrical macropores, by-passing more than 90% of the bulk soil volume during the entire experiment. From these macropores the solute diffused into the surrounding soil matrix. We illustrated the properties of the investigated solute transport by comparing it to a 1-D convective-dispersive transport and by calculating the temporal evolution of the dilution index. We furthermore showed that the tracer diffusion from one of the macropores into the surrounding soil matrix could not be exactly fitted with the cylindrical diffusion equation. We believe that similar studies will help establish links between soil structure and solute transport processes and lead to improvements in models for solute transport through undisturbed soil.

Description: A semi-analytical grid-free series solution method is presented for modeling 3-D steady-state free boundary groundwater-surface water exchange in geometrically complex stratified aquifers. Continuous solutions for pressure in the subsurface are determined semi-analytically, as is the location of the water table surface. Mass balance is satisfied exactly over the entire domain except along boundaries and interfaces between layers, where errors are shown to be acceptable. The solutions are derived and demonstrated on a number of test cases and the errors are assessed and discussed. This accurate and grid-free scheme can also be a helpful tool for providing insight into lake-aquifer and stream-aquifer interactions. Here, it is used to assess the impact of lake sediment geometry and properties on lake-aquifer interactions. Various combinations of lake sediment are considered and the appropriateness of the Dupuit-Forchheimer approximation for simulating lake bottom flux distribution is investigated. In addition, the method is applied to a test problem of surface seepage flows from a complex topographic surface; this test case demonstrated the method's efficacy for simulating physically realistic domains.

Description: This article is aimed at quantifying and discussing the relative magnitude of key terms of the equation of conservation of turbulent kinetic energy (TKE) in the inter-stem space of a flow within arrays of vertical cylinders simulating plant stems of emergent and rigid vegetation. The spatial distribution of turbulent quantities and mean flow variables are influenced by two fundamental space scales, the diameter of the stems and the local stem areal number-density. Both may vary considerably since the areal distribution of plant stems in natural systems is generally not homogeneous; they are often arranged in alternating sparse and dense patches. The magnitude of the terms of the budget of TKE in the inter-stem space has seldom been quantified experimentally and is currently not well-known. This work addresses this research need. New databases, consisting of three-component LDA velocity series and two-component PIV velocity maps, obtained in carefully controlled laboratory conditions, were used to calculate the terms of the TKE budget. The physical system comprises random arrays of rigid and emergent cylinders with longitudinally varying areal number-density. It is verified that the main source of TKE is vortex shedding from individual cylinders. The rates of production and dissipation are not in equilibrium. Regions with negative production, a previously unreported feature, are identified. Turbulent transport is particularly important along the von Kármán vortex street. Convective rate of change of TKE and pressure diffusion are most relevant in the vicinity of the cylinders.

Description: A major difficulty in modeling multiphase flow in porous media is the emergence of trapped phases. Our experiments demonstrate that gas can be trapped in either single-pores, multi-pores, or in large connected networks. These large connected clusters can comprise up to 8 grain-volumes and can contain up to 50% of the whole trapped gas volume. About 85% of the gas volume is trapped by multi-pore gas clusters. This variety of possible trapped gas clusters of different shape and volume will lead to a better process understanding of bubble-mediated mass transfer. Since multi-pore gas bubbles are in contact with the solid surface through ultra-thin adsorbed water films the interfacial area between trapped gas clusters and intergranular capillary water is only about 80% of the total gas surface. We could derive a significant (R 2 = 0.98) linear relationship between the gas-water-interface and gas saturation. We found no systematic dependency of the front velocity of the invading water phase in the velocity range from 0.1 to 0.6 cm/min corresponding to capillary numbers from 2×10 -7 to 10 -6 . Our experimental results indicate that the capillary trapping mechanism is controlled by the local pore structure and local connectivity and not by thermodynamics, i.e. by the minimum of the Free Energy , at least in the considered velocity range. Consistent with this physical picture is our finding that the trapping frequency (= bubble-size distribution) reflects the pore-size distribution for the whole range of pore radii, i.e. the capillary trapping process is determined by statistics and not by thermodynamics.

Description: Hydrological and geochemical processes controlling the pore water chemistry in a permafrost wetland, with loam overlain by sphagnum peat, were investigated. The vertical distributions of dissolved Cl, and of pore water δ 18 O, appeared unrelated to ion freeze-out and isotope ice-water fractionation processes, respectively, dismissing solute freeze-out as a main control on the water chemistry. However, concentrations of major ions, others than Cl, generally increased with depth into the active layer. A conceptual model for water and solute movement in the active layer was derived. The model indicates upwards diffusive transport of elements, released in the loam layer by mineral weathering, to the peat layer, in which lateral advective transport dominates. Active layer pore water and water of melted core sections of permafrost were of Ca-Mg-HCO 3 type (1:1:4 stoichiometry) and were subsaturated for calcite and dolomite. The results are consistent with an annual cycling of inorganic carbon species, Ca and Mg, via cryogenic carbonate precipitation during fall freeze-up and their re-dissolution following spring thaw. Similarly, elevated Fe 2+ concentrations appear to be related to cryogenic siderite formation. Pore water in the active layer showed high partial pressures CO 2 , indicating the feasibility of bubble ebullition as a greenhouse gas emission pathway from permafrost wetlands. Elevated concentrations of geogenic trace elements (Ni, Al and As) were observed, and the controlling geochemical processes are discussed. The conceptual model for water and solute movement was applied to quantify the contribution of released trace elements to a downstream lake in the permafrost catchment.

Description: Like many mountainous areas in the tropics, watersheds in the Luquillo Mountains of eastern Puerto Rico have abundant rainfall and stream discharge, and provide much of the water supply for the densely populated metropolitan areas nearby. Projected changes in regional temperature and atmospheric dynamics as a result of global warming suggest that water availability will be affected by changes in rainfall patterns. It is essential to understand the relative importance of different weather systems to water supply to determine how changes in rainfall patterns, interacting with geology and vegetation, will affect the water balance. To help determine the links between climate and water availability, stable-isotope signatures of precipitation from different weather systems were established to identify those that are most important in maintaining streamflow and groundwater recharge. Precipitation stable isotope values in the Luquillo Mountains had a large range, from fog/cloud water with δ 2 H, δ 18 O values as high as +12 ‰, −0.73 ‰, to tropical storm rain with values as low as −127 ‰, −16.8 ‰. Temporal isotope values exhibit a reverse seasonality from those observed in higher latitude continental watersheds, with higher isotopic values in the winter and lower values in the summer. Despite the higher volume of convective and low-pressure system rainfall, stable isotope analyses indicated that under the current rainfall regime, frequent trade-wind orographic showers contribute much of the groundwater recharge and stream baseflow. Analysis of rain events using 20 years of 15-minute resolution data at a mountain station (643 m) showed an increasing trend in rainfall amount, in agreement with increased precipitable water in the atmosphere, but differing from climate model projections of drying in the region. The mean intensity of rain events also showed an increasing trend. The determination of recharge sources from stable isotope tracers indicates that water supply will be affected if regional atmospheric dynamics change trade-wind orographic rainfall patterns in the Caribbean.

Description: By introducing two scalar quantities, namely, the Gini and Lorenz asymmetry coefficients, we examined their characteristics and applicability to the global analysis of changes in river flow regimes under future climate change. First, by applying these coefficients to river discharge data, we showed that various types of flow-duration curves can be interpreted quantitatively in terms of the seasonal inequality in the discharge (i.e., the unevenness of the temporal distribution of river discharge). Their statistical characteristics, based on five theoretical distribution functions frequently used in hydrological analysis, were also shown. Next, we used these coefficients to evaluate the seasonal inequality of major global rivers using the global hydrological model H08 for four 30-year time spans (1960-1989, 2010-2039, 2040-2069 and 2070-2099) under four climate-change scenarios. We used ensembles of hydrological simulation results with five general circulation models. From the analysis of the Gini coefficient, future changes in seasonal inequality show a contrasting geographical pattern: a decreasing trend at high northern latitudes and an increasing trend in most other areas. The Lorenz asymmetry coefficient shows large changes at high northern latitudes, attributable to major shifts in the flow regime accompanied by different snow-melting properties under different future climate scenarios. Although a flow-duration curve is a pictorial representation of river discharge suitable for one specific site, by depicting the geographical distribution of these two coefficients along river channels, different characteristics of flow-duration curves at different sites can be detected, even within the same river basin.

Description: Little research has examined whether forests reduce stream water eutrophication in agricultural areas during spring snowmelt periods. This study evaluated the role of forests in ameliorating deteriorated stream water quality in agricultural areas, including pasture, during snowmelt periods. Temporal variation in stream water quality at a mixed land-use basin (565 ha: pasture 13%, forestry 87%), northern Japan, was monitored for 7 years. Synoptic stream water sampling was also conducted at 16 sites across a wide range of forest and agricultural areas in a basin (18.3 km 2 ) in spring, summer, and fall. Atmospheric nitrogen (N) and phosphorus (P) deposition were measured for 4 years. The results showed that concentration pulses of nitrate (NO 3 – ), organic N (ON), and total P (TP) in stream water were observed when discharge increased during spring snowmelt. Their concentrations were high when silicate concentrations were low, suggesting surface water exported from pasture largely contributed to stream water pollution during snowmelt. Atmospheric N and P deposition (4.1 kg N ha –1 y –1 ; 0.09 kg P ha –1 y –1 , respectively) was too low to affect the background concentrations of N and P in streams from forested areas. Reduction of eutrophication caused by nutrients from pasture was mainly due to dilution by water containing low concentrations of N and P exported from forested areas, while in-stream reduction were not dominant processes. Results indicate that forests have a limited capacity to reduce the concentration pulses of N and P in stream water during snowmelt in this study basin. This article is protected by copyright. All rights reserved.

Description: A semi-analytical solution is presented for transient streaming potentials associated with flow to a pumping well in an unconfined aquifer, taking into account the effect of flow in the unsaturated zone above the water table. Flow in the unsaturated zone is modeled with a linearized form of Richards' equation using an exponential model for soil moisture retention and unsaturated hydraulic conductivity. Archie's law is invoked for unsaturated electrical conductivity. The unsaturated electrokinetic coupling coefficient is modeled with a decaying exponential, where the maximum value is at and below the water table. The coupled flow and electrokinetic problem is solved using Laplace and Hankel transforms. The results of the model predicted behavior are presented and compared to that observed in laboratory simulations of pumping tests. The early-time polarity reversal predicted the model is observable in the experiments. Other non-monotonic streaming potential behaviors predicted by the model are also evident in experimental measurements. The model is used to estimate hydraulic parameters from SP data and these compare well to those obtained from drawdown data. For example, a hydraulic conductivity of 3.6 × 10 -4 m/s is obtained from SP data compared to 3.4 × 10 -4 m/s from drawdown data.

Description: Saline conditions induce not only chemical but physical changes in swelling clays, and have a significant influence on the crack dynamics and morphology of desiccating clays. In this study, we used X-ray micro-tomography to experimentally investigate the effects of sodium chloride on the morphology and dynamics of desiccation cracks in three-dimensional mixtures of sand-bentonite slurry under varying rheological conditions. Rectangular glass containers were packed with slurries of different salt concentrations, with the top boundary exposed to air for evaporation. The growth and propagation of the cracking network that subsequently formed was visualized in 3D at multiple intervals. The characterization of cracking and branching behavior shows a high extent of localized surficial crack networks at low salinity, with a transition to less extensive but more centralize crack networks with increased salinity. The observed behavior was described in the context of the physicochemical properties of the montmorillonite clay, where shifts from an “entangled” (large platelet spacing, small pore structure) to a “stacked” (small platelet spacing, open pore structure) network influence fluid distribution and thus extent of cracking and branching behavior. This is further corroborated by vertical profiles of water distribution, which shows localized desiccation fronts that shift to uniform desaturation with increasing salt concentration. Our results provide new insights regarding the formation, dynamics, and patterns of desiccation cracks formed during evaporation from 3D saline clay structures, which will be useful in hydrological applications including water management, land surface evaporation, and subsurface contaminant transport.

Description: Climate reconstructions using tree rings and lake sediments have contributed significantly to the understanding of Holocene climates. Approaches focused specifically on reconstructing the temporal water-level response of lakes, however, are much less developed. This paper describes a statistical correlation approach based on time series with Palmer Drought Severity Index (PDSI) values derived from instrumental records or tree rings as a basis for reconstructing stage hydrographs for closed-basin lakes. We use a distributed lag correlation model to calculate a variable, ω t that represents the water level of a lake at any time t as a result of integrated climatic forcing from preceding years. The method was validated using both synthetic and measured lake-stage data and the study found that a lake's “memory” of climate fades as time passes, following an exponential-decay function at rates determined by the correlation time lag. Calculated trends in ω t for Moon Lake, Rice Lake, and Lake Mina from AD 1401 to 1860 compared well with the established chronologies (salinity, moisture, and Mg/Ca ratios) reconstructed from sediments. This method provides an independent approach for developing high-resolution information on lake behaviors in pre-instrumental times and has been able to identify problems of climate signal deterioration in sediment-based climate reconstructions in lakes with a long time lag.

Description: This study explores groundwater management policies and the effect of modeling assumptions on the projected performance of those policies. The study compares an optimal economic allocation for groundwater use subject to streamflow constraints, achieved by a central planner with perfect foresight, with a uniform tax on groundwater use and a uniform quota on groundwater use. The policies are compared with two modeling approaches, the Optimal Control Model (OCM) and the Multi-Agent System Simulation (MASS). The economic decision models are coupled with a physically based representation of the aquifer using a calibrated MODFLOW groundwater model. The results indicate that uniformly applied policies perform poorly when simulated with more realistic, heterogeneous, myopic, and self-interested agents. In particular, the effects of the physical heterogeneity of the basin and the agents undercut the perceived benefits of policy instruments assessed with simple, single-cell groundwater modeling. This study demonstrates the results of coupling realistic hydrogeology and human behavior models to assess groundwater management policies. The Republican River Basin, which overlies a portion of the Ogallala aquifer in the High Plains of the United States, is used as a case study for this analysis.