agronomy A Root Density Tradeoff in an Okra-Assisted Subsurface Pipe Drainage System for Amelioration of Saline Soil

: Subsurface pipe drainage technology can effectively improve coastal saline land in Northern China. We explored an okra ( Abelmoschus esculentus L.)-assisted subsurface pipe drainage system to improve the water and salt discharge performance and beneﬁts. In this study, the simulation box experiment was conducted to research the response of water and salt discharge performance in subsurface pipe drainage to okra root weight density (RWD). The drain outﬂow, soil salinity, and sodium adsorption ratio were determined. The results showed that okra RWD affected the vertical distribution of okra roots. Okra with an appropriate RWD (about 116 µ g · cm − 3 ) could signiﬁcantly increase the cumulative drain outﬂow. Okra with an appropriate RWD (about 136 µ g · cm − 3 ) could signiﬁcantly increase the desalting effect. Moreover, the RWD of okra also inﬂuenced the ability of subsurface pipe drainage to inhibit soil alkalization. The above results show that planting okra and installing subsurface pipe drainage to control drainage at the coastal saline land in Northern China can effectively improve the water and salt drainage effect when okra RWD is about 116–136 µ g · cm − 3 . When using subsurface pipe drainage to improve coastal saline soils, planting okra with proper density may be an appropriate choice to improve the effect and beneﬁt.


Introduction
In North China, there is a large area of coastal saline land which is a potential arable land resource which needs to be improved and utilized [1][2][3]. The subsurface pipe drainage method, has the advantage of reducing the groundwater level and increasing salt-leaching efficiency by installing underground pipes [4,5]. In recent years, relevant research showed that the use of subsurface pipe drainage can effectively improve coastal saline land in Northern China [6,7]. However, various factors constrain the application of the subsurface pipe drainage method. For example, in the early stage of saline land utilization, the soil salinity leads to the low yield of planting common crops [8,9], which will reduce farmers' enthusiasm for investing in subsurface pipe drainage to reclaim saline land. Moreover, the relatively poor soil structure, permeability, and salt-leaching efficiency of coastal saline soil constrain water and salt discharge performance under the utilization of subsurface pipe drainage [10]. In addition, the saline soil reclamation is strongly dependent on fresh water which is scarce in that area [11,12]. Rainfall is the only freshwater resource used for leaching.
Numerous studies have reported that planting plants in saline soil can improve the leaching efficiency of saline soil. It is indicated that the growth of plants has an improved effect on the saline soil through improving the structure and water permeability of soil by its roots [13,14]. Moreover, the CO 2 produced by the respiration of plant roots and the decomposition of organic matter promotes the dissolution of CaCO 3 in the soil, which further reduces the soil pH and exchangeable sodium percentage (ESP), and results in increasing the water permeability of the soil [15][16][17][18]. It is indicated that planting forage beet and barley significantly improved the stability of soil aggregates and the hydraulic conductivity of surface soil [19]. The leaching efficiency under the condition of planting crops was higher than that under the non-cropped condition [20]. The water and salt discharge performance of subsurface pipe drainage may also be enhanced by plants, but the specific conditions are not clear.
The plant selection plays an important role in saline soil utilization [13,21]. Okra (Abelmoschus esculentus (L.) Moench) is recognized as a moderate salt tolerance plant having well-developed roots. It is widely cultivated in the subtropics and tropics due to being rich in minerals and vitamins; thus, it has a high economic value [22,23]. The improving effect of growth of okra on the soil hydraulic property was more than that of wheat [24]. It may result in a great salt-leaching effect under the plantation of okra in saline soil and relatively ideal yield can be obtained in the saline land with leaching measurements [25]. Recently, okra has been successfully planted on saline soil, which achieved good profits at ShouGuang city, the biggest vegetable production base in the north of China [26].
In our present research, we explored the possibility of combining planting okra with subsurface pipe drainage to reclaim saline land. However, okra roots may have a dual effect on salt movement in the saline soil ( Figure 1). On the one hand, salt will accumulate in the root zone due to the transpiration of plants [27], and on the other hand, the root system can improve the hydraulic conductivity and leaching efficiency of saline soil [19,20]. Therefore, we set root densities as treatment variables in the present study. Our hypothesis was that installing underground pipes for subsurface pipe drainage and planting okra with an appropriate root density in a coastal saline area in Northern China will enhance the water and salt discharge performance of subsurface pipe drainage technology. The objective of the present study was to (1) determine the responses of water and salt discharge performance in subsurface pipe drainage to okra root density under the natural rainfall; and (2) find out an appropriate density of okra root density, which will effectively improve the water and salt discharge performance of subsurface pipe drainage.

