Western Iowa Missouri River Flooding -- Geo-Infrastructure Damage Assessment, Repair and Mitigation Strategies, TR-638, 2013

The 2011 Missouri River flooding caused significant damage to many geo-infrastructure systems including levees, bridge abutments/foundations, paved and unpaved roadways, culverts, and embankment slopes in western Iowa. The flooding resulted in closures of several interchanges along Interstate 29 and of more than 100 miles of secondary roads in western Iowa, causing severe inconvenience to residents and losses to local businesses. The main goals of this research project were to assist county and city engineers by deploying and using advanced technologies to rapidly assess the damage to geo-infrastructure and develop effective repair and mitigation strategies and solutions for use during future flood events in Iowa. The research team visited selected sites in western Iowa to conduct field reconnaissance, in situ testing on bridge abutment backfills that were affected by floods, flooded and non-flooded secondary roadways, and culverts. In situ testing was conducted shortly after the flood waters receded, and several months after flooding to evaluate recovery and performance. Tests included falling weight deflectometer, dynamic cone penetrometer, three-dimensional (3D) laser scanning, ground penetrating radar, and hand auger soil sampling. Field results indicated significant differences in roadway support characteristics between flooded and non-flooded areas. Support characteristics in some flooded areas recovered over time, while others did not. Voids were detected in culvert and bridge abutment backfill materials shortly after flooding and several months after flooding. A catalog of field assessment techniques and 20 potential repair/mitigation solutions are provided in this report. A flow chart relating the damages observed, assessment techniques, and potential repair/mitigation solutions is provided. These options are discussed for paved/unpaved roads, culverts, and bridge abutments, and are applicable for both primary and secondary roadways.

Zonation for 3D aquifer characterization based on joint inversions of multimethod crosshole geophysical data

Predictive groundwater modeling requires accurate information about aquifer characteristics. Geophysical imaging is a powerful tool for delineating aquifer properties at an appropriate scale and resolution, but it suffers from problems of ambiguity. One way to overcome such limitations is to adopt a simultaneous multitechnique inversion strategy. We have developed a methodology for aquifer characterization based on structural joint inversion of multiple geophysical data sets followed by clustering to form zones and subsequent inversion for zonal parameters. Joint inversions based on cross-gradient structural constraints require less restrictive assumptions than, say, applying predefined petro-physical relationships and generally yield superior results. This approach has, for the first time, been applied to three geophysical data types in three dimensions. A classification scheme using maximum likelihood estimation is used to determine the parameters of a Gaussian mixture model that defines zonal geometries from joint-inversion tomograms. The resulting zones are used to estimate representative geophysical parameters of each zone, which are then used for field-scale petrophysical analysis. A synthetic study demonstrated how joint inversion of seismic and radar traveltimes and electrical resistance tomography (ERT) data greatly reduces misclassification of zones (down from 21.3% to 3.7%) and improves the accuracy of retrieved zonal parameters (from 1.8% to 0.3%) compared to individual inversions. We applied our scheme to a data set collected in northeastern Switzerland to delineate lithologic subunits within a gravel aquifer. The inversion models resolve three principal subhorizontal units along with some important 3D heterogeneity. Petro-physical analysis of the zonal parameters indicated approximately 30% variation in porosity within the gravel aquifer and an increasing fraction of finer sediments with depth.

Evaluation of Low Cost Traffic Calming for Rural Communities – Phase II [Updated], TR-630, October, 2013

The main goal of the research described in this report was to evaluate countermeasures that agencies can use to reduce speeds as drivers enter rural communities located on high-speed roadways. The objectives of this study were as follows: * Identify and summarize countermeasures used to manage speeds in transition zones * Demonstrate the effectiveness of countermeasures that are practical for high- to low-speed transition zones * Acquire additional information about countermeasures that may show promise but lack sufficient evidence of effectiveness * Develop an application toolbox to assist small communities in selecting appropriate transition zones and effective countermeasures for entrances to small rural communities The team solicited small communities that were interested in participating in the Phase II study and several communities were also recommended. The treatments evaluated were selected by carefully considering traffic-calming treatments that have been used effectively in other countries for small rural communities, as well as the information gained from the first phase of the project. The treatments evaluated are as follows: * Transverse speed bars * Colored entrance treatment * Temporary island * Radar-activated speed limit sign * Speed feedback sign The toolbox publication and four focused tech briefs also cover the results of this work.

A new method to estimate the infilling of alluvial sediment of glacial valleys using a sloping local base level

A new method is used to estimate the volumes of sediments of glacial valleys. This method is based on the concept of sloping local base level and requires only a digital terrain model and the limits of the alluvial valleys as input data. The bedrock surface of the glacial valley is estimated by a progressive excavation of the digital elevation model (DEM) of the filled valley area. This is performed using an iterative routine that replaces the altitude of a point of the DEM by the mean value of its neighbors minus a fixed value. The result is a curved surface, quadratic in 2D. The bedrock surface of the Rhone Valley in Switzerland was estimated by this method using the free digital terrain model Shuttle Radar Topography Mission (SRTM) (~92 m resolution). The results obtained are in good agreement with the previous estimations based on seismic profiles and gravimetric modeling, with the exceptions of some particular locations. The results from the present method and those from the seismic interpretation are slightly different from the results of the gravimetric data. This discrepancy may result from the presence of large buried landslides in the bottom of the Rhone Valley.

