{"type": "FeatureCollection", "features": [{"id": "10.1016/j.geoderma.2004.02.014", "type": "Feature", "geometry": null, "properties": {"updated": "2026-05-25T16:16:13Z", "type": "Journal Article", "created": "2004-04-10", "title": "Conversion Of Grassy Cerrado Into Riparian Forest And Its Impact On Soil Organic Matter Dynamics In An Oxisol From Southeast Brazil", "description": "Abstract   The purpose of this study was to evaluate possible changes in soil organic matter (SOM) dynamics after establishing riparian forests on soils previously under Brazilian savannah (\u201ccerrado\u201d). We selected a site with a homogeneous Typic Acric Red\u2013Yellow Latosol (Anionic Acrustox). Part of this site was maintained under native vegetation (grassy cerrado C 4 -dominated), and part was planted with riparian species (C 3 ) in 1992. Litter and soil samples were collected and analysed (total organic carbon, total nitrogen,  \u03b4  13 C isotopic analysis, and SOM density fractionation). Due to the predominance of grasses, carbon input was mainly below ground in cerrado. In such a soil, the decomposition process was more efficient, and much C and N were transferred to the heavy fraction. When forest was planted, there was a change from belowground to aboveground litter input (largely superficial), leading to higher C and N stocks in the light and lower stocks in the heavy fraction (resulting in lower stocks for bulk soil). The introduction of the C 3  vegetation decreased the soil  \u03b4  13 C signature. It has occurred particularly in the topsoil (0\u20135 cm) due to the deposition of C 3  litter on the soil surface. At the same time, the presence of cerrado-remaining C below 5 cm maintained higher  \u03b4  13 C values in this layer. During the 8 years after forest plantation, the input mode influenced both the  \u03b4  13 C distribution with depth, and the C replacement: between 0 and 2.5 cm, nearly 50% of cerrado-derived C was replaced by forest-derived C, while below 5 cm, replacement was around 20%. The relatively rapid C dynamics in this Oxisol (27% replacement in the top 20 cm after 8 years of forest plantation) shows that, under tropical conditions, significant changes may occur in a short period of time.", "keywords": ["delta-c-13", "decomposition", "c-13 natural-abundance", "particle-size fractions", "turnover", "0401 agriculture", " forestry", " and fisheries", "vegetation changes", "04 agricultural and veterinary sciences", "15. Life on land", "stable carbon isotope", "density fractions", "ratios", "nitrogen"], "contacts": [{"organization": "de Alcantara, F.A., Buurman, P., Furtini Neto, A.E., Curi, N., Roscoe, R.,", "roles": ["creator"]}]}, "links": [{"href": "https://doi.org/10.1016/j.geoderma.2004.02.014"}, {"rel": "related", "href": "https://repository.soilwise-he.eu/cat/collections/metadata:main/items/Geoderma", "name": "related record", "description": "related record", "type": "application/json"}, {"rel": "self", "type": "application/geo+json", "title": "10.1016/j.geoderma.2004.02.014", "name": "item", "description": "10.1016/j.geoderma.2004.02.014", "href": "https://repository.soilwise-he.eu/cat/collections/metadata:main/items/10.1016/j.geoderma.2004.02.014"}, {"rel": "collection", "type": "application/json", "title": "Collection", "name": "collection", "description": "Collection", "href": "https://repository.soilwise-he.eu/cat/collections/metadata:main"}], "time": {"date": "2004-12-01T00:00:00Z"}}, {"id": "10.1007/s10533-010-9496-4", "type": "Feature", "geometry": null, "properties": {"updated": "2026-05-25T16:14:35Z", "type": "Journal Article", "created": "2010-07-11", "title": "Effects Of Nitrogen Additions On Above- And Belowground Carbon Dynamics In Two Tropical Forests", "description": "Anthropogenic nitrogen (N) deposition is increasing rapidly in tropical regions, adding N to ecosystems that often have high background N availability. Tropical forests play an important role in the global carbon (C) cycle, yet the effects of N deposition on C cycling in these ecosystems are poorly understood. We used a field N-fertilization experiment in lower and upper elevation tropical rain forests in Puerto Rico to explore the responses of above- and belowground C pools to N addition. As expected, tree stem growth and litterfall productivity did not respond to N fertilization in either of these N-rich forests, indicating a lack of N limitation to net primary productivity (NPP). In contrast, soil C concentrations increased significantly with N fertilization in both forests, leading to larger C stocks in fertilized plots. However, different soil C pools responded to N fertilization differently. Labile (low density) soil C fractions and live fine roots declined with fertilization, while mineral-associated soil C increased in both forests. Decreased soil CO2 fluxes in fertilized plots were correlated with smaller labile soil C pools in the lower elevation forest (R2\u00a0=\u00a00.65, p\u00a0<\u00a00.05), and with lower live fine root biomass in the upper elevation forest (R2\u00a0=\u00a00.90, p\u00a0<\u00a00.05). Our results indicate that soil C storage is sensitive to N deposition in tropical forests, even where plant productivity is not N-limited. The mineral-associated soil C pool has the potential to respond relatively quickly to N additions, and can drive increases in bulk soil C stocks in tropical forests.", "keywords": ["58 Geosciences Aboveground Biomass", "15. Life on land", "Roots", "Aboveground Biomass", "Environmental sciences", "Soil Respiration", "Dissolved Organic Carbon", "Soil Density Fractions", "Environmental Chemistry", "Nutrient Limitation", "54 Environmental Sciences", "Geosciences", "Earth-Surface Processes", "Water Science and Technology"]}, "links": [{"href": "https://escholarship.org/content/qt7ww245cp/qt7ww245cp.pdf"}, {"href": "https://doi.org/10.1007/s10533-010-9496-4"}, {"rel": "related", "href": "https://repository.soilwise-he.eu/cat/collections/metadata:main/items/Biogeochemistry", "name": "related record", "description": "related record", "type": "application/json"}, {"rel": "self", "type": "application/geo+json", "title": "10.1007/s10533-010-9496-4", "name": "item", "description": "10.1007/s10533-010-9496-4", "href": "https://repository.soilwise-he.eu/cat/collections/metadata:main/items/10.1007/s10533-010-9496-4"}, {"rel": "collection", "type": "application/json", "title": "Collection", "name": "collection", "description": "Collection", "href": "https://repository.soilwise-he.eu/cat/collections/metadata:main"}], "time": {"date": "2010-07-11T00:00:00Z"}}, {"id": "10.1016/j.geoderma.2013.06.025", "type": "Feature", "geometry": null, "properties": {"updated": "2026-05-25T16:16:16Z", "type": "Journal Article", "created": "2013-07-31", "title": "Land Use And Management Effects On Soil Organic Matter Fractions In Rhodic Ferralsols And Haplic Arenosols In Bindura And Shamva Districts Of Zimbabwe", "description": "Abstract   Soil organic carbon (SOC) is a major attribute of soil quality that responds to land management activities which is also important in the regulation of global carbon (C) cycling. This study evaluated bulk soil C and nitrogen (N) contents and C and N dynamics in three soil organic matter (SOM) fractions separated by density. The study was based on three tillage systems on farmer managed experiments (conventional tillage (CT), ripping (RP), direct seeding (DS)) and adjacent natural forest (NF) in Haplic Arenosols (sandy) and Rhodic Ferralsols (clayey) of Zimbabwe. Carbon stocks were significantly larger in forests than tillage systems, being significantly lower in sandy soils (15 and 14\u00a0Mg\u00a0C\u00a0ha\u2212\u00a01) than clayey soils (23 and 21\u00a0Mg\u00a0C\u00a0ha\u2212\u00a01) at 0\u201310 and 10\u201330\u00a0cm respectively. Nitrogen content followed the same trend. At the 0\u201310\u00a0cm depth, SOC stocks increased under CT, RP and DS by 0.10, 0.24, 0.36\u00a0Mg\u00a0ha\u2212\u00a01\u00a0yr\u2212\u00a01 and 0.76, 0.54, 0.10\u00a0Mg\u00a0ha\u2212\u00a01\u00a0yr\u2212\u00a01 on sandy and clayey soils respectively over a four year period while N stocks decreased by 0.55, 0.40, 0.56\u00a0Mg\u00a0ha\u2212\u00a01 and 0.63, 0.65, 0.55\u00a0Mg\u00a0ha\u2212\u00a01 respectively. SOM fractions were dominated by mineral associated heavy fraction (MaHF) which accounted for 86\u201393% and 94\u201398% on sandy and clayey soils respectively. Tillage systems on sandy soils had the smallest average free light fraction (fLF) and occluded light fraction (oLF) C stocks (25.3\u00a0\u00b1\u00a01.3 g m\u2212\u00a02 and 7.3\u00a0\u00b1\u00a01.2\u00a0g\u00a0m\u2212\u00a02) at 0\u201330\u00a0cm when compared with corresponding NF (58.4\u00a0\u00b1\u00a04 g\u00a0m2 and 18.5\u00a0\u00b1\u00a01.0\u00a0g\u00a0m\u2212\u00a02). Clayey soils, had the opposite, having all fLF C and N in tillage systems being higher (80.9\u00a0\u00b1\u00a012\u00a0g\u00a0C m\u2212\u00a02 and 2.7\u00a0\u00b1\u00a00.4\u00a0g\u00a0N\u00a0m\u2212\u00a02) than NF (57.4\u00a0\u00b1\u00a02.0\u00a0g\u00a0C\u00a0m\u2212\u00a02 and 2.4\u00a0\u00b1\u00a00.3\u00a0g\u00a0N\u00a0m\u2212\u00a02). Results suggest that oLF and MaHF C and N are better protected under DS and RP where they are less vulnerable to mineralisation while fLF contributes more in CT. Thus, DS and RP can be important in maintaining and improving soil quality although their practicability can be hampered by unsupportive institutional frameworks. Under prevailing climatic and management conditions, improvement of residue retention could be a major factor that can distinguish the potential of different management practices for C sequestration. The exploitation of the benefits of RP or DS and the corresponding sustainability of systems need support for surface cover retention which should also be extended to conventional tillage.", "keywords": ["2. Zero hunger", "pools", "microbial biomass", "assessment", "no-tillage", "dynamics", "04 agricultural and veterinary sciences", "15. Life on land", "term changes", "carbon sequestration", "stabilization", "soil organic carbon", "conservation agriculture", "soil organic matter", "tillage", "impact", "0401 agriculture", " forestry", " and fisheries", "climate", "density fractions", "agriculture"]}, "links": [{"href": "https://doi.org/10.1016/j.geoderma.2013.06.025"}, {"rel": "related", "href": "https://repository.soilwise-he.eu/cat/collections/metadata:main/items/Geoderma", "name": "related record", "description": "related record", "type": "application/json"}, {"rel": "self", "type": "application/geo+json", "title": "10.1016/j.geoderma.2013.06.025", "name": "item", "description": "10.1016/j.geoderma.2013.06.025", "href": "https://repository.soilwise-he.eu/cat/collections/metadata:main/items/10.1016/j.geoderma.2013.06.025"}, {"rel": "collection", "type": "application/json", "title": "Collection", "name": "collection", "description": "Collection", "href": "https://repository.soilwise-he.eu/cat/collections/metadata:main"}], "time": {"date": "2013-11-01T00:00:00Z"}}, {"id": "10.5061/dryad.0k6djhb5k", "type": "Feature", "geometry": null, "properties": {"updated": "2026-05-25T16:21:22Z", "type": "Dataset", "created": "2023-08-29", "title": "Empirical data and model simulations of the effect of repeated hurricanes on soil carbon dynamics in a humid tropical forest", "description": "unspecified<em>Site description</em> Soils  were sampled from the Bisley Experimental Watershed of the LEF, Puerto  Rico (18.3157 deg. N, 65.7487 deg W), a Long-Term Ecological Research and  Critical Zone Observatory and Network site (https://luq.lter.network). The  mean maximum daily temperature at Bisley was 27 \u00baC between 1993 and 2010  (Gonzales, 2020), with little seasonality. The mean annual precipitation  at Bisley was 3883 (\u00b1 864 s.d.) mm y<sup>-1</sup> from 1988  through 2014 (Gonz\u00e1lez, 2017; Murphy et al., 2017). Rainfall occurs all  year, though January through April experience slightly less precipitation  than other months (Heartsill-Scalley et al., 2007). The site is a humid  tropical forest with a diverse tree community of approximately 170 species  &gt; 4 cm diameter at breast height (Weaver &amp; Murphy, 1990),  and dominated by tabonuco (<em>Dacryodes excelsa</em>  Vahl<em>)</em>. Elevation of Bisley spans from 261 m a.s.l. at  the base to 450 m a.s.l. on the ridges (Scatena, 1989).  Soils in Bisley are derived from volcaniclastic sediments of  andesitic parent material (Scatena, 1989).\u00a0 Ridge soils are classified as  Ultisols (Typic Haplohumults), while slope soils are classified as Oxisols  (inceptic and Aquic Hapludox), and valley soils are classified as  Inceptisols (Typic Epiaquaepts) (Hall et al., 2015; McDowell et al., 2012;  Scatena, 1989). Detailed site descriptions can be found in Scatena (1989),  Heartsill-Scalley et al (2010), and McDowell et al (2012). Here we refer  to soil organic C (SOC) and soil C interchangeably because there is no  detectable inorganic C in these soils.  <em>Hurricane occurrence\u00a0</em>  <strong>Figure 1: Timeline of major hurricanes that have  affected Luquillo Experimental Forest between sampling dates.  </strong> Nine major hurricanes (category 3 or  higher) have impacted Puerto Rico between 1851 and 2019 (L\u00f3pez-Marrero et  al., 2019), and five of these hurricanes have impacted the LEF. Until  1998, hurricanes had historically directly impacted the LEF approximately  every 60 years (Scatena &amp; Larsen, 1991). Before the initial  sampling campaign of this study, Hurricane San Cipri\u00e1n in 1932 was the  most recent storm to cause major disturbance to the LEF (Scatena &amp;  Larsen, 1991).\u00a0 However, since sampling in 1988, four major hurricanes  have impacted the forest (Figure 1). Hurricane Hugo (Category 3-4) in  1989, Hurricane Georges (Category 3) in 1998, and Hurricanes Irma and  Maria (Categories 5 and 4, respectively) within two weeks in 2017. The  trajectory and windspeeds of all these hurricanes caused widespread  defoliation. Litterfall historically takes over five years to return to  pre-hurricane levels (Scatena et al., 1996).\u00a0  <em>Sampling</em> Sample  collection occurred in 1988 and again in 2018. In both years, samples were  collected from three depths: 0\u201310 cm (the A horizon), 10\u201335 cm (all of the  B1 horizon and part of B2), and 35\u201360 cm (B2 to C) using an 8 cm diameter  soil auger. Soils in this study were sampled at three separate sites at  least 40 m from one another for each of three topographic locations,  ridge, slope, and upland valley. Two separate cores were taken from a  fourth topographic location in the riparian valley, that characterized a  smaller proportion of the area of these watersheds (Scatena &amp;  Lugo, 1995). Riparian valley sites were ephemeral streambeds with a high  boulder presence that limited sampling to less than 25 cm depth in one  case. Sampling sites from 1988 were marked with flags, and samples from  2018 were collected from within 15 m of the same locations as the  replicates from 1988, for consistency. Samples  collected in 1988 were analyzed for bulk density, pH, soil moisture, and a  suite of soil chemical properties (see Silver <em>et al</em>.  1994). Samples were then air-dried and stored in closed Ziploc bags within  paper bags in a storage facility in Richmond, CA, USA before density  fractionation in 2018. Fresh samples collected in 2018 were also  characterized for pH, soil moisture, and soil chemistry. Approximately 3 g  subsamples from each fresh sample in 2018 were immediately extracted with  45 mL of 0.2 M sodium citrate/0.5 M ascorbate solution, shaken for 16  hours, then centrifuged and the supernatant decanted to measure  concentrations of poorly crystalline iron (Fe) oxides. Within two days of  being double-bagged in Ziploc bags, fresh samples were further subsampled  and analyzed for pH in a 1:1 soil-to-water slurry (Thomas, 1996) and for  gravimetric soil moisture by oven-drying ~10 g subsamples at 105 \u00baC until  a constant weight. Soil samples were air-dried before further processing  and analysis. Air-dried soils from both sampling years were sieved to 2 mm  and large roots were sorted out. <em>Soil Density  fractionation</em> Soil was fractionated by  density following the method of Swanston et al. (2005), as modified by  Marin-Spiotta et al., (2009). Approximately 20 g of air-dried soil was  added to centrifuge tubes. Sodium polytungstate (SPT, Na6 [H2W12O40]  TC-Tungsten Compounds, Bavaria, Germany) in solution of density 1.85 g  cm<sup>-3</sup> was added to centrifuge tubes and agitated  before centrifuging. The density of the SPT followed previous studies from  this and nearby sites to allow direct comparison (Guti\u00e9rrez del Arroyo  &amp; Silver, 2018; Hall et al., 2015). Particulate organic matter  floating at the surface after centrifugation, the free light fraction  (FLF), was aspirated and then rinsed with 100 ml of deionized water 5  times on a 0.8 \u00b5m pore polycarbonate filter (Whatman Nuclepore Track Etch  Membrane, Darmstadt, Germany). Rinsed FLF was oven-dried at 65 \u00baC until  weight had stabilized. The remainder of the sample was combined with 70 ml  of additional SPT and mixed using an electric benchtop mixer (G3U05R,  Lightning, New York, NY, USA) at 1700 rpm for 1 min and sonicated in an  ice bath for 3 min at 70% pulse (Branson 450 Sonifier, Danbury, CT, USA).  Sonication is intended to disrupt soil structure and liberate organic  matter that has been occluded in aggregates. The sonicated slurry was  centrifuged again, and the light fraction at the surface, the occluded  light fraction (OLF), was aspirated, rinsed, and dried using the same  method as for the FLF. The remaining soil pellet was considered the heavy  fraction (HF), or mineral-associated organic matter fraction. The HF was  rinsed by thoroughly mixing with 150 ml of deionized water in the  centrifuge tube, centrifuging, and removing the supernatant repeatedly  until the fraction had been rinsed 5 times. The rinsed HF was oven-dried  at 105 \u00baC until weight stabilized. The average mass recovery was  98%. <em>Soil C and N and  \u03b4<sup>13</sup>C</em> Dried bulk and  HF soils were homogenized separately using a Spex Ball mill (SPEX Sample  Prep Mixer Mill 8000D, Metuchen, NJ). The FLF and OLF were homogenized  separately by hand using a mortar and pestle. All homogenized samples were  then analyzed at U. C. Berkeley for C and N concentrations on the CE  Elantech elemental analyzer (Lakewood, NJ) and for  \u03b4<sup>13</sup>C in the Stable Isotope Laboratory at UC  Berkeley, using a CHNOS Elemental Analyzer interfaced to an IsoPrime 100  mass spectrometer (Cheadle Hulme, UK), with a long-term external precision  of 0.10 %. \u00a0Soil C stocks were calculated by multiplying the C  concentrations (%) by the oven-dry mass of bulk soil (&lt; 2 mm) and  dividing by depth and the bulk density as measured in 1988 (Silver et al.,  1994; Throop et al., 2012).  <em>Radiocarbon</em> Homogenized  soil samples were combusted to CO<sub>2</sub> in sealed glass  tubes along with silver (Ag) and copper oxide (CuO) at the Center for  Accelerator Mass Spectrometry at Lawrence Livermore National Lab. The  CO<sub>2 </sub>was then graphitized on Fe powder under  pressurized hydrogen gas (Vogel et al., 1984). Graphite was pressed into  aluminum targets and run on the Compact Accelerator Mass Spectrometer for  radiocarbon analysis (Broek et al., 2021). Radiocarbon is reported in  \u0394<sup>14</sup>C, following Stuiver &amp; Polach (1977),  and calculated based on the fraction of modern isotope composition,  corrected for the year of sampling, and corrected for mass-dependent  fractionation with observed \u03b413C values of the sample. The compact AMS had  an average \u0394<sup>14</sup>C precision of 3.2 %. We report the  corrected \u0394<sup>14</sup>C value and  \u0394\u0394<sup>14</sup>C, which is calculated as  \u0394<sup>14</sup>C of the sample minus  \u0394<sup>14</sup>C of the atmosphere, to account for rapidly  changing atmospheric \u0394<sup>14</sup>C during the study period.  Atmospheric radiocarbon has been decaying nonlinearly since the peak of  weapons testing in the 1950s. Radiocarbon signatures in the soil are  strongly influenced by the atmospheric D<sup>14</sup>C  signature, making them useful for modeling soil C age and transit time,  especially since the 1950s. To compare the contribution of modern C  between 1988 and 2018, it is useful to take the difference between soil  and atmospheric D<sup>14</sup>C values, or  DD<sup>14</sup>C, because atmospheric  D<sup>14</sup>C declined between 1988 (98 %) and 2018 (4.4 %)  in Northern Hemisphere Zone 2 (Hua et al., 2013). We note that the decline  in atmospheric D<sup>14</sup>C is nonlinear, and thus the  DD<sup>14</sup>C in 2018 soil will be less sensitive to  short-term shifts in D<sup>14</sup>C inputs than the samples  from 1988. <em>Carbon age and transit time  modeling</em> Transit times and ages of C were  modeled with the package \u201cSoilR\u201d (Sierra et al., 2012, 2014) in R, version  4.0.2. The change in C density fractions over time, termed C flow, was  modeled using a 3-pool structure with a series flow matrix, under the  simplifying assumption that C flows from the litter pool to the FLF, where  it is sequentially transferred into the OLF and HF pools (Figure 2). The  model structure is depicted in basic form in equation 1,  \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0  \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (1)\u00a0 dC(t)/dt = Inputs - k*C \u00a0in  matrix form with explicit pools in equation 2,  <em>\u00a0</em> <em>\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0  \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 </em>(2)\u00a0 dC(t)/dt = [Litter Inputs; 0; 0] +  [-<em>k</em><sub>FLF</sub>, 0, 0 ;  a<sub>21</sub>,\u00a0-<em>k</em><sub>OLF</sub>, 0; 0, a<sub>32</sub>, -<sup>k</sup><sub>HRF</sub>] * [C<sub>FLF</sub>; C<sub>OLF</sub>; C<sub>HF</sub>] where <em>k</em><strong> </strong>is the first-order decay constant for each pool, <em>a</em> is the C transfer rate between pools (<em>i.e. a<sub>21</sub> </em>is the transfer from FLF (pool 1) to OLF (pool 2) and <em>a<sub>32</sub></em> is the transfer from OLF (pool 2) to HF (pool 3)), and <em>C </em>is the C stock of each pool.<strong> </strong>The transitTime and systemAge functions within the \u201csoilR\u201d package use this model structure to solve for the distribution of ages (time since entry) of each pool, and the distribution of transit times (times between entry and exit from the bulk soil) (Sierra et al 2016). Distributions of age and transit time were time-independent and did not assume a specific distribution (Sierra et al., 2014, 2017). <strong>Figure 2: Hypothesized flow of C in soils. </strong> Free light fraction (FLF) C (pink) is either decomposed (at cycling rate -<em>k<sub>FLF </sub>* FLF</em>) or transferred to the occluded light fraction pool (OLF, blue) with the transfer proportion defined by <em>a<sub>21</sub></em>. Carbon transfer between the OLF and heavy fraction (HF, purple) is defined by transfer coefficient <em>a<sub>32</sub></em>, and is respired from these pools at cycling rates -<em>k<sub>OLF</sub>* OLF</em> and <em>-k<sub>HF</sub>* HF</em>, respectively. Figure adapted from Sierra et al. (2012). Soil D<sup>14</sup>C and C stock mean and standard deviations from each time point, depth, and fraction were used to constrain the matrix model describing the movement of C through three soil pools and losses of C from each pool. Topography was not a strong predictor of patterns in D<sup>14</sup>C, C stocks, or C fractions, so samples from all topographies were aggregated for model simulations. The model used mean observed C content in each pool for each depth in 1988 as initial conditions for SOC stocks. Above and belowground litter inputs at 0\u201310 cm were assumed to be 900 g C m<sup>-2</sup> in non-hurricane or hurricane recovery years, based on observations from the same site (Liu et al., 2018; Scatena et al., 1996; Silver et al., 1996; Vogt et al., 1996). Inputs to the 10\u201335 cm and 35\u201360 cm depths were estimated using observations of live fine roots on the surface and typical root distribution in the forest (Silver &amp; Vogt, 1993). Total root input is approximately threefold the input of fine roots alone (McCormack et al., 2015; Yaffar &amp; Norby, 2020), and live fine roots in the 0\u201310 cm depth had a mean biomass of 80 - 250 g C m<sup>-2 \u00a0</sup>(Hall et al., 2015), suggesting that total root C inputs of approximately 450 g C m<sup>-2 </sup>to the surface would be well within the expected range. Root inputs below 0\u201310 cm were estimated assuming that inputs follow the typical distribution of root biomass in Puerto Rican tropical forests, with 60\u201370% of root biomass in 0\u201310 cm, an additional 20-30% of biomass in 10\u201335 cm (~135 g C m<sup>-2\u00ad</sup>), and 5\u20138% of biomass is in the 35\u201360 cm depth (~40 g C m<sup>-2\u00ad</sup>) (Silver &amp; Vogt, 1993; Yaffar &amp; Norby, 2020). The model was parameterized under two scenarios for each depth: 1) constant inputs, assuming a steady-state undisturbed forest, and 2) hurricane inputs, which simulated the input fluxes from defoliation during the three major hurricanes, followed by a subsequent reduction in litter inputs and then litterfall increasing linearly to pre-hurricane inputs over 6 years (Scatena et al., 1996; Silver et al., 1996; Vogt et al., 1996). Hurricane inputs were imposed as an additional pulse of litter inputs to each depth interval, declining with depth. \u00a0The 0\u201310 cm interval received 100% of the surface input pulse, the 10\u201335 cm depth received a pulse of root inputs equivalent to 30% of the surface input pulse, and the 35\u201360 cm depth received root inputs equal to 10% of the surface input pulse. Surface litter pulses under hurricanes were specified according to measured litterfall values and were 42.5 g C m<sup>-2\u00ad</sup> to the surface in 1989 (Hurricane Hugo) and 1998 (Hurricane Georges) (Scatena et al., 1993; Silver et al., 1996) and 1611 g C m<sup>-2 \u00a0</sup>in 2017 (Hurricanes Irma and Maria) (Liu et al. 2018a). The same soil D<sup>14</sup>C and C stock observations were used to constrain the model under each scenario, with only the input regime varying. Parameters of the transfer matrix (<em>-k\u00ad\u00ad<sub>FLF</sub>,</em><sub> </sub><em>-k\u00ad\u00ad<sub>OLF</sub>,<sub> </sub>-k\u00ad\u00ad<sub>HF</sub>,<sub> </sub>a<sub>21</sub>, a<sub>32</sub></em>) were constrained using a cost function to accept or reject potential parameter sets over 1000 iterations, based on observed D<sup>14</sup>C and C stock means and standard errors from both time points (1988 and 2018). A Markov chain Monte Carlo (MCMC) simulation initialized with cost-optimized parameters was run to assimilate observed data and optimize parameter choices to the observations using function <em>modMCMC() </em>from R package \u201cFME\u201d (Sierra et al., 2014; Soetaert &amp; Petzoldt, 2010). The MCMC was iterated over at least 20,000 simulations or until parameter solutions converged according to the trace, which was over 100,000 iterations at the 35\u201360 cm depth. The first half of the iterations was considered the burn-in period before the chain started to converge near an equilibrium, and these iterations were discarded in calculations of optimal parameters. The model output for the surface soils of the HF pool was validated using published radiocarbon values from the mineral-associated fraction (the only fraction analyzed) of samples from the site taken in 2012 (Hall et al., 2015).\u00a0 Bulk and pool soil C age and transit time density distributions and mean values were calculated using the <em>systemAge() </em>and <em>transitTime()</em> functions from the \u201cSoilR\u201d package. Mean density distributions were calculated using the mean parameter set given from the MCMC analysis. Standard deviation from the mean was calculated using the <em>systemAge() </em>and <em>transitTime()</em> functions on 200 sets of five parameters selected randomly within one standard deviation of the mean of each parameter given as output from the MCMC. Lower and upper limits of SOC ages and transit times were calculated using the upper and lower ranges of these iterations. <em>Statistics</em> Statistics were run in R, version 4.0.2 (R Core Team, 2020). The statistical model selection followed the recommendations of Zuur et al (2009). Statistical models were chosen using a linear mixed effects model in package \u201clme4\u201d, with random slopes accounting for the influence each core, or sampling site, had on the response variable values as they varied with depth. This random effect of the core site on the depth effect was evaluated using a restricted maximum likelihood approach and was included in the initial evaluation of all model comparisons. Linear mixed effect models included year, topographic position, depth, and interactions as fixed factors, and the depth effect of each core as a random factor for each of the response variables: C concentration, N concentration, d<sup>13</sup>C, DD<sup>14</sup>C. In evaluations of some response variables with AIC and BIC criteria, the random effect no longer enhanced the model, and model comparison proceeded using ANOVAs of linear models without random effects. Topographic effects on C concentrations are discussed in the supplemental information. Model assumptions were evaluated using the check_model function in R package \u201cperformance\u201d, to check for multicollinearity, normality of residuals, homoscedasticity, homogeneity of variance, influential observations, and normality of random effects. In the cases when random effects were significant (bulk soil d<sup>13</sup>C and DD<sup>14</sup>C, FLF DD<sup>14</sup>C and HF C and N concentrations), fixed effects were chosen using ANOVA of subsequent models using maximum likelihood estimation, with the random effects held constant. Once fixed effects were established, the model was re-fitted using a restricted maximum likelihood approach to report model estimates, and an ANOVA was run to determine the significance of the response variable. In all cases, P-values were estimated using Tukey\u2019s honest significant post-hoc test to assess significant differences between variables, in the package \u201cagricolae\u201d in R, and contrasts and standard errors of contrasts were estimated using lsmeans() function in package \u201clsmeans\u201d in R. Values of\u00a0<em>P</em> &lt; 0.10 were reported as significant unless otherwise specified. The topographic position was not a significant predictor for most variables, so results are reported as means aggregated across positions.", "keywords": ["soil organic carbon", "Transit time", "Tropical forest soil", "FOS: Earth and related environmental sciences", "Soil R", "density fractions", "Radiocarbon"], "contacts": [{"organization": "Mayer, Allegra", "roles": ["creator"]}]}, "links": [{"href": "https://doi.org/10.5061/dryad.0k6djhb5k"}, {"rel": "self", "type": "application/geo+json", "title": "10.5061/dryad.0k6djhb5k", "name": "item", "description": "10.5061/dryad.0k6djhb5k", "href": "https://repository.soilwise-he.eu/cat/collections/metadata:main/items/10.5061/dryad.0k6djhb5k"}, {"rel": "collection", "type": "application/json", "title": "Collection", "name": "collection", "description": "Collection", "href": "https://repository.soilwise-he.eu/cat/collections/metadata:main"}], "time": {"date": "2024-04-01T00:00:00Z"}}, {"id": "10.5194/bg-2-159-2005", "type": "Feature", "geometry": null, "properties": {"updated": "2026-05-25T16:21:38Z", "type": "Journal Article", "created": "2010-04-29", "description": "<p>Abstract. Extreme sensitivity of soil organic carbon (SOC) to climate and land use change warrants further research in different terrestrial ecosystems. The aim of this study was to investigate the link between aggregate and SOC dynamics in a chronosequence of three different land uses of a south Chilean Andisol: a second growth Nothofagus obliqua forest (SGFOR), a grassland (GRASS) and a Pinus radiata plantation (PINUS). Total carbon content of the 0-10cm soil layer was higher for GRASS (6.7 kg C m-2) than for PINUS (4.3 kg C m-2, while TC content of SGFOR (5.8 kg C m-2) was not significantly different from either one. High extractable oxalate and pyrophosphate Al concentrations (varying from 20.3-24.4 g kg-1, and 3.9-11.1 g kg-1, respectively) were found in all sites. In this study, SOC and aggregate dynamics were studied using size and density fractionation experiments of the SOC, \uffce\uffb413C and total carbon analysis of the different SOC fractions, and C mineralization experiments. The results showed that electrostatic sorption between and among amorphous Al components and clay minerals is mainly responsible for the formation of metal-humus-clay complexes and the stabilization of soil aggregates. The process of ligand exchange between SOC and Al would be of minor importance resulting in the absence of aggregate hierarchy in this soil type. Whole soil C mineralization rate constants were highest for SGFOR and PINUS, followed by GRASS (respectively 0.495, 0.266 and 0.196 g CO2-Cm-2d-1 for the top soil layer). In contrast, incubation experiments of isolated macro organic matter fractions gave opposite results, showing that the recalcitrance of the SOC decreased in another order: PINUS&gt;SGFOR&gt;GRASS. We deduced that electrostatic sorption processes and physical protection of SOC in soil aggregates were the main processes determining SOC stabilization. As a result, high aggregate carbon concentrations, varying from 148 till 48 g kg-1, were encountered for all land use sites. Al availability and electrostatic charges are dependent on pH, resulting in an important influence of soil pH on aggregate stability. Recalcitrance of the SOC did not appear to largely affect SOC stabilization. Statistical correlations between extractable amorphous Al contents, aggregate stability and C mineralization rate constants were encountered, supporting this hypothesis. Land use changes affected SOC dynamics and aggregate stability by modifying soil pH (and thus electrostatic charges and available Al content), root SOC input and management practices (such as ploughing and accompanying drying of the soil).                     </p>", "keywords": ["DECOMPOSITION", "NEW-ZEALAND", "DENSITY FRACTIONS", "[SDU.ASTR] Sciences of the Universe [physics]/Astrophysics [astro-ph]", "HUMIC-ACID", "Life", "QH501-531", "QH540-549.5", "2. Zero hunger", "QE1-996.5", "CULTIVATED SOILS", "Ecology", "[SDU.OCEAN] Sciences of the Universe [physics]/Ocean", " Atmosphere", "Geology", "LAND-USE CHANGE", "04 agricultural and veterinary sciences", "ALUMINUM", "15. Life on land", "[SDU.ENVI] Sciences of the Universe [physics]/Continental interfaces", " environment", "MACROORGANIC MATTER", "C SEQUESTRATION", "[PHYS.ASTR.CO] Physics [physics]/Astrophysics [astro-ph]/Cosmology and Extra-Galactic Astrophysics [astro-ph.CO]", "Earth and Environmental Sciences", "FOREST SOILS", "[SDU.STU] Sciences of the Universe [physics]/Earth Sciences", "0401 agriculture", " forestry", " and fisheries"], "contacts": [{"organization": "Huygens, D., Boeckx, P., van Cleemput, O., Oyarz\u00fan, C., Godoy, R.,", "roles": ["creator"]}]}, "links": [{"href": "https://doi.org/10.5194/bg-2-159-2005"}, {"rel": "related", "href": "https://repository.soilwise-he.eu/cat/collections/metadata:main/items/Biogeosciences", "name": "related record", "description": "related record", "type": "application/json"}, {"rel": "self", "type": "application/geo+json", "title": "10.5194/bg-2-159-2005", "name": "item", "description": "10.5194/bg-2-159-2005", "href": "https://repository.soilwise-he.eu/cat/collections/metadata:main/items/10.5194/bg-2-159-2005"}, {"rel": "collection", "type": "application/json", "title": "Collection", "name": "collection", "description": "Collection", "href": "https://repository.soilwise-he.eu/cat/collections/metadata:main"}], "time": {"date": "2005-06-24T00:00:00Z"}}], "links": [{"rel": "self", "type": "application/geo+json", "title": "This document as GeoJSON", "href": "https://repository.soilwise-he.eu/cat/collections/metadata:main/items?keywords=density+fractions&f=json", "hreflang": "en-US"}, {"rel": "alternate", "type": "text/html", "title": "This document as HTML", "href": "https://repository.soilwise-he.eu/cat/collections/metadata:main/items?keywords=density+fractions&f=html", "hreflang": "en-US"}, {"rel": "collection", "type": "application/json", "title": "Collection URL", "href": "https://repository.soilwise-he.eu/cat/collections/metadata:main", "hreflang": "en-US"}, {"type": "application/geo+json", "rel": "first", "title": "items (first)", "href": "https://repository.soilwise-he.eu/cat/collections/metadata:main/items?keywords=density+fractions&", "hreflang": "en-US"}, {"rel": "last", "type": "application/geo+json", "title": "items (last)", "href": "https://repository.soilwise-he.eu/cat/collections/metadata:main/items?keywords=density+fractions&offset=5", "hreflang": "en-US"}], "numberMatched": 5, "numberReturned": 5, "distributedFeatures": [], "timeStamp": "2026-05-26T00:11:40.036699Z"}