Experimental Site
A soil box experiment was conducted at Nandagang government (117 • 05 -117 • 49 E, 38 • 09 -38 • 39 N) ( Figure 2) in coastal saline area around Bohai Sea, Hebei Province, North China from April 2019 to November 2020. The soil salt content in this area was 2-6 g·kg −1 , soil type was Salic Fluvisol. The groundwater table range was 0.3 m to 1.2 m. This area belongs to Heilonggang Basin in North China Alluvial Plain. It is a flat coastal plain with many pits and puddles. The climate of this area belongs to continental monsoon with less precipitation in spring and winter, but rich in summer. The average annual precipitation is about 660 mm, about 65% of annual precipitation concentrate in July-August. The average annual evaporation is about 1980 mm, which is 3 times the total annual average precipitation. The effective rainfall in okra growth period (May-November) was 653.1 mm during 2020 ( Figure 3).

Experiment Design and Management
In this study, the simulation boxes with volume of 1.5 × 1.5 × 1.5 m 3 were prepared by pasting 1 cm thick polyethylene plates fixed with steel pipes and the gaps between the plates were welded with plastic covered electrode. A perforated corrugated pipe with a diameter of 4 cm was installed in the middle of the box bottom as a pipe drain, the pipe was sleeved with a sealed nylon mesh cloth, and a valve was installed at the outlet of the pipe drain to control the drainage. After the pipe was installed, quartz sand as filter material was filled 30 cm deep at the bottom of the boxes. Then, the homogeneous coastal saline soil with the bulk density of 1.42 g·cm −3 (average bulk density of local saline soil) was filled in the boxes at 10 cm intervals until the whole soil depth (1 m) was reached. A 20 cm space was left from the soil surface to top of the wall of box for rainfall collection. The soil sample was taken from the local saline soil in the experiment site. After airdrying and passing through 2 mm sieve, the soil sample was mixed and made uniform in salinity during the start of experiment in April 2019. The physical and chemical properties of soil taken to fill the boxes are listed in Table 1. Soil texture was silt soil with a small hydraulic conductivity and high sodium adsorption ratio (SAR). In order to make the soil structure of the simulation box closer to the field condition, the simulation box was left untouched for one year (April 2019 to May 2020) after filling the boxes with soil. The valve of the boxes was not opened to discharge during this period.
The cultivar of okra planted in this study was Abelmoschus esculentus L. cv. Wufu.The okra with five planting densities (4, 8, 12, 16 and 20 plants per box) was planted in May 2020 to obtain five groups of root densities, with a bare simulation box as a control (Table 2). Before sowing in the beginning of May, a basal fertilizer dose of N, P 2 O 5 and K 2 O were applied in each simulation box at the rates of 120, 100 and 50 kg·ha −1 , respectively. Each planting density was replicated 4 times. A total of 24 simulation boxes were placed on the flat land using a completely random design (CRD). During the seeding period of okra, supplemental freshwater (about 50 L) was pre-showered into the simulation box to ensure the germination of okra.

