Comprehensive NMR Analysis of Pore Structures in Superabsorbing Cellulose Nanofiber Aerogels.

Highly porous cellulose nanofiber (CNF) aerogels are promising, environmentally friendly, reusable, and low-cost materials for several advanced environmental, biomedical, and electronic applications. The aerogels have a complex and hierarchical 3D porous network structure with pore sizes ranging from nanometers to hundreds of micrometers. The morphology of the network has a critical role on the performance of aerogels, but it is difficult to characterize thoroughly with traditional techniques. Here, we introduce a combination of nuclear magnetic resonance (NMR) spectroscopy techniques for comprehensive characterization of pore sizes and connectivity in the CNF aerogels. Cyclohexane absorbed in the aerogels was used as a probe fluid. NMR cryoporometry enabled us to characterize the size distribution of nanometer scale pores in between the cellulose nanofibers in the solid matrix of the aerogels. Restricted diffusion of cyclohexane revealed the size distribution of the dominant micrometer scale pores as well as the tortuosity of the pore network. T 2 relaxation filtered microscopic magnetic resonance imaging (MRI) method allowed us to determine the size distribution of the largest, submillimeter scale pores. The NMR techniques are nondestructive, and they provide information about the whole sample volume (not only surfaces). Furthermore, they show how absorbed liquids experience the complex 3D pore structure. Thorough characterization of porous structures is important for understanding the properties of the aerogels and optimizing them for various applications. The introduced comprehensive NMR analysis set is widely usable for a broad range of different kinds of aerogels used in different applications, such as catalysis, batteries, supercapacitors, hydrogen storage, etc.

Introduction

Aerogels are the lightest solid materials on Earth, and they have a complex pore structure and a high inner surface area.[1] Due to their exceptionally good electronic, mechanical, acoustical, optical, microstructural, and thermal properties, aerogels are ideal materials for various applications as adsorbents,[2−4] thermal insulators,[5,6] biomedical materials,[7,8]etc.

In recent years, a new type of aerogels derived from nanocellulose has become widely applied.[4] The nanocellulose aerogels are composed of elongated cellulose nanofibers (CNFs) or rod-shaped nanocrystals with the typical width of 2–4 nm and length reaching even 1 mm. These aerogels have a high porosity and ultralow density. Additionally, due to the base material of which they are made (cellulose) they are biodegradable and environmentally friendly. The nanocellulose aerogels can be readily tailored by surface modifications or by incorporating nanoadditives. Therefore, they are excellent materials for several applications such as catalysis and selective collection of oil and chemical spills from the environment.[9−13]

Porous structure including pore size distribution and pore interconnectivity has a critical role in the applications of aerogels, and therefore, it is highly important to characterize it accurately and thoroughly.[12] There are, however, many problems related to porous structure analysis by standard techniques. For example, scanning electron microscopy provides typically only a 2D view of the surface of the aerogel, and it is difficult to fully understand complex 3D structures or even pore sizes from the planar pictures. Furthermore, high compressive forces of mercury porosimetry may collapse the solid framework of aerogels,[14] and nitrogen adsorption technique may not be able to capture the whole range of pore sizes.[15]

NMR spectroscopy is one of the most powerful tools in chemical analysis.[16] It provides versatile chemical, dynamic, and spatial information. The radiofrequency radiation used in the technique is nondestructive and can penetrate inside even optically opaque materials. Therefore, NMR is an ideal method for investigating dynamics and adsorption phenomena of fluids in porous materials.[17−22] Furthermore, the fluids may provide detailed information about porous structures from nano- to millimeter scale. NMR cryoporometry is based on the fact that the melting point of a liquid in a small pore is lower than that of bulk liquid, and the melting point depression is inversely proportional to the pore size.[23,24] The method relies on the same principle as thermoporosimetry,[25] but it differs in the detection. While the thermoporosimetry utilizes differential scanning calorimeter (DSC), in the NMR cryoporometry the amount of unfrozen liquid in the sample is measured using a spin echo technique filtering out the signal of frozen liquid to determine the melting point distribution, and the distribution is converted into the pore size distribution using the Gibbs–Thomson equation. A thorough comparison between two methods has been provided by Gane et al.[26] The methods can be used to determine pore sizes in the nanometer scale. NMR is also one of the rare methods for measuring self-diffusion of molecules without tracers.[27] The walls of porous materials restrict the diffusion of adsorbed fluids, and therefore the diffusion experiments provide information about pore sizes in micrometer scale as well as pore connectivity.[28−31] Magnetic resonance imaging (MRI) technique, which also relies on the NMR phenomenon, can in turn be used to characterize millimeter scale pores.[32]

