X-ray Computed Tomography as a Tool for Screening Sediment Cores: An Application to the Lagoons of the Po River Delta (Italy)

: Lagoon sediments have heterogeneous structure and texture, contain shells and plants and are often highly bioturbated and disturbed by human activities. In such sediments, the selection of representative cores and the choice of a subsampling strategy are important but difﬁcult. In this study, we examine the usefulness of X-ray computed tomography (CT) for inferring sediment features that will help in making optimal decisions prior to core opening (24 cores from seven lagoons). Various algorithms (intensity projections, slice thickness, axial and sagittal images, CT number proﬁles and volumetric region of interest) are tested to visualise low- and high-density volumes or objects and to quantify the relations between the average volumetric CT number and the bulk density of the sediment matrix. The CT number is related mainly to water content and indirectly to total nitrogen and <16- µ m grain-size fraction (model R 2 = 0.94). The outliers are attributed to a weak correspondence between the fraction of sediment sampled for water content determination and the volume of sediment matrix used for CT number measurements in highly heterogeneous sediment slices. In conclusion, CT is a powerful tool for the initial screening of cores recovered from heterogeneous lagoon sediments. The adequate use of available algorithms may provide quantitative information on various sediment features, allowing the purposeful selection of cores and subsamples for further investigation. Comprehensive descriptions of the morphological, hydrological aspects of the mean and on the results and and in the


Introduction
Since the early work of Orsi et al. [1], X-ray computed tomography (CT) has been increasingly used for the characterisation of marine and freshwater sediment cores. This analytical technique allows the non-destructive examination of bulk density, porosity, multiplanar visualisation of structures and 3D images of larger objects, burrows and gas bubbles, and it appears as a valuable tool for controlling core quality and identifying lithology (e.g. [2][3][4][5][6][7][8]).
In the CT-scan analysis, the degree of attenuation of the X-radiation through the sediment core is quantified by the CT number that expresses the density range of the material on an arbitrary scale (Hounsfield units). Tomographic images are obtained by using mathematical algorithms for image reconstruction, and specialised software has advanced applications that permit detailed 2D and 3D observations of the core (e.g. [2] and references therein).
In this work, CT-scan analysis is used as an exploratory technique in the study of 24 sediment cores collected in seven lagoons of the Po River Delta (Italy). Sediments deposited in highly dynamic lagoons are characterised by strong heterogeneity, with a deposited in highly dynamic lagoons are characterised by strong heterogeneity, with a variety of sediment structures formed by materials of fluvial and marine origins, as well as authigenic components, such as shells and plants. After deposition, benthic fauna modifies the sediments by changing their geotechnical properties and erosion resistance [9]. Human activities, such as fishing and dredging, further disturb deposited sediments. Until now, lagoon sediments have not been largely studied with X-ray CT techniques. The applications include the examination of stromatolites from a hypersaline lagoon [10], the estimation of the bulk density of surface sediment in the lagoon of Venice for studies on sediment erosion [11], and the identification of sedimentary facies from the Pappas Lagoon in Greece [12]. In the present study, we explore various software algorithms for optimal observation of the structure of sediment cores prior to opening. We also examine the relations between the CT number and physical and structural properties of sediment layers, particularly the water content (WC) of the sediment matrix, in a large number of samples from different sediment types. Finally, the relations between the CT number and other properties of sediments, such as texture and organic matter (OM) content, are considered in an attempt to evaluate the limits of the X-ray CT method for predicting sediment features prior to core subsampling and analyses.

Study Site and Sediment Core Sampling
Before entering the Adriatic Sea, the Po River (the largest Italian watercourse: 650 km long, 71,000 km 2 basin surface area, 1455 m 3 /s mean discharge [13]) forms a delta with five main distributaries ( Figure 1). Two other minor streams cross the delta (Po di Levante and Po di Volano), both former branches of the Po but currently disconnected from the main watercourse. The delta includes a dozen brackish lagoons covering 18% of the 180-km 2 deltaic area [14], which are used especially for extensive aquaculture, comprising mainly mussels and Manila clams [15]. Sediments from the seven main lagoons (from north to south in Figure 1: Caleri (CL), Marinetta-Vallona (MV), Barbamarco (BA), Canarin (CN), Bonelli (BO), Scardovari (SC) and Goro (GO)) were investigated. In the spring of 2016, 24 sediment cores (35.7-59.3 cm long) in Plexiglas cylinders (internal diameter = 6.4 cm) were collected using a piston corer [16]. The investigated lagoons have diversified characteristics in terms of morphology, hydrodynamics, salinity, freshwater inputs and water exchanges with the coastal zone, although they belong to the same deltaic system. The average water depth varies from 0.9 m The investigated lagoons have diversified characteristics in terms of morphology, hydrodynamics, salinity, freshwater inputs and water exchanges with the coastal zone, although they belong to the same deltaic system. The average water depth varies from 0.9 m (CN) to 1.8 m (MV). Comprehensive descriptions of the morphological, hydrological and naturalistic aspects of the lagoons [17]; water circulation, mean salinity and renewal times, based on the results of a hydrodynamic model [18]; and surface sediment characteristics and heavy metal contamination in the seven lagoons [16] have been reported.

