Coastal wetlands are celebrated as blue carbon powerhouses, capable of storing up to roughly 200 megagrams of carbon per hectare in their top metre of soil. Crucially, about 60 to 70 percent of that carbon sits below 30 centimetres, where waterlogged, oxygen-poor conditions slow decomposition and iron minerals bind organic molecules into stable complexes. Most restoration assessments, however, monitor only vegetation recovery and the top 30 centimetres of soil, implicitly treating everything beneath as inert. The Hangzhou Bay team set out to test whether that assumption holds when management interventions dramatically alter salinity and moisture at the surface.

The researchers established three field conditions in a wetland overrun since the 1980s by the invasive cordgrass Spartina alterniflora. Untouched invaded stands served as the control. In a second set of plots, crews removed all aboveground vegetation and covered the soil with 0.1-millimetre black polyethylene film to suppress regrowth and curb evaporative salt accumulation, a technique known as plastic mulching. In a third set, vegetation was removed and the soil mechanically tilled to 160 centimetres with a rotary cultivator to boost aeration and physical disruption, a technique called deep tillage. Five replicate plots of 20 by 20 metres were established for each treatment, separated by at least 50 metres to minimise spatial interference.

Eighteen months later, the team collected soil cores from each plot and sliced them into five depth intervals spanning 0 to 100 centimetres. The samples were analysed for pH, electrical conductivity, moisture, nitrogen pools, total and organic carbon, microbial biomass, and a suite of iron fractions extracted through sequential chemical procedures. The researchers also quantified iron-bound organic carbon, the fraction of soil carbon locked to iron oxides, and sequenced bacterial 16S rRNA genes from surface layers to characterise community shifts. Statistical models, including partial least squares path modeling, were then used to explore how the measured variables related to one another across depth.

The most striking finding was the sheer reach of plastic mulching. Salinity fell by 43 to 53 percent throughout the entire profile, and electrical conductivity dropped by 45 to 50 percent in deep soils compared with untouched stands. Deep tillage, by contrast, reduced surface salinity by only about 20 percent and had limited effect below. With the osmotic stress of salt lifted, microbial biomass carbon in deep soils under mulching climbed from roughly 25 percent above control levels at 30 to 50 centimetres to more than 100 percent at 50 to 100 centimetres. Meanwhile, microbial biomass in the surface 30 centimetres declined, revealing a wholesale redistribution of microbial life from the surface toward deeper horizons.

That microbial awakening coincided with substantial carbon losses. Total carbon under mulching fell by 19 to 35 percent across the profile relative to the control, and soil organic carbon declined by 34 to 65 percent, with the strongest losses at depth. Total nitrogen dropped even more sharply, by 51 to 74 percent. The carbon declines tracked a parallel weakening of mineral protection: poorly crystalline iron oxides, the highly reactive phases that scavenge organic matter, fell by 30 to 50 percent, and iron-bound organic carbon fell by 35 to 50 percent. Integrated across the full metre of soil, the fixed-depth estimate pointed to an apparent soil carbon stock decline of 65 plus or minus 12 megagrams of carbon per hectare over just 18 months.

The mechanistic story the authors propose centres on what soil scientists call the iron gate. In anaerobic wetland soils, iron oxides form organo-mineral complexes that resist microbial attack, effectively shielding carbon for decades or centuries. Rapid desalinization can alter the solubility and mobility of reactive iron phases through changes in ionic strength and pH, potentially dissolving those protective associations. At the same time, salt-stressed deep soil microbes, once freed from osmotic constraint, may become more active in carbon turnover. The path models were consistent with this picture: mulching was strongly associated with altered soil physicochemical conditions, which in turn were negatively associated with reactive iron pools and with soil organic and iron-bound carbon. In surface soils, bacterial community attributes were linked to carbon retention, while in deep soils, increased microbial biomass covaried with reactive iron decline and lower carbon.

The bacterial sequencing revealed dramatic community reorganisation in the surface layers. Untouched stands harboured elevated abundances of sulfur- and iron-cycling taxa such as Sulfurifustis and Acidibacter, tilled soils were enriched in Nocardioides and Nitrospira, and mulched soils were characterised by Thiobacillus and other taxa. Mantel tests showed that salinity and iron-carbon variables were the strongest correlates of community composition, explaining a large share of the variation. The authors are careful, however, to note the limits of their evidence: sequencing covered only the top 20 centimetres, and no direct measurements of carbon dioxide or methane fluxes, pore-water dissolved organic carbon, iron redox speciation, or lateral hydrological export were made. The lost carbon could have been mineralised to greenhouse gases, leached away as dissolved iron-organic complexes, or simply redistributed within the profile, and the study cannot distinguish among these fates.

