Geology topics

J.C. Davis

Publications and source records attributed to J.C. Davis.

At least 19 recordsLinked to original sources

USGS2006

Two models for evaluating landslide hazards

Two alternative procedures for estimating landslide hazards were evaluated using data on topographic digital elevation models (DEMs) and bedrock lithologies in an area adjacent to the Missouri River in Atchison County, Kansas, USA. The two procedures are based on the likelihood ratio model but utilize different assumptions. The empirical likelihood ratio model is based on non-parametric empirical univariate frequency distribution functions under an assumption of conditional independence while the multivariate logistic discriminant model assumes that likelihood ratios can be expressed in terms of logistic functions. The relative hazards of occurrence of landslides were estimated by an empirical likelihood ratio model and by multivariate logistic discriminant analysis. Predictor variables consisted of grids containing topographic elevations, slope angles, and slope aspects calculated from a 30-m DEM. An integer grid of coded bedrock lithologies taken from digitized geologic maps was also used as a predictor variable. Both statistical models yield relative estimates in the form of the proportion of total map area predicted to already contain or to be the site of future landslides. The stabilities of estimates were checked by cross-validation of results from random subsamples, using each of the two procedures. Cell-by-cell comparisons of hazard maps made by the two models show that the two sets of estimates are virtually identical. This suggests that the empirical likelihood ratio and the logistic discriminant analysis models are robust with respect to the conditional independent assumption and the logistic function assumption, respectively, and that either model can be used successfully to evaluate landslide hazards. ?? 2006.

Computers & Geosciences↗

USGS2003

Using multiple logistic regression and GIS technology to predict landslide hazard in northeast Kansas, USA

Landslides in the hilly terrain along the Kansas and Missouri rivers in northeastern Kansas have caused millions of dollars in property damage during the last decade. To address this problem, a statistical method called multiple logistic regression has been used to create a landslide-hazard map for Atchison, Kansas, and surrounding areas. Data included digitized geology, slopes, and landslides, manipulated using ArcView GIS. Logistic regression relates predictor variables to the occurrence or nonoccurrence of landslides within geographic cells and uses the relationship to produce a map showing the probability of future landslides, given local slopes and geologic units. Results indicated that slope is the most important variable for estimating landslide hazard in the study area. Geologic units consisting mostly of shale, siltstone, and sandstone were most susceptible to landslides. Soil type and aspect ratio were considered but excluded from the final analysis because these variables did not significantly add to the predictive power of the logistic regression. Soil types were highly correlated with the geologic units, and no significant relationships existed between landslides and slope aspect. ?? 2003 Elsevier Science B.V. All rights reserved.

Engineering Geology↗

USGS2000

Assessing the Kansas water-level monitoring program: An example of the application of classical statistics to a geological problem

Geologists may feel that geological data are not amenable to statistical analysis, or at best require specialized approaches such as nonparametric statistics and geostatistics. However, there are many circumstances, particularly in systematic studies conducted for environmental or regulatory purposes, where traditional parametric statistical procedures can be beneficial. An example is the application of analysis of variance to data collected in an annual program of measuring groundwater levels in Kansas. Influences such as well conditions, operator effects, and use of the water can be assessed and wells that yield less reliable measurements can be identified. Such statistical studies have resulted in yearly improvements in the quality and reliability of the collected hydrologic data. Similar benefits may be achieved in other geological studies by the appropriate use of classical statistical tools.

The Compass: Earth Science Journal of Sigma Gamma ↗

USGS2000

Determination of geochemical affinities of granitic rocks from the Aue-Schwarzenberg zone (Erzgebirge, Germany) by multivariate statistics

Variscan granites of the Erzgebirge region can be effectively classified into five genetically distinct major groups by canonical analysis of geochemical variables. The same classification procedure, when applied to small plutons in the Aue-Schwarzenberg granite zone (ASGZ), shows that all ASGZ granites have compositional affinities to low-F biotite or low-F two-mica granite groups. This suggests that the ASGZ granites were emplaced during the first, late-collisional stage of silicic magmatism in the region, which occurred between about 325 and 318 Ma. The numerous biotite granite bodies in the zone are geochemically distinct from both the neighboring Kirchberg granite pluton and the spatially displaced Niederbobritzsch biotite granite massif. Instead, these bodies seem to constitute a third sub-group within the low-F biotite granite class. The ASGZ biotite granites represent three or more genetically distinct bodies, thus highlighting the enormous compositional variability within this group of granites. Least evolved samples of two-mica granites from the ASGZ apparently reflect the assimilation of low-grade metamorphic country rocks during emplacement, altering the original composition of the melts by enhancing primary Al content. The same genesis is implied for the rare "cordierite granite" facies of the Bergen massif, the type pluton for the low-F two-mica granite group in the Erzgebirge.