CT-Scan Examination
Sediment cores were analysed by X-ray CT scan with a 1-mm down-core resolution using a Philips (Amsterdam, The Netherlands) MX8000 IDT 16 System (voltage = 120 kVp, scanning dose = 211-393 mA). Tomographic intensity profiles measured along the core (CT number, expressed in HU units on the Hounsfield scale) depend on the thickness of the sediment slab considered and the software setting (i.e., the intensity of the X-ray, as well as the algorithms used) [19]. Maximum (MIP) and minimum (MinIP) intensity projection algorithms, respectively, return the highest and the lowest attenuation values encountered in a given sediment volume, highlighting the presence of high-or low-density volumes. The average intensity projection (AIP) algorithm returns the average of each component attenuation value encountered.
Intrinsic inhomogeneity in the sediment layers, as well as the presence of high-density objects (e.g., shells, shell fragments and pebbles) or low-density structures (e.g., those due to the presence of voids, water/gas "bubbles" [20][21][22], worm burrows, remains of plant roots), contribute to uncertainties in the CT number (see Supplementary Materials).
To obtain a representative value and reduce the degree of uncertainty, the CT number was measured under AIP conditions in 20-mm-thick axial slices of the core, following the procedure described in Figure 2. This procedure was chosen for optimal comparison of the CT measurement with sediment features in 20-mm-thick subsamples passed through a 1-mm mesh for water content, organic matter and grain-size measurements. (CN) to 1.8 m (MV). Comprehensive descriptions of the morphological, hydrological and naturalistic aspects of the lagoons [17]; water circulation, mean salinity and renewal times, based on the results of a hydrodynamic model [18]; and surface sediment characteristics and heavy metal contamination in the seven lagoons [16] have been reported.

