ALBERT

Description: Regional and local patterns in depth to water table, hydrochemistry and peat properties of bogs and their laggs in coastal British Columbia Hydrology and Earth System Sciences, 17, 3421-3435, 2013 Author(s): S. A. Howie and H. J. van Meerveld In restoration planning for damaged raised bogs, the lagg at the bog margin is often not given considerable weight and is sometimes disregarded entirely. However, the lagg is critical for the proper functioning of the bog, as it supports the water mound in the bog. In order to include the lagg in a restoration plan for a raised bog, it is necessary to understand the hydrological characteristics and functions of this rarely studied transition zone. We studied 13 coastal British Columbia (BC) bogs and identified two different gradients in depth to water table, hydrochemistry and peat properties: (1) a local bog expanse–bog margin gradient, and (2) a regional gradient related to climate and proximity to the ocean. Depth to water table generally increased across the transition from bog expanse to bog margin. In the bog expanse, pH was above 4.2 in the Pacific Oceanic wetland region (cooler and wetter climate) and below 4.3 in the Pacific Temperate wetland region (warmer and drier climate). Both pH and pH-corrected electrical conductivity increased significantly across the transition from bog expanse to bog margin, though not in all cases. Na + and Mg 2+ concentrations were generally highest in exposed, oceanic bogs and lower in inland bogs. Ash content in peat samples increased across the bog expanse–bog margin transition, and appears to be a useful abiotic indicator of the location of the bog margin. The observed variation in the hydrological and hydrochemical gradients across the bog expanse–bog margin transition highlights both local and regional diversity of bogs and their associated laggs.

Description: Imperfect scaling in distributions of radar-derived rainfall fields Hydrology and Earth System Sciences Discussions, 10, 11385-11422, 2013 Author(s): M. J. van den Berg, L. Delobbe, and N. E. C. Verhoest Fine scale rainfall observations for modeling exercises are often not available, but rather coarser data derived from a variety of sources are used. Effectively using these data sources in models often requires the probability distribution of the data at the applicable scale. Although numerous models for scaling distributions exist, these are often based on theoretical developments, rather than on data. In this study, we develop a model based on the α-stable distribution of rainfall fields, and tested on 5 min radar data from a Belgian weather radar. We use these data to estimate functions that describe parameters of the distribution over various scales. Moreover, we study how the mean of the distribution and the intermittency change with scale, and validate and design functions to describe the shape parameter of the distribution. This information was combined into an effective model of the distribution. Finally, the model was fitted to data from numerous storms, and the resulting parameters were compared to investigate the change in scaling behavior through time.

Description: Spatially resolved information on karst conduit flow from in-cave dye-tracing Hydrology and Earth System Sciences Discussions, 10, 11311-11335, 2013 Author(s): U. Lauber, W. Ufrecht, and N. Goldscheider Artificial tracers are powerful tools to investigate karst systems. Tracers are commonly injected into sinking streams or dolines, while springs serve as monitoring sites. The obtained flow and transport parameters represent mixed information from the vadose, epiphreatic and phreatic zones, i.e., the aquifer remains a black box. Accessible active caves constitute valuable but underexploited natural laboratories to gain detailed insights into the hydrologic functioning of the aquifer. Two multi-tracer tests in the catchment of a major karst spring (Blautopf, Germany) with injections and monitoring in two associated water caves aimed at obtaining spatially and temporally resolved information on groundwater flow in different compartments of the system. Two tracers were injected in the caves to characterize the hydraulic connections between them and with the spring. Two injections at the land surface, far from the spring, aimed at resolving the aquifer's internal drainage structure. Tracer breakthrough curves were monitored by field fluorimeters in caves and at the spring. Results demonstrate the dendritic drainage structure of the aquifer. It was possible to obtain relevant flow and transport parameters for different sections of this system. The highest mean flow velocities (275 m h −1 ) were observed in the near-spring epiphreatic section (open-channel flow), while velocities in the phreatic zone (pressurized flow) were one order of magnitude lower. Determined conduit water volumes confirm results of water balances and hydrograph analyses. In conclusion, experiments and monitoring in caves can deliver spatially resolved information on karst aquifer heterogeneity and dynamics that cannot be obtained by traditional investigative methods.

Description: On the lack of robustness of hydrologic models regarding water balance simulation – a diagnostic approach on 20 mountainous catchments using three models of increasing complexity Hydrology and Earth System Sciences Discussions, 10, 11337-11383, 2013 Author(s): L. Coron, V. Andréassian, C. Perrin, M. Bourqui, and F. Hendrickx This paper investigates the robustness of rainfall–runoff models when their parameters are transferred in time. More specifically, we studied their ability to simulate water balance on periods with different hydroclimatic characteristics. The testing procedure consisted in a series of parameter transfers between 10-yr periods and the systematic analysis of mean-volume errors. This procedure was applied to three conceptual models of different structural complexity over 20 mountainous catchments in southern France. The results showed that robustness problems are common. Errors on 10-yr-mean flows were significant for all three models and calibration periods, even when the entire record was used for calibration. Various graphical and numerical tools were used to show strong similarities between the shapes of mean flow biases calculated on a 10-yr-long sliding window when various parameter sets are used. Unexpected behavioural similarities were observed between the three models tested, considering their large differences in structural complexity. While the actual causes for robustness problems in these models remain unclear, this work stresses the limited transferability in time of the water balance adjustments made through parameter optimization. Although absolute differences between simulations obtained with different calibrated parameter sets were sometimes substantial, relative differences in simulated mean flows between time periods remained similar regardless of the calibrated parameter sets.

Description: Bayesian networks modelling in support to cross-cutting analysis of water supply and sanitation in developing countries Hydrology and Earth System Sciences, 17, 3397-3419, 2013 Author(s): C. Dondeynaz, J. López Puga, and C. Carmona Moreno Despite the efforts made towards the Millennium Development Goals targets during the last decade, improved access to water supply or basic sanitation still remains unavailable for millions of people across the world. This paper proposes a set of models that use 25 key variables and country profiles from the WatSan4Dev data set involving water supply and sanitation (Dondeynaz et al., 2012). This paper suggests the use of Bayesian network modelling methods because they are more easily adapted to deal with non-normal distributions, and integrate a qualitative approach for data analysis. They also offer the advantage of integrating preliminary knowledge into the probabilistic models. The statistical performance of the proposed models ranges between 20 and 5% error rates, which are very satisfactory taking into account the strong heterogeneity of variables. Probabilistic scenarios run from the models allow an assessment of the relationships between human development, external support, governance aspects, economic activities and water supply and sanitation (WSS) access. According to models proposed in this paper, gaining a strong poverty reduction will require the WSS access to reach 75–76% through: (1) the management of ongoing urbanisation processes to avoid slums development; and (2) the improvement of health care, for instance for children. Improving governance, such as institutional efficiency, capacities to make and apply rules, or control of corruption is positively associated with WSS sustainable development. The first condition for an increment of the HDP (human development and poverty) remains of course an improvement of the economic conditions with higher household incomes. Moreover, a significant country commitment to the environment, associated with civil society freedom of expression constitutes a favourable setting for sustainable WSS services delivery. Intensive agriculture using irrigation practises also appears as a mean for sustainable WSS thanks to multi-uses and complementarities. With a WSS sector organised at national level, irrigation practices can support the structuring and efficiency of the agriculture sector. It may then induce rural development in areas where WSS access often is set back compared to urban areas 1 . External financial support, called Official Development Assistance (ODA CI), plays a role in WSS improvement but comes last in the sensitivity analyses of models. An overall 47% of the Official Development Assistance goes first to poor countries, and is associated to governance aspects: (1) political stability and (2) country commitment to the environment and civil society degree of freedom. These governance aspects constitute a good framework for aid implementation in recipient countries. Modelling is run with the five groups of countries as defined in Dondeynaz et al. (2012). Models for profile 4 (essential external support) and profile 5 (primary material consumption) are specifically detailed and analysed in this paper. For countries in profile 4, fighting against water scarcity and progressing desertification should be the priority. However, for countries in profile 5, efforts should first concentrate on consolidation of political stability while supporting diversification of the economic activities. Nevertheless, for both profiles, reduction of poverty should remain the first priority as previously indicated. 1 JMP statistics, 2004 http://www.wssinfo.org/data-estimates/table/ , last access: 22 July 2013.

Description: Technical Note: A comparison of model and empirical measures of catchment-scale effective energy and mass transfer Hydrology and Earth System Sciences, 17, 3389-3395, 2013 Author(s): C. Rasmussen and E. L. Gallo Recent work suggests that a coupled effective energy and mass transfer (EEMT) term, which includes the energy associated with effective precipitation and primary production, may serve as a robust prediction parameter of critical zone structure and function. However, the models used to estimate EEMT have been solely based on long-term climatological data with little validation using direct empirical measures of energy, water, and carbon balances. Here we compare catchment-scale EEMT estimates generated using two distinct approaches: (1) EEMT modeled using the established methodology based on estimates of monthly effective precipitation and net primary production derived from climatological data, and (2) empirical catchment-scale EEMT estimated using data from 86 catchments of the Model Parameter Estimation Experiment (MOPEX) and MOD17A3 annual net primary production (NPP) product derived from Moderate Resolution Imaging Spectroradiometer (MODIS). Results indicated positive and significant linear correspondence ( R 2 = 0.75; P 〈 0.001) between model and empirical measures with an average root mean square error (RMSE) of 4.86 MJ m −2 yr −1 . Modeled EEMT values were consistently greater than empirical measures of EEMT. Empirical catchment estimates of the energy associated with effective precipitation ( E PPT ) were calculated using a mass balance approach that accounts for water losses to quick surface runoff not accounted for in the climatologically modeled E PPT . Similarly, local controls on primary production such as solar radiation and nutrient limitation were not explicitly included in the climatologically based estimates of energy associated with primary production ( E BIO ), whereas these were captured in the remotely sensed MODIS NPP data. These differences likely explain the greater estimate of modeled EEMT relative to the empirical measures. There was significant positive correlation between catchment aridity and the fraction of EEMT partitioned into E BIO ( F BIO ), with an increase in F BIO as a fraction of the total as aridity increases and percentage of catchment woody plant cover decreases. In summary, the data indicated strong correspondence between model and empirical measures of EEMT with limited bias that agree well with other empirical measures of catchment energy and water partitioning and plant cover.

Description: Hong-Biao Zhang, Lin Yang, and San-Xing Guo We construct so(3) algebra associated with the Morse potential and show that these operators obey so(3) commutation relations. A so(3) algebraic method is proposed in order to obtain the eigenvalues and eigenfunctions of the Morse potential. This method exhibits that Cartan operator J, the lowering ... [J. Math. Phys. 54, 092104 (2013)] published Fri Sep 6, 2013.

Description: Review Article: On the relation between the seismic activity and the Hurst exponent of the geomagnetic field at the time of the 2000 Izu swarm Natural Hazards and Earth System Science, 13, 2189-2194, 2013 Author(s): F. Masci and J. N. Thomas Many papers document the observation of earthquake-related precursory signatures in geomagnetic field data. However, the significance of these findings is ambiguous because the authors did not adequately take into account that these signals could have been generated by other sources, and the seismogenic origin of these signals have not been validated by comparison with independent datasets. Thus, they are not reliable examples of magnetic disturbances induced by the seismic activity. Hayakawa et al. (2004) claim that at the time of the 2000 Izu swarm the Hurst exponent of the Ultra-Low-Frequency (ULF: 0.001–10 Hz) band of the geomagnetic field varied in accord with the energy released by the seismicity. The present paper demonstrates that the behaviour of the Hurst exponent was insufficiently investigated and also misinterpreted by the authors. We clearly show that during the Izu swarm the changes of the Hurst exponent were strongly related to the level of global geomagnetic activity and not to the increase of the local seismic activity.

Description: Forest fire danger rating in complex topography – results from a case study in the Bavarian Alps in autumn 2011 Natural Hazards and Earth System Science, 13, 2157-2167, 2013 Author(s): C. Schunk, C. Wastl, M. Leuchner, C. Schuster, and A. Menzel Forest fire danger rating based on sparse meteorological stations is known to be potentially misleading when assigned to larger areas of complex topography. This case study examines several fire danger indices based on data from two meteorological stations at different elevations during a major drought period. This drought was caused by a persistent high pressure system, inducing a pronounced temperature inversion and its associated thermal belt with much warmer, dryer conditions in intermediate elevations. Thus, a massive drying of fuels, leading to higher fire danger levels, and multiple fire occurrences at mid-slope positions were contrasted by moderate fire danger especially in the valleys. The ability of fire danger indices to resolve this situation was studied based on a comparison with the actual fire danger as determined from expert observations, fire occurrences and fuel moisture measurements. The results revealed that, during temperature inversion, differences in daily cycles of meteorological parameters influence fire danger and that these are not resolved by standard meteorological stations and fire danger indices (calculated on a once-a-day basis). Additional stations in higher locations or high-resolution meteorological models combined with fire danger indices accepting at least hourly input data may allow reasonable fire danger calculations under these circumstances.

Description: Preface "Landslide hazard and risk assessment at different scales" Natural Hazards and Earth System Science, 13, 2169-2171, 2013 Author(s): P. Reichenbach, A. Günther, and T. Glade

Description: Tree species are predicted to track future climate by shifting their geographic distributions, but climate-mediated migrations are not apparent in a recent continental-scale analysis. To better understand the mechanisms of a possible migration lag, we analyzed relative recruitment patterns by comparing juvenile and adult tree abundances in climate space. One would expect relative recruitment to be higher in cold and dry climates as a result of tree migration with juveniles located further poleward than adults. Alternatively, relative recruitment could be higher in warm and wet climates as a result of higher tree population turnover with increased temperature and precipitation. Using the USDA Forest Service's Forest Inventory and Analysis data at regional scales, we jointly modeled juvenile and adult abundance distributions for 65 tree species in climate space of the eastern United States. We directly compared the optimal climate conditions for juveniles and adults, identified the climates where each species has high relative recruitment, and synthesized relative recruitment patterns across species. Results suggest that for 77% and 83% of the tree species, juveniles have higher optimal temperature and optimal precipitation, respectively, than adults. Across species, the relative recruitment pattern is dominated by relatively more abundant juveniles than adults in warm and wet climates. These different abundance-climate responses through life history are consistent with faster population turnover and inconsistent with the geographic trend of large-scale tree migration. Taken together, this juvenile-adult analysis suggests that tree species might respond to climate change by having faster turnover as dynamics respond to longer growing seasons and higher temperatures, before there is evidence of poleward migration at biogeographic scales. This article is protected by copyright. All rights reserved.

Description: Land use change as conversion pasture to forest produces several changes on hydrological cycle. In this paper we analyze the effects on stream discharge of afforestation of a small watershed devoted to pasture using the HBV hydrological model. Streamflow data obtained over the first ten years after planting were employed to evaluate the capacity of HBV model to simulate hydrological behavior of catchment after afforestation. Obtained results indicate that the estimation of streamflow was accurate as reflected by statistics (R 2  = 0.90, NSC = 0.89 and PBIAS = 0.34). Afterwards, streamflow under pasture land use (if afforestation had not occurred) was simulated using hydro-meteorological data collected during the period of study and model parameters optimized previously, together with two parameters, pcorr and cevpfo , that were adjusted for pasture conditions. HBV model results indicate that afforestation produced a water yield reduction around 2000 mm (22% of total stream discharge) during the first ten years of planting growth. The differences between forest and pasture land cover are increasing in all seasons year by year. The greatest streamflow reduction was observed in wet period (autumn and winter) with 76% of total reduction. In summer, streamflow reduction represents only 3% of total, however, represents 24.7% of discharge in this season. Streamflow reduction was related to increase of rainfall interception (mainly in wet periods) and the increase of evapotranspiration by plantation in dry periods. This article is protected by copyright. All rights reserved.

