Budgeting rockfall and modeling sedimentary delivery in torrent systems

This thesis is a compilation of projects to study sediment processes recharging debris flow channels. These works, conducted during my stay at the University of Lausanne, focus in the geological and morphological implications of torrent catchments to characterize debris supply, a fundamental element to predict debris flows. Other aspects of sediment dynamics are considered, e.g. the coupling headwaters - torrent, as well as the development of a modeling software that simulates sediment transfer in torrent systems. The sediment activity at Manival, an active torrent system of the northern French Alps, was investigated using terrestrial laser scanning and supplemented with geostructural investigations and a survey of sediment transferred in the main torrent. A full year of sediment flux could be observed, which coincided with two debris flows and several bedload transport events. This study revealed that both debris flows generated in the torrent and were preceded in time by recharge of material from the headwaters. Debris production occurred mostly during winter - early spring time and was caused by large slope failures. Sediment transfers were more puzzling, occurring almost exclusively in early spring subordinated to runoffconditions and in autumn during long rainfall. Intense rainstorms in summer did not affect debris storage that seems to rely on the stability of debris deposits. The morpho-geological implication in debris supply was evaluated using DEM and field surveys. A slope angle-based classification of topography could characterize the mode of debris production and transfer. A slope stability analysis derived from the structures in rock mass could assess susceptibility to failure. The modeled rockfall source areas included more than 97% of the recorded events and the sediment budgets appeared to be correlated to the density of potential slope failure. This work showed that the analysis of process-related terrain morphology and of susceptibility to slope failure document the sediment dynamics to quantitatively assess erosion zones leading to debris flow activity. The development of erosional landforms was evaluated by analyzing their geometry with the orientations of potential rock slope failure and with the direction of the maximum joint frequency. Structure in rock mass, but in particular wedge failure and the dominant discontinuities, appear as a first-order control of erosional mechanisms affecting bedrock- dominated catchment. They represent some weaknesses that are exploited primarily by mass wasting processes and erosion, promoting not only the initiation of rock couloirs and gullies, but also their propagation. Incorporating the geological control in geomorphic processes contributes to better understand the landscape evolution of active catchments. A sediment flux algorithm was implemented in a sediment cascade model that discretizes the torrent catchment in channel reaches and individual process-response systems. Each conceptual element includes in simple manner geomorphological and sediment flux information derived from GIS complemented with field mapping. This tool enables to simulate sediment transfers in channels considering evolving debris supply and conveyance, and helps reducing the uncertainty inherent to sediment budget prediction in torrent systems. Cette thèse est un recueil de projets d'études des processus de recharges sédimentaires des chenaux torrentiels. Ces travaux, réalisés lorsque j'étais employé à l'Université de Lausanne, se concentrent sur les implications géologiques et morphologiques des bassins dans l'apport de sédiments, élément fondamental dans la prédiction de laves torrentielles. D'autres aspects de dynamique sédimentaire ont été abordés, p. ex. le couplage torrent - bassin, ainsi qu'un modèle de simulation du transfert sédimentaire en milieu torrentiel. L'activité sédimentaire du Manival, un système torrentiel actif des Alpes françaises, a été étudiée par relevés au laser scanner terrestre et complétée par une étude géostructurale ainsi qu'un suivi du transfert en sédiments du torrent. Une année de flux sédimentaire a pu être observée, coïncidant avec deux laves torrentielles et plusieurs phénomènes de charriages. Cette étude a révélé que les laves s'étaient générées dans le torrent et étaient précédées par une recharge de débris depuis les versants. La production de débris s'est passée principalement en l'hiver - début du printemps, causée par de grandes ruptures de pentes. Le transfert était plus étrange, se produisant presque exclusivement au début du printemps subordonné aux conditions d'écoulement et en automne lors de longues pluies. Les orages d'été n'affectèrent guère les dépôts, qui semblent dépendre de leur stabilité. Les implications morpho-géologiques dans l'apport sédimentaire ont été évaluées à l'aide de MNT et études de terrain. Une classification de la topographie basée sur la pente a permis de charactériser le mode de production et transfert. Une analyse de stabilité de pente à partir des structures de roches a permis d'estimer la susceptibilité à la rupture. Les zones sources modélisées comprennent plus de 97% des chutes de blocs observées et les bilans sédimentaires sont corrélés à la densité de ruptures potentielles. Ce travail d'analyses des morphologies du terrain et de susceptibilité à la rupture documente la dynamique sédimentaire pour l'estimation quantitative des zones érosives induisant l'activité torrentielle. Le développement des formes d'érosion a été évalué par l'analyse de leur géométrie avec celle des ruptures potentielles et avec la direction de la fréquence maximale des joints. Les structures de roches, mais en particulier les dièdres et les discontinuités dominantes, semblent être très influents dans les mécanismes d'érosion affectant les bassins rocheux. Ils représentent des zones de faiblesse exploitées en priorité par les processus de démantèlement et d'érosion, encourageant l'initiation de ravines et couloirs, mais aussi leur propagation. L'incorporation du control géologique dans les processus de surface contribue à une meilleure compréhension de l'évolution topographique de bassins actifs. Un algorithme de flux sédimentaire a été implémenté dans un modèle en cascade, lequel divise le bassin en biefs et en systèmes individuels répondant aux processus. Chaque unité inclut de façon simple les informations géomorpologiques et celles du flux sédimentaire dérivées à partir de SIG et de cartographie de terrain. Cet outil permet la simulation des transferts de masse dans les chenaux, considérants la variabilité de l'apport et son transport, et aide à réduire l'incertitude liée à la prédiction de bilans sédimentaires torrentiels. Ce travail vise très humblement d'éclairer quelques aspects de la dynamique sédimentaire en milieu torrentiel.

