The Hydrodynamic Dispersion Characteristics of Coral Sands

Dispersion characteristics are important factors affecting groundwater solute transport in porous media. In marine environments, solute dispersion leads to the formation of freshwater aquifers under islands. In this study, a series of model tests were designed to explore the relationship between the dispersion characteristics of solute in calcareous sands and the particle size, degree of compactness, and gradation of porous media, with a discussion of the types of dispersion mechanisms in coral sands. It was found that the particle size of coral sands was an important parameter affecting the dispersion coefficient, with the dispersion coefficient increasing with particle size. Gradation was also an important factor affecting the dispersion coefficient of coral sands, with the dispersion coefficient increasing with increasing d10. The dispersion coefficient of coral sands decreased approximately linearly with increasing compactness. The rate of decrease was −0.7244 for single-grained coral sands of particle size 0.25–0.5 mm. When the solute concentrations and particle sizes increased, the limiting concentration gradients at equilibrium decreased. In this study, based on the relative weights of molecular diffusion versus mechanical dispersion under different flow velocity conditions, the dispersion mechanisms were classified into five types, and for each type, a corresponding flow velocity limit was derived.


Introduction
Coral sands are a kind of biogenic soil in marine environments, originating from the fracturing and sedimentation of coral skeletons by wind and waves [1]. The primary features of coral sands include a high content of calcium carbonate-often more than 90%-and a relatively short sediment transport distance, which results in particles frequently retaining the inner pore structure of coral bones, and are characterized by irregular shapes, high angularity, and fragility [2][3][4]. The sediments formed by coral sands have a high void ratio and display hydrogeological characteristics different from those of terrigenous sediments. Islands in the South China Sea are all composed of coral sediments, with these sediments being the only carriers in the formation of fresh groundwater aquifers in the South China Sea islands and reefs. Given that the dispersion coefficient of coral sands is a key factor affecting the conservation of fresh groundwater, uncovering the dispersion pattern of groundwater solute in coral sands will provide the basic parameters and theoretical basis for the numerical simulation of the formation and evolution of fresh groundwater aquifers in the South China Sea islands and reefs, as well as the conservation and utilization of these aquifers.
Solute dispersion in hydrated porous media has been extensively studied. The earliest work was performed by Taylor, who used a capillary tube model to investigate this topic and proposed a method vade-zone, the study of hydrodynamic dispersion is helpful in solving the problem of the effect of fertilizers and pesticides on underground water quality in islands. Finally, in terms of sewage treatment, hydrodynamic dispersion is helpful in solving the problem of the impact of sewage discharge from living and production on the underground freshwater quality in islands.

Test Scheme
The coral sands used for testing were taken from a natural reef in the South China Sea, whose gradation characteristics are shown in Figure 1. The coefficient of curvature C c was 2.12, which is within the range of 1-3. The coefficient of uniformity C u was 45.45, which is much greater than 5, indicating that the reef consisted of a type of coral sand with good gradation. The coral sand was screened to obtain six kinds of coral sand with a single particle size, as shown in Figure 2. More specifically, the test scheme was composed of three steps: (1) using a self-designed one-dimensional dispersion test device, solute dispersion tests were conducted on a total of 14 groups of coral sands of a single particle size and dry density under various conditions to explore how the dispersion coefficient varies as a function of particle size and dry density; (2) using a custom-designed pore tortuosity test device, where pore tortuosity tests were conducted on a total of five groups of coral sands of a single particle size under various conditions to explore how the pore tortuosity varies as a function of particle size; and (3) using a custom-designed molecular diffusion and mechanical dispersion test device, where for a total of six test groups, molecular diffusion tests were conducted with different concentrations of injected solutions under the condition of the pore fluid having a flow velocity. In addition, molecular diffusion and mechanical dispersion tests were conducted with pore flow velocity increasing in a stepwise manner in order to explore how the weight of molecular diffusion versus mechanical diffusion varies with a stepwise increase in pore flow velocity. The test scheme is shown in Table 1. treatment, hydrodynamic dispersion is helpful in solving the problem of the impact of sewage discharge from living and production on the underground freshwater quality in islands.

Test Scheme
The coral sands used for testing were taken from a natural reef in the South China Sea, whose gradation characteristics are shown in Figure 1. The coefficient of curvature Cc was 2.12, which is within the range of 1-3. The coefficient of uniformity Cu was 45.45, which is much greater than 5, indicating that the reef consisted of a type of coral sand with good gradation. The coral sand was screened to obtain six kinds of coral sand with a single particle size, as shown in Figure 2. More specifically, the test scheme was composed of three steps: (1) using a self-designed one-dimensional dispersion test device, solute dispersion tests were conducted on a total of 14 groups of coral sands of a single particle size and dry density under various conditions to explore how the dispersion coefficient varies as a function of particle size and dry density; (2) using a custom-designed pore tortuosity test device, where pore tortuosity tests were conducted on a total of five groups of coral sands of a single particle size under various conditions to explore how the pore tortuosity varies as a function of particle size; and (3) using a custom-designed molecular diffusion and mechanical dispersion test device, where for a total of six test groups, molecular diffusion tests were conducted with different concentrations of injected solutions under the condition of the pore fluid having a flow velocity. In addition, molecular diffusion and mechanical dispersion tests were conducted with pore flow velocity increasing in a stepwise manner in order to explore how the weight of molecular diffusion versus mechanical diffusion varies with a stepwise increase in pore flow velocity. The test scheme is shown in Table 1.

