Measuring Molecular Strain in Rewetted and Never-Dried Eucalypt

Jul 17, 2019 - of three leaning Eucalyptus bosistoana trees felled at Harewood, ... collected from each side of the samples at the same four positions...
0 downloads 0 Views 3MB Size
Subscriber access provided by BUFFALO STATE

Article

Measuring molecular strain in rewetted and neverdried eucalypt wood with Raman spectroscopy Fei Guo, and Clemens M. Altaner Biomacromolecules, Just Accepted Manuscript • DOI: 10.1021/acs.biomac.9b00808 • Publication Date (Web): 17 Jul 2019 Downloaded from pubs.acs.org on July 21, 2019

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 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 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.

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 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

1

Measuring molecular strain in rewetted and never-

2

dried eucalypt wood with Raman spectroscopy

3

Fei Guo1, 2, Clemens M. Altaner2*

4

AUTHOR ADDRESS

5

1 College

6

China

7

2 New

8

New Zealand

9

KEYWORDS: band shift, E. regnans, E. quadrangulata, growth stress, molecular deformation,

10

of Furnishings and Industrial Design, Nanjing Forestry University, 210037, Nanjing,

Zealand School of Forestry, University of Canterbury, Private Bag 4800, Christchurch,

tensile strain, PLS.

11

1 Environment ACS Paragon Plus

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

12

ABSTRACT

13

To measure growth strain in wood using Raman spectroscopy, we investigated the Raman spectra

14

of rewetted (water-saturated) E. regnans and green E. quadrangulata wood during tensile and

15

four-point bending tests. Partial least squares models to predict tensile strain were built from the

16

Raman spectra. The best model could predict tensile strain with a root mean square error of 427.5

17

>?: Apart from the widely reported band shift at 1095 cm@ upon mechanical strain, spectral

18

changes at 1420, 1120, 895, 456 cm@ were identified. The assignments of these bands were

19

discussed in relation to the molecular deformation of cellulose. The band shift rates during tensile

20

tests were @ :,. and @ : $ cm@ /% for rewetted E. regnans and green E. quadrangulata wood,

21

respectively. We successfully detected the release of molecular growth strain in green eucalyptus

22

wood with Raman spectroscopy by observing band shifts of the 1095 cm@ signal. Further, there

23

was a moderate correlation (r = 0.48) between growth-strain measured with strain gauges and the

24

1095 cm@ band position. The precision of the prediction of growth strain using Raman

25

spectroscopy was negatively affected by variation attributed to the inhomogeneity of wood on the

26

millimeter scale and instrumental instability.

27 28

1 INTRODUCTION

29

Raman and mid-infrared (MIR) spectroscopy are complementary techniques measuring

30

vibrational frequencies of chemical bonds 1. The former is based on the Raman scattering

31

phenomenon, while the latter detects the absorption of electromagnetic radiation. With

32

technological advances, Raman spectroscopy is gaining popularity in various applications

33

including material identification 2, forensic and pharmaceutical analysis 3 as well as chemical

2 Environment ACS Paragon Plus

Page 2 of 37

Page 3 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

34

distribution imaging of plant cell walls 4, 5. Portable Raman spectrometers are commercially

35

available 6, opening up the possibility of measuring wood properties in the field rapidly and non-

36

destructively. Compared to infrared spectroscopy, Raman spectroscopy is insensitive to water,

37

which interferes with the structural information of plant material 7 as it is typically water

38

saturated when alive. Challenges of using Raman spectroscopy on lignocellulosic materials are

39

the laser-excited florescence from lignin, which causes an intense background masking the

40

Raman signal, and the high energy of the laser, which leads to thermal decomposition of the

41

sample. Using 1064 nm near infrared (NIR) lasers dramatically reduces florescence from lignin 8

42

and water immersion sampling lessens thermal degradation 9.

43

Mechanical stress induces both macroscopic and molecular deformation of wood and cellulose.

44

Axial elongation and lateral contraction of the cellulose crystal lattice in wood was observed

45

upon tensile stress by X-ray diffraction 10, 11. Deformation of crystalline cellulose will cause

46

changes in vibrational frequencies of chemical bonds, which were observed as band shifts in

47

MIR 12, 13, near infrared (NIR) 14 and Raman spectra 15. IR spectroscopy detected the deformation

48

of the COC glycosidic bonds at 1160 cm@ and intramolecular O3H O5 groups at 3348 cm@

49

NIR band shifts at 6286 ± 5 cm@ and 6470 ± 10 cm@ indicated intramolecular hydrogen bonds

50

of crystalline cellulose elongated under tensile strain 14.

51

Raman spectroscopy has been a powerful tool to study the molecular deformation of polymer

52

fibers. 16 reported that the Raman frequencies assigned to the backbone structure of

53

polydiacetylene decreased with the elongation of fibers. Following this work, Raman

54

spectroscopy was used to study the molecular deformation of regenerated cellulose 17. The

55

Raman band shift at around 1095 cm@ , assigned to the stretching of the cellulose ring structure,

3 Environment ACS Paragon Plus

13.

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

56

has also been reported for wood, flax and hemp fibers 18-20. However, the effect of mechanical

57

strain on other regions in the Raman spectra is not well understood.

58

Growth stresses in trees, accumulated during cell wall maturation 21, cause problems like end-

59

splitting, heart checking and collapse and are the major cause restricting plantation-grown

60

eucalypts from in solid wood processing 22. Growth strain, similar to mechanical strain, causes

61

deformation of wood and cellulose on the molecular level 23. Methods to measure growth strain

62

in trees, the Nicholson technique 24, the CIRAD-Forêt method 25, log splitting 26 and the strain

63

gauge method 27, are destructive or time consuming. Rapid and non-destructive evaluation of

64

growth stress is needed to segregate logs with low growth strain for high-value solid wood uses.

65

Raman spectroscopy could potentially measure growth strain quickly and non-destructively by

66

accurate determination of vibrational frequencies. Previous studies on the deformation

67

mechanism of cellulosic materials focused on small areas of single fibers using a Raman

68

microscope, requiring elaborate sample preparation. However, Raman spectrometers capable of

69