NMR is also a useful tool for characterization of aerogels. For example, solid-state NMR has been exploited in the investigation of chemical structures of aerogels.[33,34] Pore structures of aerogels have been studied by MRI and diffusion measurements of gases, such as methane as well as hyperpolarized xenon and helium, adsorbed in aerogels.[35−40] Nanopore size distributions of silica aerogels have been studied with NMR cryoporometry.[41,42]

In this article, we demonstrate that the combination of three different advanced NMR techniques, namely NMR cryoporometry, diffusion NMR, and T2 relaxation weighted MRI, can be used for comprehensive determination of pore sizes of CNF aerogels from nanometer to millimeter scale, as well as the tortuosity of the pore network. The CNF aerogels were manufactured from waste cellulose fibers (recycled milk container board). The manufacturing process included a pretreatment based on a green deep eutectic solvent (DES) and mechanical grinding to form a nanofiber hydrogel, followed by freeze-drying and hydrophobic modification (silylation).[43] The resulting CNF aerogels have ultralow density (down to 0.0029 g/cm3) and high porosity (up to 99.81%).[43] Furthermore, they own outstanding absorption properties for several oils and organic solvents (absorption capacities of up to 142.9 g/g)[43] and are considered as superabsorbents. They are reusable; the absorbed liquids can be recovered by mechanical squeezing, and high absorption capacity is maintained at least for 30 cycles.[43]

Materials and Methods

Materials

Recycled milk container board (MCB) was used as raw materials for fabrication of cellulose nanofibers. MCB was pulped without any additional chemicals using a Kenwood KM020 pulper (U.K.) at a consistency of 15% using temperature of 45 °C. Pulping time was 10 min with rotor speed of 2 (i.e., ∼ 250 rpm). After pulping MCB pulp was washed and screened using a Somerville apparatus (Lorentzen & Wettre, Sweden).

A deep eutectic solvent (DES) was prepared by heating 1620 g of choline chloride and 1223 g of urea in a 5 dm3 beaker at 100 °C until the mixture melted. Next 25 g (abs) of MCB pulp (moisture content of ∼30%) was dosed to the DES for 2 h at 100 °C. Then, the pulp was washed with water using a Somerville apparatus (Lorentzen & Wettre, Sweden) until clear rinsewater was obtained.

Nanofibrillation of MCB was performed using a Masuko supermasscolloider grinder MKCA6–2J (Japan) to obtain cellulose nanofibers (CNF) using negative gap values of −20, −40, −50, −70, and −80 μm to obtain gel-like MCB CNF.

Methyltrimethoxysilane (MTMS) and hexadecyltrimethoxysilane (HDTMS) in ethanol were used for cross-linking and hydrophobization of the CNF. Freshly prepared MTMS and HDTMS solutions (20 wt %) were mixed using a magnetic stirrer for 10 min. Then, 25 wt % of freshly prepared silane solutions (MTMS and HDTMS, ratio of 50:50) against CNF amounts were added by micropipette to the CNF suspensions and stirred with a magnetic stirrer at room temperature for 2 h.

Silylated MCB CNF suspensions were placed directly into 10 mm NMR sample tubes (5 cm height pillar) with various consistencies (i.e., 0.3%, 0.5%, 0.75%, and 1.0%), frozen in liquid nitrogen, and freeze-dried using Scanvac Coolsafe 55–15 Pro (Denmark) for 72 h to obtain hydrophobized and cross-linked CNF aerogels for NMR-analysis. Similarly, aerogel samples around 30 cm3 were prepared for the FESEM analysis using 100 mL plastic beakers. The detailed characteristics of the aerogels are presented in our previous work.[43]

For the magnetic resonance analysis the samples were immersed in cyclohexane (purchased from Merck). The volume of the solvent was chosen to cover whole sample. The NMR tube was then sealed with a cap and which was tightened with Parafilm to avoid solvent evaporation during experiment.

NMR

The experiments were performed by Bruker AVANCE III spectrometer with the static field strength of 11.7 T and 1H resonance frequency of 500 MHz, using a 10 mm broad band observe (BBO) probe.

NMR Cryoporometry