A Nondestructive Method for Determining the Thickness of Sound Concrete on Older Pavements, HR-250, 1983

The primary purpose of this project was to assess the potential of a nondestructive remote sensing system, specifically, ground penetrating subsurface interface radar, for identification and evaluation of D-cracking pavement failures. A secondary purpose was to evaluate the effectiveness of this technique for locating voids under pavements and determining the location of steel reinforcement. From the data collected and the analysis performed to date, the following conclusions can be made regarding the ground penetrating radar system used for this study: (1) steel reinforcement can be accurately located; (2) pavement thickness can be determined; (3) distressed areas in pavements can be located and broadly classified as to severity of deterioration; (4) voids under pavements can be located; and (5) higher resolution recording equipment is required to accurately determine both the thickness of sound pavement remaining over distressed areas and the depth of void areas under pavements.

Geostatistical inversion of seismic and GPR reflection images: what can we actually resolve?

Estimation of the spatial statistics of subsurface velocity heterogeneity from surface-based geophysical reflection survey data is a problem of significant interest in seismic and ground-penetrating radar (GPR) research. A method to effectively address this problem has been recently presented, but our knowledge regarding the resolution of the estimated parameters is still inadequate. Here we examine this issue using an analytical approach that is based on the realistic assumption that the subsurface velocity structure can be characterized as a band-limited scale-invariant medium. Our work importantly confirms recent numerical findings that the inversion of seismic or GPR reflection data for the geostatistical properties of the probed subsurface region is sensitive to the aspect ratio of the velocity heterogeneity and to the decay of its power spectrum, but not to the individual values of the horizontal and vertical correlation lengths.

Radarsat backscattering from an agricultural scene

Orbital remote sensing in the microwave electromagnetic region has been presented as an important tool for agriculture monitoring. The satellite systems in operation have almost all-weather capability and high spatial resolution, which are features appropriated for agriculture. However, for full exploration of these data, an understanding of the relationships between the characteristics of each system and agricultural targets is necessary. This paper describes the behavior of backscattering coefficient (sigma°) derived from calibrated data of Radarsat images from an agricultural area. It is shown that in a dispersion diagram of sigma° there are three main regions in which most of the fields can be classified. The first one is characterized by low backscattering values, with pastures and bare soils; the second one has intermediate backscattering coefficients and comprises well grown crops mainly; and a third one, with high backscattering coefficients, in which there are fields with strong structures causing a kind of double bounce effect. The results of this research indicate that the use of Radarsat images is optimized when a multitemporal analysis is done making the best use of the agricultural calendar and of the dynamics of different cultures.

Joint inversion of crosshole GPR and seismic traveltime data

Joint inversion of crosshole ground-penetrating radar and seismic data can improve model resolution and fidelity of the resultant individual models. Model coupling obtained by minimizing or penalizing some measure of structural dissimilarity between models appears to be the most versatile approach because only weak assumptions about petrophysical relationships are required. Nevertheless, experimental results and petrophysical arguments suggest that when porosity variations are weak in saturated unconsolidated environments, then radar wave speed is approximately linearly related to seismic wave speed. Under such circumstances, model coupling also can be achieved by incorporating cross-covariances in the model regularization. In two case studies, structural similarity is imposed by penalizing models for which the model cross-gradients are nonzero. A first case study demonstrates improvements in model resolution by comparing the resulting models with borehole information, whereas a second case study uses point-spread functions. Although radar seismic wavespeed crossplots are very similar for the two case studies, the models plot in different portions of the graph, suggesting variances in porosity. Both examples display a close, quasilinear relationship between radar seismic wave speed in unconsolidated environments that is described rather well by the corresponding lower Hashin-Shtrikman (HS) bounds. Combining crossplots of the joint inversion models with HS bounds can constrain porosity and pore structure better than individual inversion results can.

Constraining 3-D electrical resistance tomography with GPR reflection data forimproved aquifer characterization

Surface-based ground penetrating radar (GPR) and electrical resistance tomography (ERT) are common tools for aquifer characterization, because both methods provide data that are sensitive to hydrogeologically relevant quantities. To retrieve bulk subsurface properties at high resolution, we suggest incorporating structural information derived from GPR reflection data when inverting surface ERT data. This reduces resolution limitations, which might hinder quantitative interpretations. Surface-based GPR reflection and ERT data have been recorded on an exposed gravel bar within a restored section of a previously channelized river in northeastern Switzerland to characterize an underlying gravel aquifer. The GPR reflection data acquired over an area of 240×40 m map the aquifer's thickness and two internal sub-horizontal regions with different depositional patterns. The interface between these two regions and the boundary of the aquifer with then underlying clay are incorporated in an unstructured ERT mesh. Subsequent inversions are performed without applying smoothness constraints across these boundaries. Inversion models obtained by using these structural constraints contain subtle resistivity variations within the aquifer that are hardly visible in standard inversion models as a result of strong vertical smearing in the latter. In the upper aquifer region, with high GPR coherency and horizontal layering, the resistivity is moderately high (N300 Ωm). We suggest that this region consists of sediments that were rearranged during more than a century of channelized flow. In the lower low coherency region, the GPR image reveals fluvial features (e.g., foresets) and generally more heterogeneous deposits. In this region, the resistivity is lower (~200 Ωm), which we attribute to increased amounts of fines in some of the well-sorted fluvial deposits. We also find elongated conductive anomalies that correspond to the location of river embankments that were removed in 2002.

