Detection, analysis and monitoring of slope movements by high-resolution digital elevation models

Summary Detection, analysis and monitoring of slope movements by high-resolution digital elevation modelsSlope movements, such as rockfalls, rockslides, shallow landslides or debris flows, are frequent in many mountainous areas. These natural hazards endanger the inhabitants and infrastructures making it necessary to assess the hazard and risk caused by these phenomena. This PhD thesis explores various approaches using digital elevation models (DEMs) - and particularly high-resolution DEMs created by aerial or terrestrial laser scanning (TLS) - that contribute to the assessment of slope movement hazard at regional and local scales.The regional detection of areas prone to rockfalls and large rockslides uses different morphologic criteria or geometric instability factors derived from DEMs, i.e. the steepness of the slope, the presence of discontinuities, which enable a sliding mechanism, and the denudation potential. The combination of these factors leads to a map of susceptibility to rockfall initiation that is in good agreement with field studies as shown with the example of the Little Mill Campground area (Utah, USA). Another case study in the Illgraben catchment in the Swiss Alps highlighted the link between areas with a high denudation potential and actual rockfall areas.Techniques for a detailed analysis and characterization of slope movements based on high-resolution DEMs have been developed for specific, localized sites, i.e. ancient slide scars, present active instabilities or potential slope instabilities. The analysis of the site's characteristics mainly focuses on rock slopes and includes structural analyses (orientation of discontinuities); estimation of spacing, persistence and roughness of discontinuities; failure mechanisms based on the structural setting; and volume calculations. For the volume estimation a new 3D approach was tested to reconstruct the topography before a landslide or to construct the basal failure surface of an active or potential instability. The rockslides at Åknes, Tafjord and Rundefjellet in western Norway were principally used as study sites to develop and test the different techniques.The monitoring of slope instabilities investigated in this PhD thesis is essentially based on multitemporal (or sequential) high-resolution DEMs, in particular sequential point clouds acquired by TLS. The changes in the topography due to slope movements can be detected and quantified by sequential TLS datasets, notably by shortest distance comparisons revealing the 3D slope movements over the entire region of interest. A detailed analysis of rock slope movements is based on the affine transformation between an initial and a final state of the rock mass and its decomposition into translational and rotational movements. Monitoring using TLS was very successful on the fast-moving Eiger rockslide in the Swiss Alps, but also on the active rockslides of Åknes and Nordnesfjellet (northern Norway). One of the main achievements on the Eiger and Aknes rockslides is to combine the site's morphology and structural setting with the measured slope movements to produce coherent instability models. Both case studies also highlighted a strong control of the structures in the rock mass on the sliding directions. TLS was also used to monitor slope movements in soils, such as landslides in sensitive clays in Québec (Canada), shallow landslides on river banks (Sorge River, Switzerland) and a debris flow channel (Illgraben).The PhD thesis underlines the broad uses of high-resolution DEMs and especially of TLS in the detection, analysis and monitoring of slope movements. Future studies should explore in more depth the different techniques and approaches developed and used in this PhD, improve them and better integrate the findings in current hazard assessment practices and in slope stability models.Résumé Détection, analyse et surveillance de mouvements de versant à l'aide de modèles numériques de terrain de haute résolutionDes mouvements de versant, tels que des chutes de blocs, glissements de terrain ou laves torrentielles, sont fréquents dans des régions montagneuses et mettent en danger les habitants et les infrastructures ce qui rend nécessaire d'évaluer le danger et le risque causé par ces phénomènes naturels. Ce travail de thèse explore diverses approches qui utilisent des modèles numériques de terrain (MNT) et surtout des MNT de haute résolution créés par scanner laser terrestre (SLT) ou aérien - et qui contribuent à l'évaluation du danger de mouvements de versant à l'échelle régionale et locale.La détection régionale de zones propices aux chutes de blocs ou aux éboulements utilise plusieurs critères morphologiques dérivés d'un MNT, tels que la pente, la présence de discontinuités qui permettent un mécanisme de glissement ou le potentiel de dénudation. La combinaison de ces facteurs d'instabilité mène vers une carte de susceptibilité aux chutes de blocs qui est en accord avec des travaux de terrain comme démontré avec l'exemple du Little Mill Campground (Utah, États-Unis). Un autre cas d'étude - l'Illgraben dans les Alpes valaisannes - a mis en évidence le lien entre les zones à fort potentiel de dénudation et les sources effectives de chutes de blocs et d'éboulements.Des techniques pour l'analyse et la caractérisation détaillée de mouvements de versant basées sur des MNT de haute résolution ont été développées pour des sites spécifiques et localisés, comme par exemple des cicatrices d'anciens éboulements et des instabilités actives ou potentielles. Cette analyse se focalise principalement sur des pentes rocheuses et comprend l'analyse structurale (orientation des discontinuités); l'estimation de l'espacement, la persistance et la rugosité des discontinuités; l'établissement des mécanismes de rupture; et le calcul de volumes. Pour cela une nouvelle approche a été testée en rétablissant la topographie antérieure au glissement ou en construisant la surface de rupture d'instabilités actuelles ou potentielles. Les glissements rocheux d'Åknes, Tafjord et Rundefjellet en Norvège ont été surtout utilisés comme cas d'étude pour développer et tester les diverses approches. La surveillance d'instabilités de versant effectuée dans cette thèse de doctorat est essentiellement basée sur des MNT de haute résolution multi-temporels (ou séquentiels), en particulier des nuages de points séquentiels acquis par SLT. Les changements topographiques dus aux mouvements de versant peuvent être détectés et quantifiés sur l'ensemble d'un glissement, notamment par comparaisons des distances les plus courtes entre deux nuages de points. L'analyse détaillée des mouvements est basée sur la transformation affine entre la position initiale et finale d'un bloc et sa décomposition en mouvements translationnels et rotationnels. La surveillance par SLT a démontré son potentiel avec l'effondrement d'un pan de l'Eiger dans les Alpes suisses, mais aussi aux glissements rocheux d'Aknes et Nordnesfjellet en Norvège. Une des principales avancées à l'Eiger et à Aknes est la création de modèles d'instabilité cohérents en combinant la morphologie et l'agencement structural des sites avec les mesures de déplacements. Ces deux cas d'étude ont aussi démontré le fort contrôle des structures existantes dans le massif rocheux sur les directions de glissement. Le SLT a également été utilisé pour surveiller des glissements dans des terrains meubles comme dans les argiles sensibles au Québec (Canada), sur les berges de la rivière Sorge en Suisse et dans le chenal à laves torrentielles de l'Illgraben.Cette thèse de doctorat souligne le vaste champ d'applications des MNT de haute résolution et particulièrement du SLT dans la détection, l'analyse et la surveillance des mouvements de versant. Des études futures devraient explorer plus en profondeur les différentes techniques et approches développées, les améliorer et mieux les intégrer dans des pratiques actuelles d'analyse de danger et surtout dans la modélisation de stabilité des versants.