Those caveats matter because the numbers are sobering when set against restoration gains. Aboveground biomass in restored coastal wetlands typically accumulates at only about 1 to 3 megagrams of carbon per hectare per year after native vegetation returns. A potential subsoil loss of 65 megagrams over 18 months could therefore dwarf the early carbon benefits of replanting, at least in the short term. Three trajectories remain possible for mulched sites: stabilisation as labile carbon pools are exhausted, deceleration as recovering vegetation rebuilds carbon inputs, or continued decline if weakened iron protection persists. Only multi-year monitoring with direct flux measurements can determine which path these sites actually follow, and whether they function as net atmospheric carbon sources.

The study also carries an important design lesson. Because both managed treatments involved vegetation removal before mulching or tillage, the contrasts reflect responses to integrated field management states rather than isolated single-factor effects, and the single-site results from Hangzhou Bay may not generalise to wetlands with different salinity, tides, or sediment textures. Still, the central message is hard to ignore: the magnitude and spatial extent of hydro-geochemical change, not disturbance intensity alone, appears to govern the fate of deep soil carbon after invasive plant removal. The authors argue that blue carbon accounting frameworks should incorporate full-profile soil monitoring to at least one metre, iron-associated carbon metrics as indicators of mineral protection, and multi-year measurements of gaseous and dissolved carbon fluxes. As coastal restoration accelerates worldwide in the fight against climate change, the invisible carbon beneath the marsh surface may prove to be the most important number nobody was counting.

Subject of Research: Effects of plastic mulching and deep tillage restoration interventions on deep-soil salinity, microbial biomass, and iron-bound carbon in a Spartina alterniflora-invaded coastal wetland

Article Title: Profile-wide desalinization is associated with increased deep-soil microbial biomass and reduced iron-bound carbon in coastal wetland restoration

Article References: Profile-wide desalinization is associated with increased deep-soil microbial biomass and reduced iron-bound carbon in coastal wetland restoration. (n.d.). https://doi.org/10.5194/bg-23-6671-2026

Image Credits: AI Generated

DOI: 10.5194/bg-23-6671-2026

Keywords: coastal wetlands, blue carbon, soil organic carbon, iron-bound carbon, Spartina alterniflora, soil salinity, microbial biomass, wetland restoration, deep tillage, plastic mulching, soil biogeochemistry, Hangzhou Bay

The stakes are considerable. Tillage is among the most energy-intensive operations in crop production, and improper machinery use during land preparation can irreversibly damage soil structure. Deep ploughing demands enormous power to fracture compacted horizons, which is why manufacturers have developed machines that pair a passive subsoiler shank with an active, powered rotary tiller. The rotary component pulverizes the soil the subsoiler cracks open, mixes in fertilizer, crop residues and organic amendments, and reduces the number of field passes needed. Fewer passes mean less compaction, less fuel, less labor and less time in the field. But until now, designers have lacked a reliable way to predict exactly what such combined implements do to the soil profile before the steel ever touches the ground.

DEM offers that window. The technique, pioneered in the geotechnical literature since the late 1970s, treats soil not as a continuous medium but as millions of individual particles, solving the equations of motion for each one as they collide, slide and cohere. In this study, the team modeled the soil bed in Altair EDEM 2022 software using the hysteretic spring with linear cohesion contact model, which adds a cohesion resistance to the normal contact forces between particles to mimic the behavior of real clay. The implement itself was designed in PTC Creo Parametric 3D modeling software and imported into the simulation environment, allowing the researchers to replicate the exact geometry of the subsoiler-cum-rotary tiller used in physical tests.

Calibration was the critical step. Because no universally accepted procedure exists for selecting the micromechanical parameters that govern particle interactions, the team anchored their model to laboratory measurements. They worked with clay soil from the experimental soil bin at Nanjing Agricultural University, containing 47.0 percent clay, 36.5 percent silt and 15.5 percent sand, with an average bulk density of 1.5 grams per cubic centimeter. Direct shear tests under unconsolidated undrained conditions yielded cohesion, internal friction angle and shear strength via the Mohr-Coulomb equation, while a digital penetrometer measured penetration resistance at ten locations. Using the angle of repose method and reverse parametrization, the researchers iteratively adjusted restitution and friction coefficients until the simulated soil formed a natural pile at 36.87 degrees, matching the physical material’s behavior.

