Identification of Distinct Functional Microstructural ... - ACS Publications

Sep 28, 2017 - (NanoSIMS) and a novel digital image processing routine adapted from remote ..... goethite was synthesized following the procedure of S...
1 downloads 9 Views 548KB Size
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