Symmetries and fixed point stability of stochastic differential equations modeling self-organized criticality

A stochastic nonlinear partial differential equation is constructed for two different models exhibiting self-organized criticality: the Bak-Tang-Wiesenfeld (BTW) sandpile model [Phys. Rev. Lett. 59, 381 (1987); Phys. Rev. A 38, 364 (1988)] and the Zhang model [Phys. Rev. Lett. 63, 470 (1989)]. The dynamic renormalization group (DRG) enables one to compute the critical exponents. However, the nontrivial stable fixed point of the DRG transformation is unreachable for the original parameters of the models. We introduce an alternative regularization of the step function involved in the threshold condition, which breaks the symmetry of the BTW model. Although the symmetry properties of the two models are different, it is shown that they both belong to the same universality class. In this case the DRG procedure leads to a symmetric behavior for both models, restoring the broken symmetry, and makes accessible the nontrivial fixed point. This technique could also be applied to other problems with threshold dynamics.

Agreed upon procedures report on the Lee County Narcotics Task Force for the period July 1, 2010 through May 31, 2012

Agreed upon procedures report on the Lee County Narcotics Task Force for the period July 1, 2010 through May 31, 2012

MODELING OF SOIL LOAD-BEARING CAPACITY AS A FUNCTION OF SOIL MECHANICAL RESISTANCE TO PENETRATION

Estimation of soil load-bearing capacity from mathematical models that relate preconsolidation pressure (σp) to mechanical resistance to penetration (PR) and gravimetric soil water content (U) is important for defining strategies to prevent compaction of agricultural soils. Our objective was therefore to model the σp and compression index (CI) according to the PR (with an impact penetrometer in the field and a static penetrometer inserted at a constant rate in the laboratory) and U in a Rhodic Eutrudox. The experiment consisted of six treatments: no-tillage system (NT); NT with chiseling; and NT with additional compaction by combine traffic (passing 4, 8, 10, and 20 times). Soil bulk density, total porosity, PR (in field and laboratory measurements), U, σp, and CI values were determined in the 5.5-10.5 cm and 13.5-18.5 cm layers. Preconsolidation pressure (σp) and CI were modeled according to PR in different U. The σp increased and the CI decreased linearly with increases in the PR values. The correlations between σp and PR and PR and CI are influenced by U. From these correlations, the soil load-bearing capacity and compaction susceptibility can be estimated by PR readings evaluated in different U.

