Subscriber access provided by PEPPERDINE UNIV
Article
Identification of distinct functional microstructural domains controlling C storage in soil Markus Steffens, Derek M Rogge, Carsten W. Mueller, Carmen Höschen, Johann Lugmeier, Angelika Kölbl, and Ingrid Kogel-Knabner Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b03715 • Publication Date (Web): 28 Sep 2017 Downloaded from http://pubs.acs.org on September 28, 2017
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 20
Environmental Science & Technology
TOC/graphical abstract
ACS Paragon Plus Environment
Environmental Science & Technology
1
Identification of distinct functional microstructural domains controlling C storage in soil
2 3
Markus Steffens1,2*, Derek Rogge3, Carsten W. Mueller1, Carmen Höschen1, Johann Lugmeier1, Angelika Kölbl1, Ingrid Kögel-Knabner1,4
4 1
5 6
Lehrstuhl für Bodenkunde, Wissenschaftszentrum Weihenstephan, Technische Universität München, Emil-Ramann-Str. 2, D-85354 Freising, Germany
7 8
2
9 10
3
11 12
4
13
*
[email protected] Department of Soil Sciences, Research Institute of Organic Agriculture (FiBL), Ackerstrasse 113, CH-5070 Frick, Switzerland Applied spectroscopy group, Deutsche Forschungsanstalt für Luft- und Raumfahrt Oberpfaffenhofen, D-82234 Wessling, Germany Institute of Advanced Study, Technische Universität München, Lichtenbergstraße 2a, D85748 Garching, Germany
14 15
The physical, chemical, and biological processes forming the backbone of important soil
16
functions (e.g. carbon sequestration, nutrient and contaminant storage, and water transport)
17
take place at reactive interfaces of soil particles and pores. The accessibility of these
18
interfaces is determined by the spatial arrangement of the solid mineral and organic soil
19
components, and the resulting pore system. Despite the development and application of novel
20
imaging techniques operating at the micro- and even nanometer scale, the microstructure of
21
soils is still considered as a random arrangement of mineral and organic components. Using
22
nanoscale secondary ion mass spectroscopy (NanoSIMS) and a novel digital image processing
23
routine adapted from remote sensing (consisting of image preprocessing, endmember
24
extraction and a supervised classification), we extensively analyzed the spatial distribution of
25
secondary ions characteristic for mineral and organic soil components on the sub-micrometer-
26
scale in an intact soil aggregate (40 measurements, each covering 30×30 µm2 with a lateral
27
resolution of 100×100 nm). We were surprised that the 40 spatially independent
28
measurements clustered in just two, complimentary types of micrometer-sized domains. Each
29
domain is characterized by a micro-architecture built of a definite mineral assemblage with
ACS Paragon Plus Environment
Page 2 of 20
Page 3 of 20
Environmental Science & Technology
30
various organic matter forms and a specific pore system, each fulfilling different functions in
31
soil. Our results demonstrate that these micro-architectures form due to self-organization of
32
the manifold mineral and organic soil components to distinct mineral assemblages, which are
33
in turn stabilized by biophysical feedback mechanisms acting through pore characteristics and
34
microbial accessibility. These microdomains are the smallest units in soil that fulfill specific
35
functionalities.
ACS Paragon Plus Environment
Environmental Science & Technology
36
Introduction
37
Soil functions are an essential component of ecosystem services, and cover complex processes
38
as provision of water and nutrients for plant production, purification of water and air, and
39
regulation of the global climate via carbon sequestration.1, 2 These functions emerge from the
40
soil architecture - the three-dimensional spatial arrangement of a soil’s manifold mineral and
41
organic components, and pores. This architecture determines the reaction space for all basic
42
physical, chemical, and biological processes3 that are the backbone of soil functions. Up to
43
now soil functions are deduced from studies investigating either properties of isolated mineral
44
and organic soil components4, or the pore system with aggregates as a key unit of soil
45
structure5. However, to mechanistically understand how soil functions emerge from these
46
components and their spatial arrangement, it is necessary to study the interplay between solid
47
components and the pore system within intact soil samples5-10. Prevailing concepts postulate
48
that different soil minerals like phyllosilicates and oxy-hydroxides form definite spatial
49
arrangements with distinct organic components on the micrometer scale11-13. Until today,
50
direct evidence for these concepts is missing. Contrastingly, current spectromicroscopic
51
investigations describe soil microstructures as heterogeneous mixtures of mineral and organic
52
components in random arrangement14-16. In order to elucidate the spatial heterogeneity of
53
microstructures in an intact soil macroaggregate, we extensively investigated the spatial
54
arrangement of mineral and organic components as well as pores within a millimeter-sized
55
intact soil aggregate using nanoscale secondary ion mass spectroscopy (NanoSIMS; Figure 1).
56
To account for the spatial heterogeneity of the soil aggregate, we recorded 40 images (each
57
30 µm × 30 µm) in a regular grid on an intact soil aggregate. We combined digital image
58
processing, supervised classification and cluster analyses in a new image data analyses routine
59
in order to synoptically analyze the multivariate information in this large, spatially resolved
60
data set. Using this routine, we were able to group the various spatial arrangements of the 40
61
images in three distinct, recurring types of micro-architectures. These micro-architectures ACS Paragon Plus Environment
Page 4 of 20
Page 5 of 20
Environmental Science & Technology
62
showed clearly different soil properties like local pore size distribution, microbial
63
accessibility, and sorption potential and could be identified as either storage or transport
64
domain. Here we show that the mineral composition of the mineral grains in combination with
65
the organic matter are the nucleus for the formation of distinct micro-architectures that allow
66
specific functionalities in soils - functional microdomains.
67
Materials and Methods
68
The sampling site is located at the Research Station Scheyern of the Technische Universität
69
München, 40 km north of München in southern Germany (445-498 m above sea level), with a
70
mean annual temperature of 7.4 °C and a mean annual precipitation of 803 mm. The intact
71
soil macroaggregate was taken from the Ap-horizon (0-20 cm) of a conventionally managed,
72
loamy Cambisol (27% sand, 52% silt, 21% clay), derived from coarse and fine-grained
73
tertiary sediments covered by Quaternary Loess17 - a soil material typical and representative
74
for many agriculturally managed areas in Europe. Mineralogy of the bulk material and the
75
clay-sized fraction was assessed by X-ray diffraction (XRD). Prior to XRD analyses, the
76
samples were treated with hydrogen peroxide to remove organic matter (H2O2). XRD was
77
performed on random powder samples and on oriented samples (Co-Kα; Philips
78
diffractometer PW 1830) after saturation with Ca2+ and glycerol (room temperature) and K+
79
(room temperature and stepwise heated to 560°). The bulk sample consists of 50% quartz,
80
30% micas, and 20% feldspars, while the clay-sized fraction contained 70% illite, 20%
81
smectite, and 10% kaolinite.
82
The millimeter-sized, air-dried intact soil aggregate was embedded using epoxy resin
83
(Araldite kit 502, electron microscope sciences, Hatfield, USA). For larger specimen as the
84
analyzed macroaggregate (>10 mm in diameter) consisting of a prominent amount of quartz
85
grains, the chemical embedding and sectioning approach is a secure approach to obtain intact
86
microscale structures with no distortion due to drying18. We used an epoxy resin that ACS Paragon Plus Environment
Environmental Science & Technology
Page 6 of 20
87
demonstrates a unique matrix effect with clearly distinguishable differences between resin and
88
organic matter derived signals in the NanoSIMS19. The embedded sample was cured for 48
89
hours at 60°C, cut, and subsequently polished in order to obtain a sample surface with low
90
topography20. Prior to NanoSIMS analysis, the sample was investigated using a reflectance
91
light microscope (Axio Imager Z2m, Zeiss, Germany). We designed a regular, staggered
92
sampling grid with 400 µm sampling distance yielding 40 locations for NanoSIMS imaging
93
(Figure 1). In total 1.2% of the polished sample surface were analyzed with NanoSIMS.
94
NanoSIMS measurements were carried out with a NanoSIMS 50L instrument (Cameca,
95
Gennevilliers, France) at the Lehrstuhl für Bodenkunde, Technische Universität München,
96
Freising, Germany18. To compensate for charging on the non-conductive mineral soil
97
particles, we sputter-coated (SCD 005 sputter coater, Bal-Tec GmbH, Germany) our sample
98
with a conductive layer of Au/Pd (∼30 nm) prior to analysis and used the electron flood gun
99
during the NanoSIMS measurements. The Cs+ primary beam was used at an approximate
100
lateral resolution of 100 nm. The primary beam current was ∼1-2 pA, with 2 ms pix-1 dwell
101
time, and 5 planes were recorded. The size of each image was 30 µm × 30 µm at 512 × 512
102
pixels. In preparation, an area larger than the spot to be measured was pre-sputtered with a
103
high primary beam current (∼500 pA). This eliminated the Au-coating layer, cleaned the
104
sample of air adsorbates, and promoted ionization of the elements by Cs implantation. When
105
signal equilibrium was reached, the NanoSIMS measurement was initiated. Electron
106
multiplier secondary ion collectors were used for
107
56
108
organic components. Nitrogen does not ionize, and was detected as
109
secondary ions
110
components. The secondary ions are retrieved from the sample surface with a depth resolution
111
of approximately 10 nm. We calibrated the machine using a Si grid sample and regularly
Fe16O-. The
12
C- and
28
Si-,
12
C-,
16
O-,
12
C14N-,
28
Si-,
27
Al16O- and
12 14
C N- secondary ions were recorded to monitor the distribution of
27
Al16O-, and
56
12 14
C N- cluster ion. The
Fe16O- were recorded as representative of the mineral
ACS Paragon Plus Environment
Page 7 of 20
Environmental Science & Technology
112
checked for optimal beam centering using the Cameca integrated tuning software. We
113
suppose that ratios between secondary ions remain constant for one matrix if measured at
114
different locations or dates, even if the absolute counts can vary.
115
To allow a direct comparison of a large number of NanoSIMS images measured at several
116
days and different sample locations, we developed a new routine in IDL/ENVI 5.2 (Exelis
117
Visual Information Solutions, Boulder, Colorado, USA) for NanoSIMS image analyses. The
118
key advantage of this routine for NanoSIMS image analysis is the ability to evaluate the
119
whole data set simultaneously across all measured secondary ions at all measured locations.
120
This routine makes it possible to sort every pixel into a class with a similar or comparable
121
“mass signature” and map these classes across all images. It is an important prerequisite for
122
the synoptical analyses of spectromicroscopic data in geoscience. In order to analyse the
123
whole data set simultaneously we first need to concatenate the 40 images, each containing the
124
images of six secondary ions, to produce an L×S×E image cube (L = number of lines; S =
125
number of samples; E = number of ion species) where each pixel of the contiguous image
126
represents a vector, or mass signature, comprising all measured secondary ions. This image
127
cube being analogous to a reflectance image in imaging spectroscopy comprising spectral
128
measurements at specific wavelengths, either contiguous or not, resulting in each pixel having
129
a reflectance profile. With such an image cube array it is now possible to analyse all pixels as
130
vectors comprising all measured secondary ions.
131
The developed routine consists of three main steps: Step 1 is the data preparation including
132
filtering to reduce random noise and intra-class variability, and normalization to allow for
133
transferability across different measurement dates or locations. Data preparation is followed
134
by an image-based endmember extraction (step 2) and a supervised classification of the whole
135
image using a spectral angle approach on the extracted endmembers (step 3).
ACS Paragon Plus Environment
Environmental Science & Technology
136
The high sensitivity (ppm to ppb range) of the NanoSIMS can lead to an imaging “noise” in
137
the recorded measurements that can hamper image segmentation. Reduction of this noise can
138
be applied per ion-species image using simple approaches such as median or mean filtering.
139
In this study we make use of the iterative adaptive spatial-spectral filter (IAS filter) that is
140
designed to reduce random noise and intra-class variability, while at the same time being
141
effective at retaining edges and boundaries across materials (step 1). The advantage of this
142
particular filter is that it smooths the image data by considering the local spatial variance of
143
each pixel compared to its surroundings using the vector based standard spectral similarity
144
measures spectral angle21 and root-mean-squared error. The IAS filter is a simplified variant
145
of a vector based anisotropic diffusion approach22 and is characterized by minimal
146
algorithmic complexity and computational load.
147
Once filtering is complete the next data preparation step is to apply a per-pixel spectral mean
148
normalization23 to normalize varying count rates between the six secondary ion
149
measurements. The result of this normalization is to retain the ratios between secondary ions,
150
while removing overall amplitude. The approach allows for transferability across the different
151
measurement locations and is vital for the synoptical analysis of the six secondary ion images.
152
A mean normalized secondary ion-profile is calculated per-pixel with the following equation:
153
, =
154
where c is the ion species identifier, x,y is the pixel vector location in the image, Nx,y(c) is the
155
normalized secondary ion profile for a specific pixel, Rx,y(c) is the original secondary ion
156
profile, and Mx,y is the mean value across all secondary ions for the pixel. By applying this
157
normalization, the ratio of all secondary ions to each other is kept constant for one pixel but
158
the absolute count numbers are normalized and made comparable between different pixels
159
and therefore between matrices and measurement locations. In this step we suppose that the
, ,
(1)
ACS Paragon Plus Environment
Page 8 of 20
Page 9 of 20
Environmental Science & Technology
160
absolute counts vary between sampling locations, matrices and sampling dates but the ratios
161
between the measured secondary ions remain constant (Figure S5, S6, and S7).
162
Endmember extraction (step 2) is an approach commonly used in remote sensing to extract
163
pure spectra of unique materials24. We use the approach to identify pixels, which are
164
characterised by a unique mass signature, so called endmembers from the full L×S×E image
165
cube and use their signature for further steps. These endmembers define the extreme points
166
(vertices) of the six-dimensional data space (simplex) spanned by the six analysed secondary
167
ions, and all other pixels are considered mixtures of these few identified endmembers. In this
168
study we use the Spatial Spectral Endmember Extraction (SSEE) method25 which is
169
commonly used and established in the analysis of hyperspectral remote sensing imagery. The
170
advantage of this specific approach over many other endmember extraction methods is its
171
ability to use both spectral and spatial information to define unique signatures, resulting in a
172
larger set of potential endmembers that may include spectrally similar, but spatially
173
independent materials. This approach allows for a good representation of the intra-class
174
variability within a given class type and allows for expert knowledge to be integrated in the
175
selection process to select a final endmember set appropriate for a specific application. We
176
identified 30 endmembers with unique mass signatures in the 40 NanoSIMS images (Figure
177
S7). We analysed the likeliness of the mass signatures of the 30 endmembers using a
178
hierarchical cluster analysis based on the squared Euclidean distance (SPSS 23; IBM SPSS
179
Statistics for Windows, Armonk, NY). This analysis in combination with the reference spectra
180
from test data sets 1 and 2 showed that the 30 endmembers could be clustered into 10
181
compound groups with comparable mass signatures (Figure S9, table S3).
182
For the classification of the 11,010,048 pixels in the NanoSIMS images into these 30 mass
183
signature classes we used the spectral angle mapper (SAM) algorithm26 to assign each pixel to
184
one of the unique endmember classes (step 3; Figure S8). SAM is a nonparametric classifier
ACS Paragon Plus Environment
Environmental Science & Technology
185
that determines the signature similarity between each pixel’s signature with the set of
186
endmembers. Each pixel is assigned to the endmember class with the highest similarity based
187
on the smallest angle to the reference signature. In this study, we use the ion profile
188
endmembers derived from SSEE as our reference to classify the image cube. The angle (ө)
189
between two ion-signatures si = si1, …, siC and sj = sj1, …, sjC, where c is the number of bands
190
(or ion-species), is calculated as:
191
, = ө , =
∑!" # "⁄% "⁄% % %) $∑!" & (∑!" #
* (2)
192
193
Based on the contribution of these 30 mass signature classes to each image, we clustered the
194
40 images in three types of micro-architectures using a second cluster analyses based on their
195
squared Euclidean distance (Figure S10, table S4).
196
We used two test data sets with increasing spatial complexity to evaluate the power of our
197
NanoSIMS image analyses routine. Test data set 1: We applied the technique on seven pure
198
standard materials representing pedogenic aluminium and iron oxides, and phyllosilicates in
199
order to evaluate the power of our routine to identify homogeneous and pure minerals. These
200
materials were measured one by one so that no mixing of materials could hamper their
201
identification. The iron oxide goethite was synthesized following the procedure of
202
Schwertmann & Cornell27 in three varieties with increasing Al substitution (Figure S1 a-c).
203
Additionally we purchased four commercial minerals – Fe oxide Goethite (‘Bayferrox 920’;
204
Lanxess Deutschland GmbH, Leverkusen, Germany), Al oxide Boehmite (‘APYRAL AOH
205
20’; Nabaltec AG, Schwandorf, Germany), and the phyllosilicates Illite (Inter-ILI Engineering
206
Co. Ltd., Kosd, Hungary) and Montmorillonite (‘Ceratosil® WG’, Süd-Chemie AG,
207
Moosburg, Germany). All standard materials were attached to adhesive carbon foil
208
(polycarbonate graphite mixture; PLANO Leit-Tabs; Plano, Marburg, Germany) and ACS Paragon Plus Environment
Page 10 of 20
Page 11 of 20
Environmental Science & Technology
209
measured individually on the NanoSIMS 50L at the Lehrstuhl für Bodenkunde at the TU
210
München. All seven measurements were conducted under the same settings as the main
211
experiment and the following six masses were recorded:
212
adhesive carbon film, and 16O, 28Si, 27Al16O, 56Fe16O, and 56Fe16O2 for the identification of the
213
different mineral species (Figures S1). In test data set 1, the routine identified the seven pure
214
mineral grains on the adhesive carbon film and clearly discriminated the different mineral
215
species (Figures S1 and S2, table S1). Even the increasing Al substitution in the synthesized
216
goethite minerals could be visualized with this digital image processing routine.
217
Test data set 2: We prepared an artificial soil to evaluate the power of our routine to identify
218
pure materials in more complex soil-like mixtures, and especially to test the discrimination
219
between organic matter and resin. We built the artificial soil by mixing the two pure minerals
220
Quartz (Carl Roth, Karlsruhe, Germany) and Illite (Inter-ILI Engineering Co. Ltd., Kosd,
221
Hungary) with organic matter (needle from Norway Spruce (Picea abies L. Karst)), and
222
embedded this mixture in resin (Araldite kit 502, electron microscope sciences, Hatfield,
223
USA). The embedded sample (10 mm in diameter) was cured for 48 hours at 60°C, cut, and
224
subsequently polished in order to obtain a sample surface with low topography20. Prior to
225
NanoSIMS analysis, the sample was investigated using a reflectance light microscope (Axio
226
Imager Z2m, Zeiss, Germany) in order to identify 18 locations for the more detailed
227
NanoSIMS measurements. We used the same settings as the main experiment and the
228
following six masses were recorded: 12C and 12C14N for the identification of resin and organic
229
matter, and
230
species (Figures S3). In test data set 2, we could show that the discrimination of different
231
mineral species is possible in artificial and resin-embedded soil mixtures, too (Figures S3 and
232
S4, table S2). Furthermore, this test data set showed that resin and organic matter can be
233
discriminated based on their unique mass signatures, even in very complex cases in which
16
O,
28
Si,
27
Al16O, and
56
12
C for the identification of the
Fe16O for the identification of the different mineral
ACS Paragon Plus Environment
Environmental Science & Technology
234
resin infiltrated the internal structure of the spruce needle and lead to pixels with a mixed
235
mass signature showing characteristics of both organic matter and resin.
236
Results and discussion
237
By using digital image processing consisting of data preparation, endmember extraction, and
238
supervised classification, we were able to analyze all 40 images in one approach. We were
239
surprised that these 40 images did not show various spatial arrangements14, 15, but clustered in
240
just three distinct, recurring types of micro-architectures (Figure 2, table 1, figure S8). These
241
three types of micro-architectures can be distinguished by their specific mineral composition
242
(Figure S11), quantity and quality of organic matter, and their pore system. Our findings show
243
for the first time that mineral grains differing in size and mineral composition are not
244
randomly distributed in the soil matrix14-16 but are arranged in well-defined mineral
245
assemblages. Therefore, we define the combination of specific mineral assemblages with
246
different quantities and qualities of organic matter as arranged in distinct micro-architectures
247
as microdomains.
248
Microdomain A is characterized by weathering-resistant, coarse, Si-rich mineral grains (mass
249
signature classes 23-27; figure S7), glued together by thin layers of phyllosilicates (mass
250
signatures classes 15-22; for the separation of different types of phyllosilicates see figures S1
251
and S2 for NanoSIMS, and S11 for SEM-EDX experiments). This microdomain is low in
252
organic matter (mass signature classes 5 and 6), but high in porosity (mass signature classes
253
1-4; Figure 3, table 1). Microdomain B contains coarse, Si-rich mineral grains in a dense
254
matrix of phyllosilicates, and shows low porosity and intermediate amounts of organic matter.
255
In contrast to Microdomain A, here the coarse grains are connected to the dense phyllosilicate
256
matrix via thin Fe oxide coatings (mass signature classes 28-30). Microdomain C has low
257
porosity as it is built up by densely arranged Al-rich phyllosilicates (mass signature classes 7-
258
10) in association with high proportions of organic matter. Single, coarse, Al-rich mineral ACS Paragon Plus Environment
Page 12 of 20
Page 13 of 20
Environmental Science & Technology
259
grains (mass signature classes 11-14) remain as parent minerals embedded in the
260
phyllosilicate matrix. Overall, we classified 19 out of 40 images as Microdomain A and 13 as
261
Microdomain C. Four images showed the characteristics of Microdomain B and four images
262
were taken on large pure mineral grains without any micro-architecture (Table S4). It is
263
important to note that the hierarchical cluster analyses showed a clear discrimination between
264
the three microdomains and no smooth transition from one to the other (Figure S10).
265
We assume these microdomains with their unique micro-architectures to be the result of their
266
specific mineral assemblages that in turn induced different feedback mechanisms leading to
267
the formation and stabilization of the respective microdomain28. Large, Si-rich mineral grains,
268
most likely quartz crystals, are the main mineral components of microdomain A. These
269
smooth grains are weathering slowly and produce little, Si-rich phyllosilicates that bind the
270
coarse mineral grains into the soil structure via mineral-mineral interactions. Here, attractive
271
forces of physical origin (water menisci, swelling pressure) pull the particles together29 and
272
stabilize this porous micro-architecture. Binding via organic matter plays a minor role due to
273
low amounts and the low sorptive capacities of the mineral grains13. In Microdomain C, we
274
assume coarse, Al-rich mineral grains to have higher weathering rates producing plenty of Al-
275
rich phyllosilicates (Table 1 and figure 2) with high isomorphic substitution, surface charge,
276
and sorption capacity. This results in a microstructure dominated by a phyllosilicate matrix
277
that promotes accumulation of organic matter28, which in turn increases the association to
278
other reactive mineral surfaces and results in a dense packing with low porosity29, 30. This low
279
accessibility hampers the decomposition of the organic matter, keeping the structure intact,
280
and the organic matter protected. In Microdomain B similar mechanisms are active as in
281
Microdomain C with the exception that coarse, Si-rich mineral grains are embedded in the
282
dense phyllosilicate matrix. This is possible because thin Fe oxide coatings add sorption
283
potential to the otherwise low sorptive capacity of the Si-rich mineral grains and enable the
284
binding of these minerals in the phyllosilicate matrix31. The different feedback mechanisms ACS Paragon Plus Environment
Environmental Science & Technology
285
described above keep the three microdomains stable and support the concept of soil self-
286
organization forwarded by Young and Crawford32, and Smucker et al.33. Going beyond this
287
assessment of organic and mineral soil components and their interactions, it is now possible to
288
link their spatial arrangement with distinct pore characteristics.
289
The three microdomains are characterized by different pore systems, and in consequence
290
contribute differently to important soil functions, i.e. carbon sequestration, microbial habitat,
291
nutrient and contaminant exchange and storage, and water transport (Table 1). The pore
292
system of Microdomain A is characterized by large pore sizes between 0.2 and 14.7 µm
293
(mean = 7 µm), and high pore volumes (Figure 3). It allows the fast and advective transport
294
and exchange of water, nutrients, contaminants, and dissolved organic matter. We assume that
295
this pore system provides the living space for soil microorganisms7, 16, and its soil solution to
296
be available to plants. We did not see particulate organic matter, indicative for fresh or partly
297
decomposed litter material (Table 1, figure 2, and figure S8). We assume the particulate
298
organic matter in Microdomain A to be easily accessible and thus prone to microbial
299
degradation8, 34. Organic matter is found only in thin layers intimately associated with the
300
surfaces of phyllosilicates and may in this configuration be protected from degradation35-37.
301
Overall, the carbon sequestration in this microdomain is comparably low (Figure 3 and
302
table 1). We define microdomain A as the supply domain, holding low but readily available
303
amounts of water, nutrients, contaminants, and living biomass. In contrast, the microdomains
304
B and C are characterized by small pore sizes (0.2 to 7.2 µm, mean = 2 µm) and a generally
305
low porosity (Figure 3). The accessibility of the phyllosilicate surfaces with their high
306
sorption capacity is limited for nutrients, contaminants, and dissolved organic matter because
307
of restriction to diffusive transport38. Therefore, their exchange and cycling can be considered
308
as strongly retarded. The small pore sizes also hamper the access of microorganisms39, 40. This
309
reduced microbial accessibility promotes the accumulation of both, particulate and mineral
ACS Paragon Plus Environment
Page 14 of 20
Page 15 of 20
Environmental Science & Technology
310
bound organic matter. We define the microdomains B and C as the storage domains with high
311
but not readily available amounts of organic matter, nutrients and contaminants.
312
We were surprised to find recurring spatial arrangements of mineral and organic components,
313
and pores in 40 spatially independent NanoSIMS images on an intact soil aggregate.
314
Moreover, these micro-architectures were always associated with the same specific mineral
315
assemblages, characteristic amounts and forms of organic matter, and distinct pore systems.
316
Based on the extensive measurements we were able to statistically cluster these manifold
317
components in just three clearly different microdomains based on the contribution of the
318
different mineral and organic components. We assume that self-organization and feedback
319
mechanisms enable their development and keep these microdomains stable. It is striking that
320
these microdomains are associated with clearly different functions (Table 1). We suppose
321
these microdomains as the smallest spatial units of a soil that fulfil specific functions. Based
322
on these findings we propose these functional microdomains to be an adequate level for
323
generalization in soil modelling. Functional microdomains are the threshold where basic
324
physical, chemical, and biological properties reach a maximum level of complexity. On the
325
scale of functional microdomains, these basic properties can be standardized and used as
326
characteristic for a given type of microdomain. We suggest functional microdomains to be the
327
smallest building block of soil functionality and to define the lowest level of generalization as
328
proposed by Vereecken et al.9. We assume that the contribution of the different types of
329
functional microdomains is different in other soil materials and their variability represents the
330
variety of soil functionality. However, the surprising results steming from the extensive
331
measurement and the applied statistical methods used for their identification encourage us to
332
propose these functional microdomains as the smallest functional unit of soils and as an
333
important spatial scale in soil structure. We call for more studies on this scale that combine
334
different disciplines and data analyses techniques in order to proof or challenge our findings
335
of functional microdomains and their underlying micro-architectures. Our findings do not ACS Paragon Plus Environment
Environmental Science & Technology
336
redundantise studies on the smaller, molecular or coarser, aggregate-sized scales. We want to
337
set the threshold where generalization of soil materials can start. The characterization of
338
physico-chemical processes on the one and the understanding of the temporal more dynamic
339
aggregate-scale processes are vital for the understanding of soil systems. The functional
340
microdomains define the threshold where basic physico-chemical processes become functions
341
and where heterogeneity of soil functionality begins.
ACS Paragon Plus Environment
Page 16 of 20
Page 17 of 20
Environmental Science & Technology
342
References
343 344
1. Amundson, R.; Berhe, A. A.; Hopmans, J. W.; Olson, C.; Sztein, A. E.; Sparks, D. L., Soil and human security in the 21st century. Science 2015, 348, (6235).
345 346
2. Paustian, K.; Lehmann, J.; Ogle, S.; Reay, D.; Robertson, G. P.; Smith, P., Climatesmart soils. Nature 2016, 532, (7597), 49-57.
347 348 349
3. Totsche, K. U.; Rennert, T.; Gerzabek, M. H.; Koegel-Knabner, I.; Smalla, K.; Spiteller, M.; Vogel, H.-J., Biogeochemical interfaces in soil: The interdisciplinary challenge for soil science. Journal of Plant Nutrition and Soil Science 2010, 173, (1), 88-99.
350 351
4. Mikutta, R.; Kaiser, K., Organic matter bound to mineral surfaces: Resistance to chemical and biological oxidation. Soil Biology & Biochemistry 2011, 43, (8), 1738-1741.
352 353 354
5. Peth, S.; Chenu, C.; Leblond, N.; Mordhorst, A.; Garnier, P.; Nunan, N.; Pot, V.; Ogurreck, M.; Beckmann, F., Localization of soil organic matter in soil aggregates using synchrotron-based X-ray microtomography. Soil Biology & Biochemistry 2014, 78, 189-194.
355 356 357
6. Baveye, P. C.; Laba, M., Moving away from the geostatistical lamppost: Why, where, and how does the spatial heterogeneity of soils matter? Ecological Modelling 2015, 298, 2438.
358 359
7. Hallett, P. D.; Karim, K. H.; Bengough, A. G.; Otten, W., Biophysics of the Vadose Zone: From Reality to Model Systems and Back Again. Vadose Zone Journal 2013, 12, (4).
360 361 362
8. Kravchenko, A. N.; Negassa, W. C.; Guber, A. K.; Rivers, M. L., Protection of soil carbon within macro-aggregates depends on intra-aggregate pore characteristics. Scientific Reports 2015, 5.
363 364 365 366 367 368 369 370 371
9. Vereecken, H.; Schnepf, A.; Hopmans, J. W.; Javaux, M.; Or, D.; Roose, D. O. T.; Vanderborght, J.; Young, M. H.; Amelung, W.; Aitkenhead, M.; Allison, S. D.; Assouline, S.; Baveye, P.; Berli, M.; Bruggemann, N.; Finke, P.; Flury, M.; Gaiser, T.; Govers, G.; Ghezzehei, T.; Hallett, P.; Franssen, H. J. H.; Heppell, J.; Horn, R.; Huisman, J. A.; Jacques, D.; Jonard, F.; Kollet, S.; Lafolie, F.; Lamorski, K.; Leitner, D.; McBratney, A.; Minasny, B.; Montzka, C.; Nowak, W.; Pachepsky, Y.; Padarian, J.; Romano, N.; Roth, K.; Rothfuss, Y.; Rowe, E. C.; Schwen, A.; Simunek, J.; Tiktak, A.; Van Dam, J.; van der Zee, S.; Vogel, H. J.; Vrugt, J. A.; Wohling, T.; Young, I. M., Modeling Soil Processes: Review, Key Challenges, and New Perspectives. Vadose Zone Journal 2016, 15, (5).
372 373
10. Six, J.; Paustian, K., Aggregate-associated soil organic matter as an ecosystem property and a measurement tool. Soil Biology & Biochemistry 2014, 68, A4-A9.
374 375 376
11. Kleber, M.; Sollins, P.; Sutton, R., A conceptual model of organo-mineral interactions in soils: self-assembly of organic molecular fragments into zonal structures on mineral surfaces. Biogeochemistry 2007, 85, (1), 9-24.
377 378 379
12. Basile-Doelsch, I.; Balesdent, J.; Rose, J., Are Interactions between Organic Compounds and Nanoscale Weathering Minerals the Key Drivers of Carbon Storage in Soils? Environmental Science & Technology 2015, 49, (7), 3997-3998.
ACS Paragon Plus Environment
Environmental Science & Technology
380 381 382
13. Barre, P.; Fernandez-Ugalde, O.; Virto, I.; Velde, B.; Chenu, C., Impact of phyllosilicate mineralogy on organic carbon stabilization in soils: incomplete knowledge and exciting prospects. Geoderma 2014, 235, 382-395.
383 384 385
14. Lehmann, J.; Solomon, D.; Kinyangi, J.; Dathe, L.; Wirick, S.; Jacobsen, C., Spatial complexity of soil organic matter forms at nanometre scales. Nature Geoscience 2008, 1, (4), 238-242.
386 387 388 389
15. Solomon, D.; Lehmann, J.; Harden, J.; Wang, J.; Kinyangi, J.; Heymann, K.; Karunakaran, C.; Lu, Y.; Wirick, S.; Jacobsen, C., Micro- and nano-environments of carbon sequestration: Multi-element STXM-NEXAFS spectromicroscopy assessment of microbial carbon and mineral associations. Chemical Geology 2012, 329, 53-73.
390 391
16. Young, I. M.; Crawford, J. W.; Nunan, N.; Otten, W.; Spiers, A., Microbial distribution in soils: physics and scaling. Advances in Agronomy 2008, 100, 81-121.
392 393 394
17. Kölbl, A.; Kögel-Knabner, I., Content and composition of free and occluded particulate organic matter in a differently textured arable Cambisol as revealed by solid-state C-13 NMR spectroscopy. Journal of Plant Nutrition and Soil Science 2004, 167, (1), 45-53.
395 396 397
18. Mueller, C. W.; Weber, P. K.; Kilburn, M. R.; Hoeschen, C.; Kleber, M.; Pett-Ridge, J., Advances in the Analysis of Biogeochemical Interfaces: NanoSIMS to Investigate Soil Microenvironments. Advances in Agronomy, Vol 121 2013, 121, 1-46.
398 399 400 401
19. Mueller, C. W.; Remusat, L.; Rumpel, C., Characterization of Biogeochemical Processes at the Microscale: Concepts and Applications of NanoSIMS. In Terrestrial Ecosystem Research Infrastructures: Challenges and Opportunities, Chabbi, A.; Loescher, H. W., Eds. CRC Press: 2017; pp 193-212.
402 403 404
20. Mueller, C. W.; Höschen, C.; Kölbl, A.; Hillion, F.; Herrmann, A.; Heister, K.; KögelKnabner, I., Submicron scale imaging of soil organic matter dynamics using NanoSIMS from single particles to intact aggregates. Organic Geochemistry 2012, 42, (12), 1476-1488.
405 406
21. Price, J. C., How unique are spectral signatures. Remote Sensing of Environment 1994, 49, (3), 181-186.
407 408
22. Perona, P.; Malik, J., Scale-space and edge-detection using anisotropic diffusion. Ieee Transactions on Pattern Analysis and Machine Intelligence 1990, 12, (7), 629-639.
409 410
23. Wu, C. S., Normalized spectral mixture analysis for monitoring urban composition using ETM plus imagery. Remote Sensing of Environment 2004, 93, (4), 480-492.
411 412 413
24. Plaza, J.; Hendrix, E. M. T.; Garcia, I.; Martin, G.; Plaza, A., On Endmember Identification in Hyperspectral Images Without Pure Pixels: A Comparison of Algorithms. Journal of Mathematical Imaging and Vision 2012, 42, (2-3), 163-175.
414 415 416
25. Rogge, D. M.; Rivard, B.; Zhang, J.; Sanchez, A.; Harris, J.; Feng, J., Integration of spatial-spectral information for the improved extraction of endmembers. Remote Sensing of Environment 2007, 110, (3), 287-303.
417 418
26. Kruse, F. A.; Lefkoff, A. B.; Boardman, J. W.; Heidebrecht, K. B.; Shapiro, A. T.; Barloon, P. J.; Goetz, A. F. H., The spectral image-processing system (SIPS) - interactive
ACS Paragon Plus Environment
Page 18 of 20
Page 19 of 20
Environmental Science & Technology
419 420
visualization and analysis of imaging spectrometer data. Remote Sensing of Environment 1993, 44, (2-3), 145-163.
421 422
27. Schwertmann, U.; Cornell, R. M., Iron oxides in the laboratory: preparation and characterization. Second Edition ed.; John Wiley & Sons: 2000.
423 424 425 426
28. Fernandez-Ugalde, O.; Barre, P.; Hubert, F.; Virto, I.; Girardin, C.; Ferrage, E.; Caner, L.; Chenu, C., Clay mineralogy differs qualitatively in aggregate-size classes: clay-mineralbased evidence for aggregate hierarchy in temperate soils. European Journal of Soil Science 2013, 64, (4), 410-422.
427 428
29. Dexter, A. R., Advances in characterization of soil structure. Soil & Tillage Research 1988, 11, (3-4), 199-238.
429 430 431
30. Denef, K.; Six, J., Clay mineralogy determines the importance of biological versus abiotic processes for macroaggregate formation and stabilization. European Journal of Soil Science 2005, 56, (4), 469-479.
432 433
31. Tisdall, J. M.; Oades, J. M., Organic-Matter and Water-Stable Aggregates in Soils. Journal of Soil Science 1982, 33, (2), 141-163.
434 435
32. Young, I. M.; Crawford, J. W., Interactions and self-organization in the soil-microbe complex. Science 2004, 304, (5677), 1634-1637.
436 437 438
33. Smucker, A. J. M.; Park, E.-J.; Dorner, J.; Horn, R., Soil micropore development and contributions to soluble carbon transport within macroaggregates. Vadose Zone Journal 2007, 6, (2), 282-290.
439 440 441
34. Or, D.; Smets, B. F.; Wraith, J. M.; Dechesne, A.; Friedman, S. P., Physical constraints affecting bacterial habitats and activity in unsaturated porous media - a review. Advances in Water Resources 2007, 30, (6-7), 1505-1527.
442 443 444
35. Kleber, M.; Eusterhues, K.; Keiluweit, M.; Mikutta, C.; Mikutta, R.; Nico, P. S., Mineral–Organic Associations: Formation, Properties, and Relevance in Soil Environments. Advances in Agronomy 2015, 130, 1-140.
445 446 447 448
36. Koegel-Knabner, I.; Guggenberger, G.; Kleber, M.; Kandeler, E.; Kalbitz, K.; Scheu, S.; Eusterhues, K.; Leinweber, P., Organo-mineral associations in temperate soils: Integrating biology, mineralogy, and organic matter chemistry. Journal of Plant Nutrition and Soil Science 2008, 171, (1), 61-82.
449 450 451
37. Mikutta, R.; Mikutta, C.; Kalbitz, K.; Scheel, T.; Kaiser, K.; Jahn, R., Biodegradation of forest floor organic matter bound to minerals via different binding mechanisms. Geochimica Et Cosmochimica Acta 2007, 71, (10), 2569-2590.
452 453
38. Hartmann, A.; Gräsle, W.; Horn, R., Cation exchange processes in structured soils at various hydraulic properties. Soil & Tillage Research 1998, 47, (1-2), 67-72.
454 455
39. Sollins, P.; Homann, P.; Caldwell, B. A., Stabilization and destabilization of soil organic matter: Mechanisms and controls. Geoderma 1996, 74, (1-2), 65-105.
ACS Paragon Plus Environment
Environmental Science & Technology
456 457 458
Page 20 of 20
40. Dungait, J. A. J.; Hopkins, D. W.; Gregory, A. S.; Whitmore, A. P., Soil organic matter turnover is governed by accessibility not recalcitrance. Global Change Biology 2012, 18, (6), 1781-1796.
459
460
Acknowledgements
461
We are
462
Forschungsgemeinschaft (KO 1035/38-1). We kindly acknowledge financial support by the
463
Deutsche Forschungsgemeinschaft within the framework of the research unit “MAD Soil -
464
Microaggregates: Formation and turnover of the structural building blocks of soils” (DFG RU
465
2179) through project KO 1035/48-1. We thank Lars R. Zimmermann for designing figure 1,
466
Stefan Elgeti for SEM-EDX measurements, and Werner Häusler for XRD measurements.
467
Correspondence and requests should be addressed to M.S. Supplementary Information
468
accompanies this paper.
grateful for funding
of
the
NanoSIMS instrument by
ACS Paragon Plus Environment
the
Deutsche