USGS Science⌕ Search

SEARCH · USGS Science

Results for “Scientific Data - Nature”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,153 records · Page 64Linked to original sources

Simulation of zones of contribution to wells at site GM–38, Naval Weapons Industrial Reserve Plant, Bethpage, New York

A three-dimensional groundwater-flow model is coupled with the particle-tracking program MODPATH to delineate zones of contribution to wells pumping from the Magothy aquifer and supplying water to a chlorinated volatile organic compound removal plant at site GM–38, Naval Weapons Industrial Reserve Plant, Bethpage, New York. By use of driller’s logs, a transitional probability approach generated three alternative realizations of heterogeneity within the Magothy aquifer to assess uncertainty in model representation. Finer-grained sediments with low hydraulic conductivity were realized as laterally discontinuous, thickening towards the south, and comprising about 17 percent of the total aquifer volume. Particle-tracking evaluations of a steady state present conditions model with alternative heterogeneity realizations were used to develop zones of contribution of remedial pumping wells. Because of heterogeneity and high rates of advection within the coarse-grained sediments, transport by dispersion and (or) diffusion was assumed to be negligible. Resulting zones of contribution of existing remedial wells are complex shapes, influenced by heterogeneity of each realization and other nearby hydrologic stresses. The use of two particle tracking techniques helped identify zones of contribution to wells. Backtracking techniques and observations of points of intersection of backward-tracked particles at shells of the GM–38 Hot Spot, as defined by surfaces of equal total volatile organic compound concentration, identified the source of water within the GM–38 Hot Spot to simulated wells. Forward-tracking techniques identified the fate of water within the GM–38 Hot Spot, including well capture and discharge to model constant head and drain boundaries. The percentage of backward-tracked particles, started at GM–38 wells that were sourced from within the Hot Spot, varied from 72.0 to 98.2, depending on the Hot Spot delineation used (present steady state model and Magothy aquifer heterogeneity realization A). The percentage of forward-tracked particles that were captured by GM–38 wells varied from 81.1 to 94.6, depending on the Hot Spot delineation used, with the remainder primarily captured by Bethpage Water District Plant 4 production wells (present steady state model and Magothy aquifer heterogeneity realization A). Less than 1 percent of forward-tracked particles ultimately discharge at model constant head and drain boundaries. The differences between forward- and backward-tracked particle percentage ranges are due to some forward-tracked particles not being captured by GM–38 wells, and some backward-tracked particles not intersecting specific regions of the Hot Spot. During 2013, an aquifer test generated detailed time series of well pumping rates and corresponding water-level responses were recorded at numerous locations. These data were used to verify the present conditions steady state model and demonstrate the sensitivity of model results to transient-state changes.

New York↗

Hydroclimate of the Spring Mountains and Sheep Range, Clark County, Nevada

Precipitation, potential evapotranspiration, and actual evapotranspiration often are used to characterize the hydroclimate of a region. Quantification of these parameters in mountainous terrains is difficult because limited access often hampers the collection of representative ground data. To fulfill a need to characterize ecological zones in the Spring Mountains and Sheep Range of southern Nevada, spatially and temporally explicit estimates of these hydroclimatic parameters are determined from remote-sensing and model-based methodologies. Parameter-elevation Regressions on Independent Slopes Model (PRISM) precipitation estimates for this area ranges from about 100 millimeters (mm) in the low elevations of the study area (700 meters [m]) to more than 700 mm in the high elevations of the Spring Mountains (> 2,800 m). The PRISM model underestimates precipitation by 7–15 percent based on a comparison with four high‑elevation precipitation gages having more than 20 years of record. Precipitation at 3,000-m elevation is 50 percent greater in the Spring Mountains than in the Sheep Range. The lesser amount of precipitation in the Sheep Range is attributed to partial moisture depletion by the Spring Mountains of eastward-moving, cool-season (October–April) storms. Cool-season storms account for 66–76 percent of annual precipitation. Potential evapotranspiration estimates by the Basin Characterization Model range from about 700 mm in the high elevations of the Spring Mountains to 1,600 mm in the low elevations of the study area. The model realistically simulates lower potential evapotranspiration on northeast-to-northwest facing slopes compared to adjacent southeast-to-southwest facing slopes. Actual evapotranspiration, estimated using a Moderate Resolution Imaging Spectroradiometer based water-balance model, ranges from about 100 to 600 mm. The magnitude and spatial variation of simulated, actual evapotranspiration was validated by comparison to PRISM precipitation. Estimated groundwater recharge, computed as the residual of precipitation depleted by actual evapotranspiration, is within the range of previous estimates. A climatic water deficit dataset and aridity-index-based climate zones are derived from precipitation and evapotranspiration datasets. Climate zones range from arid in the lower elevations of the study area to humid in small pockets on north- to northeast-facing slopes in the high elevations of the Spring Mountains. Correlative analyses between hydroclimatic variables and mean ecosystem elevations indicate that the climatic water deficit is the best predictor of ecosystem distribution (R 2 = 0.92). Computed water balances indicate that substantially more recharge is generated in the Spring Mountains than in the Sheep Range. A geospatial database containing compiled and developed hydroclimatic data and other pertinent information accompanies this report.