Observations and Analyses
The initial soil samples were taken before sowing in May 2020 while the last soil samples were taken at the beginning of November 2020. Soil was sampled by layers (20 cm) with an auger (5.0 cm diameter, 100 cm high) in the soil profile (up to 1 m depth). Three points were randomly selected in each simulation box to take samples and the soil in the same soil layer was mixed into a sample. The soil samples were air-dried, passed through a 1 mm sieve, and then water-soluble salts content of 1:5 soil-water extract was measured. The contents of cations including ions of Ca 2+ , Mg 2+ , K + , Na + and anions of Cl − , SO 4 2− , NO 3 − were determined by ion chromatography (ICS-2100, DIONEX, Sunnyvale, CA, USA). The contents of CO 3 2− and HCO 3 − were determined by potentiometric titration [28]. The salinity was calculated as the sum of all the above-measured ions. The average salinity value (Sa) in the soil profile (0-100 cm) of each group was also calculated: where Sp is the average salinity value of the soil profile (0-100 cm) in each simulation box, m is the number of soil layers, n is the number of replicates, S is the salinity of each soil layer in the simulation box. Sodium adsorption ratio (SAR) (mmol L) −1/2 was used to analyze the potential risk of soil alkalization. The formula for calculating SAR (mmol·L) −1/2 of effluent was as follows [29]: Here, the unit of concentration is mmol·L −1 .
The calculation method of average SAR in soil profile of each group was the same as average salinity value in soil profile of each group.
Soil desalination rate and SAR reduction rate were used to analyze the improvement effect in subsurface pipe drainage of coastal saline soil. The calculation formula was as follows: where DR is desalination rate (%), S1 is the salinity (g·kg −1 ) of soil which was sampled before okra was planted (before leaching), and S2 is the salinity (g·kg −1 ) of soil after okra harvest (after leaching). RR represents reduction rate (%) of SAR. SAR1 represents the SAR (mmol·L) −1/2 of soil which was sampled before okra was planted (before leaching), and SAR2 represents SAR (mmol·L) −1/2 of soil after okra harvest (after leaching). During the growth season of okra, high rainfall occurred four times. The drainage water was collected through opening the valve of the simulation box within 48 h after occurrence of above four rainfall events. Then, the volume of collected discharged water was measured by a cylinder. The cumulative drain outflow was the sum of four times of drain outflow (L).
After okra stopped growth in the beginning of November, the soil in the simulation box was dug out from every 10 cm layer and the roots of okra were washed out to determine the dry root weight by using the oven-dry method at 70 • C for 48 h. Root weight density (µg·cm −3 ) was used to describe and analyze the density and distribution of okra roots. The RWD of each simulation box was calculated using the following formula: DW (µg) is the total root dry mass of the sum of root dry mass in all soil layers within a simulation box while V (cm 3 ) is the volume of saline soil in the simulation box (2,700,000 cm 3 ). Using a vertical root distribution model to calculate root distribution coefficient (β, dimensionless), the calculation formula was as follows [30]: Y is the cumulative percentage (%) of root biomass with a depth of d (cm) from surface to the bottom of soil.
Some studies indicated that salt stress affects the development and growth of okra root system [31], and the heterogeneity of the salts present in the soil will increase the variation in okra RWD under the same planting density. In order to eliminate this effect, the RWD of 20 simulation boxes was arranged from small to large and the RWD of every 4 simulation boxes was divided into 5 groups (Figure 4). Each group from S1 to S5 was analyzed as the treatment object of this study ( Figure 5). The difference of RWD among five groups reached a significant level (p < 0.05). There were highly significant differences between S5 and S1, S2 and S3, and between S4 and S1 and S2 (p < 0.01).

Statistical Analyses
Paired t-test was used to evaluate the difference of salinity and SAR in soil sampled before and after leaching (p < 0.05). The normality of data was tested by Shapiro-Wilk method, and the homogeneity of data was tested based on the mean. One-way analysis of variance (ANOVA) and multiple comparisons (LSD test) were used to analyze the differences of RWD, the cumulative drain outflow, soil salt content, and SAR between different groups. Regression analysis was used to study the relationship between the cumulative drain outflow, the desalination rate, the decreasing rate of SAR, and the RWD of okra. A computer software SPSS 26.0 (SPSS Inc., Chicago, IL, USA) was used for statistical analysis. The software Origin Pro 2018 (Origin Lab Inc., Northampton, MA, USA) was used for graphical representation.

Root Vertical Distribution
The RWD significantly affected the vertical distribution characteristics of okra root dry mass in the soil profile ( Figure 6). The percentage of root dry mass in 0-20 cm soil layer decreased with the increase in RWD, but increased in 20-50 cm soil layer. For 40-50 cm soil depth, the percentage of root dry mass in all groups except for S5 was almost reduced to 0. Group S1 intersected with all other groups at the soil depth of 20 cm.
The roots of the five groups were fitted with the vertical distribution model of roots (Table 3). The curve of five groups fitted very well, for which R2 was 0.96 or above. RWD had great influence on root distribution coefficient (β) of okra, which decreased with the increase in RWD. Root distribution coefficient (β) increased to 115% of itself from 0.85 of S1 to 0.98 of S5. From S1 to S2, β increased slightly by 0.03. From S2 to S4, β increased very slightly by 0.01, and there was no significant change. From S4 to S5, β significantly increased by 0.09. Figure 6. Percentage of okra root dry mass for each soil layer in relation to root weight density.