Description: Air flows from the atmosphere into an unconfined aquifer when the water table falls during pumping tests. Pumping test results in unconfined aquifers may be significantly affected by low-permeability zones (LPZs) near the initial water table position because they restrict the downward movement of air. A transient, three-dimensional (3D) air-water two-phase flow model is employed to investigate numerically the effects of local heterogeneity on pumping test results in unconfined aquifers. Two cases of local heterogeneities are considered herein: a LPZ around the pumping well and on one side of the pumping well. Results show that the drawdown with the LPZ is significantly greater than that of the homogeneous aquifer. The differences in drawdown are the most significant at intermediate times and gradually diminish at later times. The LPZ significantly reduces air flow from the atmosphere to the aquifer. The pore air velocity in the LPZ is very low. The air pressure at the observation point under the LPZ when air begins to enter is significantly lower than the air pressure of the homogeneous aquifer at the same point. After that, the air pressure increases quickly and then increases slowly. The time for the air pressure to reach the atmospheric pressure is significantly longer. This article is protected by copyright. All rights reserved.

Description: Over the last half century, since logging for timber production became intensive, Borneo has lost much of its pristine tropical forests. The long-term consequences of associated decline in precipitation are evident, and might, in turn, cause much more severe deforestation. This article is protected by copyright. All rights reserved.

Description: Antecedent flow conditions and nitrate concentrations in the Mississippi River Basin Hydrology and Earth System Sciences Discussions, 10, 11451-11484, 2013 Author(s): J. C. Murphy, R. M. Hirsch, and L. A. Sprague The influence of antecedent flow conditions on nitrate concentrations was explored at eight sites in the Mississippi River Basin, USA. Antecedent moisture conditions have been shown to influence nutrient export from small, relatively homogenous basins, but this influence has not been observed at a regional or continental scale. Antecedent flow conditions were quantified as the ratio between the mean daily flow of the previous year and the mean daily flow from the period of record ( Q ratio), and the Q ratio was statistically related to nitrate anomalies (the unexplained variability in nitrate concentration after filtering out season, long-term trend, and contemporaneous flow effects) at each site. Nitrate anomaly and Q ratio were negatively related at three of the four major tributary sites and upstream in the Mississippi River, indicating that when the previous year was drier than average, at these sites, nitrate concentrations were higher than expected. The strength of these relationships increased when data were subdivided by contemporaneous flow conditions. Five of the eight sites had significant negative relationships ( p ≤ 0.05) at high or moderately high contemporaneous flows, suggesting nitrate that accumulates in these basins during a drought is flushed during subsequent storm events. At half of the sites, when flow during the previous year was 50% drier than average, nitrate concentration can be from 9 and 27% higher than nitrate concentrations that follow a year with average daily flow. Conversely, nitrate concentration can be from 8 and 21% lower than expected when the previous year was 50% wetter than average. These relationships between nitrate concentration and Q ratio serve as the basis for future studies that can better define specific hydrologic processes occurring during and after a drought, which influence nitrate concentration, such as the duration or magnitude of low flows, and the timing of low and high flows.

Description: Tidally induced velocity variations of the Beardmore Glacier, Antarctica, and their representation in satellite measurements of ice velocity The Cryosphere, 7, 1375-1384, 2013 Author(s): O. J. Marsh, W. Rack, D. Floricioiu, N. R. Golledge, and W. Lawson Ocean tides close to the grounding line of outlet glaciers around Antarctica have been shown to directly influence ice velocity, both linearly and non-linearly. These fluctuations can be significant and have the potential to affect satellite measurements of ice discharge, which assume displacement between satellite passes to be consistent and representative of annual means. Satellite observations of horizontal velocity variation in the grounding zone are also contaminated by vertical tidal effects, the importance of which is highlighted here in speckle tracking measurements. Eight TerraSAR-X scenes from the grounding zone of the Beardmore Glacier are analysed in conjunction with GPS measurements to determine short-term and decadal trends in ice velocity. Diurnal tides produce horizontal velocity fluctuations of 〉50% on the ice shelf, recorded in the GPS data 4 km downstream of the grounding line. This variability decreases rapidly to 〈5% only 15 km upstream of the grounding line. Daily fluctuations are smoothed to 〈1% in the 11-day repeat pass TerraSAR-X imagery, but fortnightly variations over this period are still visible and show that satellite-velocity measurements can be affected by tides over longer periods. The measured tidal displacement observed in radar look direction over floating ice also allows the grounding line to be identified, using differential speckle tracking where phase information cannot be easily unwrapped.

Description: The 170 km long river course of the Guadalete River (western Andalucía) provides an excellent record of Late Pleistocene and Holocene fluvial sedimentation dynamics. Furthermore, its floodplain sediments are very well suited to describe geomorphic changes forced by climate fluctuations, sea-level changes, tectonic influences and human activity. Multiproxy investigations were based on field mapping and the study of 18 profile sections, mainly including sedimentological characterisation, soil-chemical analyses and radiocarbon dating of 34 samples. Findings were complemented by drillings and electrical resistivity tomography. The lowermost 50 km of the river section are divided into an upper and lower part (each 25 km long), based on different sediment preservation conditions. The boundary corresponds to the disappearance of the Late Pleistocene river terrace. Significant floodplain aggradation occurred at around 10 000 cal. years BP, while dynamics were strongly affected by sea-level fluctuations until the early Holocene. Furthermore, sedimentation starting at 8000, 6100, 4600, 2200, 900 and 400 cal. years BP is related to enhanced fluvial dynamics due to environmental stress that presumably was triggered by climate fluctuations, that is, aridification. However, the strongest intensity of sedimentation at 400 cal. years BP points to climate anomalies in the course of the Little Ice Age. In contrast, several periods of stability associated with alluvial soil formation took place during the Bølling and Allerød interstadials, prior to 8000, 6100 and 5100, and after 4300 and 2000 cal. years BP. The anthropogenic signal in floodplain evolution is not clearly distinguishable from that of climate. However, human land use had the potential to amplify geomorphic processes, especially during periods of climate deteriorations that caused increasing stress on the environment.