Nevada↗

Hydrologic characteristics of low-impact stormwater control measures at two sites in northeastern Ohio, 2008-13

This report updates and examines hydrologic data gathered to characterize the performance of two stormwater-control measure (SCM) sites in the Chagrin River watershed, Ohio. At the Sterncrest Drive site, roadside bioswales and rain gardens were used to alleviate drainage problems in this residential neighborhood area. At the Washington Street site, a treatment train (including a pervious-paver system, rain garden, and bioswales) was used to reduce and delay stormwater runoff at a small business development. Selected metrics were used to demonstrate SCM system performance with regard to stormwater-management objectives at each site. Rain-garden overflow-frequency data collected at the Sterncrest Drive site during 2008–13 were used to characterize system sensitivity to rainfall characteristics. Approximately 70 percent of storms exceeding 0.75 inches during 3 hours or more resulted in overflows. Drainage-design features that may restrict flow through the system were identified. Overall, the data and local observations confirmed the continued success of the SCM at the Sterncrest Drive site in preventing roadway closure due to flooding. The additional years of data collected at the Washington Street site indicated that a previous analysis of increased runoff removal, based on only the first 2 years (2009–10) of data, provided premature conclusions. With 5 years of data (2009–13) and adjusting for changes in rainfall characteristics, it appears that the percentage of runoff removed by the system is decreasing; however, the lag time (time from onset of rainfall to runoff) has remained nearly constant. The annual mean percent removal for 2010–13 ranged from 55 to 37 percent with an overall mean of 45 percent, and this does meet the project objective of reducing runoff from the business complex. One possible explanation for the combination of increased volume of runoff and no change in the timing of runoff is the preferential flow paths developed in the SCM, increasing the capacity for internal drainage. Data indicated that the SCM system at the Washington Street site had reduced functionality over time.

Ohio↗

Hydrology of the Claiborne aquifer and interconnection with the Upper Floridan aquifer in southwest Georgia

The U.S. Geological Survey conducted a study, in cooperation with the Georgia Environmental Protection Division, to define the hydrologic properties of the Claiborne aquifer and evaluate its connection with the Upper Floridan aquifer in southwest Georgia. The effort involved collecting and compiling hydrologic data from the aquifer in subarea 4 of southwestern Georgia. Data collected for this study include borehole geophysical logs in 7 wells, and two 72-hour aquifer tests to determine aquifer properties. The top of the Claiborne aquifer extends from an altitude of about 200 feet above the North American Vertical Datum of 1988 (NAVD 88) in Terrell County to 402 feet below NAVD 88 in Decatur County, Georgia. The base of the aquifer extends from an altitude of about 60 feet above NAVD 88 in eastern Sumter County to about 750 feet below NAVD 88 in Decatur County. Aquifer thickness ranges from about 70 feet in eastern Early County to 400 feet in Decatur County. The transmissivity of the Claiborne aquifer, determined from two 72-hour aquifer tests, was estimated to be 1,500 and 700 feet squared per day in Mitchell and Early Counties, respectively. The storage coefficient was estimated to be 0.0006 and 0.0004 for the same sites, respectively. Aquifer test data from Mitchell County indicate a small amount of leakage occurred during the test. Groundwater-flow models suggest that the source of the leakage was the underlying Clayton aquifer, which produced about 2.5 feet of drawdown in response to pumping in the Claiborne aquifer. The vertical hydraulic conductivity of the confining unit between the Claiborne and Clayton aquifers was simulated to be about 0.02 foot per day. Results from the 72-hour aquifer tests run for this study indicated no interconnection between the Claiborne and overlying Upper Floridan aquifers at the two test sites. Additional data are needed to monitor the effects that increased withdrawals from the Claiborne aquifer may have on future water resources.