One-Dimensional Dispersion Test
The hydrodynamic dispersion coefficient is a tensor related to the average pore velocity and the characteristics of the porous medium. Studies show that the dispersion is directional even in an isotropic medium and is more complex when the medium is anisotropic [15]. Given this information, only the one-dimensional dispersion characteristics of a solute in coral sands were explored in this study. The test was conducted with a custom-designed one-dimensional dispersion test device, as shown in Figure 3. The complete test device consisted of four parts (denoted by numbers in the figure): dispersion columns (01), a freshwater supply tank (02), a tracer supply tank (03), and a data acquisition system (04). The dispersion columns were composed of several organic glass columns; each column was 8 cm in inner diameter and 50 cm in height. The preparation and tests of the porous medium samples were performed in the dispersion columns, with a 5-cm-thick buffer layer of glass beads placed separately at the top and bottom of the samples. The sensors used in the data acquisition system were CS655 multi-parameter sensors (Campbell Scientific, State of California, America), which can simultaneously measure volumetric moisture content, temperature, conductivity, dielectric constant, signal propagation time, and signal attenuation. The tracer was a 20 g/L NaCl solution and the samples were single-grained coral sands. The test was conducted by continuously injecting the tracer and collecting the data at a fixed location. The samples in the dispersion columns were first saturated with fresh water, then the valve of the freshwater supply tank was shut, followed by opening the valve of the tracer supply tank and starting data acquisition in a synchronous manner. The relative concentration of NaCl was calculated according to C = (C cj − C 0 )/(C max − C 0 ), with C cj denoting the NaCl concentration at time j. C 0 was the NaCl concentration at the initial time, and C max was the maximum concentration (final steady concentration) in the test.

One-Dimensional Dispersion Test
The hydrodynamic dispersion coefficient is a tensor related to the average pore velocity and the characteristics of the porous medium. Studies show that the dispersion is directional even in an isotropic medium and is more complex when the medium is anisotropic [15]. Given this information, only the one-dimensional dispersion characteristics of a solute in coral sands were explored in this study. The test was conducted with a custom-designed one-dimensional dispersion test device, as shown in Figure 3. The complete test device consisted of four parts (denoted by numbers in the figure): dispersion columns (01), a freshwater supply tank (02), a tracer supply tank (03), and a data acquisition system (04). The dispersion columns were composed of several organic glass columns; each column was 8 cm in inner diameter and 50 cm in height. The preparation and tests of the porous medium samples were performed in the dispersion columns, with a 5-cm-thick buffer layer of glass beads placed separately at the top and bottom of the samples. The sensors used in the data acquisition system were CS655 multi-parameter sensors (Campbell Scientific, State of California, America), which can simultaneously measure volumetric moisture content, temperature, conductivity, dielectric constant, signal propagation time, and signal attenuation. The tracer was a 20 g/L NaCl solution and the samples were single-grained coral sands. The test was conducted by continuously injecting the tracer and collecting the data at a fixed location. The samples in the dispersion columns were first saturated with fresh water, then the valve of the freshwater supply tank was shut, followed by opening the valve of the tracer supply tank and starting data acquisition in a synchronous manner. The relative concentration of NaCl was calculated according to C = (Ccj − C0)/(Cmax − C0), with Ccj denoting the NaCl concentration at time j. C0 was the NaCl concentration at the initial time, and Cmax was the maximum concentration (final steady concentration) in the test.  The concentration distribution was derived by solving the one-dimensional steady flow problem in a semi-infinite soil column with a constant concentration at one end. With t (s) as the time since the tracer has been turned on, x (cm) as the distance of the multi-parameter sensor from the tracer injection port, E(x, t) (dS/m) as the conductivity measured at time t, E0 (dS/m) as the conductivity of the tracer, and u (cm/s) as the average flow velocity of the pore fluid in the soil column, the longitudinal dispersion coefficient DL (cm 2 /s) is calculated as given below, and the coordinate system is shown in Figure 4. The concentration distribution was derived by solving the one-dimensional steady flow problem in a semi-infinite soil column with a constant concentration at one end. With t (s) as the time since the tracer has been turned on, x (cm) as the distance of the multi-parameter sensor from the tracer injection port, E(x, t) (dS/m) as the conductivity measured at time t, E 0 (dS/m) as the conductivity of the tracer, and u (cm/s) as the average flow velocity of the pore fluid in the soil column, the longitudinal dispersion coefficient D L (cm 2 /s) is calculated as given below, and the coordinate system is shown in Figure 4. The above problem is converted into a mathematical equation as follows: The Laplace transform of Equation (1) with respect to t gives is a function of x, and p is a parameter. The original problem is converted into the following definite solution problem as an ordinary differential equation: The solution of the second-order homogeneous linear ordinary differential equation is: where erfc is the complementary error function, and ξ in When x or t is large: The above problem is converted into a mathematical equation as follows: The Laplace transform of Equation (1) with respect to t gives where E(x, p) is a function of x, and p is a parameter. The original problem is converted into the following definite solution problem as an ordinary differential equation: The solution of the second-order homogeneous linear ordinary differential equation is: where erfc is the complementary error function, and ξ in er f c(u) = 1 − 2 √ π u 0 e −ξ 2 dξ is the mathematical expectation value.
When x or t is large: Let ξ 2 = η 2 /2 and dξ = 1/ √ 2dη (where η is the substitute for the mathematical expectation value ξ and does not have practical meaning), then we have: Then, according to Equation (6): According to the normal distribution table, N(−1) = 0.1587, N(1) = 0.8413, thereby leading to: Solving this equation leads to the dispersion coefficient as follows: Here, t 0.1587 is the time when

