Analysis, modelling and visualisation of spatiotemporal data for urban studies

The proportion of population living in or around cites is more important than ever. Urban sprawl and car dependence have taken over the pedestrian-friendly compact city. Environmental problems like air pollution, land waste or noise, and health problems are the result of this still continuing process. The urban planners have to find solutions to these complex problems, and at the same time insure the economic performance of the city and its surroundings. At the same time, an increasing quantity of socio-economic and environmental data is acquired. In order to get a better understanding of the processes and phenomena taking place in the complex urban environment, these data should be analysed. Numerous methods for modelling and simulating such a system exist and are still under development and can be exploited by the urban geographers for improving our understanding of the urban metabolism. Modern and innovative visualisation techniques help in communicating the results of such models and simulations. This thesis covers several methods for analysis, modelling, simulation and visualisation of problems related to urban geography. The analysis of high dimensional socio-economic data using artificial neural network techniques, especially self-organising maps, is showed using two examples at different scales. The problem of spatiotemporal modelling and data representation is treated and some possible solutions are shown. The simulation of urban dynamics and more specifically the traffic due to commuting to work is illustrated using multi-agent micro-simulation techniques. A section on visualisation methods presents cartograms for transforming the geographic space into a feature space, and the distance circle map, a centre-based map representation particularly useful for urban agglomerations. Some issues on the importance of scale in urban analysis and clustering of urban phenomena are exposed. A new approach on how to define urban areas at different scales is developed, and the link with percolation theory established. Fractal statistics, especially the lacunarity measure, and scale laws are used for characterising urban clusters. In a last section, the population evolution is modelled using a model close to the well-established gravity model. The work covers quite a wide range of methods useful in urban geography. Methods should still be developed further and at the same time find their way into the daily work and decision process of urban planners. La part de personnes vivant dans une région urbaine est plus élevé que jamais et continue à croître. L'étalement urbain et la dépendance automobile ont supplanté la ville compacte adaptée aux piétons. La pollution de l'air, le gaspillage du sol, le bruit, et des problèmes de santé pour les habitants en sont la conséquence. Les urbanistes doivent trouver, ensemble avec toute la société, des solutions à ces problèmes complexes. En même temps, il faut assurer la performance économique de la ville et de sa région. Actuellement, une quantité grandissante de données socio-économiques et environnementales est récoltée. Pour mieux comprendre les processus et phénomènes du système complexe "ville", ces données doivent être traitées et analysées. Des nombreuses méthodes pour modéliser et simuler un tel système existent et sont continuellement en développement. Elles peuvent être exploitées par le géographe urbain pour améliorer sa connaissance du métabolisme urbain. Des techniques modernes et innovatrices de visualisation aident dans la communication des résultats de tels modèles et simulations. Cette thèse décrit plusieurs méthodes permettant d'analyser, de modéliser, de simuler et de visualiser des phénomènes urbains. L'analyse de données socio-économiques à très haute dimension à l'aide de réseaux de neurones artificiels, notamment des cartes auto-organisatrices, est montré à travers deux exemples aux échelles différentes. Le problème de modélisation spatio-temporelle et de représentation des données est discuté et quelques ébauches de solutions esquissées. La simulation de la dynamique urbaine, et plus spécifiquement du trafic automobile engendré par les pendulaires est illustrée à l'aide d'une simulation multi-agents. Une section sur les méthodes de visualisation montre des cartes en anamorphoses permettant de transformer l'espace géographique en espace fonctionnel. Un autre type de carte, les cartes circulaires, est présenté. Ce type de carte est particulièrement utile pour les agglomérations urbaines. Quelques questions liées à l'importance de l'échelle dans l'analyse urbaine sont également discutées. Une nouvelle approche pour définir des clusters urbains à des échelles différentes est développée, et le lien avec la théorie de la percolation est établi. Des statistiques fractales, notamment la lacunarité, sont utilisées pour caractériser ces clusters urbains. L'évolution de la population est modélisée à l'aide d'un modèle proche du modèle gravitaire bien connu. Le travail couvre une large panoplie de méthodes utiles en géographie urbaine. Toutefois, il est toujours nécessaire de développer plus loin ces méthodes et en même temps, elles doivent trouver leur chemin dans la vie quotidienne des urbanistes et planificateurs.

Percolation thresholds in chemical disordered excitable media

The behavior of chemical waves advancing through a disordered excitable medium is investigated in terms of percolation theory and autowave properties in the framework of the light-sensitive Belousov-Zhabotinsky reaction. By controlling the number of sites with a given illumination, different percolation thresholds for propagation are observed, which depend on the relative wave transmittances of the two-state medium considered.