Georgia↗

Simulation of potential groundwater recharge for the glacial aquifer system east of the Rocky Mountains, 1980–2011, using the Soil-Water-Balance Model

An understanding of the spatial and temporal extent of groundwater recharge is critical for many types of hydrologic assessments involving water quality, contaminant transport, ecosystem health, and sustainable use of groundwater. Annual potential groundwater recharge was simulated at a 1-kilometer resolution with the Soil-Water-Balance (SWB) model for the glacial aquifer system east of the Rocky Mountains, from central Montana east to Maine, for calendar years 1980–2011. The SWB model used high resolution meteorological, land cover, and soil hydrology datasets that are nationally consistent and publicly available. The SWB model computed daily potential groundwater recharge as precipitation in excess of interception, runoff, evapotranspiration, and soil-water storage capacity. Daily potential recharge values within each year of the simulation were summed to produce annual potential recharge rates. Potential recharge as described in this report is water that infiltrates vertically below the plant rooting zone and is assumed to reach the water table. The calibrated SWB model in this report is called the glacial SWB model. Model calibration assumed that the area contributing to groundwater discharge equaled the surface watershed. The model was calibrated to stream base flows from 39 watersheds throughout the model domain that had hydrologic conditions appropriate for hydrograph separation. Base flows were calculated from daily streamflow records with the HYSEP local minimum hydrograph separation method The glacial SWB model reproduced the mean annual base-flow calibration targets well; the Nash-Sutcliffe efficiency coefficient was 0.94, and the root mean squared error was 1.28 inches per year. The glacial SWB model provides insight into the spatial and temporal variability in potential annual recharge across the glacial aquifer system. About 20 percent of the active model area had an average potential recharge rate of less than 1 inch per year. Total precipitation, total recharge, and recharge as a percentage of precipitation increased from west to east. A substantial amount of the recharge water (39 percent) entering the glacial aquifer system travels through developed (urbanized) and agricultural landscapes, which are known to cause water-quality impairments. Regional climatic events, such as the 1988 to 1989 drought, are apparent in the potential recharge time series. Potential recharge generally increased across the glacial aquifer system between 2001 and 2011. A comparison of the potential recharge from the glacial SWB model to previous broad-scale recharge estimates reveals several important considerations for future SWB modeling applications. Shifts in the overall distribution of potential recharge between separate models can be explained by methods used to generate base-flow calibration target datasets. Spatial patterns in potential recharge simulated by SWB models are strongly dependent on the data and assumptions used to assign model cells to hydrologic soil groups. A review of several SWB models used to estimate groundwater recharge (and not surface runoff) revealed that model results are most sensitive to input climatic data, followed by surface runoff (curve number) and root-zone depth parameters.

Scientific Investigations Report↗

Status of water-level altitudes and long-term water-level changes in the Chicot and Evangeline (undifferentiated) and Jasper aquifers, greater Houston area, Texas, 2021

Since the early 1900s, groundwater withdrawn from the primary aquifers that compose the Gulf Coast aquifer system—the Chicot and Evangeline (undifferentiated) and Jasper aquifers—has been the primary source of water in the greater Houston area, Texas. This report, prepared by the U.S. Geological Survey in cooperation with the Harris-Galveston Subsidence District, City of Houston, Fort Bend Subsidence District, Lone Star Groundwater Conservation District, and Brazoria County Groundwater Conservation District, is one in an annual series of reports depicting the status of water-level altitudes and water-level changes in aquifers in the greater Houston area. In contrast to previous reports, the Chicot and Evangeline aquifers are treated as a single hydrogeologic unit in this report. In 2021, shaded depictions of water-level altitudes for the Chicot and Evangeline aquifers (undifferentiated) ranged from 300 feet (ft) below the North American Vertical Datum of 1988 (NAVD 88) to 300 ft above NAVD 88. The largest decline in water-level altitudes indicated by the 1977–2021 long-term water-level-change map for the Chicot and Evangeline aquifers (undifferentiated) was in the north-central part of The Woodlands, Tex., whereas the 1990–2021 long-term water-level-change map for the Chicot and Evangeline aquifers (undifferentiated) depicts a large area of decline in water-level altitudes in northwestern Harris County, northwest of Jersey Village, Tex. The largest rise in water-level altitudes in the Chicot and Evangeline aquifers (undifferentiated) was observed in a relatively large area in southeastern Harris County for 1977–2021, whereas the largest rise in water-level altitudes for 1990–2021 was in a relatively large area in central Harris County. In 2021, shaded depictions of water-level altitudes for the Jasper aquifer ranged from 250 ft below NAVD 88 to 300 ft above NAVD 88. The 2000–21 long-term water-level-change map for the Jasper aquifer depicts water-level declines throughout most of the study area where water-level-altitude data from the Jasper aquifer were collected, with the largest decline in northern Harris County southwest of The Woodlands.