The Drain Outflow
In the process of four times of drainage, the drain outflow between different groups had no significant difference due to the large difference within groups (Table 4). In the first and second drainage process, the drain outflow of S1 was the highest, while that of S5 was the lowest. The drain outflow of S1 and S5 was 189% and 192%, about 140% and 150% of the control, respectively. In the third drainage process, the drain outflow of S2 was the highest, while that of the CK was the least. The drain outflow of S2 was 242% of the control. In the fourth drainage process, the drain outflow of S3 was the highest, while that of the S1 was the least. The drain outflow of S3 was 161% of S5, and it was about 126% of the CK.
Regression analyses were performed on the RWD with the cumulative drain outflow (Figure 7). With the increase in RWD of okra, the cumulative drain outflow increased first and then decreased. The cumulative drain outflow of S2 was the highest and about 165% of that of S5. The fitting curve was a downward parabola and reached its maximum at the RWD of 116 µg·cm −3 , which was 140% of the cumulative drain outflow in CK. Except for group S5, the cumulative drain outflow of all groups was greater than that of the control.

Salt Content and Sodium Adsorption Ratio of Different Soil Layers
The salt content and SAR of soil sampled before and after leaching were analyzed ( Figure 8). The salinity of each soil layer before leaching ranged between 2.10-4.58 g·kg −1 , and the results of ANOVA in groups were not significant. Except for S2 and S5, the topsoil (0-20 cm) salt content of all groups was the highest among all soil layers. The salt contents of each soil layer after leaching ranged between 1.49-3.66 g·kg −1 , and the result of ANOVA was significant in salt contents of 40-60 cm, 60-80 cm, and 80-100 cm soil layers (p < 0.05). Except S1 and S5, the topsoil salt content of different groups decreased significantly (p < 0.05) after leaching. The SAR of each soil layer before leaching was 4.07-12.57 (mmol·L) −1/2 . The results of variance analysis among groups in the 60-80 cm soil layer reached a significant level (p < 0.05). The SAR of topsoil (0-20 cm) of all groups was the highest or a higher level in all soil layers. After leaching, the SAR of each soil layer ranged 4.63-9.09 (mmol·L) −1/2 , and the variance analysis results of SAR in the 0-20 cm and 80-100 cm soil layers were significant among the groups. Except S3, the SAR of topsoil and bottom soil layer in all groups was the lowest or a lower level in soil layers.
The desalination rate and SAR reduction rate of all groups in each soil layer were analyzed (Figure 9). In groups S2, S3, and S4, all soil layers at 0-100 cm soil depth were in a desalting state. The desalination rate of the top and sub-topsoil layer (20-40 cm) was not significantly different from that of the control, while the desalination rate of the middle soil layer (40-60 cm) and lower soil layers (60-100 cm) was significantly higher than that of the control. The overall desalting effect of S1, S5, and CK was poor, with one or more layers of salt accumulation after leaching.
The SAR of soil in S1, S2, S4, and S5 decreased effectively, with the SAR of all layers decreased after leaching during the rainy season. The reduction rate of SAR in topsoil at S2 was the highest among all groups. In the control group, the reduction rate of SAR in topsoil was higher, decreasing from the top layer to the bottom layer in turn. In S3, the reduction rate of SAR increased from the top layer to the lower layer, but it was negative in the top layer. The effects of okra RWD on soil desalination rate and reduction rate of SAR at different depths were analyzed ( Table 5). The regression analysis results showed that the soil desalination rate in all soil layers first increased and then decreased with the increase in okra RWD, showing a parabolic trend with downward opening. All values of R 2 of the regression analyses reached 0.59 and above, but none of them were significant. There was no obvious correlation between the reduction rate of SAR in each soil layer and the RWD of okra.