Water percolation estimated with time domain reflectometry (TDR) in drainage lysimeters

Due to the difficulty of estimating water percolation in unsaturated soils, the purpose of this study was to estimate water percolation based on time-domain reflectometry (TDR). In two drainage lysimeters with different soil textures TDR probes were installed, forming a water monitoring system consisting of different numbers of probes. The soils were saturated and covered with plastic to prevent evaporation. Tests of internal drainage were carried out using a TDR 100 unit with constant dielectric readings (every 15 min). To test the consistency of TDR-estimated percolation levels in comparison with the observed leachate levels in the drainage lysimeters, the combined null hypothesis was tested at 5 % probability. A higher number of probes in the water monitoring system resulted in an approximation of the percolation levels estimated from TDR - based moisture data to the levels measured by lysimeters. The definition of the number of probes required for water monitoring to estimate water percolation by TDR depends on the soil physical properties. For sandy clay soils, three batteries with four probes installed at depths of 0.20, 0.40, 0.60, and 0.80 m, at a distance of 0.20, 0.40 and 0.6 m from the center of lysimeters were sufficient to estimate percolation levels equivalent to the observed. In the sandy loam soils, the observed and predicted percolation levels were not equivalent even when using four batteries with four probes each, at depths of 0.20, 0.40, 0.60, and 0.80 m.

Estimation of water percolation by different methods using TDR

Detailed knowledge on water percolation into the soil in irrigated areas is fundamental for solving problems of drainage, pollution and the recharge of underground aquifers. The aim of this study was to evaluate the percolation estimated by time-domain-reflectometry (TDR) in a drainage lysimeter. We used Darcy's law with K(θ) functions determined by field and laboratory methods and by the change in water storage in the soil profile at 16 points of moisture measurement at different time intervals. A sandy clay soil was saturated and covered with plastic sheet to prevent evaporation and an internal drainage trial in a drainage lysimeter was installed. The relationship between the observed and estimated percolation values was evaluated by linear regression analysis. The results suggest that percolation in the field or laboratory can be estimated based on continuous monitoring with TDR, and at short time intervals, of the variations in soil water storage. The precision and accuracy of this approach are similar to those of the lysimeter and it has advantages over the other evaluated methods, of which the most relevant are the possibility of estimating percolation in short time intervals and exemption from the predetermination of soil hydraulic properties such as water retention and hydraulic conductivity. The estimates obtained by the Darcy-Buckingham equation for percolation levels using function K(θ) predicted by the method of Hillel et al. (1972) provided compatible water percolation estimates with those obtained in the lysimeter at time intervals greater than 1 h. The methods of Libardi et al. (1980), Sisson et al. (1980) and van Genuchten (1980) underestimated water percolation.

Magnetic microstructures from magnetic force microscopy and Monte Carlo simulation in CoFe-Ag-Cu granular films

CoFe-Ag-Cu granular films, prepared by rf sputtering, displayed magnetic domain microstructures for ferromagnetic concentrations above about 32% at, and below the percolation threshold. All samples have a fcc structure with an (111) texture perpendicular to the film plane. Magnetic force microscopy (MFM) showed a variety of magnetic domain microstructures, extremely sensitive to the magnetic history of the sample, which arise from the balance of the ferromagnetic exchange, the dipolar interactions and perpendicular magnetocrystalline anisotropy, MFM images indicate that in virgin samples, magnetic bubble domains with an out-of-plane component of the magnetization are surrounded by a quasicontinuous background of opposite magnetization domains. The application of a magnetic field in different geometries drastically modifies the microstructure of the system in the remanent state: i) for an in-plane field, the MFM images show that most of the magnetic moments are aligned along the film plane, ii) for an out-of-plane field, the MFM signal increases about one order of magnitude, and out-of-plane striped domains with alternating up and down magnetization are stabilized. Numerical simulations show that a variety of metastable domain structures (similar to those observed experimentally) can be reached, depending on magnetic history, in systems with competing perpendicular anisotropy, exchange and dipolar interactions.

Hydrodynamic induced deformation of and orientation of a microscopic elastic filament

We describe simulations of an elastic filament immersed in a fluid and subjected to a body force. The coupling between the fluid flow and the friction that the filament experiences induces bending and alignment perpendicular to the force. With increasing force there are four shape regimes, ranging from slight distortion to an unsteady tumbling motion. We also find marginally stable structures. The instability of these shapes and the alignment are explained by induced bending and nonlocal hydrodynamic interactions. These effects are experimentally relevant for stiff microfilaments.

Generalized percolation in random directed networks