Molecular Diffusion and Pore Tortuosity Tests
Hydrodynamic dispersion is a process that includes molecular diffusion and mechanical dispersion. Molecular diffusion is the process of molecular transport associated with the stochastic movement of molecules due to a concentration gradient. Mechanical dispersion is the process of mechanical mixing that takes place in porous media as a result of the movement of fluid through the pore space. When the average pore flow velocity is 0, dispersion takes the form of molecular diffusion. In this study, each sample consisted of single-grained coral sands, and there was a total of six test groups, with particle sizes ranging from 0-0.1 mm, 0.1-0.25 mm, 0.25-0.5 mm, 0.5-1 mm, and 1-2 mm, while the dry density of each sample was a fixed value: ρ d = 1.3 g/cm 3 . The tracer and sensor types were the same as those in the one-dimensional dispersion test. The test instrument was a custom-designed and built molecular diffusion device, as shown in Figure 5. There was a control partition in the middle of the device, and after the partition was lifted, the fluids in the left and right sample box compartments could flow back and forth between the compartments. The sizes of the left and right sample boxes were both 10 cm × 10 cm × 10 cm, and the left sample was saturated with NaCl solution while the right sample was saturated with fresh water. Timing and data reading were started simultaneously with the lifting of the partition. Let EN-t and EW-t be the measured conductivity at time t on the left and right sides, respectively, and EM-t = (EN-t + EW-t)/2 be the mean value of the two measured conductivities at time t.

Molecular Diffusion and Mechanical Dispersion Tests
When the pore flow velocity is greater than 0, molecular diffusion and mechanical dispersion Let EN-t and EW-t be the measured conductivity at time t on the left and right sides, respectively, and EM-t = (EN-t + EW-t)/2 be the mean value of the two measured conductivities at time t.

Molecular Diffusion and Mechanical Dispersion Tests
When the pore flow velocity is greater than 0, molecular diffusion and mechanical dispersion processes coexist. With an increase in average pore flow velocity, the weight of molecular diffusion versus mechanical dispersion changes. In order to further study the relationship between the weight and the flow velocity, a custom-designed molecular diffusion and mechanical dispersion test device was employed in this study, as shown in Figure 6. The device consisted of three components, namely, a dispersion body, a water supply part, and a data acquisition system. The dispersion body was made up of organic glass tubes and samples, with each tube being 8 cm in inner diameter and 80 cm in length. The samples were coral sands with particle sizes of 0.25-0.5 mm and a dry density of 1.3 g/cm 3 . The data acquisition system was equipped with CS655 multi-parameter sensors; the tracer was a 20 g/cm 3 NaCl solution.