Description: Contribution of snow and glacier melt to discharge for highly glacierised catchments in Norway Hydrology and Earth System Sciences Discussions, 10, 11485-11517, 2013 Author(s): M. Engelhardt, T. V. Schuler, and L. M. Andreassen Glacierised catchments significantly alter the streamflow regime due to snow and glacier meltwater contribution to discharge. In this study, we modelled the mass balance and discharge rates for three highly glacierised catchments (〉50% glacier cover) in western Norway over the period 1961–2012. The spatial pattern of the catchments follows a gradient in climate continentality from west to east. The model uses gridded temperature and precipitation values from seNorge ( http://senorge.no ) as input which are available at a daily resolution. It accounts for accumulation of snow, transformation of snow to firn and ice, evaporation and melt. The model was calibrated for each catchment based on measurements of seasonal glacier mass-balances and daily discharge rates. For validation, daily melt rates were compared with measurements from sonic rangers located in the ablation zones of two of the glaciers and an uncertainty analysis was performed for the third catchment. The discharge contributions from snowmelt, glacier melt and rain were analysed with respect to spatial variations and temporal evolution. The model simulations reveal an increase of the relative contribution from glacier melt for the three catchments from less than 10% in the early 1990s to 15–30% in the late 2000s. The decline in precipitation by 10–20% in the same period was therefore overcompensated resulting in an increase of the annual discharge by 5–20%. Annual discharge sums and annual glacier melt are strongest correlated with annual and winter precipitation at the most maritime glacier and, with increased climate continentality, variations in both glacier melt contribution and annual discharge are becoming stronger correlated with variations in summer temperatures.

Description: N. Aizawa, Z. Kuznetsova, and F. Toppan Inequivalent [script N]=2 supersymmetrizations of the [script-l]-conformal Galilei algebra in d-spatial dimensions are constructed from the chiral (2, 2) and the real (1, 2, 1) basic supermultiplets of the [script N]=2 supersymmetry. For non-negative integer and half-integer [script-l], both superal ... [J. Math. Phys. 54, 093506 (2013)] published Wed Sep 11, 2013.

Description: Decadal trends in the Antarctic sea ice extent ultimately controlled by ice-ocean feedback The Cryosphere Discussions, 7, 4585-4632, 2013 Author(s): H. Goosse and V. Zunz The large natural variability of the Antarctic sea ice is a key characteristic of the system that might be responsible for the small positive trend in sea ice extent observed since 1979. In order to gain insight in the processes responsible for this variability, we have analysed in a control simulation performed with a coupled climate model a strong positive ice-ocean feedback that amplifies sea ice variations. When sea ice concentration increases in a region, in particular close to the ice edge, the mixed layer depth tends to decrease. This can be caused by a net inflow of ice and thus of freshwater that stabilizes the water column. Another stabilizing mechanism at interannual time scales that appears more widespread in our simulation is associated with the downward salt transport due to the seasonal cycle of ice formation: brine is released in winter when ice is formed and mixed over a deep layer while the freshwater flux caused by ice melting is included in a shallow layer, resulting in a net vertical transport of salt. Because of this stronger stratification due to the presence of sea ice, more heat is stored at depth in the ocean and the vertical oceanic heat flux is reduced, which contributes to maintain a higher ice extent. This positive feedback is not associated with a particular spatial pattern. Consequently, the spatial distribution of the trend in ice concentration is largely imposed by the wind changes that can provide the initial perturbation. A positive freshwater flux could alternatively be the initial trigger but the amplitude of the final response of the sea ice extent is finally set up by the amplification related to ice-ocean feedback. Initial conditions have also an influence as the chance to have a large increase in ice extent is higher if starting from a state characterized by a low value.

Description: Accepting the concept of standardization introduced by the standardized precipitation index (SPI), similar methodologies have been developed to construct some other standardized drought indices such as the standardized precipitation evapotranspiration index (SPEI). In this study, the authors provided deep insight into the SPEI and recognized potential deficiencies/limitations in relating to the climatic water balance it used. By coupling another well-known Palmer drought severity index (PDSI), we proposed a new standardized Palmer drought index (SPDI) through a moisture departure probabilistic approach, which allows multiscalar calculation for accurate temporal and spatial comparison of the hydro-meteorological conditions of different locations. Using datasets of monthly precipitation, temperature and soil available water capacity, the moisture deficit/surplus was calculated at multiple temporal scales and a couple of techniques were adopted to adjust corresponding time series to a generalized extreme value distribution out of several candidates. Results of the historical records (1900–2012) for diverse climates by multiple indices showed that the SPDI was highly consistent and correlated with the SPEI and self-calibrated PDSI (SC-PDSI) at most analyzed time scales. Furthermore, a simple experiment of hypothetical temperature and/or precipitation change scenarios also verified the effectiveness of this newly-derived SPDI index in response to climate change impacts. Being more robust and preferable in spatial consistency and comparability as well as combining the simplicity of calculation with sufficient accounting of the physical nature of water supply and demand relating to droughts, the SPDI is promising to serve as a competent reference and alternative for drought assessment and monitoring. This article is protected by copyright. All rights reserved.

Description: Contribution of directly connected and isolated impervious areas to urban drainage network hydrographs Hydrology and Earth System Sciences, 17, 3473-3483, 2013 Author(s): Y. Seo, N.-J. Choi, and A. R. Schmidt This paper addresses the mass balance error observed in runoff hydrographs in urban watersheds by introducing assumptions regarding the contribution of infiltrated rainfall from pervious areas and isolated impervious area (IIA) to the runoff hydrograph. Rainfall infiltrating into pervious areas has been assumed not to contribute to the runoff hydrograph until Hortonian excess rainfall occurs. However, mass balance analysis in an urban watershed indicates that rainfall infiltrated to pervious areas can contribute directly to the runoff hydrograph, thereby offering an explanation for the long hydrograph tail commonly observed in runoff from urban storm sewers. In this study, a hydrologic analysis based on the width function is introduced, with two types of width functions obtained from both pervious and impervious areas, respectively. The width function can be regarded as the direct interpretation of the network response. These two width functions are derived to obtain distinct response functions for directly connected impervious areas (DCIA), IIA, and pervious areas. The results show significant improvement in the estimation of runoff hydrographs and suggest the need to consider the flow contribution from pervious areas to the runoff hydrograph. It also implies that additional contribution from flow paths through joints and cracks in sewer pipes needs to be taken into account to improve the estimation of runoff hydrographs in urban catchments.

Description: Modeled hydrologic processes are represented in a set of numerical equations, the complexity of which can be measured by the total number of variables needed. A single dominant hydrologic process could control the hydrologic response of a watershed, and so the identification of the corresponding dominant variable(s) would aid in identifying a parsimonious model and in collecting more reliable data. By accounting for both model complexity and serial correlation in the variables, a model is used to identify the dominant variables for representing watershed scale streamflow, sediment transport, and phosphorus yields. Long-term water quantity and quality data was used to show that rainfall and non-linear soil water storage were the dominant variables for weekly streamflow, suspended sediment, and particulate phosphorus. Model accuracy did not consistently improve when other statistically significant variables were included. The results suggest that improved model performance may not justify the added model complexity. As such, identification of dominant variables would be the priority for developing parsimonious hydrologic models, especially at watershed scales. This article is protected by copyright. All rights reserved.

Description: The development of a glacial lake impounded along the retreating, northeastern ice margin of the Fennoscandian Ice Sheet during the last deglaciation and environmental conditions directly following the early Holocene deglaciation have been studied in NE Finland. This so-called Sokli Ice Lake has been reconstructed previously using topographic and geomorphologic evidence. In this paper a multiproxy approach is employed to study a 3-m-thick sediment succession consisting of laminated silts grading into gyttja cored in Lake Loitsana, a remnant of the Sokli Ice Lake. Variations in the sediment and siliceous microfossil records indicate distinct changes in water depth and lake size in the Loitsana basin as the Sokli Ice Lake was drained through various spillways opening up along the retreating ice front. Geochemical data (XRF core-scanning) show changes in the influence of regional catchment geochemistry (Precambrian crystalline rocks) in the glacial lake drainage area versus local catchment geochemistry (Sokli Carbonatite Massif) within the Lake Loitsana drainage area during the lake evolution. Principal component analysis on the geochemical data further suggests that grain-size is an additional factor responsible for the variability of the sediment geochemistry record. The trophic state of the lake changed drastically as a result of morphometric eutrophication once the glacial lake developed into Lake Loitsana. The AMS radiocarbon dating on tree birch seeds found in the glaciolacustrine sediment indicates that Lake Loitsana was deglaciated sometime prior to 10 700 cal. a BP showing that tree Betula was present on the deglaciated land surrounding the glacial lake. Although glacial lakes covered large areas of northern Finland during the last deglaciation, only few glaciolacustrine sediment successions have been studied in any detail. Our study shows the potential of these sediments for multiproxy analysis and contributes to the reconstruction of environmental conditions in NE Finland directly following deglaciation in the early Holocene.

Description: Bias correction can modify climate model-simulated precipitation changes without adverse affect on the ensemble mean Hydrology and Earth System Sciences Discussions, 10, 11585-11611, 2013 Author(s): E. P. Maurer and D. W. Pierce When applied to remove climate model biases in precipitation, quantile mapping can in some settings modify the simulated trends. This has important implications when the precipitation will be used to drive an impacts model that is sensitive to changes in precipitation. We use daily precipitation output from 12 general circulation models (GCMs) over the conterminous United States interpolated to a common 1° grid, and gridded observations aggregated to the same scale, to compare precipitation differences before and after quantile mapping bias correction. The change in seasonal mean (winter, DJF, and summer, JJA) precipitation between different 30-yr historical periods is compared to examine (1) the consensus among GCMs as to whether the bias correction tends to amplify or diminish their simulated precipitation trends, and (2) whether the modification of the change in precipitation tends to improve or degrade the correspondence to observed changes in precipitation for the same periods. In some cases, for a particular GCM, the trend modification can be as large as the original simulated change, though the areas where this occurs varies among GCMs so the ensemble median shows smaller trend modification. In specific locations and seasons the trend modification by quantile mapping improves correspondence with observed trends, and in others it degrades it. In the majority of the domain the ensemble median is for little effect on the correspondence of simulated precipitation trends with observed. This highlights the need to use an ensemble of GCMs rather than relying on a small number of models to estimate impacts.

Description: Marcus A. Khuri and Gilbert Weinstein In this paper, we show how a natural coupling of the Dirac equation with the generalized Jang equation leads to a proof of the rigidity statement in the positive mass theorem with charge, without the maximal slicing condition, provided a solution to the coupled system exists. ... [J. Math. Phys. 54, 092501 (2013)] published Fri Sep 13, 2013.

Description: Andre Unterberger We rely on an analysis of initial states of the wave function associated with the simplest relativistic particlethe Weyl neutrinoto construct two series of representations of SL(2,[openface R]), or of the twofold cover of this group, by operators acting on scalar functions defined on the real line. ... [J. Math. Phys. 54, 093509 (2013)] published Fri Sep 13, 2013.

Description: Linking ICT and society in early warning and adaptation to hydrological extremes in mountains Natural Hazards and Earth System Science, 13, 2253-2270, 2013 Author(s): C. de Jong The assessment of the societal impact of hydrological extremes is particularly important in mountain regions, since mountains can be considered both as the generators and victims of extreme events. ICT can provide a powerful tool for transmitting hydro-meteorological information to predict, prepare and adapt to such events. However, in remote regions, such as mountains, the poles, deserts and islands, preventive and adaptive measures are often restricted by data availability and lack and/or incoherence of data networks. This paper distinguishes between early warning of floods and droughts, emphasising the latter in particular in mountains and explores the possibilities of enhancing the role of society in data collection, the identification, activation and application of stakeholder knowledge and transferral of data from gauged to ungauged catchments.

Description: This paper proposes a method to identify blocking onset and decay by means of two stability indicators: enstrophy advection and its integral. The key to this technique is the use of local Lyapunov exponents for the barotropic vorticity equation which can be approximated by the integral of enstrophy (IRE) over a fixed, finite, region. The IRE can then be viewed as a measure of stability. However, by differentiating the IRE with respect to time, two measures of stability can be derived to assess blocking onset and decay: (i) the integral of enstrophy advection (DIRE), for which a time series is used to assess stability; and (ii) enstrophy advection, for which contours are plotted in conjunction with 500 hPa heights to locate blocking. One year of Northern Hemisphere blocking events from July 2011-July 2012 are studied to demonstrate that the integral of enstrophy advection is a useful diagnostic. In particular, time series of IRE and DIRE for four of the blocking cases are presented, while contour plots of enstrophy advection for one case are presented. In all cases studied, the diagnostics were seen to detect the instability in an incipient blocking event and in its decay.

Description: [1]  Although they are key components of the surface ocean carbon budget, physical processes inducing carbon fluxes across the mixed-layer base, i.e. subduction and obduction, have received much less attention than biological processes. Using a global model analysis of the pre-industrial ocean, physical carbon fluxes are quantified and compared to the other carbon fluxes in and out of the surface mixed-layer, i.e. air-sea CO 2 gas exchange and sedimentation of biogenic material. Model-based carbon obduction and subduction are evaluated against independent data-based estimates to the extent that was possible. We find that climatological physical fluxes of DIC are two orders of magnitude larger than the other carbon fluxes and vary over the globe at smaller spatial scale. At temperate latitudes, the subduction of DIC and to a much lesser extent (〈10%) the sinking of particles maintain CO 2 undersaturation, whereas DIC is obducted back to the surface in the tropical band (75%) and Southern Ocean (25%).At the global scale, these two large counter-balancing fluxes of DIC amount to +275.5 PgC y -1 for the supply by obduction and -264.5 PgC y -1 for the removal by subduction which is ~ 3 to 5 times larger than previous estimates. Moreover, we find that subduction of organic carbon (dissolved and particulate) represents ~ 20% of the total export of organic carbon: at the global scale, we evaluate that, of the 11 PgC y–1 of organic material lost from the surface every year, 2.1 PgC y -1 are lost through subduction of organic carbon. Our results emphasis the strong sensitivity of the oceanic carbon cycle to changes in mixed-layer depth, ocean currents and wind.

Description: Climate warming threatens to increase mass coral bleaching events, and several studies have projected the demise of tropical coral reefs this century. However, recent evidence indicates corals may be able to respond to thermal stress though adaptive processes (e.g., genetic adaptation, acclimatization, and symbiont shuffling). How these mechanisms might influence warming induced bleaching is largely unknown. This study compared how different adaptive processes could affect coral bleaching projections. We used the latest bias-corrected global sea surface temperature (SST) output from the NOAA/GFDL Earth System Model 2 (ESM2M) for the pre-industrial period though 2100 to project coral bleaching trajectories. Initial results showed that, in the absence of adaptive processes, application of a pre-industrial climatology to the NOAA Coral Reef Watch bleaching prediction method over-predicts the present day bleaching frequency. This suggests that corals may have already responded adaptively to some warming over the industrial period. We then modified the prediction method so that the bleaching threshold either permanently increased in response to thermal history (e.g., simulating directional genetic selection) or temporarily increased for 2-10 years in response to a bleaching event (e.g., simulating symbiont shuffling). A bleaching threshold that changes relative to the preceding 60 years of thermal history reduced the frequency of mass bleaching events by 20-80% compared with the ‘no adaptive response’ prediction model by 2100, depending on the emissions scenario. When both types of adaptive responses were applied, up to 14% more reef cells avoided high frequency bleaching by 2100. However, temporary increases in bleaching thresholds alone only delayed the occurrence of high frequency bleaching by ~10 years in all but the lowest emissions scenario. Future research should test the rate and limit of different adaptive responses for coral species across latitudes and ocean basins to determine if and how much corals can respond to increasing thermal stress. This article is protected by copyright. All rights reserved.

Description: The snow-masking effect of vegetation exerts strong control on albedo in northern high latitude ecosystems. Large-scale changes in the distribution and stature of vegetation in this region will thus have important feedbacks to climate. The snow-albedo feedback is controlled largely by the contrast between snow-covered and snow-free albedo (Δα), which influences predictions of future warming in coupled climate models, despite being poorly constrained at seasonal and century time scales. Here we compare satellite observations and coupled climate model representations of albedo and tree cover for the boreal and Arctic region. Our analyses reveal consistent declines in albedo with increasing tree cover, occurring south of latitudinal tree line, that are poorly represented in coupled climate models. Observed relationships between albedo and tree cover differ substantially between snow-covered and snow-free periods, and among plant functional type (PFT). Tree cover in models varies widely but surprisingly does not correlate well with model albedo. Further, our results demonstrate a relationship between tree cover and snow-albedo feedback that may be used to accurately constrain high latitude albedo feedbacks in coupled climate models under current and future vegetation distributions. This article is protected by copyright. All rights reserved.

Description: Because of global land surface warming, extreme temperature events are expected to occur more often and more intensely, affecting the growth and development of the major cereal crops in several ways, thus affecting the production component of food security. In this paper, we have identified rice and maize crop responses to temperature in different, but consistent, phenological phases and development stages. A literature review and data compilation of around 140 scientific articles have determined the key temperature thresholds and response to extreme temperature effects for rice and maize, complementing an earlier study on wheat. Lethal temperatures and cardinal temperatures, together with error estimates, have been identified for phenological phases and development stages. Following the methodology of previous work, we have collected and statistically analysed temperature thresholds of the three crops for the key physiological processes such as leaf initiation, shoot growth and root growth and for the most susceptible phenological phases such as sowing to emergence, anthesis and grain filling. Our summary shows that cardinal temperatures are conservative between studies and are seemingly well-defined in all three crops. Anthesis and ripening are the most sensitive temperature stages in rice as well as in wheat and maize. We call for further experimental studies of the effects of transgressing threshold temperatures so such responses can be included into crop impact and adaptation models. This article is protected by copyright. All rights reserved.

Description: J. A. de Azcarraga and J. M. Izquierdo Two types of higher order Lie [script-l]-ple systems are introduced in this paper. They are defined by brackets with [script-l] 〉 3 arguments satisfying certain conditions, and generalize the well-known Lie triple systems. One of the generalizations uses a construction that allows us to associate a ... [J. Math. Phys. 54, 093510 (2013)] published Tue Sep 17, 2013.

Description: Snow on the Ross Ice Shelf: comparison of reanalyses and observations from automatic weather stations The Cryosphere, 7, 1399-1410, 2013 Author(s): L. Cohen and S. Dean Snow accumulation measurements from automatic weather stations (AWS) around the Ross Ice Shelf (RIS), Antarctica, are used to provide a new set of ground-based observations which are compared to precipitation from the ECMWF ERA-Interim and NCEP/NCAR Reanalysis-2 datasets. The high temporal resolution of the AWS snow accumulation measurements allow for an event-based comparison of reanalyses precipitation to the in situ observations. Snow accumulation records from nine AWS provide multiple years of accumulation data between 2008 and 2012 over a relatively large, homogeneous region of Antarctica, and also provide the basis for a statistical evaluation of accumulation and precipitation events. The complex effects of wind on snow accumulation (which can both limit and enhance accumulation) complicate the use of the accumulation measurements, but this analysis shows that they can provide a valuable source of ground-based observations for comparisons to modelled precipitation on synoptic timescales. The analysis shows that ERA-Interim reproduces more precipitation events than NCEP-2, and these events correspond to an average 8.2% more precipitation. Significant correlations between reanalyses and AWS event sizes are seen at several stations and show that ERA-Interim consistently produces larger precipitation events than NCEP-2.

Description: The Weather Research and Forecasting (WRF) model was used to simulate the evolution of Tropical Storm Ivan (2004) in the southeast (SE) US using both the Yonsei University (YSU) and Mellor-Yamada-Janjić (MYJ) boundary layer parameterizations. In contrast with tropical cyclone (TC) simulations over the ocean, the effect of surface layer becomes secondary for a dissipating hurricane along its terrestrial track. Although these two schemes can reproduce Ivan reasonably well, our results suggest that the mixing properties for damped mechanical turbulent conditions (weakly stable) are strongly underestimated by both parameterizations. This underestimation impacts the thermodynamic properties of the storm, leading to significant differences in the storm areal extent and the simulated precipitation fields. Suggestions for further improvements are provided. An evaluation of the impact of using or not using a convective parameterization, specifically the Kain-Fritsch (KF) scheme, at 3 km grid spacing shows marginal impact on storm coverage, intensity and precipitation, except for the presence of widespread light rainfall in the Piedmont east of the mountains when the KF is employed. Analysis of the thermal structure of the simulated storm indicates that, in the inner-storm region, the KF is either not activated or primarily produces ( parameterized ) shallow convection. As a result, the net heating tendency associated with adiabatic and diabatic processes is almost unaltered inside the storm, together with a nearly equivalent surface momentum sink, leading to similar storm areal extent and intensity. Light rainfall to the east of the mountains can be due to the trigger mechanism of KF, which depends on boundary layer convergence, forcing parameterized deep convection near the coast, where surface roughness changes enhance convergence.

Description: Will warming lead to an increased use of older soil organic carbon (SOC) by microbial communities, thereby inducing C losses from C-rich alpine soils? We studied soil microbial community composition, activity and substrate use after three and four years of soil warming (+4°C, 2007-2010) at the alpine treeline in Switzerland. The warming experiment was nested in a free air CO 2 enrichment experiment using depleted 13 CO 2 (δ 13 C = –30‰, 2001-2009). We traced this depleted 13 C label in phospholipid fatty acids (PLFA) of the organic layer (0-5 cm soil depth) and in C mineralized from root-free soils to distinguish substrate ages used by soil microorganisms: fixed before 2001 (“old”), from 2001 to 2009 (“new”) or in 2010 (“recent”). Warming induced a sustained stimulation of soil respiration (+38%) without decline in mineralizable SOC. PLFA concentrations did not reveal changes in microbial community composition due to soil warming, but soil microbial metabolic activity was stimulated (+66%). Warming decreased the amount of new and recent C in the fungal biomarker 18:2ω6,9 and the amount of new C mineralized from root-free soils, implying a shift in microbial substrate use towards a greater use of old SOC. This shift in substrate use could indicate an imbalance between C inputs and outputs, which could eventually decrease SOC storage in this alpine ecosystem. This article is protected by copyright. All rights reserved.

Description: Asif Shakeel and Peter J. Love Quantum cellular automata (QCA) are models of quantum computation of particular interest from the point of view of quantum simulation. Quantum lattice gas automata (QLGA - equivalently partitioned quantum cellular automata) represent an interesting subclass of QCA. QLGA have been more deeply analyze ... [J. Math. Phys. 54, 092203 (2013)] published Thu Sep 19, 2013.

Description: Heavy precipitation events (HPEs) affect the south-eastern area of France frequently during the months of September to November. Very high amounts of rain can fall during these events, with the ensuing flash-floods causing widespread damage. The cases of the 6th of September 2010 and the 1st to the 4th of November 2011 represent the different large-scale conditions in which these episodes can occur. These HPEs are forecast with differing levels of skill by the Méso-NH model at a 2.5 km resolution. The case of the 6th of September is used to test different methods of addressing cloud physics parameterisation uncertainties. Three ensembles are constructed, where the warm process microphysical time tendencies are perturbed by different methods. Results are compared by examining the spatio-temporal distribution of the precipitation field as well as looking at ensemble statistics. The ensemble methodology which induces the most dispersion in the rainfall field is deemed the most suitable. This method is then used to examine the sensitivity of four cases from November 2011 to errors in the microphysical and turbulent parameterisations. It appears that according to the model skill for the HPE, the sensitivity to microphysical perturbations varies. Events where the model skill is high (low) show low (moderate) sensitivity. These cases show a stronger sensitivity to perturbations performed upon the turbulent tendencies, while perturbing the microphysical and turbulent tendencies together produces even further dispersion. The results show the importance and the usefulness of ensembles with perturbed physical parameterisations in the forecasting of HPEs.

Description: This paper analyses the annual mean vertical and latitudinal structure of the Brewer-Dobson circulation in the CMIP5 models. The strength of the tropical mass upwelling is found to increase at all altitudes throughout the stratosphere due to climate change. However, the width of the tropical upwelling region narrows below about 20hPa, and widens above 20hPa, suggesting different physical mechanisms may play a role in this change above and below 20hPa. In the lower stratosphere, an equatorward shift in the stationary wave critical line allows waves to propagate further into the tropics. However, in the upper stratosphere, where the behaviour is dominated by what happens during the winter, an increase in the extratropical zonal mean westerly jet leads to a reduced equatorward refraction of planetary waves. The seasonal cycle of the change in the Brewer-Dobson circulation is also considered, and differences are found in the latitudinal structure of the increased extratropical downwelling between the Northern and Southern Hemispheres in winter.

Description: Sub-seasonal forecasts have been routinely produced at ECMWF since 2002 with re-forecasts produced "on the fly" to calibrate the real-time sub-seasonal forecasts. In this study, the skill of the re-forecasts from April 2002 to March 2012 and covering a common set of years (1995 to 2001) has been evaluated. Results indicate that the skill of the ECMWF re-forecasts to predict the Madden Julian Oscillation has improved significantly since 2002, with an average gain of about 1 day of prediction skill per year. The amplitude of the MJO has also become more realistic, although the model still tends to produce MJOs which are weaker than in the ECMWF re-analysis. As a consequence, the ability of the ECMWF model to simulate realistic MJO teleconnections over the northern and southern Extratropics has improved dramatically over the 10-year period. Forecast skill scores have also improved in the Extratropics. For instance, weekly mean forecasts of the North Atlantic Oscillation Index are more skillful in recent years than ten years ago. A large part of this improvement seems to be linked to the improvements in the representation of the Madden Julian Oscillation. Skill to predict 2-metre temperature anomalies over the northern Extratropics has also improved almost continuously since 2002. Changes in the horizontal and vertical resolutions of the atmospheric model had only a small impact on the skill scores, suggesting that most of the improvements in the ECMWF sub-seasonal forecasts were due to changes in model physics which were primarily designed to improve the model climate and medium-range forecasts. The impact of changes in the data assimilation system and in the observing data has not been considered in this study, since all the re-forecasts used for this study were initialized from the same re-analysis over a common set of years.

Description: At present continental to global scale flood forecasting predicts at a point discharge, with little attention to detail and accuracy of local scale inundation predictions. Yet, inundation variables are of interest and all flood impacts are inherently local in nature. This paper proposes a large scale flood inundation ensemble forecasting model that uses best available data and modeling approaches in data scarce areas. The model was built for the Lower Zambezi River to demonstrate current flood inundation forecasting capabilities in large data-scarce regions. ECMWF ensemble forecast (ENS) data were used to force the VIC (Variable Infiltration Capacity) hydrological model, which simulated and routed daily flows to the input boundary locations of a 2-D hydrodynamic model. Efficient hydrodynamic modeling over large areas still requires model grid resolutions that are typically larger than the width of channels that play a key a role in flood wave propagation. We therefore employed a novel sub-grid channel scheme to describe the river network in detail whilst representing the floodplain at an appropriate scale. The modeling system was calibrated using channel water levels from satellite laser altimetry and then applied to predict the February 2007 Mozambique floods. Model evaluation showed that simulated flood edge cells were within a distance of between one and two model resolutions compared to an observed flood edge and inundation area agreement was on average 86%. Our study highlights that physically plausible parameter values and satisfactory performance can be achieved at spatial scales ranging from tens to several hundreds of thousands of km 2 and at model grid resolutions up to several km 2 .

Description: A class of capillary flows in unsaturated porous media is characterized by quasi-steady viscous flow confined behind curved air-water interfaces and within liquid bodies held by capillary forces along crevices and grain contacts. The geometry of the connected capillary liquid network within the pore space resembles channels that form between adjacent bubbles in foam (Plateau borders) with solid grains representing gas bubbles in foam. For simplified channel geometry we combine expressions for viscous flow with continuity considerations to describe the evolution of the channels cross-sectional area during gravity drainage. This formulation enables modeling of unsaturated flow without invoking the Richards equation and associated hydraulic functions. We adapt a formalism originally developed for foam “free drainage” (drainage under gravity) or “forced drainage” (infiltration front motion) to a class of unsaturated flows in porous media that require a few input parameters only (mean channel corner angle, air entry value and porosity) for certain initial and boundary conditions. We demonstrate that the reduction in capillary channel cross section yields a consistent description of self-regulating internal fluxes towards attainment of the so-called “field capacity” in soil and provides an alternative method for interpretation of outflow experiments for prescribed pressure boundary conditions. Additionally, the geometrically-explicit formulation provides a more intuitive picture of capillary flows across textural boundaries (changes in channel cross-section and number of channels). The foam drainage methodology expands the range of tools available for analyses of unsaturated flow processes and offers more realistic links between liquid configuration and flow dynamics in unsaturated porous media.

Description: Synthetic streamflows at different sites in a river basin are needed for planning, operation and management of water resources projects. Modeling the temporal and spatial dependence structure of monthly streamflow at different sites is generally required. In this study, the maximum entropy copula method is proposed for multisite monthly streamflow simulation, in which the temporal and spatial dependence structure is imposed as constraints to derive the maximum entropy copula. The monthly streamflows at different sites are then generated by sampling from the conditional distribution. A case study for the generation of monthly streamflow at three sites in the Colorado River basin illustrates the application of the proposed method. Simulated streamflow from the maximum entropy copula is in satisfactory agreement with observed streamflow.

Description: The Krycklan Catchment Study (KCS) provides a unique field infrastructure for hillslope to landscape-scale research on short and long-term ecosystem dynamics in boreal landscapes. The site is designed for process-based research assessing the role of external drivers including forest management, climate change, and long-range pollutant transport on forests, mires, soils, streams, lakes and groundwater. The over-arching objectives of KCS are to (1) provide a state-of-the-art infrastructure for experimental and hypothesis driven research, (2) maintain a collection of high quality, long-term climatic, biogeochemical, hydrological and environmental data, and (3) support the development of models and guidelines for research, policy and management.

Description: 3-D Hydraulic Tomography (3-D HT) is a method for aquifer characterization whereby the 3-D spatial distribution of aquifer flow parameters (primarily hydraulic conductivity, K) is estimated by joint inversion of head change data from multiple partially-penetrating pumping tests. While performance of 3-D HT has been studied extensively in numerical experiments, few field studies have demonstrated the real-world performance of 3-D HT. Here we report on a 3-D transient hydraulic tomography (3-D THT) field experiment at the Boise Hydrogeophysical Research Site which is different from prior approaches in that it represents a “baseline” analysis of 3-D THT performance using only a single arrangement of a central pumping well and 5 observation wells with nearly complete pumping and observation coverage at 1m intervals. We jointly analyze all pumping tests using a geostatistical approach based on the quasi-linear estimator of kitanidis [1995]. We reanalyze the system after progressively removing pumping and/or observation intervals; significant progressive loss of information about heterogeneity is quantified as reduced variance of the K field overall, reduced correlation with slug test K estimates at wells, and reduced ability to accurately predict independent pumping tests. We verify that imaging accuracy is strongly improved by pumping and observational densities comparable to the aquifer heterogeneity geostatistical correlation lengths. Discrepancies between K profiles at wells, as obtained from HT and slug tests, are greatest at the tops and bottoms of wells where HT observation coverage was lacking.

Description: ABSTRACT The value of a formally defined Anthropocene for geomorphologists is discussed. Human impacts have been diachronistic, multifaceted and episodic, as demonstrated by the record of alluvial deposition in the UK. Rather than boxing time into discrete eras or periods, modern research uses calendar dates and multiple dating techniques to explore co-trajectories for a range of human impacts. Despite the value of ‘The Anthropocene’ as an informal concept and as a prompt to useful debate, arriving at a single, generally acceptable  formal definition is impractical, and has some disadvantages. This article is protected by copyright. All rights reserved.

Description: Circumpolar expansion of tall shrubs and trees into Arctic tundra is widely thought to be occurring as a result of recent climate warming, but little quantitative evidence exists for northern Siberia, which encompasses the world's largest forest-tundra ecotonal belt. We quantified changes in tall shrub and tree canopy cover in eleven, widely-distributed Siberian ecotonal landscapes by comparing very-high-resolution photography from the Cold War-era “Gambit” and “Corona” satellite surveillance systems (1965-1969) with modern imagery. We also analyzed within-landscape patterns of vegetation change to evaluate the susceptibility of different landscape components to tall shrub and tree increase. The total cover of tall shrubs and trees increased in nine of eleven ecotones. In northwest Siberia, alder ( Alnus ) shrubland cover increased 5.3 – 25.9% in five ecotones. In Taymyr and Yakutia, larch ( Larix ) cover increased 3.0 – 6.7% within three ecotones, but declined 16.8% at a fourth ecotone due to thaw of ice-rich permafrost. In Chukotka, the total cover of alder and dwarf pine ( Pinus ) increased 6.1% within one ecotone and was little-changed at a second ecotone. Within most landscapes, shrub and tree increase was linked to specific geomorphic settings, especially those with active disturbance regimes such as permafrost patterned-ground, floodplains, and colluvial hillslopes. Mean summer temperatures increased at most ecotones since the mid-1960s, but rates of shrub and tree canopy cover expansion were not strongly correlated with temperature trends and were better correlated with mean annual precipitation. We conclude that shrub and tree cover is increasing in tundra ecotones across most of northern Siberia, but rates of increase vary widely regionally and at the landscape-scale. Our results indicate that extensive changes can occur within decades in moist, shrub-dominated ecotones, as in northwest Siberia, while changes are likely to occur much more slowly in the highly continental, larch-dominated ecotones of central and eastern Siberia. This article is protected by copyright. All rights reserved.

Description: The adaptation of different species to warming temperatures has been increasingly studied. Moose ( Alces alces ) is the largest of the ungulate species occupying the northern latitudes across the globe, and in Finland it is the most important game species. It is very well adapted to severe cold temperatures, but has a relatively low tolerance to warm temperatures. Previous studies have documented changes in habitat use by moose due to high temperatures. In many of these studies the used areas have been classified according to how much thermal cover they were assumed to offer based on satellite/aerial imagery data. Here, we identified the vegetation structure in the areas used by moose under different thermal conditions. For this purpose we used airborne laser scanning (ALS) data extracted from the locations of GPS-collared moose. This provided us with detailed information about the relationships between moose and the structure of forests it uses in different thermal conditions and we were therefore able to determine and differentiate between the canopy structures at locations occupied by moose during different thermal conditions. We also discovered a threshold beyond which moose behaviour began to change significantly: as day temperatures began to reach 20 ○ C and higher, the search for areas with higher and denser canopies during daytime became evident. The difference was clear when compared to habitat use at lower temperatures, and was so strong that it provides supporting evidence to previous studies, suggesting that moose are able to modify their behaviour to cope with high temperatures, but also that the species is likely to be affected by warming climate. This article is protected by copyright. All rights reserved.

Description: Mountain ecosystems are particularly susceptible to climate change. Characterizing intraspecific variation of alpine plants along elevational gradients is crucial for estimating their vulnerability to predicted changes. Environmental conditions vary with elevation, which might influence plastic responses and affect selection pressures that lead to local adaptation. Thus, local adaptation and phenotypic plasticity among low and high elevation plant populations in response to climate, soil and other factors associated with elevational gradients might underlie different responses of these populations to climate warming. Using a transplant experiment along an elevational gradient, we investigated reproductive phenology, growth and reproduction of the nutrient-poor grassland species Ranunculus bulbosus , Trifolium montanum , and Briza media . Seeds were collected from low and high elevation source populations across the Swiss Alps and grown in nine common gardens at three different elevations with two different soil depths. Despite genetic differentiation in some traits, the results revealed no indication of local adaptation to the elevation of population origin. Reproductive phenology was advanced at lower elevation in low and high elevation populations of all three species. Growth and reproduction of T. montanum and B. media were rarely affected by garden elevation and soil depth. In R. bulbosus , however, growth decreased and reproductive investment increased at higher elevation. Furthermore, soil depth influenced growth and reproduction of low elevation R. bulbosus populations. We found no evidence for local adaptation to elevation of origin and hardly any differences in the responses of low and high elevation populations. However, the consistent advanced reproductive phenology observed in all three species shows that they have the potential to plastically respond to environmental variation. We conclude that populations might not be forced to migrate to higher elevations as a consequence of climate warming, as plasticity will buffer the detrimental effects of climate change in the three investigated nutrient-poor grassland species. This article is protected by copyright. All rights reserved.

Description: The impact of climate change on the stability of soil organic carbon (SOC) remains a major source of uncertainty in predicting future changes in atmospheric CO 2 levels. One unsettled issue is whether the mineralization response to temperature depends on SOC mineralization rate. Long-term (〉25 years) bare fallow experiments (LTBF) in which the soil is kept free of any vegetation and organic inputs, and their associated archives of soil samples represent a unique research platform to examine this issue as with increasing duration of fallow, the lability of remaining total SOC decreases. We retrieved soils from LTBF experiments situated at Askov (Denmark), Grignon (France), Ultuna (Sweden) and Versailles (France) and sampled at the start of the experiments and after 25, 50, 52, and 79 years of bare fallow, respectively. Soils were incubated at 4, 12, 20 and 35 °C and the evolved CO 2 monitored. The apparent activation energy ( Ea ) of SOC was then calculated for similar loss of CO 2 at the different temperatures. The Ea was always higher for samples taken at the end of the bare-fallow period, implying a higher temperature sensitivity of stable C than of labile C. Our results provide strong evidence for a general relationship between temperature sensitivity and SOC stability upon which significant improvements in predictive models could be based. This article is protected by copyright. All rights reserved.

Description: [1]  In recent decades, it has been found useful to partition the pelagic environment using the concept of biogeochemical provinces, or BGCPs, within each of which it is assumed that environmental conditions are distinguishable and unique at global scale. The boundaries between provinces respond to features of physical oceanography and, ideally, should follow seasonal and inter-annual changes in ocean dynamics. But this ideal has not been fulfilled except for small regions of the oceans. Moreover, BGCPs have been used only as static entities having boundaries that were originally established to compute global primary production. In the present study, a new statistical methodology based on non-parametric procedures is implemented to capture the environmental characteristics within 56 BGCPs. Four main environmental parameters (bathymetry, chlorophyll-a concentration, surface temperature and salinity) are used to infer the spatial distribution of each BGCP over 1997–2007. The resulting dynamic partition allows us to integrate changes in the distribution of BGCPs at seasonal and inter-annual timescales, and so introduces the possibility of detecting spatial shifts in environmental conditions.

Description: The temperature dependence of aerobic scope has been suggested to be a major determinant of how marine animals will cope with future rises in environmental temperature. Here we present data suggesting that in some animals, the temperature dependence of anaerobic scope (i.e. the capacity for surviving severe hypoxia) may determine present-day latitudinal distributions and potential for persistence in a warmer future. As a model for investigating the role of anaerobic scope, we studied two sibling species of coral-dwelling gobies, Gobiodon histrio and G. erythrospilus , with different latitudinal distributions, but which overlap in equal abundance at Lizard Island (14°40'S) on the Great Barrier Reef. These species did not differ in the temperature dependence of resting oxygen consumption or critical oxygen concentration (the lowest oxygen level where resting oxygen consumption can be maintained). By contrast, the more equatorial species ( G. histrio ) had a better capacity to endure anaerobic conditions at oxygen levels below the critical oxygen concentration at the high temperatures (32 – 33 °C) more likely to occur near the equator, or in a warmer future. These results suggest that anaerobic scope, in addition to aerobic scope, could be important in determining the impacts of global warming on some marine animals. This article is protected by copyright. All rights reserved.

Description: Most North American forests are at some stage of post-disturbance regrowth, subject to a changing climate, and exhibit growth and mortality patterns that may not be closely coupled to annual environmental conditions. Distinguishing the possibly interacting effects of these processes is necessary to put short-term studies in a longer-term context, and particularly important for the carbon-dense, fire-prone boreal forest. The goals of this study were to combine dendrochronological sampling, inventory records, and machine-learning algorithms to understand how tree growth and death have changed at one highly studied site (Northern Old Black Spruce, NOBS) in the central Canadian boreal forest. Over the 1999-2012 inventory period, mean tree diameter increased even as stand density and basal area declined significantly. Tree mortality averaged 1.4±0.6% yr −1 , with most mortality occurring in medium-sized trees; new recruitment was minimal. There have been at least two, and probably three, significant influxes of new trees since stand initiation, but none in recent decades. A combined tree ring chronology constructed from sampling in 2001, 2004, and 2012 showed several periods of extreme growth depression, with increased mortality lagging depressed growth by ~5 years. Higher minimum and maximum air temperatures exerted a negative influence on tree growth, while precipitation and climate moisture index had a positive effect; both current- and previous-year data exerted significant effects. Models based on these variables explained 23-44% of the ring-width variability. We suggest that past climate extremes led to significant mortality still visible in the current forest structure, with decadal dynamics superimposed on slower patterns of fire and succession. These results have significant implications for our understanding of previous work at NOBS, the carbon sequestration capability of old-growth stands in a disturbance-prone landscape, and the sustainable management of regional forests in a changing climate. This article is protected by copyright. All rights reserved.

Description: The 20th century was a pivotal period at high northern latitudes as it marked the onset of a rapid climatic warming brought on by major anthropogenic changes in global atmospheric composition. In parallel, Arctic sea ice extent has been decreasing over the period of available satellite data record. Here we document how these changes influenced vegetation productivity in adjacent eastern boreal North America. To do this, we used normalized difference vegetation index (NDVI) data, model simulations of net primary productivity (NPP), and tree-ring width measurements covering the last 300 years. Climatic and proxy-climatic datasets were used to explore the relationships between vegetation productivity and Arctic sea ice concentration and extent, and temperatures. Results indicate that an unusually large amount of black spruce ( Picea mariana ) trees entered into a period of growth decline during the late 20th century (68% of sampled trees; n = 724 cross-sections of age 〉 70 years). This finding is coherent with evidence encoded in NDVI and simulated NPP data. Analyses of climatic and vegetation productivity relationships indicate that the influence of recent climatic changes in the studied forests has been via the enhanced moisture stress (i.e. greater water demands) and autotrophic respiration amplified by the declining sea ice concentration in the Hudson Bay and Hudson Strait. The recent decline strongly contrasts with other growth reduction events that occurred during the 19 th century, which were associated with cooling and high sea ice severity. The recent decline of vegetation productivity is the first one to occur under circumstances related to excess heat in a 300-year period, and further culminates with an intensifying wildfire regime in the region. Our results concur with observations from other forest ecosystems about intensifying temperature-driven drought stress and tree mortality with ongoing climatic changes. This article is protected by copyright. All rights reserved.

Description: We applied graphical methods and multivariate statistics to understand impacts of an unsewered slum catchment on nutrients and hydrochemistry of groundwater in Kampala, Uganda. Data was collected from 56 springs (groundwater), 22 surface water sites and 13 rain samples. Groundwater was acidic and dominated by Na, Cl and NO 3 . These ions were strongly correlated indicating pollution originating from wastewater infiltration from on-site sanitation systems. Results also showed that rain, which was acidic, impacted on groundwater chemistry. Using Q-mode hierarchical cluster analysis, we identified three distinct water quality groups. The first group had springs dominated by Ca-Cl-NO 3 , low values of EC, pH and cations, and relatively high NO 3 values. These springs were shown to have originated from the acidic rains because their chemistry closely corresponded to ion concentrations that would occur from rainfall recharge, which was around 3.3 times concentrated by evaporation. The second group had springs dominated by Na-K-Cl-NO 3 and Ca-Cl-NO 3 , low pH but with higher values of EC, NO 3 and cations. We interpreted these as groundwater affected by both acid rain and infiltration of wastewater from urban areas. The third group had the highest EC values (average of 688 μS/cm), low pH and very high concentrations of NO 3 (average of 2.15 mmol/L) and cations. These springs were exclusively located in slum areas and we interpreted these springs as groundwater affected by infiltration of wastewater from poorly sanitized slums areas. Surface water was slightly reducing and eutrophic due to wastewater effluents, but the contribution of groundwater to nutrients in surface water was minimal because o-PO 4 was absent whereas NO 3 was lost by denitification. Our findings suggest that groundwater chemistry in the catchment is strongly influenced by anthropogenic inputs derived from nitrogen-containing rains and domestic wastewater. This article is protected by copyright. All rights reserved.

Description: The environmental impact of the Puyehue–Cordon Caulle 2011 volcanic eruption on Buenos Aires Natural Hazards and Earth System Science, 13, 2319-2330, 2013 Author(s): G. B. Raga, D. Baumgardner, A. G. Ulke, M. Torres Brizuela, and B. Kucienska On 4 June 2011, the volcanic complex Puyehue–Cordon Caulle located in the Chilean Andes erupted, producing a plume of gases and particles that eventually circled the Southern Hemisphere, disrupting air travel and depositing ash in large quantities. On eight occasions, the plume passed over the city of Buenos Aires, Argentina, leading local authorities to close the two international airports. The eruption occurred during an on-going field campaign when measurements of the properties of atmospheric aerosol particles were being made in Buenos Aires as part of a year-long study of the concentration and optical properties of aerosol at one site in the city. The suite of instruments deployed in Buenos Aires were not tailored to measurements of volcanic ash, but were designed to characterize urban conditions. Nevertheless, these measurements were analysed for periods when vertical profiles of aerosol backscatter, made with a ceilometer, clearly showed the presence of the volcano plume over the research site and resulted in airport closure. Aerosol optical thickness derived from AERONET, MODIS and a ceilometer at our research site, all show enhanced values clearly indicating that the three platforms identified the volcanic plume simultaneously. However, a quantitative comparison of the different estimates proves difficult, suggesting large spatial and temporal variability of the plume. Our results indicate that the number concentration of condensation nuclei (CN), the mass concentration of particle-bound polycyclic aromatic hydrocarbons (PPAH) and the light absorption coefficient exceeded the average background values by more than one standard deviation during the events of volcanic plume. The anomalous concentrations of CN suggest new particle formation, presumably from the conversion of SO 2 , while the anomalous concentrations of PPAH may come from the uptake of PAHs on the plume particles or from chemical reactions on the surface of plume particles. The anomalous absorption coefficients indicate that plume particles may contain certain compounds that can absorb radiation at 550 nm. Another possible explanation consistent with the observations is the scavenging of black carbon from urban sources as the plume descends through the boundary layer to the surface. In addition, the volcanic plume influenced the local meteorology resulting in a decrease of the temperature when compared to the average temperature during days with no plume present.

Description: Contribution of land use changes to future flood damage along the river Meuse in the Walloon region Natural Hazards and Earth System Science, 13, 2301-2318, 2013 Author(s): A. Beckers, B. Dewals, S. Erpicum, S. Dujardin, S. Detrembleur, J. Teller, M. Pirotton, and P. Archambeau Managing flood risk in Europe is a critical issue because climate change is expected to increase flood hazard in many european countries. Beside climate change, land use evolution is also a key factor influencing future flood risk. The core contribution of this paper is a new methodology to model residential land use evolution. Based on two climate scenarios ("dry" and "wet"), the method is applied to study the evolution of flood damage by 2100 along the river Meuse. Nine urbanization scenarios were developed: three of them assume a "current trend" land use evolution, leading to a significant urban sprawl, while six others assume a dense urban development, characterized by a higher density and a higher diversity of urban functions in the urbanized areas. Using damage curves, the damage estimation was performed by combining inundation maps for the present and future 100 yr flood with present and future land use maps and specific prices. According to the dry scenario, the flood discharge is expected not to increase. In this case, land use changes increase flood damages by 1–40%, to €334–462 million in 2100. In the wet scenario, the relative increase in flood damage is 540–630%, corresponding to total damages of €2.1–2.4 billion. In this extreme scenario, the influence of climate on the overall damage is 3–8 times higher than the effect of land use change. However, for seven municipalities along the river Meuse, these two factors have a comparable influence. Consequently, in the "wet" scenario and at the level of the whole Meuse valley in the Walloon region, careful spatial planning would reduce the increase in flood damage by no more than 11–23%; but, at the level of several municipalities, more sustainable spatial planning would reduce future flood damage to a much greater degree.

Description: Dual-polarisation radar measurements provide valuable information about the shapes and orientations of atmospheric ice particles. For quantitative interpretation of these data in the Rayleigh regime, common practice is to approximate the true ice crystal shape with that of a spheroid. Calculations using the discrete dipole approximation for a wide range of crystal aspect ratios demonstrate that approximating hexagonal plates as spheroids leads to significant errors in the predicted differential reflectivity, by as much as 1.5dB. An empirical modification of the shape factors in Gans's spheroid theory was made using the numerical data. The resulting simple expressions, like Gans's theory, can be applied to crystals in any desired orientation, illuminated by an arbitrarily polarised wave, but are much more accurate for hexagonal particles. Calculations of the scattering from more complex branched and dendritic crystals indicate that these may be accurately modelled using the new expression, but with a reduced permittivity dependent on the volume of ice relative to an enclosing hexagonal prism.

Description: The parameterisation of diabatic processes in numerical models is critical for the accuracy of weather forecasts and for climate projections. A novel approach to the evaluation of these processes in models is introduced in this contribution. The approach combines a suite of on-line tracer diagnostics with off-line trajectory calculations. Each tracer tracks accumulative changes in potential temperature associated with a particular parameterised diabatic process in the model. A comparison of tracers therefore allows the identification of the most active diabatic processes and their downstream impacts. The tracers are combined with trajectories computed using model-resolved winds, allowing the various diabatic contributions to be tracked back to their time and location of occurrence. We have used this approach to investigate diabatic processes within a simulated extratropical cyclone. We focus on the warm conveyor belt, in which the dominant diabatic contributions come from large-scale latent heating and parameterised convection. By contrasting two simulations, one with standard convection parameterisation settings and another with reduced parameterised convection, the effects of parameterised convection on the structure of the cyclone have been determined. Under reduced parameterised convection conditions, the large-scale latent heating is forced to release convective instability that would otherwise have been released by the convection parameterisation. Although the spatial distribution of precipitation depends on the details of the split between parameterised convection and large-scale latent heating, the total precipitation amount associated with the cyclone remains largely unchanged. For reduced parameterised convection, a more rapid and stronger latent heating episode takes place as air ascends within the warm conveyor belt.

Description: Probabilistic estimates of future water levels and river discharge are usually simulated with hydrologic models using ensemble weather forecasts as main inputs. As hydrologic models are imperfect and the meteorological ensembles tend to be biased and underdispersed, the ensemble forecasts for river runoff typically are biased and underdispersed, too. Thus, in order to achieve both reliable and sharp predictions statistical post-processing is required. In this work Bayesian model averaging (BMA) is applied to statistically post-process ensemble runoff raw forecasts for a catchment in Switzerland, at lead-times ranging from 1 to 240 hours. The raw forecasts have been obtained using deterministic and ensemble forcing meteorological models with different forecast lead-time ranges. First, BMA is applied based on mixtures of univariate normal distributions, subject to the assumption of independence between distinct lead-times. Then, the independence assumption is relaxed in order to estimate multivariate runoff forecasts over the entire range of lead-times simultaneously, based on a BMA version that uses multivariate normal distributions. Since river runoff is a highly skewed variable, Box-Cox transformations are applied in order to achieve approximate normality. Both univariate and multivariate BMA approaches are able to generate well calibrated probabilistic forecasts that are considerably sharper than climatological forecasts. Additionally, multivariate BMA provides a promising approach for incorporating temporal dependencies into the post-processed forecasts. Its major advantage against univariate BMA is an increase in reliability when the forecast system is changing due to model availability.

Description: A general treatment of snow microstructure exemplified by an improved relation for thermal conductivity The Cryosphere, 7, 1473-1480, 2013 Author(s): H. Löwe, F. Riche, and M. Schneebeli Finding relevant microstructural parameters beyond density is a longstanding problem which hinders the formulation of accurate parameterizations of physical properties of snow. Towards a remedy, we address the effective thermal conductivity tensor of snow via anisotropic, second-order bounds. The bound provides an explicit expression for the thermal conductivity and predicts the relevance of a microstructural anisotropy parameter Q , which is given by an integral over the two-point correlation function and unambiguously defined for arbitrary snow structures. For validation we compiled a comprehensive data set of 167 snow samples. The set comprises individual samples of various snow types and entire time series of metamorphism experiments under isothermal and temperature gradient conditions. All samples were digitally reconstructed by micro-computed tomography to perform microstructure-based simulations of heat transport. The incorporation of anisotropy via Q considerably reduces the root mean square error over the usual density-based parameterization. The systematic quantification of anisotropy via the two-point correlation function suggests a generalizable route to incorporate microstructure into snowpack models. We indicate the inter-relation of the conductivity to other properties and outline a potential impact of Q on dielectric constant, permeability and adsorption rate of diffusing species in the pore space.

Description: We derive a series solution for the nonlinear Boussinesq equation in terms of the similarity variable of the Boltzmann transformation in a semi-infinite domain. The first few coefficients of the series have been known for a long time, having been obtained by a truncated inversion of the series solution of the Blasius equation, but no direct recurrence relation was known for the complete series representing the solution of the Boussinesq equation. The series turns out to have a finite radius of convergence, which we estimate with a numerical complex-plane integration method that identifies the singularities of the solution when the equation is extended to the complex plane. The homogeneous condition at the origin produces a singularity which complicates numerical solutions with Runge-Kutta methods. We present two variable transformations that circumvent the problem and that are best suited to the complex-variable and the real-variable versions of the equation, respectively. Using those tools, an approximate solution accurate to 1.75 × 10 -10 and valid for the entire positive real axis is then developed by matching a Padé approximant of the exact series and an asymptotic solution (to overcome the restriction imposed by the finite radius of convergence of the series), along the same lines of the expression proposed by Hogarth and Parlange [1999]. The accuracies of all of the existing and the newly proposed solutions are obtained.

Description: Field sampling in unwadeable and flashy flood events encounters the problem that lateral variability of flow hydraulics and sediment transport cannot be captured adequately, and there is also an accuracy problem because parameters change while being measured. Moreover, event based gravel-sand mixed transport data in rapidly changing conditions are largely missing, in particular for gravel-bed rivers in small catchments. In this study, field measurements of bed load, suspended load, flow velocities, water depths and cross section geometry were collected during flood events at a monitoring station near the mouth of the Versilia river, Italy. Since the observed hydrographs are characterized by short durations, to the order of a few hours, an analysis of the lateral and temporal flow variability was carried out to enable the design of a sampling strategy and to minimize the errors created by the time variations of discharge associated with unsteady flow conditions. The measurements were interpreted using a 1D hydro-morphodynamic numerical model simulating the dynamics of flow and sediment discharges during a flood event for a given return period. The flow and sediment rating curves were then developed through an integrated approach combining different methodologies: field measurements, laboratory analyses and mathematical modeling. The developed approach allows one to capture the main physical mechanisms associated to the transport of sand–gravel mixtures, such as selective transport, and the hysteretic behaviour of sediment transport produced by rapid and intense flood events.

Description: Influence of high-order mechanics on simulation of glacier response to climate change: insights from Haig Glacier, Canadian Rocky Mountains The Cryosphere, 7, 1527-1541, 2013 Author(s): S. Adhikari and S. J. Marshall Evolution of glaciers in response to climate change has mostly been simulated using simplified dynamical models. Because these models do not account for the influence of high-order physics, corresponding results may exhibit some biases. For Haig Glacier in the Canadian Rocky Mountains, we test this hypothesis by comparing simulation results obtained from 3-D numerical models that deal with different assumptions concerning physics, ranging from simple shear deformation to comprehensive Stokes flow. In glacier retreat scenarios, we find a minimal role of high-order mechanics in glacier evolution, as geometric effects at our site (the presence of an overdeepened bed) result in limited horizontal movement of ice (flow speed on the order of a few meters per year). Consequently, high-order and reduced models all predict that Haig Glacier ceases to exist by ca. 2080 under ongoing climate warming. The influence of high-order mechanics is evident, however, in glacier advance scenarios, where ice speeds are greater and ice dynamical effects become more important. Although similar studies on other glaciers are essential to generalize such findings, we advise that high-order mechanics are important and therefore should be considered while modeling the evolution of active glaciers. Reduced model predictions may be adequate for other glaciologic and topographic settings, particularly where flow speeds are low and where mass balance changes dominate over ice dynamics in determining glacier geometry.

Description: Meteorological drivers of ablation processes on a cold glacier in the semi-arid Andes of Chile The Cryosphere, 7, 1513-1526, 2013 Author(s): S. MacDonell, C. Kinnard, T. Mölg, L. Nicholson, and J. Abermann Meteorological and surface change measurements collected during a 2.5 yr period are used to calculate surface mass and energy balances at 5324 m a.s.l. on Guanaco Glacier, a cold-based glacier in the semi-arid Andes of Chile. Meteorological conditions are marked by extremely low vapour pressures (annual mean of 1.1 hPa), strong winds (annual mean of 10 m s −1 ), shortwave radiation receipt persistently close to the theoretical site maximum during cloud-free days (mean annual 295 W m −2 ; summer hourly maximum 1354 W m −2 ) and low precipitation rates (mean annual 45 mm w.e.). Snowfall occurs sporadically throughout the year and is related to frontal events in the winter and convective storms during the summer months. Net shortwave radiation provides the greatest source of energy to the glacier surface, and net longwave radiation dominates energy losses. The turbulent latent heat flux is always negative, which means that the surface is always losing mass via sublimation, which is the main form of ablation at the site. Sublimation rates are most strongly correlated with net shortwave radiation, incoming shortwave radiation, albedo and vapour pressure. Low glacier surface temperatures restrict melting for much of the period, however episodic melting occurs during the austral summer, when warm, humid, calm and high pressure conditions restrict sublimation and make more energy available for melting. Low accumulation (131 mm w.e. over the period) and relatively high ablation (1435 mm w.e.) means that mass change over the period was negative (−1304 mm w.e.), which continued the negative trend recorded in the region over the last few decades.

Description: Antarctic ice-mass balance 2003 to 2012: regional reanalysis of GRACE satellite gravimetry measurements with improved estimate of glacial-isostatic adjustment based on GPS uplift rates The Cryosphere, 7, 1499-1512, 2013 Author(s): I. Sasgen, H. Konrad, E. R. Ivins, M. R. Van den Broeke, J. L. Bamber, Z. Martinec, and V. Klemann We present regional-scale mass balances for 25 drainage basins of the Antarctic Ice Sheet (AIS) from satellite observations of the Gravity and Climate Experiment (GRACE) for time period January 2003 to September 2012. Satellite gravimetry estimates of the AIS mass balance are strongly influenced by mass movement in the Earth interior caused by ice advance and retreat during the last glacial cycle. Here, we develop an improved glacial-isostatic adjustment (GIA) estimate for Antarctica using newly available GPS uplift rates, allowing us to more accurately separate GIA-induced trends in the GRACE gravity fields from those caused by current imbalances of the AIS. Our revised GIA estimate is considerably lower than previous predictions, yielding an estimate of apparent mass change of 53 ± 18 Gt yr −1 . Therefore, our AIS mass balance of −114 ± 23 Gt yr −1 is less negative than previous GRACE estimates. The northern Antarctic Peninsula and the Amundsen Sea sector exhibit the largest mass loss (−26 ± 3 Gt yr −1 and −127 ± 7 Gt yr −1 , respectively). In contrast, East Antarctica exhibits a slightly positive mass balance (26 ± 13 Gt yr −1 ), which is, however, mostly the consequence of compensating mass anomalies in Dronning Maud and Enderby Land (positive) and Wilkes and George V Land (negative) due to interannual accumulation variations. In total, 6% of the area constitutes about half the AIS imbalance, contributing 151 ± 7 Gt yr −1 (ca. 0.4 mm yr −1 ) to global mean sea-level change. Most of this imbalance is caused by ice-dynamic speed-up expected to prevail in the near future.

Description: Successful species interactions require that both partners share a similar cue. For many species, spring warming acts as a shared signal to synchronize mutualist behaviors. Spring flowering plants and the ants that disperse their seeds respond to warming temperatures so that ants forage when plants drop seeds. However, where warm-adapted ants replace cold-adapted ants, changes in this timing might leave early seeds stranded without a disperser. We investigate plant seed dispersal south and north of a distinct boundary between warm- and cold-adapted ants to determine if changes in the ant species influence local plant dispersal. The warm-adapted ants forage much later than the cold-adapted ants, and so we first assess natural populations of early and late blooming plants. We then transplant these plants south and north of the ant boundary to test whether distinct ant climate requirements disrupt the ant-plant mutualism. Whereas the early blooming plant's inability to synchronize with the warm-adapted ant leaves its populations clumped and patchy and its seedlings clustered around the parents in natural populations, when transplanted into the range of the cold-adapted ant, effective seed dispersal recovers. In contrast, the mutualism persists for the later blooming plant regardless of location because it sets seed later in spring when both warm- and cold-adapted ant species forage, resulting in effective seed dispersal. These results indicate that the climate response of species interactions, not just the species themselves, is integral in understanding ecological responses to a changing climate. Data linking phenological synchrony and dispersal are rare, and these results suggest a viable mechanism by which a species’ range is limited more by biotic than abiotic interactions – despite the general assumption that biotic influences are buried within larger climate drivers. These results show that biotic partner can be as fundamental a niche requirement as abiotic resources. This article is protected by copyright. All rights reserved.

Description: The effect of training image and secondary data integration with multiple-point geostatistics in groundwater modeling Hydrology and Earth System Sciences Discussions, 10, 11829-11860, 2013 Author(s): X. He, T. O. Sonnenborg, F. Jørgensen, and K. H. Jensen Multiple-point geostatistic simulation (MPS) has recently become popular in stochastic hydrogeology, primarily because of its capability to derive multivariate distributions from the training image (TI). However, its application in three dimensional simulations has been constrained by the difficulty of constructing 3-D TI. The object-based TiGenerator may be a useful tool in this regard; yet the sensitivity of model predictions to the training image has not been documented. Another issue in MPS is the integration of multiple geophysical data. The best way to retrieve and incorporate information from high resolution geophysical data is still under discussion. This work shows that TI from TiGenerator delivers acceptable results when used for groundwater modeling, although the TI directly converted from high resolution geophysical data leads to better simulation. The model results also indicate that soft conditioning in MPS is a convenient and efficient way of integrating secondary data such as 3-D airborne electromagnetic data, but over conditioning has to be avoided.

Description: Multi-decadal river flows variations in France Hydrology and Earth System Sciences Discussions, 10, 11861-11900, 2013 Author(s): J. Boé and F. Habets In this article, multi-decadal variations in French hydroclimate are investigated, with a specific focus on river flows. Based on long observed series, it is shown that river flows in France generally exhibit large multi-decadal variations on the historical period, especially in spring. Differences of means between two 21 yr periods of the 20th century as large as 40% are indeed found for many gauging stations. Multi-decadal spring river flows variations are associated with variations in spring precipitation and temperature. These multi-decadal variations in precipitation are themselves found to be driven by large-scale atmospheric circulation, more precisely by a multi-decadal oscillation in a sea level pressure dipole between western Europe and the East Atlantic. It is suggested that the Atlantic Multidecadal Variability, the main mode of decadal variability in the North Atlantic/Europe sector, controls those variations in large-scale circulation and is therefore the main ultimate driver of multi-decadal variations in spring river flows. Multi-decadal variations in river flows in other seasons, and in particular summer, are also noted. As they are not associated with significant surface climate anomalies (i.e. temperature, precipitation) in summer, other mechanisms are investigated based on hydrological simulations. The impact of climate variations in spring on summer soil moisture, and the impact of soil moisture in summer on the runoff to precipitation ratio, could potentially play a role in multi-decadal summer river flows variations. The large amplitude of the multi-decadal variations in French river flows suggests that internal variability may play a very important role in the evolution of river flows during the next decades, potentially temporarily limiting, reversing or seriously aggravating the long-term impacts of anthropogenic climate change.

Description: Sylwia Kondej In this paper, we study a two-dimensional quantum system governed by the Hamiltonian corresponding to the formal expression Delta(alpha+omega())delta(Gamma), where alpha 〉 0, omega() is a continuous compactly supported function and Gamma is a curve in [openface R]. We derive the asymptotics of eigen ... [J. Math. Phys. 54, 093511 (2013)] published Thu Sep 26, 2013.

Description: Estimating soil suction from electrical resistivity Natural Hazards and Earth System Science, 13, 2369-2379, 2013 Author(s): E. Piegari and R. Di Maio Soil suction and resistivity strongly depend on the degree of soil saturation and, therefore, both are used for estimating water content variations. The main difference between them is that soil suction is measured using tensiometers, which give point information, while resistivity is obtained by tomography surveys, which provide distributions of resistivity values in large volumes, although with less accuracy. In this paper, we have related soil suction to electrical resistivity with the aim of obtaining information about soil suction changes in large volumes, and not only for small areas around soil suction probes. We derived analytical relationships between soil matric suction and electrical resistivity by combining the empirical laws of van Genuchten and Archie. The obtained relationships were used to evaluate maps of soil suction values in different ashy layers originating in the explosive activity of the Mt Somma-Vesuvius volcano (southern Italy). Our findings provided a further example of the high potential of geophysical methods in contributing to more effective monitoring of soil stress conditions; this is of primary importance in areas where rainfall-induced landslides occur periodically.

Description: J. Sadeghi, M. Rostami, and Z. Sadeghi In this paper we consider modified wave equations for spinless particles in an external magnetic field. We consider 4-potentials which guarantee Lorentz' and Coulomb's conditions. The new variable for modified wave equation leads us to consider the associated Laguerre differential equation. We take ... [J. Math. Phys. 54, 093512 (2013)] published Fri Sep 27, 2013.

Description: O. L. de Lange and R. E. Raab Multipole expansions for the macroscopic charge and current densities in a dielectric half-space involve a hierarchy of singular functions comprising the Dirac delta function and its derivatives. For these, Maxwell's differential equations yield corresponding singular expansions of the macroscopic e ... [J. Math. Phys. 54, 093513 (2013)] published Fri Sep 27, 2013.

Description: A data set of world-wide glacier length fluctuations The Cryosphere Discussions, 7, 4775-4811, 2013 Author(s): P. W. Leclercq, J. Oerlemans, H. J. Basagic, I. Bushueva, A. J. Cook, and R. Le Bris Glacier fluctuations contribute to variations in sea level and historical glacier length fluctuations are natural indicators of climate change. To study these subjects, long-term information of glacier change is needed. In this paper we present a~data set of global long-term glacier length fluctuations. The data set is a compilation of available information on changes in glacier length world-wide, including both measured and reconstructed glacier length fluctuations. All 471 length series start before 1950 and cover at least four decades. The longest record starts in 1534, but the majority of time series start after 1850. The number of available records decreases again after 1962. The data set has global coverage including records from all continents. However, the Canadian Arctic is not represented in the data set. The glacier length series show relatively small fluctuations until the mid-19th century followed by a global retreat that was strongest in the first half of the 20th century, although large variability in the length change of the different glaciers is observed. During the 20th century, calving glaciers retreated more than land terminating glaciers, but their relative length change was approximately equal. Besides calving, the glacier slope is the most important glacier property determining length change: steep glaciers have retreated less than glaciers with a gentle slope.

Description: Increasing ocean temperatures and strengthening boundary currents have caused the poleward migration of many marine species. Cubozoan jellyfish known to cause Irukandji syndrome have historically been confined to tropical waters but may be expanding into sub-tropical regions. Here we examine the interactive effects of warming and acidification on the population dynamics of polyps of an Irukandji jellyfish, Alatina nr mordens, and the formation of statoliths in newly metamorphosed medusae, to determine if this jellyfish could tolerate future conditions predicted for southeast Queensland (SEQ), Australia. Two experiments, examining the orthogonal factors of temperature and pH were undertaken. Experiment 1 mimicked the current, ca. 2050 and ca. 2100 summer temperature and pH conditions predicted for SEQ using A1F1 scenarios (temperature: 25, 27, 29°C; pH: 7.9, 7.8, 7.6) and Experiment 2 mimicked current and future winter conditions (18 and 22°C, pH 7.9, 7.8, 7.6). All polyps in Experiment 1 survived and budded. Fewer polyps budded in the lower pH treatments but patterns varied slightly among temperature treatments. Statoliths at pH 7.6 were 24% narrower than those at pH 7.8 and 7.9. Most polyps survived the winter conditions mimicked by Experiment 2 but only polyps in the 22°C, pH 7.9 treatment increased significantly. The current absence of A . nr mordens medusae in SEQ, despite the polyps’ ability to tolerate the current temperature and pH conditions, suggests that ecological, rather than abiotic factors currently limit their distribution. Observations that budding was lower under low pH treatments suggest that rates of asexual reproduction will likely be much slower in the future. We consider that A . nr mordens polyps are likely to tolerate future conditions but are unlikely to thrive in the long term. However, if polyps can overcome potential ecological boundaries and acidification proceeds slowly A . nr mordens could expand polewards in the short-term. This article is protected by copyright. All rights reserved.

Description: Soil CO 2 efflux ( F soil ) is the largest source of carbon from forests and reflects primary productivity as well as how carbon is allocated within forest ecosystems. Through early stages of stand development, both elevated [CO 2 ] and availability of soil nitrogen (N; sum of mineralization, deposition, and fixation) have been shown to increase gross primary productivity, but the long-term effects of these factors on F soil are less clear. Expanding on previous studies at the Duke Free Air CO 2 Enrichment (FACE) site, we quantified the effects of elevated [CO 2 ] and N fertilization on F soil using daily measurements from automated chambers over 10 years. Consistent with previous results, compared to ambient-unfertilized plots, annual F soil increased under elevated [CO 2 ] (~17%) and decreased with N (~21%). N fertilization under elevated [CO 2 ] reduced F soil to values similar to untreated plots. Over the study period, base respiration rates increased with leaf productivity but declined after productivity saturated. Despite treatment-induced differences in aboveground biomass, soil temperature and water content were similar among treatments. Inter-annually, low soil water content decreased annual F soil from potential values – estimated based on temperature alone assuming non-limiting soil water content – by ~0.7% per 1.0% reduction in relative extractable water. This effect was only slightly ameliorated by elevated [CO 2 ]. Variability of soil N availability among plots accounted for the spatial variability of F soil , showing a decrease of ~114 g C m -2 y -1 per 1 g m -2 increase in soil N availability, with consistently higher F soil in elevated [CO 2 ] plots ~127 g C per 100 ppm [CO 2 ] over the +200 ppm enrichment. Altogether, reflecting increased belowground carbon partitioning in response to greater plant nutritional needs, the effects of elevated [CO 2 ] and N fertilization on F soil in this stand are sustained beyond the early stages of stand development and through stabilization of annual foliage production. This article is protected by copyright. All rights reserved.

Description: To meet growing global food demand with limited land and reduced environmental impact, agricultural greenhouse gas (GHG) emissions are increasingly evaluated with respect to crop productivity, i.e. on a yield-scaled as opposed to area basis. Here, we compiled available field data on CH 4 and N 2 O emissions from rice production systems to test the hypothesis that in response to fertilizer nitrogen (N) addition, yield-scaled global warming potential (GWP) will be minimized at N rates that maximize yields. Within each study, yield N surplus was calculated to estimate deficit or excess N application rates with respect to the optimal N rate (defined as the N rate at which maximum yield was achieved). Relationships between yield N surplus and GHG emissions were assessed using linear and nonlinear mixed-effects models. Results indicate that yields increased moving from deficit to optimal N rates. At N rates contributing to a yield N surplus, N 2 O and yield-scaled N 2 O emissions increased exponentially. In contrast, CH 4 emissions were not impacted by N inputs. Accordingly, yield-scaled CH 4 emissions decreased with N addition. Overall, yield-scaled GWP was minimized at optimal N rates, decreasing by 21% compared to treatments without N addition. These results are unique compared to aerobic cropping systems in which N 2 O emissions are the primary contributor to GWP, meaning yield-scaled GWP may not necessarily decrease for aerobic crops when yields are optimized by N fertilizer additions. Balancing gains in agricultural productivity with climate change concerns, this work supports the concept that high rice yields can be achieved with minimal yield-scaled GWP through optimal N application rates. Moreover, additional improvements in N use efficiency may further reduce yield-scaled GWP, thereby strengthening the economic and environmental sustainability of rice systems. This article is protected by copyright. All rights reserved.

Description: Climate change is projected to push the limits of cropping systems and has the potential to disrupt the agricultural sector from local to global scales. This article introduces the Coordinated Climate-Crop Modeling Project (C3MP), an initiative of the Agricultural Model Intercomparison and Improvement Project (AgMIP) to engage a global network of crop modelers to explore the impacts of climate change via an investigation of crop responses to changes in carbon dioxide concentration ([CO 2 ]), temperature, and water. As a demonstration of the C3MP protocols and enabled analyses, we apply the Decision Support System for Agrotechnology Transfer (DSSAT) CROPGRO-Peanut crop model for Henry County, Alabama, to evaluate responses to the range of plausible [CO 2 ], temperature changes, and precipitation changes projected by climate models out to the end of the 21 st century. These sensitivity tests are used to derive crop model emulators that estimate changes in mean yield and the coefficient of variation for seasonal yields across a broad range of climate conditions, reproducing mean yields from sensitivity test simulations with deviations of ~2% for rainfed conditions. We apply these statistical emulators to investigate how peanuts respond to projections from various global climate models, time periods, and emissions scenarios, finding a robust projection of modest (〈10%) median yield losses in the middle of the 21 st century accelerating to more severe (〉20%) losses and larger uncertainty at the end of the century under the more severe representative concentration pathway 8.5. This projection is not substantially altered by the selection of the AgMERRA global gridded climate dataset rather than the local historical observations, differences between the Third and Fifth Coupled Model Intercomparison Project (CMIP3 and CMIP5), or the use of the delta method of climate impacts analysis rather than the C3MP impacts response surface and emulator approach. This article is protected by copyright. All rights reserved.

Description: Large-scale, long-term FACE (Free Air CO 2 -enrichment) experiments indicate that increases in atmospheric CO 2 concentrations will influence forest C cycling in unpredictable ways. It has been recently suggested that responses of mycorrhizal fungi could determine whether forest NPP (net primary production) is increased by elevated CO 2 over long time periods and if forests soils will function as sources or sinks of C in the future. We studied the dynamic responses of ectomycorrhizae to N fertilization and atmospheric CO 2 -enrichment at the Duke FACE experiment using minirhizotrons over a six year period (2005-2010). Stimulation of mycorrhizal production by elevated CO 2 was observed during only one (2007) of six years. This increased the standing crop of mycorrhizal tips during 2007 and 2008; during 2008, significantly higher mortality returned standing crop to ambient levels for the remainder of the experiment. It is therefore unlikely that increased production of mycorrhizal tips can explain the lack of progressive nitrogen limitations and associated increases in N uptake observed in CO 2 -enriched plots at this site. Fertilization generally decreased tree reliance on mycorrhizae as tip production declined with the addition of nitrogen as has been shown in many other studies. Annual NPP of mycorrhizal tips was greatest during years with warm January temperatures and during years with cool spring temperatures. A 2° C increase in average late spring temperatures (May and June) decreased annual production of mycorrhizal root tip length by 50%. This has important implications for ecosystem function in a warmer world in addition to potential for forest soils to sequester atmospheric C. This article is protected by copyright. All rights reserved.

Description: An original interpretation of the wet edge of the surface temperature–albedo space to estimate crop evapotranspiration (SEB-1S), and its validation over an irrigated area in northwestern Mexico Hydrology and Earth System Sciences, 17, 3623-3637, 2013 Author(s): O. Merlin The space defined by the pair surface temperature ( T ) and surface albedo (α), and the space defined by the pair T and fractional green vegetation cover ( f vg ) have been extensively used to estimate evaporative fraction (EF) from solar/thermal remote sensing data. In both space-based approaches, evapotranspiration (ET) is estimated as remotely sensed EF times the available energy. For a given data point in the T -α space or in the T - f vg space, EF is derived as the ratio of the distance separating the point from the line identified as the dry edge to the distance separating the dry edge and the line identified as the wet edge. The dry and wet edges are classically defined as the upper and lower limit of the spaces, respectively. When investigating side by side the T -α and the T - f vg spaces, one observes that the range covered by T values on the (classically determined) wet edge is different for both spaces. In addition, when extending the wet and dry lines of the T -α space, both lines cross at α ≈ 0.4 although the wet and dry edges of the T - f vg space never cross for 0 ≤ f vg 〈 1. In this paper, a new ET (EF) model (SEB-1S) is derived by revisiting the classical physical interpretation of the T -α space to make its wet edge consistent with that of the T - f vg space. SEB-1S is tested over a 16 km by 10 km irrigated area in northwestern Mexico during the 2007–2008 agricultural season. The classical T -α space-based model is implemented as benchmark to evaluate the performance of SEB-1S. Input data are composed of ASTER (Advanced Spaceborne Thermal Emission and Reflection radiometer) thermal infrared, Formosat-2 shortwave, and station-based meteorological data. The fluxes simulated by SEB-1S and the classical T -α space-based model are compared on seven ASTER overpass dates with the in situ measurements collected at six locations within the study domain. The ET simulated by SEB-1S is significantly more accurate and robust than that predicted by the classical T -α space-based model. The correlation coefficient and slope of the linear regression between simulated and observed ET is improved from 0.82 to 0.93, and from 0.63 to 0.90, respectively. Moreover, constraining the wet edge using air temperature data improves the slope of the linear regression between simulated and observed ET.

Description: Teleconnection analysis of runoff and soil moisture over the Pearl River basin in South China Hydrology and Earth System Sciences Discussions, 10, 11943-11982, 2013 Author(s): J. Niu, J. Chen, and B. Sivakumar This study explores the teleconnection of two climatic patterns, namely the El Niño-Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD), with hydrological processes over the Pearl River basin in South China. The Variable Infiltration Capacity (VIC) model is used to simulate the daily hydrological processes over the basin for the study period 1952–2000, and then, using the simulation results, the time series of the monthly runoff and soil moisture anomalies for its ten sub-basins are aggregated. Wavelet analysis is performed to explore the variability properties of these time series at 49 timescales ranging from 2 months to 9 yr. Use of wavelet coherence and rank correlation method reveals that the dominant variabilities of the time series of runoff and soil moisture are basically correlated with IOD. The influences of ENSO on the terrestrial hydrological processes are mainly found in the eastern sub-basins. The teleconnections between climatic patterns and hydrological variability also serve as a reference basis for inferences on the occurrence of extreme hydrological events (e.g. floods and droughts).

Description: On an improved sub-regional water resources management representation for integration into earth system models Hydrology and Earth System Sciences, 17, 3605-3622, 2013 Author(s): N. Voisin, H. Li, D. Ward, M. Huang, M. Wigmosta, and L. R. Leung Human influence on the hydrologic cycle includes regulation and storage, consumptive use and overall redistribution of water resources in space and time. Representing these processes is essential for applications of earth system models in hydrologic and climate predictions, as well as impact studies at regional to global scales. Emerging large-scale research reservoir models use generic operating rules that are flexible for coupling with earth system models. Those generic operating rules have been successful in reproducing the overall regulated flow at large basin scales. This study investigates the uncertainties of the reservoir models from different implementations of the generic operating rules using the complex multi-objective Columbia River Regulation System in northwestern United States as an example to understand their effects on not only regulated flow but also reservoir storage and fraction of the demand that is met. Numerical experiments are designed to test new generic operating rules that combine storage and releases targets for multi-purpose reservoirs and to compare the use of reservoir usage priorities and predictors (withdrawals vs. consumptive demands, as well as natural vs. regulated mean flow) for configuring operating rules. Overall the best performing implementation is with combined priorities rules (flood control storage targets and irrigation release targets) set up with mean annual natural flow and mean monthly withdrawals. The options of not accounting for groundwater withdrawals, or on the contrary, of assuming that all remaining demand is met through groundwater extractions, are discussed.

Description: Permafrost thaw in the Arctic driven by climate change is mobilizing ancient terrigenous organic carbon (OC) into fluvial networks. Understanding the controls on metabolism of this OC is imperative for assessing its role with respect to climate feedbacks. In this study we examined the effect of inorganic nutrient supply and dissolved organic matter (DOM) composition on aquatic extracellular enzyme activities (EEAs) in waters draining the Kolyma River Basin (Siberia), including permafrost derived OC. Reducing the phenolic content of the DOM pool resulted in dramatic increases in hydrolase EEAs (e.g. phosphatase activity increased 〉 28 fold) supporting the idea that high concentrations of polyphenolic compounds in DOM (e.g. plant structural tissues) inhibit enzyme synthesis or activity, limiting OC degradation. EEAs were significantly more responsive to inorganic nutrient additions only after phenolic inhibition was experimentally removed. In controlled mixtures of modern OC and thawed permafrost endmember OC sources, respiration rates per unit dissolved OC were 1.3 – 1.6 times higher in waters containing ancient carbon, suggesting that permafrost derived OC was more available for microbial mineralization. In addition, waters containing ancient permafrost derived OC supported elevated phosphatase and glucosidase activities. Based on these combined results, we propose that both composition and nutrient availability regulates DOM metabolism in Arctic aquatic ecosystems. Our empirical findings are incorporated into a mechanistic conceptual model highlighting two key enzymatic processes in the mineralization of riverine OM: 1) the role of phenol oxidase activity in reducing inhibitory phenolic compounds; and 2) the role of phosphatase in mobilizing organic P. Permafrost derived DOM degradation was less constrained by this initial “phenolic-OM” inhibition; thus, informing reports of high biological availability of ancient, permafrost derived DOM with clear ramifications for its metabolism in fluvial networks and feedbacks to climate. This article is protected by copyright. All rights reserved.