CT-Scan Examination
Sediment cores were analysed by X-ray CT scan with a 1-mm down-core resolution using a Philips (Amsterdam, The Netherlands) MX8000 IDT 16 System (voltage = 120 kVp, scanning dose = 211-393 mA). Tomographic intensity profiles measured along the core (CT number, expressed in HU units on the Hounsfield scale) depend on the thickness of the sediment slab considered and the software setting (i.e., the intensity of the X-ray, as well as the algorithms used) [19]. Maximum (MIP) and minimum (MinIP) intensity projection algorithms, respectively, return the highest and the lowest attenuation values encountered in a given sediment volume, highlighting the presence of high-or low-density volumes. The average intensity projection (AIP) algorithm returns the average of each component attenuation value encountered.
Intrinsic inhomogeneity in the sediment layers, as well as the presence of high-density objects (e.g., shells, shell fragments and pebbles) or low-density structures (e.g., those due to the presence of voids, water/gas "bubbles" [20][21][22], worm burrows, remains of plant roots), contribute to uncertainties in the CT number (see Supplementary Materials).
To obtain a representative value and reduce the degree of uncertainty, the CT number was measured under AIP conditions in 20-mm-thick axial slices of the core, following the procedure described in Figure 2. This procedure was chosen for optimal comparison of the CT measurement with sediment features in 20-mm-thick subsamples passed through a 1-mm mesh for water content, organic matter and grain-size measurements. High-density volumes (Figure 2a-d) were identified under MIP conditions as the region of interest (ROI), isolated with contours and then, using the ROI editing procedure, were discarded in the recalculation of the CT number under AIP conditions. A similar procedure was used under MinIP conditions to discard large low-density volumes (Figure High-density volumes (Figure 2a-d) were identified under MIP conditions as the region of interest (ROI), isolated with contours and then, using the ROI editing procedure, were discarded in the recalculation of the CT number under AIP conditions. A similar procedure was used under MinIP conditions to discard large low-density volumes (Figure 2e-h). Both kinds of ROIs were carefully observed before being excluded using the 2D and 3D tools made available by the software. Unless otherwise stated, from here onwards, the term CT number is meant for measurements obtained under AIP conditions using the procedure illustrated in Figure 2.
The CT numbers measured in this manner were better correlated with the W C of sediments than the CT numbers measured under MIP or MinIP algorithms. More details on the procedure are described in Supplementary Materials.
In each core, the procedure described above was applied to all the 2-cm-thick axial slices from which sediment aliquots were taken for the analyses described in Section 2.3. The final CT number was then compared with various features of the sediment subsample sieved through a 1-mm mesh. Obviously, after the core opening, free water and gases (from bubbles) were lost.

Core Subsampling and Sediment Analyses
The Plexiglas cylinders were split lengthwise, and one-half sediment cores were sampled by selecting 2-cm-thick layers at different depths along the profile; the top (0-10 cm) sections were continuously sampled, producing 5 samples, and the layers below were selected on the basis of CT images. Wet sediment samples were sieved through a 1-mm Teflon mesh to remove shells and other debris and thoroughly homogenised. W C was determined from an aliquot of a sample dried in an oven at 105 • C until it reached a constant weight [23]. Based on variation coefficients of 7 replicate measurements in 2 sediment samples, the analytical error of W C determination was less than 0.65%.
Other variables were derived from W C , following Hobbs' assumptions and formulas [24]. The assumptions were as follows: (i) The pore space in the sediment is completely water saturated (no gases); (ii) The density of pore water ( w ) is 1.0 g cm −3 ; (iii) The mean density of dry sediment particles ( p ) is correctly estimated. As the terminology in soil and sediment studies is not uniform, we specify the definitions and the equations used in this study, as follows: W c is the weight fraction of water in the total weight of wet sediment, where W w denotes water weight and W s means dry sediment weight. For a unit weight of water sediment mix (W w + W s = 1), W c in g equals W w . Porosity (Φ) is the volume fraction of water in the total volume of wet sediment, where V w signifies water volume and V s refers to dry sediment volume.
V w = W c / w = W c , as water density w = 1 where p is the mean sediment particle density in g cm −3 ; thus, Bulk density ( bw ) is the weight of wet sediment in a unit volume (g cm −3 ) or bw = Φ/W c It should be noted that the bulk density used to match the values of the CT number was calculated from W C in a fraction of sediment matrix after removing all larger shells or shell fragments and other debris. More precisely, it was thus the bulk density of the sediment matrix rather than the bulk density of the total sediment.
The mean sediment particle density ( p ) was estimated by assuming constant densities of mineral fraction ( m ) and OM fraction ( o ). Following Avnimelech et al. [25], we assumed the conversion factor of organic carbon (OC) to OM to be 1.7, and thus: where OM = 1.7 × OC in %, and o = 1.25 and m = 2.65 g cm −3 were assumed for the density of organic and mineral fraction, respectively [25].
W C was not corrected for salt, as this correction is considered negligible ( [26], cited after [24]). Based on the modelled mean surface water salinity on various sites in lagoons [18], we estimated that introducing a correction for salt would change the calculated bulk density value by a maximum of 0.33%.
Grain-size analysis was carried out on a wet aliquot sample of about 2 g, using a Mastersizer 2000 analyser (Mastersizer 2000, Malvern Instruments, Malvern, UK), which provided the volumetric percentages of particles within 100 contiguous dimension classes covering the 0.02-2000-µm range.
Total nitrogen (TN) and total carbon (TC) were determined on duplicate samples using a ThermoFisher Flash 2000 IRMS Elemental Analyzer (EA, ThermoFisher Scientific Inc., Aurora, OH, USA). Subsamples for OC measurements were first decarbonated with HCl 1.5 N. Determined by replicate analyses of the same sample, each measurement had an average standard deviation (STD) of ±0.07% for OC and ±0.009% for TN. Inorganic carbon (IC) was calculated as the difference between TC and OC.
About 0.4 g of the sample dried at 55 • C was digested with 8 mL of HNO 3 8N in a microwave oven [27], and after filtration, the iron (Fe) concentration was determined by Inductively Coupled Plasma-Atomic Emission Spectrometry (ICP-AES-Optima 2100DV, PerkinElmer, Waltham, MA, USA) [28]. Determined using six replicates of homogenised samples, the analytical precision was 6%. A certified reference material sample (PACS-2-Marine Sediment Reference Materials for Trace Metals and other Constituents, National Research Council of Canada) was digested and analysed five times in accordance with the sediment sample analysis procedure. The mean recovery was 75%, indicating that a part of Fe was bound to refractory minerals (e.g., alumina silicates and/or some resistant sulphides). The measured content of Fe is thus considered an extractable fraction. The calibration for ICP-AES analysis was achieved using external standards via the multiplepoint curve approach. Full calibration was performed after each set of 48 samples.
The CT-scan data were processed with Horos v2.0.2 (Horos Project, Free DICOM Medical Image Viewer), an open-source code software distributed free of charge under the Lesser General Public License at Horosproject.org and sponsored by Nimble Co LLC d/b/a Purview in Annapolis, MD, USA.
Statistical calculations were performed using PAST 3.4 for Windows (University of Oslo, Oslo, Norway) and StatPlus for Mac (AnalystSoft Inc., Walnut, Canada).

Sediment Characteristics
The statistical data for the distribution of the relevant variables measured in the sediment samples are shown in Table 1. Distributions of water, fine particles (<16 µm) in the clay plus fine silt dimensional classes and OC percentages with sediment depth are shown in Figure 3. Table 1. Statistical description of the relevant variables measured in the 252 sediment samples collected from the 24 cores. (1) calculated from water content (W C ) using Equations (3) and (4); (2)  It is worth pointing out that the set of 252 samples shows large ranges in the values of the main physico-chemical variables considered, as evidenced by the coefficient of variation in percent (CV%) values (Table 1), particularly for water, fine particles and TN and OC contents. Such a well-diversified set of sediment samples is suitable for the evaluation of the usefulness of the CT-scan investigation in reflecting sediment features prior to core opening. Table 1. Statistical description of the relevant variables measured in the 252 sediment samples collected from the 24 cores. (1) calculated from water content (WC) using Equations (3) and (4); (2) percentage of particles with diameter <16 μm; (3) Fe concentration extractable with 8N HNO3; (4) coefficient of variation.

CT Number, Bulk Density and Water Content
The mean CT numbers obtained from the selected volume of the examined slices, as shown in Figure 2 (Section 2.2), are compared with the matrix bulk density calculated from the measured WC and the assumed density of sediment particles (Equations (4) and (5)). A relatively narrow range of mean particle densities (2.59-2.65 g cm −3 ) is estimated from the proportions of mineral and organic particles. The calculated matrix bulk density varies largely (1.3-2.0 g cm −3 ) because of the considerable range of WC.
In Figure 4, the CT numbers are plotted against the calculated bulk densities of 252 samples from 24 cores recovered from the 7 lagoons. For comparison, calibration fits Fine particle distribution depends on the hydrodynamic characteristics, sediment transport from the Po River distributaries and processes that occur along the coastline. A higher content of fine particles in the surface sediment layer (0-10 cm) is observed in the lagoons located south of the main distributary (Po di Pila), with the exception of the sites close to the tidal inlets (CN11, SC15) or near the river branch (BO13), where strong currents prevent their sedimentation [16]. On a few sites, sediments in the deeper part of the cores are depleted in fine fraction (e.g., CL3, MV5, GO18), suggesting more energetic hydrodynamics in the past.
OC content generally decreases with depth, but in some cores, clear minima or maxima are observed at intermediate depths, indicating modifications of the environmental setting or the sediment delivery.
It is worth pointing out that the set of 252 samples shows large ranges in the values of the main physico-chemical variables considered, as evidenced by the coefficient of variation in percent (CV%) values (Table 1), particularly for water, fine particles and TN and OC contents. Such a well-diversified set of sediment samples is suitable for the evaluation of the usefulness of the CT-scan investigation in reflecting sediment features prior to core opening.

CT Number, Bulk Density and Water Content
The mean CT numbers obtained from the selected volume of the examined slices, as shown in Figure 2 (Section 2.2), are compared with the matrix bulk density calculated from the measured W C and the assumed density of sediment particles (Equations (4) and (5)). A relatively narrow range of mean particle densities (2.59-2.65 g cm −3 ) is estimated from the proportions of mineral and organic particles. The calculated matrix bulk density varies largely (1.3-2.0 g cm −3 ) because of the considerable range of W C .
In Figure 4, the CT numbers are plotted against the calculated bulk densities of 252 samples from 24 cores recovered from the 7 lagoons. For comparison, calibration fits obtained by other authors [29][30][31] who used various density standards and a similar setting of the CT scanner, particularly an X-ray energy of 120 keV, are also shown. Although the determination coefficient for our data is lower than those obtained for calibration fits using various standards, the differences in regression coefficients are minor ( obtained by other authors [29][30][31] who used various density standards and a similar setting of the CT scanner, particularly an X-ray energy of 120 keV, are also shown. Although the determination coefficient for our data is lower than those obtained for calibration fits using various standards, the differences in regression coefficients are minor ( Table 2).  Table 2). Table 2. Statistics of regressions shown in Figure 4. Regressions 2-4 are obtained using various homogeneous and heterogeneous density standards, such as water, bentonites, wood and rocks. All CT numbers are obtained using X-ray energy of 120 keV. Considering a large set of samples with broad ranges of WC, sediment texture and some variations in chemical composition (Table 1), the fact that the calculated bulk density explains nearly 90% of the variability of the CT number is a quite satisfactory result, better than that reported in a study [32] from a lake sediment core using direct measurements of bulk density (R 2 = 0.59).

Regression
Although the methods for deriving the bulk density value from the CT number seem reliable, the underlying principles are not always well constrained for various types of sediments. Orsi et al. [1] studied shallow-water Louisiana shelf sediments and found a generally good correlation between the CT number and the bulk density, measured with a traditional method. They noticed occasional deviations from expected values, explained by the presence of shell fragments and gas bubbles, as well as imperfect adequacy between volumes measured with a CT scan and sampled for bulk density determinations. Consequently, the CT number value derived from the total volume of a sediment slab corresponds to a "true" bulk density, while the selective measurement that avoids the obvious heterogeneities due to shells or bubbles reflects a "matrix" density. The variability of the bulk density of recent, water-saturated sediments depends essentially on WC, OM content and the mean density of the mineral phase. These parameters are strongly interrelated in soils and sediments [25,33,34].
As shown in Figure 5a, the CT number depends inversely and linearly on porosity (Φ), with the determination coefficient similar to the coefficient found for the CT numberbulk density relation. This is because porosity is calculated from the measured WC and the  Table 2). Figure 4. Regressions 2-4 are obtained using various homogeneous and heterogeneous density standards, such as water, bentonites, wood and rocks. All CT numbers are obtained using X-ray energy of 120 keV. Considering a large set of samples with broad ranges of W C , sediment texture and some variations in chemical composition (Table 1), the fact that the calculated bulk density explains nearly 90% of the variability of the CT number is a quite satisfactory result, better than that reported in a study [32] from a lake sediment core using direct measurements of bulk density (R 2 = 0.59).