Surface Complexation Modeling in Variable Charge Soils: Charge Characterization by Potentiometric Titration

ABSTRACT Intrinsic equilibrium constants of 17 representative Brazilian Oxisols were estimated from potentiometric titration measuring the adsorption of H+ and OH− on amphoteric surfaces in suspensions of varying ionic strength. Equilibrium constants were fitted to two surface complexation models: diffuse layer and constant capacitance. The former was fitted by calculating total site concentration from curve fitting estimates and pH-extrapolation of the intrinsic equilibrium constants to the PZNPC (hand calculation), considering one and two reactive sites, and by the FITEQL software. The latter was fitted only by FITEQL, with one reactive site. Soil chemical and physical properties were correlated to the intrinsic equilibrium constants. Both surface complexation models satisfactorily fit our experimental data, but for results at low ionic strength, optimization did not converge in FITEQL. Data were incorporated in Visual MINTEQ and they provide a modeling system that can predict protonation-dissociation reactions in the soil surface under changing environmental conditions.

Surface Complexation Modeling in Variable Charge Soils: Prediction of Cadmium Adsorption

ABSTRACT Intrinsic equilibrium constants for 22 representative Brazilian Oxisols were estimated from a cadmium adsorption experiment. Equilibrium constants were fitted to two surface complexation models: diffuse layer and constant capacitance. Intrinsic equilibrium constants were optimized by FITEQL and by hand calculation using Visual MINTEQ in sweep mode, and Excel spreadsheets. Data from both models were incorporated into Visual MINTEQ. Constants estimated by FITEQL and incorporated in Visual MINTEQ software failed to predict observed data accurately. However, FITEQL raw output data rendered good results when predicted values were directly compared with observed values, instead of incorporating the estimated constants into Visual MINTEQ. Intrinsic equilibrium constants optimized by hand calculation and incorporated in Visual MINTEQ reliably predicted Cd adsorption reactions on soil surfaces under changing environmental conditions.

Task 4: Testing Iowa Portland Cement Concrete Mixtures for the AASHTO Mechanistic-Empirical Pavement Design Procedure, May 2008

The present research project was designed to identify the typical Iowa material input values that are required by the Mechanistic-Empirical Pavement Design Guide (MEPDG) for the Level 3 concrete pavement design. It was also designed to investigate the existing equations that might be used to predict Iowa pavement concrete for the Level 2 pavement design. In this project, over 20,000 data were collected from the Iowa Department of Transportation (DOT) and other sources. These data, most of which were concrete compressive strength, slump, air content, and unit weight data, were synthesized and their statistical parameters (such as the mean values and standard variations) were analyzed. Based on the analyses, the typical input values of Iowa pavement concrete, such as 28-day compressive strength (f’c), splitting tensile strength (fsp), elastic modulus (Ec), and modulus of rupture (MOR), were evaluated. The study indicates that the 28-day MOR of Iowa concrete is 646 + 51 psi, very close to the MEPDG default value (650 psi). The 28-day Ec of Iowa concrete (based only on two available data of the Iowa Curling and Warping project) is 4.82 + 0.28x106 psi, which is quite different from the MEPDG default value (3.93 x106 psi); therefore, the researchers recommend re-evaluating after more Iowa test data become available. The drying shrinkage (εc) of a typical Iowa concrete (C-3WR-C20 mix) was tested at Concrete Technology Laboratory (CTL). The test results show that the ultimate shrinkage of the concrete is about 454 microstrain and the time for the concrete to reach 50% of ultimate shrinkage is at 32 days; both of these values are very close to the MEPDG default values. The comparison of the Iowa test data and the MEPDG default values, as well as the recommendations on the input values to be used in MEPDG for Iowa PCC pavement design, are summarized in Table 20 of this report. The available equations for predicting the above-mentioned concrete properties were also assembled. The validity of these equations for Iowa concrete materials was examined. Multiple-parameters nonlinear regression analyses, along with the artificial neural network (ANN) method, were employed to investigate the relationships among Iowa concrete material properties and to modify the existing equations so as to be suitable for Iowa concrete materials. However, due to lack of necessary data sets, the relationships between Iowa concrete properties were established based on the limited data from CP Tech Center’s projects and ISU classes only. The researchers suggest that the resulting relationships be used by Iowa pavement design engineers as references only. The present study furthermore indicates that appropriately documenting concrete properties, including flexural strength, elastic modulus, and information on concrete mix design, is essential for updating the typical Iowa material input values and providing rational prediction equations for concrete pavement design in the future.