Inversion of multiple intersecting high-resolution crosshole GPR profiles for hydrological characterization at the Boise Hydrogeophysical Research Site

The integration of geophysical data into the subsurface characterization problem has been shown in many cases to significantly improve hydrological knowledge by providing information at spatial scales and locations that is unattainable using conventional hydrological measurement techniques. In particular, crosshole ground-penetrating radar (GPR) tomography has shown much promise in hydrology because of its ability to provide highly detailed images of subsurface radar wave velocity, which is strongly linked to soil water content. Here, we develop and demonstrate a procedure for inverting together multiple crosshole GPR data sets in order to characterize the spatial distribution of radar wave velocity below the water table at the Boise Hydrogeophysical Research Site (BHRS) near Boise, Idaho, USA. Specifically, we jointly invert 31 intersecting crosshole GPR profiles to obtain a highly resolved and consistent radar velocity model along the various profile directions. The model is found to be strongly correlated with complementary neutron porosity-log data and is further corroborated by larger-scale structural information at the BHRS. This work is an important prerequisite to using crosshole GPR data together with existing hydrological measurements for improved groundwater flow and contaminant transport modeling.

Kernel-based Mapping of Meteorological Fields in Complex Orography

Due to the advances in sensor networks and remote sensing technologies, the acquisition and storage rates of meteorological and climatological data increases every day and ask for novel and efficient processing algorithms. A fundamental problem of data analysis and modeling is the spatial prediction of meteorological variables in complex orography, which serves among others to extended climatological analyses, for the assimilation of data into numerical weather prediction models, for preparing inputs to hydrological models and for real time monitoring and short-term forecasting of weather.In this thesis, a new framework for spatial estimation is proposed by taking advantage of a class of algorithms emerging from the statistical learning theory. Nonparametric kernel-based methods for nonlinear data classification, regression and target detection, known as support vector machines (SVM), are adapted for mapping of meteorological variables in complex orography.With the advent of high resolution digital elevation models, the field of spatial prediction met new horizons. In fact, by exploiting image processing tools along with physical heuristics, an incredible number of terrain features which account for the topographic conditions at multiple spatial scales can be extracted. Such features are highly relevant for the mapping of meteorological variables because they control a considerable part of the spatial variability of meteorological fields in the complex Alpine orography. For instance, patterns of orographic rainfall, wind speed and cold air pools are known to be correlated with particular terrain forms, e.g. convex/concave surfaces and upwind sides of mountain slopes.Kernel-based methods are employed to learn the nonlinear statistical dependence which links the multidimensional space of geographical and topographic explanatory variables to the variable of interest, that is the wind speed as measured at the weather stations or the occurrence of orographic rainfall patterns as extracted from sequences of radar images. Compared to low dimensional models integrating only the geographical coordinates, the proposed framework opens a way to regionalize meteorological variables which are multidimensional in nature and rarely show spatial auto-correlation in the original space making the use of classical geostatistics tangled.The challenges which are explored during the thesis are manifolds. First, the complexity of models is optimized to impose appropriate smoothness properties and reduce the impact of noisy measurements. Secondly, a multiple kernel extension of SVM is considered to select the multiscale features which explain most of the spatial variability of wind speed. Then, SVM target detection methods are implemented to describe the orographic conditions which cause persistent and stationary rainfall patterns. Finally, the optimal splitting of the data is studied to estimate realistic performances and confidence intervals characterizing the uncertainty of predictions.The resulting maps of average wind speeds find applications within renewable resources assessment and opens a route to decrease the temporal scale of analysis to meet hydrological requirements. Furthermore, the maps depicting the susceptibility to orographic rainfall enhancement can be used to improve current radar-based quantitative precipitation estimation and forecasting systems and to generate stochastic ensembles of precipitation fields conditioned upon the orography.

A la recerca del temple d’August a Tarragona: una experiència entre arqueologia i geofísica

Aplicació de diferents tècniques de prospecció geofísica aplicades al subsòl de la Catedral de Tarragona: tomografia de resistivitat elèctrica (ERT), cartografia de conductivitat (EM) i radar de subsòl (GPR).

On-Pavement Signing, 2014

Pavement marking legends are placed on the roadway to remind drivers of the speed limit or to slow down. Use of wording on the pavement surface is more dramatic than use of signing only, which can get lost in the clutter of a streetscape. On-pavement speed limit markings have been used to reinforce speed limits or to indicate a transition zone. For more on this topic by these authors, see also "Evaluation of Dynamic Speed Feedback Signs on Curves: A National Demonstration Project": http://www.trb.org/main/blurbs/172092.aspx