Linear relations between writhe and minimal crossing number in Conway families of ideal knots and links

We showed earlier how to predict the writhe of any rational knot or link in its ideal geometric configuration, or equivalently the average of the 3D writhe over statistical ensembles of random configurations of a given knot or link (Cerf and Stasiak 2000 Proc. Natl Acad. Sci. USA 97 3795). There is no general relation between the minimal crossing number of a knot and the writhe of its ideal geometric configuration. However, within individual families of knots linear relations between minimal crossing number and writhe were observed (Katritch et al 1996 Nature 384 142). Here we present a method that allows us to express the writhe as a linear function of the minimal crossing number within Conway families of knots and links in their ideal configuration. The slope of the lines and the shift between any two lines with the same

Origin of a restraining bend in an evolving strike-slip system: The Cenozoic As Pontes basin (NW Spain)

The As Pontes basin (12 km2), NW Iberian Peninsula, is bounded by a double restraining bend of a dextral strike-slip fault, which is related to the western onshore end of the Pyrenean belt. Surface and subsurface data obtained from intensive coal exploration and mining in the basin since the 1960s together with additional structural and stratigraphic sequence analysis allowed us to determine the geometric relationships between tectonic structures and stratigraphic markers. The small size of the basin and the large amount of quality data make the As Pontes basin a unique natural laboratory for improving our understanding of the origin and evolution of restraining bends. The double restraining bend is the end stage of the structural evolution of a compressive underlapping stepover, where the basin was formed. During the first stage (stepover stage), which began ca. 30 Ma ago (latest Rupelian) and lasted 3.4 My, two small isolated basins bounded by thrusts and normal faults were formed. For 1.3 My, the strike-slip faults, which defined the stepover, grew towards each other until joining and forming the double restraining bend, which bounds one large As Pontes basin (transition stage). The history of the basin was controlled by the activity of the double restraining bend for a further 3.4 My (restraining bend stage) and ended in mid-Aquitanian times (ca. 22 Ma).

