Quantitative Prediction and Evaluation of Mineral Resources Based on GIS: A Case Study in Sanjiang Region,Southwestern China

Quantitative prediction and evaluation of mineral resources are one of the important topics of mathematical geology. On the basis of GIS technologies and weights of evidence modeling, MapGIS is integrated with GIS and mineral-resource prediction and evaluation. The final product is a predictor map of posterior probabilities of occurrence of the discrete event within a small unit cell. Predictor layers were created on a digital database that includes 1:200,000 scale geological, and geochemical, and geophysical maps, and remote-sensing images in study area. According to metallogenetic factors extractiont and weights of evidence modeling, there are four main metal ore belts in the study area: (1) the Batang belt; (2) the Lei Wuqi belt; (3) the Basu-Chayu belt; and (4) the Ganzi-Litang belt. The predictor map of posterior probabilities show that 29% of study area as zones with potential for porphyry copper, and 81% known mineral occurrences success rate is circled in the metallogenetic posterior probabilities map. The results demonstrate plausibility of weights-of-evidence modeling of mineral potential in large areas with small number of mineral prospects.     

Evaluation of Weights of Evidence to Predict Epithermal-Gold Deposits in the Great Basin of the Western United States

The weights-of-evidence method provides a simple approach to the integration of diverse geologic information. The application addressed is to construct a model that predicts the locations of epithermal-gold mineral deposits in the Great Basin of the western United States. Weights of evidence is a data-driven method requiring known deposits and occurrences that are used as training sites in the evaluated area. Four hundred and fifteen known hot spring gold–silver, Comstock vein, hot spring mercury, epithermal manganese, and volcanogenic uranium deposits and occurrences in Nevada were used to define an area of 327.4 km² as training sites to develop the model. The model consists of nine weighted-map patterns that are combined to produce a favorability map predicting the distribution of epithermal-gold deposits. Using a measure of the association of training sites with predictor features (or patterns), the patterns can be ranked from best to worst predictors. Based on proximity analysis, the strongest predictor is the area within 8 km of volcanic rocks younger than 43 Ma. Being close to volcanic rocks is not highly weighted, but being far from volcanic rocks causes a strong negative weight. These weights suggest that proximity to volcanic rocks define where deposits do not occur. The second best pattern is the area within 1 km of hydrothermally altered areas. The next best pattern is the area within 1 km of known placer-gold sites. The proximity analysis for gold placers weights this pattern as useful when close to known placer sites, but unimportant where placers do not exist. The remaining patterns are significantly weaker predictors. In order of decreasing correlation, they are: proximity to volcanic vents, proximity to east-west to northwest faults, elevated airborne radiometric uranium, proximity to northwest to west and north-northwest linear features, elevated aeromagnetics, and anomalous geochemistry. This ordering of the patterns is a function of the quality, applicability, and use of the data. The nine-pattern favorability map can be evaluated by comparison with the USGS National Assessment for hot spring gold–silver deposits. The Spearman's ranked correlation coefficient between the favorability and the National Assessment permissive tracts is 0.5. Tabulations of the areas of agreement and disagreement between the two maps show 74% agreement for the Great Basin. The posterior probabilities for 51 significant deposits in the Great Basin, both used and not used in the model, show the following: 26 classified as favorable; 25 classified as permissive; and 1, Florida Canyon, classified as nonpermissive.The Florida Canyon deposit has a low favorability because there are no volcanic rocks near the deposit on the Nevada geologic map used. The largest areas of disagreement are caused by the USGS National Assessment team concluding that volcanic rocks older than 27 Ma in Nevada are not permissive, which was not assumed in this model. The weights-of-evidence model is evaluated as reasonable and delineates permissive areas for epithermal deposits comparable to expert's delineation. The weights-of-evidence model has the additional characteristics that it is well defined, reproducible, objective, and provides a quantitative measure of confidence.     

A Modified Weights-of-Evidence Method for Regional Mineral Resource Estimation