Task 6: Material Thermal Input for Iowa Materials, May 2008

The present research project was designed to determine thermal properties, such as coefficient of thermal expansion (CTE) and thermal conductivity, of Iowa concrete pavement materials. These properties are required as input values by the Mechanistic-Empirical Pavement Design Guide (MEPDG). In this project, a literature review was conducted to determine the factors that affect thermal properties of concrete and the existing prediction equations for CTE and thermal conductivity of concrete. CTE tests were performed on various lab and field samples of portland cement concrete (PCC) at the Iowa Department of Transportation and Iowa State University. The variations due to the test procedure, the equipment used, and the consistency of field batch materials were evaluated. The test results showed that the CTE variations due to test procedure and batch consistency were less than 5%, and the variation due to the different equipment was less than 15%. Concrete CTE values were significantly affected by different types of coarse aggregate. The CTE values of Iowa concrete made with limestone+graval, quartzite, dolomite, limestone+dolomite, and limestone were 7.27, 6.86, 6.68, 5.83, and 5.69 microstrain/oF (13.08, 12.35, 12.03, 10.50, and 10.25 microstrain/oC), respectively, which were all higher than the default value of 5.50 microstrain/oF in the MEPDG program. The thermal conductivity of a typical Iowa PCC mix and an asphalt cement concrete (ACC) mix (both with limestone as coarse aggregate) were tested at Concrete Technology Laboratory in Skokie, Illinois. The thermal conductivity was 0.77 Btu/hr•ft•oF (1.33 W/m•K) for PCC and 1.21 Btu/hr•ft•oF (2.09 W/m•K) for ACC, which are different from the default values (1.25 Btu/hr•ft•oF or 2.16 W/m•K for PCC and 0.67 Btu/hr•ft•oF or 1.16 W/m•K for ACC) in the MEPDG program. The investigations onto the CTE of ACC and the effects of concrete materials (such as cementitious material and aggregate types) and mix proportions on concrete thermal conductivity are recommended to be considered in future studies.

EFNS guidelines on cognitive rehabilitation: report of an EFNS task force.

Disorders of language, spatial perception, attention, memory, calculation and praxis are a frequent consequence of acquired brain damage [in particular, stroke and traumatic brain injury (TBI)] and a major determinant of disability. The rehabilitation of aphasia and, more recently, of other cognitive disorders is an important area of neurological rehabilitation. We report here a review of the available evidence about effectiveness of cognitive rehabilitation. Given the limited number and generally low quality of randomized clinical trials (RCTs) in this area of therapeutic intervention, the Task Force considered, besides the available Cochrane reviews, evidence of lower classes which was critically analysed until a consensus was reached. In particular, we considered evidence from small group or single cases studies including an appropriate statistical evaluation of effect sizes. The general conclusion is that there is evidence to award a grade A, B or C recommendation to some forms of cognitive rehabilitation in patients with neuropsychological deficits in the post-acute stage after a focal brain lesion (stroke, TBI). These include aphasia therapy, rehabilitation of unilateral spatial neglect (ULN), attentional training in the post-acute stage after TBI, the use of electronic memory aids in memory disorders, and the treatment of apraxia with compensatory strategies. There is clearly a need for adequately designed studies in this area, which should take into account specific problems such as patient heterogeneity and treatment standardization.

Concrete Overlay Field Application Program Iowa Task Report US 18 Concrete Overlay Construction Under Traffic, May 2012