Regression
Although the methods for deriving the bulk density value from the CT number seem reliable, the underlying principles are not always well constrained for various types of sediments. Orsi et al. [1] studied shallow-water Louisiana shelf sediments and found a generally good correlation between the CT number and the bulk density, measured with a traditional method. They noticed occasional deviations from expected values, explained by the presence of shell fragments and gas bubbles, as well as imperfect adequacy between volumes measured with a CT scan and sampled for bulk density determinations. Consequently, the CT number value derived from the total volume of a sediment slab corresponds to a "true" bulk density, while the selective measurement that avoids the obvious heterogeneities due to shells or bubbles reflects a "matrix" density. The variability of the bulk density of recent, water-saturated sediments depends essentially on W C , OM content and the mean density of the mineral phase. These parameters are strongly interrelated in soils and sediments [25,33,34]. Figure 5a, the CT number depends inversely and linearly on porosity (Φ), with the determination coefficient similar to the coefficient found for the CT number-bulk density relation. This is because porosity is calculated from the measured W C and the estimated mean sediment particle density ( p ), combining Equations (3) and (5). Although calculated porosity explains 89% of the CT number variability, other factors must significantly contribute to some data points' considerable deviations from the best fit. It is evident that the uncertainty of p cannot explain these deviations, as shown by the blue and the red regression lines depicting the linear relations, assuming the mean specific density of sediment particles fixed at 2.5 and 2.8 g cm −3 , respectively. This range of the mean density of particles is likely valid for all samples examined in this study. Fortin et al. [32] reported a determination coefficient (R 2 ) of 0.83 for the linear regression CT number versus porosity (calculated from wet and dry bulk density), comparable to R 2 = 0.89 in this study. They observed that the correlation was weaker for lower porosity and attributed this tendency to the errors in the discrete sampling for bulk density measurements; in fact, these errors correspond to a percentage of the mass per cubic centimetre and hence they increase in denser samples. generally good correlation between the CT number and the bulk density, measured with a traditional method. They noticed occasional deviations from expected values, explained by the presence of shell fragments and gas bubbles, as well as imperfect adequacy between volumes measured with a CT scan and sampled for bulk density determinations. Consequently, the CT number value derived from the total volume of a sediment slab corresponds to a "true" bulk density, while the selective measurement that avoids the obvious heterogeneities due to shells or bubbles reflects a "matrix" density. The variability of the bulk density of recent, water-saturated sediments depends essentially on WC, OM content and the mean density of the mineral phase. These parameters are strongly interrelated in soils and sediments [25,33,34]. As shown in Figure 5a, the CT number depends inversely and linearly on porosity (Φ), with the determination coefficient similar to the coefficient found for the CT numberbulk density relation. This is because porosity is calculated from the measured WC and the estimated mean sediment particle density (ρp), combining Equations (3) and (5). Although calculated porosity explains 89% of the CT number variability, other factors must significantly contribute to some data points' considerable deviations from the best fit. It is evident that the uncertainty of ρp cannot explain these deviations, as shown by the blue and the red regression lines depicting the linear relations, assuming the mean specific density of sediment particles fixed at 2.5 and 2.8 g cm −3 , respectively. This range of the mean density of particles is likely valid for all samples examined in this study. Fortin et al. [32] reported a determination coefficient (R 2 ) of 0.83 for the linear regression CT number versus porosity (calculated from wet and dry bulk density), comparable to R 2 = 0.89 in this study. They observed that the correlation was weaker for lower porosity and attributed this tendency to the errors in the discrete sampling for bulk density measurements; in fact, these errors correspond to a percentage of the mass per cubic centimetre and hence they increase in denser samples. and linear (not shown) fits to the data points are very similar. Blue and red lines show the expected relations between water content and CT number, calculated using Equation (4) with two particle densities of 2.5 and 2.8 g cm −3 , respectively, and coefficients resulting from CT versus bulk density regression (Figure 5c). (c) Calculated bulk density against CT number plotted as dependent variable. Error bars show 1 STD of the mean CT number measured on selected volumes, as described in Section 2.2. Samples with residuals larger than 3 STDs of CT number are filled with red and blue for positive and negative residuals, respectively. Samples with residuals larger than 2 STDs of CT number are filled with orange and light blue for positive and negative residuals, respectively.  linear (not shown) fits to the data points are very similar. Blue and red lines show the expected relations between water content and CT number, calculated using Equation (4) with two particle densities of 2.5 and 2.8 g cm −3 , respectively, and coefficients resulting from CT versus bulk density regression (Figure 5c). (c) Calculated bulk density against CT number plotted as dependent variable. Error bars show 1 STD of the mean CT number measured on selected volumes, as described in Section 2.2. Samples with residuals larger than 3 STDs of CT number are filled with red and blue for positive and negative residuals, respectively. Samples with residuals larger than 2 STDs of CT number are filled with orange and light blue for positive and negative residuals, respectively.