Weights-of-evidence (WofE) modeling and weighted logistic regression (WLR) are two methods of regional mineral resource estimation, which are closely related: For example, if all the map layers selected for further analysis are binary and conditionally independent of the mineral occurrences, expected WofE contrast parameters are equal to WLR coefficients except for the constant term that depends on unit area size. Although a good WofE strategy is supposed to achieve approximate conditional independence, a common problem is that the final estimated probabilities are biased. If there are N deposits in a study area and the sum of all estimated probabilities is written as S, then WofE generally results in S > N. The difference S − N can be tested for statistical significance. Although WLR yields S = N, WLR coefficients generally have relatively large variances. Recently, several methods have been developed to obtain WofE weights that either result in S = N, or become approximately unbiased. A method that has not been applied before consists of first performing WofE modeling and following this by WLR applied to the weights. This method results in modified weights with unbiased probabilities satisfying S = N. An additional advantage of this approach is that it automatically copes with missing data on some layers because weights of unit areas with missing data can be set equal to zero as is generally practiced in WofE applications. Some practical examples of application are provided.     

Geochemistry of waters and brines from the Salinas Grandes basin,Córdoba,Argentina. I. Geomorphology and hydrochemical characteristics

This paper reports on hydrochemical features of diluted waters in the source areas and the brine end-members dominant in the playa of the Salinas Grandes Basin, Córdoba, Argentina. Special emphasis was placed on the study of the relationship between geomorphology and the resulting hydrochemical fractionation. Inflow is from springs and mountain streams which disappear before reaching the saline complex. The playa and intermittent saline lakes are mainly fed by groundwater flow and a few moderately saline and perennial springs. Conversely, ephemeral lakes are fed by atmospheric precipitation and groundwater, whereas small ponds are only fed by atmospheric precipitation. The absence of a clear linkage between geomorphological units and water types was evident in the source areas. Up to four types of water were recognized in a given geomorphological unit. From the sandflat downward towards the playa, the correspondence between geomorphological units and water types was clear, coinciding with a decreased hydrochemical heterogeneity. In this subenvironment, I have recognized two dominant types of water (SO ₄ ^2– –Cl^––HCO ₃ ^– –Na⁺ and Cl^––SO ₄ ^2– –HCO ₃ ^– –Na⁺), which can be considered the original members of the neutral brine in the playa (Cl^––SO ₄ ^2– –Na⁺ where Ca²⁺+Mg²⁺ do not surpass 5 meq per cent), and an intermediate type in the mudflat (Cl^––SO ₄ ^2– –Na⁺ where the contribution of Ca²⁺+Mg²⁺ reaches up to 15 meq per cent). It seems evident that in the zone between the source areas and the sandflat, hydrochemistry is governed by chemical weathering. In stream floodplains and in the distal alluvial plain, close to the saline complex, the increase in water types as well as the increase of HCO ₃ ^– with respect to SO ₄ ^2– , were explained by the mixture of aquifers controlled by the Salinas Grandes-Salinas de Ambargasta fracture. In the saline complex, the more concentrated end-members are the result of evaporation of the two more frequent water-types in the sandflat subenvironment, and salt dissolution of ancient evaporite deposit.     

Growth, ionic and osmotic relations of an Allenrolfea occidentalis population in an inland salt playa of the Great Basin Desert

Variation in growth, physiology and ionic relations patterns of Allenrolfea occidentalis, a perennial halophyte of dry habitats, was studied under field conditions from May 1996 to November 1997. An A. occidentalis community has a characteristic soil pH of 7·3–8·3. During the two years, the population was exposed to great variations in soil salinity, from 29 to 146 dS m⁻¹, and soil moisture, ranging from drought (9·2%) to wet (19%). The salt concentrations were significantly higher in the surface soil layers than in the subsurface layers. Seasonal changes in dry weight are directly related to soil salinity stress. Allenrolfea occidentalis had greater growth and biomass production under saline conditions. Na⁺and Cl⁻ions were accumulated in plant tissues in much greater amounts than K⁺, Ca²⁺, and Mg²⁺. Soil salinities were significantly reduced at the end of the growing season. Water potentials of the shoots decreased significantly with increasing salinity. The plant (F_v/F_mratio) was more affected by salinity and irradiation levels during the summer period.     

Extended weights of evidence modeling for the pseudo-estimation of metal grades