For the T2 relaxation and NMR cryoporometry analysis, a series of CMPG experiments was performed over the temperature range of 239.7–287.4 K, from low to high temperatures. From 239.70 to 266.2 K, the temperature step was about 2 K, from 267.2 to 283.2 K, the step was about 1 K, and from 285.2 to 287.4 K, the step was about 2 K. The experiment series was set using modified TReNDS acquisition script.[44] After each temperature change system was let to stabilize for 7.5 min. Then the probe was automatically tuned and a CPMG experiment was acquired with 8 scans, 512 echoes, 1 ms echo time, 30 s relaxation delay and experiment time of 4 min. The overall heating rate was about 10 and 5 K/h for the steps of 2 and 1 K, respectively, which are slow enough for keeping the sample well in thermal equilibrium. We note that, although the air-filled aerogels are excellent thermal insulators (as the immobilized air is a great thermal insulator), the aerogel filled by cyclohexane is not, and thermal equilibrium is reached relatively fast.

According to the Gibbs–Thomson equation,[45] the melting point depression isHere, T0 and T are the melting points of bulk and confined liquids, respectively, K is a constant characteristic of the liquid, and R is the pore radius. In the NMR cryoporometry analysis, the K of 190.1 K nm was used for cyclohexane.[46]

By assuming that the intensity of NMR signal, I, is proportional to the volume of the pores containing molten liquid at a given temperature, the pore size distribution can be calculated by the following equation:[24,47]

The T2 relaxation time distributions were reconstructed from CPMG data using ITAMeD[48] in MATLAB. After the reconstruction, the region corresponding to the liquid phase (T2 > 6 ms) was integrated for the NMR cryoporometry analysis. The effect of the temperature dependence of the thermal polarization on the intensity was compensated out according to the Curie law,[49] and the NMR cryoporometry pore distributions were calculated using eq .

NMR Diffusometry

The apparent diffusion coefficients D were measured by pulsed gradient stimulated echo (PGSTE) NMR. The diffusion times Δ used in the experiments were 50, 100, 200, 400, and 800 ms and the corresponding gradient pulse lengths were 2, 2, 1.5, 1.06, and 0.75 ms. The experiments were performed at temperature 285 K with 8 scans and 12 s relaxation delay. The gradient strengths were linearly increased from 0.963 to 47.187 G/cm (total 32 values).

The apparent diffusion coefficients were calculated using ITAMeD program.[48] Surface to volume ratios were determined by fitting the following equation to the initial part of the D vs Δ data:[21]

Here, S is the surface of the pores, V is the volume of the pores, and D0 is the diffusion coefficient of bulk liquid.

The tortuosity, α, was determined by fitting the following equation to the whole D vs Δ data:[28]Here, c is defined asand θ is a scaling constant proportional to the pore size. Further, the pore diameter was calculated using the following equation introduced by Lawrence et al.:[24,28−30]

The fittings were done by using Origin 2018b (S/V and tortuosity) and Microsoft Excel 2018 for the diameter.

MRI

MRI experiments were performed by the Bruker Avance III 300 MHz spectrometer equipped with Micro2.5 microimaging unit with x, y, and z gradients, using a 10 mm RF insert. The multi-slice multi-echo (MSME) sequence was used in the experiments. The number of echoes was 16, and the echo time was 20 ms. Altogether eight slices were measured from each sample with 500 μm slice thickness and a field of view of 1 cm. The image was acquired on the grid of 300 × 300 points, and the size of individual pixel was 33 μm. The image was averaged 16 times with the repetition time of 48 s. The axial images were acquired for all four samples, and for the 1% consistency sample, an additional coronal image was acquired.

After the acquisition, the analysis was performed in MATLAB software. For each pixel the T2 value was calculated. Then all the pixels with T2 within the range 350–650 ms where assigned as macro-pores and selected for further analysis. The selected pixels where grouped according to connections with other pixels to separate shapes corresponding to single pore. For each pore the ellipsis was fitted to determine the length, width and angle of main axis. The analysis was possible by using MATLAB functions “bwlabel” and “regionprops” from Image processing toolbox. The pore diameter was defined to be a mean of the longer and shorter axis of the ellipsis.

FESEM

A FEI Quanta 200 FEG ESEM was used for FESEM analysis of the aerogels. The samples were Au-coated using a Quorum Q150T ES (Quorum Technologies, Ashford, Kent, UK) sputter coater for 140 s with a current of 20 mA. A deposition rate of 20 nm/min is typically achieved, however, the porosity and resulting high surface area of the aerogel lead to a much thinner coating (<10 nm thickness on average) throughout the sample. Samples were imaged with an Everhardt Thornley detector (ETD) for secondary electrons at 10 kV in high vacuum. The pore size analysis followed the MRI method. Each FESEM picture was imported to MATLAB. The pixels were classified as part of a pore if they had the intensity lower than 40% of the most intense pixel. Then pixels were grouped and fitted with ellipsis identically to that in MRI analysis