As shown in
At a W C exceeding 40%, p has very little influence on the CT number, but at a lower W C , the importance of p increases (Figure 5b). For example, at 60% W C , the difference between the CT numbers calculated with p of 2.5 and 2.8 is 40, while it is about 150 for 20% W C . It should be noted that linear and polynomial fits to our experimental data are very similar. It might be that at a very high W C , part of the water is lost during subsampling; thus, W C is underestimated in these samples, increasing linearity.
A closer examination of the plot of the calculated matrix bulk density against the mean CT number (Figure 5c) reveals strong differences in STDs of the mean CT number and some distinct outliers. The STD increases with the heterogeneity of sediments. An inspection of residuals of the regression indicates that 32 out of the 252 (12.7%) samples deviate from the best fit by more than 2 STDs of the mean CT number. Among them, 17 (6.7%) samples deviate by more than 3 STDs. These samples are marked with different colours in Figure 5c.
Most probably, some deviations are due to poor correspondence between the measured mean CT number of the sediment matrix and the measured W C . On one hand, the W C is determined with a small aliquot of sediment. On the other hand, in some sediment layers with abundant shells or numerous low-density volumes, CT numbers are determined using only a part of the total layer volume. Another reason could be the underestimation of the W C due to a loss of a fraction of water during subsampling, which may occur especially at a high W C .
It is unlikely that strong deviations are due to variations in the mean sediment particle density, unless heavy minerals are strongly enriched in some grain-size classes. Such enrichment can occur due to size-density sorting [35], but no data on heavy minerals in the recent sediment from lagoons in the Po Delta are available.
It is interesting that the majority of the samples showing high positive residuals of regression (Figure 5c) are found in cores from the CL Lagoon. In contrast, most of the samples that have negative residuals are collected from the GO Lagoon. These two lagoons are situated at the northern and the southern extremities of the modern Po Delta (Figure 1), respectively, and have accumulated sediments of different origins and experienced different environmental conditions [16,18].
X-ray attenuation expressed as the CT number is sensitive not only to the Compton effect, which is related to bulk density, but also to the photoelectric effect, which increases with the atomic number (Z) of the elements present in sediments. Duchense et al. [36] demonstrated that the photoelectric effect increased for aragonite and calcite with respect to quartz, the minerals common in lagoon sediments. Orsi and Anderson [37] showed that CT numbers could be precisely calibrated in terms of the bulk density using an artificial mixture of SiO 2 and seawater, although they noticed a difference in calibration when using quartz or carbonates as sediment surrogates. This is because the photoelectric effect is more important for Ca (atomic number = 20) than for Si (atomic number = 14). However, at an X-ray energy of >100 keV, variations in the effective atomic numbers of sediments forming major elements have little impact on the CT number due to the small photoelectric effect [34]. The difference in CT numbers increases with decreasing W C or increasing bulk density. For example, at bulk densities of 1.3 and 2.0 g cm −3 , the CT number for calcite compared with silica increases by 12% and 20%, respectively, at an X-ray energy of 120 kVp [38].
However, in our study, the samples showing positive residuals of regression ( Figure 5c) higher than three STDs of the CT number, all from the CL Lagoon, have relatively elevated W C (mean = 50%, range = 36-61%), corresponding to a mean bulk density of 1.6 g cm −3 (range = 1.4-1.8). In these samples, the measured CT numbers are 15-50% higher than those predicted from regression. Neither carbonates nor extractable Fe contents are significantly higher in these samples compared with the remaining samples from the CL Lagoon.
Consequently, it seems that Ca from carbonates or Fe oxides/sulphides are not responsible for the increased CT number. At present, the CL Lagoon is not directly connected to the Po River distributaries, but deltaic deposits in this region are alternatively alimented by the Po and the Adige Rivers. The composition of major elements of sediments carried by these two rivers is not distinct enough to provide a fingerprint of provenance [39]. More detailed studies would be necessary to confirm and possibly explain the observed deviations, but some hints can be obtained from the available data (Sections 3.3 and 3.4).