The weights of evidence modeling (WEM) for binary patterns is extended to take account of general categorical variables. The extension makes it possible to use the weights of evidence model in estimating the conditional probability distributions of metal grades. First, the target feature is converted into a set of binary target indicators. Second, the posterior probabilities are estimated for each of the target categories. Third, the estimates are combined to yield the posterior probability distribution of the target feature. Finally, the pseudometal estimates are derived from the probability distribution. The metal grade estimates are prefixed with pseudo, because the estimates are created from indirect evidence (explanatory variables). The pseudo-estimates provide a unique quantitative means to the delineation of exploration targets. This advantage reduces the ambiguities of target selection based solely on probability estimates. In order to use the generalized WEM, continuous geoscience attributes must be converted into categorical variables by means of optimal segmentation based on the target attribute of interest. The segmentation may be viewed as a process of defining evidence of the target feature. The extended weights of evidence model is demonstrated on a case study to select gold targets of Carlin type. The dataset used in the modeling includes apparent resistivity fields, soil geochemical samples, lithological and alteration information, and structural data.     

Diatom-based models for reconstructing past water quality and productivity in New Zealand lakes

The trophic status of lakes in New Zealand is, on average, low compared to more densely populated areas of the globe. Despite this, trends of eutrophication are currently widespread due to recent intensification in agriculture. In order to better identify baseline productivity and establish long-term trends in lake trophic status, diatom-based transfer functions for productivity-related parameters were developed. Water quality data and surface sediment diatom assemblages from 53 lakes across the North and South Islands of New Zealand were analysed to determine species responses to the principal environmental gradients in the data set. Repeat sampling of water chemistry over a 12-month period enabled examination of species responses to annual means as well as means calculated for stratified and mixed periods. Variables found to be most strongly correlated with diatom species distributions were chlorophyll a (Chl a), total phosphorus (TP), dissolved reactive phosphorus (DRP), ionic concentration (measured as electrical conductivity (EC)) and pH. These variables were used to develop diatom-based transfer functions using weighted averaging regression and calibration (simple, tolerance down-weighted and with partial least squares algorithm applied). Overall, models derived for stratified means were weaker than those using annual or isothermal means. For specific variables, the models derived for the isothermal mean of EC (WA-tol r²_jack = 0.79; RMSEP = 0.15 log₁₀ S cm^–1),the annual mean of pH (WA r²_jack = 0.72; RMSEP = 0.25 pH units) and the isothermal mean of Chl a (WA r²_jack = 0.71; RMSEP = 0.18 log₁₀ mg m^–3 Chl a) performed best. The models derived for TP were weak in comparison (for the annual mean of TP: WA r²_jack = 0.50; RMSEP = 0.24 log₁₀ mg m^–3 TP) and residuals on estimates for this model were correlated with several other water quality variables, suggesting confounding of species responses to TP concentrations. The model derived for the isothermal mean of DRP was relatively strong (WA-tol r²_jack = 0.78; RMSEP = 0.17 log₁₀ mg m^–3 DRP); however, residual values for this model were also found to be strongly correlated with several other water quality variables. It is concluded that the poor performance of the TP and DRP transfer functions relative to that of the Chl a model reflects the coexistence of nitrogen and phosphorus limitation within the lakes in the data set. In spite of this, the suite of transfer functions developed from the training set is regarded as a valuable addition to palaeolimnological studies in NewZealand.     

Fractal pattern integration for mineral potential estimation

Concepts of fractal/multifractal dimensions and fractal measure were used to derive the prior and posterior probabilities that a small unit cell on a geological map contains one or more mineral deposits. This has led to a new version of the weights of evidence technique which is proposed for integrating spatial datasets that exhibit nonfractal and fractal patterns to predict mineral potential. The method is demonstrated with a case study of gold mineral potential estimation in the Iskut River area, northwestern British Columbia. Several geological, geophysical, and geochemical patterns (Paleozoic-Mesozoic sedimentary and volcanic clastic rocks; buffer zones around the contacts between sedimentary rocks and Mesozoic intrusive rocks; a linear magnetic anomaly; and geochemical anomalies for Au and associated elements in stream sediments) were integrated with the gold mineral occurrences which have fractal and multifractal properties with a box-counting dimension of 1.335±0.077 and cluster dimension of 1.219±0.037.     

Fuzzy Weights of Evidence Method and Its Application in Mineral Potential Mapping