We develop a general theory for percolation in directed random networks with arbitrary two-point correlations and bidirectional edgesthat is, edges pointing in both directions simultaneously. These two ingredients alter the previously known scenario and open new views and perspectives on percolation phenomena. Equations for the percolation threshold and the sizes of the giant components are derived in the most general case. We also present simulation results for a particular example of uncorrelated network with bidirectional edges confirming the theoretical predictions.

Clustering in complex networks. II. Percolation properties

The percolation properties of clustered networks are analyzed in detail. In the case of weak clustering, we present an analytical approach that allows us to find the critical threshold and the size of the giant component. Numerical simulations confirm the accuracy of our results. In more general terms, we show that weak clustering hinders the onset of the giant component whereas strong clustering favors its appearance. This is a direct consequence of the differences in the k-core structure of the networks, which are found to be totally different depending on the level of clustering. An empirical analysis of a real social network confirms our predictions.

Effects of Pavement Surface Texture on Noise and Frictional Characteristics HR-281, February 1987

An experimental modification of the transverse groove surface texture of a section of an urban interstate highway was performed by the Iowa Department of Transportation. Transverse groove texturing i s a design feature required by the Federal Highway Administration t o reduce skidding under wet pavement conditions. Adjacent residents claimed the texturing was the cause of especially annoying tonal characteristics within the traffic noise. A research proposal to modify the existing texture pattern by surface grinding and to study the noise and friction effects was approved for funding by the Iowa Highway Research Board. Results i n the form of a comparison between traffic noise before modification and traffic noise immediately after and 15 months after modification indicate that the change in surface texture has lowered overall traffic noise levels by reducing a high frequency component of the traffic noise spectrum. Fraffic testing data show reduced capacity of the roadway to inhibit wet pavement skidding as a result of the surface modification.

Percolation thresholds in chemical disordered excitable media

The behavior of chemical waves advancing through a disordered excitable medium is investigated in terms of percolation theory and autowave properties in the framework of the light-sensitive Belousov-Zhabotinsky reaction. By controlling the number of sites with a given illumination, different percolation thresholds for propagation are observed, which depend on the relative wave transmittances of the two-state medium considered.

The Potential of Friction as a Tool for Winter Maintenance, TR-400, 1998

It is intuitively obvious that snow or ice on a road surface will make that surface more slippery and thus more hazardous. However, quantifying this slipperiness by measuring the friction between the road surface and a vehicle is rather difficult. If such friction readings could be easily made, they might provide a means to control winter maintenance activities more efficiently than at present. This study is a preliminary examination of the possibility of using friction as an operational tool in winter maintenance. In particular, the relationship of friction to traffic volume and speed, and accident rates is examined, and the current lack of knowledge in this area is outlined. The state of the art of friction measuring techniques is reviewed. A series of experiments whereby greater knowledge of how friction deteriorates during a storm and is restored by treatment is proposed. The relationship between plowing forces and the ice-pavement bond strength is discussed. The challenge of integrating all these potential sources of information into a useful final product is presented together with a potential approach. A preliminary cost-benefit analysis of friction measuring devices is performed and suggests that considerable savings might be realized if certain assumptions should hold true. The steps required to bring friction from its current state as a research tool to full deployment as an operational tool are presented and discussed. While much remains to be done in this regard, it is apparent that friction could be an extremely effective operational tool in winter maintenance activities of the future.

Pile Design and Tests for Integral Abutment Bridges, HR-273, 1987

Expansion joints increase both the initial cost and the maintenance cost of bridges. Integral abutment bridges provide an attractive design alternative because expansion joints are eliminated from the bridge itself. However, the piles in these bridges are subjected to horizontal movement as the bridge expands and contracts during temperature changes. The objective of this research was to develop a method of designing piles for these conditions. Separate field tests simulating a pile and a bridge girder were conducted for three loading cases: (1) vertical load only, (2) horizontal displacement of pile head only, and (3) combined horizontal displacement of pile head with subsequent vertical load. Both tests (1) and (3) reached the same ultimate vertical load, that is, the horizontal displacement had no effect on the vertical load capacity. Several model tests were conducted in sand with a scale factor of about 1:10. Experimental results from both the field and model tests were used to develop the vertical and horizontal load-displacement properties of the soil. These properties were input into the finite element computer program Integral Abutment Bridge Two-Dimensional (IAB2D), which was developed under a previous research contract. Experimental and analytical results compared well for the test cases. Two alternative design methods, both based upon the American Association of State Highway and Transportation Officials (AASHTO) Specification, were developed. Alternative One is quite conservative relative to IAB2D results and does not permit plastic redistribution of forces. Alternative Two is also conservative when compared to IAB2D, but plastic redistribution is permitted. To use Alternative Two, the pile cross section must have sufficient inelastic rotation capacity before local buckling occurs. A design example for a friction pile and an end-bearing pile illustrates both alternatives.