Relations between Water Content and Sediment Characteristics
The W C of sediments is well correlated with the percentage of fine silt and clay fractions (<16 µm), with OC and TN contents. The regressions of W C against these variables are shown in Figure 6a, and the corresponding statistics are listed in Table 3. Generally, fine-grained and organic-rich unconsolidated sediments retain more water than sandy and organic-poor sediments (e.g. [25]). A comprehensive discussion of the relations between bulk density, W C and sediment composition in surface shelf and lagoon sediments is presented by Flemming and Delafontaine [40]. The higher determination coefficient for the regression of W C against TN as compared with OC can probably be explained by concomitant preferential nitrogen mineralisation of OM and sediment compaction with depth. The linearity of the relations between W C and both TN or OC seems lost for the highest values of independent variables ( Figure 6a); thus, second-degree polynomial fits represent these relations better (Table 3). It seems that when TN and OC contents exceed 0.2 and 1.8%, respectively, water saturation capacity is reached. On the contrary, the relation between WC and <16-µm fraction remains linear for the full range of measured values.  Table 3.
Stepwise multiple linear regression of WC against the percentage of TN and <16-μm μm fraction explains 84% of WC variability (Figure 7a), whereas an addition of OC as the third variable does not improve R 2 adj. anymore (Table 3). In 10 out of the 252 samples, measured values are higher than predicted values by more than 20%, all from the cores collected in the CL Lagoon. Among these 10 samples, 8 also show over 20% deviations from the best fit of the CT number versus bulk density (Figure 5c).  Table 3.