Texas↗

Flood-frequency estimates for Kentucky streamgages based on data through water year 2021 and results of updating the fundamental layers in Kentucky StreamStats

The U.S. Geological Survey, in cooperation with the Kentucky Transportation Cabinet, analyzed flood-frequency statistics for streamgages in Kentucky. Using annual peak-flow data through water year 2021, flood-frequency estimates were computed for 261 streamgages, including unregulated and regulated sites as well as sites with mixed regulation records. Methods followed those outlined in “Guidelines for Determining Flood Flow Frequency—Bulletin 17C” (U.S. Geological Survey Techniques and Methods 4–B5). These estimates included flows corresponding to annual exceedance probabilities of 50, 20, 10, 4, 2, 1, and 0.2 percent. Temporal trend analyses using the Mann-Kendall test indicated that 18 percent of unregulated streamgages with (1) at least 30 years of peak-flow record and (2) peak-flow record at least as recent as water year 2000 showed statistically significant trends, most of which were weak to moderate increases in peak flows. Concurrently, the fundamental geospatial datasets that support the Kentucky StreamStats application were updated by using high-resolution digital elevation models and hydrography datasets to derive flow direction, flow accumulation, and stream definition rasters. Comparisons of regression-based flood-frequency models using the old and new layers demonstrated consistent results, with a statewide root-mean-square error of 0.019, in the base-10 logarithm of cubic feet per second. Furthermore, to assess model performance, flood-frequency estimates made by using the updated layers and previously published regression-based models were compared to flood-frequency estimates newly computed by following Bulletin 17C. This analysis showed the models performed adequately for most Kentucky stream locations. The updated statistics and geospatial layers provide stakeholders with more accurate, current data for flood-risk assessment, infrastructure design, and water-resource management.

Kentucky↗

Spatio-temporal ecological models via physics-informed neural networks for studying chronic wasting disease

To mitigate the negative effects of emerging wildlife diseases in biodiversity and public health it is critical to accurately forecast pathogen dissemination while incorporating relevant spatio-temporal covariates. Forecasting spatio-temporal processes can often be improved by incorporating scientific knowledge about the dynamics of the process using physical models. Ecological diffusion equations are often used to model epidemiological processes of wildlife diseases where environmental factors play a role in disease spread. Physics-informed neural networks (PINN) are deep learning algorithms that constrain neural network predictions based on physical laws and therefore are powerful forecasting models useful even in cases of limited and imperfect training data. In this paper, we develop a novel ecological modeling tool using PINNs, which fits a feedforward neural network and simultaneously performs parameter identification in a partial differential equation (PDE) with varying coefficients. We demonstrate the applicability of our model by comparing it with the commonly used Bayesian stochastic partial differential equation method and traditional machine learning approaches, showing that our proposed model exhibits superior prediction and forecasting performance when modeling chronic wasting disease in deer in Wisconsin. Furthermore, our model provides the opportunity to obtain scientific insights into spatiotemporal covariates affecting spread and growth of diseases. This work contributes to future machine learning and statistical methodology development by studying spatio-temporal processes enhanced by prior physical knowledge.

Spatial Statistics↗

Assessment of environments for Mars Science Laboratory entry, descent, and surface operations