Optical Laser Microscopy

The CNF aerogel samples were examined under a 3D Laser Scanning Confocal Microscope, VK-X200 by Keyence Ltd. To acquire the high-resolution images, the aerogel samples were prepared in circular Petri dishes (90 × 15 mm). Random, representative, and macroscopically flat areas (1.4 × 1.06 mm) were used to obtain the images. The photoreceptor detects the reflected light at each pixel by dividing the planar field of view into 2048 × 1536 pixels to scan a layer of an image. The intensity of each pixel at every position on the Z-axis was obtained with a pitch of 0.1 μm by changing the position of the objective lens. With the Z axis position of the highest intensity being the focal point, the height information and the laser intensity were detected through the aerogel to a depth where no pores could be longer detected.

Such information provides three types of images: a deep field color image, a laser intensity image, and a height map. The final maps were created in the way that the latter two images were overlaid creating a laser intensity image with height information encoded with color.

Results and Discussion

Four hydrophobic CNF aerogel samples with the initial consistencies (mass of cellulose in solvent before freeze-drying during manufacturing) of 0.3%, 0.5%, 0.75%, and 1% were prepared for the NMR analysis. Reference data was collected by field emission scanning electron microscopy (FESEM) and laser microscopy.

Nanometer Scale Pores

These inside the solid matrix of the CNF aerogel were investigated by NMR cryoporometry. Before the analysis, the samples were submerged in cyclohexane in an NMR tube at room temperature. No visible changes in the size and shape of the samples after addition of solvent were observed. 1H CPMG temperature series of the adsorbed cyclohexane was performed over the temperature range of 239.7–328.4 K. Figure A shows the T2 distributions resulting from the Laplace inversion[48] of the CPMG data of the aerogel sample with a consistency of 0.3%. Above the melting point of bulk cyclohexane (280 K), there is only a single component visible with T2 over 100 ms. Below the melting point, cyclohexane forms cubic plastic crystal[50] with a long T2 (slightly below 5 and 1 ms at the higher and lower temperatures, respectively) as compared to typical solids. Below 265 K, the distributions include another, minor component with T2 around 10 ms and intensity smaller than 10% of the main component. That minor component was observed from all the aerogel samples and a neat cyclohexane sample as well, and it was tentatively assigned to arise from crystal defects of cyclohexane.[51]

Figure 1

(A) 1H T2 distributions of cyclohexane adsorbed in the CNF aerogel with the consistency of 0.3% as the function of temperature. (B) Integrals of liquid cyclohexane signals (T2 > 6 ms) of the CNF aerogels with variable consistency (0.3–1%) as well as of the bulk cyclohexane sample. (C) Nanopore size distribution derived from the integrals using the Gibbs–Thomson equation. (D) Mean and median pore sizes of the aerogels.

Figure B plots the integral of the liquid component (T2 > 6 ms) with temperature for all the aerogel samples as well as for bulk cyclohexane. The aerogel samples show significantly higher liquid intensities below 280 K than the bulk cyclohexane sample due to the lowered melting point of cyclohexane in the nanopores of the aerogels. The changes of the integrals revealed the melting point distributions, which were converted into pore size distributions by using the Gibbs–Thomson equation as explained in Materials and Methods. The distributions (see Figure C) include a dominant peak around 0–200 nm, which was interpreted to arise from interfibrous nanopores between the cross-linked cellulose nanofibers, as the size range matches with the widths of the CNF observed by transmission electron microscopy (TEM, 2–80 nm).[43] Previously, a notably smaller pore size of 8 nm with very narrow size distribution obtained from Barrett–Joyner–Halenda analysis of BET measurements was reported for the cellulose nanocrystal (CNC) based aerogels.[52] However, these CNCs had significantly smaller diameter of <10 nm compared to the CNF used in the present work. The nanopore size distributions of CNF aerogels were similar for all four samples, and the median and mean pore sizes (Figure D) varied between 70 and 110 nm. It is likely that interfibrous voids are more affected by cross-linking chemistry and colloidal forces between the CNFs, and a small variation in CNF consistency has only a minor effect on nanoporosity.

Micrometer Scale Pores