Impacts on Safety of Left-Turn Treatment at High Speed Signalized Intersections, HR-347, 1994

Left-turning traffic is a major source of conflicts at intersections. Though an average of only 10% to 15% of all approach traffic turns left, these vehicles are involved in approximately 45% of all accidents. This report presents the results of research conducted to develop models which estimate approach accident rates at high speed signalized intersections. The objective of the research was to quantify the relationship between traffic and intersection characteristics, and accident potential of different left turn treatments. Geometric, turning movement counts, and traffic signal phasing data were collected at 100 intersections in Iowa using a questionnaire sent to municipalities. Not all questionnaires resulted in complete data and ultimately complete data were derived for 63 intersections providing a database of 248 approaches. Accident data for the same approaches were obtained from the Iowa Department of Transportation Accident Location and Analysis System (ALAS). Regression models were developed for two different dependent variables: 1) the ratio of the number of left turn accidents per approach to million left turning vehicles per approach, and 2) the ratio of accidents per approach to million traffic movements per approach. A number of regression models were developed for both dependent variables. One model using each dependent variable was developed for intersections with low, medium, and high left turning traffic volumes. As expected, the research indicates that protected left turn phasing has a lower accident potential than protected/permitted or permitted phasing. Left turn lanes and multiple lane approaches are beneficial for reducing accident rates, while raised medians increase the likelihood of accidents. Signals that are part of a signal system tend to have lower accident rates than isolated signals. The resulting regression models may be used to determine the likely impact of various left turn treatments on intersection accident rates. When designing an intersection approach, a traffic engineer may use the models to estimate the accident rate reduction as a result of improved lane configurations and left turn treatments. The safety benefits may then be compared to any costs associated with operational effects to the intersection (i.e., increased delay) to determine the benefits and costs of making intersection safety improvements.

Common variants in UMOD associate with urinary uromodulin levels: a meta-analysis.

Uromodulin is expressed exclusively in the thick ascending limb and is the most abundant protein excreted in normal urine. Variants in UMOD, which encodes uromodulin, are associated with renal function, and urinary uromodulin levels may be a biomarker for kidney disease. However, the genetic factors regulating uromodulin excretion are unknown. We conducted a meta-analysis of urinary uromodulin levels to identify associated common genetic variants in the general population. We included 10,884 individuals of European descent from three genetic isolates and three urban cohorts. Each study measured uromodulin indexed to creatinine and conducted linear regression analysis of approximately 2.5 million single nucleotide polymorphisms using an additive model. We also tested whether variants in genes expressed in the thick ascending limb associate with uromodulin levels. rs12917707, located near UMOD and previously associated with renal function and CKD, had the strongest association with urinary uromodulin levels (P<0.001). In all cohorts, carriers of a G allele of this variant had higher uromodulin levels than noncarriers did (geometric means 10.24, 14.05, and 17.67 μg/g creatinine for zero, one, or two copies of the G allele). rs12446492 in the adjacent gene PDILT (protein disulfide isomerase-like, testis expressed) also reached genome-wide significance (P<0.001). Regarding genes expressed in the thick ascending limb, variants in KCNJ1, SORL1, and CAB39 associated with urinary uromodulin levels. These data indicate that common variants in the UMOD promoter region may influence urinary uromodulin levels. They also provide insights into uromodulin biology and the association of UMOD variants with renal function.