direct surface sampling are available, enabling the collection of spectra without the limitation of

70

sample size. This is also advantageous as mechanical properties of wood have been known to

71

depend on sample size 28. Thin wood samples were shown to have different load-displacement

72

curves compared to bulk wood 29 and the modulus of elasticity of microtomed wood sections

73

were lower than that of normal sized specimens 30.

74

In this study, thick wood strips were mechanically stretched to confirm the effect of strain on

75

Raman spectra, previously only reported for microscopic samples, differing in their mechanical

76

behavior. Band shift rates were calculated to quantify molecular strain in rewetted and green

77

wood. Following on from this, the release of growth strain was monitored by Raman

4 Environment ACS Paragon Plus

Page 4 of 37

Page 5 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

78

spectroscopy on green, never-dried stems before finally evaluating the possibility of measuring

79

growth strain in green stems using Raman spectroscopy.

80 81

2 MATERIALS AND METHODS

82

2.1 Tensile tests of rewetted and green wood monitored by Raman spectroscopy

83

An air-dry Eucalyptus regnans board with an oven-dry density of 0.51 g/cm3 was sourced from

84

Canterbury, New Zealand. The microfibril angle (MFA) of the E. regnans wood was 7°

85

(standard deviation of 2.9°), as previously measured by X-ray diffraction. Ten wood strips with a

86

thickness of 0.5 to 1 mm (radial) were cut with a circular saw from the sapwood of the air-dry E.

87

regnans board and sanded before being cut into the shape shown in Figure 1. Wood sample used

88

for tensile tests while collecting Raman spectra. The grain (longitudinal direction) runs along the

89

longest dimension of the sample. Dimensions are labelled in millimeters.

90

. The samples were re-saturated through immersion in water for at least 7 days and were referred

91

to as “rewetted” wood in this study. A 2-year old E. quadrangulata tree with an oven-dry density

92

of 0.68 g/cm3 was harvested from a site in Woodville, New Zealand. The tree was kept in water

93

during transport to the lab. Eight wood strips were prepared from the never-dried green tree in

94

accordance to the procedure described above for the E. regnans samples. E. quadrangulata

95

samples were put in a sealed plastic bag and stored in a freezer until tested.

5 Environment ACS Paragon Plus

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

96 97

Figure 1. Wood sample used for tensile tests while collecting Raman spectra. The grain

98

(longitudinal direction) runs along the longest dimension of the sample. Dimensions are labelled

99

in millimeters.

100

A purpose-built tensile test rig 14 was placed under a dispersive Raman spectrometer (Sierra 2.0,

101

Metrohm Raman, USA) equipped with a 1064 nm NIR laser and an InGaAs array detector. The

102

wood samples were mechanically stretched while simultaneously collecting Raman spectra at

103

twelve strain levels, increasing progressively from “Strain 0” (approx. 0 >?< to “Strain 11”

104

(maximum strain of approx. 4500 >?< with increments of approximately 410 >?) before returning

105

to the relaxed (RL) state. Applied macroscopic strain was measured by strain gauges (FLA-5-11,

106

Tokyo Sokki Kenkyujo, Japan) glued (Loctite 454, Australia) onto the samples and connected to

107

a TC-31K strain-meter (Tokyo Sokki Kenkyujo, Japan). Strain was defined as the amount of

108

deformation per unit length and 1 >? equals to 1×10@. strain. The wood surface was kept wet

109

during spectra collection to restrict thermal decomposition 31 and improve spectral quality 9. Two

110

Raman spectra were taken at each strain level in the region from 2300 to 200 cm@ with a laser

111

power of 500 mW and an integration time of 30 to 35 s. In total 26 spectra were acquired for

112

each sample. The spectrometer featured orbital raster scan (ORS), a technique to reduce sample

113

damage by moving the laser over a measuring area of 1 to 2 mm in diameter. The polarization

114

direction of the incident light was parallel with the longitudinal direction of the wood samples.

6 Environment ACS Paragon Plus

Page 6 of 37

Page 7 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

ACS Paragon Plus Environment

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

Page 8 of 37

Page 9 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

151

2.4 Data analysis

152

In order to visualize the band shift caused by mechanical strain at around 1095 cm@ upon tensile

153

stress, Raman spectra were processed and plotted in R 33. Raman spectra were smoothed using

154

the Savitzky-Golay filter with a window length of 9 in the “signal” package 34 before being

155

normalized using the standard normal variate (SNV) method 35 in the range of 1180 to 950 cm@ .

156

For E. regnans wood upon tensile tests, PLS models were built to predict the applied mechanical

157

strain based on Raman spectra in the region between 2300 and 200 cm@ . The 258 spectra from

158

10 E. regnans samples were split randomly into a training group with 180 spectra and a

159

validation group with 78 spectra using the “sample” function in R. Five PLS models were

160

developed using leave-one-out (LOO) cross-validation with the “pls” package 36 based on

161

different spectra pre-processing, including baseline correction, SNV normalization, second

162

derivatives and spectra truncation (1180 ~ 950 cm@ ). Baseline correction was conducted with a

163

modified polynomial algorithm in the “baseline” package 37. Coefficients of determination (R2)

164

and root mean squared errors (RMSE) of the PLS models for the training and the validation

165

group were extracted.

166

Spectral changes caused by mechanical strain were also analyzed by the regression coefficients

167

plot of the PLS model and ‘slope’ spectra. ‘Slope’ spectra were calculated to show the rate of

168

change in spectral intensity at each wavenumber induced by mechanical strain as described

169

previously 14, 38. It was calculated from baseline corrected and normalized spectra. The average

170

‘slope’ spectrum was calculated as the average of 10 ‘slope’ spectra of 10 rewetted E. regnans

171

samples. To quantify the band shift at around 1095 cm@ , a Cauchy-Lorentzian function was used

9 Environment ACS Paragon Plus

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

172

to fit the spectra in the range from 1103 to 1084 cm@ using the nonlinear least squares method in

173

R. Band positions were extracted from the fitted model.

174 175

3 RESULTS AND DISCUSSIONS

176

3.1 Raman spectra of rewetted wood under tensile strain

177

3.1.1 Visualization of Raman band shifts

178

The signal intensity of Raman bands mainly depends on laser power and its frequency,

179

polarizability and concentration of the chemical group, as well as the sample orientation 39. For

180

biological materials, fluorescence and various physical effects such as surface roughness, particle

181

size and sample inhomogeneity have noticeable influence on the Raman scattering intensity as

182

well 40, 41. Due to these interfering factors, the intensities of the Raman spectra of the rewetted

183

wood fluctuated and it was difficult to identify the small spectral changes caused by mechanical

184

strain from the original spectra (Figure 4A). The fluctuations in the spectra could be removed by

185

SNV normalization in range from 1180 to 950 cm@ , allowing to determine the effect of

186

mechanical tensile strain on the position of the band centered at around 1095 cm@ (Figure 4A

187

inset). The 1095 cm@ band was assigned to CC and CO stretching of the cellulose ring structure

188

17, 42,

189

cellulose backbone. This band shift had been reported for regenerated cellulose as well as single

190

wood and other natural fibers using Raman microscopy 18, 19, 43-45. Our results showed that this

191

tensile strain induced Raman band shift can also be observed in relatively thick (>0.5 mm) E.

192

regnans wood strips, with mechanical behavior similar to solid wood 30. A strong linear

indicating the elongation of cellulose chains and weakening of the covalent bonds in the

10 Environment ACS Paragon Plus

Page 10 of 37

Page 11 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

193

relationship (r = 0.997) was found between the peak position and the applied strain (Figure 4B).

194

As strain increased from @$ >? to 4649 >?) the band position shifted from 1095.35 to 1093.57