Effect of Particle Size on the Dispersion Coefficient
In order to uncover the effect of particle size on the one-dimensional dispersion coefficient, dispersion tests were performed on a total of six groups of single-grained coral sand samples, with the sample groups having a particle sizes of < 0.1 mm, 0.1-0.25 mm, 0.25-0.5 mm, 0.5-1 mm, 1-2 mm, and 2-5 mm; the groups had particle size ranges within the categories of silt, fine sand, medium sand, coarse sand, gravelly sand, and crushed stone (angular gravel), respectively. These groupings were determined according to the soil classification method in the Code for the Investigation of Geotechnical Engineering (GB50021-2001) (Ministry of Construction of the People's Republic of China, 2009) [16]. The dry density of all samples was 1.3 g/cm 3 . Figure 7 presents a one-dimensional dispersion curve for single-grained coral sands, which indicates that when the particle size was 0-0.1 mm, the dispersion process consisted of three stages: a drainage stage, a displacement stage, and a stabilization stage. In the drainage stage, the original saturated fluid in the soil column was discharged outward under the displacement of the tracer, a circumstance in which the discharged fluid had the same concentration as the original saturated fluid, thereby leading to a flat dispersion curve. In the displacement stage, with the continuous injection of the tracer, the tracer underwent diffusion in the original saturated solution, in addition to producing the displacement effect. The diffusion "interface" existed in the form of a concentration transition The required coral sand samples were first prepared in the dispersion tubes. The samples were then saturated with fresh water, after which the inlet and outlet valves of the dispersion tube were shut. The NaCl solution was then injected into the tracer injection port and the injection port was shut. This was completed while simultaneously starting data acquisition, with the pore flow velocity in the device being 0. The above steps were repeated using different concentrations of NaCl solution.
Next, dispersion tests with pore flow velocities greater than 0 were conducted according to the above steps, with the exception that the inlet and outlet valves of the dispersion tube were reopened after the injection of the tracer to allow the pore fluid at a certain flow velocity to pass through the samples. The pore flow velocity was controlled by adjusting the output power of the water pump. By repeating the above test procedure, dispersion tests were performed at several different pore flow velocities.

Effect of Particle Size on the Dispersion Coefficient
In order to uncover the effect of particle size on the one-dimensional dispersion coefficient, dispersion tests were performed on a total of six groups of single-grained coral sand samples, with the sample groups having a particle sizes of < 0.1 mm, 0.1-0.25 mm, 0.25-0.5 mm, 0.5-1 mm, 1-2 mm, and 2-5 mm; the groups had particle size ranges within the categories of silt, fine sand, medium sand, coarse sand, gravelly sand, and crushed stone (angular gravel), respectively. These groupings were determined according to the soil classification method in the Code for the Investigation of Geotechnical Engineering (GB50021-2001) (Ministry of Construction of the People's Republic of China, 2009) [16]. The dry density of all samples was 1.3 g/cm 3 . Figure 7 presents a one-dimensional dispersion curve for single-grained coral sands, which indicates that when the particle size was 0-0.1 mm, the dispersion process consisted of three stages: a drainage stage, a displacement stage, and a stabilization stage. In the drainage stage, the original saturated fluid in the soil column was discharged outward under the displacement of the tracer, a circumstance in which the discharged fluid had the same concentration as the original saturated fluid, thereby leading to a flat dispersion curve. In the displacement stage, with the continuous injection of the tracer, the tracer underwent diffusion in the original saturated solution, in addition to producing the displacement effect. The diffusion "interface" existed in the form of a concentration transition zone between the original saturated fluid and the tracer. When the edge of the transition zone reached the sensor position, the measured concentration of the discharged fluid began to increase, indicative of the onset of the displacement stage, as manifested by the curve slope starting to increase. In the stabilization stage, when the tracer had completely displaced the original saturated fluid, the discharged fluid had a concentration similar to that of the tracer, with the curve reaching the highest value and tending to flatten.  As revealed by the above observations, different sample groups of different particle sizes exhibited dispersion pattern similarities, although the drainage stage gradually shortened with increasing particle size. The curve's slope in the displacement stage, i.e., the diffusion rate, increased with increasing particle size, and the time taken for the concentration to reach the stabilization stage gradually shortened with increasing particle size. With the exception of the 0-0.1 mm group, the other three groups were relatively similar to each other in terms of the parameters of the three stages, indicating a significant difference in dispersion characteristics between the particle size 0-0.1 mm and the other particle sizes. Figure 8 shows the variation pattern of the one-dimensional dispersion coefficient with respect to particle size. The dispersion coefficient of coral sands increased with increasing particle size. When the particle size exceeded 0.25 mm, the dispersion coefficient increased rapidly, with the dispersion As revealed by the above observations, different sample groups of different particle sizes exhibited dispersion pattern similarities, although the drainage stage gradually shortened with increasing particle size. The curve's slope in the displacement stage, i.e., the diffusion rate, increased with increasing particle size, and the time taken for the concentration to reach the stabilization stage gradually shortened with increasing particle size. With the exception of the 0-0.1 mm group, the other three groups were relatively similar to each other in terms of the parameters of the three stages, indicating a significant difference in dispersion characteristics between the particle size 0-0.1 mm and the other particle sizes. Figure 8 shows the variation pattern of the one-dimensional dispersion coefficient with respect to particle size. The dispersion coefficient of coral sands increased with increasing particle size. When the particle size exceeded 0.25 mm, the dispersion coefficient increased rapidly, with the dispersion coefficient differing by a factor of 202. When the particle size was 0.1-2 mm (i.e., fine sands-gravels), the dispersion coefficient was between 0.152-2.37 cm 2 /s, with the difference being a factor of approximately 15. However, when the particle size was greater than 2 mm, the change in the dispersion coefficient with respect to particle size became smaller, with the dispersion coefficient tending to reach a constant value as the particle size increased, finally reaching a value of 2.5 cm 2 /s in this test condition. These observations suggest that 0.25 mm and 2 mm could be considered the characteristic particle sizes of coral sands.

Effect of Dry Density on the Dispersion Coefficient
The effect of the degree of compactness on the one-dimensional dispersion coefficient was investigated using the group of coral sands having a particle size range of 0.25-0.5 mm and dry densities of 1.2 g/cm 3 (relative degree of compactness of 1.021), 1.3 g/cm 3 (relative degree of compactness of 1.212), and 1.4 g/cm 3 (relative degree of compactness of 1.375). The test results, which are shown in Figure 9, indicate that with the increase in density, the displacement stage became longer, the diffusion rate decreased, and the time spent before reaching a steady state increased. Figure 10 represents the variation of the dispersion coefficient with the change in dry density, which reveals that the dispersion coefficient decreased linearly with increasing dry density.

Effect of Dry Density on the Dispersion Coefficient
The effect of the degree of compactness on the one-dimensional dispersion coefficient was investigated using the group of coral sands having a particle size range of 0.25-0.5 mm and dry densities of 1.2 g/cm 3 (relative degree of compactness of 1.021), 1.3 g/cm 3 (relative degree of compactness of 1.212), and 1.4 g/cm 3 (relative degree of compactness of 1.375). The test results, which are shown in Figure 9, indicate that with the increase in density, the displacement stage became longer, the diffusion rate decreased, and the time spent before reaching a steady state increased. Figure 10 represents the variation of the dispersion coefficient with the change in dry density, which reveals that the dispersion coefficient decreased linearly with increasing dry density. are shown in Figure 9, indicate that with the increase in density, the displacement stage became longer, the diffusion rate decreased, and the time spent before reaching a steady state increased. Figure 10 represents the variation of the dispersion coefficient with the change in dry density, which reveals that the dispersion coefficient decreased linearly with increasing dry density.

Effect of Particle Gradation on the Dispersion Coefficient
To simulate natural-gradation sands, stepwise removal of the particles smaller than a certain size was conducted in order to change the gradation [17]. This was followed by one-dimensional dispersion tests on coral sands of different gradations to investigate the effect of particle gradation on the one-dimensional dispersion coefficient. The sample dry density was 1.3 g/cm 3 and the gradation parameters of each sample are listed in Table 2 [18]. As shown by the data, all the samples were poorly graded coral sands, with the exception of the natural gradation sands. The gradation curves of the samples are displayed in Figure 11 and the dispersion coefficient curves are illustrated in Figure 12.

Effect of Particle Gradation on the Dispersion Coefficient
To simulate natural-gradation sands, stepwise removal of the particles smaller than a certain size was conducted in order to change the gradation [17]. This was followed by one-dimensional dispersion tests on coral sands of different gradations to investigate the effect of particle gradation on the one-dimensional dispersion coefficient. The sample dry density was 1.3 g/cm 3 and the gradation parameters of each sample are listed in Table 2 [18]. As shown by the data, all the samples were poorly graded coral sands, with the exception of the natural gradation sands. The gradation curves of the samples are displayed in Figure 11 and the dispersion coefficient curves are illustrated in Figure 12.   As shown in Figure 12, the dispersion coefficient of the coral sands gradually increased with increasing d10. Following the removal of soil particles smaller than 0.25 mm (d10 increased to 0.31 mm from the initial 0.05 mm), the dispersion coefficient of the coral sands increased dramatically. After removal of the soil particles smaller than 0.5 mm (d10 exceeded 0.6 mm), the dispersion coefficient of the coral sands tended to reach a constant value. These results suggest that 0.25 mm and 2 mm were the characteristic particle sizes for the dispersion properties of graded coral sands.

Analysis of the Dispersion Mechanisms in Coral Sands
There are many significant factors affecting solute transport, such as convection, mechanical dispersion, molecular diffusion, interactions between the solid phase and solute (such as dissolution and adsorption), chemical reactions within the solution, and other source-sink solute interactions (such as attenuation of radioactive elements and absorption of certain solutes by crop roots) [19]. Among these, mechanical dispersion and molecular diffusion are collectively referred to as dispersion. Studies have shown that mechanical dispersion and molecular diffusion generally occur simultaneously during solute transport, although they exhibit different variation patterns with changing average pore flow velocity. When the flow velocity is low, molecular diffusion is stronger than mechanical dispersion. In contrast, when the flow velocity is high, mechanical dispersion is stronger than molecular diffusion. In order to reveal the roles of these two mechanisms in the solute dispersion process of coral sands, pore tortuosity tests, as well as molecular and mechanical dispersion tests, were performed on the coral sands in this study. As shown in Figure 12, the dispersion coefficient of the coral sands gradually increased with increasing d 10 . Following the removal of soil particles smaller than 0.25 mm (d 10 increased to 0.31 mm from the initial 0.05 mm), the dispersion coefficient of the coral sands increased dramatically. After removal of the soil particles smaller than 0.5 mm (d 10 exceeded 0.6 mm), the dispersion coefficient of the coral sands tended to reach a constant value. These results suggest that 0.25 mm and 2 mm were the characteristic particle sizes for the dispersion properties of graded coral sands.

Analysis of the Dispersion Mechanisms in Coral Sands
There are many significant factors affecting solute transport, such as convection, mechanical dispersion, molecular diffusion, interactions between the solid phase and solute (such as dissolution and adsorption), chemical reactions within the solution, and other source-sink solute interactions (such as attenuation of radioactive elements and absorption of certain solutes by crop roots) [19]. Among these, mechanical dispersion and molecular diffusion are collectively referred to as dispersion. Studies have shown that mechanical dispersion and molecular diffusion generally occur simultaneously during solute transport, although they exhibit different variation patterns with changing average pore flow velocity. When the flow velocity is low, molecular diffusion is stronger than mechanical dispersion. In contrast, when the flow velocity is high, mechanical dispersion is stronger than molecular diffusion. In order to reveal the roles of these two mechanisms in the solute dispersion process of coral sands, pore tortuosity tests, as well as molecular and mechanical dispersion tests, were performed on the coral sands in this study. Figure 13 depicts a curve showing the diffusion concentrations of solute molecules in coral sands having different particle sizes. E NaCl is the conductivity measured by the left sensor and E water is measured by the right sensor, with E Med = (E NaCl + E water )/2. When there was a concentration gradient in the saturated porous medium and the pore flow velocity was 0, the NaCl solution concentrations gradually decreased and increased in the left and right samples, respectively, until the two concentrations reached a relative equilibrium. It is noteworthy that the concentrations of the NaCl solution on both sides were not necessarily the same when equilibrium was reached; this was due to the fact that when the concentration gradient was low (defined as the limiting concentration gradient at equilibrium, ranging from 0-1, with larger values representing larger concentration gradients at equilibrium), the "obstruction" effect of the porous medium occurred, and the concentration gradient would not continue to decrease. Porous media with different particle sizes (pore sizes) had different limiting concentration gradients at equilibrium. As shown by Figure 13, when the particle size was larger than 0.1 mm, the time spent prior to reaching the limiting concentration gradient at equilibrium gradually decreased with increasing particle size, concomitant with a gradual decrease in the limiting concentration gradient at equilibrium. When the particle size was larger than 1 mm, the limiting concentration gradient at equilibrium was infinitely close to 0. When the particle size was smaller than 0.1 mm, the molecular driving force generated by the concentration gradient was less than the "obstruction" effect of the porous medium; thus, the limiting concentration gradient at equilibrium was rapidly achieved and relatively large (approaching 1). Figure 14 illustrates how the amount of time taken to reach the limiting concentration gradient at equilibrium varied with particle size. As shown in Figure 14, 0.1-0.25 mm was still the characteristic particle size of coral sands for molecular dispersion, and the group of coral sands with this particle size range took the longest time to reach the limiting concentration gradient at equilibrium. This finding is in agreement with the conclusion drawn from the one-dimensional dispersion test conducted in this study. equilibrium), the "obstruction" effect of the porous medium occurred, and the concentration gradient would not continue to decrease. Porous media with different particle sizes (pore sizes) had different limiting concentration gradients at equilibrium. As shown by Figure 13, when the particle size was larger than 0.1 mm, the time spent prior to reaching the limiting concentration gradient at equilibrium gradually decreased with increasing particle size, concomitant with a gradual decrease in the limiting concentration gradient at equilibrium. When the particle size was larger than 1 mm, the limiting concentration gradient at equilibrium was infinitely close to 0. When the particle size was smaller than 0.1 mm, the molecular driving force generated by the concentration gradient was less than the "obstruction" effect of the porous medium; thus, the limiting concentration gradient at equilibrium was rapidly achieved and relatively large (approaching 1). Figure 14 illustrates how the amount of time taken to reach the limiting concentration gradient at equilibrium varied with particle size. As shown in Figure 14, 0.1-0.25 mm was still the characteristic particle size of coral sands for molecular dispersion, and the group of coral sands with this particle size range took the longest time to reach the limiting concentration gradient at equilibrium. This finding is in agreement with the conclusion drawn from the one-dimensional dispersion test conducted in this study.   In order to quantitatively characterize the "obstruction" effect of porous media on molecular diffusion, earlier studies introduced the concept of pore tortuosity (θ), i.e., the ratio of the length of the porous medium sample to the actual path traveled by the fluid particles through a sample of that length [20]. Thus, θ = * ⁄ , where is the dispersion coefficient in the open water body and * is the dispersion coefficient in the porous medium. Given the assumption that the limiting concentration gradient at equilibrium infinitely approaches 0 for molecular diffusion in an open and still water body, the θ of a porous medium should range between 0 and 1. According to the definition of pore tortuosity, θ = S/(v × t), S (cm) is the length of the porous medium sample, v (cm/s, taken as scalar without considering direction) is the velocity of fluid particle movement, and t (s) is the time duration of fluid particle movement, namely, the time spent by the fluid particles prior to reaching the limiting concentration gradient at equilibrium. Based on the data in Figure 13, the time t in each particle size group of coral sands could be obtained, with t set to positive infinity in the case of particle sizes less than 0.1 mm. In addition, based on the time t and the known distance of 4.9 cm between the sensor and the central partition, it was possible to obtain the relationship curve of pore tortuosity versus particle size, as shown in Figure 15, where v is the velocity of fluid particle movement at a certain temperature.

Molecular Diffusion in Coral Sands
A previous study on the dispersion characteristics of terrigenous sediments found that θ ranges from 0.01-0.5 [21]. In addition, this parameter has been proposed to be 0.1 for clays and 0.7 for sandy soils [22]. In order to quantitatively characterize the "obstruction" effect of porous media on molecular diffusion, earlier studies introduced the concept of pore tortuosity (θ), i.e., the ratio of the length of the porous medium sample to the actual path traveled by the fluid particles through a sample of that length [20]. Thus, θ = D 0 /D * , where D 0 is the dispersion coefficient in the open water body and D * is the dispersion coefficient in the porous medium. Given the assumption that the limiting concentration gradient at equilibrium infinitely approaches 0 for molecular diffusion in an open and still water body, the θ of a porous medium should range between 0 and 1. According to the definition of pore tortuosity, θ = S/(v × t), S (cm) is the length of the porous medium sample, v (cm/s, taken as scalar without considering direction) is the velocity of fluid particle movement, and t (s) is the time duration of fluid particle movement, namely, the time spent by the fluid particles prior to reaching the limiting concentration gradient at equilibrium. Based on the data in Figure 13, the time t in each particle size group of coral sands could be obtained, with t set to positive infinity in the case of particle sizes less than 0.1 mm. In addition, based on the time t and the known distance of 4.9 cm between the sensor and the central partition, it was possible to obtain the relationship curve of pore tortuosity versus particle size, as shown in Figure 15, where v is the velocity of fluid particle movement at a certain temperature.  Figure 15. Relationship curve between pore tortuosity and particle size of coral sands.

Mechanical Dispersion in Coral Sands
Dispersion consists of molecular diffusion and mechanical dispersion. When the flow velocity is greater than 0, the two processes usually coexist, albeit with different weights on the overall Figure 15. Relationship curve between pore tortuosity and particle size of coral sands.
A previous study on the dispersion characteristics of terrigenous sediments found that θ ranges from 0.01-0.5 [21]. In addition, this parameter has been proposed to be 0.1 for clays and 0.7 for sandy soils [22].

Mechanical Dispersion in Coral Sands
Dispersion consists of molecular diffusion and mechanical dispersion. When the flow velocity is greater than 0, the two processes usually coexist, albeit with different weights on the overall dispersion coefficient under different conditions (i.e., concentration and flow velocity). Figure 16 presents the variation curves of relative conductivity at a fixed pore flow velocity of 0 for two different concentrations of tracer NaCl solution (20 g/L versus 60 g/L). Here, E C is the sensor-measured conductivity, E 0 is the initial conductivity before the injection of the tracer, and E max is the maximum conductivity measured by each sensor. The rate of increase (slope) of relative conductivity was greater for the 60 g/L NaCl solution than for its 20 g/L counterpart, suggesting that the higher the concentration gradient, the greater the molecular diffusion rate. Figure 15. Relationship curve between pore tortuosity and particle size of coral sands.

Mechanical Dispersion in Coral Sands
Dispersion consists of molecular diffusion and mechanical dispersion. When the flow velocity is greater than 0, the two processes usually coexist, albeit with different weights on the overall dispersion coefficient under different conditions (i.e., concentration and flow velocity). Figure 16 presents the variation curves of relative conductivity at a fixed pore flow velocity of 0 for two different concentrations of tracer NaCl solution (20 g/L versus 60 g/L). Here, EC is the sensormeasured conductivity, E0 is the initial conductivity before the injection of the tracer, and Emax is the maximum conductivity measured by each sensor. The rate of increase (slope) of relative conductivity was greater for the 60 g/L NaCl solution than for its 20 g/L counterpart, suggesting that the higher the concentration gradient, the greater the molecular diffusion rate.   Figure 17 shows the variation curves of relative conductivity at different flow velocities, with Figure 17A-D depicting the variation curves of measured conductivities from the sensors EC_1 to EC_3. Figure 17a-d depict the variation curves of measured conductivities from the sensors EC_4 to EC_6. In relation to the tracer injection port, the sensors EC_1 to EC_3 were downstream, while the sensors EC_4 to EC_6 were upstream. The dispersion process downstream from the port would involve both mechanical dispersion and molecular diffusion, with both taking place in the same direction. In contrast, in the dispersion process occurring against the flow direction upstream from the port, mechanical dispersion would take place in a different direction than molecular diffusion. Therefore, the conductivity change detected by the sensor would be a net conductivity change of the solutes with their molecular diffusion overcoming their mechanical dispersion. Generally, when mechanical dispersion dominates, the shape of concentration-time curve will be as shown in Figure 18a. In this situation, concentration attenuation occurs after the maximum concentration has been reached, with t p denoting the time taken to reach the peak point of the curve. When molecular diffusion dominates, the shape of the concentration-time curve is as shown in Figure 18b. In this case, the concentration remains constant after reaching the maximum value, with t s denoting the time taken to reach the stable point of the curve.
18-a. In this situation, concentration attenuation occurs after the maximum concentration has been reached, with tp denoting the time taken to reach the peak point of the curve. When molecular diffusion dominates, the shape of the concentration-time curve is as shown in Figure 18-b. In this case, the concentration remains constant after reaching the maximum value, with ts denoting the time taken to reach the stable point of the curve. It can be seen that when the flow velocity was 1.36 × 10 −4 cm/s ( Figures 17A and 17a), the curve characteristics of the left figure fell between those of Figure 18-a and Figure 18-b, indicating the simultaneous presence of the two mechanisms of molecular diffusion and mechanical dispersion downstream from the port. Data comparison between the left and right figures indicates that the time to reach the maximum conductivity was similar in the two figures, suggesting that the two mechanisms were comparable in terms of the dispersion process under this condition. The relative conductivity decreased after the peak value in the right figure, which was likely the result of the displacement effect of the pore fluids.
With the increase in flow velocity, the curve shape in the left figure (Fig.17-A,B,C and D)becomes gradually similar to that in Figure 18-a, indicating that mechanical dispersion was gradually enhanced. This was manifested by the gradual shortening of the displacement stage and the It can be seen that when the flow velocity was 1.36 × 10 −4 cm/s ( Figures 17A and 17a), the curve characteristics of the left figure fell between those of Figures 18a and 18b, indicating the simultaneous presence of the two mechanisms of molecular diffusion and mechanical dispersion downstream from the port. Data comparison between the left and right figures indicates that the time to reach the maximum conductivity was similar in the two figures, suggesting that the two mechanisms were comparable in terms of the dispersion process under this condition. The relative conductivity decreased after the peak value in the right figure, which was likely the result of the displacement effect of the pore fluids.
With the increase in flow velocity, the curve shape in the left figure ( Figure 17A-D) becomes gradually similar to that in Figure 18a, indicating that mechanical dispersion was gradually enhanced. This was manifested by the gradual shortening of the displacement stage and the continuous increase in the diffusion rate. In contrast, the curve shape in the right figure (Figure 17a-d) became increasingly atypical, and the time taken to reach the maximum value gradually increased and was much longer than the corresponding time in the left figure, indicating that molecular diffusion was increasingly subject to the flow velocity. By the time the flow velocity reached 6.16 × 10 −3 cm/s, the molecular diffusion upstream from the tracer injection port was basically negligible.
Based on a comprehensive analysis of the above test results, the velocity boundaries of the pore fluids that control the dispersion mechanisms were derived, as shown in Table 3.

Conclusions
The solute diffusion coefficient of coral sands is a key factor affecting fresh groundwater conservation, and knowledge of the freshwater dispersion pattern in coral sands will provide both the basic parameters and a theoretical basis for numerical simulation of the formation and evolution of freshwater aquifers under islands, as well as the conservation and utilization of these aquifers. Using coral sands collected from a reef in the South China Sea, a series of tests and theoretical analyses were conducted in this study, from which the following conclusions were drawn: 1) The particle size of coral sands was an important parameter affecting the dispersion coefficient, which increased with increasing particle size. The reason is that when the volume and porosity of the soil are constant, the number of pores decreased and the pore size increased with the increase in particle size. Because the pore fluid was more likely to mechanically flow in larger pores, the diffusion coefficient, therefore, could increase. The diameters of 0.25 mm and 2 mm were discovered to be the characteristic particle sizes of coral sands due to their dispersion characteristics. The dispersion coefficient could vary by a factor of more than 200 between the 0.1-0.25 mm size group and the 0.25-0.5 mm size group, while the dispersion coefficient increased at a relatively small rate in the group with particle sizes larger than 2 mm. 2) Gradation was also an important factor affecting the dispersion coefficient of coral sands, with the diffusion coefficient increasing with increasing d 10 . When d 10 was larger than 0.31 mm, the dispersion coefficient increased dramatically. In contrast, when d 10 was larger than 0.6 mm, the dispersion coefficient tended to be constant. This was because an increase in d 10 reflected that the filling degree of the fine particles in the soil became worse. Similarly, the pore size and diffusion coefficient increased with an increase in d 10 . However, when d 10 was 0.31 mm or less, the diffusion coefficient did not increase with increasing d 10 . The reason was the increase in pore size was minimal with respect to the macroscopic scale. Alternatively, the importance of pore size for the macro scale grew with an increase in d 10 , and the diffusion coefficient also rose. Finally, when the pore size reached a critical value, the influence of pore size on the solute diffusion gradually reduced, and the diffusion coefficient also approached a stable state.
3) The dispersion coefficient of coral sands decreased linearly with an increasing degree of compactness. The dispersion coefficient decreased at a rate of −0.7244 in the group of coral sands with particle sizes of 0.25-0.5 mm. The reason was that when the particle size distribution remained stable, the increase in compactness directly resulted in a decrease in porosity, which could restrict the flow of the pore liquid. The diffusion coefficient, therefore, was reduced. By analyzing the experimental dates, the diffusion coefficient was reduced when the dry density was increased from 1.3 g/cm 3 to 1.5 g/cm 3 . 4) Both the concentration of solute in the coral sands and the particle size of the porous medium affected the limiting concentration gradient at equilibrium. With an increase in solute concentration and particle size, the limiting concentration gradient at equilibrium decreased. As particle size increased, the pore tortuosity decreased. This was because when the soil volume and porosity were constant, the number of pores decreased and the pore size increased with an increase in the particle size, and the actual distance that the solution molecules flowed from point A to point B in the soil was shortened. Therefore, the calculated porosity curvature was reduced. 5) The dispersion mechanisms in coral sands could be determined by the flow velocity of the pore fluid. Based on the different weights of molecular diffusion versus mechanical dispersion for different conditions of flow velocity, the dispersion mechanisms could be classified into five types.