The National Concrete Pavement Technology Center, Iowa Department of Transportation, and Federal Highway Administration set out to demonstrate and document the design and construction of portland cement concrete (PCC) overlays on two-lane roadways while maintaining two-way traffic. An 18.82 mile project was selected for 2011 construction in northeast Iowa on US 18 between Fredericksburg and West Union. This report documents planning, design, and construction of the project and lessons learned. The work included the addition of subdrains, full-depth patching, bridge approach replacement, and drainage structural repair and cleaning prior to overlay construction. The paving involved surface preparation by milling to grade and the placement of a 4.5 inch PCC overlay and 4 foot of widening to the existing pavement. In addition, the report makes recommendations on ways to improve the process for future concrete overlays.

High-throughput SELEX SAGE method for quantitative modeling of transcription-factor binding sites.

The ability to determine the location and relative strength of all transcription-factor binding sites in a genome is important both for a comprehensive understanding of gene regulation and for effective promoter engineering in biotechnological applications. Here we present a bioinformatically driven experimental method to accurately define the DNA-binding sequence specificity of transcription factors. A generalized profile was used as a predictive quantitative model for binding sites, and its parameters were estimated from in vitro-selected ligands using standard hidden Markov model training algorithms. Computer simulations showed that several thousand low- to medium-affinity sequences are required to generate a profile of desired accuracy. To produce data on this scale, we applied high-throughput genomics methods to the biochemical problem addressed here. A method combining systematic evolution of ligands by exponential enrichment (SELEX) and serial analysis of gene expression (SAGE) protocols was coupled to an automated quality-controlled sequence extraction procedure based on Phred quality scores. This allowed the sequencing of a database of more than 10,000 potential DNA ligands for the CTF/NFI transcription factor. The resulting binding-site model defines the sequence specificity of this protein with a high degree of accuracy not achieved earlier and thereby makes it possible to identify previously unknown regulatory sequences in genomic DNA. A covariance analysis of the selected sites revealed non-independent base preferences at different nucleotide positions, providing insight into the binding mechanism.

Outcomes from studies of Antineutrophil Cytoplasm Antibody Associated Vasculitis: a systematic review by the European League Against Rheumatism systemic vasculitis task force.

Objectives: We undertook a systematic literature review as a background to the European League Against Rheumatism (EULAR) recommendations for conducting clinical trials in anti-neutrophil cytoplasm antibody associated vasculitis (AAV), and to assess the quality of evidence for outcome measures in AAV. Methods: Using a systematic Medline search, we categorised the identified studies according to diagnoses. Factors affecting remission, relapse, renal function and overall survival were identified. Results: A total of 44 papers were reviewed from 502 identified by our search criteria. There was considerable inconsistency in definitions of end points. Remission rates varied from 30% to 93% in Wegener granulomatosis (WG), 75% to 89% in microscopic polyangiitis (MPA) and 81% to 91% in Churg¿Strauss syndrome (CSS). The 5-year survival for WG, MPA and CSS was 74¿91%, 45¿76% and 60¿97%. Relapse (variably defined) was common in the first 2 years but the frequency varied: 18% to 60% in WG, 8% in MPA, and 35% in CSS. The rate of renal survival in WG varied from 23% at 15 months to 23% at 120 months. Methods used to assess morbidity varied between studies. Ignoring the variations in definitions of the stage of disease, factors influencing remission, relapse, renal and overall survival included immunosuppressive therapy used, type of organ involvement, presence of ANCA, older age and male ender. Conclusions: Factors influencing remission, relapse, renal and overall survival include the type of immunosuppressive therapy used, pattern of organ involvement, presence of ANCA, older age and male gender. Methodological variations between studies highlight the need for a consensus on terminology and definitions for future conduct of clinical studies in AAV.

Missouri River Flood Coordination Task Force Report, September 25, 2011