These are the dominant pores in the CNF aerogels and were studied via diffusion of adsorbed cyclohexane. The self-diffusion coefficients (D) were determined by the pulsed-field-gradient stimulated-echo (PGSTE) sequence[53] including two magnetic field gradient pulses separated by a diffusion delay Δ, and the results are shown in Figure A. The observed apparent D become shorter with increasing diffusion time because the molecules experience more the restrictions caused by pore walls, and the initial slope of D vs Δ is proportional to the surface-to-volume ratio (S/V) of the material (see Materials and Methods).[21] The S/V values derived from the initial slopes increases with increasing sample consistency from 1.2 × 104 to 2.6 × 104 m–1 (see Figure B). The values were converted into pore diameters by assuming a cylindrical pore geometry (see Figure C). Naturally, the pore geometry of the CNF aerogels is very complex and irregular, but based on FESEM pictures, the cylindrical model was the most appropriate approximation. Opposite to nanopores, here a clear correlation between pore size and consistency of aerogels was noted as the average pore size decreased with increasing consistency from 4.9 to 2.5 μm.

Figure 2

(A) Apparent diffusion coefficients of cyclohexane adsorbed in CNF aerogels as a function of diffusion times. The solid lines show the fits of eq to the data. (B) Surface-to-volume ratios derived from the initial slope of the D vs. Δ curves. (C) Pore diameters calculated from the S/V values assuming a cylindrical pore geometry. (D) Tortuosities of the pore networks resulting from the fits.

The tortuosity factor α is defined as a ratio between the diffusion coefficients corresponding to 0 and infinite diffusion delay (D0 and D∞), and it was determined by fitting eq to the D vs Δ data (see Materials and Methods).[21] The tortuosity of the CNF aerogels was high (see Figure D) due to open and interconnected structure of the porous networks. The tortuosity decreased from 0.95 to 0.92 when the consistency of the samples increased from 0.3% to 1%.

For comparison, the micrometer scale pore structures were investigated by FESEM as well. The FESEM images of the CNF aerogel samples are shown in parts A–D of Figure . The mean and median pore sizes resulting from the automatized analysis of the FESEM images (see Materials and Methods) are plotted in Figure E. Similar to diffusion NMR analysis, the FESEM analysis shows the decreasing trend of the pore size with increasing consistency. However, the pore size values given by FESEM (about 1- 3 μm) are slightly smaller than the values measured by NMR (about 2–5 μm) presumably because the diffusion NMR experiments explore the real 3D pore space, while the FESEM analysis is based on 2D projections.

Figure 3

(A–D) FESEM pictures of the aerogel samples with the consistencies of 0.3%, 0.5%, 0.75%, and 1%. (E) Mean and median pore diameters estimated from the FESEM images. (F) Pore diameter distributions.

The micrometer scale pores of the aerogels are attributed to hierarchical ordering of nanofibers to larger aggregates due to cross-linking and hydrogen bonding. These pores exist between and within the aggregated, sheet-like areas as illustrated by FESEM images (Figure A–D).[43] In addition, the ice crystal formation during aerogel fabrication creates assembled structures of the CNFs within the aerogel.[52] Decreasing trend of micropore size against aerogel consistency is logical due to density increase of the aerogels (presented in our previous work[43]) which decreases the available open space. Overall, the NMR and FESEM analysis are in good qualitative agreement and even quantitatively close to each other, confirming the reliability of the NMR analysis.

Submillimeter Pores

The CNF aerogel samples included also some larger pores with the pore size close to 1 mm, which existed as few inhomogeneous patches between the large sheet-like aggregated areas (Figure ). These pores were investigated by high-resolution microscopic MRI. Figure A shows the spin-density image (slice thickness 500 μm, resolution 33 μm) of cyclohexane adsorbed in the CNF aerogel with the consistency of 1%. There are significant differences in the amounts of adsorbed cyclohexane in different parts of the sample due to the millimeter range heterogeneities of the sample. Corresponding T2 map is shown in Figure B. The T2 maps were utilized to identify the largest, submillimeter pores. T2 relaxation time of bulk cyclohexane at room temperature was determined to be 420 ms. The values in the T2 maps in the range of 350–650 ms were interpreted to arise from bulk like cyclohexane in the large pores. Based on this criterion, the T2 maps were converted into submillimeter pore maps (see Figure C), and the size distributions as well as mean and median pore sizes were calculated (Figure , parts E and F). The mean and median pore sizes were about 200 and 160 μm and quite similar for all the aerogel samples.