Average Salt Content and Sodium Adsorption Ratio in the Soil Profile
The average salt content and average SAR in the soil profile (0-100 cm) before and after leaching were analyzed ( Table 6). The results of variance analysis of soil salt content among groups were not significant. There were significant differences between the salt content of soil sampled before leaching and after leaching in the S3 and S4 groups (p < 0.05). Except S1, the desalination rate of all groups was higher than that of the control. The desalination rate of S2 was the highest, about 257% of the control. The soil desalination rate in groups from high to low was in the following order: S2, S3, S4, S5, CK, and S1. The soil SAR of different groups before leaching was 7.90-10.82 (mmol·L) −1/2 , and the results of variance analysis were significant (p < 0.05). After leaching, the soil SAR was 6.53-7.92 (mmol·L) −1/2 with no significant differences among groups. There were significant differences between SAR of soil sampled before and after leaching in S1 and S4 groups (p < 0.05). Except S3, the SAR reduction rate of groups was greater than that of the control. The reduction rate of SAR in S4 was the highest (161.81 µg·cm −3 ), which was about 180% of that of the control. The reduction rate of SAR from large to small was in the following order: S4, S2, S5, S1, CK and S3.
The effects of okra RWD on average soil desalination rate and average reduction rate of SAR in the soil profile were also analyzed ( Figure 10). The soil desalination rate first increased and then decreased with the increase in okra RWD. The fitting curve reached the maximum value (28.83%) when the RWD was 136 µg·cm −3 , which was about 229% of the control. There was no significant correlation between the reduction rate of soil SAR and okra RWD. Figure 10. The average soil desalination rate and average SAR reduction rate in relation to okra root weight density. CK represents the value corresponding to bare land. Data are the means of four replicates. SAR refers to sodium adsorption ratio.

Effect of Okra Root Weight Density on Root Vertical Distribution
Soil salinity, planting density, soil water, fertilizer distribution, and some other factors strongly affect the vertical distribution of roots in the soil profile [32][33][34]. The increase in planting density increases the root biomass, aggravates the competition among roots, and changes the vertical distribution of roots [35]. In the present study, the results showed that with the increase in RWD, more roots of okra grew into the deep soil. The percentage of root dry mass in the 0-20 cm soil layer decreased with the increase in RWD, while it increased with the increase in RWD in the 20-50 cm soil layer. These results were similar to those reported by an article studying the effect of planting density on the distribution of maize roots [36].
The fitting results of the root vertical distribution model also showed this phenomenon, as β decreased with the increase in RWD. However, β increased less at low RWD but more at high RWD, especially for the S4 to S5 groups, which reached 69% of the increase from group S1 to S5. This indicated that there was a threshold (161.81-181.07µg·cm −3 ) between S4 and S5. The okra RWD exceeded the threshold, the vertical distribution of okra roots was significantly affected, and more root penetrated beyond the surface layer to the deep soil layer. This may be one of the reasons why okra with different RWD has different effects on water and salt discharge performance in subsurface pipe.

Effect of Okra Root Weight Density on Drain Outflow
The drain outflow of subsurface pipe drainage was affected by many factors. Beside the depth and the spacing of the pipe drain, plants could influence the soil permeability [37], which affected the drainage capacity of subsurface pipe drainage. With the increase in root density, the more macropores were produced by roots, which further improved the structure and permeability of soil [38]. This would have a greater impact on the soil water and salt transport and discharge performance in subsurface pipe drainage, and would have increased the drain outflow from simulation boxes in the rainy season, but it could also increase the evapotranspiration per unit area of the soil at the same time [39], and reduce the soil water content and drain outflow.
The present study found that the cumulative drain outflow firstly increased and then decreased with the increase in RWD, and had a significant parabolic relationship with okra RWD. The fitting curve reached the maximum at the RWD about 116 µg·cm −3 , which was about 140% of that of the control. This indicated okra with an appropriate density (116 µg·cm −3 ) in coastal saline land could significantly increase the drain outflow of subsurface pipe drainage in the rainy season and further reduce ponding, surface runoff, and soil erosion caused by strong rainfall.
The results of this study were different from some studies. Their results showed that some forage grasses with deeper roots significantly reduced soil water content and groundwater level [40,41], which should lead to the reduction in the drain outflow. The possible reason was that the drain outflow in this experiment was measured within 48 h after heavy rainfall. Deep-rooted crops may reduce soil water content over a long period of time. However, the existing results of this study could not show the changes of soil water content over a long interval, so further research is needed in the future.

