Spatial Modeling of Iron Transformations Within Artificial Soil Aggregates

Structured soils exhibit significant variations in transport characteristics at the aggregate scale. Preferential flow occurs through macropores while predominantly diffusive exchange takes place in intra-aggregate micropores. Such environments characterized by mass transfer limitations are conducive to the formation of small-scale chemical gradients and promote strong spatial variation in processes controlling the fate of redox-sensitive elements such as Fe. In this study, we present a reactive transport model used to spatially resolve iron bioreductive processes occurring within a spherical aggregate at the interface between advective and diffusive domains. The model is derived from current conceptual models of iron(hydr)oxide (HFO) transformations and constrained by literature and experimental data. Data were obtained from flow-through experiments on artificial soil aggregates inoculated with Shewanella putrefaciens strain CN32, and include the temporal evolution of the bulk solution composition, as well as spatial information on the final solid phase distribution within aggregates. With all iron initially in the form of ferrihydrite, spatially heterogeneous formation of goethite/lepidocrocite, magnetite and siderite was observed during the course of the experiments. These transformations were reproduced by the model, which ascribes a central role to divalent iron as a driver of HFO transformations and master variable in the rate laws of the considered reaction network. The predicted dissolved iron breakthrough curves also match the experimental ones closely. Thus, the computed chemical concentration fields help identify factors governing the observed trends in the solid phase distribution patterns inside the aggregate. Building on a mechanistic description of transformation reactions, fluid flow and solute transport, the model was able to describe the observations and hence illustrates the importance of small-scale gradients and dynamics of bioreductive processes for assessing bulk iron cycling. As HFOs are ubiquitous in soils, such process-level understanding of aggregate-scale iron dynamics has broad implications for the prediction of the subsurface fate of nutrients and contaminants that interact strongly with HFO surfaces.