Description: On selection of the optimal data time interval for real-time hydrological forecasting Hydrology and Earth System Sciences, 17, 3639-3659, 2013 Author(s): J. Liu and D. Han With the advancement in modern telemetry and communication technologies, hydrological data can be collected with an increasingly higher sampling rate. An important issue deserving attention from the hydrological community is which suitable time interval of the model input data should be chosen in hydrological forecasting. Such a problem has long been recognised in the control engineering community but is a largely ignored topic in operational applications of hydrological forecasting. In this study, the intrinsic properties of rainfall–runoff data with different time intervals are first investigated from the perspectives of the sampling theorem and the information loss using the discrete wavelet transform tool. It is found that rainfall signals with very high sampling rates may not always improve the accuracy of rainfall–runoff modelling due to the catchment low-pass-filtering effect. To further investigate the impact of a data time interval in real-time forecasting, a real-time forecasting system is constructed by incorporating the probability distributed model (PDM) with a real-time updating scheme, the autoregressive moving-average (ARMA) model. Case studies are then carried out on four UK catchments with different concentration times for real-time flow forecasting using data with different time intervals of 15, 30, 45, 60, 90 and 120 min. A positive relation is found between the forecast lead time and the optimal choice of the data time interval, which is also highly dependent on the catchment concentration time. Finally, based on the conclusions from the case studies, a hypothetical pattern is proposed in three-dimensional coordinates to describe the general impact of the data time interval and to provide implications of the selection of the optimal time interval in real-time hydrological forecasting. Although nowadays most operational hydrological systems still have low data sampling rates (daily or hourly), the future is that higher sampling rates will become more widespread, and there is an urgent need for hydrologists both in academia and in the field to realise the significance of the data time interval issue. It is important that more case studies in different catchments with various hydrological forecasting systems are explored in the future to further verify and improve the proposed hypothetical pattern.

