Palaeosecular variation over the past 5 Myr based on a new generalized database

A generalized database of global palaeomagnetic data from 3719 lava flows and thin dykes of age 0–5 Ma has been constructed for use with a relational database management system. The database includes all data whose virtual geomagnetic poles (VGP) lie within 45 of the spin axis and can be used for studies of palaeosecular variation and for geomagnetic field modelling. Because many of these data were collected and processed more than 15–20 years ago, each result has been characterized according to the demagnetization procedures carried out. Analysis of these data in terms of the latitude variation of the angular dispersion of VGPs (palaeosecular variation from lavas) strongly suggests that careful data selection is required and that many of the older studies may need to be redone using more modern methods. Differences between the angular dispersions for separate normal- and reverse-polarity data sets confirm that many older studies have not been adequately cleaned magnetically. Therefore, the use of the database for geomagnetic field modelling should be carried out with some caution. Using a VGP cut-off angle that varies with latitude, the best data set consists of 2636 results that show a smooth increase of VGP angular dispersion with latitude. Model G for palaeosecular variation, which is based on modelling of the antisymmetric (dipole) and symmetric (quadrupole) dynamo families, provides a good fit to these results.     

High-resolution climate simulation of the last glacial maximum

The climate of the last glacial maximum (LGM) is simulated with a high-resolution atmospheric general circulation model, the NCAR CCM3 at spectral truncation of T170, corresponding to a grid cell size of roughly 75 km. The purpose of the study is to assess whether there are significant benefits from the higher resolution simulation compared to the lower resolution simulation associated with the role of topography. The LGM simulations were forced with modified CLIMAP sea ice distribution and sea surface temperatures (SST) reduced by 1°C, ice sheet topography, reduced CO₂, and 21,000 BP orbital parameters. The high-resolution model captures modern climate reasonably well, in particular the distribution of heavy precipitation in the tropical Pacific. For the ice age case, surface temperature simulated by the high-resolution model agrees better with those of proxy estimates than does the low-resolution model. Despite the fact that tropical SSTs were only 2.1°C less than the control run, there are many lowland tropical land areas 4–6°C colder than present. Comparison of T170 model results with the best constrained proxy temperature estimates (noble gas concentrations in groundwater) now yield no significant differences between model and observations. There are also significant upland temperature changes in the best resolved tropical mountain belt (the Andes). We provisionally attribute this result in part as resulting from decreased lateral mixing between ocean and land in a model with more model grid cells. A longstanding model-data discrepancy therefore appears to be resolved without invoking any unusual model physics. The response of the Asian summer monsoon can also be more clearly linked to local geography in the high-resolution model than in the low-resolution model; this distinction should enable more confident validation of climate proxy data with the high-resolution model. Elsewhere, an inferred salinity increase in the subtropical North Atlantic may have significant implications for ocean circulation changes during the LGM. A large part of the Amazon and Congo Basins are simulated to be substantially drier in the ice age—consistent with many (but not all) paleo data. These results suggest that there are considerable benefits derived from high-resolution model regarding regional climate responses, and that observationalists can now compare their results with models that resolve geography at a resolution comparable to that which the proxy data represent.     

A time-distribution study of doline collapse: Framework for prediction

A temporal analysis of doline collapse on the Western Highland Rim (Tennessee) indicated one approach to estimating the probability of collapse in areas where the geologic and hydrologic criteria associated with collapse have been identified. The distribution of collapse was examined for trend, autocorrelation, and goodness of fit. The distribution of doline collapse during one 12-month period conformed to a Poisson distribution with a mean occurrence rate of λ=0.346 collapses per week and with the interoccurrence times being exponentially distributed (0.01 level). Although the proposed model is spatially and temporally restricted, it may provide an initial framework for estimating the probability of doline collapse in other karst terrains of similar geologic and hydrologic settings.     

The karst contagion model: Synopsis and environmental implications

