SELFLIMITING ECOSYSTEM ENGINEERING IN NONKARST REGOLITH
This is a followup to my previous post on emergent ecosystem engineering in epikarst, so I won't repeat much of the background or analytical details. There I argued that interactions among rock weathering, moisture flux, biological effects (particularly roots and their symbionts) and soil operate such that if weathering is moisturelimited, and biota are limited by water availability and belowground space, the system is dynamically unstable. Positive feedbacks dominate so as to reinforce or accelerate dissolution, joint/fracture widening, root growth, and soil accumulation. The net effect is to develop the epikarst as increasingly hospitable habitat. This continues, according to the analysis, until weathering becomes reactionlimited and subsurface space and moisture are no longer significant limiting factors for plant growth. Under the latter circumstances the system is dynamically stable, implying resilience to relatively small changes or disturbances and slower change.
But what about noncarbonate parent rocks? There are some key differences. In epikarst, ground water flow is dominated by fissures and conduits. In nonkarst, matrix and macropore flow may be more important, and chemical weathering may facilitate moisture penetration independently of widening of joints and fractures. Also, though this is a bit oversimplified, dissolution of carbonates results in removal of material in solution and the creation or widening of fissures and cavities. In other rocks, chemical weathering is often more selective, removing or modifying more soluble and less resistant minerals first, and increasingly concentrating more resistant ones. Thus the rate of weathering slows down as weatherable minerals are depleteda selflimitation effect independent of moisture penetration. Also, because conduits are less common (and smaller), plugging or flushing of sediments therein is rarely an issue.
Granitic saprolite from upland South Carolina with fine roots in a weathered parting >2 m below the surface.
Figure 1 below is a similar system model to that of the epikarst post, to which you can refer for a more complete explanation of the links. To cut to the chase, applying the same type of analysis to the nonepikarst model produces the same general result.
Figure 1. Interaction system. Dotted lines indicate links that are positive when present, but that may become inapplicable as regolith develops.
When selflimitations (depletion of weatherable minerals, biological saturation, external moisture supply) are not approached, biological production is stimulated by additional space and subsurface moisture, and chemical weathering is moisturelimited, the system is dynamically unstable. Positive feedbacks dominate, and absent disturbances that radically change any component, weathering, soil, biota, and moisture flux are all mutually reinforcing. As this creates ever more suitable habitat, this can be viewed as ecosystem engineering (EE). Dynamical stability obtains when selflimitations are approached, biota are not generally limited by water or belowground space, and chemical weathering is reaction rather than moisture limited most of the time.
As with the epikarst case, though the processes by which EE occurs may continue, they no longer operate so as to create more favorable habitat for the surface organisms (woody plants) primarily responsible. Thus in both cases, the EE is emergent in that it derives from interactions that occur in certain exposed rock or thinsoil settings and is not directly associated with traits of the engineer organisms. And in both epikarst and noncarbonate settings, the EE is selflimiting as factors other than water and subsurface space become constraining.
Table 1. Interaction matrix for the system shown in Figure 1.

Joints, fractures^{1} 
Weathering 
Biota 
Water flux 
Soil 
JF 
0 
0 
a_{13} 
a_{14} 
0 
W 
a_{21} 
a_{22} 
0 
a_{24} 
a_{25} 
B 
0 
a_{32} 
a_{33} 
a_{34} 
a_{35} 
WF 
0 
a_{42} 
a_{43} 
a_{44} 
+a_{45} 
S 
0 
a_{52} 
a_{53} 
0 
a_{55} 
Table 2. Feedbacks for the system shown in Figure 1 and Table 1. One of the criteria for dynamical stability is that all FN < 0.
FN = feedback at level N 
F1 = (a_{22}) + (a_{33}) + (a_{44}) + (a_{55}) 
F2 = a_{24}a_{42 }+ a_{43}a_{34} + a_{52} a_{25} + a_{53}a_{35} +  (a_{33})(a_{44})  (a_{33})(a_{55})  (a_{44})(a_{55}) . . . 
F3 = a_{13}a_{32}a_{21 }+ a_{14}a_{42}a_{21 }+ a_{24}a_{43}a_{32 }+ a_{24}(+a_{45})a_{52} a_{25}a_{53}a_{32 }+ a_{35}(a_{54})a_{42} +a_{45}a_{53}a_{34 } a_{43}a_{34}(a_{22})  a_{53}a_{35}(a_{22}) + (a_{33})(a_{44})(a_{55}) 
F4 = a_{21}a_{14}a_{43}a_{32} + a_{24}a_{43}a_{35}a_{52}_{ +} a_{52}a_{21}a_{14}(+a_{45}) + a_{35}a_{52}a_{21}a_{13} + (a_{54})a_{43}a_{32}a_{25}  a_{13}a_{32}a_{21}(a_{44})  a_{25}a_{53}a_{32}(a_{44})  a_{13}a_{32}a_{21}(a_{55})  a_{14}a_{42}a_{21}(a_{33})  a_{14}a_{42}a_{21}(a_{55})  (a_{43} a_{34})(a_{52} a_{25})  a_{13}a_{32}a_{21}(a_{44})  a_{13}a_{32}a_{21}(a_{55}) 
F5 = a_{52}a_{21}a_{13}a_{34}(+a_{45}) + a_{53}a_{32}a_{21}a_{14}(+a_{45})  a_{21}a_{14}a_{43}a_{32}(a_{55})  + a_{52}a_{21}a_{14}(+a_{45})(a_{33})  a_{35}a_{52}a_{21}a_{13}(a_{44})  (a_{14}a_{42}a_{21})(a_{53}a_{35}) + (a_{13}a_{32}a_{21})(a_{44})(a_{55}) + a_{14}a_{42}a_{21}(a_{33})(a_{55}) 
(Posted 23 September 2017)