This paper proposes a new approach of weights of evidence method based on fuzzy sets and fuzzy probabilities for mineral potential mapping. It can be considered as a generalization of the ordinary weights of evidence method, which is based on binary or ternary patterns of evidence and has been used in conjunction with geographic information systems for mineral potential mapping during the past few years. In the newly proposed method, instead of separating evidence into binary or ternary form, fuzzy sets containing more subjective genetic elements are created; fuzzy probabilities are defined to construct a model for calculating the posterior probability of a unit area containing mineral deposits on the basis of the fuzzy evidence for the unit area. The method can be treated as a hybrid method, which allows objective or subjective definition of a fuzzy membership function of evidence augmented by objective definition of fuzzy or conditional probabilities. Posterior probabilities calculated by this method would depend on existing data in a totally data-driven approach method, but depend partly on expert's knowledge when the hybrid method is used. A case study for demonstration purposes consists of application of the method to gold deposits in Meguma Terrane, Nova Scotia, Canada.     

Development of a diatom-based model for inferring total phosphorus in southeastern Australian water storages

Diatom-based transfer functions for inferring epilimnetic total phosphorus (TP) have been developed from a data set of 33 southeastern Australian water storages. Regular institutional monitoring of these sites has allowed comparison of models developed from TP data covering different time periods. A model based on mean annual TP performs better than models derived from winter maximum TP, spring minimum TP or TP nearest the time of diatom sampling. A mean annual TP model (WA-PLS 2 component) has a jack-knifed diatom-inferred versus measured TP correlation coefficient (r ² _jack) of 0.69 and a root-mean-square-error of prediction (RMSEP) of 0.246 log₁₀g TP l^–1, while alternative models have RMSEP > 0.27. Deletion of two samples with uncharacteristic species composition and environmental conditions improved performance of the mean annual TP model (r ² _jack= 0.74; RMSEP = 0.233 log₁₀g TP l^–1). Comparison with other published diatom-TP calibration models indicates that this model performs relatively well, with possible contributing factors including the extensive characterisation of TP (with an average 15 determinations making up the annual mean) and the dominance of planktonic diatoms in most sites. Downcore application of the model will allow the reconstruction of reservoir nutrient histories since commissioning, and thus provide a basis for understanding and management of reservoirs.     

A Hybrid Fuzzy Weights-of-Evidence Model for Mineral Potential Mapping

This paper describes a hybrid fuzzy weights-of-evidence (WofE) model for mineral potential mapping that generates fuzzy predictor patterns based on (a) knowledge-based fuzzy membership values and (b) data-based conditional probabilities. The fuzzy membership values are calculated using a knowledge-driven logistic membership function, which provides a framework for treating systemic uncertainty and also facilitates the use of multiclass predictor maps in the modeling procedure. The fuzzy predictor patterns are combined using Bayes’ rule in a log-linear form (under an assumption of conditional independence) to update the prior probability of target deposit-type occurrence in every unique combination of predictor patterns. The hybrid fuzzy WofE model is applied to a regional-scale mapping of base-metal deposit potential in the south-central part of the Aravalli metallogenic province (western India). The output map of fuzzy posterior probabilities of base-metal deposit occurrence is classified subsequently to delineate zones with high-favorability, moderate favorability, and low-favorability for occurrence of base-metal deposits. An analysis of the favorability map indicates (a) significant improvement of probability of base-metal deposit occurrence in the high-favorability and moderate-favorability zones and (b) significant deterioration of probability of base-metal deposit occurrence in the low-favorability zones. The results demonstrate usefulness of the hybrid fuzzy WofE model in representation and in integration of evidential features to map relative potential for mineral deposit occurrence.     

Identification of steps and pools from stream longitudinal profile data