Effect of Okra Root Weight Density on Soil Desalination Rate and Reduction Rate of Sodium Adsorption Ratio
The results of leaching in saline soil were related to many factors such as initial soil salt content, soil permeability, water volume, water quality, and so on [42]. Plant roots can improve the permeability of soil and then influence the desalination rate [43].
In this study, the leaching volume of water and water quality were the same for all the RWD groups. The average initial soil salt content in the soil profile of each group was 2.33-4.09 g·kg −1 where all groups belong to mild saline soil and there was no significant difference in groups. However, the variation in salt content had some influence on the results of the experiment, so needs to be reduced in future research to avoid experimental error.
The present study found that okra RWD has a significant effect on the soil desalination rate. The soil desalination rate of average salt content in the soil profile first increased and then decreased with the increase in okra RWD. The fitting curve reached the maximum value (28.83%) when the RWD was 136 µg·cm −3 , which was about 229% of the control. In S2, S3, and S4, all soil layers at 0-100 cm soil depth were in a desalting state. The desalination rate of the top and sub-topsoil layer was not significantly different from that of the control, while the desalination rate of the middle soil layer (40-60 cm) and lower soil layers (60-100 cm) was significantly higher than that of the control. This means that when the RWD is 128-161 µg·cm −3 , the desalting effect of each soil layer is better than that of CK or as good as it is. The above results indicate that okra with an appropriate root (136 µg·cm −3 ) density in coastal saline land could significantly improve the desalting effect of subsurface pipe drainage in the rainy season and shorten the time needed for the improvement in coastal saline land.
The results of this study contradicted the views of an article, which suggested that growing crops increased leaching efficiency because it reduced soil water content, macropore bypass, and drainage volume [13]. However, in this study, the cumulative drain outflow of the group with the highest soil desalination rate was also very high. The possible reason was that the existence of the subsurface pipe drainage system generally increased the drainage rate, which led to higher drainage within a short time, and it was too late for crops to show the effect of reducing soil moisture.
The soil SAR reflected the relative content of Na + , Mg 2+ , and Ca 2+ in soil solution and can indicate the exchangeable sodium percentage to a certain extent [44]. In the present study, the soil SAR before leaching was 7.90-10.82 (mmol·L) −1/2 and the result of variance analysis was significant (p < 0.05). After leaching, it decreased to 6.53-7.92 (mmol·L) −1/2 with no significant difference among groups. It showed that subsurface pipe drainage can reduce the soil SAR after leaching.
The RWD of okra had a great effect on the reduction rate of SAR. Except S3, the reduction rate of average SAR in the soil profile in other groups was greater than the control. Group S4 was the highest, which was about 180% of the control. The above results indicate that okra with an appropriate root density in coastal saline soil could improve the ability of subsurface pipe drainage to inhibit soil alkalization that will change the relative content of ions and improve the physiochemical properties of the soil. However, the relationship between RWD and reduction rate of SAR needs a further study.
The results were similar to the results of an article on the improvement in saline soil by subsurface pipe drainage technology through field experiments. Its results showed that subsurface pipe drainage technology can effectively restrain the trend in soil alkalization, but that soil SAR increases after the rainy season [45]. The effect of an okra-assisted subsurface pipe drainage system on soil SAR needs to be further studied under field conditions.

Conclusions
The okra root density affected the vertical distribution of okra roots in soil layers. When the okra RWD exceeded a threshold, the vertical distribution of okra roots was significantly affected, and more root penetrated beyond the surface layer to the deep soil layers. The root density of okra influenced the water discharge performance of subsurface pipe drainage under natural precipitation. Okra with an appropriate RWD (about 116 µg·cm −3 ) in coastal saline soil could significantly increase the cumulative drain outflow of subsurface pipe drainage in the rainy season. The root density of okra influenced the salt discharge performance of subsurface pipe drainage under natural precipitation. In addition, okra with an appropriate RWD (about 136 µg·cm −3 ) in coastal saline land could significantly improve the desalting effect of subsurface pipe drainage in the rainy season. The root density of okra also influenced the ability of subsurface pipe drainage to inhibit soil alkalization. Planting okra and installing subsurface pipe drainage to control drainage at the coastal saline land in Northern China can effectively improve the water and salt drainage effect of subsurface pipe drainage technology when okra RWD is about 116-136 µg·cm −3 ; it will also enhance the ability of subsurface pipe drainage to inhibit soil alkalization, and shorten the time needed to improve the coastal saline soil. The results of this study provided a scientific basis for an okra-assisted subsurface pipe drainage system to reclaim coastal saline land in Northern China. When using subsurface pipe drainage to improve coastal saline soils, planting okra with proper density is an appropriate choice to improve the effect and benefit.

Data Availability Statement:
The authors confirm that the data supporting the findings of this study are available within the article.