Composite Precast Prestressed Concrete Bridge Slabs, HR-310, 1991

Precast prestressed concrete panels have been used as subdecks in bridge construction in Iowa and other states. To investigate the performance of these types of composite slabs at locations adjacent to abutment and pier diaphragms in skewed bridges, a research prcject which involved surveys of design agencies and precast producers, field inspections of existing bridges, analytical studies, and experimental testing was conducted. The survey results from the design agencies and panel producers showed that standardization of precast panel construction would be desirable, that additional inspections at the precast plant and at the bridge site would be beneficial, and that some form of economical study should be undertaken to determine actual cost savings associated with composite slab construction. Three bridges in Hardin County, Iowa were inspected to observe general geometric relationships, construction details, and to note the visual condition of the bridges. Hairline cracks beneath several of the prestressing strands in many of the precast panels were observed, and a slight discoloration of the concrete was seen beneath most of the strands. Also, some rust staining was visible at isolated locations on several panels. Based on the findings of these inspections, future inspections are recommended to monitor the condition of these and other bridges constructed with precast panel subdecks. Five full-scale composite slab specimens were constructed in the Structural Engineering Laboratory at Iowa State University. One specimen modeled bridge deck conditions which are not adjacent to abutment or pier diaphragms, and the other four specimens represented the geometric conditions which occur for skewed diaphragms of 0, 15, 30, and 40 degrees. The specimens were subjected to wheel loads of service and factored level magnitudes at many locations on the slab surface and to concentrated loads which produced failure of the composite slab. The measured slab deflections and bending strains at both service and factored load levels compared reasonably well with the results predicted by simplified Finite element analyses of the specimens. To analytically evaluate the nominal strength for a composite slab specimen, yield-line and punching shear theories were applied. Yield-line limit loads were computed using the crack patterns generated during an ultimate strength test. In most cases, these analyses indicated that the failure mode was not flexural. Since the punching shear limit loads in most instances were close to the failure loads, and since the failure surfaces immediately adjacent to the wheel load footprint appeared to be a truncated prism shape, the probable failure mode for all of the specimens was punching shear. The development lengths for the prestressing strands in the rectangular and trapezoidal shaped panels was qualitatively investigated by monitoring strand slippage at the ends of selected prestressing strands. The initial strand transfer length was established experimentally by monitoring concrete strains during strand detensioning, and this length was verified analytically by a finite element analysis. Even though the computed strand embedment lengths in the panels were not sufficient to fully develop the ultimate strand stress, sufficient stab strength existed. Composite behavior for the slab specimens was evaluated by monitoring slippage between a panel and the topping slab and by computation of the difference in the flexural strains between the top of the precast panel and the underside of the topping slab at various locations. Prior to the failure of a composite slab specimen, a localized loss of composite behavior was detected. The static load strength performance of the composite slab specimens significantly exceeded the design load requirements. Even with skew angles of up to 40 degrees, the nominal strength of the slabs did not appear to be affected when the ultimate strength test load was positioned on the portion of each slab containing the trapezoidal-shaped panel. At service and factored level loads, the joint between precast panels did not appear to influence the load distribution along the length of the specimens. Based on the static load strength of the composite slab specimens, the continued use of precast panels as subdecks in bridge deck construction is recommended.

Fast Geodesic Active Fields for Image Registration Based on Splitting and Augmented Lagrangian Approaches.