Neues Jahrbuch fur Mineralogie, Monatshefte↗

USGS1999

Multivariate analyses of Erzgebirge granite and rhyolite composition: Implications for classification of granites and their genetic relations

High-precision major, minor and trace element analyses for 44 elements have been made of 329 Late Variscan granitic and rhyolitic rocks from the Erzgebirge metallogenic province of Germany. The intrusive histories of some of these granites are not completely understood and exposures of rock are not adequate to resolve relationships between what apparently are different plutons. Therefore, it is necessary to turn to chemical analyses to decipher the evolution of the plutons and their relationships. A new classification of Erzgebirge plutons into five major groups of granites, based on petrologic interpretations of geochemical and mineralogical relationships (low-F biotite granites; low-F two-mica granites; high-F, high-P2O5 Li-mica granites; high-F, low-P2O5 Li-mica granites; high-F, low-P2O5 biotite granites) was tested by multivariate techniques. Canonical analyses of major elements, minor elements, trace elements and ratio variables all distinguish the groups with differing amounts of success. Univariate ANOVA's, in combination with forward-stepwise and backward-elimination canonical analyses, were used to select ten variables which were most effective in distinguishing groups. In a biplot, groups form distinct clusters roughly arranged along a quadratic path. Within groups, individual plutons tend to be arranged in patterns possibly reflecting granitic evolution. Canonical functions were used to classify samples of rhyolites of unknown association into the five groups. Another canonical analysis was based on ten elements traditionally used in petrology and which were important in the new classification of granites. Their biplot pattern is similar to that from statistically chosen variables but less effective at distinguishing the five groups of granites. This study shows that multivariate statistical techniques can provide significant insight into problems of granitic petrogenesis and may be superior to conventional procedures for petrological interpretation.

Computers & Geosciences↗

USGS1999

Subdivision of Holocene Baltic sea sediments by their physical properties [Gliederung holozaner ostseesedimente nach physikalischen Eigenschaften]

The Holocene sediment sequence of a core taken within the centre of the Eastern Gotland Basin was subdivided into 12 lithostratigraphic units based on MSCL-data (sound velocity, wet bulk density, magnetic susceptibility) using a multivariate classification method. The lower 6 units embrace the sediments until the Litorina transgression, and the upper 6 units subdivide the brackish-marine Litorina- and post-Litorina sediments. The upper lithostratigraphic units reflect a change of anoxic (laminated) and oxic (non-laminated) sediments. By application of a numerical stratigraphic correlation method the zonation was extended laterally onto contiguous sediment cores within the central basin. Consequently the change of anoxic and oxic sediments can be used for a general lithostratigraphic subdivision of sediments of the Gotland Basin. A quantitative criterion based on the sediment-physical lithofacies is added to existing subdivisions of the Holocene in the Baltic Sea.

Petermanns Geographische Mitteilungen↗

USGS1998

Singularity and Nonnormality in the Classification of Compositional Data

Geologists may want to classify compositional data and express the classification as a map. Regionalized classification is a tool that can be used for this purpose, but it incorporates discriminant analysis, which requires the computation and inversion of a covariance matrix. Covariance matrices of compositional data always will be singular (noninvertible) because of the unit-sum constraint. Fortunately, discriminant analyses can be calculated using a pseudo-inverse of the singular covariance matrix; this is done automatically by some statistical packages such as SAS. Granulometric data from the Darss Sill region of the Baltic Sea is used to explore how the pseudo-inversion procedure influences discriminant analysis results, comparing the algorithm used by SAS to the more conventional Moore-Penrose algorithm. Logratio transforms have been recommended to overcome problems associated with analysis of compositional data, including singularity. A regionalized classification of the Darss Sill data after logratio transformation is different only slightly from one based on raw granulometric data, suggesting that closure problems do not influence severely regionalized classification of compositional data.

Mathematical Geology↗

USGS1997

Geochemical patterns in soils of the karst region, Croatia

Soil samples were collected at 420 locations in a 5-km grid pattern in the Istria and Gorski Kotar areas of Croatia, and on the Croatian islands of Cres, Rab and Krk, in order to relate geochemical variation in the soils to underlying differences in geology, bedrock lithology, soil type, environment and natural versus anthropogenic influences. Specific objectives included assessment of possible agricultural and industrial sources of contamination, especially from airborne effluent emitted by a local power plant. The study also tested the adequacy of a fixed-depth soil sampling procedure developed for meager karstic soils. Although 40 geochemical variables were analyzed, only 15 elements and 5 radionuclides are common to all the sample locations. These elements can be divided into three groups: (1) those of mostly anthropogenic origin — Pb, V, Cu and Cr; (2) those of mixed origin — radionuclides and Zn; and (3) those of mostly geogene origin — Ba, Sr, Ti, Al, Na, Ca, Mg, Fe, Mn, Ni and Co. Variation in Pb shows a strong correlation with the pattern of road traffic in Istria. The distributions of Ca, Na and Mg in the flysch basins of southern Istria and Slovenia are clearly distinguishable from the distributions of these elements in the surrounding carbonate terrains, a consequence of differences in bedrock permeability, type of drainage and pH. The spatial pattern of 137 Cs from the Chernobyl nuclear power plant accident reflects almost exclusively the precipitation in Istria during the days immediately after the explosion.

