Experimental techniques for characterising water in wood covering the range from dry to fully water-saturated

Water plays a central role in wood research, since it affects all material properties relevant to the performance of wood materials. Therefore, experimental techniques for characterising water within wood are an essential part of nearly all scientific investigations of wood materials. This review focuses on selected experimental techniques that can give deeper insights into various aspects of water in wood in the entire moisture domain from dry to fully water-saturated. These techniques fall into three broad categories: (1) gravimetric techniques that determine how much water is absorbed, (2) fibre saturation techniques that determine the amount of water within cell walls, and (3) spectroscopic techniques that provide insights into chemical wood–water interactions as well as yield information on water distribution in the macro-void wood structure. For all techniques, the general measurement concept is explained, its history in wood science as well as advantages and limitations.


Introduction
Wood is a hygroscopic material which means that it absorbs and exchanges water molecules with the surroundings.Water molecules found within cell walls interfere with internal secondary wood-wood bonds, for example hydrogen bonds.Therefore, it is no surprise that absorbed water, often termed moisture, has a profound effect on wood properties dependent on these bonds, for example strength and stiffness (Tiemann 1906;Wagner et al. 2015).However, moisture also influences the performance of wood and wood products in other ways important for the structural application of these materials such as durability, i.e. resistance to biodegradation (Ammer 1963a;Meyer and Brischke 2015;Stamm and Baechler 1960;Stienen et al. 2014;Thybring 2013).The central role of moisture in the performance of wood materials makes the wood-water interaction an obvious target for enhancing material performance to specific applications.Thus, a wide range of physicochemical modification processes exist aimed at improving the performance of wood materials.Common targets for wood modification include improved durability and dimensional stability (Stamm and Baechler 1960).However, recent years have seen an increase in advanced modifications adding novel functionalities to wood such as stimuli-responsive properties, for example ''sweating wood'' for cooling purposes or stimuli-responsive filtration (Keplinger et al. 2016).For all these types of modifications, wood-water interactions within the modified material are important for the performance.For instance, the resistance to biological decay is directly linked to the water capacity of modified cell walls (Stamm and Baechler 1960;Thybring 2013), while advanced modifications such as stimuli-responsive hydrogel formation within the wood structure use water uptake or release for cooling purposes by water evaporation or for filtration (Keplinger et al. 2016).To advance the understanding of the fundamental mechanisms of wood performance as affected by water as well as for tailoring wood modifications to specific applications, characterisation of the amount, distribution, and physicochemical interactions of water within wood materials is important.
This review covers selected experimental techniques for characterising moisture in wood in the entire moisture content range from absolutely dry to fully watersaturated, i.e.where water takes up all available space in cell walls and voids (pits, lumina, vessels) in the wood structure.Furthermore, techniques which elucidate other fundamental aspects of wood-water interactions in untreated and modified wood are included.Experimental techniques which concern effects of water on, for example, mechanics or dimensions will not be discussed.The schematic sorption isotherm shown in Fig. 1 covering the entire relative humidity (RH) range provides an overview of the characterisation techniques covered in this review as well as the RH range to which each technique can be applied.

Terminology
One of the most fundamental and commonly used descriptors for moisture in wood is the moisture content (MC).In this review, it is described by water mass divided by dry wood mass.When dealing with physicochemically modified wood materials, it is often necessary to adjust MC for added dry mass due to the modification process (Hill 2008;Thybring 2013) in order to correctly compare a given modified material with untreated or differently modified materials.This is done by multiplying MC of the modified material, MC mod [%], with (1 ?R mod ) as shown in (1): where MC R [%] is adjusted moisture content also termed ''reduced moisture content'', m water [kg] is water mass, and m wood,dry [kg] and m mod,dry [kg] are dry masses of the wood itself and the modified wood material, respectively.R mod [-] is the modification ratio defined as relative mass gain due to modification R mod ¼ m mod;dry À m wood;dry m wood;dry Mass change is commonly used in the literature to describe the intensity of modification either by the term ''weight per cent gain'' (WPG) [% WPG] in cases of increased mass due to modification, for example acetylation, furfurylation, etc., or by the term ''mass loss'' (ML) [% ML] in cases of decreased mass due to modification, for example thermal modification.It was chosen to describe the mass change by the common term ''modification ratio'' (R mod ) and define it as a dimensionless fraction to avoid confusion with the relative mass gain given in per cent (WPG).Correction of MC in (1) and following equations should, however, only be done in cases where R mod is positive.Thus, MC should not be corrected for thermally modified wood where R mod \ 0.
As water molecules absorbed by wood are continuously exchanged with the surroundings as well as meandering around within the material, MC is often characterised in the equilibrium state where there is no net change in water with the surroundings.This alleviates issues of comparing various samples with differences in moisture gradients depending on sample geometry, wood density, etc.In case of equilibrium, it is common to use the terms ''equilibrium moisture content'' (EMC) Fig. 1 Schematic illustration of moisture uptake in wood in absorption from absolute dry state (lower black line) and desorption from fully water-saturated state (upper black line) in the entire relative humidity (RH) range from dry to water-saturated.Characterisation (gravimetric, fibre saturation, and spectroscopic) techniques described in the present review are indicated along with the RH range they cover

Gravimetric techniques
Wood can absorb significant amounts of water in its structure, i.e. both within cell walls and in the void structure.The maximum amount depends on the wood density, which varies widely between species.As an example, a low-density wood species like balsa (Ochroma lagopus Sw.) with a bulk density around 100 kg/m 3 can fit a total amount of water of around 960-980 kg/m 3 in its bulk volume, equal to a maximum MC of 960-980%.On the other hand, a high-density wood species such as ironwood (Lophira alata Banks ex.Gaertn.) with a bulk density around 1000 kg/ m 3 can only fit about 630-780 kg/m 3 water inside its bulk volume, yielding a maximum MC of just 63-78%.The noticeable change in mass associated with absorption of water makes gravimetric techniques ideal for determining the amount of water inside the bulk material.