Description: Due to the intrinsic complexities associated with modeling land-atmosphere interactions, global models typically use elaborate land surface and boundary layer physics parameterizations. Unfortunately, it is difficult to use elaborate models, by themselves, to develop a deeper understanding of how land surface parameters affect the coupled land-atmosphere system. At the same time, it is also increasingly important to gain a deeper understanding of the role of changes in land cover, land use, and ecosystem function as forcings and feedbacks in past and future climate change. To improve the foundation of our understanding, we outline a framework for boundary layer climate sensitivity based on surface energy balance; just as global climate sensitivity is based on top-of-atmosphere energy balance. We develop an analytic theory for the boundary layer climate sensitivity of an idealized model of a diurnally-averaged well-mixed boundary layer over land [ Betts, 2000]. This analytic sensitivity theory identifies changes in the properties of the land surface — including moisture availability, albedo, and aerodynamic roughness — as forcings, and identifies strong negative feedbacks associated with the surface fluxes of latent and sensible heat. We show that our theory can explain nearly all of the sensitivity of the Betts [2000] full system of equations. Favorable comparison of the theory and the simulation results from a two-column radiative convective model suggests that the theory may be broadly useful for unifying our understanding of how changes in land use or ecosystem function may affect climate change.

Description: Ecosystem functioning is simultaneously affected by changes in community composition and environmental change such as increasing atmospheric carbon dioxide (CO 2 ) and subsequent ocean acidification. However, it largely remains uncertain how the effects of these factors compare to each other. Addressing this question, we experimentally tested the hypothesis that initial community composition and elevated CO 2 are equally important to the regulation of phytoplankton biomass. We full-factorially exposed three compositionally different marine phytoplankton communities to two different CO 2 levels and examined the effects and relative importance (ω 2 ) of the two factors and their interaction on phytoplankton biomass at bloom peak. The results showed that initial community composition had a significantly greater impact than elevated CO 2 on phytoplankton biomass, which varied largely among communities. We suggest that the different initial ratios between cyanobacteria, diatoms, and dinoflagellates might be the key for the varying competitive and thus functional outcome among communities. Furthermore, the results showed that depending on initial community composition elevated CO 2 selected for larger sized diatoms, which led to increased total phytoplankton biomass. Our study highlights the relevance of initial community composition, which strongly drives the functional outcome, when assessing impacts of climate change on ecosystem functioning. In particular, the increase in phytoplankton biomass driven by the gain of larger sized diatoms in response to elevated CO 2 potentially has strong implications for nutrient cycling and carbon export in future oceans. This article is protected by copyright. All rights reserved.