Relations between CT Number and Sediment Characteristics
Boespflug et al. [41] used CT number profiles for the detection of lamination in St. Laurence estuary sediments and observed an inverse correlation of the CT number with sediment grain size and OM measured in light and dark laminae. As these variables control the WC of sediments, the relations were probably indirect.
Here, the relations between the CT number and the sediments' TN, OC and <16-μm fraction are tested by simple (Figure 6b) and stepwise multiple regressions (Figure 7b,c  Table 3. Stepwise multiple linear regression of W C against the percentage of TN and <16-µm fraction explains 84% of W C variability (Figure 7a), whereas an addition of OC as the third variable does not improve R 2 adj. anymore (Table 3). In 10 out of the 252 samples, measured values are higher than predicted values by more than 20%, all from the cores collected in the CL Lagoon. Among these 10 samples, 8 also show over 20% deviations from the best fit of the CT number versus bulk density (Figure 5c).  Table 3.
Stepwise multiple linear regression of WC against the percentage of TN and <16-μm μm fraction explains 84% of WC variability (Figure 7a), whereas an addition of OC as the third variable does not improve R 2 adj. anymore (Table 3). In 10 out of the 252 samples, measured values are higher than predicted values by more than 20%, all from the cores collected in the CL Lagoon. Among these 10 samples, 8 also show over 20% deviations from the best fit of the CT number versus bulk density (Figure 5c).  Table 3.