In this paper, we present an efficient numerical scheme for the recently introduced geodesic active fields (GAF) framework for geometric image registration. This framework considers the registration task as a weighted minimal surface problem. Hence, the data-term and the regularization-term are combined through multiplication in a single, parametrization invariant and geometric cost functional. The multiplicative coupling provides an intrinsic, spatially varying and data-dependent tuning of the regularization strength, and the parametrization invariance allows working with images of nonflat geometry, generally defined on any smoothly parametrizable manifold. The resulting energy-minimizing flow, however, has poor numerical properties. Here, we provide an efficient numerical scheme that uses a splitting approach; data and regularity terms are optimized over two distinct deformation fields that are constrained to be equal via an augmented Lagrangian approach. Our approach is more flexible than standard Gaussian regularization, since one can interpolate freely between isotropic Gaussian and anisotropic TV-like smoothing. In this paper, we compare the geodesic active fields method with the popular Demons method and three more recent state-of-the-art algorithms: NL-optical flow, MRF image registration, and landmark-enhanced large displacement optical flow. Thus, we can show the advantages of the proposed FastGAF method. It compares favorably against Demons, both in terms of registration speed and quality. Over the range of example applications, it also consistently produces results not far from more dedicated state-of-the-art methods, illustrating the flexibility of the proposed framework.

A multimodality vascular imaging phantom with fiducial markers visible in DSA, CTA, MRA, and ultrasound.

The objective was to design a vascular phantom compatible with digital subtraction angiography, computerized tomography angiography, ultrasound and magnetic resonance angiography (MRA). Fiducial markers were implanted at precise known locations in the phantom to facilitate identification and orientation of plane views from three-dimensional (3-D) reconstructed images. A vascular conduit connected to tubing at the extremities of the phantom ran through an agar-based gel filling it. A vessel wall in latex was included around the conduit to avoid diffusion of contrast agents. Using a lost-material casting technique based on a low melting point metal, geometries of pathological vessels were modeled. During the experimental testing, fiducial markers were detectable in all modalities without distortion. No leak of gadolinium through the vascular wall was observed on MRA after 5 hours. Moreover, no significant deformation of the vascular conduit was noted during the fabrication process (confirmed by microtome slicing along the vessel). The potential use of the phantom for calibration, rescaling, and fusion of 3-D images obtained from the different modalities as well as its use for the evaluation of intra- and inter-modality comparative studies of imaging systems are discussed. In conclusion, the vascular phantom can allow accurate calibration of radiological imaging devices based on x-ray, magnetic resonance and ultrasound and quantitative comparisons of the geometric accuracy of the vessel lumen obtained with each of these methods on a given well defined 3-D geometry.

Shielding requirements in helical tomotherapy.

Helical tomotherapy is a relatively new intensity-modulated radiation therapy (IMRT) treatment for which room shielding has to be reassessed for the following reasons. The beam-on-time needed to deliver a given target dose is increased and leads to a weekly workload of typically one order of magnitude higher than that for conventional radiation therapy. The special configuration of tomotherapy units does not allow the use of standard shielding calculation methods. A conventional linear accelerator must be shielded for primary, leakage and scatter photon radiations. For tomotherapy, primary radiation is no longer the main shielding issue since a beam stop is mounted on the gantry directly opposite the source. On the other hand, due to the longer irradiation time, the accelerator head leakage becomes a major concern. An analytical model based on geometric considerations has been developed to determine leakage radiation levels throughout the room for continuous gantry rotation. Compared to leakage radiation, scatter radiation is a minor contribution. Since tomotherapy units operate at a nominal energy of 6 MV, neutron production is negligible. This work proposes a synthetic and conservative model for calculating shielding requirements for the Hi-Art II TomoTherapy unit. Finally, the required concrete shielding thickness is given for different positions of interest.

Effects of Road Construction Intensity and Operations on Rural Freeway Work Zone Capacity, TPF-5(081), 2014

Capacity is affected by construction type and its intensity on adjacent open traffic lanes. The effect on capacity is a function of vehicles moving in and out of the closed lanes of the work zone, and the presence of heavy construction vehicles. Construction activity and its intensity, however, are not commonly considered in estimating capacity of a highway lane. The main purpose of this project was to attempt to quantify the effects of construction type and intensity (e.g. maintenance, rehabilitation, reconstruction, and milling) on work zone capacity. The objective of this project is to quantify the effects of construction type and its intensity on work zone capacity and to develop guidelines for MoDOT to estimate the specific operation type and intensity that will improve the traffic flow by reducing the traffic flow and queue length commonly associated with work zones. Despite the effort put into field data collection, the data collected did not show a full speed-flow chart therefore extracting a reliable capacity value was difficult. A statistical comparison between the capacity values found in this study using either methodologies indicates that there is an effect of construction activity on the values work zone capacity. It was found that the heavy construction activity reduces the capacity. It is very beneficial to conduct similar studies on the capacity of work zone with different lane closure barriers, which is also directly related to the type of work zone being short-term or long-term work zones. Also, the effect of different geometric and environmental characteristics of the roadway should be considered in future studies.