Description: Less than half of anthropogenic carbon emissions are accumulating in the atmosphere, due to large net fluxes into both the oceans and the land (Le Queré et al., 2012). The land sink in particular has increased markedly, doubling in strength since the 1960's, to reach 26 petagrams of carbon in the latest decade. However, the location and drivers of this large terrestrial sink are still relatively poorly constrained by atmospheric measurements (Ciais et al. 2013). Pan et al. (2011) recently utilised 〉1 million forest inventory plots to provide summaries of forest carbon stocks, and the first global bottom-up estimates of carbon fluxes for the world's forest biomes for the period 1990-2007. One key result was that almost all the residual global terrestrial carbon sink (i.e. carbon uptake after accounting for land use change), some 2.4 ± 0.4 Pg of carbon per year, is located in the world's established forests (Pan et al., 2011). The sink is distributed worldwide, with globally significant net fluxes into boreal and temperate forests, and a large sink in intact tropical forest, albeit with large uncertainty. Furthermore, Pan et al. (2011) showed that this tropical intact forest sink - may have faded from an estimated annual 1.3 ± 0.4 Pg C in the 1990's to 1.0 ± 0.5 Pg C for 2000-2007. The tropical intact sink is offset by a net land-use emission (1.5 Pg C yr −1 [1990-1999]) declining to 1.1 Pg C yr −1 [2000-2007]), and as a consequence aircraft measurements and inverse modelling studies indicate the tropics to be close to neutral in terms of net carbon fluxes (reviewed by Ciais et al. 2013). While the intact tropical forest sink values represent updates from similar values published previously (e.g., Lewis et al., 2009a), the fact that almost the entire residual terrestrial carbon sink is accounted for by the forests of the world was a notable discovery. Evidence from the ground now points to established forests being a net sink across almost every major forest region, including all extra-tropical forest regions analysed. This article is protected by copyright. All rights reserved.

