Statistical modelling and analysis of soil physicochemical properties influencing gully development across contrasting geomorphic settings in Anambra State, Nigeria.
# Quantitative Geochemical and Geotechnical Assessment of Gully Erosion Susceptibility in the Niger Delta Basin, Southeastern Nigeria
## Overview
This repository contains the dataset, statistical analysis workflow, and machine learning models for the geoscientific research study evaluating the intrinsic geochemical, mineralogical, and geomechanical controls driving severe gully erosion in the Niger Delta Basin, Southeastern Nigeria.
By contrasting baseline control terrains in the Imo Shale Formation (**Awka**) with active gully corridors in the Nanka Sands Formation (**Idemili**), this study establishes that gully vulnerability is fundamentally governed by a geochemical trade-off between unbonded quartzic framework grains and cohesive clay-sesquioxide binders ($\text{Fe}_2\text{O}_3$, $\text{Al}_2\text{O}_3$, $\text{MgO}$).
---
## Research Workflow & Methodology
The analytical pipeline follows a multi-stage quantitative workflow implemented entirely in Python:
1. **Laboratory Data Extraction & Structuring**:
* **Geochemical Screening**: Major elemental oxides ($\text{SiO}_2$, $\text{Al}_2\text{O}_3$, $\text{Fe}_2\text{O}_3$, $\text{MgO}$, $\text{TiO}_2$, $\text{K}_2\text{O}$, $\text{Na}_2\text{O}$, $\text{CaO}$) and Silica-to-Sesquioxide Ratios ($\text{SSR}$).
* **Mineralogical Scoring**: Quantitative Mineral Scores for Quartz, Kaolinite, Smectite, Goethite, and Hematite.
* **Hydro-Textural Testing**: Particle Size Distribution (Sand, Silt, Clay fractions) and Atterberg Limits (Liquid Limit, Plastic Limit, Plasticity Index).
2. **Multivariate Pattern Recognition**:
* Principal Component Analysis (PCA) for variance decomposition and dimensional reduction.
* Unsupervised Ward's Hierarchical Clustering evaluated via Adjusted Rand Index (ARI).
3. **Predictive Diagnostic Modeling**:
* $L_2$-Regularized Binary Logistic Regression to classify stable vs. active gully terrains based on core geochemical predictors.
* Performance validation using Receiver Operatin …