The Mars Science Laboratory mission aims to land a car-sized rover on Mars' surface and operate it for at least one Mars year in order to assess whether its field area was ever capable of supporting microbial life. Here we describe the approach used to identify, characterize, and assess environmental risks to the landing and rover surface operations. Novel entry, descent, and landing approaches will be used to accurately deliver the 900-kg rover, including the ability to sense and "fly out" deviations from a best-estimate atmospheric state. A joint engineering and science team developed methods to estimate the range of potential atmospheric states at the time of arrival and to quantitatively assess the spacecraft's performance and risk given its particular sensitivities to atmospheric conditions. Numerical models are used to calculate the atmospheric parameters, with observations used to define model cases, tune model parameters, and validate results. This joint program has resulted in a spacecraft capable of accessing, with minimal risk, the four finalist sites chosen for their scientific merit. The capability to operate the landed rover over the latitude range of candidate landing sites, and for all seasons, was verified against an analysis of surface environmental conditions described here. These results, from orbital and model data sets, also drive engineering simulations of the rover's thermal state that are used to plan surface operations.

Space Science Reviews↗

Slow slip phenomena in Cascadia from 2007 and beyond: a review

Recent technological advances combined with more detailed analyses of seismologic and geodetic observations have fundamentally changed our understanding of the ways in which tectonic stresses arising from plate motions are accommodated by slip on faults. The traditional view that relative plate motions are accommodated by a simple cycle of stress accumulation and release on “locked” plate-boundary faults has been revolutionized by the serendipitous discovery and recognition of the significance of slow-slip phenomena, mostly in the deeper reaches of subduction zones. The Cascadia subduction zone, located in the Pacific Northwest of the conterminous United States and adjacent Canada, is an archetype of exploration and learning about slow-slip phenomena. These phenomena are manifest as geodetically observed aseismic transient deformations accompanied by a previously unrecognized class of seismic signals. Although secondary failure processes may be involved in generating the seismic signals, the primary origins of both aseismic and seismic phenomena appear to be episodic fault slip, probably facilitated by fluids, on a plate interface that is critically stressed or weakened. In Cascadia, this transient slip evolves more slowly and over more prolonged durations relative to the slip in earthquakes, and it occurs between the 30- and 40-km-depth contours of the plate interface where information was previously elusive. Although there is some underlying organization that relaxes nearly all the accrued plate-motion stresses along the entirety of Cascadia, we now infer that slow slip evolves in complex patterns indicative of propagating stress fronts. Our new understanding provides key constraints not only on the region where the slow slip originates, but also on the probable characteristics of future megathrust earthquakes in Cascadia. Herein, we review the most significant scientific issues and progress related to understanding slow-slip phenomena in Cascadia and highlight some of their societal implications. We provide a comprehensive review, from the big picture as inferred from studies of regional-scale monitoring data to the details revealed by innovative, focused experiments and new instrumentation. We focus on what has been learned largely since 2007, when several major investments in monitoring and temporary deployments dramatically increased the quality and quantity of available data.

Oregon↗

South Florida Ecosystem Program web site

Introduction The South Florida Ecosystem Restoration Program is an intergovernmental effort, involving a number of agencies, to reestablish and maintain the ecosystem of south Florida. One element of the restoration effort is the development of a firm scientific basis for resource decision making. The U.S. Geological Survey (USGS), one of the agencies, provides scientific information a part of the South Florida Ecosystem Restoration Program (SFEP). The USGS began their own program, called the South Florida Ecosystem Program, in fiscal year 1995 for the purpose of gathering hydrologic, cartographic, geologic, and biologic data that relate to the mainland of south Florida, Florida Bay, and the Florida Keys and Reef ecosystems.

Florida↗

Reconstructing historical changes in the environmental health of watersheds by using sediment cores from lakes and reservoirs in Salt Lake Valley, Utah

The Great Salt Lake Basins study area of the National Water-Quality Assessment (NAWQA) program, which began in 1997, is increasing the scientific understanding of factors that affect surface-water quality within the study-area boundaries (fig. 1). One way to improve the understanding of these factors is to look at historical trends in existing water-quality data. Unfortunately, short record lengths, in- consistent analytical methods, numerous measurements at less than detection levels, and questionable accuracy limit the usefulness of historical monitoring data for most trace inorganic and organic contaminants found in streams, rivers, and lakes in the study area.

Idaho, Utah, Wyoming↗

