Grade and Cross Slope Estimation from LIDAR-based Surface Models, 2003

Many transportation agencies maintain grade as an attribute in roadway inventory databases; however, the information is often in an aggregated format. Cross slope is rarely included in large roadway inventories. Accurate methods available to collect grade and cross slope include global positioning systems, traditional surveying, and mobile mapping systems. However, most agencies do not have the resources to utilize these methods to collect grade and cross slope on a large scale. This report discusses the use of LIDAR to extract roadway grade and cross slope for large-scale inventories. Current data collection methods and their advantages and disadvantages are discussed. A pilot study to extract grade and cross slope from a LIDAR data set, including methodology, results, and conclusions, is presented. This report describes the regression methodology used to extract and evaluate the accuracy of grade and cross slope from three dimensional surfaces created from LIDAR data. The use of LIDAR data to extract grade and cross slope on tangent highway segments was evaluated and compared against grade and cross slope collected using an automatic level for 10 test segments along Iowa Highway 1. Grade and cross slope were measured from a surface model created from LIDAR data points collected for the study area. While grade could be estimated to within 1%, study results indicate that cross slope cannot practically be estimated using a LIDAR derived surface model.

Simulation of surface waves in porous media

We present a novel numerical algorithm for the simulation of seismic wave propagation in porous media, which is particularly suitable for the accurate modelling of surface wave-type phenomena. The differential equations of motion are based on Biot's theory of poro-elasticity and solved with a pseudospectral approach using Fourier and Chebyshev methods to compute the spatial derivatives along the horizontal and vertical directions, respectively. The time solver is a splitting algorithm that accounts for the stiffness of the differential equations. Due to the Chebyshev operator the grid spacing in the vertical direction is non-uniform and characterized by a denser spatial sampling in the vicinity of interfaces, which allows for a numerically stable and accurate evaluation of higher order surface wave modes. We stretch the grid in the vertical direction to increase the minimum grid spacing and reduce the computational cost. The free-surface boundary conditions are implemented with a characteristics approach, where the characteristic variables are evaluated at zero viscosity. The same procedure is used to model seismic wave propagation at the interface between a fluid and porous medium. In this case, each medium is represented by a different grid and the two grids are combined through a domain-decomposition method. This wavefield decomposition method accounts for the discontinuity of variables and is crucial for an accurate interface treatment. We simulate seismic wave propagation with open-pore and sealed-pore boundary conditions and verify the validity and accuracy of the algorithm by comparing the numerical simulations to analytical solutions based on zero viscosity obtained with the Cagniard-de Hoop method. Finally, we illustrate the suitability of our algorithm for more complex models of porous media involving viscous pore fluids and strongly heterogeneous distributions of the elastic and hydraulic material properties.

Le Briançonnais, terrain exotique dans les Alpes ?

