Here’s a detailed description for updating and improving your crop recommendation system based on soil data:

Description of a Crop Recommendation System with Soil Data

A crop recommendation system helps farmers select the best crops to grow based on the specific properties of their soil. This system uses soil characteristics and environmental factors to determine the crops that are most likely to thrive. Recommendations are provided to improve crop yield, optimize resource use, and ensure sustainable farming practices.

Core Components for Recommendations

The system should consider the following soil parameters and external factors to make accurate recommendations:

1. Soil Nutrients:

   • Nitrogen (N): Promotes leafy growth; ideal for crops like spinach, lettuce, and wheat.
   • Phosphorus (P): Essential for root development; crucial for legumes, peas, and root vegetables like carrots.
   • Potassium (K): Enhances disease resistance and fruit quality; important for fruiting plants like tomatoes, bananas, and potatoes.

2. Soil pH:

   • Indicates soil acidity or alkalinity.
   • Neutral pH (6.5-7.5) supports most crops like rice, wheat, and maize.
   • Acidic soil (<6.5) favors crops like tea and coffee.
   • Alkaline soil (>7.5) supports crops like barley and asparagus.

3. Organic Matter:

   • High organic content improves water retention and nutrient availability.
   • Crops like vegetables and fruits benefit from rich organic matter.

4. Moisture Level:

   • Determines irrigation needs and crop suitability.
   • High moisture crops: Paddy, sugarcane.
   • Low moisture crops: Millet, sunflower.

5. Temperature:

   • Warm crops: Maize, rice, and cotton.
   • Cool crops: Wheat, barley, and cabbage.

6. Rainfall:

   • Rain-fed crops (e.g., rice) thrive in high rainfall areas.
   • Drought-resistant crops (e.g., millets) perform well in low-rainfall zones.

7. Geographical Factors:

   • Altitude, latitude, and local climate conditions.
   • Example: Coffee grows well in high altitudes, while coconut thrives in coastal regions.

How to Update Recommendations

1. Dynamic Soil Profiles:

   • Use real-time soil testing data to determine nutrient levels, pH, and moisture.
   • Example: If the nitrogen level is low, recommend nitrogen-fixing crops like legumes.

2. Crop Rotation Insights:

   • Suggest crop rotations to maintain soil health.
   • Example: After a nitrogen-depleting crop like wheat, recommend a nitrogen-fixing crop like lentils.

3. Fertilizer Suggestions:

   • Provide recommendations for fertilizers based on deficiencies.
   • Example: If phosphorus is low, suggest adding rock phosphate.

4. Weather and Climate Integration:

   • Include real-time weather data like rainfall forecasts and temperature trends.
   • Example: Recommend drought-tolerant crops during dry seasons.

5. Regional Crop Suitability:

   • Use regional data to match crops with local soil and climate.
   • Example: Recommend paddy in water-rich regions like Punjab, and millet in arid regions like Rajasthan.

Sample Output for Crop Recommendations

Based on soil and environmental data: - Soil Parameters: - pH: 6.8 (neutral) - Nitrogen: Medium - Phosphorus: Low - Potassium: High - Moisture: Moderate - Recommendations: - Primary Crops: Wheat, Maize, Barley. - Secondary Crops (Improving Soil Health): Lentils, Chickpeas (for nitrogen fixation). - Fertilizer Recommendation: Use phosphorus-rich fertilizers (e.g., DAP).

How to Present Recommendations

• Use a dashboard or mobile app for farmers.
• Show clear visualizations of soil test results and matched crops.
• Include:
  • Top recommended crops.
  • Fertilizer and irrigation tips. -``

Soil Classification Datasets

Please ensure to cite the paper when utilizing the dataset in a research study. Refer to the paper link or BibTeX provided below.

This repository contains comprehensive datasets for soil classification and recognition research. The Original Dataset comprises soil images sourced from various online repositories, which have been meticulously cleaned and preprocessed to ensure data quality and consistency. To enhance the dataset's size and diversity, we employed Generative Adversarial Networks (GANs), specifically the CycleGAN architecture, to generate synthetic soil images. This augmented collection is referred to as the CyAUG Dataset. Both datasets are specifically designed to advance research in soil classification and recognition using state-of-the-art deep learning methodologies.