Description: Forecasting how global warming will affect onset of the growing season is essential for predicting terrestrial productivity, but suffers from conflicting evidence. We show that accurate estimates require ways to connect discrete observations of changing tree status (e.g., pre- vs. post-budbreak) with continuous responses to fluctuating temperatures. By coherently synthesizing discrete observations with continuous responses to temperature variation, we accurately quantify how increasing temperature variation accelerates onset of growth. Application to warming experiments at two latitudes demonstrates that maximum responses to warming are concentrated in late winter, weeks ahead of the main budbreak period. Given that warming will not occur uniformly over the year, knowledge of when temperature variation has the most impact can guide prediction. Responses are large and heterogeneous, yet predictable. The approach has immediate application to forecasting effects of warming on growing season length, requiring only information that is readily available from weather stations and generated in climate models. This article is protected by copyright. All rights reserved.

Description: Climate change is expected to cause geographic redistributions of species. To the extent that species within assemblages have different niche requirements, assemblages may no longer remain intact and dis- and reassemble at current or new geographic locations. We explored how climate change projected by 2100 may transform the world's avian assemblages (characterized at a 110 km spatial grain) by modelling environmental niche-based changes to their dietary guild structure under 0 km, 500 km, and 2000 km dispersal distances. We examined guild structure changes at coarse (primary, high-level, and mixed consumers) and fine (frugivores, nectarivores, insectivores, herbivores, granivores, scavengers, omnivores, and carnivores) ecological resolutions to determine whether or not geographic co-occurrence patterns among guilds were associated with the magnitude to which guilds are functionally resolved. Dietary guilds vary considerably in their global geographic prevalence, and under broad-scale niche-based redistribution of species, these are projected to change very heterogeneously. A non-dispersal assumption results in the smallest projected changes to guild assemblages, but with significant losses for some regions and guilds, such as South American insectivores. Longer dispersal distances are projected to cause greater degrees of disassembly, and lead to greater homogenization of guild composition, especially in northern Asia and Africa. This arises because projected range gains and losses result in geographically heterogeneous patterns of guild compensation. Projected decreases especially of primary and mixed consumers most often are compensated by increases in high-level consumers, with increasing uncertainty about these outcomes as dispersal distance and degree of guild functional resolution increases. Further exploration into the consequences of these significant broad-scale ecological functional changes at the community or ecosystem level should be increasingly on the agenda for conservation science. This article is protected by copyright. All rights reserved.