Field research on step–pool channels has largely focused on the dimensions and sequence of steps and pools and how these features vary with slope, grain sizes and other governing variables. Measurements by different investigators are frequently compared, yet no means to identify steps and pools objectively have been used. Automated surveying instruments record the morphology of streams in unprecedented detail making it possible to objectively identify steps and pools, provided an appropriate classification procedure can be developed.To achieve objective identification of steps and pools from long profile survey data, we applied a series of scale-free geometric rules that include minimum step length (2.25% of bankfull width (W_b)), minimum pool length (10% of W_b), minimum residual depth (0.23% of W_b), minimum drop height (3.3% of W_b), and minimum step slope (10° greater than the mean slope). The rules perform as well as the mean response of 11 step–pool researchers who were asked to classify four long profiles, and the results correspond well with the channel morphologies identified during the field surveys from which the long profiles were generated. The method outperforms four other techniques that have been proposed. Sensitivity analysis shows that the method is most sensitive to the choice of minimum pool length and minimum drop height.Detailed bed scans of a step–pool channel created in a flume reveal that a single long profile with a fixed sampling interval poorly characterizes the steps and pools; five or more long profiles spread across the channel are required if a fixed sampling interval is used and the data suggest that survey points should be located more frequently than the diameter of the step-forming material. A single long profile collected by a surveyor who chooses breaks in slope and representative survey points was found to adequately characterize the mean bed profile.     

Salinity tolerance in Schinopsis quebracho colorado: Seed germination, growth, ion relations and metabolic responses

Soil salinity is a major abiotic stress influencing plant productivity worldwide. Schinopsis quebracho colorado is one of the most important woody species in the Gran Chaco, an arid and salt-prone subtropical biome of South America. To gain a better understanding of the physiological mechanisms that allow plant establishment under salt conditions, germination and seedling growth of S. quebracho colorado were examined under treatment with a range of NaCl solutions (germination: 0–300 mmol l⁻¹ NaCl; seedling growth: 0–200 mmol l⁻¹ NaCl). The aim was to test the hypothesis that S. quebracho colorado is a glycophite that shows different salt tolerance responses with development stage. Proline content, total soluble carbohydrates and Na⁺, K⁺ and Cl⁻ concentrations in leaves and roots of seedlings, and the chlorophyll concentration and relative water content of leaves were measured. Germination was not affected by 100 mmol l⁻¹ NaCl, but decreased at a concentration of 200 mmol l⁻¹. At 300 mmol l⁻¹ NaCl, germination did not occur. Seedling growth decreased drastically with increasing salinity. An increase in NaCl from 0 to 100 mmol l⁻¹ also significantly reduced the leaf relative water content by 22% and increased the proline concentration by 60% in roots. In contrast, total soluble carbohydrates were not significantly affected by salinity. Seedlings showed a sodium exclusion capacity, and there was an inverse correlation between Cl⁻ concentration and the total chlorophyll concentration. S. quebracho colorado was more tolerant to salinity during germination than in the seedling phase. The results suggest that this increased tolerance during germination might, in part, be the result of lower sensitivity to high tissue Na⁺ concentrations. The significant increment of proline in the roots suggests the positive role of this amino acid as a compatible solute in balancing the accumulation of Na⁺ and Cl⁻ as a result of salinity. These results help clarify the physiological mechanisms that allow establishment of S. quebracho colorado on salt-affected soils in arid and semi-arid Gran Chaco.     

Influence of shrubs on soil characteristics and their function in Sahelian agro-ecosystems in semi-arid Niger

Soil beneath shrubs form ‘fertile islands’ in fallow sites and millet fields in semi-arid Niger. To gain more information about this phenomena different shrub species in fallow sites following a gradient from 350–650 mm precipitation were examined. For each shrub two different areas were distinguished: an area under the canopy of the shrubs and an area in the nearby open land. Soil samples were taken from a depth of 0–10 cm and analysed for C_org, N_total, P_Bray, pH(H₂O), exchangeable cations, effective cation exchange capacity (ECEC) and soil texture. Significantly higher concentrations between 38–51% for C, N, P and 22% on ECEC for K⁺were found in the soil under the shrubs. The pH showed only slight but significant differences, whereas Al³⁺and H⁺rates on ECEC under the shrubs were increased by 44–55%. For Guiera senegalensis, the most common shrub of the studied area, enrichment ratios of most soil properties increased relatively more with increasing aridity. In general, enrichment ratios decreased with the age of the fallows, whereas concentrations showed no clear evolution. The chemical composition of the shrub litter seems to influence the degree of soil enrichment. The main step of fertile island formation takes place during the cultivation period by trapping wind-blown sediment. This work shows that shrubs are of vital importance for the accumulation of nutrients and maintenance of soil fertility within agro-ecosystems of Niger.     

The spatial characterization of turbulence around large roughness elements in a gravel-bed river