Social network architecture and the maintenance of deleterious cultural traits.

How have changes in communications technology affected the way that misinformation spreads through a population and persists? To what extent do differences in the architecture of social networks affect the spread of misinformation, relative to the rates and rules by which individuals transmit or eliminate different pieces of information (cultural traits)? Here, we use analytical models and individual-based simulations to study how a 'cultural load' of misinformation can be maintained in a population under a balance between social transmission and selective elimination of cultural traits with low intrinsic value. While considerable research has explored how network architecture affects percolation processes, we find that the relative rates at which individuals transmit or eliminate traits can have much more profound impacts on the cultural load than differences in network architecture. In particular, the cultural load is insensitive to correlations between an individual's network degree and rate of elimination when these quantities vary among individuals. Taken together, these results suggest that changes in communications technology may have influenced cultural evolution more strongly through changes in the amount of information flow, rather than the details of who is connected to whom.

Regional Approach to Landslide Interpretation and Repair, TR-430, 2001

The objective of this report is to provide Iowa county engineers and highway maintenance personnel with procedures that will allow them to efficiently and effectively interpret and repair or avoid landslides. The research provides an overview of basic slope stability analyses that can be used to diagnose the cause and effect associated with a slope failure. Field evidence for identifying active or potential slope stability problems is outlined. A survey of county engineers provided data for presenting a slope stability risk map for the state of Iowa. Areas of high risk are along the western border and southeastern portion of the state. These regions contain deep to moderately deep loess. The central portion of the state is a low risk area where the surficial soils are glacial till or thin loess over till. In this region, the landslides appear to occur predominately in backslopes along deeply incised major rivers, such as the Des Moines River, or in foreslopes. The south-central portion of the state is an area of medium risk where failures are associated with steep backslopes and improperly compacted foreslopes. Soil shear strength data compiled from the Iowa DOT and consulting engineers files are correlated with geologic parent materials and mean values of shear strength parameters and unit weights were computed for glacial till, friable loess, plastic loess and local alluvium. Statistical tests demonstrate that friction angles and unit weights differ significantly but in some cases effective stress cohesion intercept and undrained shear strength data do not. Moreover, effective stress cohesion intercept and undrained shear strength data show a high degree of variability. The shear strength and unit weight data are used in slope stability analyses for both drained and undrained conditions to generate curves that can be used for a preliminary evaluation of the relative stability of slopes within the four materials. Reconnaissance trips to over fifty active and repaired landslides in Iowa suggest that, in general, landslides in Iowa are relatively shallow [i.e., failure surfaces less than 6 ft (2 m) deep] and are either translational or shallow rational. Two foreslope and two backslope failure case histories provide additional insights into slope stability problems and repair in Iowa. These include the observation that embankment soils compacted to less than 95% relative density show a marked strength decrease from soils at or above that density. Foreslopes constructed of soils derived from shale exhibit loss of strength as a result of weathering. In some situations, multiple causes of instability can be discerned from back analyses with the slope stability program XSTABL. In areas where the stratigraphy consists of loess over till or till over bedrock, the geologic contracts act as surfaces of groundwater accumulation that contribute to slope instability.

Longitudinal Grinding and Transverse Grooving to Improve Frictional and Profile Properties, HR-1035, 1983

There are many miles of portland cement concrete pavement in Iowa that due to normal wear, and in some cases accelerated wear from studded tires, the surface has become polished resulting in less than desirable friction values. Retexturing the surface may be an economical way to re-establish desirable friction values. Retexturing by grinding with diamond blades and transverse grooving with diamond blades are two methods of rehabilitating p.c.c. surfaces. MU Inc. of Lebanan, Tennessee proposed to provide without charge to the Iowa Department of Transportation, one 1500 ft x 12 ft section each of three methods of texturing. They are longitudinal grinding, transverse grooving and longitudinal grinding followed by transverse grooving. A section of 1500 feet is needed to properly evaluate a texturing method. It was decided by Iowa DOT personnel that due to possible differential friction it would be undesirable to texture only one lane. The decision was made to do test sections of 1500 ft x 24 ft with the cost of the additional texturing paid by the Iowa DOT. Iowa also has areas where the p.c.c. pavement has faulted at the joints and cracks which results in poor riding quality. Methods of correcting the faulting are to underseal the pavement where needed and/or grinding the surface to eliminate the faulted areas. It was decided to include in this research project a section for profiling by grinding.