Description: Climate change scenarios predict increases in the frequency and duration of ENSO-related droughts for parts of South-East Asia until the end of this century exposing the remaining rainforests to increasing drought risk. A pan-tropical review of recorded drought-related tree mortalities in more than 100 monitoring plots before, during and after drought events suggested a higher drought-vulnerability of trees in South-East Asian than in Amazonian forests. Here, we present the results of a replicated (n=3 plots) throughfall exclusion experiment in a perhumid tropical rainforest in Sulawesi, Indonesia. In this first large-scale roof experiment outside semi-humid eastern Amazonia, 60% of the throughfall was displaced during the first 8 months and 80% during the subsequent 17 months, exposing the forest to severe soil desiccation for about 17 months. In the experiment's second year, wood production decreased on average by 40% with largely different responses of the tree families (ranging from -100 to +100% change). Most sensitive were trees with high radial growth rates under moist conditions. In contrast, tree height was only a secondary factor and wood specific gravity had no influence on growth sensitivity. Fine root biomass was reduced by 35% after 25 months of soil desiccation while fine root necromass increased by 250% indicating elevated fine root mortality. Cumulative aboveground litter production was not significantly reduced in this period. The trees from this Indonesian perhumid rainforest revealed similar responses of wood and litter production and root dynamics as those in two semi-humid Amazonian forests subjected to experimental drought. We conclude that trees from paleo- or neotropical forests growing in semi-humid or perhumid climates may not differ systematically in their growth sensitivity and vitality under sub-lethal drought stress. Drought vulnerability may depend more on stem cambial activity in moist periods than on tree height or wood specific gravity. This article is protected by copyright. All rights reserved.

Description: We examined the water balance a forested ombrotrophic peatland and adjacent burned peatland in the boreal plain of western Canada over a three-year period. Complete combustion of foliage and fine branches dramatically increased shortwave radiation inputs to the peat surface while halting all tree transpiration at the burned site. End-of-winter snowpack was 7-25% higher at the burned site likely due to decreased ablation from the tree canopy at the unburned site. Shrub regrowth at the burned site was rapid post-fire, and shading by the shrub canopy in the burned site approached that of the unburned site within three years after fire. Site-averaged surface resistance to evaporation was not different between sites, though surface resistance in hollows was lower in the burned site. Water loss at both burned and unburned sites is largely driven by surface evaporative losses. Evaporation at the burned site marginally exceeded the sum of pre-fire transpiration and interception at the unburned site, suggesting that ET during the growing season was 2 0–40  mm greater at the the burned peatland. While the net change in water storage during the growing season was largely unchanged by fire, the lack of low-density surface peat in the burned site appears to have decreased specific yield, leading to greater water table decline at the burned site despite similar net change in storage. This article is protected by copyright. All rights reserved.

Description: Soil microbial communities in Chihuahuan Desert grasslands generally experience highly variable spatiotemporal rainfall patterns. Changes in precipitation regimes can affect belowground ecosystem processes such as decomposition and nutrient cycling by altering soil microbial community structure and function. The objective of this study was to determine if increased seasonal precipitation frequency and magnitude over a seven-year period would generate a persistent shift in microbial community characteristics and soil nutrient availability. We supplemented natural rainfall with large events (one/winter and three/summer) to simulate increased precipitation based on climate model predictions for this region. We observed a two year delay in microbial responses to supplemental precipitation treatments. In Years 3-5, higher microbial biomass, arbuscular mycorrhizae abundance, and soil enzyme C and P acquisition activities were observed in the supplemental water plots even during extended drought periods. In Years 5-7, available soil P was consistently lower in the watered plots compared to control plots. Shifts in soil P corresponded to higher fungal abundances, microbial C utilization activity, and soil pH. This study demonstrated that 25% shifts in seasonal rainfall can significantly influence soil microbial and nutrient properties, which in turn may have long-term effects on nutrient cycling and plant P uptake in this desert grassland. This article is protected by copyright. All rights reserved.

Description: It is proposed that increases in anthropogenic reactive nitrogen (N r )-deposition may cause temperate and boreal forests to sequester a globally significant quantity of carbon (C); however, long-term data from boreal forests describing how C sequestration responds to realistic levels of chronic N r -deposition are scarce. Using a long term (14-year) stand scale (0.1 ha) N-addition experiment (three levels: 0, 12.5, and 50 kg N ha −1 yr −1 ) in the boreal zone of northern Sweden, we evaluated how chronic N additions altered N uptake and biomass of understory communities, and whether changes in understory communities explained N uptake and C sequestration by trees. We hypothesized that understory communities (i.e. mosses and shrubs) serve as important sinks for low-level N additions, with the strength of these sinks weakening as chronic N addition rates increase, due to shifts in species composition. We further hypothesized that trees would exhibit non-linear increases in N acquisition, and subsequent C sequestration as N addition rates increased, due to a weakening understory N sink. Our data showed that understory biomass was reduced by 50% in response to the high N addition treatment, mainly due to reduced moss biomass. A 15 N labelling experiment showed that feather mosses acquired the largest fraction of applied label, with this fraction decreasing as the chronic N addition level increased. Contrary to our hypothesis, the proportion of label taken up by trees was equal (~8%) across all three N addition treatments. The relationship between N addition and C sequestration in all vegetation pools combined was linear, and had a slope of 16 kg C kg −1 N. While canopy retention of N r deposition may cause C sequestration rates to be slightly different than this estimate, our data suggests that a minor quantity of annual anthropogenic CO 2 emissions are sequestered into boreal forests as a result of N r deposition. This article is protected by copyright. All rights reserved.

Description: Benchmark products for land evapotranspiration: LandFlux-EVAL multi-data set synthesis Hydrology and Earth System Sciences, 17, 3707-3720, 2013 Author(s): B. Mueller, M. Hirschi, C. Jimenez, P. Ciais, P. A. Dirmeyer, A. J. Dolman, J. B. Fisher, M. Jung, F. Ludwig, F. Maignan, D. G. Miralles, M. F. McCabe, M. Reichstein, J. Sheffield, K. Wang, E. F. Wood, Y. Zhang, and S. I. Seneviratne Land evapotranspiration (ET) estimates are available from several global data sets. Here, monthly global land ET synthesis products, merged from these individual data sets over the time periods 1989–1995 (7 yr) and 1989–2005 (17 yr), are presented. The merged synthesis products over the shorter period are based on a total of 40 distinct data sets while those over the longer period are based on a total of 14 data sets. In the individual data sets, ET is derived from satellite and/or in situ observations (diagnostic data sets) or calculated via land-surface models (LSMs) driven with observations-based forcing or output from atmospheric reanalyses. Statistics for four merged synthesis products are provided, one including all data sets and three including only data sets from one category each (diagnostic, LSMs, and reanalyses). The multi-annual variations of ET in the merged synthesis products display realistic responses. They are also consistent with previous findings of a global increase in ET between 1989 and 1997 (0.13 mm yr −2 in our merged product) followed by a significant decrease in this trend (−0.18 mm yr −2 ), although these trends are relatively small compared to the uncertainty of absolute ET values. The global mean ET from the merged synthesis products (based on all data sets) is 493 mm yr −1 (1.35 mm d −1 ) for both the 1989–1995 and 1989–2005 products, which is relatively low compared to previously published estimates. We estimate global runoff (precipitation minus ET) to 263 mm yr −1 (34 406 km 3 yr −1 ) for a total land area of 130 922 000 km 2 . Precipitation, being an important driving factor and input to most simulated ET data sets, presents uncertainties between single data sets as large as those in the ET estimates. In order to reduce uncertainties in current ET products, improving the accuracy of the input variables, especially precipitation, as well as the parameterizations of ET, are crucial.

Description: Urbanization is one of the most extensive and ecologically significant changes happening to terrestrial environments, as it strongly affects the distribution of biodiversity. It is well established that native species richness is reduced in urban and suburban areas, but the species traits that predict tolerance to urbanization are yet little understood. In birds, one of the most studied groups in this respect, evidence is appearing that acoustic traits influence urban living, but it is unknown how this compares to the effects of more obvious ecological traits that facilitate urban living. Therefore, it remains unclear whether acoustic communication is an important predictor of urban tolerance among species. Here, with a comparative study across 140 European and North American passerines, I show that high song frequency, which is less masked by the low-frequency anthropogenic noise, is associated with urban tolerance, with an effect size over half that of the most important ecological trait studied: off-ground nesting. Other nesting and foraging traits accepted to facilitate urban living did not differ for species occurring in urban environments. Thus, the contribution of acoustic traits for passerine urban tolerance approximates that of more obvious ecological traits. Nonetheless, effect sizes of the biological predictors of urban tolerance were low and the phylogenetic signal for urban tolerance was null, both of which suggest that factors other than phenotypic traits have major effects on urban tolerance. A simple possibility is exposure to urbanization, since there was a higher proportion of urban-tolerant species in Europe, which is more urbanized than North America. This article is protected by copyright. All rights reserved.

Description: Spatial patterns in timing of the diurnal temperature cycle Hydrology and Earth System Sciences, 17, 3695-3706, 2013 Author(s): T. R. H. Holmes, W. T. Crow, and C. Hain This paper investigates the structural difference in timing of the diurnal temperature cycle (DTC) over land resulting from choice of measuring device or model framework. It is shown that the timing can be reliably estimated from temporally sparse observations acquired from a constellation of low Earth-orbiting satellites given record lengths of at least three months. Based on a year of data, the spatial patterns of mean DTC timing are compared between temperature estimates from microwave Ka-band, geostationary thermal infrared (TIR), and numerical weather prediction model output from the Global Modeling and Assimilation Office (GMAO). It is found that the spatial patterns can be explained by vegetation effects, sensing depth differences and more speculatively the orientation of orographic relief features. In absolute terms, the GMAO model puts the peak of the DTC on average at 12:50 local solar time, 23 min before TIR with a peak temperature at 13:13 (both averaged over Africa and Europe). Since TIR is the shallowest observation of the land surface, this small difference represents a structural error that possibly affects the model's ability to assimilate observations that are closely tied to the DTC. The equivalent average timing for Ka-band is 13:44, which is influenced by the effect of increased sensing depth in desert areas. For non-desert areas, the Ka-band observations lag the TIR observations by only 15 min, which is in agreement with their respective theoretical sensing depth. The results of this comparison provide insights into the structural differences between temperature measurements and models, and can be used as a first step to account for these differences in a coherent way.

Description: The eddy covariance (EC) method was used in a 30-month study to quantify evapotranspiration (ET) and vegetation coefficient (K CW ) for a wetland on a ranch in subtropical south Florida. To evaluate the errors in ET estimates, the EC-based ET (ET C.-EC ) and the Food and Agricultural Organization (FAO) Penman-Monteith (PM) based ET (ET C.-PM ) estimates (with literature crop coefficient K C ) were compared to each other. The ET C.-EC and FAO-PM reference ET were used to develop K CW . Regression models were developed to estimate K CW using climatic and hydrologic variables. Annual and daily ET C.-EC values were 1152 mm and 3.27 mm, respectively. The FAO-PM model underestimated ET by 25% with ET C.-EC being statistically higher than ET C.-PM . The K CW varied from 0.79 (December) to 1.06 (November). The mean K CW for dry (November-April) season (0.95) was much higher than values reported for wetlands in literature, while wet (May-October) season K CW (0.97) was closer to literature values. Higher than expected K CW values during dry season were due to higher temperature, lower humidity and perennial wetland vegetation. Regression analyses showed that factors affecting the K CW were different during the dry (soil moisture, temperature, and relative humidity) and wet (net radiation, inundation, and wind speed) seasons. Separate regression models for the dry and wet seasons were developed. ET and K CW from this study, one of the first for the agricultural wetlands in subtropical environment, will help improve the ET estimates for similar wetlands. This article is protected by copyright. All rights reserved.

Description: J. Bakosi and J. R. Ristorcelli The method of potential solutions of Fokker-Planck equations is used to develop a transport equation for the joint probability of N stochastic variables with Lochner's generalized Dirichlet distribution [R. H. Lochner, A generalized Dirichlet distribution in Bayesian life testing, J. R. Stat. Soc. S ... [J. Math. Phys. 54, 102701 (2013)] published Tue Oct 1, 2013.

Description: M. V. Altaisky and N. E. Kaputkina Nondispersive wave packets are widely used in optics and acoustics. We found it interesting that such packets could be also a subject of quantum field theory. Canonical commutation relations for the nondispersive wave packets are constructed. ... [J. Math. Phys. 54, 102101 (2013)] published Tue Oct 1, 2013.