This study characterizes the flow field above and around multiple instream submerged cobbles, boulders, and pebble clusters in order to obtain a better understanding of the hydrodynamics associated with large roughness elements (LREs) in gravel-bed rivers. Spatially distributed high frequency, three-dimensional velocity measurements were recorded in situ using acoustic Doppler velocimeters at different flow stages. The spatial distributions of turbulent kinetic energy, k_e, longitudinal component integral timescales, ITS_u, and Reynolds shear stresses were characterized and are presented for selected sites. The longitudinal–vertical Reynolds shear stress increased with flow stage more strongly than the longitudinal–lateral or lateral–vertical Reynolds shear stresses and dominate at the highest measured flows. Canonical redundancy analysis was used to relate LRE morphometrics and mean flow conditions to the turbulence parameters estimated in the LRE wakes (i.e., k_e, ITS_u, and Reynolds shear stresses). LRE size and mean unobstructed velocity explained the highest proportion of the variance in the turbulent wake statistics. Multivariate regression models based on LRE width, mean unobstructed longitudinal velocity and flow depth are presented offering a tool to predict LRE wake turbulence.     

Downstream hydraulic geometry of rivers in Victoria, Australia

Downstream variation of hydraulic geometry in rivers, characterized by fine textured banks and low width to depth ratios (7–25), is investigated in Victoria, Australia, with the aim of developing predictive models of channel geometry for large-scale spatial modeling applications. A one-dimensional hydraulic model is used to determine the mean bank-full geometry and discharge (Q_bf) for 93 sites which are investigated in relation to discharge of fixed average recurrence interval (ARI). The median bank-full ARI is estimated at 0.8 years with 75% of sites between 0.5 and 2.5 years. Exponents in the downstream hydraulic geometry relations for width, depth and velocity are respectively 0.43, 0.40 and 0.18 (Q = Q_bf) and 0.44, 0.38 and 0.03 (Q = Q₂, i.e., 2-year ARI), falling near the mode of global values. Q₂ and slope explain 66% of variance in Q_bf, while Q₂ explains 73% and 69% of the variance in width and depth relations, respectively: Q₂ provides a reliable substitute for Q_bf in spatial modeling applications. Spatial variation in hydraulic geometry relations within and between river basins remains largely unexplained. The W/D ratio characteristically decreases with increasing distance along the lower reaches of most rivers and this has contributed to the lower than expected value for the width exponent.     

Estimating Petroleum Resources Using Geo-Anchored Method--A Sensitivity Study

The Geo-anchored method, based on a moment-type estimator, has been developed for estimating parent population properties from a successive sample of discoveries. By substituting the expectation of the waiting time z _(n+1) of the (n + 1) ^th discovery to occurrence for an unknown parameter in the anchored method, the Geo-anchored method allows estimation of inclusion probabilities directly from observed data, thus eliminating the need for a priori selection of a value of N, R, or some other feature of the parent population. Because direct estimation of N and R requires an ordered sample, the Geo-anchored method is more sensitive to the data-generating process than the anchored method. This paper presents a sensitivity study on the Geo-anchored method. The test is based on simulated discovery sequences with different assumptions regarding discovery efficiency, exploration maturity, and the shape of the parent field size distribution. As a reference for comparison, estimates from the Horvitz–Thompson estimator also are presented.     

Prescribed fire in a Great Basin sagebrush ecosystem: Dynamics of soil extractable nitrogen and phosphorus

Pinyon and juniper have been expanding into sagebrush (Artemisia tridentata) ecosystems since settlement of the Great Basin around 1860. Herbaceous understory vegetation is eliminated as stand densities increase and the potential for catastrophic fires increases. Prescribed fire is increasingly used to remove trees and promote recovery of sagebrush ecosystems. We quantified the effects of prescribed fire, vegetation type, and time following fire on soil KCl extractable nitrogen and NaHCO₃ extractable phosphorus in a pinyon–juniper woodland and its associated sagebrush ecosystem immediately before and for 4 years after a spring prescribed burn. Potassium chloride extractable NH₄⁺ and total inorganic-N increased immediately following prescribed fire, and extractable NO₃⁻ decreased immediately after the burn. In the surface layer (top 8 cm), extractable NH₄⁺ remained elevated compared to the control through year 2 after the burn. By the first fall post-burn extractable NO₃⁻ and total extractable inorganic-N increased and remained elevated over the control through year 3 after the burn in the surface layer. For the entire soil profile (52 cm), the burn had no effect on NH₄⁺, and the effects on total extractable inorganic-N were no longer significant after year 1. However, NO₃⁻ remained elevated over the control through year 2 post-fire for the soil profile. Near surface NaHCO₃ extractable ortho-P increased immediately following fire, and remained elevated through year 2 post-fire. No fire effects were observed for extractable ortho-P in deeper horizons. Our data show that plant available nitrogen can remain elevated for extended periods following prescribed fire. This can influence regrowth and seedling establishment of native plant species, invasion of exotic plant species and, ultimately, site recovery potential.     