Topographic Map of the Northwest Ascraeus Mons Region of Mars - MTM 500k 15/252E OMKT

This map is part of a series of topographic maps of areas of special scientific interest on Mars. The topography was compiled photogrammetically using Viking Orbiter stereo image pairs. The contour interval is 250 meters. Horizontal and vertical control was established using the USGS Mars Digital Image Model 2.0 (MDIM 2.0) and data from the Mars Orbiter Laser Altimeter (MOLA).

IMAP↗

Topographic Map of the Northeast Ascraeus Mons Region of Mars - MTM 500k 15/257E OMKT

This map is part of a series of topographic maps of areas of special scientific interest on Mars. The topography was compiled photogrammetically using Viking Orbiter stereo image pairs. The contour interval is 250 meters. Horizontal and vertical control was established using the USGS Mars Digital Image Model 2.0 (MDIM 2.0) and data from the Mars Orbiter Laser Altimeter (MOLA).

IMAP↗

Topographic Map of the Southeast Ascraeus Mons Region of Mars - MTM 500k 10/257E OMKT

This map is part of a series of topographic maps of areas of special scientific interest on Mars. The topography was compiled photogrammetically using Viking Orbiter stereo image pairs. The contour interval is 250 meters. Horizontal and vertical control was established using the USGS Mars Digital Image Model 2.0 (MDIM 2.0) and data from the Mars Orbiter Laser Altimeter (MOLA).

IMAP↗

Topographic Map of the Southwest Ascraeus Mons Region of Mars - MTM 500k 10/252E OMKT

This map is part of a series of topographic maps of areas of special scientific interest on Mars. The topography was compiled photogrammetically using Viking Orbiter stereo image pairs. The contour interval is 250 meters. Horizontal and vertical control was established using the USGS Mars Digital Image Model 2.0 (MDIM 2.0) and data from the Mars Orbiter Laser Altimeter (MOLA).

IMAP↗

A model to locate potential areas for lake sturgeon spawning habitat construction in the St. Clair–Detroit River System

In response to a need for objective scientific information that could be used to help remediate loss of fish spawning habitat in the St. Clair River and Detroit River International Areas of Concern, this paper summarizes a large-scale geographic mapping investigation. Our study integrates data on two variables that many riverine fishes respond to in selecting where to spawn in these waters (water flow velocity and water depth) with available maps of the St. Clair–Detroit River System (SC–DRS). Our objectives were to locate and map these two physical components of fish habitat in the St. Clair and Detroit rivers and Lake St. Clair using a geographic information system (GIS) and to identify where, theoretically, fish spawning habitat could be remediated in these rivers. The target fish species to which this model applies is lake sturgeon (Acipenser fulvescens), but spawning reefs constructed for lake sturgeon in this system have been used for spawning by 17 species of fish. Our analysis revealed areas in each river that possessed suitable water velocity and depth for fish spawning and therefore could theoretically be remediated by the addition of rock-rubble substrate like that used at two previously remediated sites in the Detroit River at Belle Isle and Fighting Island. Results of our analysis revealed that only 3% of the total area of the SC–DRS possesses the necessary combination of water depth and high flow velocity to be indicated by the model as potential spawning habitat for lake sturgeon.

Michigan, Ontario↗

Dynamics and spatio-temporal variability of environmental factors in Eastern Australia using functional principal component analysis

This paper introduces a new technique in ecology to analyze spatial and temporal variability in environmental variables. By using simple statistics, we explore the relations between abiotic and biotic variables that influence animal distributions. However, spatial and temporal variability in rainfall, a key variable in ecological studies, can cause difficulties to any basic model including time evolution. The study was of a landscape scale (three million square kilometers in eastern Australia), mainly over the period of 19982004. We simultaneously considered qualitative spatial (soil and habitat types) and quantitative temporal (rainfall) variables in a Geographical Information System environment. In addition to some techniques commonly used in ecology, we applied a new method, Functional Principal Component Analysis, which proved to be very suitable for this case, as it explained more than 97% of the total variance of the rainfall data, providing us with substitute variables that are easier to manage and are even able to explain rainfall patterns. The main variable came from a habitat classification that showed strong correlations with rainfall values and soil types. ?? 2010 World Scientific Publishing Company.

Journal of Biological Systems↗