Journal of Geochemical Exploration↗

USGS1997

Spatial analysis of temperature (BHT/DST) data and consequences for heat-flow determination in sedimentary basins

Large numbers of bottom-hole temperatures (BHTs) and temperatures measured during drill-stem tests (DSTs) are available in areas explored for hydrocarbons, but their usefulness for estimating geothermal gradients and heat-flow density is limited. We investigated a large data set of BHT and DST measurements taken in boreholes in the American Midcontinent, a geologically uniform stable cratonic area, and propose an empirical correction for BHTs based on relationships between BHTs, DSTs, and thermal logs. This empirical correction is compared with similar approaches determined for other areas. The data were analyzed by multivariate statistics prior to the BHT correction to identify anomalous measurements and quantify external influences. Spatial patterns in temperature measurements for major stratigraphic units outline relations to regional structure. Comparision of temperature and structure trend-surface residuals reveals a relationship between temperature highs and local structure highs. The anticlines, developed by continuous but intermittent movement of basement fault blocks in the Late Paleozoic, are subtle features having closures of 10-30 m and contain relatively small hydrocarbon reservoirs. The temperature anomalies of the order of 5-7 ??C may reflect fluids moving upward along fractures and faults, rather than changes in thermal conductivity resulting from different pore fluids. ?? Springer-Verlag 1997.

Conference Paper↗

USGS1995

Adjusting stream-sediment geochemical maps in the Austrian Bohemian Massif by analysis of variance

The Austrian portion of the Bohemian Massif is a Precambrian terrane composed mostly of highly metamorphosed rocks intruded by a series of granitoids that are petrographically similar. Rocks are exposed poorly and the subtle variations in rock type are difficult to map in the field. A detailed geochemical survey of stream sediments in this region has been conducted and included as part of the Geochemischer Atlas der Republik O??sterreich, and the variations in stream sediment composition may help refine the geological interpretation. In an earlier study, multivariate analysis of variance (MANOVA) was applied to the stream-sediment data in order to minimize unwanted sampling variation and emphasize relationships between stream sediments and rock types in sample catchment areas. The estimated coefficients were used successfully to correct for the sampling effects throughout most of the region, but also introduced an overcorrection in some areas that seems to result from consistent but subtle differences in composition of specific rock types. By expanding the model to include an additional factor reflecting the presence of a major tectonic unit, the Rohrbach block, the overcorrection is removed. This iterative process simultaneously refines both the geochemical map by removing extraneous variation and the geological map by suggesting a more detailed classification of rock types. ?? 1995 International Association for Mathematical Geology.

Mathematical Geology↗

USGS1995

Estimation of regionalized compositions: A comparison of three methods

A regionalized composition is a random vector function whose components are positive and sum to a constant at every point of the sampling region. Consequently, the components of a regionalized composition are necessarily spatially correlated. This spatial dependence-induced by the constant sum constraint-is a spurious spatial correlation and may lead to misinterpretations of statistical analyses. Furthermore, the cross-covariance matrices of the regionalized composition are singular, as is the coefficient matrix of the cokriging system of equations. Three methods of performing estimation or prediction of a regionalized composition at unsampled points are discussed: (1) the direct approach of estimating each variable separately; (2) the basis method, which is applicable only when a random function is available that can he regarded as the size of the regionalized composition under study; (3) the logratio approach, using the additive-log-ratio transformation proposed by J. Aitchison, which allows statistical analysis of compositional data. We present a brief theoretical review of these three methods and compare them using compositional data from the Lyons West Oil Field in Kansas (USA). It is shown that, although there are no important numerical differences, the direct approach leads to invalid results, whereas the basis method and the additive-log-ratio approach are comparable. ?? 1995 International Association for Mathematical Geology.

Mathematical Geology↗

USGS1993

Prediction of hydrocarbons in sedimentary basins

To estimate the undiscovered hydrocarbon potential of sedimentary basins, quantitative play assessments specific for each location in a region may be obtained using geostatistical methods combined with the theory of classification of geological objects, a methodology referred to as regionalization. The technique relies on process modeling and measured borehole data as well as probabilistic methods to exploit the relationship between geology (the "predictor") and known hydrocarbon productivity (the "target") to define prospective stratigraphic intervals within a basin. It is demonstrated in case studies from the oil-producing region of the western Kansas Pennsylvanian Shelf and the gas-bearing Rotliegend sediments of the Northeast German Basin. ?? 1993 International Association for Mathematical Geology.