Plant water relations of thornscrub shrub species, north-eastern Mexico

As an approach to understand how diurnal and seasonal plant water potentials (Ψ) are related to soil water-content and evaporative demand components, the responses of six thornscrub shrubs species (Havardia pallens, Acacia rigidula, Eysenhardtia texana, Diospyros texana, Randia rhagocarpa, and Bernardia myricaefolia) of the north-eastern region of Mexico to drought stress were investigated during the course of 1 year. All study species showed the typical diurnal pattern of variation in Ψ. That is, Ψ decreased gradually from predawn (Ψ_pd) maximal values to reach minima at midday (Ψ_md) and began to recover in the late afternoon. On a diurnal basis and with adequate soil water-content (>0.20 kg kg⁻¹), diurnal Ψ values differed among shrub species and were negatively and significantly (p<0.001) correlated with air temperature (r=−0.741 to −0.883) and vapor pressure deficit (r=−0.750 to −0.817); in contrast, a positive and significant (p<0.001) relationship was found to exist with relative humidity (r=0.758–0.842). On a seasonal basis, during the wettest period (soil water-content>0.20 kg kg⁻¹), higher Ψ_pd (−0.10 MPa) and Ψ_md (−1.16 MPa) values were observed in R. rhagocarpa, whereas lower figures (−0.26 and −2.73 MPa, respectively) were detected in A. rigidula. On the other hand, during the driest period (soil water-content<0.15 kg kg⁻¹), Ψ_pd and Ψ_md values were below −7.3 MPa; i.e. when shrubs species faced severe water deficit. Soil water-content at different soil layers, monthly mean relative humidity and monthly precipitation were significantly correlated with both Ψ_pd (r=0.538–0.953; p<0.01) and Ψ_md (r=0.431–0.906; p<0.05). Average soil water-content in the 0–50 cm soil depth profile explained between 70% and 87% of the variation in Ψ_pd. Results have shown that when gravimetric soil water-content values were above 0.15 kg kg⁻¹, Ψ_pd values were high and constant; below this threshold value, Ψ declined gradually. Among all shrub species, A. rigidula appeared to be the most drought tolerant of the six species since during dry periods it tends to sustain significantly higher Ψ_pd in relation to B. myricaefolia. The remaining species showed an intermediate pattern. It is concluded that the ability of shrub species to cope with drought stress depends on the pattern of water uptake and the extent to control water loss through the transpirational flux.     

Conditional estimates of the number of podiform chromite deposits

A desirable guide for estimating the number of undiscovered mineral deposits is the number of known deposits per unit area from another well-explored permissive terrain. An analysis of the distribution of 805 podiform chromite deposits among ultramafic rocks in 12 subareas of Oregon and 27 counties of California is used to examine and extend this guide. The average number of deposits in this sample of 39 areas is 0.225 deposits per km² of ultramafic rock; the frequency distribution is significantly skewed to the right. Probabilistic estimates can be made by using the observation that the lognormal distribution fits the distribution of deposits per unit area. A further improvement in the estimates is available by using the relationship between the area of ultramafic rock and the number of deposits.The number (N) of exposed podiform chromite deposits can be estimated by the following relationship: log₁₀(N)=–0.194+0.577 log₁₀(area of ultramafic rock). The slope is significantly different from both 0.0 and 1.0. Because the slope is less than 1.0, the ratio of deposits to area of permissive rock is a biased estimator when the area of ultramafic rock is different from the median 93 km². Unbiased estimates of the number of podiform chromite deposits can be made with the regression equation and 80 percent confidence limits presented herein.