195

cm@ . The band shift rate for this sample was @ :-4 cm@ /%, which was greater compared to the

196

values ;@ :

197

to 3.2.

to @ :. cm@ /%) reported for single wood fibers 19. For a detailed discussion refer

198 199

Figure 4 Unmodified Raman spectra at 12 tensile strain levels of one rewetted E. regnans wood

200

sample (A) and the corresponding band shifts of the 1095 cm@ band (B). Inset (A) shows the

201

band at 1095 cm@ after SNV normalization in the region of 1180 to 950 cm@ . “Strain 0” to

202

“Strain 3” were represented with solid lines, “Strain 4” to “Strain 7” with dashed lines and

203

“Strain 8” to “Strain 11” with dotted lines, each in the color sequence blue, red, pink and green.

204

The band at 1095 cm@ shifted 1.78 cm@ to lower frequencies with a strain increase of 4702 >?:

205

3.1.2 Predicting strain based on Raman spectra using PLS modelling

206

Table 1 shows the effects of spectral processing on the PLS models to predict mechanical strain

207

in wood. Baseline correction aimed to remove the effects of lignin fluorescence, while

208

normalization attempted to compensate multiplicative effects such as laser fluctuation and

11 Environment ACS Paragon Plus

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 37

209

uncertainties in focusing 40. Second derivative is another way to reduce baseline effects and

210

truncation to the region from 1180 to 950 cm@ eliminated many variables unaffected by

211

mechanical strain. The root mean square error of cross-validation (RMSECV) ranged from 354.7

212

to 493.1 >?) while the coefficient of determination of cross-validation (R2CV) varied between

213

0.89 and 0.94 (Table 1). Using the models with the validation data resulted in comparable

214

statistics, indicating no over-fitting. Baseline correction followed by normalization resulted in

215

the best performing PLS model (RMSEP = 397.2 >?< using 4 components. The predicted and

216

measured strain for the training group for this PLS model was shown in Figure 5A. The

217

distribution of residuals suggested an unbiased model (Figure 5B). The models based on

218

truncated spectra in the region of 1180 to 950 cm@ and unmodified spectra achieved slightly

219

reduced model statistics with RMSEP of 426.7 and 443.6 >?) respectively. However, more

220

components were required for the unmodified spectra (Table 1). PLS modelling coupled with

221

reflection NIR spectroscopy was successfully used to predict bending (R2 V 0.90) and tensile

222

load (R2 = 0.96) on small air-dry wood samples 46, 47. As NIR is sensitive to moisture, this

223

technique is of limited use for wood with high moisture content. Non-destructive reflection NIR

224

measurements of mechanical strain in green wood were found to be less accurate (R2 = 0.36),

225

due to the high moisture content interfering with structural information and a low signal to noise

226

ratio 7.

227

Table 1 Metrics of PLS models utilizing unmodified and processed Raman spectra to predict

228

mechanical strain in rewetted E. regnans wood. The spectral region was 2300 to 200 cm@ , except

229

for the truncated spectra, for which ranged from 1180 to 950 cm@ .

Pre-processing method

Training group Number of 2 components R CV RMSECV

12 Environment ACS Paragon Plus

Validation group R2P RMSEP

Page 13 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

Raw data Baseline correction Baseline correction + SNV Second derivative Truncation + SNV

0.91 0.91 0.94 0.89 0.93

440.5 432.8 354.7 493.1 385.9

6 4 4 3 3

0.92 0.92 0.94 0.90 0.93

443.6 456.1 397.2 510.5 426.7

230 231

SNV: standard normal variate (SNV) normalization; CV: Leave-one-out cross-validation;

232

RMSECV: root mean squared error of cross-validation using the training data; RMSEP: root mean

233

squared error of prediction for the validation data; R2: coefficient of determination.

234 235

Figure 5 Measured versus predicted tensile strain of rewetted E. regnans wood (A) and distribution

236

of the residuals (B) for the PLS model based on baseline corrected and SNV normalized Raman

237

spectra of the training data.

238

3.1.3 Spectral changes analyzed by regression coefficients and the average ‘slope’ spectrum

239

Both, regression coefficients of the PLS model and ‘slope’ spectra reflect spectral changes

240

caused by mechanical strain, but were calculated with different algorithms. The regression

241

coefficients of the PLS model represent the change in the response variable (i.e. tensile strain)

13 Environment ACS Paragon Plus

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

242

induced by one unit change in the predictor variables (i.e. the Raman intensities at each

243

wavenumber). As shown in Figure 6, the PLS regression coefficients and the average ‘slope’