Figure 4

(A) Axial spin-density MR image of cyclohexane in the CNF aerogel of 1% consistency. (B) Corresponding T2 map. (C) Image of the submillimeter pores derived from the T2 map using the T2 = 500 ± 150 μs selection criterion. (D) Coronal plane T2 map of the same sample. (E) Distributions of the submillimeter pore diameters in all the aerogel samples. (F) Mean and median values of the pore diameters.

For the comparison, the submillimeter pores were also investigated using a 3D laser scanning confocal microscope. In the laser microscopy images shown in parts A–D of Figure , a clear transition from low to high consistency of aerogels is seen, the reddish and blueish colors indicating higher and lower depths, respectively. The bluish regions highlight the largest pores, and the mean and median sizes of these submillimeter pores were determined from the images (see Figure E). The values vary between 140 and 200 μm, and they are in good agreement with the values given by the MRI analysis (Figure F), which supports its reliability. In the MRI analysis, the pore size distribution was actually derived from eight slices measured from different heights of the sample, resulting in the overall surface area of the analysis of 50 mm2, while in the laser microscopy images the analyzed are was only 1.5 mm2 to ensure sufficient resolution. Therefore, in this regard, the MRI analysis was statistically more comprehensive, and the small deviations between the MRI and laser microscopy values are addressed to the low coverage of the laser microscopy analysis.

Figure 5

(A−D) Superimposed laser microscopy image and height map of the CNF aerogels with consistencies of 0.3%, 0.5%, 0.75%, and 1%, respectively. Red and blue colors correspond to high and low levels, and the large continuous blue regions represent the submillimeter pores. (E) Mean and median pore diameters derived from the laser microscopy images.

Overall Porous Structure

The results of the mean pore sizes measured by different methods are summarized in Figures and 7. The complex and hierarchical 3D structure of aerogels is reflected as trimodal pore size distribution and can be provided by the combination of three different independent NMR techniques, namely NMR cryoporometry, diffusion NMR, and T2 relaxation weighted MRI. The interfibrous small voids between the cross-linked nanofibers are illustrated as nanometer size pores, the size of which is relative independent of aerogel density. The larger micrometer and submillimeter pores are attributed to patches between the aggregated nanofibers, and the size of micrometer pores clearly decrease as a function of aerogel density, while the few submillimeter voids have a relatively constant size. The comparison of results from microscopic visualization with NMR based techniques show good similarity and indicates good reliability of the obtained results.

Figure 6

Summary of mean pore diameters determined by different methods.

Figure 7

Pore size distributions combined from all NMR methods. The intensity of each method is normalized separately. The micrometer scale pores were visualized by assuming a Gaussian pore size distribution with the mean and standard deviation corresponding to the mean pore diameter and statistical error measured by diffusion NMR, respectively.

The interaction between the aerogel material and cyclohexane solvent during the adsorption process may change the structure of the material, and therefore, the NMR analysis may not represent precisely the structure of intact, dry aerogel. However, the changes are expected to be minor, because we did not observe visually any changes in the dimensions of the aerogel piece during the adsorption process. On the other hand, NMR provides the information about the material structure in the conditions for which it was designed, which can be considered to be a strength of the method. As discussed earlier, the CNF aerogel is an organic specific superabsorbent, and therefore, it is highly important to know its structure during the interaction between the solvent and material.

Conclusions

We have proved that combination of a set of advanced NMR techniques can be used for comprehensive analysis of porous networks of novel, environmentally friendly, recyclable, superabsorbing CNF aerogels. NMR cryoporometry provided information about the size distribution of nanometer scale pores in between the cross-linked cellulose nanofibers. The observed pore sizes (0–200 nm) were similar to the widths of the nanofibers and did not vary much with the consistency of the samples. Diffusion NMR experiments revealed the mean pore diameters of the dominant, micrometer pores as well as the tortuosity of the pore network. The mean pore diameter decreased from 4.9 to 2.5 μm when the consistency of the sample increased from 0.3% to 1%. The observed tortuosity was very high due to the open structure of the pore network, and it decreased from 0.95 to 0.92 with increasing consistency. Relaxation filtered MRI showed that the sample included also some larger submillimeter pores with the pore size of about 200 μm. FESEM and laser microscopy confirmed the reliability of the NMR analysis. The NMR toolbox introduced in this article can be broadly applied in the analysis of the pore networks in various aerogels and other porous media including very heterogeneous pore structures.