New reconstructions of the Western Alps from late Early Jurassic till early Tertiary are proposed. These reconstructions use deep lithospheric data gathered through recent seismic surveys and tomographic studies carried out in the Alps. The present day position, under the Po plain, of the southern limit of the European plate (fig. 1), allows to define the former geometry of the Brianconnais peninsula. The Brianconnais domain is regarded as an exotic terrane formerly belonging to the European margin until Late Jurassic, then transported eastward during the drift of Iberia (fig. 5). Therefore, on a present day Western Alps cross section, a duplication of the European continental margin can be recognized (fig. 10). Stratigraphic and sedimentological data along a zone linking the Pyrenean fracture zone to the Brianconnais, can be related to a rifting event starting in Oxfordian time. This event is responsible for the Late Jurassic till mid-Cretaceous drift of Iberia opening, first the northern Atlantic, then the Gulf of Biscay. Simultaneously, the drift of the Brianconnais will open the Valais ocean and close the Piemontese ocean. The resulting oblique collision zone between the Brianconnais and the Apulian margin generates HP/LT metamorphism starting in Early Cretaceous. The eastward drift of the Brianconnais peninsula will eventually bring it in front of a more northerly segment of the former European margin. The thrusting of the Brianconnais unto that margin takes place in early Tertiary (fig. 9), following the subduction of the Valais ocean. The present nappe pile results not only from continent/continent frontal collision, but also from important lateral displacement of terranes, the most important one being the Brianconnais. The dilemma of `'en echelon'' oceanic domains in the Alps is an outcome of these translations. A solution is found when considering the opening of a Cretaceous Valais ocean across the European margin, running out eastward into the Piemontese ocean, where the drift is taken up along a former transform fault and compensated by subduction under the Apulian margin (fig. 8). In the Western Alps we are then dealing with two oceans, the Piemontese and the Valaisan and a duplicated European margin. In the Eastern Alps the single Piemontese ocean is cut by newly created oceanic crust. All these elements will be incorporated into the Penninic structural domain which does not represent a former unique paleogeographic area, it is a composite accretionary domain squeezed between Europe and Apulia.

Développement de l'interprétation quantitative de la sismique réflexion lacustre : application de l'AVO aux sédiments quaternaires dans le lac Léman

Résumé Des développements antérieurs, au sein de l'Institut de Géophysique de Lausanne, ont permis de développer des techniques d'acquisition sismique et de réaliser l'interprétation des données sismique 2D et 3D pour étudier la géologie de la région et notamment les différentes séquences sédimentaires du Lac Léman. Pour permettre un interprétation quantitative de la sismique en déterminant des paramètres physiques des sédiments la méthode AVO (Amplitude Versus Offset) a été appliquée. Deux campagnes sismiques lacustres, 2D et 3D, ont été acquises afin de tester la méthode AVO dans le Grand Lac sur les deltas des rivières. La géométrie d'acquisition a été repensée afin de pouvoir enregistrer les données à grands déports. Les flûtes sismiques, mises bout à bout, ont permis d'atteindre des angles d'incidence d'environ 40˚ . Des récepteurs GPS spécialement développés à cet effet, et disposés le long de la flûte, ont permis, après post-traitement des données, de déterminer la position de la flûte avec précision (± 0.5 m). L'étalonnage de nos hydrophones, réalisé dans une chambre anéchoïque, a permis de connaître leur réponse en amplitude en fonction de la fréquence. Une variation maximale de 10 dB a été mis en évidence entre les capteurs des flûtes et le signal de référence. Un traitement sismique dont l'amplitude a été conservée a été appliqué sur les données du lac. L'utilisation de l'algorithme en surface en consistante a permis de corriger les variations d'amplitude des tirs du canon à air. Les sections interceptes et gradients obtenues sur les deltas de l'Aubonne et de la Dranse ont permis de produire des cross-plots. Cette représentation permet de classer les anomalies d'amplitude en fonction du type de sédiments et de leur contenu potentiel en gaz. L'un des attributs qui peut être extrait des données 3D, est l'amplitude de la réflectivité d'une interface sismique. Ceci ajoute une composante quantitative à l'interprétation géologique d'une interface. Le fond d'eau sur le delta de l'Aubonne présente des anomalies en amplitude qui caractérisent les chenaux. L'inversion de l'équation de Zoeppritz par l'algorithme de Levenberg-Marquardt a été programmée afin d'extraire les paramètres physiques des sédiments sur ce delta. Une étude statistique des résultats de l'inversion permet de simuler la variation de l'amplitude en fonction du déport. On a obtenu un modèle dont la première couche est l'eau et dont la seconde est une couche pour laquelle V P = 1461 m∕s, ρ = 1182 kg∕m3 et V S = 383 m∕s. Abstract A system to record very high resolution (VHR) seismic data on lakes in 2D and 3D was developed at the Institute of Geophysics, University of Lausanne. Several seismic surveys carried out on Lake Geneva helped us to better understand the geology of the area and to identify sedimentary sequences. However, more sophisticated analysis of the data such as the AVO (Amplitude Versus Offset) method provides means of deciphering the detailed structure of the complex Quaternary sedimentary fill of the Lake Geneva trough. To study the physical parameters we applied the AVO method at some selected places of sediments. These areas are the Aubonne and Dranse River deltas where the configurations of the strata are relatively smooth and the discontinuities between them easy to pick. A specific layout was developed to acquire large incidence angle. 2D and 3D seismic data were acquired with streamers, deployed end to end, providing incidence angle up to 40˚ . One or more GPS antennas attached to the streamer enabled us to calculate individual hydrophone positions with an accuracy of 50 cm after post-processing of the navigation data. To ensure that our system provides correct amplitude information, our streamer sensors were calibrated in an anechoic chamber using a loudspeaker as a source. Amplitude variations between the each hydrophone were of the order of 10 dB. An amplitude correction for each hydrophone was computed and applied before processing. Amplitude preserving processing was then carried out. Intercept vs. gradient cross-plots enable us to determine that both geological discontinuities (lacustrine sediments/moraine and moraine/molasse) have well defined trends. A 3D volume collected on the Aubonne river delta was processed in order ro obtain AVO attributes. Quantitative interpretation using amplitude maps were produced and amplitude maps revealed high reflectivity in channels. Inversion of the water bottom of the Zoeppritz equation using the Levenberg-Marquadt algorithm was carried out to estimate V P , V S and ρ of sediments immediately under the lake bottom. Real-data inversion gave, under the water layer, a mud layer with V P = 1461 m∕s, ρ = 1182 kg∕m3 et V S = 383 m∕s.

Estudio de alta resolución de las fallas activas y procesos sedimentarios asociados al Cañón de São Vicente, margen Sudoeste de Iberia

The Southwest Iberian Margin is caracterized by an intense and diffuse seismic activity due to the convergence between Eurasian and African plates...

Deep incision in an Aptian carbonate succession indicates major sea-level fall in the Cretaceous

Long-term relative sea-level cycles (0 5 to 6 Myr) have yet to be fully understood for the Cretaceous. During the Aptian, in the northern Maestrat Basin (Eastern Iberian Peninsula), fault-controlled subsidence created depositional space, but eustasy governed changes in depositional trends. Relative sea-level history was reconstructed by sequence stratigraphic analysis. Two forced regressive stages of relative sea-level were recognized within three depositional sequences. The first stage is late Early Aptian age (intra Dufrenoyia furcata Zone) and is characterized by foreshore to upper shoreface sedimentary wedges, which occur detached from a highstand carbonate platform, and were deposited above basin marls. The amplitude of relative sea-level drop was in the order of tens of metres, with a duration of <1 Myr. The second stage of relative sea-level fall occurred within the Late Aptian and is recorded by an incised valley that, when restored to its pre-contractional attitude, was >2 km wide and cut 115 m down into the underlying Aptian succession. With the subsequent transgression, the incision was back-filled with peritidal to shallow subtidal deposits. The changes in depositional trends, lithofacies evolution and geometric relation of the stratigraphic units characterized are similar to those observed in coeval rocks within the Maestrat Basin, as well as in other correlative basins elsewhere. The pace and magnitude of the two relative sea-level drops identified fall within the glacio-eustatic domain. In the Maestrat Basin, terrestrial palynological studies provide evidence that the late Early and Late Aptian climate was cooler than the earliest part of the Early Aptian and the Albian Stage, which were characterized by warmer environmental conditions. The outcrops documented here are significant because they preserve the results of Aptian long-term sea-level trends that are often only recognizable on larger scales (i.e. seismic) such as for the Arabian Plate.

Deep incision in an Aptian carbonate succession indicates major sea-level fall in the Cretaceous

Long-term relative sea-level cycles (0 5 to 6 Myr) have yet to be fully understood for the Cretaceous. During the Aptian, in the northern Maestrat Basin (Eastern Iberian Peninsula), fault-controlled subsidence created depositional space, but eustasy governed changes in depositional trends. Relative sea-level history was reconstructed by sequence stratigraphic analysis. Two forced regressive stages of relative sea-level were recognized within three depositional sequences. The first stage is late Early Aptian age (intra Dufrenoyia furcata Zone) and is characterized by foreshore to upper shoreface sedimentary wedges, which occur detached from a highstand carbonate platform, and were deposited above basin marls. The amplitude of relative sea-level drop was in the order of tens of metres, with a duration of <1 Myr. The second stage of relative sea-level fall occurred within the Late Aptian and is recorded by an incised valley that, when restored to its pre-contractional attitude, was >2 km wide and cut 115 m down into the underlying Aptian succession. With the subsequent transgression, the incision was back-filled with peritidal to shallow subtidal deposits. The changes in depositional trends, lithofacies evolution and geometric relation of the stratigraphic units characterized are similar to those observed in coeval rocks within the Maestrat Basin, as well as in other correlative basins elsewhere. The pace and magnitude of the two relative sea-level drops identified fall within the glacio-eustatic domain. In the Maestrat Basin, terrestrial palynological studies provide evidence that the late Early and Late Aptian climate was cooler than the earliest part of the Early Aptian and the Albian Stage, which were characterized by warmer environmental conditions. The outcrops documented here are significant because they preserve the results of Aptian long-term sea-level trends that are often only recognizable on larger scales (i.e. seismic) such as for the Arabian Plate.

3D crosshole ERT for aquifer characterization and monitoring of infiltrating river water

The hydrogeological properties and responses of a productive aquifer in northeastern Switzerland are investigated. For this purpose, 3D crosshole electrical resistivity tomography (ERT) is used to define the main lithological structures within the aquifer (through static inversion) and to monitor the water infiltration from an adjacent river. During precipitation events and subsequent river flooding, the river water resistivity increases. As a consequence, the electrical characteristics of the infiltrating water can be used as a natural tracer to delineate preferential flow paths and flow velocities. The focus is primarily on the experiment installation, data collection strategy, and the structural characterization of the site and a brief overview of the ERT monitoring results. The monitoring system comprises 18 boreholes each equipped with 10 electrodes straddling the entire thickness of the gravel aquifer. A multi-channel resistivity system programmed to cycle through various four-point electrode configurations of the 180 electrodes in a rolling sequence allows for the measurement of approximately 15,500 apparent resistivity values every 7 h on a continuous basis. The 3D static ERT inversion of data acquired under stable hydrological conditions provides a base model for future time-lapse inversion studies and the means to investigate the resolving capability of our acquisition scheme. In particular, it enables definition of the main lithological structures within the aquifer. The final ERT static model delineates a relatively high-resistivity, low-porosity, intermediate-depth layer throughout the investigated aquifer volume that is consistent with results from well logging and seismic and radar tomography models. The next step will be to define and implement an appropriate time-lapse ERT inversion scheme using the river water as a natural tracer. The main challenge will be to separate the superposed time-varying effects of water table height, temperature, and salinity variations associated with the infiltrating water.

The Quaternary Active Faults Database of Iberia (QAFI v.2.0)

The Quaternary Active Faults Database of Iberia (QAFI) is an initiative lead by the Institute of Geology and Mines of Spain (IGME) for building a public repository of scientific data regarding faults having documented activity during the last 2.59 Ma (Quaternary). QAFI also addresses a need to transfer geologic knowledge to practitioners of seismic hazard and risk in Iberia by identifying and characterizing seismogenic fault-sources. QAFI is populated by the information freely provided by more than 40 Earth science researchers, storing to date a total of 262 records. In this article we describe the development and evolution of the database, as well as its internal architecture. Aditionally, a first global analysis of the data is provided with a special focus on length and slip-rate fault parameters. Finally, the database completeness and the internal consistency of the data are discussed. Even though QAFI v.2.0 is the most current resource for calculating fault-related seismic hazard in Iberia, the database is still incomplete and requires further review.

Inferring the mode of origin of polyploid species from next-generation sequence data.

Many eukaryote organisms are polyploid. However, despite their importance, evolutionary inference of polyploid origins and modes of inheritance has been limited by a need for analyses of allele segregation at multiple loci using crosses. The increasing availability of sequence data for nonmodel species now allows the application of established approaches for the analysis of genomic data in polyploids. Here, we ask whether approximate Bayesian computation (ABC), applied to realistic traditional and next-generation sequence data, allows correct inference of the evolutionary and demographic history of polyploids. Using simulations, we evaluate the robustness of evolutionary inference by ABC for tetraploid species as a function of the number of individuals and loci sampled, and the presence or absence of an outgroup. We find that ABC adequately retrieves the recent evolutionary history of polyploid species on the basis of both old and new sequencing technologies. The application of ABC to sequence data from diploid and polyploid species of the plant genus Capsella confirms its utility. Our analysis strongly supports an allopolyploid origin of C. bursa-pastoris about 80 000 years ago. This conclusion runs contrary to previous findings based on the same data set but using an alternative approach and is in agreement with recent findings based on whole-genome sequencing. Our results indicate that ABC is a promising and powerful method for revealing the evolution of polyploid species, without the need to attribute alleles to a homeologous chromosome pair. The approach can readily be extended to more complex scenarios involving higher ploidy levels.

Connexins: key mediators of endocrine function.

The appearance of multicellular organisms imposed the development of several mechanisms for cell-to-cell communication, whereby different types of cells coordinate their function. Some of these mechanisms depend on the intercellular diffusion of signal molecules in the extracellular spaces, whereas others require cell-to-cell contact. Among the latter mechanisms, those provided by the proteins of the connexin family are widespread in most tissues. Connexin signaling is achieved via direct exchanges of cytosolic molecules between adjacent cells at gap junctions, for cell-to-cell coupling, and possibly also involves the formation of membrane "hemi-channels," for the extracellular release of cytosolic signals, direct interactions between connexins and other cell proteins, and coordinated influence on the expression of multiple genes. Connexin signaling appears to be an obligatory attribute of all multicellular exocrine and endocrine glands. Specifically, the experimental evidence we review here points to a direct participation of the Cx36 isoform in the function of the insulin-producing β-cells of the endocrine pancreas, and of the Cx40 isoform in the function of the renin-producing juxtaglomerular epithelioid cells of the kidney cortex.

Electromyographic, cerebral, and muscle hemodynamic responses during intermittent, isometric contractions of the biceps brachii at three submaximal intensities.

This study examined the electromyographic, cerebral and muscle hemodynamic responses during intermittent isometric contractions of biceps brachii at 20, 40, and 60% of maximal voluntary contraction (MVC). Eleven volunteers completed 2 min of intermittent isometric contractions (12/min) at an elbow angle of 90° interspersed with 3 min rest between intensities in systematic order. Surface electromyography (EMG) was recorded from the right biceps brachii and near infrared spectroscopy (NIRS) was used to simultaneously measure left prefrontal and right biceps brachii oxyhemoglobin (HbO2), deoxyhemoglobin (HHb), and total hemoglobin (Hbtot). Transcranial Doppler ultrasound was used to measure middle cerebral artery velocity (MCAv) bilaterally. Finger photoplethysmography was used to record beat-to-beat blood pressure and heart rate. EMG increased with force output from 20 to 60% MVC (P < 0.05). Cerebral HbO2 and Hbtot increased while HHb decreased during contractions with differences observed between 60% vs. 40% and 20% MVC (P < 0.05). Muscle HbO2 decreased while HHb increased during contractions with differences being observed among intensities (P < 0.05). Muscle Hbtot increased from rest at 20% MVC (P < 0.05), while no further change was observed at 40 and 60% MVC (P > 0.05). MCAv increased from rest to exercise but was not different among intensities (P > 0.05). Force output correlated with the root mean square EMG and changes in muscle HbO2 (P < 0.05), but not changes in cerebral HbO2 (P > 0.05) at all three intensities. Force output declined by 8% from the 1st to the 24th contraction only at 60% MVC and was accompanied by systematic increases in RMS, cerebral HbO2 and Hbtot with a leveling off in muscle HbO2 and Hbtot. These changes were independent of alterations in mean arterial pressure. Since cerebral blood flow and oxygenation were elevated at 60% MVC, we attribute the development of fatigue to reduced muscle oxygen availability rather than impaired central neuronal activation.

New directions in Stochastic Multicriteria Acceptability Analysis

Decisions taken in modern organizations are often multi-dimensional, involving multiple decision makers and several criteria measured on different scales. Multiple Criteria Decision Making (MCDM) methods are designed to analyze and to give recommendations in this kind of situations. Among the numerous MCDM methods, two large families of methods are the multi-attribute utility theory based methods and the outranking methods. Traditionally both method families require exact values for technical parameters and criteria measurements, as well as for preferences expressed as weights. Often it is hard, if not impossible, to obtain exact values. Stochastic Multicriteria Acceptability Analysis (SMAA) is a family of methods designed to help in this type of situations where exact values are not available. Different variants of SMAA allow handling all types of MCDM problems. They support defining the model through uncertain, imprecise, or completely missing values. The methods are based on simulation that is applied to obtain descriptive indices characterizing the problem. In this thesis we present new advances in the SMAA methodology. We present and analyze algorithms for the SMAA-2 method and its extension to handle ordinal preferences. We then present an application of SMAA-2 to an area where MCDM models have not been applied before: planning elevator groups for high-rise buildings. Following this, we introduce two new methods to the family: SMAA-TRI that extends ELECTRE TRI for sorting problems with uncertain parameter values, and SMAA-III that extends ELECTRE III in a similar way. An efficient software implementing these two methods has been developed in conjunction with this work, and is briefly presented in this thesis. The thesis is closed with a comprehensive survey of SMAA methodology including a definition of a unified framework.

Fuzzy Principal Components Analysis

Due to the large number of characteristics, there is a need to extract the most relevant characteristicsfrom the input data, so that the amount of information lost in this way is minimal, and the classification realized with the projected data set is relevant with respect to the original data. In order to achieve this feature extraction, different statistical techniques, as well as the principal components analysis (PCA) may be used. This thesis describes an extension of principal components analysis (PCA) allowing the extraction ofa finite number of relevant features from high-dimensional fuzzy data and noisy data. PCA finds linear combinations of the original measurement variables that describe the significant variation in the data. The comparisonof the two proposed methods was produced by using postoperative patient data. Experiment results demonstrate the ability of using the proposed two methods in complex data. Fuzzy PCA was used in the classificationproblem. The classification was applied by using the similarity classifier algorithm where total similarity measures weights are optimized with differential evolution algorithm. This thesis presents the comparison of the classification results based on the obtained data from the fuzzy PCA.

Nappe geometry in the Western Swiss Alps

Detailed geological mapping during the last 20 years in the Western Swiss Alps has shown clearly that most of the lower basement nappes are fold nappes possessing normal and inverted limbs. Moreover their cores are made of strongly deformed gneisses indicating that important ductile strain took place during the formation of the fold nappes. It is therefore probably wrong to imagine deep basement nappes as rigid slices as often actually claimed, especially when interpreting seismic profiles. True `brittle type' thrust nappes involving basement rocks only occur in the internal and upper parts of the belt. Cover nappes, on the contrary, are in most parts of the Alpine belt thrust sheets following more or less the rules of thin-skinned tectonics. Many basement fold nappes lost part of their sedimentary cover during or just before their formation, by decollement along ductile horizons. The result is that many cover thrust nappes in the external part of the Alps are directly related to their original basement fold nappes.