The contagion model of karst terrane evolution focuses on the environmental implications for a large karst depression population on the Pennyroyal Plain (southern Kentucky) and the adjacent Western Highland Rim (Tennessee) immediately south of the Mammoth Cave Plateau. In karst terranes where the contagion model applies, there is a well-defined infrastructure comprised of hydrologic, structural geologic and geomorphic interacting elements that result in clustered depressions underlain by a radial conduit system. Clusters tend to be randomly distributed and typically contain a parent depression surrounded by numerous daughters.Groundwater flow is assumed to be turbulent and confined largely to conduits that are 3-dimensionally configured between clusters in a dendritic to trellis network. Parent depressions serve as conduit nodes for collecting groundwater migrating from beneath daughter depressions.Flow velocities in the 3-dimensional cluster-cell conduits exceed those in granular media by several orders of magnitude making pathogen and chemical contaminant migration rapid. Groundwater quality assessment in karst conduit hydrogeologic settings is difficult because monitoring wells are inappropriate. Monitoring wells may have a low probability of intercepting a major conduit and therefore the sampling regime must take into consideration the pulse discharge of pollutants in karst conduits. Representative water quality data must come from springs located near the local base level.     

Deep crustal structure and rheology of the Gascoyne volcanic margin,western Australia

Interpretation of deep seismic reflection data across the Gascoyne Margin reveals six distinct seismic facies units related to the tectono-magmatic breakup history. On the outer Exmouth Plateau four large scale units are identified: (1) an extensively block-faulted upper crust; (2) a middle-crustal unit of discontinuous, undulatory reflectors; (3) a reflection-free deep crustal unit; and (4) a lower-crustal band of low-frequency, high-amplitude reflectors. Two additional units are found near the continent-ocean boundary (COB); (5) seaward-dipping reflectors (SDR); and (6) landward-dipping reflectors in the lower crust below the SDR. The lower-crustal high-reflectivity band, located near the top of a high-velocity unit (V_p > 7 kms^–1), is interpreted as magmatic underplating. There is a spatial correlation between the underplated area and the presence of extensive upper-crustal block-faulting and intrusive rocks in the shallow crust. The undulatory middle-crustal reflector unit is also only identified in the outer plateau area, and is interpreted as a zone in which the upper-crustal faults terminate. The inner parts of the margin consist of a deep basin showing little upper-crustal faulting and no evidence of middle crustal deformation or underplating. Theoretical modeling of the effect of rifting and magmatic underplating on crustal strength profiles suggests that the brittle-ductile transition may migrate at least 5 km upwards during several million years after the underplating event. Based on the seismic interpretation and crustal strength modeling we propose that the seismic structure of the outer Exmouth Plateau is severely modified by a transient change in the crustal rheological structure associated with magmatic underplating.     

Assessing spatial and attribute errors in large national datasets for population distribution models: a case study of Philadelphia county schools

Geospatial technologies and digital data have developed and disseminated rapidly in conjunction with increasing computing efficiency and Internet availability. The ability to store and transmit large datasets has encouraged the development of national infrastructure datasets in geospatial formats. National datasets are used by numerous agencies for analysis and modeling purposes because these datasets are standardized and considered to be of acceptable accuracy for national scale applications. At Oak Ridge National Laboratory a population model has been developed that incorporates national schools data as one of the model inputs. This paper evaluates spatial and attribute inaccuracies present within two national school datasets, Tele Atlas North America and National Center of Education Statistics (NCES). Schools are an important component of the population model, because they are spatially dense clusters of vulnerable populations. It is therefore essential to validate the quality of school input data. Schools were also chosen since a validated schools dataset was produced in geospatial format for Philadelphia County; thereby enabling a comparison between a local dataset and the national datasets. Analyses found the national datasets are not standardized and incomplete, containing 76 to 90 percent of existing schools. The temporal accuracy of updating annual enrollment values resulted in 89 percent inaccuracy for 2003. Spatial rectification was required for 87 percent of NCES points, of which 58 percent of the errors were attributed to the geocoding process. Lastly, it was found that by combining the two national datasets, the resultant dataset provided a more useful and accurate solution.     

Stresses in column-supported hyperboloidal shells subject to seismic loading

The computation of stresses within a finite element displacement method analysis of a shell of revolution is considered. The common procedure of applying the kinematic and constitutive laws to the displacement functions is examined and justified for models where the displacements are represented by high-order polynomial expansions. Also, two alternative computational formats within this technique are explored. The influence of the column-supported base condition on a hyperboloidal shell of revolution is studied with respect to the stresses calculated from a response spectrum analysis. These studies emphasize the importance of accurately modelling the base region of a column-supported shell such as a hyperbolic cooling tower.