This dataset was curated as part of the research study titled "An advanced artificial intelligence framework integrating ensembled convolutional neural networks and Vision Transformers for precise soil classification with adaptive fuzzy logic-based crop recommendations" by Farhan Sheth, Priya Mathur, Amit Kumar Gupta, and Sandeep Chaurasia, published in Engineering Applications of Artificial Intelligence.

Links

• Paper: An advanced artificial intelligence framework integrating ensembled convolutional neural networks and Vision Transformers for precise soil classification with adaptive fuzzy logic-based crop recommendations
• Code: Development Code (model and recommendation system)

Application produced by this research is available at:

• Agro Companion: Agro Companion
• Server Code: Agro Companion Server

Note: If you are using any part of this project; dataset, code, application, then please cite the work as mentioned in the Citation section below.

Dataset

Both dataset consists of images of 7 different soil types.

The Soil Classification Dataset is structured to facilitate the classification of various soil types based on images. The dataset includes images of the following soil types:

• Alluvial Soil
• Black Soil
• Laterite Soil
• Red Soil
• Yellow Soil
• Arid Soil
• Mountain Soil

The dataset is organized into folders, each named after a specific soil type, containing images of that soil type. The images vary in resolution and quality, providing a diverse set of examples for training and testing classification models.

Original Dataset Details

• Total Images: 1189 images
• Image Format: JPG/JPEG
• Image Size: Varies
• Source: Collected from various online repositories and cleaned for consistency.

CyAUG Dataset Details

• Total Images: 5097 images
• Image Format: JPG/JPEG
• Image Size: Varies
• Source: Generated using CycleGAN to augment the original dataset, enhancing its size and diversity.

Input and Output Parameters

• Input Parameters:
  • Image: The images of the soils (JPG/JPEG format).
  • Label: The labels are in the format 'soil types' (folder names).
• Output Parameter:
  • Classification: The predicted class (soil type) based on the input image.

Citation

If you are using any of the derived dataset, please cite the following paper:

@article{SHETH2025111425,
  title = {An advanced artificial intelligence framework integrating ensembled convolutional neural networks and Vision Transformers for precise soil classification with adaptive fuzzy logic-based crop recommendations},
  journal = {Engineering Applications of Artificial Intelligence},
  volume = {158},
  pages = {111425},
  year = {2025},
  issn = {0952-1976},
  doi = {https://doi.org/10.1016/j.engappai.2025.111425},
  url = {https://www.sciencedirect.com/science/article/pii/S0952197625014277},
  author = {Farhan Sheth and Priya Mathur and Amit Kumar Gupta and Sandeep Chaurasia},
  keywords = {Soil classification, Crop recommendation, Vision transformers, Convolutional neural network, Transfer learning, Fuzzy logic}
}

We developed a comprehensive, gridded Global Soil Dataset for use in Earth System Models (GSDE) and other applications as well. GSDE provides soil information including soil particle-size distribution, organic carbon, and nutrients, etc. and quality control information in terms of confidence level. GSDE is based on the Soil Map of the World and various regional and national soil databases, including soil attribute data and soil maps. We used a standardized data structure and data processing procedures to harmonize the data collected from various sources. We then used a soil type linkage method (i.e. taxotransfer rules) and the polygon linkage method to derive the spatial distribution of soil properties. To aggregate the attributes of different compositions of a mapping unit, we used three mapping approaches: area-weighting method, the dominant soil type method and the dominant binned soil attribute method. In the released gridded dataset, we used the area-weighting method as it will meet the demands of most applications. The dataset can be also aggregate to a lower resolution. The resolution is 30 arc-seconds (about 1 km at the equator). The vertical variation of soil property was captured by eight layers to the depth of 2.3 m (i.e. 0- 0.045, 0.045- 0.091, 0.091- 0.166, 0.166- 0.289, 0.289- 0.493, 0.493- 0.829, 0.829- 1.383 and 1.383- 2.296 m).

This dataset is the cleaned up version of the "Soil Image Dataset"(https://www.kaggle.com/datasets/jayaprakashpondy/soil-image-dataset), which had lots of corrupted images. The dataset contains 1555 images divided into two subsets (train and test set) of 4 classes each.

(Link to Metadata) This data set is a digital soil survey and generally is the most detailed level of soil geographic data developed by the National Cooperative Soil Survey. The information was prepared by digitizing maps, by compiling information onto a planimetric correct base and digitizing, or by revising digitized maps using remotely sensed and other information. This data set consists of georeferenced digital map data and computerized attribute data. The map data are in a soil survey area extent format and include a detailed, field verified inventory of soils and miscellaneous areas that normally occur in a repeatable pattern on the landscape and that can be cartographically shown at the scale mapped. A special soil features layer (point and line features) is optional. This layer displays the location of features too small to delineate at the mapping scale, but they are large enough and contrasting enough to significantly influence use and management. The soil map units are linked to attributes in the National Soil Information System relational database, which gives the proportionate extent of the component soils and their properties. Survey Dates - https://www.nrcs.usda.gov/wps/portal/nrcs/surveylist/soils/survey/state/?stateId=VT

Description

In this dataset, there are soil data analyses with properties such as pH, organic matter (OM), salinity (EC), etc., major elements (N, P, K, Mg) as well as some microelements (Fe, Zn, Mn, Cu, B) with significant impact on plant nutrition.

Categories

Agricultural Soil

Acknowledgements & Source

Panagiotis Tziachris

Data Source

View Details

Image Source

This layer shows soil type, based on the result of a classification established from Kalimantan RePPProT data on 'SL_ORDER' field (1990, 1:250,000 scale) . This data was provided and processed by Daemeter Consulting. Soil categories from RePPProT were then re-classified by the World Resources Institute according to the FAO Digital Soil Map of the World, for use in the Suitability Mapper (2012). The FAO data is available at http://www.fao.org/geonetwork/srv/en/metadata.show?id=14116 . Data separated into categories: Inceptisol; Oxisol; Alfisol; Ultisol; Spodosol; Entisol; Histosol.

This dataset is a digital soil survey and generally is the most detailed level of soil geographic data developed by the National Cooperative Soil Survey. The information was prepared by digitizing maps, by compiling information onto a planimetric correct base and digitizing, or by revising digitized maps using remotely sensed and other information.

This dataset consists of georeferenced digital map data and computerized attribute data. The map data are in a soil survey area extent format and include a detailed, field verified inventory of soils and miscellaneous areas that normally occur in a repeatable pattern on the landscape and that can be cartographically shown at the scale mapped. A special soil features layer (point and line features) is optional. This layer displays the location of features too small to delineate at the mapping scale, but they are large enough and contrasting enough to significantly influence use and management. The soil map units are linked to attributes in the National Soil Information System relational database, which gives the proportionate extent of the component soils and their properties.

Note: This metadata record was created by MnGeo to serve as a generic record for all SSURGO data sets within Minnesota. See the individual county metadata records created by NRCS for county-specific information; these records are included in the data set download files.

This data set provides gridded data for selected soil parameters derived from data and methods developed by the Global Soil Data Task, an international collaborative project with the objective of making accurate and appropriate data relating to soil properties accessible to the global change research community. The task was coordinated by the International Geosphere-Biosphere Programme (IGBP-DIS). The data in this data set were produced by the International Satellite Land-Surface Climatology Project, Initiative II (ISLSCP II) staff from data obtained from the Oak Ridge National Laboratory Distributed Active Archive Center (ORNL DAAC, http://daac.ornl.gov/). See the related data sets section below. Two-dimensional gridded maps of selected soil parameters, including soil texture, at a 1.0 by 1.0 degree spatial resolution and for two soil depths are provided. All data layers have been adjusted to match the ISLSCP II land/water mask. There are 36 data files with this data set.

The data were derived from scanned soil maps. Attributes include a soil number (2-180), corresponding to runoff coefficient values in a Hydrology Manual, provided by the Los Angeles County Department of Public Works, Water Resources Division.Purpose: For use in DPW’s Modified Rational Method Hydrology Model.Supplemental Information:Stormwater Engineering is a Division of the Los Angeles County Department of Public Works. Please visit their website for posted publications, including the above mentioned Hydrology Manual.

The Harmonized World Soil Database version 2.0 (HWSD v2.0) is a unique global soil inventory providing information on the morphological, chemical and physical properties of soils at approximately 1 km resolution. Its main objective is to serve as a basis for prospective studies on agro-ecological zoning, food security and climate change. The Harmonized World Soil Database (HWSD) was established in 2008 by the International Institute for Applied Systems Analysis (IIASA) and FAO, and in partnership with International Soil Reference and Information Centre (ISRIC), the European Soil Bureau Network (ESBN) and the Institute for Soil Sciences Chinese Academy of Sciences (CAS). The data entry and harmonization within a Geographic Information System (GIS) was carried out at IIASA, with verification of the database undertaken by all partners. HWSD was then updated in 2013 (HWSD v1.2) and in 2023 (HWSD v2.0). This updated version (HWSD v2.0) is built on the previous versions of HWSD with several improvements on (i) the data source that now includes several national soil databases, (ii) an enhanced number of soil attributes available for seven soil depth layers, instead of two in HWSD v1.2, and (iii) a common soil reference for all soil units (FAO1990 and the World Reference Base for Soil Resources). This contributes to a further harmonization of the database. The GIS raster image file is linked to the soil attribute database. The HWSD v2.0 soil attribute database provides information on the soil unit composition for each of the near 30 000 soil association mapping units. The HWSD v2.0 Viewer, provided with the database, creates this link automatically and provides direct access to the soil attribute data and the soil association information. Note: - A tutorial for accessing HWSD ver. 2.0 using R (prepared by David Rossiter, June 2023) has been added as an 'associated resource' (NOTE: Needs the SQLite version of HWSD v2 as provided below). - Soil property estimates in HWSDv2 were derived from Batjes (2016), Geoderma (https://doi.org/10.1016/j.geoderma.2016.01.034).

ISRIC World Soil Information is compiling legacy soil profile data of Sub Saharan Africa, as a project activity of the AfSIS project (Globally integrated Africa Soil Information Service). http://africasoils.net/services/data/soil-databases/ Africa Soil Profiles database, version. 1.0 (April 2012) identifies less than 15700 unique soil profiles inventoried from a wide variety of data sources. From the less than 14600 profiles that are geo-referenced, soil layer attribute data are available for less than 12500 and soil analytical data for less than 10000 profiles. The database includes, but is not limited, to the soil attributes specified by GlobalSoilMap.net. Soil attribute values are standardized according to e-SOTER conventions and validated according to routine rules. Odd values are flagged. The degree of validation, and associated reliability of the data, varies because reference soil profile data, that are previously and thoroughly validated, are compiled together with non-reference soil profile data of lesser inherent representativeness. Updated milestone versions of the dataset have been posted online and made available to the project serving as input to the soil property maps generated by AfSIS. The continuously growing dataset will also be made available through the World Soil Information Service upon continuation of the project activity. The version is released here is version 1.0., the latest version is 1.1.

The International Soil Reference and Information Centre-World Inventory of Soil Emission Potentials (ISRIC-WISE) international soil profile data set consists of a homogenized, global set of 1,125 soil profiles for use by global modelers. These profiles provided the basis for the Global Pedon Database (GPDB) of the International Geosphere-Biosphere Programme (IGBP) - Data and Information System (DIS). The data set consists of a selection of 665 profiles originating from the Natural Resources Conservation Service (NRCS, Lincoln), 250 profiles obtained from the Food and Agriculture Organization (FAO, Rome), and 210 profiles from the reference collection of the International Soil Reference and Information Centre (ISRIC, Wageningen). All profiles are georeferenced and classified according to the 1974 Legend of the FAO-UNESCO Soil Map (FAC-UNESCO, 1974) of the World, as well as the 1988 Revised Legend of FAO-UNESCO (FAO, 1990). The data set includes information on soil classification, site data, soil horizon data, source of data, and methods used for determining analytical data. The data files are in a comma-delimited format. Data Citation: The data set should be cited as follows: Batjes, N. H. (ed). 2000. Global Soil Profile Data (ISRIC-WISE). Available on-line from the ORNL Distributed Active Archive Center, Oak Ridge National Laboratory, Oak Ridge, Tennessee, U.S.A.

The National Soil Database has produced a national database of soil geochemistry including point and spatial distribution maps of major nutrients, major elements, essential trace elements, trace elements of special interest and minor elements. In addition, this study has generated a National Soil Archive, comprising bulk soil samples and a nucleic acids archive each of which represent a valuable resource for future soils research in Ireland. The geographical coherence of the geochemical results was considered to be predominantly underpinned by underlying parent material and glacial geology. Other factors such as soil type, land use, anthropogenic effects and climatic effects were also evident. The coherence between elements, as displayed by multivariate analyses, was evident in this study. Examples included strong relationships between Co, Fe, As, Mn and Cu. This study applied large-scale microbiological analysis of soils for the first time in Ireland and in doing so also investigated microbial community structure in a range of soil types in order to determine the relationship between soil microbiology and chemistry. The results of the microbiological analyses were consistent with geochemical analyses and demonstrated that bacterial community populations appeared to be predominantly determined by soil parent material and soil type.

A global data set of soil types is available at 0.5-degree latitude by 0.5-degree longitude resolution. There are 106 soil units, based on Zobler?s (1986) assessment of the FAO/UNESCO Soil Map of the World. This data set is a conversion of the Zobler 1-degree resolution version to a 0.5-degree resolution. The resolution of the data set was not actually increased. Rather, the 1-degree squares were divided into four 0.5-degree squares with the necessary adjustment of continental boundaries and islands. The computer code used to convert the original 1-degree data to 0.5-degree is provided as a companion file. A JPG image of the data is provided in this document. The Zobler data (1-degree resolution) as distributed by Webb et al. (1993) [http://www.ngdc.noaa.gov/seg/eco/cdroms/gedii_a/datasets/a12/wr.htm#top] contains two columns, one column for continent and one column for soil type. The Soil Map of the World consists of 9 maps that represent parts of the world. The texture data that Webb et al.(1993) provided allowed for the fact that a soil type in one part of the world may have different properties than the same soil in a different part of the world. This continent-specific information is retained in this 0.5-degree resolution data set, as well as the soil type information which is the second column. A code was written (one2half.c) to take the file CONTIZOB.LER distributed by Webb et al. (1993) [http://www.ngdc.noaa.gov/seg/eco/cdroms/gedii_a/datasets/a12/wr.htm#top] and simply divide the 1-degree cells into quarters. This code also reads in a land/water file (land.wave) that specifies the cells that are land at 0.5 degrees. The code checks for consistency between the newly quartered map and the land/water map to which the quartered map is to be registered. If there is a discrepancy between the two, an attempt was made to make the two consistent using the following logic. If the cell is supposed to be water, it is forced to be water. If it is supposed to be land but was resolved to water at 1 degree, the code looks at the surrounding 8 cells and picks the most frequent soil type and assigns it to the cell. If there are no surrounding land cells then it is kept as water in the hopes that on the next pass one or more of the surrounding cells might be converted from water to a soil type. The whole map is iterated 5 times. The remaining cells that should be land but couldn't be determined from surrounding cells (mostly islands that are resolved at 0.5 degree but not at 1 degree) are printed out with coordinate information. A temporary map is output with -9 indicating where data is required. This is repeated for the continent code in CONTIZOB.LER as well. A separate map of the temporary continent codes is produced with -9 indicating required data. A nearly identical code (one2half.c) does the same for the continent codes. The printout allows one to consult the printed versions of the soil map and look up the soil type with the largest coverage in the 0.5-degree cell. The program manfix.c then will go through the temporary map and prompt for input to correct both the soil codes and the continent codes for the map. This can be done manually or by preparing a file of changes (new_fix.dat) and redirecting stdin. A new complete version of the map is outputted. This is in the form of the original CONTIZOB.LER file (contizob.half) but four times larger. Original documentation and computer codes prepared by Post et al. (1996) are provided as companion files with this data set. Image of 106 global soil types available at 0.5-degree by 0.5-degree resolution. Additional documentation from Zobler?s assessment of FAO soil units is available from the NASA Center for Scientific Information.

This data set consists of soil texture classification data derived from field surveys as part of the Soil Moisture Active Passive Validation Experiment 2012 (SMAPVEX12). The soil texture classification map provides information about vegetation present in the study area.

Abstract

This dataset created by the Native Lands Advocacy Project contains spatial and tabular data derived from the NRCS's 2020 national SSURGO soils database and has been formatted to include only data within current American Indian areas for the conterminous United States. Boundary data was acquired from the US Census Bureau's Tiger Database (2010) reported by the federally recognized tribal governments through the Census Bureau's Boundary and Annexation Survey (BAS) and includes a total of 549 areas associated with 398 unique American Indian groups. The boundary data includes all Indian Lands defined by CFR 25 CFR § 502.12 which includes both lands held in trust by the United States Government and fee lands within the boundaries of American Indian Lands encompassing a total of 109,465,623 acres of land. The spatial data includes each soil unit and its associated MUKEY making it possible to join this data with other GSSURGO soil characteristics tables.According to the NRCS4, "The SSURGO database provides the most detailed level of information and was designed primarily for farm and ranch, landowner/user, township, county, or parish natural resource planning and management. Using the soil attributes, this database serves as an excellent source for determining erodible areas and developing erosion control practices; reviewing site development proposals and land use potential; making land use assessments and chemical fate assessments; and identifying potential wetlands and sand and gravel aquifer areas." The NRCS's SSURGO data is used in numerous land valuation, carbon and hydrologic assessment models including the proprietary AcreValue™ valuation estimation tool5, NRCS's Rapid Carbon Assessment RaCa6, the EPA's Automated Geospatial Watershed Assessment (AGWA), and the Soil and Water Assessment Tool (SWAT), to name a few.

This data set is a digital soil survey and generally is the most detailed level of soil geographic data developed by the National Cooperative Soil Survey (SSURGO). The information was prepared by digitizing maps, by compiling information onto a planimetric correct base and digitizing, or by revising digitized maps using remotely sensed and other information. The map units delineated on the detailed soil maps in a soil survey represent the soils or miscellaneous areas in the survey area. The map unit descriptions in this report, along with the maps, can be used to determine the composition and properties of a unit. A map unit delineation on a soil map represents an area dominated by one or more major kinds of soil or miscellaneous areas. A map unit is identified and named according to the taxonomic classification of the dominant soils. Within a taxonomic class there are precisely defined limits for the properties of the soils. On the landscape, however, the soils are natural phenomena, and they have the characteristic variability of all natural phenomena. Thus, the range of some observed properties may extend beyond the limits defined for a taxonomic class. Areas of soils of a single taxonomic class rarely, if ever, can be mapped without including areas of other taxonomic classes. Consequently, every map unit is made up of the soils or miscellaneous areas for which it is named and some minor components that belong to taxonomic classes other than those of the major soils.The Map Unit Description (Brief, Generated) report displays a generated description of the major soils that occur in a map unit. Descriptions of non-soil (miscellaneous areas) and minor map unit components are not included. This description is generated from the underlying soil attribute data. To see the Non-Technical description of the soil types, click here.

For more information, contact: GIS Manager Information Technology & Innovation (ITI) Montgomery County Planning Department, MNCPPC T: 301-650-5620 U.S. Department of Agriculture USDA Natural Resources Conservation Service p: 1-833-ONE-USDA e: askusda@usda.gov

The Food and Agriculture Organization of the United Nations (FAO) and the International Institute for Applied Systems Analysis (IIASA) combined the available regional and national soil information with the data already contained within the 1:5,000,000 scale FAO-UNESCO map, into a new comprehensive Harmonized World Soil Database (HWSD_v121). This map has a resolution of about 1 km (30 arc seconds) and consists of a 30-cm topsoil layer, and a 70-cm subsoil layer. The soil variables provided in the Harmonized World Soil Database (2009) and FAO/UNESCO Soil Map of the World included soil texture (%sand, %silt, %clay), organic carbon, pH, and EC. However, from a hydrological point of view, we are in need of parameters such as bulk density, water storage capacity, and hydraulic conductivity for different soil layers. Hence, we have used various pedotransfer functions from the literature to estimate the soil parameters needed in a Soil and Water Assessment Tool (SWAT model). The associated SWAT2012.mdb and lookup table is available at 2w2e GmbH website. […]

Forty two samples were selected from the Kellogg Soil Survey Laboratory (KSSL) archive. The soils (41) were taken from the A horizon except for one sample that came from an O horizon. The samples represented 9 of the 12 US soil Orders, including Mollisols (23), Alfisols (5), Ultisols (5), Andisols (2), Entisols (2), Inceptisols (2), Aridisols (1), Histosols (1) and Vertisols (1). The soils varied widely in SOC (3.0 – 288.4 g kg-1; mean 31 g kg-1), pH (4.3 – 8.5; mean 6.2) and clay content (3.6 – 47.0%; mean 21.5%) The geographic origin of the selected samples and the distribution of SOC concentrations, clay contents and pH values are in the sample selected materials. All samples were analyzed using methods of the KSSL Manual (Soil Survey Staff 2014). Samples were dried in an open tray in a low temperature oven at 30-35 °C for 3-7 days and ground to