Relations between CT Number and Sediment Characteristics
Boespflug et al. [41] used CT number profiles for the detection of lamination in St. Laurence estuary sediments and observed an inverse correlation of the CT number with sediment grain size and OM measured in light and dark laminae. As these variables control the WC of sediments, the relations were probably indirect.
Here, the relations between the CT number and the sediments' TN, OC and <16-μm fraction are tested by simple (Figure 6b) and stepwise multiple regressions (Figure 7b,c  Table 3.

Relations between CT Number and Sediment Characteristics
Boespflug et al. [41] used CT number profiles for the detection of lamination in St. Laurence estuary sediments and observed an inverse correlation of the CT number with sediment grain size and OM measured in light and dark laminae. As these variables control the W C of sediments, the relations were probably indirect.
Here, the relations between the CT number and the sediments' TN, OC and <16-µm fraction are tested by simple (Figure 6b) and stepwise multiple regressions (Figure 7b,c and Table 3). These relations are similar to those for W C , but all determination coefficients are higher. Multiple linear regression with just two independent variables (TN and <16-µm fraction) explains 90% of the variability of the CT number compared with 84% for W C (Table 3). Moreover, none of the 10 samples from the CL Lagoon cores showing more than 20% deviations from predicted values in W C multiple regression (Figure 7a) against the same independent variables strongly deviates from predicted CT number values (Figure 7b). This suggests that W C values, rather than CT number measurements, are not representative for these samples.
Although the reason for inadequate W C values is not clear, it should be noted that some CL Lagoon samples show considerable heterogeneity due to their high content of shell fragments, and possibly measured CT number values reflect true W C better than direct measurements in a small subsample. Five samples show higher measured CT number values than the predicted ones by more than 20%. These samples, all from the top layers of three GO Lagoon cores, have TN contents among the highest measured (0.23-0.32%). At these contents, the relations between the CT number or W C and TN lose their linearity ( Figure 6). When W C is added to the independent variables, this problem is partly corrected, increasing the R 2 adj. to 0.94 (multiple regression 4 in Table 3 and Figure 7c).

Conclusions
Various software algorithms available for the CT-scan investigation of sediment cores have been explored to optimise conditions for the observation of sediment structure along vertical profiles and the visualisation of high-density objects (Maximal Intensity Projection) or low-density volumes (Minimal Intensity Projection). CT numbers obtained with Average Intensity Projection reflect the bulk density of the sediment matrix throughout the lagoons of the Po Delta, and the segmentation procedure with the selection of regions of interest reduces measurement uncertainty.
The water content of sediments and the CT number are well correlated. Multiple linear regression shows that two variables, total nitrogen and <16-µm fraction, explain 90% of the CT number variability.
The measurement of the CT number prior to core opening provides not only a proxy of sediment water content but also an indirect estimation of fine fraction and organic matter (organic carbon and total nitrogen) contents.