244

spectrum showed almost identical features, with spectral changes in the regions of 1500 to 850

245

cm@ and 500 to 280 cm@ . Split peaks at 1420, 1120, 1095, 895, 456 cm@ with negative peaks

246

immediately followed by positive peaks indicated band shifts to lower frequencies 38, 48.

247

The most dominant Raman band shift induced by mechanical strain at around 1095 cm@ (Figure

248

4 and Figure 6) was reported for various (ligno)cellulosic materials 18-20. The band at 1127 cm@

249

was also previously reported to shift to lower wavenumbers for wood fibers under tensile strain

250

19.

251

deformation of the cellulose skeleton 5. Additionally, in regenerated cellulose fibers the band at

252

895 cm@ was also reported to shift upon stress, but with a lower sensitivity to strain than the

253

band at 1095 cm@

254

HCC at C6 of cellulose 9, 42, indicating that HCO and HCC at C6 were deformed upon tensile

255

strain. NIR spectroscopy observed the deformation of 2OH, involved in intramolecular hydrogen

256

bonds to O6, upon stretching 14, likely affecting the other covalent bonds connected to this group.

257

The region from 1430 to 1350 cm@ was assigned to the bending vibrations of COH groups 42.

258

Therefore, the Raman band shift at 1420 cm@ suggested the molecular deformation of COH

259

groups, which was consistent with MIR and NIR evidence of tensile strain induced deformation

260

of OH groups involved in hydrogen bonds 14, 49-51. Consistent with our observations, Raman band

261

shifts at 1414, 1095 and 895 cm@ were reported for high-modulus cellulose fibers under tensile

262

stress 52 and additionally at around 459 cm@ in single ramie fibers during tensile tests 53. The

263

region between 250 and 550 cm@ was mainly attributed to the skeletal bending motions of

264

cellulose, including CCC, COC, OCC and OCO bonds 42, 54.

Both bands were assigned to the stretching of cellulose CC and CO bonds, confirming the

17.

The band at 895 cm@ was tentatively assigned to the bending of HCO and

14 Environment ACS Paragon Plus

Page 14 of 37

Page 15 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

265 266

Figure 6 Regression coefficients (grey) of the PLS model based on baseline corrected and

267

normalized Raman spectra and the average ‘slope’ spectrum (black) for 10 rewetted E. regnans

268

wood samples upon stretching in the region from 2300 to 200 cm@ . The dashed line shows a

269

representative unmodified Raman spectrum of rewetted E. regnans wood.

270

3.2 Quantification of Raman band shifts caused by stretching

271

Apart from PLS modelling, another approach to predict strain is to use linear models based on

272

the positions of the Raman band at around 1095 cm@ . Raman band positions at around 1095

273

cm@ were determined by fitting a Cauchy-Lorentzian function to the signal. The dependence of

274

the Raman band position on tensile strain was quantified for rewetted E. regnans (Figure 7A)

275

and green E. quadrangulata (Figure 7B) wood (Figure 7). As tensile strain increased from the

276

relaxed state to the maximal applied strain level, the band at 1095 cm@ shifted to lower

277

wavenumbers for both, previously dried water-saturated and never-dried green wood, indicating

278

elongation of the cellulose ring structure. With the release of the applied tensile stress, the

15 Environment ACS Paragon Plus

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

279

Raman band returned to its original position. No significant differences were found between the

280

band positions before (“Strain 0”) and after the tensile tests (RL) (p = 0.31 and 0.27 for rewetted

281

and green wood, respectively). Similar behavior has been reported for single wood fibers 19 but

282

no similar studies have been done for bulk wood. This study showed that the Raman band at

283

around 1095 cm@ was also sensitive to strain for wood samples with a thickness of 0.5 to 1 mm.

284

Thin wood samples, only a few cells thick, behave differently from bulk wood during tensile

285

tests due to unrestricted cellular and supramolecular deformations caused by removing the

286

supporting structure of neighboring tissue 29. With eucalyptus fibers having diameters of 10 to 20

287

µm 55, this effect was minimized here, as the wood strips used in this study were comprised of

288

more than 25 cells in thickness. For even larger wood samples MOE was shown to increase with

289

sample thickness 28, 30.

290

16 Environment ACS Paragon Plus

Page 16 of 37

Page 17 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

291

Figure 7 Dependence of the position of the Raman band at around 1095 cm@ on tensile strain for

292

rewetted E. regnans (A) wood (n = 10) and green E. quadrangulata (B) wood (n = 8). Raman

293