The validation results are striking. Across eighteen soil bin experiments, the mean horizontal tillage resistance measured by sensors on the three-point linkage of a motorized trolley was 1,950.474 newtons, while vertical resistance averaged 295.92 newtons. The DEM model reproduced these values with a relative error of just 0.4 percent for horizontal resistance and 4.9 percent for vertical resistance. Horizontal resistance predictions achieved a coefficient of determination of 0.9997 with a normalized root mean square error of 0.04, while vertical resistance reached an R-squared of 0.9. An unpaired t-test found no statistically significant difference between simulated and measured values for either component, meaning the virtual experiment was statistically indistinguishable from the real one.

The model also predicted the shape of the furrow left behind. Soil profilometer measurements of the tilled profile matched the simulation with a relative error of 4.3 percent, an R-squared of 0.9936 and a normalized RMSE of 0.23. Both experiment and simulation produced U-shaped furrows with loose soil at the bottom, a profile considered favorable for seed coverage. The average top width of the soil profile was 0.405 meters in the soil bin versus 0.3875 meters in the simulation, and bottom widths were 0.26 and 0.25 meters respectively, a close correspondence that underscores the model’s fidelity.

Perhaps the most consequential finding concerns soil mixing. By coloring the simulated topsoil and subsoil particles differently, the researchers could watch, in effect, inside the soil as the machine passed through. The L-shaped rotary blades promoted predominantly horizontal and downward soil movement with limited upward throw, allowing topsoil particles to migrate toward deeper layers without wholesale inversion of the profile. Even at a tillage depth of 30 centimeters, the natural layering remained discernible, a critical advantage over conventional rotary systems that churn the profile indiscriminately. The mixing of the two layers increased with cutting depth, as expected given the larger soil volume involved, but the vertical-axis rotation design avoided dragging subsoil to the surface, preserving stratification that underpins fertility and moisture retention.

The study also mapped how operating parameters drive energy demand. Both horizontal and vertical resistance rose as tillage depth increased from 0.15 to 0.30 meters, because deeper operation cuts, disperses and moves a greater volume of soil. Resistance likewise grew when forward speed increased from 1.5 to 2.5 kilometers per hour, since faster-moving soil particles gain acceleration, raising normal loads on the tool and thus frictional resistance. The substantial gap between horizontal and vertical forces, roughly 1,800 to 2,060 newtons versus 264 to 323 newtons, shows that far more energy is spent dragging the implement forward than penetrating downward. Notably, the rotary tiller’s rotational motion generated a forward thrust component that offset part of the draft requirement, one reason combined active-passive implements outperform conventional systems in energy efficiency.

The practical implications reach from the design office to the farm field. For manufacturers, the validated model provides a virtual testing platform to optimize blade shape, operating depth and rotational speed across diverse soil conditions without costly physical prototyping. The evidence points toward vertical-axis rotary tillers equipped with L-shaped blades as the configuration of choice for deep tillage that respects soil architecture. For farmers, the combined implement approach promises reduced operational time and fuel consumption while maintaining soil quality, addressing economic and environmental goals simultaneously. As the authors conclude, DEM can serve as an accurate, consistent and fast method for predicting the final soil condition and the resistances required for tillage operations, and with proper blade selection, deep tillage and soil structure conservation need no longer be opposing goals.

Subject of Research: Discrete element simulation of rotary mixing effects on soil structure and tillage resistance for a combined subsoiler and rotary tiller

Article Title: Evaluation of rotary mixing effects on soil structure and tillage resistance using discrete element method

Article References: Makange, N. R., & Ji, C. (2026). Evaluation of rotary mixing effects on soil structure and tillage resistance using discrete element method. Discover Soil, 3(1), Article 148. https://doi.org/10.1007/s44378-026-00308-8

Image Credits: AI Generated

DOI: 10.1007/s44378-026-00308-8

Keywords: discrete element method, tillage resistance, rotary tiller, subsoiler, soil structure, soil bin experiments, EDEM simulation, blade geometry, deep tillage, soil profile, precision agriculture, soil-tool interaction