Early-stage sorption balances, salt solutions, and climate chambers
In the early scientific investigations of wood materials, it is recognised that the MC of wood exposed to climatic conditions changes with changes in temperature and RH (Forest Products Laboratory 1919;Glass et al. 2014;Zeller 1920).In the early days of wood research, it was, however, a common practice to achieve a specific MC in wood samples by drying the material from green or water-saturated conditions (Carrington 1922;Record 1914;Tiemann 1906;Volbehr 1896) and then maybe leaving the wood for moisture to redistribute within the material before further testing, for example mechanical.Alternatively, wood samples were kept under laboratory conditions for a period of time to obtain an MC around 12% (Carrington 1921).The first published study in a scientific journal on wood MC in equilibrium with a controlled, stable climate appears to be the work of Zeller (1920), who used sulphuric acid solutions of various concentrations to generate various levels of RH in a specially designed climate chamber.Similar investigations performed at the Forest Products Laboratory in the USA might, however, predate this publication (Glass et al. 2014).In these early investigations, the now wellknown effect of sorption history, i.e. that conditioning wood to a specific RH through moisture uptake (absorption) or through moisture release (desorption) will not produce similar MC, was not considered.Therefore, the relation between wood MC and RH at constant temperature was described by a single curve, i.e. sorption isotherm.The seemingly first absorption and desorption isotherms were reported by Larian and co-workers (Larian et al. 1930;Lavine and Gauger 1930) who conditioned birch shavings (about 1.5 g each) in desiccators over saturated salt solutions generating different RH levels.Initially, all samples were in the green condition, and after equilibrium with the specific RH had been obtained, they were dried over pure sulphuric acid (0% RH) followed by conditioning by absorption to the same level of RH as before.These data can therefore be regarded as the first desorption and absorption isotherms for wood, where the previous is initiated from a condition close to fully water-saturated.The same year, Pidgeon and Maass (1930) reported on absorption and desorption isotherms for white spruce and jack pine found by exposing small samples (150 mg) suspended from a quartz helix spring inside a vacuum tube system to various levels of RH.Desorption was, however, initiated from hygroscopic conditions (\ 98% RH), the effect of which will be discussed later.RH in the vacuum tube system was controlled by tuning the vapour pressure above a water bath via controlling the bath temperature while keeping the rest of the system at constant temperature.Moisture uptake was measured by the spring elongation, and the system is a precursor to modern-day automated sorption balances.Several later studies on wood use this kind of technique for studying absorption and desorption in wood (Christensen 1959(Christensen , 1960;;Christensen and Kelsey 1959;Grace and Maass 1932;Kelly and Hart 1970;Spalt 1958).For modified wood, the first sorption isotherm appears to be due to Mo ¨rath (1931) who steamed beech at 150 and 190 °C and used sulphuric acid solutions for conditioning the material at different RH levels.The second reported sorption isotherm for modified wood seems to be the absorption and desorption isotherms reported by Spalt (1958) for acetylated (R mod = 0.32) white pine, although studies from the intervening period of time mention reduced hygroscopicities of modified wood (Buro 1954;Seborg et al. 1953) without showing sorption isotherms.
Of the different methods for generating specific RH levels detailed above, one of the most common methods is by saturated salt solutions (Ammer 1963b;Eriksson and Noren 1965;Forster 1998;Higgins 1957;Nearn 1955;Papadopoulos and Hill 2003;Wangaard and Granados 1967;Yasuda et al. 1995) since they are relatively simple to use and cover a wide range of RH levels (see, for example, Greenspan 1977) for solutions in the range 3-98% RH.Their use limits, however, which RH levels can be achieved.This problem is avoided when RH is generated by mixing of wet and dry air streams to produce any RH level desired in the range 0-97% RH.Such method has been employed in climate chambers (Strømdahl 2000;Thygesen et al. 2010) and has even be used below the freezing temperature of water (Hedlin 1968).
Whether using early-stage sorption balances, saturated salt solutions, or climate chambers for conditioning wooden samples, it is important to select appropriate sample masses taking into consideration the balance uncertainty and the stability of the conditioning climate in order to minimise the uncertainty in determined MC.The ASTM C1498 (ASTM 2016) and ISO 12571 (ISO 2013) standards which concern the determination of hygroscopic moisture sorption isotherms require the balance accuracy to be at least 0.01-0.1% of the sample mass.However, these standards also require a sample mass of 10 g or more, which is often not feasible or wanted in scientific studies, for example on sorption kinetics.To determine EMC, it is of course important that equilibrium has been obtained as marked by a constant sample mass after conditioning at constant RH for a given time.Many studies employ stability criteria adopted from the ASTM C1498 (ASTM 2016) and ISO 12571 (ISO 2013) standards to specify when the variation in sample mass is low enough for equilibrium to be assumed obtained.These standards require mass variations of \ 0.1% over a period of 5 days after daily mass readings or over 7 days in two consecutive mass readings, respectively.Other stability criteria have been proposed in the literature, for example the one by Wadso ¨et al. (2004) which also takes into account the change in RH from the previous to the current climate.
One important thing to note about the use of early-stage sorption balances, saturated salt solutions, and climate chambers is that these methods are unsuitable for conditioning wood samples to saturation.Nonetheless, several studies condition wood samples in desiccators over liquid water as a way of saturating the material (Boonstra et al. 2007;Feist and Sell 1987;Harju et al. 2002;Kamdem et al. 2002;Meyer and Brischke 2015;Williams and Hale 2003), but the resultant EMC after such conditioning is questionable.This is due to the great difficulty with obtaining saturated water vapour as discussed by Strømdahl (2000).Differences in temperature between vapour and surrounding surfaces will result in either water condensation on surfaces if these are cooler or a non-saturated vapour if the surfaces are warmer.Even slight temperature fluctuations can give noticeable effects.Thus, fluctuations of 0.5 °C around room temperature cause a ''saturated RH'' around 97%, while even small changes of 0.2 °C limit the maximum RH to 99%.As shown by Hoffmeyer et al. (2011), prolonged conditioning of untreated wood in such presumably saturated water vapour will not give the same isotherm upon subsequent desorption as when the material is desorbing from the fully watersaturated state.Moreover, as discussed later under ''Fibre saturation techniques'', it is highly uncertain that conditioning wood in presumably saturated water vapour will actually saturate cell walls.

Automated sorption balances
This type of equipment, often referred to as dynamic vapour sorption (DVS), combines automated climate control in the range 0-97% RH with continuous highresolution mass measurements.RH is generated by mixing wet and dry gas streams, for example air or nitrogen.As only a small sample mass (typically in the 5-20 mg range) is needed like in the quartz helix spring systems, a faster approach to equilibrium is observed compared with conditioning larger samples over saturated salt solutions or in climate chambers.As automated sorption balances are becoming a common research tool in wood research laboratories, the literature is abundant with sorption isotherms determined with such equipment, also for various types of modified wood (Hill et al. 2013;Jalaludin et al. 2010a, b;Popescu and Hill 2013;Xie et al. 2010Xie et al. , 2011)).While the low sample mass makes possible detailed studies of differences in water sorption between earlywood and latewood (Fredriksson and Thygesen 2017), even at the growth ring level (Hill et al. 2015;Song et al. 2014), it also requires the user to consider how to collect representative samples for the study in question.For generating representative samples of larger wood volumes, several studies have used milled wood (Hill et al. 2009;Himmel and Mai 2015;Xie et al. 2010;Zaihan et al. 2009Zaihan et al. , 2011)), but it could be speculated if this gives a representative material or might produce unwanted side effects, for example a change in sample chemistry.Milling is not as intense a method for subdividing wood as ball milling which has been shown to change the material chemistry (Mao et al. 2006) and internal chemical bonding (Schwanninger et al. 2004).Nonetheless, different particle size fractions generated from milling might not have identical chemical composition if some anatomical fractions (e.g.lignin-rich middle lamella fragments) are over-represented in smaller size fractions that are more easily lost upon handling the milled material.Instead of milling, several studies have used wood slices cut with a microtome or razor blade to produce a large sample surface area while preserving the microscopic wood structure (Fredriksson and Thygesen 2017;Fredriksson et al. 2010;Glass et al. 2017;Hill et al. 2015;Hosseinpourpia et al. 2016;Song et al. 2014), while a few use small wood blocks (Fredriksson and Johansson 2016).All of these studies use the automated humidity control by which RH is automatically changed upon reaching a specified mass stability (dm/dt) criterion.If such dm/dt criterion is not selected carefully, however, it can lead to inaccurate results.If the stability criteria of the ASTM C1498 (ASTM 2016) and ISO 12571 (ISO 2013) standards are converted into %/min, i.e. per cent mass change per minute, employed in most automated sorption balances, the dm/dt is around 0.00001%/min.While the stability criteria of the standards are based on mass change relative to the previously observed mass, the dm/dt of automated sorption balances is typically based on mass change relative to the dry sample mass which is often determined by initial drying.Thus, it would seem that the latter gives stricter stability requirements.However, the dm/dt suggested by the producer of one of the most common types of automated sorption balances is 0.002%/min for 10-20 mg samples and 0.0005%/min for 100 mg samples, i.e. about 50-200 times higher than required by the standards for larger ([ 10 g) samples.The other important parameter for establishing if equilibrium is obtained is the ''stability window'', i.e. the time period over which the dm/dt needs to meet the stability criterion.In the ASTM C1498 (ASTM 2016) and ISO 12571 (ISO 2013) standards, the mass change is calculated over 5 and 7 days, respectively, whereas the default stability window in the most popular automated sorption balance in wood research is 15 min.The 0.002%/min over 15-min stability criterion is commonly employed for wood materials, but Glass et al. (2017) recently demonstrated that it gives too low EMC values in absorption and too high in desorption.The error was up to 0.7% MC for southern pine wood and up to 2.9% MC for isolated lignin for changes in the high-RH regime.Moreover, Glass et al. (2017) were unable to find literature data to support the claim that a dm/dt of 0.002%/min is sufficient for accurate determinations of EMC, even in publications often cited to document this claim.The use of a dm/dt criterion becomes critical when considering sample drying which in many studies is performed initially.The results of Christensen and Hergt (1969) show that the longer the time under constant climate before a given change in RH, the slower the rate of sorption.Thus, if a sample is initially dried out using a dm/dt criterion after experiencing stable conditions for a long while, it will exhibit a slower rate of desorption than in a subsequent second drying stage after re-wetting.This will affect the determined dry mass in the two drying stages as the dm/dt criterion is met at higher residual MC in the first stage, thus causing an overestimation of the dry mass.Examining the mass curves in the literature when a dm/dt criterion of 0.002%/min is used, it is often clearly seen that the sample material is far from having reached equilibrium with 0% RH before subsequent rewetting is initiated.The overestimated dry mass is problematic as it shifts the entire sorption isotherm downwards, partly due to the increased dry mass but mainly due to decreased estimated water masses.Moreover, a second drying stage will yield a lower dry mass, perhaps leading to speculations about ''trapped water'' in the first drying stage (Hill et al. 2015).It is recognised that Hergt and Christensen (1965) did find that up to 1% MC is retained in the wood after vacuum-drying if the material is not rapidly dried from the water-saturated state, i.e. conditioning to any level of RH before drying results in incomplete drying even at a drying temperature of 65 °C.However, their experiments used thin, microtomed wood sections and drying times of 1-5 days which is considerably longer than the drying duration if using a dm/dt criterion of 0.002%/min.In the authors' experience, drying for 6 h at 60 °C in dry nitrogen gas of a 20-mg wood block sample gives a reproducible dry mass within 0.05% even when several cycles of re-wetting to various RH (10-90%) and drying are used.For samples cut into thin slices, this drying protocol (6 h at 60 °C, 0% RH) should be more than enough to secure complete drying.Therefore, it is recommended that longer, fixed drying times, preferably at slightly elevated temperatures, are used in automated sorption balances in order to secure accurate dry mass determinations.For getting closer to the ''true'' EMC upon conditioning at a specific RH, it might be useful to lower the dm/dt as well as increase the stability window.However, further investigations are needed for establishing guidelines for the appropriate dm/dt settings for wood.Until such have been established, a degree of uncertainty exists as to the accuracy of reported EMC and sorption isotherms by the use of automated sorption balances.
The difference in MC in equilibrium with a given RH depends on whether equilibrium is reached by absorption or desorption (Urquhart 1929).This phenomenon is termed sorption hysteresis and can be seen as a difference between the absorption isotherm from dry state and desorption isotherm from water-saturated state (see Fig. 2).It has been speculated that it is linked to non-elastic mechanical responses during swelling, but the exact mechanism is not yet understood (Engelund et al. 2013).Nonetheless, wood modification can have different effects on sorption hysteresis.For instance, acetylation decreases the absolute hysteresis (Popescu et al. 2014), i.e.MC difference between the absorption and desorption isotherms, but this result appears to be correlated with the general reduction in MC at given RH.A similar general reduction in MC is found for thermal modification of wood; however, this treatment appears to increase the absolute sorption hysteresis (Olek et al. 2013).Automated sorption balances are useful for studies of sorption hysteresis, and several published studies have focused on this phenomenon.In most studies, absolute sorption hysteresis has a maximum around 70-75% RH and decreases for higher RH (Hill et al. 2009(Hill et al. , 2012b;;Himmel and Mai 2015;Hosseinpourpia et al. 2016;Popescu et al. 2014;Simon et al. 2017;Xie et al. 2011).The decreasing sorption hysteresis at high RH might, however, be an artefact due to desorption being initiated from a non-saturated state (often by conditioning to 95% RH).Thus, in recent results by Fredriksson and co-workers (Fredriksson and Johansson 2016;Fredriksson and Thygesen 2017), sorption hysteresis in untreated Norway spruce gradually increases up to at least 95% RH (see Fig. 2).Furthermore, as shown by Hoffmeyer et al. (2011), desorption from a non-saturated state will yield a lower desorption isotherm than for water-saturated wood until around 75% RH.Therefore, most desorption isotherms using automated sorption balances reported in the literature can be regarded as scanning curves above 70-75% RH Based on data from Fredriksson and Johansson (2016), Fredriksson et al. (2013) at 20 °C and from authors at 25 °C (dm/dt = 0.0005%/min).Note that the desorption and absorption isotherms do not coalesce at 97% RH; b absolute sorption hysteresis of untreated wood from moisture content data acquired with automated sorption balances.Based on data for Norway spruce (Fredriksson and Johansson 2016;Fredriksson and Thygesen 2017;Ka ¨llbom 2015;Sedighi-Gilani and Schwarze 2015) and Sitka spruce (Hill et al. 2009).Note that in all studies except (Fredriksson and Johansson 2016;Fredriksson and Thygesen 2017) desorption is initiated from a non-saturated state (conditioning to 95% RH) produced by a reversal in sorption direction (absorption/desorption) at intermediate RH, i.e. away from the dry or saturated states.
Besides deriving MC in equilibrium with various climates, automated sorption balances have been used for studying the kinetics of sorption.In many polymeric materials, the sorption process exhibits two-stage behaviour (Christensen and Kelsey 1959;Crank 1953;Kelly and Hart 1970;Newns 1956) and cannot be explained by Fickian diffusion into the material.Thus, sorption is often a much slower process than expected for diffusion through the solid material which is also the case for wood (Christensen 1959;Kelly and Hart 1970).Instead, sorption might be controlled by relaxation of stresses arising from swelling of the dense and stiff wood cell walls (Christensen 1959(Christensen , 1960;;Christensen and Kelsey 1959).This also explains why samples conditioned for longer times after a change in RH exhibit slower rate of sorption as observed by Christensen and Hergt (1969).Due to the two-stage nature of the sorption process, it has often been analysed by the parallel exponential kinetics (PEK) model featuring two moisture sorption processes (fast and slow exponentials) using data for both untreated (Hill et al. 2010) and modified wood (Hill et al. 2012a;Himmel and Mai 2016;Jalaludin et al. 2010a;Popescu et al. 2014;Xie et al. 2010).It may be difficult to compare the results across studies due to differences in sample mass and geometry as diffusion into the sample at some point becomes relevant, at least in the low-RH regime (Christensen 1965).More importantly, however, the recent results by Glass et al. (2017) demonstrate that the PEK model is not capable of capturing the long-term kinetics of untreated wood if the data acquisition at a given RH is disrupted when a dm/dt stability criterion of 0.002%/min is met.This seriously questions whether results derived using PEK model fittings on sorption data with a 0.002%/min dm/dt criterion are correct.Adding to this uncertainty, the PEK model is quite sensitive to initial guesses of parameter values of the two moisture sorption processes (Himmel and Mai 2016).
More details about wood-water interactions can be gained from automated sorption balances using heavy water, i.e. deuterium oxide ( 2 H 2 O or D 2 O), instead of normal water ( 1 H 2 O or H 2 O) as conditioning vapour.Hereby, hydrogen of the hydroxyl (OH) groups interacting with heavy water will be exchanged for deuterium if the material is continuously exposed to 2 H 2 O vapour for sufficient time.As deuterium ( 2 H) is about 1 g/mol heavier than protium ( 1 H), the material dry mass increases slightly upon deuteration (* 1% for untreated wood), which can only be measured with enough accuracy using very sensitive balances.Therefore, vapourphase deuteration has been used on untreated wood in an early-stage sorption balance (Taniguchi et al. 1978) as well as recently on untreated (Thybring et al. 2017) and modified wood (Popescu et al. 2014;Rautkari et al. 2013) using modern automated sorption balances, yielding information about how many hydroxyls are accessible for water.It should be noted, however, that not all accessible hydroxyls will have their hydrogen exchanged.For instance, Lindh et al. (2016) found that only two out of three surface hydroxyls on cellulose microfibrils are accessible for deuteration as O(3)H hydroxyls act solely as hydrogen bond acceptors, i.e. interacting exclusively with water via the O(3) (Lindh et al. 2016).Therefore, out of around 6 mmol/g accessible hydroxyls in wood cellulose (Engelund et al. 2013), only 4 mmol/g will contribute to the increased dry mass after deuteration.In order to determine the hydroxyl accessibility, the sample has to be dried initially and after deuteration for determining the change in dry mass.Hereafter, the hydroxyl accessibility, OH access [mol/g], can be calculated as where Dm dry [g] is the change in dry mass due to deuteration, m 0,dry [g] is initial dry mass, DM hydrogen [g/mol] is molar mass difference between deuterium and protium (1.006 g/mol), and (1 ?R mod ) -1 [-] corrects for mass gain due to modification if present.The change in dry mass is dependent on not only the added deuterium on hydroxyls, but also the residual normal and heavy water in the first and last drying stages, respectively, i.e.
where Dm dry,true [g] is the ''true'' dry mass increase due to hydroxyl deuteration, m 1,water and m 2,water [g] are residual water masses after first and last drying stages, respectively, c [mol/g] is molar concentration of residual water per gram of wood, M [g/mol] is molar mass, and indices ''D 2 O'' and ''H 2 O'' refer to heavy water and normal water, respectively.Even if the concentrations are similar after both drying stages, the dry mass is still increased due to deuterium associated with the heavy water.Thus, a residual water content corresponding to 1% MC of normal water will add 0.11% mass.Compared with an expected dry mass increase due to deuteration in the range of 0.5-1.0%,this will yield an error of 11-22%.If only 0.1% MC residual water is left after drying, the error is one order of magnitude lower.However, this is only for similar concentrations of normal and heavy water after drying, which is an unlikely situation, even after similar drying conditions and durations in both first and last drying stages, since transport of heavy water is about 20% slower than that of normal water (Mills 1973).Therefore, if assuming a conservative estimate of 50% higher concentration of heavy water than of normal water after the drying stage, the measurement errors in determined hydroxyl accessibility for residual water contents of 1 and 0.1% MC would be 66-133 and 7-13%, respectively.This clearly demonstrates the need for sufficient drying of the material to eliminate residual water as a significant source of error.Therefore, the use of a dm/dt stability criterion for determining hydroxyl accessibility is not recommended due to the uncertainties it involves as discussed previously for generating sorption isotherm data.As mass differences involved in the determination of hydroxyl accessibility are much lower than those for sorption isotherms, the potential errors introduced are much higher.
The following is a short guideline for using heavy water in automated sorption balances.In general, it is advisable to use at least three but preferably more replicates for each sample batch due to measurement uncertainties.Only one conditioning stage in D 2 O vapour is needed as the determined accessibility depends on time of exposure and not on number of repeated conditioning cycles (Po ¨nni et al. 2013).Moreover, the carrier gas (often nitrogen) should be dried out before the gas

Pressure plate technique
As previously discussed, it is very difficult to condition samples close to saturation due to problems with keeping the temperature stable enough to achieve equilibrium in the over-hygroscopic RH range (98-100% RH).Therefore, in this range the pressure plate technique is often employed for conditioning, a method in which one or several samples are placed within an extractor where pressures up to 100 bar (10 MPa) can be applied.Typically, samples are cuboids or cylindrical discs with side length/diameter in the range 10-60 mm and a thickness of 2-20 mm (Almeida andHernandez 2006, 2007;Cloutier and Fortin 1991;Fredriksson and Johansson 2016;Fredriksson et al. 2013;Thygesen et al. 2010;Zauer et al. 2016); however, Stone and Scallan (1967) used the pressure plate technique to condition microtome sections of 100 lm thickness and a cross-sectional area of 2 cm 2 .After watersaturating samples and the supporting substrate (ceramic plate or cellulose membrane, dependent on pressure range), pressure is applied to the extractor using dry gas, whereby water is pressed out of pores in which capillary forces are lower than the applied pressure.From a combination of the Young-Laplace and Kelvin equations, a correlation between applied pressure in the extractor and RH can be found: where DP [Pa] is applied pressure (above atmospheric pressure), M w [kg/mol] is molar mass of water (0.018 kg/mol), R [J/mol K] is the universal gas constant (8.3145J/mol K), T [K] is temperature, and q w [kg/m 3 ] is density of water.By adjusting the applied pressure, it is possible to tune RH much more precisely in the 98-100% RH range than with other methods.For instance, an applied pressure of 15 bar corresponds to 98.90% RH, while 9 bar is equivalent to 99.34% RH.When equilibrium is reached as monitored on the amount of water pressed out over time, samples are taken out and weighed on a balance of sufficient mass resolution.
According to standard ISO 11274 (ISO 1998) equilibrium is obtained when two consecutive weekly mass readings show \ 0.1% variation in sample mass.The technique has primarily been used to condition untreated wood (Almeida andHernandez 2006, 2007;Cloutier and Fortin 1991;Fortin 1979;Fredriksson and Johansson 2016;Fredriksson et al. 2013;Griffin 1977;Stone and Scallan 1967), but some studies of modified wood have been published (Thygesen et al. 2010;Zauer et al. 2016) (see Fig. 3).The traditional use of the technique has been for conditioning in desorption from water-saturated state; however, Fredriksson and Johansson (2016) used a rebuilt extractor for absorption measurements which avoids problems of previous modifications of the technique to handle absorption (Cloutier and Fortin 1991;Fortin 1979).
EMC in the over-hygroscopic RH range depends on the size of voids in the wood structure such as pits, lumina, vessels, etc., and the sorption isotherms therefore depend heavily on wood species due to the high variability of these anatomical characteristics.Moreover, while the over-hygroscopic absorption isotherms depend on water condensation in voids and pores of gradually increasing void sizes (Engelund et al. 2010), desorption isotherms are governed by the size of bottleneck pits connecting lumina/vessels to the surrounding climate (Fredriksson and Johansson 2016;Salin 2008).During absorption, lumina/vessels take up moisture when an RH corresponding to equilibrium with their lumen/vessel radius is reached, but during desorption they will empty at a lower RH corresponding to radius of their largest pits facing the surroundings.This is known as the ''ink bottle effect'' (Ravikovitch and Neimark 2002) and might explain the large sorption hysteresis in the over-hygroscopic regime.In an interesting study by Zelinka et al. (2016), mercury intrusion porosimetry is explored as a complementary technique for studying over-hygroscopic moisture sorption in wood.The technique is based on applying pressure to push liquid mercury into thinly cut transverse wood sections.The applied pressure correlates with the radius of voids filled with mercury as described by the Kelvin equation.Therefore, after adjusting for differences in contact angle and surface tension between mercury and liquid water, the volumetric Fig. 3 a Desorption isotherms in the over-hygroscopic range determined after conditioning with the pressure plate technique and saturated salt solutions for untreated and modified Norway spruce (Thygesen et al. 2010;Zauer et al. 2016).Data points in the range of 94-98% RH from climate chamber measurements have been omitted from the results of Thygesen et al. (2010) for clarity.b Zoom of (a) in the 99.5-100.0%RH range Wood Sci Technol (2018) 52:297-329 uptake of mercury as a function of applied pressure can be used to calculate the capillary water in the wood void structure as a function of RH.The order of mercury intrusion into the wood structure, i.e. vessels-lumina-pits, makes the calculated capillary water curves correspond to the desorption isotherms of thin wood sections.For thicker wood sections without direct access to cell lumina, the order of emptying of voids as well as filling of voids with mercury changes.Mercury intrusion porosimetry would therefore need to be performed on such thicker samples in order to correctly represent over-hygroscopic desorption isotherms of capillary water.
Although wood modifications rarely cause significant changes to the void geometry unless modification agent fills up lumina or pits, they may considerably affect the over-hygroscopic sorption isotherms.This is due to potential changes after modification in the contact angle between wood surface and wetting liquid.Such change will affect EMC in the over-hygroscopic range as predicted by Engelund et al. (2010).The contact angle of liquids with modified wood has been investigated in several studies by the Wilhelmy and sessile drop methods.The first determines the force required during immersion and withdrawal from a probe liquid, for example water, while in the latter the contact angle is determined over time after dropping a liquid droplet onto the surface.Despite the shortcomings of both methods as discussed by Wa ˚linder and co-workers (Wa ˚linder and Johansson 2001; Wa ˚linder and Stro ¨m 2001), they both show that modification can alter the contact angle with water.For instance, both thermal modification (Bakar et al. 2013;Hakkou et al. 2006;Kocaefe et al. 2008) and acetylation (Bryne and Wa ˚linder 2010) appear to increase the contact angle with water, while furfurylation (Bryne and Wa ˚linder 2010) slightly decreases it.The effect of these changes on the overhygroscopic sorption isotherms has not been investigated.

Fibre saturation techniques
The original definition of the fibre saturation point (FSP) was suggested by Tiemann (1906) as the MC reached upon drying where lumina are empty of liquid water, cell walls begin to dry, and the mechanical properties start increasing.As shown by Stamm (1971), however, these three phenomena do not occur at the same MC.FSP was later suggested to be determined as MC at the first shrinkage upon drying from green or water-saturated states (Koehler and Thelen 1926).However, due to inconsistencies in the previous definitions, Stone and Scallan (1967) focused on the emptying of lumina in the over-hygroscopic RH range and defined FSP as the MC where lumina are empty, but cell walls are still water-saturated.Using the pressure plate technique, they found FSP for black spruce at 99.75% RH, while Griffin (1977) added his own experimental results and adjusted it to 99.93% RH.This definition is also problematic as liquid water may still be present in some cell lumina even at equilibrium conditions where cell walls in other parts of the same sample are not saturated (Passarini et al. 2015).
Although the traditional FSP concept originated more than a century ago, it is still widely used to describe the transition point imagined between cell wall saturation and (significant) moisture uptake in the wood void structure in the overhygroscopic regime.For practical purposes, for example for structural engineering applications of wood with a strong focus on mechanical performance, the traditional concept of FSP works well.However, from a scientific standpoint the concept is problematic (Babiak and Kudela 1995;Feist and Tarkow 1967), and it mischaracterises the wood-water relations.Instead, in this review, FSP is regarded as the maximum cell wall water capacity measured in the fully water-saturated state (Hill 2008).The first part of this definition is in line with the intension behind the original definition by Tiemann (1906), i.e. the amount of water found within cell walls when these are water-saturated.By determining FSP in the fully water-saturated state, i.e.where water takes up all available space in cell walls and voids (pits, lumina, vessels) in the wood structure, uncertainties about partial drying of any cell walls are avoided.However, investigating FSP in the fully water-saturated state requires experimental techniques that can differentiate between water within and outside cell walls.

Differential scanning calorimetry
One technique which is capable of separating cell wall MC from the total MC is differential scanning calorimetry (DSC).This type of equipment measures the energy involved in changing the temperature of a small sample over a wide temperature range.The amount of sample depends on the geometry of the crucible of the specific DSC instrument, but typical sample masses and sizes are in the range 5-15 mg (Park et al. 2006;Weise et al. 1996) and 1-mm discs of 3-4 mm diameter (Simpson and Barton 1991;Zauer et al. 2014), respectively.The technique is particularly useful for determining the energy involved in exothermic or endothermic reactions occurring during an experiment, for example phase transitions such as freezing or thawing of water.As cell wall water is tightly associated with the chemical constituents, it will not exhibit a phase change at least down to -70 °C (Berthold et al. 1994).This was first used by Magne et al. (1947) to determine the amount of cell wall water in water-saturated textile fibres (cotton, rayon, nylon, glass fibres) by measuring the energy for freezing water at -4 °C.Since then, DSC techniques have been used to determine cell wall water within, for example, pulp and paper fibres (Maloney 2000;Nelson 1977;Park et al. 2006;Weise et al. 1996), untreated wood (Simpson and Barton 1991;Zauer et al. 2014;Zelinka et al. 2012), and modified wood (Repellin and Guyonnet 2005;Zauer et al. 2014).By cooling or quenching a wet sample to a given sub-zero temperature and ramping the temperature at constant rate to a specified temperature above the freezing point of water, the energy input (enthalpy) required to melt the ice formed can be measured.Temperature ranges and ramping rates employed vary between studies, but are most often found within the -40 to 40 °C range (Park et al. 2006;Repellin and Guyonnet 2005;Simpson and Barton 1991;Weise et al. 1996;Zauer et al. 2014) and from 1 to 3 °C/min (Park et al. 2006;Repellin and Guyonnet 2005;Weise et al. 1996;Zauer et al. 2014) to 5 °C/min (Zelinka et al. 2012) or even 10 °C/min (Nelson 1977;Simpson and Barton 1991), respectively.The cell wall moisture content, MC cell-wall [%], can then be found as where m water [g] is total mass of water, m tot [g] is total wet sample mass, m 0,dry [g] is dry sample mass, DH tot [J/g] is measured specific enthalpy of melting per total wet sample mass, H f [J/g] is melting enthalpy (often termed enthalpy of fusion) of pure water (333.6J/g) (Lide 2013), and (1 ?R mod ) -1 [-] corrects for mass gain due to modification if present.If the sample is fully water-saturated, MC cell-wall corresponds to FSP, and using such samples avoids problems of non-saturation of any parts of the wood volume examined.The technique can, however, also be used to differentiate cell wall water from liquid water outside cell walls in the over-hygroscopic RH range as long as a sufficient amount of freezable water is present.
Although the DSC technique removes the sample from its conditioning temperature and FSP has been reported to change with temperature (Babiak and Kudela 1995;Stamm and Loughborough 1935), it would be unreasonable to relate the determined FSP to a temperature around 0 °C since the experimental run in the DSC is much shorter than the time to equilibrate MC.Therefore, it can be assumed that the determined FSP relates to the conditioning temperature.Besides the above-described method for determining FSP, an additional one has been reported in several studies (Simpson and Barton 1991;Zauer et al. 2014;Zelinka et al. 2012) in which samples within a wide range of MC have been used and a linear regression is constructed between MC and specific enthalpy measured.The inflection point of this linear regression is then FSP, and the slope of the line is the enthalpy of water associated with the wood structure (Hager and Macrury 1980).In some studies (Simpson and Barton 1991;Zauer et al. 2014;Zelinka et al. 2012), it is stated that this method is more correct as it does not assume any specific value of the melting enthalpy.However, all studies show that the melting enthalpies determined are close to the generally accepted value of 333.6 J/g (Magne et al. 1947;Nelson 1977;Zauer et al. 2014) as shown in Fig. 4, considering the experimental uncertainties (small sample masses, high/low temperature ramping rates, heat flow accuracy, etc.) involved in the measurements.The only exception to this is the study by Zelinka et al. (2012) in which a very low melting enthalpy of 240 J/g is found, leading the Fig. 4 Melting enthalpy as a function of reduced moisture content (corrected for added mass of modification) for untreated and variously modified Norway spruce samples.Some of the samples had excess water on surfaces, resulting in very high calculated moisture contents, even though this water was not physically inside the material.Authors' own data authors to speculate about temperature buffering from the sample heat capacity.However, this effect should be mostly caught by baseline correction of the recorded energy flow.That the melting enthalpy can differ from that of bulk water is clear for water found in nano-pores (Ja ¨hnert et al. 2008;Shimizu et al. 2015).A vast majority of liquid water in wood is, however, found in voids larger than 0.1 lm.This is, for instance, seen in the negligible capillary water uptake at 98% RH (Engelund et al. 2010) corresponding to a void size of 0.1 lm.For voids larger than this size, the melting enthalpy reduction due to confinement is \ 1% (Ja ¨hnert et al. 2008).If freezable water was found in wood nano-pores in detectable amounts, it would exhibit a significant freezing/melting point depression (Shimizu et al. 2015); however, no peak of appreciable lower melting point has been found for solid wood.
Besides the effect of heavy confinement on melting enthalpy of liquids, there are no results found in the literature to suggest that surface chemistry or other factors should affect the melting enthalpy.Therefore, determining FSP by the DSC method using fully water-saturated samples is deemed appropriate, whereas conditioning to lower MC might not result in fully saturated cell walls and as a result yield a lower determined FSP.

Solute exclusion technique
The amount of different pools of water in very small voids can be investigated by the solute exclusion technique (SET).The method was originally developed by Aggebrandt and Samuelson (1964) and modified by Stone and co-workers (Stone et al. 1969;Stone andScallan 1967, 1968a, b) to characterise pore volumes and pore size distributions in pulps in the water-swollen state.Since then, SET has been used to characterise both untreated (Ahlgren et al. 1972;Farahani 2003;Feist and Tarkow 1967;Flournoy et al. 1991Flournoy et al. , 1993;;Forster 1998) and modified wood (Farahani 2003;Forster 1998).In SET, a water-saturated sample with known mass of water is placed in an aqueous solution of molecular probe solutes of a known size and concentration.The typical sample is a cuboid of 20-30 mm side length and a thickness of 3-5 mm (Farahani 2003;Flournoy et al. 1991Flournoy et al. , 1993;;Forster 1998), although several older studies have used shavings, pulp, or milled wood (Ahlgren et al. 1972;Feist and Tarkow 1967;Stone and Scallan 1968b).As solutes diffuse into the water within the sample, their concentration in the surrounding solution decreases.The change in concentration is related to the cell wall water volume that can be penetrated by solutes, and from this change in concentration the inaccessible moisture content, MC inaccess [% g/g], can be calculated: where m water [g] is total mass of water in the sample before immersion, m sol [g] is mass of solution added, m dry [g] is dry mass of the sample, c initial [mol/g] is initial concentration of solutes in the solution, Dc [mol/g] is change in solute concentration after immersion of the sample, and (1 ?R mod ) -1 [-] corrects for mass gain due to modification if present.Concentration differences are small, and a HPLC or calibrated refractometer is needed for determining the concentration before and after immersion of the wet sample.To maximise the concentration difference, solution volumes used to immerse samples in are typically not more than 3-10 times the sample volume (Farahani 2003;Flournoy et al. 1991;Forster 1998;Stone and Scallan 1968b).It is important that equilibrium has been achieved before the change in concentration is determined.However, the incubation time employed in the literature varies dramatically from a couple of minutes for pulp samples (Stone and Scallan 1968b) to 1-14 days for solid wood samples (Farahani 2003;Flournoy et al. 1991;Forster 1998).Clearly, more research is needed concerning the approach to equilibrium and appropriate incubation times as a function of diffusion pathways into samples.
When selecting solutes for SET, it is important that they are chemically inert towards the sample and do not adsorb to surfaces within cell walls (Stone and Scallan 1967) as this would shift the equilibrium between pore solution and surrounding solution.The solutes most often used are small sugars, PEG, and dextrans.Their size is determined from the hydrodynamic diameter calculated from the Stokes-Einstein equation based on their diffusivity.The smallest solutes are saccharides such as glucose (0.8 nm) (Flournoy et al. 1991;Kuga 1981;Lin et al. 1987;Stone and Scallan 1968a;Wang et al. 2011), mannitol (0.8 nm) (Carpita et al. 1979), sucrose (1.0 nm) (Carpita et al. 1979;Kuga 1981;Van Dyke 1972), cellobiose (1.0 nm) (Lin et al. 1987), maltose (1.0 nm) (Flournoy et al. 1991;Stone and Scallan 1968a), raffinose (1.2 nm) (Flournoy et al. 1991;Kuga 1981;Stone and Scallan 1968a;Van Dyke 1972), and stachyose (1.4 nm) (Carpita et al. 1979;Stone and Scallan 1968a).PEGs range in size from 0.44 nm for ethylene glycol (Kuga 1981) to 13.0 nm (Lin et al. 1987;Van Dyke 1972), while dextrans have been reported in the literature from 1.8 nm in diameter (Ishizawa et al. 2007) up to 586 nm (Arond and Frank 1954). Figure 5 shows how the solute diameter correlates with molecular mass for both PEGs and dextrans.Arond and Frank 1954;Carpita et al. 1979;Day et al. 1978;Flournoy et al. 1991;Granath and Kvist 1967;Hill et al. 2005;Ishizawa et al. 2007;Kuga 1981;Lin et al. 1987;Peters 1986;Stone et al. 1969;Stone and Scallan 1968a;Van Dyke 1972;Wang et al. 2011;Wong et al. 1988).b Inaccessible cell wall water volume as a function of size of different molecular probes.Data for untreated and acetylated (R mod = 0.21) Corsican pine from Forster (1998).Note that the fibre saturation point (FSP) can be determined from the plateau in inaccessible cell wall volume for probe sizes too large to enter cell walls The inaccessible MC describes the mass of water inaccessible for solutes of a given size.Thus, by using several differently sized solutes it is in principle possible to determine the pore size distribution within cell walls as indicated in Fig. 5.However, the underlying assumption of (6) which allows calculation of probed water mass is that solute concentration within accessible pores is the same as in the surrounding solution (Lin et al. 1985).This assumption has been shown to be incorrect, since solute concentration within accessible pores decreases as the sizes of probe solutes and accessible pores get closer (Alince 1991;Day et al. 1978Day et al. , 1979)).Thus, for solutes about half as big as the accessible pore, the solute concentration within the pore is around 30% of that in the surrounding solution (Alince 1991), while it decreases to nearly zero for solute sizes close to the pore size.For this reason, SET does not provide absolute pore size distributions, but can only hint at it.However, for solutes too large to enter cell walls, the voids they occupy (pits, lumina, vessels, etc.) are much larger than the solute size.Therefore, the assumption of equal solute concentration in these voids to the ones in the surrounding solution is valid, and SET can offer a good estimate of total pore volume, i.e.FSP (Alince 1991).If several sizes of solutes are used which are too big to enter cell walls, a plateau corresponding to FSP is found for MC inaccess as a function of solute size (see Fig. 5).

Spectroscopic techniques
While the gravimetric and fibre saturation techniques are useful for determining total amount of water in wood or amount of cell wall water in the over-hygroscopic range, they are not able to provide a deeper insight into wood-water interactions or the distribution of water within the wood structure in the entire RH range.For this, spectroscopic techniques are needed that can distinguish water in different environments or provide more detail about molecular interactions and chemistry within cell walls.A wide range of spectroscopic techniques are available, but here two of the most important ones are highlighted.

Low-field nuclear magnetic resonance spectroscopy and magnetic resonance imaging
Low-field (LF) nuclear magnetic resonance (NMR) spectroscopy, also known as time-domain NMR (Thygesen and Elder 2008) or NMR relaxometry (Kekkonen et al. 2014), can distinguish between water-bound hydrogen located in different physical and chemical environments.Consequently, it is a good technique for determining the distribution of water within different parts of the wood structure (Araujo et al. 1992).The technique uses the very, very slight magnetic polarisation in a material caused by a permanent magnet within the spectrometer (typically up to a few tesla in strength).The polarisation arises because the permanent magnetic field B 0 gives rise to two energy states for the spins (along and against the field) of hydrogen nuclei within the sample with a tiny difference in abundance given by the Boltzmann distribution.This difference is what causes a polarisation vector M of the material aligned with B 0 .The spectrometer is also equipped with an electromagnetic coil which can generate short-lived radiofrequency pulses that result in a magnetic field B 1 () B 0 ) oriented perpendicular to B 0 .These pulses will change the orientation of M, and the spectrometer will afterwards record the decay or relaxation of M back to its equilibrium value, realigned along B 0 .This relaxation is characterised by two different relaxation times termed T 1 (longitudinal or spinlattice relaxation) and T 2 (transverse or spin-spin relaxation).T 1 relaxation is studied as the increase in M in the direction along the permanent field B 0 , while T 2 relaxation is studied in the transverse plane, i.e. perpendicular to B 0 , as the decrease in the component of M in this plane.This latter process does not only depend on the realignment of M along B 0 , but also on the dephasing of spins in the transverse plane due to transfer of spins between neighbouring nuclei.T 1 and T 2 relaxations are both first-order rate processes and can therefore be described by exponential decay functions with characteristic relaxation times T 1 and T 2 , respectively.For fluids in porous materials, the relaxation times depend on the size of pores confining the fluid as well as interactions between pore wall and fluid (Mitchell et al. 2005).If a material contains fluid in micro-and macro-pores of various sizes, the fluid will exhibit several relaxation times characteristic of pore sizes and chemistry.By measuring the decay of the low-field (LF) NMR signal and deconvoluting it, the relaxation times of the different components, i.e. fluid in variously sized pores, can be determined and the distribution of fluid in different environments (pore sizes and chemistry) can be analysed.The deconvolution can be done by fitting discrete exponential functions to the decaying signal, for example to add components until residuals are randomly distributed.Another way is to fit continuous relaxation time distributions to the signal by non-negative least square fitting (Whittall et al. 1991) using algorithms often found in computer software for LFNMR data analysis, for example the CONTIN algorithm by Provencher (1982a, b).Before decay of the LFNMR signal can be measured, the spins of the sample need to be perturbed by one or more B 1 pulses.For this, several pulse sequences have been developed.For determining T 1 relaxation times, one of the pulse sequences termed ''inversion recovery'' and ''saturation recovery'' is used, while T 2 relaxation times are determined after free induction decay (FID) or the Carr-Purcell-Meiboom-Gill (CPMG) sequence (Carr and Purcell 1954;Meiboom and Gill 1958).In particular, CPMG is a commonly used sequence for determining T 2 relaxation as it reduces effects of magnet inhomogeneity by multiple refocusing of the spins in the transverse plane.
Several studies have employed LFNMR for characterising water in both untreated (Almeida et al. 2007;Araujo et al. 1992Araujo et al. , 1994;;Cox et al. 2010;Flibotte et al. 1990;Fredriksson and Thygesen 2017;Labbe ´et al. 2002Labbe ´et al. , 2006;;Menon et al. 1987;Passarini et al. 2015;Telkki et al. 2013;Thygesen andElder 2008, 2009) and modified wood (Elder et al. 2006;Hietala et al. 2002;Javed et al. 2015;Kekkonen et al. 2014;Thygesen andElder 2008, 2009).Based on a deconvolution of the decay curves, signals from water in different compartments within the wood structure have been analysed (Almeida et al. 2007;Araujo et al. 1992Araujo et al. , 1993;;Fredriksson and Thygesen 2017;Passarini et al. 2015) (see Fig. 6).A comparison of void geometries and volumes shows good correspondence between pore volumes of pits and lumina and the areas of respective T 2 peaks assigned to these voids for water-saturated samples (Fredriksson and Thygesen 2017).Even more details about the chemical environment of water in pores within wood can be gained from T 1 to T 2 2D LFNMR correlation spectroscopy as reported by Cox et al. (2010).In this method, T 1 and T 2 relaxations are determined simultaneously allowing an analysis of their correlation (Song 2009;Song et al. 2002).This is useful since T 2 relaxation occurs both as a result of realignment of M, like T 1 relaxation, and due to local interactions independent of T 1 relaxation.For a material with a spatially homogenous chemistry, relaxation times depend only on pore sizes which would result in similar T 1 /T 2 ratios for the peaks observed (Song et al. 2002).On the other hand, if surface interactions vary between differently sized pores, this would show up as a markedly different T 1 / T 2 ratio.In the results of Cox et al. (2010), the T 2 peak assigned to water in cell walls by using a regular CPMG splits up into two components in the T 1 T 2 correlation diagram implying two different environments for this water.
There are a number of important issues to remember when using LFNMR to characterise water in wood.First of all, recorded decay curves should all be fully decayed within the duration of a scan, i.e. one decay measurement, in order to retrieve meaningful results from the data analysis algorithm.Secondly, temporal resolution should be adequate for capturing the shortest expected relaxation time.It is recommended that at least five data points are acquired before the shortest expected characteristic relaxation time in order to yield reasonable curve fits.Thus, the shortest expected relaxation time should be divided by five to find the maximum temporal spacing between equidistant pulses in a CPMG sequence.Since water in wood exhibits typical T 2 relaxation times in the range from milliseconds to seconds as shown in Fig. 6, each scan needs to capture a decaying signal over a range of [ 1 s with a resolution of * 1/5 ms.As several scans are needed for adequate signal-to-noise ratio, each wood sample is exposed to a large number of Fig. 6 Low-field NMR characterisation of water in untreated and acetylated (R mod = 0.20) Norway spruce in a fully water-saturated state: a zoom of CPMG decay curves for three replicates of untreated and acetylated spruce, b average T 2 relaxation distribution curves for the three replicates derived from nonnegative least square analysis with regularisation using Prospa version 3.1 (Magritek, NZ) software.Note that upon acetylation, the T 2 relaxation times increase for water in pits and lumina due to weakened wood-water interactions.The cell wall water appears to split into two distinct pools after acetylation.The amount of moisture (MC R ) found in lumina remains the same, while the amounts in pits, small cavities, and cell walls are reduced.Original LFNMR data from Thygesen and Elder (2008) radiofrequency pulses.Although the amount of energy put into the sample by each pulse is tiny, the accumulated number of pulses can heat up the sample considerably in extreme cases.This is problematic as isothermal conditions are required for LFNMR since the magnitude of M is temperature dependent (Dinesh and Rogers 1971).To limit this excessive heating, the waiting time between scans can be increased.
Magnetic resonance (MR) imaging uses the same principle as LFNMR although typically with considerably stronger permanent magnetic (2-14 tesla) and electromagnetic fields.By employing CPMG sequences, relaxation spectra from different parts of the wood structure can be analysed.This has been done on both annual ring level (Dvinskikh et al. 2011;Quick et al. 1990) and cell level (Almeida et al. 2008;Hernandez and Caceres 2010;Meder et al. 2003;Passarini et al. 2015), allowing a visualisation of water distributions in both the hygroscopic (Hernandez and Caceres 2010;Passarini et al. 2015) and over-hygroscopic RH range (Almeida et al. 2008).MR imaging is mainly appropriate for water with longer relaxation times than bound water (Hernandez and Caceres 2010).

Vibrational spectroscopy
Molecular bonds in a material vibrate at certain frequencies dependent on involved atomic masses and stiffness (elastic force constant) of the bond in question.For historic reasons, vibrational frequencies are often given in wavenumbers (cm -1 ) determined as frequency divided by speed of light (3 9 10 8 m/s).As an example, the covalent bonds C=C, C-C, and C-H all have different vibrational wavenumbers due to differences in either bond stiffness (single or double bond) or in masses of the atoms involved.In vibrational spectroscopy, this is used to probe the chemistry of materials as different chemical bonds are seen to vibrate at different wavenumbers.For instance, infrared (IR) spectroscopy detects the absorption of photon (light) energy at different wavelengths from an incoming infrared beam containing all wavelengths in the range 4000-400 cm -1 , while near-infrared (NIR) spectroscopy uses a similar principle, but a different wavenumber range (14.000-4000 cm -1 ) in which mainly vibrational overtones and coupling of vibrations are found.Unlike these techniques, Raman spectroscopy detects the shift in scattered photon energy from a material which occurs for 10 -6 -10 -8 of scattered photons from an incoming laser beam of fixed wavelength.This shift corresponds to a change in vibrational energy of molecular bonds in the material.Thus, Raman spectroscopy probes the same molecular bonds as IR and NIR spectroscopy, but is based on a different physical phenomenon.
Vibrational spectroscopy has been widely used to characterise chemical changes during wood modification (see, for example, Gonza ´lez-Pen ˜a and Hale 2011; Kotilainen et al. 2000;Mitsui et al. 2008;Nuopponen et al. 2005;Stefke et al. 2008;Sun and Sun 2002;Temiz et al. 2007;Tserki et al. 2005;Vena ˚s et al. 2006;and Fig. 7).Moreover, spatial chemical distributions in advanced modified materials have been mapped on the cell wall level using confocal Raman microscopy (Cabane et al. 2014(Cabane et al. , 2016;;Keplinger et al. 2015Keplinger et al. , 2016)).For studying the relation with water, vibrational spectroscopy has been used to examine the presence of water or the sorption of water from the molecular vibrations of water (Christensen et al. 2006;Inagaki et al. 2008;Tsuchikawa and Tsutsumi 1998).One of the advantages of vibrational spectroscopy is, however, the insights it can provide into chemical interactions of wood materials with water.By deuteration of a material, i.e. exchanging protium ( 1 H) for deuterium ( 2 H) in chemical groups hydrogen-bonding with liquid D 2 O or D 2 O vapour (Lindh et al. 2016), a marked shift in vibrational wavenumbers of affected molecular bonds is seen.For instance, deuteration changes the fundamental stretching vibrations of O-H, N-H, and C-H by around 1000, 900, and 800 cm -1 , respectively, which is due to the near-doubling of the vibrating hydrogen mass.If the extra neutron of deuterium instead was added to the heavier vibrating atom of the O-H, N-H, and C-H bonds, for example exchanging 16 O for 17 O, the shift in vibrational wavenumber would hardly be detectable with vibrational spectroscopy as the relative mass increase in O, N, and C would be significantly lower.Compared with gravimetric techniques using deuteration, the advantage of vibrational spectroscopy is that deuteration can be achieved without drying the material prior to exposure to D 2 O, unlike gravimetric techniques which require determination of initial dry mass.
From vibrational spectroscopic studies, it can be seen that wood-water interactions are dominated by interactions with O-H groups (hydroxyls), since deuteration only affects O-H vibrations, while C-H stretching vibrations remain unchanged and no C-D stretching vibrations are observed (Hofstetter et al. 2006;Mann and Marrinan 1956a;Schmidt et al. 2006;Taniguchi et al. 1966;Watanabe et al. 2006).After deuteration and drying, the relative accessibility of hydroxyls to water can be determined as the ratio of integrated areas of O-D to the sum of O-H and O-D vibrations (Suchy et al. 2010b;Taniguchi et al. 1966) or alternatively as the ratio of weighted intensities of O-D to the sum of O-H and O-D vibrations (Mann and Marrinan 1956b;Sepall and Mason 1961).These calculations are based on the Beer-Lambert law which states that the absorbance of a given vibration is given by the product of the three following parameters: abundance of the vibration ) vibrations decrease, while the C=O stretching (* 1750 cm -1 ) and CH 3 deformations (* 1350-1250 cm -1 ) increase.Authors' own data in the material, molar absorptivity of the vibration, and path length of the light beam through the material.For calculating the relative accessibility of hydroxyls to water, it is implicitly assumed that molar absorptivities of O-H and O-D vibrations are equal.This is not necessarily correct as O-D stretching vibrations in cellulose, starch, and wood have been reported to be about 10% more intense than O-H stretching vibrations (Mann and Marrinan 1956b;Nara et al. 1981;Thybring et al. 2017).From theoretical calculations, O-D stretching vibrations are expected to be less intense than O-H stretching vibrations (Crawford 1952;Swenson 1965), in line with studies on liquid normal and heavy water showing 29% lower intensity for O-D stretching vibrations (Venyaminov and Prendergast 1997).Moreover, it should be noted, as discussed for gravimetric techniques, that deuteration will not exchange hydrogens on all hydroxyls interacting with water as shown by Lindh et al. (2016), only those capable of acting as hydrogen bond donor, i.e. in which the hydrogen of hydroxyls is involved in a hydrogen bond with a water molecule.Water-interacting hydroxyls restricted to act as hydrogen bond acceptors, i.e.only interacting with water through their oxygen, are presumably only found on the surfaces of the stiff, aggregated cellulose microfibrils.Before the effect of these issues on the determined hydroxyl accessibility has been resolved for wood and plant materials, it is uncertain whether vibrational spectroscopy can yield quantitative results.Nonetheless, vibrational spectroscopy can be a useful technique for getting qualitative results for interactions between water and wood hydroxyls (see Fig. 7).Both IR and NIR spectroscopy have been applied to study untreated wood (Altaner et al. 2006;Chow 1972;Fackler and Schwanninger 2011;Fernandes et al. 2011;Suchy et al. 2010a, b;Taniguchi et al. 1966;Tsuchikawa and Siesler 2003a, b), while no published data on modified wood are available.Similarly, no known study uses Raman spectroscopy to probe the wood-water interactions by deuteration, whereas a few have been published on deuteration of proteins (Bolton and Scherer 1989;Sa ´nchez-Gonza ´lez et al. 2008) and lactose (Whiteside et al. 2008).
Even more information about hydroxyls can be gained from 2D correlation infrared spectroscopy (Noda 1990(Noda , 2006;;Noda et al. 1993).By this technique, the IR spectrometer is coupled with a dynamic excitation which could be a change in climate (RH or temperature) or a mechanical excitation.This has, for instance, been done using dynamic mechanical analysis (Hinterstoisser and Salme ´n 2000;Hofstetter et al. 2006;Stevanic and Salme ´n 2006;A ˚kerholm and Salme ´n 2003).From a correlation analysis between obtained spectra and external excitation, information about which vibrations in the rather broad and unstructured O-H stretching region (3700-3000 cm -1 ) are affected by external stimulus can be gained, for example the location and orientation of hydroxyls in cellulose.

Conclusion
The interaction of water with wood is highly important for the performance of wood materials.This review describes a range of experimental techniques for characterising water within wood in the entire moisture domain from dry to fully watersaturated.The techniques fall into three broad categories: (1) gravimetric techniques that determine how much water is absorbed, (2) fibre saturation techniques that determine the amount of water within cell walls, and (3) spectroscopic techniques that provide insights into chemical wood-water interactions as well as yield information on water distribution in the macro-void wood structure.For each technique, advantages and limitations are discussed.The experimental data obtained from the various techniques show the many different aspects of wood-water interactions which can be probed.Thus, by combination of several experimental techniques, a wealth of information can be extracted about the state and distribution of water in wood, if the techniques are properly employed.The aim of this review is that more researchers become aware of the wide suite of available experimental techniques and by taking advantage of these techniques, will provide new insights into water interacting with wood materials.

Fig. 2 a
Fig.2a Sorption isotherms for desorption from water-saturated state and absorption from dry state for untreated and acetylated (R mod = 0.18) Norway spruce in the range of 0-97% relative humidity (RH).Based on data fromFredriksson and Johansson (2016),Fredriksson et al. (2013) at 20 °C and from authors at 25 °C (dm/dt = 0.0005%/min).Note that the desorption and absorption isotherms do not coalesce at 97% RH; b absolute sorption hysteresis of untreated wood from moisture content data acquired with automated sorption balances.Based on data for Norway spruce(Fredriksson and Johansson 2016;Fredriksson and Thygesen 2017;Ka ¨llbom 2015;Sedighi-Gilani and Schwarze 2015) and Sitka spruce(Hill et al. 2009).Note that in all studies except(Fredriksson and Johansson 2016;Fredriksson and Thygesen 2017) desorption is initiated from a non-saturated state (conditioning to 95% RH)

Fig. 5 a
Fig. 5 a Diameter of dextran and PEG probe molecules as a function of molecular weight (based on data fromArond and Frank 1954;Carpita et al. 1979;Day et al. 1978;Flournoy et al. 1991;Granath and Kvist 1967;Hill et al. 2005;Ishizawa et al. 2007;Kuga 1981;Lin et al. 1987;Peters 1986;Stone et al. 1969;Stone and Scallan 1968a;Van Dyke 1972;Wang et al. 2011;Wong et al. 1988).b Inaccessible cell wall water volume as a function of size of different molecular probes.Data for untreated and acetylated (R mod = 0.21) Corsican pine fromForster (1998).Note that the fibre saturation point (FSP) can be determined from the plateau in inaccessible cell wall volume for probe sizes too large to enter cell walls

Fig. 7
Fig. 7 Infrared vibrational spectra determined with ATR-FTIR equipment on untreated and acetylated (R mod = 0.13) Norway spruce after deuteration in liquid D 2 O. Average curves based on 4-5 replicates are shown.Note that for acetylated wood both O-H stretching (* 3400 cm -1 ) and O-D stretching (* 2400 cm -1 ) vibrations decrease, while the C=O stretching (* 1750 cm -1 ) and CH 3 deformations (* 1350-1250 cm -1 ) increase.Authors' own data (Becker et al. 196918) 52:297-329  307inlet to the automated sorption balance, e.g. with molecular sieves to get rid of tracer amounts of normal water.Even high-purity gas (like grade 5) has enough water to slightly reverse deuteration in the last drying stage.Lastly, one important issue should be mentioned concerning the RH level(s) used for D 2 O vapour conditioning or their evaluation by RH sensors.For automated sorption balances, the calibration of generated RH or built-in RH sensors is commonly done using saturated salt solutions which generate known equilibrium RH at a given temperature.However, two issues complicate the use of D 2 O in a similar way: the saturated vapour pressure of D 2 O is lower than of H 2 O, and saturated salt solutions with D 2 O generate a slightly (1-3% RH) lower RH than their H 2 O counterparts(Becker et al. 1969).These differences are important to remember if the deuteration technique is used to generate D 2 O sorption isotherms or deuterate a material at a specific RH or vapour pressure.