spectra were collected as mechanical strain increased stepwise from “Strain 0” (approx. 0 >?< to

294

“Strain 11” (approx. 4500 >?< with a step of approx. 410 >? and then returned to the relaxed (RL)

295

state.

296

The mean Raman band position in the relaxed state for rewetted E. regnans and green E.

297

quadrangulata was 1095.41 (±0.15) and 1095.26 (±0.17) cm@ , respectively (Table 2). The

298

difference in band position between E. regnans and E. quadrangulata was not statistically

299

significant (p = 0.08). The observed variation could be caused by instrumentation like laser

300

fluctuation or variation in chemical composition and wood anatomy features such as microfibril

301

angle and grain angle 56. Additionally, the drying history was reported to affect the cellulose

302

structure 57, 58.

303

The band shift rate with respect to strain for rewetted E. regnans and green E. quadrangulata

304

during the tensile tests was @ :,. and @ : $ cm@ /% (Table 2). The significantly lower (p ?:

cm@ /% for air-dried hemp fibers 20 and @ :

Species

Sample condition

Sample size

E. regnans

rewetted

10

E. quadrangulata

green

8

to @ :. cm@ /% for air-dried single wood

Band position (cm@ ) 1095.41 (0.15) 1095.26 (0.17)

Band shift rate (cm@ /%)

R2

@ :,. (0.48)

0.95 (0.02)

@ : $ (0.19)

0.89 (0.05)

327 328

3.3 Monitoring the release of growth stress using Raman spectroscopy

329

Wood formed in the upper side of leaning hardwood trees, referred to as tension wood, is usually

330

associated with high tensile stress 62. Correspondingly, growth strain ranged from 1294 to 1446

18 Environment ACS Paragon Plus

Page 19 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

331

>? in the tension sides and from @ -3 to 64 >? in the opposite sides of the e investigated E.

332

bosistoana stems (Figure 8). Raman spectra were taken from these stems before and after

333

releasing growth stress. For the tension sides, the band position shifted to higher wavenumbers (p

334

< 0.016) after cutting the grooves, suggesting the release of tensile strain (Figure 8). In contrast,

335

no significant differences were found for the opposite sides (p > 0.140), attributed to the

336

relatively low stress levels in opposite wood and the detection limit of the Raman method.

337 338

Figure 8 Raman band positions of three leaning E. bosistoana trees before (open boxes) and after

339

(filled boxes) releasing growth stress. Positive strain values stand for tensile strain (contraction)

340

and negative values for compressive strain (expansion). “Side T” represents tension wood and

341

“Side O” opposite wood.

342

The variation in the Raman band positions was greater between the samples than within the

343

samples. The reasons are not fully understood, but could be related to the variation of chemical

344

composition 63 and anatomical structures 55 at the millimeter scale. The variation between

345

samples is a challenge for the prediction of growth strain based on Raman band positions.

346

3.4 Correlation of growth strain to the position of the 1095 cm@ Raman band

19 Environment ACS Paragon Plus

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

347

Methods to measure growth strain in trees, including the Nicholson technique 24, the CIRAD-

348

Forêt method 25 and the strain gauge method, had been reviewed by 64. The splitting test

349

proposed by 26 estimates growth strain of the stem based on the opening at the ends after splitting

350

the sample along the pith. This estimation will be affected by the three-dimensional distribution

351

of growth strain throughout the whole sample. In contrast, the other methods are a localized

352

measurement of surface growth strain. Although linear correlations were found between the

353

methods 26, 27, it is still challenging to build quantitative relationships between the different

354

methods 22.

355

As shown in Figure 9A, a strong linear relationship (r = 0.78) was found between the average

356

growth strain measured by strain gauges and the EGS using the splitting method. But the growth

357

strain measured by strain gauges was significantly lower compared to the EGS using the splitting

358

method. Additionally, the growth strain measured with the strain gauges at each side of the same

359

stem differed by up to more than four times (data not shown), confirming asymmetric SGS

360

around the stem 24.

361

The positions of the Raman band at around 1095 cm@ were moderately correlated (r = 0.48) with

362

growth strain measured using strain gauges (Figure 9B), but only a lower correlation (r = 0.27)

363

was found between the Raman band positions and the EGS (data not shown). This agreed with

364

the understanding that the Raman measurements reflected the localized molecular strain in the

365

stem and that this strain is variable throughout the stem 24.

366

The intercept of the linear fit was 1095.24 cm@ , consistent with the band position in the relaxed

367

state for green E. quadrangulata wood strips (Table 2). However, the slope of the linear fit was

368

@-: cm@ /%, more than twice the band shift rate of mechanically deformed samples (Table 2).

20 Environment ACS Paragon Plus

Page 20 of 37

Page 21 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

369

This suggested that the growth strain measured by strain gauges were underestimated, while the

370

band shift rate based on EGS calculated from the splitting test was @ :$ cm@ /%, similar to

371

those obtained by mechanical strain (Table 2). Similar discrepancies were reported between the

372

CIRAD-Forêt method and strain gauge measurements (Yoshida & Okuyama, 2002). The groove

373

depth and the distance to the strain gauge can affect the measured growth strain 27. The

374

surrounding wood tissue will impede the release of growth stresses 21. Additionally, according to

375

an X-ray diffraction study (Toba et al., 2013), molecular growth stress remained in a wood block

376

after being removed from the living stem. It is likely that not all growth stresses were released in

377

this experiment.

378 379

Figure 9 Relationship between growth strain measured by strain gauges and the estimated growth

380

strain (EGS) by the splitting method (A) and the dependence of the Raman band at around 1095

21 Environment ACS Paragon Plus

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

381

cm@ on growth strain measured by strain gauges (B). Each band position was calculated as the

382

average of 9 Raman spectra collected on the wood surface at the location of the strain gauge.

383 384

While this study has shown that molecular strain causing surface growth strain in stems can be

385

detected with Raman spectroscopy, technological advances are needed for Raman spectroscopy

386

to become commercially feasible for non-destructive segregation of eucalyptus logs for solid

387

wood processing. The relatively low correlation (r = 0.48) in the present study implied that

388

estimations of growth strain by Raman spectroscopy will result in falsely rejected usable logs.

389

This might be acceptable for a pulp chip business, aiming to increase revenue by identifying top

390

logs suitable for solid wood production, but is likely to be uneconomic for high value solid wood

391

plantations, which will lose revenue by falsely rejecting quality logs. Technological challenges

392

like measurement time and the focusing procedure of the dispersive system renders industrial

393

applications inefficient. Furthermore, assuming a band shift rate of @ cm@ /%, 1000 >? only

394

induces a band shift of 0.3 cm@ , making it necessary to address factors like instrumental

395

uncertainties, sample surface conditions, fluorescence, as well as variations in wood properties

396

and growth strain 24, 25, 65. Better prediction of growth strain can be expected with a Raman

397

instrument that can cover a larger surface area and averaging out local variations in chemical

398

composition and growth strain.

399 400 401

4 CONCLUSION

22 Environment ACS Paragon Plus

Page 22 of 37

Page 23 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

402

A partial least squares model based on pre-processed spectra was able to predict mechanical

403

strain with a root mean square error of 397.2 >? and R2 of 0.94. The major spectral changes

404

caused by mechanical strain were in the regions of 1500 to 850 cm@ and 500 to 280 cm@ . Apart

405

from the predominant band shift at 1095 cm@ , detailed analysis suggested that mechanical strain

406

also caused Raman band shifts at 1420, 1120, 895, 456 cm@ . The spectral changes could be

407

interpreted as the molecular deformation of the cellulose skeleton, HCO and HCC groups at C6

408

as well as COH structures. Further, studies with green wood indicated that the Raman band shift

409

rate with respect to strain was lower in never-dried green wood than previously dried rewetted

410

wood. The release of growth strain was successfully detected by Raman spectroscopy based on

411

the band positions at around 1095 cm@ . A moderate correlation was found between the band

412

positions at around 1095 cm@ and the growth strain measured by strain gauges, indicating that

413

the technology can potentially be used to segregate logs for solid wood production in certain

414

circumstances. However, improvements in Raman instrumentation, allowing more precise

415

measurements and spatial averaging, are needed.

416 417

AUTHOR INFORMATION

418

Corresponding Author

419

*New Zealand School of Forestry, University of Canterbury, Private Bag 4800, Christchurch,

420

New Zealand.

421

Tel: +64 33695949. Email: [email protected].

422

Author Contributions

23 Environment ACS Paragon Plus

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

423

The manuscript was written through contributions of both authors. FG conducted the

424

experiments and drafted the manuscript. CA provided guidance to the experimental design and

425

manuscript revision. All authors have given approval to the final version of the manuscript.

426

ACKNOWLEDGMENT

427

This work was financially supported by the Ministry of Business, Innovation and Employment,

428

New Zealand, through the Specialty Wood Products Partnership (FFRX1501). Nigel Pink and

429

Gert Hendriks (University of Canterbury) helped with the preparation of the wood samples.

430 431

REFERENCES

432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458

1. Larkin, P., Infrared and raman spectroscopy: principles and spectral interpretation. Elsevier: 2011. 2. Amer, M., Raman spectroscopy for soft matter applications. John Wiley & Sons: 2009. 3. Vankeirsbilck, T.; Vercauteren, A.; Baeyens, W.; Van der Weken, G.; Verpoort, F.; Vergote, G.; Remon, J. P., Applications of Raman spectroscopy in pharmaceutical analysis. TrAC, Trends Anal. Chem. 2002, 21 (12), 869-877. 4. Richter, S.; Müssig, J.; Gierlinger, N., Functional plant cell wall design revealed by the Raman imaging approach. Planta 2011, 233 (4), 763-772. 5. Gierlinger, N.; Schwanninger, M., Chemical imaging of poplar wood cell walls by confocal Raman microscopy. Plant Physiol. 2006, 140 (4), 1246-1254. 6. Carron, K.; Cox, R., Qualitative analysis and the answer box: A perspective on portable Raman spectroscopy. Anal. Chem. 2010, 82 (9), 3419-3425. 7. Guo, F. Molecular deformation of wood and cellulose studied by near infrared and Raman spectroscopy. University of Canterbury, Christchurch, New Zealand, 2019. 8. Agarwal, U. P., 1064 nm FT-Raman spectroscopy for investigations of plant cell walls and other biomass materials. Front. Plant Sci. 2014, 5 (490). 9. Agarwal, U. P., An overview of Raman spectroscopy as applied to lignocellulosic materials. In Advances in lignocellulosics characterization, Argyropoulos, D. S., Ed. TAPPI Press: 1999; pp 201-225. 10. Peura, M.; Kölln, K.; Grotkopp, I.; Saranpää, P.; Müller, M.; Serimaa, R., The effect of axial strain on crystalline cellulose in Norway spruce. Wood Sci. Technol. 2007, 41 (7), 565-583. 11. Zabler, S.; Paris, O.; Burgert, I.; Fratzl, P., Moisture changes in the plant cell wall force cellulose crystallites to deform. J. Struct. Biol. 2010, 171 (2), 133-141. 12. Åkerholm, M.; Salmén, L., Interactions between wood polymers studied by dynamic FTIR spectroscopy. Polymer 2001, 42 (3), 963-969. 13. Salmén, L.; Bergström, E., Cellulose structural arrangement in relation to spectral changes in tensile loading FTIR. Cellulose 2009, 16 (6), 975-982.

24 Environment ACS Paragon Plus

Page 24 of 37

Page 25 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504

Biomacromolecules

14. Guo, F.; Altaner, C. M., Molecular deformation of wood and cellulose studied by near infrared spectroscopy. Carbohydr. Polym. 2018, 197, 1-8. 15. Eichhorn, S. J.; Young, R. J., The Young's modulus of a microcrystalline cellulose. Cellulose 2001, 8 (3), 197-207. 16. Mitra, V. K.; Risen Jr, W. M.; Baughman, R. H., A laser Raman study of the stress dependence of vibrational frequencies of a monocrystalline polydiacetylene. J. Chem. Phys. 1977, 66 (6), 2731-2736. 17. Hamad, W. Y.; Eichhorn, S., Deformation micromechanics of regenerated cellulose fibers using Raman spectroscopy. J. Eng. Mater. Technol. 1997, 119 (3), 309-313. 18. Eichhorn, S. J.; Sirichaisit, J.; Young, R. J., Deformation mechanisms in cellulose fibres, paper and wood. J. Mater. Sci. 2001, 36 (13), 3129-3135. 19. Gierlinger, N.; Schwanninger, M.; Reinecke, A.; Burgert, I., Molecular changes during tensile deformation of single wood fibers followed by Raman microscopy. Biomacromolecules 2006, 7 (7), 2077-2081. 20. Peetla, P.; Schenzel, K. C.; Diepenbrock, W., Determination of mechanical strength properties of hemp fibers using near-infrared Fourier transform Raman microspectroscopy. Appl. Spectrosc. 2006, 60 (6), 682-691. 21. Alméras, T.; Clair, B., Critical review on the mechanisms of maturation stress generation in trees. J. R. Soc. Interface 2016, 13 (122), 20160550. 22. Yang, J. L.; Waugh, G., Growth stress, its measurement and effects. Aust. For. 2001, 64 (2), 127-135. 23. Clair, B.; Alméras, T.; Pilate, G.; Jullien, D.; Sugiyama, J.; Riekel, C., Maturation stress generation in poplar tension wood studied by synchrotron radiation microdiffraction. Plant Physiol. 2011, 155 (1), 562-570. 24. Nicholson, J. E., A rapid method for estimating longitudinal growth stresses in logs. Wood Sci. Technol. 1971, 5 (1), 40-48. 25. Gerard, J.; Bailleres, H.; Fournier, M.; Thibaut, B., Wood quality in plantation Eucalyptus - a study of variation in three reference properties. Bois For. Trop. 1995, 245, 101110. 26. Chauhan, S. S.; Entwistle, K., Measurement of surface growth stress in Eucalyptus nitens Maiden by splitting a log along its axis. Holzforschung 2010, 64 (2), 267-272. 27. Yoshida, M.; Okuyama, T., Techniques for measuring growth stress on the xylem surface using strain and dial gauges. Holzforschung 2002, 56 (5), 461. 28. Buchelt, B.; Pfriem, A., Influence of wood specimen thickness on its mechanical properties by tensile testing: solid wood versus veneer. Holzforschung 2011, 65 (2), 249-252. 29. Navi, P.; Rastogi, P. K.; Gresse, V.; Tolou, A., Micromechanics of wood subjected to axial tension. Wood Sci. Technol. 1995, 29 (6), 411-429. 30. Yu, Y.; Jiang, Z.-h.; Tian, G.-l., Size effect on longitudinal MOE of microtomed wood sections and relevant theoretical explanation. For. Stud. China 2009, 11 (4), 243-248. 31. Conti, C.; Botteon, A.; Bertasa, M.; Colombo, C.; Realini, M.; Sali, D., Portable sequentially shifted excitation Raman spectroscopy as an innovative tool for in situ chemical interrogation of painted surfaces. Analyst 2016, 141 (15), 4599-4607. 32. Davies, N. T.; Apiolaza, L. A.; Sharma, M., Heritability of growth strain in Eucalyptus bosistoana: a Bayesian approach with left-censored data. N. Z. J. For. Sci. 2017, 47, 5. 33. R Core Team, A language and environment for statistical computing. In R foundation for statistical computing, Vienna, Austria, 2017.

25 Environment ACS Paragon Plus

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549

34. Signal developers, Signal: Signal processing. 2013. 35. Stevens, A.; Ramirez-Lopez, L., An introduction to the prospectr package. In R package version 0.1.3, 2013. 36. Mevik, B.-H.; Wehrens, R.; Liland, K. H., Pls: Partial least squares and principal component regression. In R package version 2.6-0, 2016. 37. Liland, K. H.; Mevik, B.-H., Baseline: Baseline xorrection of spectra. In R package version 1.2-1, 2015. 38. Altaner, C. M.; Horikawa, Y.; Sugiyama, J.; Jarvis, M. C., Cellulose (c investigated by IR-spectroscopy at low temperatures. Cellulose 2014, 21 (5), 3171-3179. 39. Smith, E.; Dent, G., Modern Raman spectroscopy: a practical approach. John Wiley & Sons, Ltd: 2005. 40. Liland, K. H.; Kohler, A.; Afseth, N. K., Model-based pre-processing in Raman spectroscopy of biological samples. J. Raman Spectrosc. 2016, 47 (6), 643-650. 41. Huang, J.; Romero-Torres, S.; Moshgbar, M., Practical considerations in data pretreatment for NIR and Raman spectroscopy. American Pharmaceutical Review 2010, 13 (6), 116. 42. Wiley, J. H.; Atalla, R. H., Band assignments in the raman spectra of celluloses. Carbohydr. Res. 1987, 160 (Supplement C), 113-129. 43. Hsieh, Y.-C.; Yano, H.; Nogi, M.; Eichhorn, S. J., An estimation of the Young’s modulus of bacterial cellulose filaments. Cellulose 2008, 15 (4), 507-513. 44. Šturcová, A.; Davies, G. R.; Eichhorn, S. J., Elastic modulus and stress-transfer properties of tunicate cellulose whiskers. Biomacromolecules 2005, 6 (2), 1055-1061. 45. Eichhorn, S. J.; Young, R. J.; Davies, G. R., Modeling crystal and molecular deformation in regenerated cellulose fibers. Biomacromolecules 2005, 6 (1), 507-513. 46. André, N.; Labbé, N.; Rials, T. G.; Kelley, S. S., Assessment of wood load condition by Near Infrared (NIR) spectroscopy. J. Mater. Sci. 2006, 41 (7), 1879-1886. 47. Sandak, J.; Sandak, A.; Pauliny, D.; Krasnoshlyk, V.; Hagman, O., Near infrared spectroscopy as a tool for estimation of mechanical stresses in wood. Adv. Mater. Res. (DurntenZurich, Switz.) 2013, 778, 448-453. 48. Åkerholm, M.; Hinterstoisser, B.; Salmén, L., Characterization of the crystalline structure of cellulose using static and dynamic FT-IR spectroscopy. Carbohydr. Res. 2004, 339 (3), 569-578. 49. Hinterstoisser, B.; Åkerholm, M.; Salmén, L., Load distribution in native cellulose. Biomacromolecules 2003, 4 (5), 1232-1237. 50. Altaner, C. M.; Thomas, L. H.; Fernandes, A. N.; Jarvis, M. C., How cellulose stretches: synergism between covalent and hydrogen bonding. Biomacromolecules 2014, 15 (3), 791-798. 51. Kong, K.; Eichhorn, S. J., The influence of hydrogen bonding on the deformation micromechanics of cellulose fibers. J. Macromol. Sci., Part B: Phys. 2005, 44 (6), 1123-1136. 52. Eichhorn, S. J.; Young, R. J.; Davies, R. J.; Riekel, C., Characterisation of the microstructure and deformation of high modulus cellulose fibres. Polymer 2003, 44 (19), 59015908. 53. Gierlinger, N., New insights into plant cell walls by vibrational microspectroscopy. Appl. Spectrosc. Rev. 2018, 53 (7), 517-551. 54. Österberg, M.; Schmidt, U.; Jääskeläinen, A.-S., Combining confocal Raman spectroscopy and atomic force microscopy to study wood extractives on cellulose surfaces. Colloids and Surfaces A: Physicochemical and Engineering Aspects 2006, 291 (1), 197-201.

26 Environment ACS Paragon Plus

Page 26 of 37

Page 27 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577

Biomacromolecules

55. Wilkes, J., Variations in wood anatomy within species of Eucalyptus. IAWA J 1988, 9 (1), 13-23. 56. Agarwal, U. P., Raman imaging to investigate ultrastructure and composition of plant cell walls: distribution of lignin and cellulose in black spruce wood (Picea mariana). Planta 2006, 224 (5), 1141. 57. Hill, S. J.; Kirby, N. M.; Mudie, S. T.; Hawley, A. M.; Ingham, B.; Franich, R. A.; Newman, R. H., Effect of drying and rewetting of wood on cellulose molecular packing. Holzforschung 2010, 64 (4), 421-427. 58. Toba, K.; Yamamoto, H.; Yoshida, M., Mechanical interaction between cellulose microfibrils and matrix substances in wood cell walls induced by repeated wet-and-dry treatment. Cellulose 2012, 19 (4), 1405-1412. 59. Agarwal, U. P.; Ralph, S. A.; Reiner, R. S.; Baez, C., Probing crystallinity of neverdried wood cellulose with Raman spectroscopy. Cellulose 2016, 23 (1), 125-144. 60. Keckes, J.; Burgert, I.; Fruhmann, K.; Muller, M.; Kolln, K.; Hamilton, M.; Burghammer, M.; Roth, S. V.; Stanzl-Tschegg, S.; Fratzl, P., Cell-wall recovery after irreversible deformation of wood. Nat. Mater. 2003, 2 (12), 810-813. 61. Green, D. W.; Winandy, J. E.; Kretschmann, D. E., Mechanical properties of wood. In Wood handbook: wood as an engineering material, Ross, R. J., Ed. Department of Agriculture, Forest Service, Forest Products Laboratory: Madison, Wisconsin, 2010. 62. Gardiner, B.; Barnett, J.; Saranpää, P.; Gril, J., The biology of reaction wood. Springer: 2014. 63. Downes, G. M.; Touza, M.; Harwood, C.; Wentzel-Vietheer, M., NIR detection of nonrecoverable collapse in sawn boards of Eucalyptus globulus. Eur. J. Wood Wood Prod. 2014, 72 (5), 563-570. 64. Yang, J. L.; Baillères, H.; Okuyama, T.; Muneri, A.; Downes, G., Measurement methods for longitudinal surface strain in trees: a review. Aust. For. 2005, 68 (1), 34-43. 65. Yang, J. L.; Fife, D.; Matheson, A. C., Growth strain in three provenances of plantationgrown Eucalyptus globulus Labill. Aust. For. 2001, 64 (4), 248-256.

578

27 Environment ACS Paragon Plus

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

159x80mm (600 x 600 DPI)

ACS Paragon Plus Environment

Page 28 of 37

Page 29 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

Figure 1. Wood sample used for tensile tests while collecting Raman spectra. The grain (longitudinal direction) runs along the longest dimension of the sample. Dimensions are labelled in millimeters. 83x37mm (600 x 600 DPI)

ACS Paragon Plus Environment

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 2. Schematic drawing of the experimental setup for monitoring the release of growth strain by Raman spectroscopy on green stems. The dashed red lines indicate the grooves cut to release growth stress. T: tension wood side; O: opposite wood side. 85x37mm (600 x 600 DPI)

ACS Paragon Plus Environment

Page 30 of 37

Page 31 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

Figure 3. Illustration of evaluating growth strain of green stems using Raman spectroscopy. In total, 18 spectra were collected from each sample. The dashed red line indicates the cut for the splitting test. 85x34mm (600 x 600 DPI)

ACS Paragon Plus Environment

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4 Unmodified Raman spectra at 12 tensile strain levels of one rewetted E. regnans wood sample (A) and the corresponding band shifts of the 1095 cm−1 band (B). Inset (A) shows the band at 1095 cm−1 after SNV normalisation in the region of 1180 to 950 cm−1. “Strain 0” to “Strain 3” were represented with solid lines, “Strain 4” to “Strain 7” with dashed lines and “Strain 8” to “Strain 11” with dotted lines, each in the color sequence blue, red, pink and green. The band at 1095 cm−1 shifted 1.78 cm−1 to lower frequencies with a strain increase of 4702 µε. 169x67mm (600 x 600 DPI)

ACS Paragon Plus Environment

Page 32 of 37

Page 33 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

Figure 5 Measured versus predicted tensile strain of rewetted E. regnans wood (A) and distribution of the residuals (B) for the PLS model based on baseline corrected and SNV normalized Raman spectra of the training data 219x92mm (600 x 600 DPI)

ACS Paragon Plus Environment

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 6 Regression coefficients (grey) of the PLS model based on baseline corrected and normalized Raman spectra and the average ‘slope’ spectrum (black) for 10 strained rewetted E. regnans wood samples in the region from 2300 to 200 cm−1. The dashed line shows a representative unmodified Raman spectrum of rewetted E. regnans wood. 85x85mm (600 x 600 DPI)

ACS Paragon Plus Environment

Page 34 of 37

Page 35 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

Figure 7 Dependence of the position of the Raman band at around 1095 cm−1 on tensile strain for rewetted E. regnans (A) wood (n = 10) and green E. quadrangulata (B) wood (n = 8). Raman spectra were collected as mechanical strain increased stepwise from “Strain 0” (approx. 0 µε) to “Strain 11” (approx. 4500 µε) with a step of approx. 410 µε and then returned to the relaxed (RL) state. 85x103mm (600 x 600 DPI)

ACS Paragon Plus Environment

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 8 Raman band positions of three leaning E. bosistoana trees before (open boxes) and after (filled boxes) releasing growth stress. Positive strain values stand for tensile strain (contraction) and negative values for compressive strain (expansion). “Side T” represents tension wood and “Side O” opposite wood. 169x66mm (600 x 600 DPI)

ACS Paragon Plus Environment

Page 36 of 37

Page 37 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

Figure 9 Relationship between growth strain measured by strain gauges and the estimated growth strain (EGS) by the splitting method (A) and the dependence of the Raman band at around 1095 cm−1 on growth strain measured by strain gauges (B). Each band position was calculated as the average of 9 Raman spectra collected on the wood surface at the location of the strain gauge. 85x103mm (600 x 600 DPI)

ACS Paragon Plus Environment