Mathematical Geology↗

USGS1993

Boundary assessment under uncertainty: A case study

Estimating certain attributes within a geological body whose exact boundary is not known presents problems because of the lack of information. Estimation may result in values that are inadmissible from a geological point of view, especially with attributes which necessarily must be zero outside the boundary, such as the thickness of the oil column outside a reservoir. A simple but effective way to define the boundary is to use indicator kriging in two steps, the first for the purpose of extrapolating control points outside the body, the second to obtain a weighting function which expresses the uncertainty attached to estimations obtained in the boundary region. ?? 1993 International Association for Mathematical Geology.

Mathematical Geology↗

USGS1993

Volumetric calculations in an oil field: The basis method

The basis method for estimating oil reserves in place is compared to a traditional procedure that uses ordinary kriging. In the basis method, auxiliary variables that sum to the net thickness of pay are estimated by cokriging. In theory, the procedure should be more powerful because it makes full use of the cross-correlation between variables and forces the original variables to honor interval constraints. However, at least in our case study, the practical advantages of cokriging for estimating oil in place are marginal. ?? 1993.

Computers & Geosciences↗

USGS1993

A fortran program for Monte Carlo simulation of oil-field discovery sequences

We have developed a program for performing Monte Carlo simulation of oil-field discovery histories. A synthetic parent population of fields is generated as a finite sample from a distribution of specified form. The discovery sequence then is simulated by sampling without replacement from this parent population in accordance with a probabilistic discovery process model. The program computes a chi-squared deviation between synthetic and actual discovery sequences as a function of the parameters of the discovery process model, the number of fields in the parent population, and the distributional parameters of the parent population. The program employs the three-parameter log gamma model for the distribution of field sizes and employs a two-parameter discovery process model, allowing the simulation of a wide range of scenarios. ?? 1993.

Computers & Geosciences↗

USGS1992

Quantitative assessment of mineral resources with an application to petroleum geology

The probability of occurrence of natural resources, such as petroleum deposits, can be assessed by a combination of multivariate statistical and geostatistical techniques. The area of study is partitioned into regions that are as homogeneous as possible internally while simultaneously as distinct as possible. Fisher's discriminant criterion is used to select geological variables that best distinguish productive from nonproductive localities, based on a sample of previously drilled exploratory wells. On the basis of these geological variables, each wildcat well is assigned to the production class (dry or producer in the two-class case) for which the Mahalanobis' distance from the observation to the class centroid is a minimum. Universal kriging is used to interpolate values of the Mahalanobis' distances to all locations not yet drilled. The probability that an undrilled locality belongs to the productive class can be found, using the kriging estimation variances to assess the probability of misclassification. Finally, Bayes' relationship can be used to determine the probability that an undrilled location will be a discovery, regardless of the production class in which it is placed. The method is illustrated with a study of oil prospects in the Lansing/Kansas City interval of western Kansas, using geological variables derived from well logs. ?? 1992 Oxford University Press.

Nonrenewable Resources↗

USGS1990

Multivariate statistical analysis of stream-sediment geochemistry in the Grazer Paläozoikum, Austria

The Austrian reconnaissance study of stream-sediment composition — more than 30000 clay-fraction samples collected over an area of 40000 km 2 — is summarized in an atlas of regional maps that show the distributions of 35 elements. These maps, rich in information, reveal complicated patterns of element abundance that are difficult to compare on more than a small number of maps at one time. In such a study, multivariate procedures such as simultaneous R-Q mode components analysis may be helpful. They can compress a large number of variables into a much smaller number of independent linear combinations. These composite variables may be mapped and relationships sought between them and geological properties. As an example, R-Q mode components analysis is applied here to the Grazer Paläozoikum, a tectonic unit northeast of the city of Graz, which is composed of diverse lithologies and contains many mineral deposits.

Mineralium Deposita↗

USGS1990

Regionalization in geology by multivariate classification

The concept of multivariate classification of "geological objects" can be combined with the concept of regionalized variables to yield a procedure for typification of geological objects, such as rock units, well records, or samples. Numerical classification is followed by subdivision of the area of investigation, and culminates in a regionalization or mapping of the classification onto the plane. Regions are subdivisions of the map area which are spatially contiguous and relatively homogeneous in their geological properties. The probability of correct classification of each point within a region as being part of that region can be assessed in terms of Bayesian probability as a space-dependent function. The procedure is applied to subsurface data from western Kansas. The geologic properties used are quantitative variables, and relationships are expressed by Mahalanobis' distances. These functions could be replaced by other metrics if qualitative or binary data derived from geological descriptions or appraisals were included in the analysis. ?? 1990 International Association for Mathematical Geology.

Mathematical Geology↗