Uniform estimates in the Poincare-Aronszajn theorem on the separation of singularities of analytic functions

We study the possibility of splitting any bounded analytic function $f$ with singularities in a closed set $E\cup F$ as a sum of two bounded analytic functions with singularities in $E$ and $F$ respectively. We obtain some results under geometric restrictions on the sets $E$ and $F$ and we provide some examples showing the sharpness of the positive results.

Enhanced compressed sensing recovery with level set normals.

We propose a compressive sensing algorithm that exploits geometric properties of images to recover images of high quality from few measurements. The image reconstruction is done by iterating the two following steps: 1) estimation of normal vectors of the image level curves, and 2) reconstruction of an image fitting the normal vectors, the compressed sensing measurements, and the sparsity constraint. The proposed technique can naturally extend to nonlocal operators and graphs to exploit the repetitive nature of textured images to recover fine detail structures. In both cases, the problem is reduced to a series of convex minimization problems that can be efficiently solved with a combination of variable splitting and augmented Lagrangian methods, leading to fast and easy-to-code algorithms. Extended experiments show a clear improvement over related state-of-the-art algorithms in the quality of the reconstructed images and the robustness of the proposed method to noise, different kind of images, and reduced measurements.

Evaluation of the Confirmation Light System at Signalized Intersections in Medium to Large Iowa Communities to Reduce Red Light Running Violations”, TR-657, 2015

Red light running continues to be a serious safety concern for many communities in the United States. The Federal Highway Administration reported that in 2011, red light running accounted for 676 fatalities nationwide. Red light running crashes at a signalized intersections are more serious, especially in high speed corridors where speeds are above 35 mph. Many communities have invested in red light countermeasures including low-cost strategies (e.g. signal backplates, targeted enforcement, signal timing adjustments and improvement with signage) to high-cost strategies (e.g. automated enforcement and intersection geometric improvements). This research study investigated intersection confirmation lights as a low-cost strategy to reduce red light running violations. Two intersections in Altoona and Waterloo, Iowa were equipped with confirmation lights which targeted the through and left turning movements. Confirmation lights enable a single police officer to monitor a specific lane of traffic downstream of the intersection. A before-after analysis was conducted in which a change in red light running violations prior to- and 1 and 3 months after installation were evaluated. A test of proportions was used to determine if the change in red light running violation rates were statistically significant at the 90 and 95 percent levels of confidence. The two treatment intersections were then compared to the changes of red light running violation rates at spillover intersections (directly adjacent to the treatment intersections) and control intersections. The results of the analysis indicated a 10 percent reduction of red light running violations in Altoona and a 299 percent increase in Waterloo at the treatment locations. Finally, the research team investigated the time into red for each observed red light running violation. The analysis indicated that many of the violations occurred less than one second into the red phase and that most of the violation occurred during or shortly after the all-red phase.

Muscle moment-arms: a key element in muscle-force estimation.

A clear and rigorous definition of muscle moment-arms in the context of musculoskeletal systems modelling is presented, using classical mechanics and screw theory. The definition provides an alternative to the tendon excursion method, which can lead to incorrect moment-arms if used inappropriately due to its dependency on the choice of joint coordinates. The definition of moment-arms, and the presented construction method, apply to musculoskeletal models in which the bones are modelled as rigid bodies, the joints are modelled as ideal mechanical joints and the muscles are modelled as massless, frictionless cables wrapping over the bony protrusions, approximated using geometric surfaces. In this context, the definition is independent of any coordinate choice. It is then used to solve a muscle-force estimation problem for a simple 2D conceptual model and compared with an incorrect application of the tendon excursion method. The relative errors between the two solutions vary between 0% and 100%.