Processes and rates of formation defined by modelling in alkaline to acidic soil systems in Brazilian Pantanal wetland

Within a soil system, it is generally difficult to relate the observed secondary phases and their spatial distribution to past or current processes. Geochemical modelling can help to identify present-days processes and to quantify formation rates if the system is sufficiently constrained. Here we used the PHREEQC model for this purpose, applying it to soil alkaline to acidic soil systems in the Nhecolandia area. A large amount of data was available, but scientific questions remained opened. We carried out a complementary field study to clarify a key point necessary for modelling, confirming the disconnection between two types of groundwater, perched, alkaline and deep, acidic. Our simulations showed that the observed soil system corresponds to the current biogeochemistry and that paleoclimatic genetic hypotheses are unnecessary. We found that the ridge soil formation rate ranges from 396 to 638 kg ha-1 y-1, depending on hypotheses. We showed that the downslope soils act as buffers to store labile species during the dry season, allowing alkalinity to be maintained from one year to the next. At least 125 y were necessary to obtain downslope alkaline lakes, and this alkaline nature can quickly disappear after changing drainage conditions. We explained how hyperacidic horizons (pH < 4) can neighbor alkaline horizons (pH > 10) with a sharp contact. When correctly informed, geochemical modelling appeared as an effective tool to answer questions about soil processes.

Automated summary

Geochemical modelling is an effective tool for studying soil processes and answering scientific questions. In this study, we used the PHREEQC model to investigate alkaline to acidic soil systems in the Nhecolandia area of Brazil. Our simulations provided insights into soil formation rates, the stability of alkaline soils, and the interactions between different soil layers.