The Missouri River floods of 2011 will go down in history as the longest duration flooding event this state has seen to date. The combination of above normal snowfall in the upper Missouri River basin followed by the equivalent of nearly one year’s worth of rainfall in May created an above normal runoff situation which filled the Missouri River and the six main reservoirs within the basin. Compounding this problem was colder than normal temperatures which kept much of the snowpack in the upper basin on the ground longer into the spring, setting the stage for this historic event. The U.S. Army Corps of Engineers (USACE) began increasing the outflow at Gavin’s Point, near Yankton, South Dakota in May. On June 14, 2011, the outflow reached a record rate of over 160,000 cubic feet per second (cfs), over twice the previous record outflow set in 1997. This increased output from Gavin’s Point caused the Missouri River to flow out of its banks covering over 283,000 acres of land in Iowa, forcing hundreds of evacuations, damaging 255,000 acres of cropland and significantly impacting the levee system on the Missouri River basin. Over the course of the summer, approximately 64 miles of primary roads closed due to Missouri River flooding, including 54 miles of Interstate Highway. Many county secondary roads were closed by high water or overburdened due to the numerous detours and road closures in this area. As the Missouri River levels began to increase, municipalities and counties aided by State and Federal agencies began preparing for a sustained flood event. Citizens, businesses, state agencies, local governments and non‐profits made substantial preparations, in some cases expending millions of dollars on emergency protective measures to protect their facilities from the impending flood. Levee monitors detected weak spots in the levee system in all affected counties, with several levees being identified as at risk levees that could potentially fail. Of particular concern was the 28 miles of levees protecting Council Bluffs. Based on this concern, Council Bluffs prepared an evacuation plan for the approximately 30,000 residents that resided in the protected area. On May 25, 2011, Governor Branstad directed the execution of the Iowa Emergency Response Plan in accordance with Section 401 of the Stafford Act. On May 31, 2011, HSEMD Administrator, Brigadier General J. Derek Hill, formally requested the USACE to provide technical assistance and advanced measures for the communities along the Missouri River basin. On June 2, 2011 Governor Branstad issued a State of Iowa Proclamation of Disaster Emergency for Fremont, Harrison, Mills, Monona, Pottawattamie, and Woodbury counties. The length of this flood event created a unique set of challenges for Federal, State and local entities. In many cases, these organizations were conducting response and recovery operations simultaneously. Due to the length of this entire event, the State Emergency Operations Center and the local Emergency Operations Centers remained open for an extended period of time, putting additional strain on many organizations and resources. In response to this disaster, Governor Branstad created the Missouri River Recovery Coordination Task Force to oversee the State’s recovery efforts. The Governor announced the creation of this Task Force on October 17, 2011 and appointed Brigadier General J. Derek Hill, HSEMD Administrator as the chairman. This Task Force would be a temporary group of State agency representatives and interested stakeholders brought together to support the recovery efforts of the Iowa communities impacted by the Missouri River Flood. Collectively, this group would analyze and share damage assessment data, coordinate assistance across various stakeholders, monitor progress, capture best practices and identify lessons learned.

Neural networks as mechanisms to regulate division of labor.

Abstract In social insects, workers perform a multitude of tasks, such as foraging, nest construction, and brood rearing, without central control of how work is allocated among individuals. It has been suggested that workers choose a task by responding to stimuli gathered from the environment. Response-threshold models assume that individuals in a colony vary in the stimulus intensity (response threshold) at which they begin to perform the corresponding task. Here we highlight the limitations of these models with respect to colony performance in task allocation. First, we show with analysis and quantitative simulations that the deterministic response-threshold model constrains the workers' behavioral flexibility under some stimulus conditions. Next, we show that the probabilistic response-threshold model fails to explain precise colony responses to varying stimuli. Both of these limitations would be detrimental to colony performance when dynamic and precise task allocation is needed. To address these problems, we propose extensions of the response-threshold model by adding variables that weigh stimuli. We test the extended response-threshold model in a foraging scenario and show in simulations that it results in an efficient task allocation. Finally, we show that response-threshold models can be formulated as artificial neural networks, which consequently provide a comprehensive framework for modeling task allocation in social insects.