ARTICLE   Open Access    

Salt tolerance of seven genotypes of zoysiagrass (Zoysia spp.)

More Information
  • Seven zoysiagrass genotypes were evaluated for salt tolerance in a greenhouse study. The plant materials included Zoysia matrella 'Diamond', Z. japonica 'Palisades', three Z. matrella × Z. japonica hybrids DALZ 1701, DALZ 1713, and 'Innovation', and two Z. minima × Z. matrella hybrids (DALZ 1309 and 'Lazer'). Treatments included a control (nutrient solution) and two saline treatments representing moderate and high salt levels. The electrical conductivity (EC) was 1.3 dS m−1 for control and moderate (EC5) and high salinity (EC10) were 5.0 and 10.0 dS m−1, respectively. At the end of eight-weeks of treatment, the relative (percent control) shoot dry weight (DW) was greatest in 'Diamond' in EC10, and the relative root DW was greatest in DALZ 1309 in EC5. A cluster analysis based on the relative tissue dry weight identified 'Diamond', DALZ 1309, and DALZ 1713 as the most salt tolerant genotypes. Additionally, the green leaf area (GLA) index of 'Diamond' and DALZ 1713 were 98.8% and 100%, respectively, indicating excellent visual appearance under high salt levels. Bi-weekly clipping DW showed that 'Diamond' continued to produce biomass throughout the duration of the study under the EC10 treatment. Sodium (Na+) and chloride (Cl) content in the shoot tissue of the seven turfgrass genotypes indicated that lower concentrations corresponded to greater salt tolerance indicating exclusion of Na+ and Cl from the shoot tissue. Taken together, the genotypes 'Diamond' and DALZ 1713 were determined to be the most salt tolerant and recommended for use in areas with high soil or water salinity.
  • In a recent report on Latin America's next petroleum boom, The Economist refers to the current and future situation in oil producing countries in the region. In the case of Argentina, the increase in oil and gas output 'have led to an increase in production in Vaca Muerta, a mammoth field in Argentina's far west. It holds the world's second-largest shale gas deposits and its fourth-largest shale oil reserves… Rystad Energy expects shell-oil production in Argentina will more than double by the end of the decade, to over a million barrels per day'[1].

    Oil production in Argentina is currently dominated by three Patagonian areas: Neuquén, San Jorge, and Austral. Based on 2021 information, 49% of oil reserves of Argentina are located in Neuquén, whereas San Jorge has 46%. Neuquén is also the largest source of oil (57%) and gas (37%) in the country. According to 2018 data, conventional oil produced in Argentina amounts to 87%, whereas non-conventional, shale production represents 13%; however, non-conventional oil is increasing due to Vaca Muerta shale oil exploitation.

    This increase of production in the Patagonian fields requires the use of a large fluid storage capacity by means of vertical oil storage tanks having different sizes and configurations. Tanks are required to store not just oil but also water. The exploitation of nonconventional reservoirs, such as Vaca Muerta, involves massive water storage to carry out hydraulic stimulation in low-permeability fields, and for managing the return fluid and production water at different stages of the process (storage, treatment, and final disposal).

    Storage tanks in the oil industry are large steel structures; they may have different sizes, and also different shell configurations, such as vertical cylinders with a fixed roof or with a floating roof and opened at the top[2]. It is now clear that such oil infrastructure is vulnerable to accidents caused by extreme weather events[35].

    Data from emergencies occurring in oil fields shows that accidents due to regional winds, with wind speed between 150 and 240 km/h, may cause severe tank damage. Seismic activity in the region, on the other hand, is of less concern to tank designers in Patagonia.

    Damage and failure mechanisms of these tanks largely depend on tank size and configuration, and their structural response should be considered from the perspective of shell mechanics and their consequences. In a report on damage observed in tanks following hurricanes Katrina and Rita in 2005[6,7], several types of damage were identified. The most common damage initiation process is due to shell buckling[811], which may progress into plasticity at higher wind speeds. In open-top tanks, a floating roof does not properly slide on a buckled cylindrical shell, and this situation may lead to different failure mechanisms. Further, damage and loss of integrity have the potential to induce oil spills, with direct consequences of soil contamination and also of fire initiation.

    Concern about an emergency caused by such wind-induced hazards involves several stakeholders, because the consequences may affect the operation of oil plants, the local and regional economies, the safety of the population living in the area of a refinery or storage farm, and the environment[6]. In view of the importance of preserving the shell integrity and avoiding tank damage, there is a need to evaluate risk of existing tanks at a regional level, such as in the Neuquén and San Jorge areas. This information may help decision makers in adopting strategies (such as structural reinforcement of tanks to withstand expected wind loads) or post-event actions (like damaged infrastructure repair or replacement).

    The studies leading to the evaluation of risk in the oil infrastructure are known as vulnerability studies, and the most common techniques currently used are fragility curves[12]. These curves evaluate the probability of reaching or exceeding a given damage level as a function of a load parameter (such as wind speed in this case).

    Early studies in the field of fragility of tanks were published[13] from post-event earthquake damage observations. Studies based on computational simulation of tank behavior under seismic loads were reported[14]. The Federal Emergency Management Administration in the US developed fragility curves for tanks under seismic loads for regions in the United States, and more recently, this has been extended to hurricane and flood events in coastal areas[15]. Seismic fragility in Europe has been reviewed by Pitilakis et al.[16], in which general concepts of fragility are discussed. Bernier & Padgett[17] evaluated the failure of tanks due to hurricane Harvey using data from aerial images and government databases. Fragility curves were developed based on finite element analyses and damage of the tank population was identified in the Houston Ship Channel. Flood and wind due to hurricane Harvey were also considered[18] to develop fragility curves.

    Because fragility curves for tanks under wind depend on the wind source (either hurricane or regional winds), and the type and size of tanks identified in a region, fragility curves developed for one area are not possible to be directly used in other areas under very different inventory and wind conditions.

    This paper addresses problems of shell buckling and loss of integrity of open top tanks, with wind-girders and floating roof and it focuses on the development of fragility curves as a way to estimate damage states under a given wind pressure. The region of interest in this work covers the oil producing areas in Patagonia, Argentina. Damage of tanks under several wind pressures are evaluated by finite element analyses together with methodologies to evaluate the structural stability.

    The construction of fragility curves requires information from the following areas: First, an inventory of tanks to be included in the analysis; second, data about the loads in the region considered; third, data about structural damage, either observed or computed via modeling; and fourth, a statistical model that links damage and load/structure data. This section describes the main features of the tank population considered in the study.

    The construction of an inventory at a regional level is a very complex task, which is largely due to a lack of cooperation from oil companies to share information about their infrastructure. Thus, to understand the type of tanks in an oil producing region, one is left collecting a limited number of structural drawings and aerial photography. A detailed inventory of the Houston Ship Channel was carried out by Bernier et al.[19], who identified 390 floating roof tanks. An inventory for Puerto Rico[20] identified 82 floating roof tanks. Although both inventories used different methodologies and addressed very different tank populations, some common features were found in both cases.

    An alternative strategy to carry out fragility studies is to develop a database using a small number of tanks, for which a detailed structural behavior is investigated using finite element analysis. This is a time-consuming task, but it allows identification of buckling pressures, buckling modes, and shell plasticity. This information serves to build approximate fragility curves, and it can also be used to develop what are known as meta-models, which predict structural damage based on tank/load characteristics. Such meta-models take the form of equations that include the tank geometry and wind speed to estimate damage. Meta-models were used, for example, in the work of Kameshwar & Padgett[18].

    This work employs a simplified strategy, and addresses the first part of the procedure described above. The use of a limited number of tanks in a database, for which a finite element structural analysis is carried out. This leads to fragility curves based on a simplified tank population (reported in this work) and the development of a meta-model together with enhanced fragility results will be reported and compared in a future work.

    Partial information of tanks in the Patagonian region was obtained from government sources, and this was supplemented by aerial photography showing details of tank farms in the region. As a result of that, it was possible to establish ranges of tank dimensions from which an artificial database was constructed.

    The present study is restricted to open-top tanks with a wind girder at the top. They are assumed to have floating roofs, which are designed and fabricated to allow the normal operation of the roof without the need of human intervention. The main characteristics of tanks investigated in this paper, are illustrated in Fig. 1.

    Figure 1.  Geometric characteristics of open-topped oil storage considered in this paper.

    The range of interest in terms of tank diameter D was established between 35 m < D < 60 m. Based on observation of tanks in the region, the ratios D/H were found to be in the range 0.20 < D/H < 0.60, leading to cylinder height H in the range 12 m < H < 20 m. These tanks were next designed using API 650[21] regulations to compute their shell thickness and wind girder dimensions. A variable thickness was adopted in elevation, assuming 3 m height shell courses. The geometries considered are listed in Table 1, with a total of 30 tanks having combinations of five values of H and six values of D. The volume of these tanks range between 55,640 and 272,520 m3.

    Table 1.  Geometry and course thickness of 30 tanks considered in this work.
    H
    (m)
    CoursesThickness t (m)
    D = 35 mD = 40 mD = 45 mD = 50 mD = 55 mD = 60 m
    12V10.0140.0160.0180.0180.0200.022
    V20.0120.0120.0140.0160.0160.018
    V30.0080.0100.0100.0100.0120.012
    V40.0060.0080.0080.0080.0080.008
    14V10.0160.0180.0200.0220.0250.025
    V20.0140.0140.0160.0180.0200.020
    V30.0100.0120.0120.0140.0140.016
    V40.0080.0080.0080.0100.0100.010
    V50.0060.0080.0080.0080.0080.008
    16V10.0180.0200.0220.0250.0280.028
    V20.0160.0180.0180.0200.0220.025
    V30.0120.0140.0160.0160.0180.020
    V40.0100.0100.0120.0120.0140.014
    V50.0060.0080.0080.0080.0080.010
    V60.0060.0080.0080.0080.0080.010
    18V10.0200.0220.0250.0280.0300.032
    V20.0180.0200.0220.0250.0250.028
    V30.0140.0160.0180.0200.0200.022
    V40.0120.0120.0140.0160.0160.018
    V50.0080.0100.0100.0100.0120.012
    V60.0080.0100.0100.0100.0120.012
    20V10.0220.0250.0280.0300.0320.035
    V20.0200.0220.0250.0280.0280.032
    V30.0160.0180.0200.0220.0250.028
    V40.0140.0140.0160.0180.0200.020
    V50.0100.0120.0120.0140.0140.016
    V60.0100.0120.0120.0140.0140.016
    V70.0100.0120.0120.0140.0140.016
     | Show Table
    DownLoad: CSV

    The material assumed in the computations was A36 steel, with modulus of elasticity E = 201 GPa and Poisson's ratio ν = 0.3.

    For each tank, a ring stiffener was designed as established by API 650[21], in order to prevent buckling modes at the top of the tank. The minimum modulus Z to avoid ovalization at the top of the tank is given by

    Z=D2H17(V190)2 (1)

    where V is the wind speed, in this case taken as V = 172.8 km/h for the Patagonian region. Intermediate ring stiffeners were not observed in oil tanks in Patagonia, so they were not included in the present inventory.

    Because a large number of tanks need to be investigated in fragility studies, it is customary to accept some simplifications in modeling the structure to reduce the computational effort. The geometry of a typical ring stiffener at the top is shown in Fig. 2a, as designed by API 650. A simplified version was included in this research in the finite element model, in which the ring stiffener is replaced by an equivalent thickness at the top, as suggested in API Standard 650[21]. This approach has been followed by most researchers in the field. The equivalent model is shown in Fig. 2b.

    Figure 2.  Ring stiffener, (a) design according to API 650, (b) equivalent section[22].

    The pressure distribution due to wind around a short cylindrical shell has been investigated in the past using wind tunnels and computational fluid dynamics, and a summary of results has been included in design regulations.

    There is a vast number of investigations on the pressures in storage tanks due to wind, even if one is limited to isolated tanks, as in the present paper. For a summary of results, see, for example, Godoy[11], and only a couple of studies are mentioned here to illustrate the type of research carried out in various countries. Wind tunnel tests were performed in Australia[23], which have been the basis of most subsequent studies. Recent tests in Japan on small scale open top tanks were reported[24,25]. In China, Lin & Zhao[26] reported tests on fixed roof tanks. CFD models, on the other hand, were computed[27] for open top tanks with an internal floating roof under wind flow. Although there are differences between pressures obtained in different wind tunnels, the results show an overall agreement.

    The largest positive pressures occur in the windward meridian covering an angle between 30° and 45° from windward. Negative pressures (suction), on the other hand, reach a maximum at meridians located between 80° and 90° from windward. An evaluation of US and European design recommendations has been reported[28,29], who also considered the influence of fuel stored in the tank.

    The circumferential variation of pressures is usually written in terms of a cosine Fourier series. The present authors adopted the series coefficients proposed by ASCE regulations[30], following the analytical expression:

    q=λinCicos(iφ) (2)

    in which λ is the amplification factor; the angle φ is measured from the windward meridian; and coefficients Ci represent the contribution of each term in the series. The following coefficients were adopted in this work (ASCE): C0 = −0.2765, C1 = 0.3419, C2 = 0.5418, C3 = 0.3872, C4 = 0.0525, C5 = 0.0771, C6 = −0.0039 and C7 = 0.0341. For short tanks, such as those considered in this paper, previous research reported[31] that for D/H = 0.5 the variation of the pressure coefficients in elevation is small and may be neglected to simplify computations. Thus, the present work assumes a uniform pressure distribution in elevation at each shell meridian.

    In fragility studies, wind speed, rather than wind pressures, are considered, so that the following relation from ASCE is adopted in this work:

    qz=0.613KztKdV2IV=qz0.613KztKdI (3)

    in which I is the importance factor; Kd is the directionality factor; and Kzt is the topographic factor. Values of I = 1.15, Kd = 0.95 and Kzt = 1, were adopted for the computations reported in this paper.

    Because shell buckling was primarily investigated in this work using a bifurcation analysis, the scalar λ was increased in the analysis until the finite element analysis detected a singularity.

    Fragility curves are functions that describe the probability of failure of a structural system (oil tanks in the present case) for a range of loads (wind pressures) to which the system could be exposed. In cases with low uncertainty in the structural capacity and acting loads, fragility curves take the form of a step-function showing a sudden jump (see Fig. 3a). Zero probability occurs before the jump and probability equals to one is assumed after the jump. But in most cases, in which there is uncertainty about the structural capacity to withstand the load, fragility curves form an 'S' shape, as shown in Fig. 3a and b probabilistic study is required to evaluate fragility.

    Figure 3.  Examples of fragility curves, (a) step-function, (b) 'S' shape function.

    The construction of fragility curves is often achieved by use of a log-normal distribution. In this case, the probability of reaching a certain damage level is obtained by use of an exponential function applied to a variable having a normal distribution with mean value μ and standard deviation σ. If a variable x follows a log-normal distribution, then the variable log(x) has a normal distribution, with the following properties:

    • For x < 0, a probability equal to 0 is assigned. Thus, the probability of failure for this range is zero.

    • It can be used for variables that are computed by means of a number of random variables.

    • The expected value in a log-normal distribution is higher than its mean value, thus assigning more importance to large values of failure rates than would be obtained in a normal distribution.

    The probability density function for a log-normal distribution may be written in the form[32]:

    f(xi)=12πσ21xexp[(lnxµ)2/(2σ2)] (4)

    in which f(xi) depends on the load level considered, and is evaluated for a range of interest of variable x; and μ* is the mean value of the logarithm of variable x associated with each damage level. Damage levels xi are given by Eqn (5).

    µ(xi)=1NNn=1ln(xin) (5)

    where the mean value is computed for a damage level xi, corresponding to I = DSi; summation in n extends to the number of tanks considered in the computation of the mean value. Damage levels in this work are evaluated using computational modeling and are defined in the next section. Variance is the discrete variable xi 2), computed from:

    σ2(xi,µ)=1NNn=1(ln(xin)µ)2=1NNn=1ln(xin)2µ2 (6)

    The probability of reaching or exceeding a damage level DSi is computed by the integral of the density function using Eqn (7), for a load level considered (the wind speed in this case):

    P[DS/x]=x=V0x=0f(x)dx (7)

    where V0 is the wind speed at which computations are carried out, and x is represented by wind speed V.

    Various forms of structural damage may occur as a consequence of wind loads, including elastic or plastic deflections, causing deviations from the initial perfect geometry; crack initiation or crack extension; localized or extended plastic material behavior; and structural collapse under extreme conditions. For the tanks considered in this work, there are also operational consequences of structural damage, such as blocking of a floating roof due to buckling under wind loads that are much lower than the collapse load. For this reason, a damage study is interested in several structural consequences but also in questions of normal operation of the infrastructure. Several authors pointed out that there is no direct relation between structural damage and economic losses caused by an interruption of normal operation of the infrastructure.

    Types of damage are usually identified through reconnaissance post-event missions, for example following Hurricanes Katrina and Rita[6,7]. Damage states reported in Godoy[7] include shell buckling, roof buckling, loss of thermal insulation, tank displacement as a rigid body, and failure of tank/pipe connections. These are qualitative studies, in which damage states previously reported in other events are identified and new damage mechanisms are of great interest in order to understand damage and failure modes not taken into account by current design codes.

    In this work, in which interest is restricted to open top tanks having a wind girder at the top, four damage states were explored, as shown in Table 2. Regarding the loss of functionality of a tank, several conditions may occur: (1) No consequences for the normal operation of a tank; (2) Partial loss of operation capacity; (3) Complete loss of operation.

    Table 2.  Damage states under wind for open-top tanks with a wind girder.
    Damage states (DS)Description
    DS0No damage
    DS1Large deflections on the cylindrical shell
    DS2Buckling of the cylindrical shell
    DS3Large deflections on the stiffening ring
     | Show Table
    DownLoad: CSV

    DS1 involves displacements in some area of the cylindrical body of the tank, and this may block the free vertical displacement of the floating roof. Notice that this part of the tank operation is vital to prevent the accumulation of inflammable gases on top of the fluid stored. Blocking of the floating roof may cause a separation between the fuel and the floating roof, which in turn may be the initial cause of fire or explosion.

    DS2 is associated with large shell deflections, which may cause failure of pipe/tank connections. High local stresses may also arise in the support of helicoidal ladders or inspection doors, with the possibility of having oil spills.

    DS3 is identified for a loss of circularity of the wind girder. The consequences include new deflections being transferred to the cylindrical shell in the form of geometrical imperfections.

    In summary, DS1 and DS3 may affect the normal operation of a floating roof due to large shell or wind-girder deflections caused by buckling.

    Tank modeling was carried out in this work using a finite element discretization within the ABAQUS environment[33] using rectangular elements with quadratic interpolation functions and reduced integration (S8R5 in the ABAQUS nomenclature). Two types of shell analysis were performed: Linear Bifurcation Analysis (LBA), and Geometrically Nonlinear Analysis with Imperfections (GNIA). The tank perimeter was divided into equal 0.35 m segments, leading to between 315 and 550 elements around the circumference, depending on tank size. Convergence studies were performed and errors in LBA eigenvalues were found to be lower than 0.1%.

    The aim of an LBA study is to identify a first critical buckling state and buckling mode by means of an eigenvalue problem. The following expression is employed:

    (K0+λCKG)ΦC=0 (8)

    where K0 is the linear stiffness matrix of the system; KG is the load-geometry matrix, which includes the non-linear terms of the kinematic relations; λC is the eigenvalue (buckling load); and ΦC is the critical mode (eigenvector). For a reference wind state, λ is a scalar parameter. One of the consequences of shell buckling is that geometric deviations from a perfect geometry are introduced in the shell, so that, due to imperfection sensitivity, there is a reduced shell capacity for any future events.

    The aim of the GNIA study is to follow a static (non-linear) equilibrium path for increasing load levels. The GNIA study is implemented in this work using the Riks method[34,35], which can follow paths in which the load or the displacement decrease. The geometric imperfection was assumed with the shape of the first eigenvector at the critical state in the LBA study, and the amplitude of the imperfection was defined by means of a scalar ξ [10]. To illustrate this amplitude, for a tank with D = 45 m and H = 12 m, the amplitude of imperfection is equal to half the minimum shell thickness (ξ = 4 mm in this case).

    It was assumed that a damage level DS1 is reached when the displacement amplitudes do not allow the free vertical displacement of the floating roof. Based on information from tanks in the Patagonian region, the limit displacement was taken as 10 mm. This state was detected by GNIA, and the associated load level is identified as λ = λDS1.

    The load at which damage state DS2 occurs was obtained by LBA, leading to a critical load factor λC and a buckling mode. An example of damage levels is shown in Fig. 4.

    Figure 4.  Damage computed for a tank with D = 45 m and H = 12 m. (a) Deflected shape for damage DS1; (b) Equilibrium path for node A (DS1); (c) Deflected shape for damage DS2 (critical mode).

    An LBA study does not account for geometric imperfections. It is well known that the elastic buckling of shells is sensitive to imperfections, so that a reduction in the order of 20% should be expected for cylindrical shells under lateral pressure. This consideration allows to estimate DS0 (a state without damage) as a lower bound of the LBA study. An approach to establish lower bounds for steel storage tanks is the Reduced Stiffness Method (RSM)[3640]. Results for tanks using the RSM to estimate safe loads show that λDS0 = 0.5λDS2 provides a conservative estimate for present purposes.

    DS3 was computed using a linear elastic analysis to evaluate the wind pressure at which a 10 mm displacement of the wind girder is obtained.

    In a similar study for tanks with a fixed conical roof, Muñoz et al.[41] estimated a collapse load based on plastic behavior. However, in the present case the top ring has a significant stiffness, and this leads to extremely high wind speeds before reaching collapse (higher than 500 km/h). For this reason, the most severe damage level considered here was that of excessive out-of-plane displacements of the wind girder, and not shell collapse.

    The methodology to construct fragility curves has been presented by several authors[42,43]. The following procedure was adapted here[44]: (1) Establish qualitative damage categories (Table 2). (2) Compute a data base for different tanks, using LBA and GNIA. In each case, the damage category based on step (1) was identified (Table 3). (3) Approximate data obtained from step (2) using a log-normal distribution. (4) Plot the probabilities computed in step (3) with respect to wind speed x.

    Table 3.  Wind speed for each tank considered reaching a damage level.
    HDIDDS0DS1DS2DS3
    H12D351137.76162.06194.82336.02
    D402160.62181.31227.16360.73
    D453153.32174.19216.82374.04
    D504145.23165.27205.39373.76
    D555152.76180.83216.03374.75
    D606145.11170.75205.22370.98
    H14D357145.57162.05205.87295.03
    D408148.55166.20210.08311.24
    D459136.42153.72192.92334.54
    D5010155.36177.51219.71339.86
    D5511145.24165.34205.39343,17
    D6012141.89167.26200.67338.77
    H16D3513131.32161.94185.71262.20
    D4014146.95163.99207.82277.08
    D4515150.58170.90212.95293.37
    D5016138.97161.05196.54303.62
    D5517138.51174.17195.88313.97
    D6018156.34182.78221.10326.83
    H18D3519146.80160.79207.60223.18
    D4020159.01177.71224.87243.63
    D4521157.10179.51222.17265.32
    D5022152.54172.17215.72293.32
    D5523164.93188.10233.25305.94
    D6024163.69180.32231.49315.63
    H20D3525163.64199.59231.42195.03
    D4026171.24195.14242.18216.47
    D4527171.58203.68242.64293.32
    D5028182.46209.43258.03259.41
    D5529178.95208.23253.07272.48
    D6030174.47196.11246.74290.86
     | Show Table
    DownLoad: CSV

    Wind speeds for each tank, obtained via Eqn (3), are shown in Table 3 for the pressure level associated with each damage level DSi. A scalar ID was included in the table to identify each tank of the population in the random selection process. Wind speed was also taken as a random variable, so that wind speed in the range between 130 and 350 km/h have been considered at 5 km/h increase, with intervals of −2.5 and +2.5 km/h.

    Out of the 30-tank population considered, a sample of 15 tanks were chosen at random and were subjected to random wind forces. The random algorithm allowed for the same tank geometry to be chosen more than once as part of the sample.

    The type of damage obtained in each case for wind speed lower or equal to the upper bound of the interval were identified. Table 4 shows a random selection of tanks, together with the wind speed required to reach each damage level. For example, for a wind speed of 165 km/h, the wind interval is [162.5 km/h, 167.5 km/h]. This allows computation of a damage matrix (shown in Table 5). A value 1 indicates that a damage level was reached, whereas a value 0 shows that a damage level was not reached. In this example, 13 tanks reached DS0; six tanks reached DS1; and there were no tanks reaching DS2 or DS3. The ratio between the number of tanks with a given damage DSi and the total number of tanks selected is h, the relative accumulated frequency. The process was repeated for each wind speed and tank selection considered.

    Table 4.  Random tank selection for V = 165 km/h, assuming wind interval [162.5 km/h, 167.5 km/h].
    IDDS0DS1DS2DS3
    11145.2165.3205.4343.2
    6145.1170.7205.2371.0
    3153.3174.2216.8374.0
    9136.4153.7192.9334.5
    28182.5209.4258.0259.4
    22152.5172.2215.7293.3
    13131.3161.9185.7262.2
    19146.8160.8207.6223.2
    3153.3174.2216.8374.0
    12141.9167.3200.7338.8
    30174.5196.1246.7290.9
    23164.9188.1233.2305.9
    2160.6181.3227.2360.7
    17138.5174.2195.9314.0
    11145.2165.3205.4343.2
     | Show Table
    DownLoad: CSV
    Table 5.  Damage matrix for random tank selection (V = 165 km/h), assuming wind interval [162.5 km/h, 167.5 km/h].
    DS0DS1DS2DS3
    1100
    1000
    1000
    1100
    0000
    1000
    1100
    1100
    1000
    1100
    0000
    1000
    1000
    1000
    1100
    Total13600
    hi0.870.400
     | Show Table
    DownLoad: CSV

    Table 6 shows the evaluation of the fragility curve for damage level DS0. This requires obtaining the number of tanks for each wind speed (fi), the cumulative number as wind speed is increased (Fi), and the frequency with respect to the total number of the sample of 15 tanks is written on the right-hand side of Table 6, for relative frequency (hi) and accumulated frequency (Hi).

    Table 6.  Damage DS0: Wind speed intervals [km/h] shown on the left; logarithm of wind speed; and relative and absolute frequencies (shown on the right).
    V inf
    (km/h)
    V m
    (km/h)
    V sup
    (km/h)
    Ln
    (Vm)
    fiFihiHi
    127.5130132.54.87000.0000
    132.5135137.54.91220.1330.133
    137.5140142.54.94130.0670.200
    142.5145147.54.98360.2000.400
    147.5150152.55.01170.0670.467
    152.5155157.55.044110.2670.733
    157.5160162.55.080110.0000.733
    162.5165167.55.112130.1330.867
    167.5170172.55.140130.0000.867
    172.5175177.55.160130.0000.867
    177.5180182.55.192150.1331.000
     | Show Table
    DownLoad: CSV

    With the values of mean and deviation computed with Eqns (5) & (6), it is possible to establish the log normal distribution of variable V for damage level DS0, usually denoted as P[DS0/V]. Values obtained in discrete form and the log-normal distribution are shown in Fig. 5a for DS0. For the selection shown in Table 6, the media is μ* = 5.03 and the deviation is σ = 0.09.

    Figure 5.  Probability of reaching a damage level P[DSi/V], (a) DS0, (b) DS0, DS1, DS2 and DS3.

    The process is repeated for each damage level to obtain fragility curves for DS1, DS2, and DS3 (Fig. 5b). Notice that the wind speeds required to reach DS3 are much higher than those obtained for the other damage levels. Such values should be compared with the regional wind speeds in Patagonia, and this is done in the next section.

    The oil producing regions in Argentina having the largest oil reserves are the Neuquén and the San Jorge regions, both located in Patagonia. This needs to be placed side by side with wind loads to understand the risk associated with such oil production.

    Figure 6 shows the geographical location of these regions. The Neuquén region includes large areas of four provinces in Argentina (Neuquén, south of Mendoza, west of La Pampa, and Río Negro). The San Jorge region is in the central Patagonia area, including two provinces (south of Chubut, north of Santa Cruz). Another area is the Austral region covering part of a Patagonian province (Santa Cruz).

    Figure 6.  Oil producing regions in Argentina. (Adapted from IAPG[47]).

    A map of basic wind speed for Argentina is available in the Argentinian code CIRSOC 102[45], which is shown in Fig. 7. Notice that the highest wind speeds are found in Patagonia, and affect the oil-producing regions mentioned in this work. For the Neuquén region, wind speeds range from 42 to 48 m/s (151.2 to 172.8 km/h), whereas for San Jorge Gulf region they range between 52 and 66 m/s (187.2 and 237.6 km/h).

    Figure 7.  Wind speed map of Argentina. (Adapted from CIRSOC 102[45]).

    The wind values provided by CIRSOC 102[45] were next used to estimate potential shell damage due to wind. Considering the fragility curves presented in Fig. 4, for damage levels DS0, DS1, DS2 and DS3 based on a log-normal distribution, it may be seen that it would be possible to have some form of damage in tanks located in almost any region of Argentina because CIRSOC specifies wind speeds higher than 36 m/s (129.6 km/h). The fragility curve DS0 represents the onset of damage for wind speeds higher than 130 km/h, so that only winds lower than that would not cause tank damage.

    Based on the fragility curves shown in Fig. 8, it is possible to estimate probable damage levels for the wind speed defined by CIRSOC. Because design winds in Patagonia are higher than 165.6 km/h (46 m/s), it is possible to conclude that there is 81% probability to reach DS0 and 25% to reach DS1.

    Figure 8.  Probability P[DSi/V] to reach damage levels DS1, DS2 and DS3 in tanks located in the Patagonia region of Argentina.

    For the geographical area of the Neuquén region in Fig. 6, together with the wind map of Fig. 7, the expected winds range from 150 to 172.8 km/h (42 to 48 m/s). Such wind range is associated with a DS0 probability between 41% and 92%, whereas the DS1 probability is in the order of 48%.

    A similar analysis was carried out for the San Jorge region, in which winds between 187.2 and 237 km/h (52 and 66 m/s). The probability of reaching DS1 is 87%, and the probability of DS2 is 88%. Wind girder damage DS3 could only occur in this region, with a lower probability of 18%.

    This work focuses on open top tanks having a floating roof, and explores the probability of reaching damage levels for wind loads, using the methodology of fragility curves. A population of 30 tanks was defined with H/D ratios between 0.2 and 0.6; such aspect ratios were found to be the most common in the oil producing regions of Patagonia. The data employed assumed diameters D between 35 and 60 m, together with height between 12 and 20 m. The tanks were designed using current API 650 regulations which are used in the region, in order to define the shell thickness and wind girder. All tanks were assumed to be empty, which is the worst condition for shell stability because a fluid stored in a tank has a stabilizing effect and causes the buckling load to be higher.

    Both structural damage (shell buckling) and operational damage (blocking of the floating roof due to deflections of the cylindrical shell) were considered in the analysis. The qualitative definition of damage levels in this work was as follows: The condition of no damage was obtained from a lower bound of buckling loads. This accounts for geometric imperfections and mode coupling of the shell. Shell buckling was evaluated using linear bifurcation analysis to identify damage level DS2. A geometrically non-linear analysis with imperfections was used to identify deflection levels that would block a floating roof, a damage level identified as DS1. Finally, deflections in the wind girder were investigated using a linear elastic analysis to define damage DS3.

    The present results were compared with the wind conditions of Patagonia, to show that several damage levels may occur as a consequence of wind speeds higher than 130 km/h, which is the expected base value identified for the region. The most frequent expected damage is due to the loss of vertical displacements of the floating roof due to large displacements in the cylindrical shell of the tank, and this may occur for wind speed up to 200 km/h. Damage caused by shell buckling may occur for wind speeds higher than 190 km/h, and for that wind speed, further damage due to displacements in the wind girder may also occur, but with a lower probability. This latter damage form requires much higher wind speed to reach a probability of 20%, and would be more representative of regions subjected to hurricanes.

    The number of tanks considered in the present analysis was relatively low, mainly because the aim of this work was to collect data to build a meta-model, i.e. a simple model that may estimate damage based on shell and load characteristics[46]. In future work, the authors expect to develop and apply such meta-models to a larger number of tank/wind configurations, in order to obtain more reliable fragility curves.

    Fragility studies for an oil producing region, like those reported in this work, may be important to several stakeholders in this problem. The fragility information links wind speed levels to expected infrastructure damage, and may be of great use to government agencies, engineering companies, and society at large, regarding the risk associated with regional oil facilities. At a government level, this helps decision makers in allocating funding to address potential oil-related emergencies cause by wind. This can also serve as a guide to develop further modifications of design codes relevant to the oil infrastructure. The engineering consequences may emphasize the need to strengthen the present regional infrastructure to reduce risk of structural damage and its consequences. The impact of damage in the oil infrastructure on society was illustrated in the case of Hurricane Katrina in 2005, in which a large number of residents had to be relocated due to the conditions created by the consequences of infrastructure failure.

    The authors confirm contribution to the paper as follows: study conception and design: Jaca RC, Godoy LA; data collection: Grill J, Pareti N; analysis and interpretation of results: Jaca RC, Bramardi S, Godoy LA; draft manuscript preparation: Jaca RC, Godoy LA. All authors reviewed the results and approved the final version of the manuscript.

    All data generated or analyzed during this study are included in this published article.

    The authors are thankful for the support of a grant received from the National Agency for the Promotion of Research, Technological Development and Innovation of Argentina and the YPF Foundation. Luis A. Godoy thanks Prof. Ali Saffar (University of Puerto Rico at Mayaguez) for introducing him to the field of fragility studies.

  • The authors declare that they have no conflict of interest.

  • Supplemental Fig. S1 Photos of the seven genotypes that were treated for 8 weeks: control (left), EC5 (middle), and EC10 (right). The genotype arrangement in each treatment was randomized. Photos were taken right after trimming.
  • [1]

    Chavarria M, Wherley B, Thomas J, Chandra A, Raymer P. 2019. Salinity tolerance and recovery attributes of warm-season turfgrass cultivars. HortScience 54:1625−31

    doi: 10.21273/hortsci13963-19

    CrossRef   Google Scholar

    [2]

    Chaudhary J. 2019. Application of reclaimed water for irrigation: A review. (Accessed 9 Nov. 2022) https://soils.ifas.ufl.edu/media/soilsifasufledu/sws-main-site/pdf/technical-papers/Chaudhary_Juhi_Immediate_Release.pdf

    [3]

    United States Environmental Protection Agency (EPA). 2012. Guidelines for water reuse. EPA/600/R-12/618. (Accessed 10 Oct. 2022) https://www.epa.gov/sites/default/files/2019-08/documents/2012-guidelines-water-reuse.pdf

    [4]

    Chen W, Lu S, Jiao W, Wang M, Chang AC. 2013. Reclaimed water: A safe irrigation water source. Environmental Development 8:74−83

    doi: 10.1016/j.envdev.2013.04.003

    CrossRef   Google Scholar

    [5]

    Wu S, Sun Y, Niu G, Altland J, Cabrera R. 2016. Response of 10 aster species to saline water irrigation. HortScience 51:197−201

    doi: 10.21273/HORTSCI.51.2.197

    CrossRef   Google Scholar

    [6]

    Qian YL, Fu JM, Wilhelm SJ, Christensen D, Koski AJ. 2007. Relative salinity tolerance of turf-type saltgrass selections. HortScience 42:205−9

    doi: 10.21273/HORTSCI.42.2.205

    CrossRef   Google Scholar

    [7]

    Kronzucker HJ, Britto DT. 2011. Sodium transport in plants: a critical review. New Phytologist 189:54−81

    doi: 10.1111/j.1469-8137.2010.03540.x

    CrossRef   Google Scholar

    [8]

    Parida AK, Das AB. 2005. Salt tolerance and salinity effects on plants: a review. Ecotoxicology and Environmental Safety 60:324−49

    doi: 10.1016/j.ecoenv.2004.06.010

    CrossRef   Google Scholar

    [9]

    Flowers TJ, Munns R, Colmer TD. 2015. Sodium chloride toxicity and the cellular basis of salt tolerance in halophytes. Annals of Botany 115:419−31

    doi: 10.1093/aob/mcu217

    CrossRef   Google Scholar

    [10]

    Munns R, Tester M. 2008. Mechanisms of salinity tolerance. Annual Review of Plant Biology 59:651−81

    doi: 10.1146/annurev.arplant.59.032607.092911

    CrossRef   Google Scholar

    [11]

    Bernstein L, François LE, Clark RA. 1972. Salt tolerance of ornamental shrubs and ground covers. Journal of the American Society for Horticultural Science 97:550−56

    doi: 10.21273/JASHS.97.4.550

    CrossRef   Google Scholar

    [12]

    Niu G, Cabrera RI. 2010. Growth and physiological responses of landscape plants to saline water irrigation: A review. HortScience 45:1605−9

    doi: 10.21273/HORTSCI.45.11.1605

    CrossRef   Google Scholar

    [13]

    Marcum KB, Pessarakli M, Kopec DM. 2005. Relative salinity tolerance of 21 turf-type desert saltgrasses compared to bermudagrass. HortScience 40:827−29

    doi: 10.21273/HORTSCI.40.3.827

    CrossRef   Google Scholar

    [14]

    Wu Y, Chandra A, Genovesi AD, Meeks M. 2017. Freezing tolerance of field-and growth chamber-acclimated zoysiagrass hybrids. International Turfgrass Society Research Journal 13:421−26

    doi: 10.2134/itsrj2016.05.0427

    CrossRef   Google Scholar

    [15]

    Patton AJ, Schwartz BM, Kenworthy KE. 2017. Zoysiagrass (Zoysia spp.) history, utilization, and improvement in the United States: A review. Crop Science 57:S37−S72

    doi: 10.2135/cropsci2017.02.0074

    CrossRef   Google Scholar

    [16]

    Meeks M, Kong ST, Genovesi AD, Smith B, Chandra A. 2022. Low-input golf course putting green performance of fine-textured inter- and intra-specific zoysiagrass (Zoysia spp.). International Turfgrass Society Research Journal 14:610−21

    doi: 10.1002/its2.63

    CrossRef   Google Scholar

    [17]

    Marcum KB, Anderson SJ, Engelke MC. 1998. Salt gland ion secretion: A salinity tolerance mechanism among five zoysiagrass species. Crop Science 38:806−10

    doi: 10.2135/cropsci1998.0011183X003800030031x

    CrossRef   Google Scholar

    [18]

    Maeda Y. 2019. Effects of calcium application on the salt tolerance and sodium excretion from salt glands in zoysiagrass (Zoysia japonica). Grassland Science 65:189−96

    doi: 10.1111/grs.12234

    CrossRef   Google Scholar

    [19]

    Marcum KB, Murdoch CL. 1994. Salinity tolerance mechanisms of six C4 turfgrasses. Journal of the American Society for Horticultural Science 119:779−84

    doi: 10.21273/JASHS.119.4.779

    CrossRef   Google Scholar

    [20]

    Qian YL, Engelke MC, Foster MJV. 2000. Salinity effects on zoysiagrass cultivars and experimental lines. Crop Science 40:488−92

    doi: 10.2135/cropsci2000.402488x

    CrossRef   Google Scholar

    [21]

    Cavins TJ, Whipker BE, Fonteno WC. 2008. PourThru: A method for monitoring nutrition in the greenhouse. Acta Horticulturae 779:289−98

    doi: 10.17660/actahortic.2008.779.35

    CrossRef   Google Scholar

    [22]

    Havlin JL, Soltanpour PN. 1980. A nitric acid plant tissue digest method for use with inductively coupled plasma spectrometry. Communications in Soil Science and Plant Analysis 11:969−980

    doi: 10.1080/00103628009367096

    CrossRef   Google Scholar

    [23]

    Isaac RA, Johnson WC. 1975. Collaborative study of wet and dry ashing techniques for the elemental analysis of plant tissue by atomic absorption spectrophotometry. Journal of Association of Official Analytical Chemists 58:436−40

    Google Scholar

    [24]

    Gavlak RG, Horneck DA, Miller RO. 1994. Plant, Soil, and Water Reference Methods for the Western Region. Corvallis, OR, USA: Western Regional Extension Publication.

    [25]

    Grismer ME. 1990. Leaching fraction, soil salinity, and drainage efficiency. California Agriculture 44:24−26

    doi: 10.3733/ca.v044n06p24

    CrossRef   Google Scholar

    [26]

    Franco JA, Martínez-Sánchez, JJ, Fernández JA, Bañón S. 2006. Selection and nursery production of ornamental plants for landscaping and xerogardening in semi-arid environments. The Journal of Horticultural Science and Biotechnology 81:3−17

    doi: 10.1080/14620316.2006.11512022

    CrossRef   Google Scholar

    [27]

    Sainju UM, Allen BL, Lenssen AW, Ghimire RP. 2017. Root biomass, root/shoot ratio, and soil water content under perennial grasses with different nitrogen rates. Field Crops Research 210:183−191

    doi: 10.1016/j.fcr.2017.05.029

    CrossRef   Google Scholar

    [28]

    Beard JB, Green RL. 1994. The role of turfgrasses in environmental protection and their benefits to humans. Journal of Environmental Quality 23:452−60

    doi: 10.2134/jeq1994.00472425002300030007x

    CrossRef   Google Scholar

    [29]

    Marcum KB, Pessarakli M. 2006. Salinity tolerance and salt gland excretion efficiency of bermudagrass turf cultivars. Crop Science 46:2571−74

    doi: 10.2135/cropsci2006.01.0027

    CrossRef   Google Scholar

    [30]

    Fu JM, Koski AJ, Qian YL. 2005. Responses of creeping Bentgrass to salinity and mowing management: growth and turf quality. HortScience 40:463−67

    doi: 10.21273/HORTSCI.40.2.463

    CrossRef   Google Scholar

  • Cite this article

    Hooks T, Masabni J, Ganjegunte G, Sun L, Chandra A, et al. 2022. Salt tolerance of seven genotypes of zoysiagrass (Zoysia spp.). Technology in Horticulture 2:8 doi: 10.48130/TIH-2022-0008
    Hooks T, Masabni J, Ganjegunte G, Sun L, Chandra A, et al. 2022. Salt tolerance of seven genotypes of zoysiagrass (Zoysia spp.). Technology in Horticulture 2:8 doi: 10.48130/TIH-2022-0008

Figures(4)  /  Tables(3)

Article Metrics

Article views(5025) PDF downloads(753)

ARTICLE   Open Access    

Salt tolerance of seven genotypes of zoysiagrass (Zoysia spp.)

Technology in Horticulture  2 Article number: 8  (2022)  |  Cite this article

Abstract: Seven zoysiagrass genotypes were evaluated for salt tolerance in a greenhouse study. The plant materials included Zoysia matrella 'Diamond', Z. japonica 'Palisades', three Z. matrella × Z. japonica hybrids DALZ 1701, DALZ 1713, and 'Innovation', and two Z. minima × Z. matrella hybrids (DALZ 1309 and 'Lazer'). Treatments included a control (nutrient solution) and two saline treatments representing moderate and high salt levels. The electrical conductivity (EC) was 1.3 dS m−1 for control and moderate (EC5) and high salinity (EC10) were 5.0 and 10.0 dS m−1, respectively. At the end of eight-weeks of treatment, the relative (percent control) shoot dry weight (DW) was greatest in 'Diamond' in EC10, and the relative root DW was greatest in DALZ 1309 in EC5. A cluster analysis based on the relative tissue dry weight identified 'Diamond', DALZ 1309, and DALZ 1713 as the most salt tolerant genotypes. Additionally, the green leaf area (GLA) index of 'Diamond' and DALZ 1713 were 98.8% and 100%, respectively, indicating excellent visual appearance under high salt levels. Bi-weekly clipping DW showed that 'Diamond' continued to produce biomass throughout the duration of the study under the EC10 treatment. Sodium (Na+) and chloride (Cl) content in the shoot tissue of the seven turfgrass genotypes indicated that lower concentrations corresponded to greater salt tolerance indicating exclusion of Na+ and Cl from the shoot tissue. Taken together, the genotypes 'Diamond' and DALZ 1713 were determined to be the most salt tolerant and recommended for use in areas with high soil or water salinity.

    • Turfgrass is an important landscape groundcover that is widely used, not only for its aesthetics, but also for its function, such as in lawns, parks, athletic fields, and golf courses. However, in arid and semi-arid regions of the United States where fresh water sources are limited, reclaimed water sources that typically have elevated salt levels are increasingly being used for landscape irrigation[1,2]. Reclaimed water, also known as recycled or reused water, is non-potable wastewater from a variety of sources, including residential, industrial, or stormwater runoff[2]. Reclaimed water can be treated to varying degrees to be suitable for a wide range of uses, including agricultural and landscape irrigation, that effectually relieves the demand from freshwater sources[3]. Because of these benefits, reclaimed water is increasingly being used for irrigation in arid and semi-arid regions, particularly for turfgrass areas[4]. However, reclaimed water sources are typically saline, with high concentrations of sodium and chloride, which can be detrimental to the growth and aesthetic quality of salt-sensitive plants[5]. Therefore, there is a continued need for salt tolerant turfgrasses for sustainable landscaping in arid or semi-arid regions in order to utilize saline, reclaimed water sources and conserve fresh water sources, especially in the Southwestern U.S.[6].

      Sodium and chloride are the two major soluble salts that can be detrimental to glycophytes at high concentrations[7]. When these salts accumulate in the rhizosphere, they can impose osmotic stress on the plant, which leads to the inhibition of water uptake and can rapidly reduce plant growth and even lead to mortality[8]. If salts are taken up by the roots and translocated to the shoots, then ionic stress can occur which can result in metabolic disruption in the cytosol of cells, as well as damage to chloroplasts by reactive oxygen species (ROS)[9]. Ionic stress can lead to leaf burn, which appears as brown and necrotic tissue[10] and can significantly degrade the quality of landscape and ornamental plants[11]. For turfgrasses, it is important that salt tolerant cultivars are not susceptible to ionic stress and leaf burn and can maintain an appearance that is aesthetically appealing under saline conditions[12].

      Traditionally, bermudagrasses (Cynodon spp.) have been used as warm season turfgrasses for landscaping in arid and semi-arid regions, although they are considered high water consumption plants[13]. Alternatively, zoysiagrasses (Zoysia spp.) are warm season turfgrasses that are moderately tolerant to salinity and have good potential for the selection and development of new salt tolerant cultivars that can be irrigated with saline, reclaimed water sources[14,15]. Both bermudagrasses and zoysiagrasses are being used as breeding material for new cultivars with the aim of having warm season, salt tolerant varieties that can grow well and maintain high visual quality, or greenness, when irrigated with saline water.

      Zoysiagrass was originally introduced to the U.S. in 1892, from East Asia, and have since been very influential in the turfgrass industry with more than 50 cultivars that have been developed, particularly for stress and pest tolerances. Two major species, Z. japonica and Z. matrella, readily hybridize with each other and are known for high quality turfgrasses primarily used in residential and commercial lawns, and golf courses[15]. Another species, Z. minima, is native to New Zealand and has a diminutive growth habit that has potential for use in golf course putting greens[16]. The recognized salt tolerance of zoysiagrasses is largely based on their ability to excrete salt out of their leaves through specialized salt glands[17,18]. However, differences in salt tolerance have been noted among species and cultivars. Marcuum & Murdoch[19] reported greater salt tolerance in Z. matrella compared to Z. japonica under solution culture up to 400 mM (approximately 30 dS m−1) NaCl. Qian et al.[20] reported differences in relative salt tolerance among 29 zoysiagrass experimental lines and cultivars under solution culture up to 42.5 dS m−1.

      In the present study, seven genotypes of zoysiagrass (Zoysia spp.), including several hybrids, were selected for evaluation for salt tolerance in a greenhouse study. The seven genotypes were developed by the turfgrass breeding program at Texas A&M AgriLife (Texas, USA). The parent species for these seven genotypes represent different leaf morphologies and traits: Z. japonica, wide leaf blade with drought and cold tolerance; Z. matrella, fine leaf blade with salt and shade tolerance; and Z. minima, very fine leaf blade with shade tolerance and good visual quality. The objectives of the study were to identify salt tolerant genotypes for the continued improvement of salt tolerance in turfgrass breeding programs, and for the potential use of these genotypes for landscaping under saline conditions in arid and semi-arid regions of the US.

    • Seven turfgrass genotypes were acquired from the Texas A&M AgriLife Turfgrass Breeding Program and used in this study, including a Z. matrella cultivar 'Diamond', a Z. japonica 'Palisades', three Z. matrella × Z. japonica hybrids DALZ 1701 ([(Z. matrella × Z. matrella) × Z. japonica] × Z. japonica), DALZ 1713 ([Z. japonica × (Z. matrella × Z. matrella)], and 'Innovation' (Z. matrella × Z. japonica), and two Z. minima × Z. matrella hybrids (DALZ 1309 and 'Lazer'). Approximately 10 rhizomes were transplanted into 10-cm (top diameter) round plastic pots (volume: 450 mL, height: 8.5 cm) filled with potting mix (Sun Gro, Agawam, MA, USA) and fertigated through the surface of the pot with 20-10-20 (N-P2O5-K2O) Peters Excel fertilizer (ICL, Sommerville, SC, USA) at a rate of 150 mg L−1 N, on an as-needed basis (when the substrate surface became dry). The nutrient solution was made by mixing 1.0 g of the above fertilizer to 1 L of tap water. The final electrical conductivity (EC) and pH was 1.3 dS m−1 and 6.5, respectively. Turfgrass cuttings were established and rooted for four weeks in a greenhouse at the Texas A&M AgriLife Research Center in Dallas, Texas (USA). A total of 24 pots with uniform growth of each genotype were selected. The genotypes were randomized on a greenhouse bench for the initiation of the saline treatments.

    • Two saline treatments were used, in addition to a non-saline control, in this experiment. In the control group, plants were irrigated with the nutrient solution as mentioned above. The two saline treatments were prepared by the addition of NaCl to the nutrient solution to achieve EC levels of 5.0 dS m−1 (EC5) and 10.0 dS m−1 (EC10). These two salinity levels were chosen based on available information in the literature on salt tolerance of other turfgrasses. The treatment solution for EC5 or EC10 was prepared by adding 230 g or 550 g of NaCl to 100 L nutrient solution. The actual EC and pH were recorded each time. Treatments were applied to the plants overhead on an as-needed basis. Approximately 200 mL was applied to each plant/pot per treatment application which provided a leaching fraction of approximately 35% to reduce the accumulation of salts in the substrate throughout the experiment. The experiment was arranged in a split-plot design with treatments randomized in greenhouse benches and genotypes randomized within treatments. There were eight pots (replicates) per treatment. Weekly measurements of the leachate EC and pH were recorded to track the salinity level in the substrate and rhizosphere. The leachate was collected via the 'Pourthru' method as described by Cavins et al.[21]. The treatments were initiated on 06 May 2020 and lasted eight weeks and was terminated on 02 July 2020.

    • The experiment was conducted in a greenhouse at the Texas A&M AgriLife Research Center in Dallas, TX, USA (32°59'13.2" N 96°45'59.8" W; elevation 131 m). The greenhouse air temperature was controlled by an evaporative cooling wall and two exhaust fans. A 50% shade fabric was used throughout the experiment to reduce sunlight and heat in the greenhouse. Throughout the experiment, greenhouse air temperature and photosynthetic active radiation (PAR) were recorded by a datalogger (Campbell Scientific, Logan, UT, USA). The air temperature and quantum sensor (for PAR measurement) were installed right above the bench to capture the actual air temperature and light intensity near the plant canopy. The daily average air temperature during the experiment was 26.0 ± 3.56 °C (mean ± standard deviation) and the average daily light integral (DLI) was 12.0 ± 2.98 mol m−2 d−1.

    • Throughout the experiment, the plants were clipped on a biweekly schedule and the clippings were collected and dried in a drying oven at 70 °C for dry weight determination. Clipping was accomplished by hand with scissors and a ruler, following a treatment application. The turfgrass was clipped to a 2-cm height and the perimeter of the pots were also trimmed. Additionally, after each biweekly trimming, the percent canopy green leaf area (GLA) was determined visually by two persons to assess the quality of the plants under the saline treatments. At harvest, shoot and root tissue were separated and dried in a drying oven at 70 °C for biomass determination. Roots were washed of substrate and rinsed briefly in reverse osmosis water before being placed in paper bags and dried in the drying oven. Following dry weight determination, three shoot samples from each treatment were ground in a Wiley mill (Thomas Scientific, Swedesboro, NH, USA) to pass a 40-mesh screen. Shoot tissue mineral contents were analyzed using inductively coupled plasma mass spectrometry (ICP-MS) using the methods described by Havlin & Soltanpour[22] and Isaac & Johnson[23]. Shoot tissue chloride content was determined by extraction with 2% acetic acid and analyzed using an M926 Chloride Analyzer (Cole Parmer Instrument Company, Vernon Hills, IL, USA) according to the methods described by Gavlak et al.[24]

    • There was a total of three treatments and seven genotypes with eight replications each (N = 168). Data were analyzed as a two-way analysis of variance (ANOVA) with an alpha of 0.05 using JMP 15 (SAS, Cary, NC, USA). Means were separated using Tukey's Honest Significant Difference (HSD) test with an alpha of 0.05. Relative shoot dry weight (DW) was calculated as the shoot DW in the saline treatment/average shoot DW in control × 100%. Relative root DW and relative total DW in percentage were calculated in a similar fashion compared to control. Student's t-test was used for comparing the relative growth parameters between the two saline treatments.

    • Throughout the duration of the study, leachate EC of the control, moderate, and high salt treatments averaged 1.7, 7.4, and 14.2 mS cm−1, respectively (Fig. 1). Although leachate EC increased steadily in the moderate and high salt treatments throughout most of the study due to a buildup of salts in the substrate, the averages of the treatments were significantly different, as expected. The maximum EC of the salt treatments peaked during week six, at 11.2 and 18.8 dS m−1 for the moderate and high salt treatments, respectively. During weeks seven and eight, the EC of the salt treatments started to decline which was attributed to the retention of moisture in the substrate due to reduced water uptake by the osmotically stressed grasses, which ultimately lead to increased leaching fractions during irrigation.

      Figure 1. 

      Electrical conductivity (EC) of leachate collected from seven turfgrass genotypes treated with control or saline solutions (EC5 or EC10) for a total of eight weeks. Vertical bars indicate standard error (n = 5).

    • For relative (percent control) shoot DW, there were no treatment differences but there were significant genotype differences in the EC10 treatment (Table 1), as expected, with Diamond showing the greatest increase of 130% compared to the control (Fig. 2). For relative root DW, there were significant treatment differences in Lazer and DALZ 1713, with reductions of 20% and 40%, respectively, in the EC10 treatment. There were genotypic differences in the EC5 treatment only, with DALZ 1309 showing the greatest increase of 130% compared to the control, while DALZ 1701 and Palisades decreased by 15% and 17%, respectively, compared to the control. Overall, for total DW, there were significant treatment differences in Lazer and DALZ 1713, with reductions of 20% and 22%, respectively, in the EC10 treatment. Both salt treatments had significant genotype differences, with DALZ 1309, DALZ 1713, and Diamond showing the greatest increases of 120%, 122%, 118%, respectively, in the EC5 treatment compared to the control, while Diamond showed the greatest increase of 121% in the EC10 treatment.

      Table 1.  ANOVA summary of the response variables of the seven zoysiagrass genotypes irrigated with a nutrient solution (control) or saline solution at electrical conductivity (EC) of 5 dS m−1 or 10 dS m−1 for eight weeks. The response variables are shoot DW (dry weight), root DW, total DW, relative shoot DW (R. shoot DW), relative root DW (R. root DW), relative total DW (R. total DW), green leaf area (GLA), cumulative clipping DW, shoot sodium (Na) and chloride (Cl) concentration.

      SourceShoot DWRoot DWTotal DWR. Shoot DWR. Root DWR. Total DWGLAClipping DWShoot Na+Shoot Cl
      Model0.0004< 0.0001< 0.00010.00030.0001<.0001< 0.0001< 0.0001< 0.0001< 0.0001
      Treatment (T)0.0080.00010.001NS0.00020.0013< 0.0001< 0.0001< 0.0001< 0.0001
      Genotype (G)0.0006< 0.0001< 0.0001< 0.00010.0014< 0.0001< 0.0001< 0.0001< 0.0001< 0.0001
      T × GNS0.04240.0109NSNSNS< 0.0001< 0.0001< 0.0001< 0.0001

      Figure 2. 

      Relative (percent of control) dry weight (DW) of shoot and root tissue, and the total (shoot + root) of the seven turfgrass genotypes treated with control or saline solutions (EC5 or EC10) for a total of eight weeks. Bars represent standard error (n = 8). Different letters indicate significant differences among genotypes for the same treatment according to Tukey's HSD test (P < 0.05). That is, the comparison was made for EC5 (red bars) or EC10 (green bars) separately. For those without any letters such as EC5 for shoot DW, no difference was observed. Asterisks indicate significant differences between treatments (EC5 and EC10) according to Student's t-test (P < 0.05). No asterisks mean no differences.

    • A cluster analysis was performed on the relative shoot and root DW of the seven turfgrass genotypes treated with moderate and high salinity and found two distinct clusters as indicated by the distance graph (Fig. 3). Cluster 1 (red) indicates the least salt tolerant genotypes (based on lowest relative tissue DW) and included Lazer, DALZ 1701, Innovation, and Palisades. Cluster 2 (green) indicates the most salt tolerant genotypes (based on greatest relative tissue DW) and included DALZ 1309, DALZ 1713, and Diamond.

      Figure 3. 

      Hierarchal cluster analysis based on relative (percent of control) tissue dry weight (DW) of the seven turfgrass genotypes treated with control or saline solutions (EC5 or EC10) for a total of eight weeks. Cluster 1 (red) indicates the least salt tolerant genotypes and Cluster 2 (green) indicates the most salt tolerant genotypes. The two clusters were determined by the distance graph at the bottom of the figure that shows the best separation between clusters.

    • The GLA Index averaged 98.7 in the control, 96.7 in the EC5 treatment, and 92.4 in the EC10 treatment (Table 2, Supplemental Fig. S1). There were significant treatment, genotype, and treatment × genotype interactions in GLA (Table 1). The interactions were attributed to DALZ 1309 showing substantial reductions in GLA in the EC10 treatment, while other genotypes such as DALZ 1713, showing no reductions. In fact, DALZ 1309 showed the greatest reductions in GLA in all treatments, specifically 95.6, 88.8, and 69.4 in the control, EC5, and EC10 treatments, respectively. In contrast, DALZ 1701, DALZ 1713, Diamond, and Palisades showed no significant differences in GLA among the treatments and maintained excellent scores under the saline irrigation treatments.

      Table 2.  Green Leaf Area (GLA) index of the seven turfgrass genotypes treated with control or saline solutions (EC5 or EC10: electrical conductivity at 5 or 10 dS m−1) for a total of eight weeks. Means and standard errors are presented (n = 8). The GLA was assessed visually following a clipping.

      GenotypeControlEC5EC10
      Lazer100.0 ± 0.0Aa99.4 ± 0.6Aab95 ± 2.1Ab
      DALZ 130995.6 ± 2.0Ba88.8 ± 3.0Ca69.4 ± 6.4Bb
      DALZ 170198.1 ± 0.9ABa99.4 ± 0.6Aa97.5 ± 0.9Aa
      DALZ 1713100.0 ± 0.0Aa99.4 ± 0.6Aa100.0 ± 0.0Aa
      Diamond99.4 ± 0.6ABa100.0 ± 0.0Aa98.8 ± 1.3Aa
      Innovation98.1 ± 1.3ABa91.9 ± 2.5BCab88.1 ± 3.3Ab
      Palisades100.0 ± 0.0Aa98.1 ± 0.9ABa98.1 ± 0.9Aa
      Different letters indicate significant differences Tukey's HSD test; uppercase among genotypes and lowercase among treatments.
    • For Clipping DW, there were significant treatment, genotype, and treatment × genotype interactions (Table 1). The interactions can be explained by some genotypes showing an increase in clipping DW in all treatments throughout the study, while other genotypes showed a decrease, particularly in the EC10 treatment during the final weeks of the study (Fig. 4). There were significant treatment differences as early as week 2 in DALZ 1701 and Palisades, and in all genotypes for the remaining weeks of the study. Overall, clipping DW was greatest in the control, followed by the EC5 and then EC10 treatment. By the end of the study, clipping DW in EC5 and EC10 plateaued or declined in all genotypes except for Lazer, DALZ 1309, DALZ 1701 (EC5), and Diamond (EC10), which still showed increases despite the high saline conditions as indicated by the leachate EC. Declines in the control treatment towards the end of the study in DALZ 1713, Innovation, and Palisades can be attributed to the plants exceeding the growth capacity of the containers.

      Figure 4. 

      Bi-weekly clipping dry weight (DW) of the seven turfgrass genotypes treated with control or saline solutions (EC5 or EC10: electrical conductivity at 5 or 10 dS m−1) for a total of eight weeks. The plants were clipped to a height of 2-cm. Bars represent standard error (n = 8). Significant differences among treatments per week are indicated by asterisks (*, P < 0.05; **, P < 0.01; and ***, P < 0.001).

    • There were significant treatment, genotype, and treatment x genotype interactions for both sodium (Na+) and chloride (Cl) concentrations in the shoot tissue (Table 1). Overall, the Na+ and Cl concentrations in the shoot tissue increased in the salt treatments compared to the control, as expected due to the higher amount of Na+ and Cl ions in the salt treatments (Table 3). The average Na+ concentration in the tissue of plants treated with control, EC5, and EC10 was 2.00, 8.49, and 12.04 mg g−1, respectively. For Cl, the average amount was 6.41, 12.53, and 18.54 mg g−1 in the tissue of plants treated with control, EC5, and EC10, respectively. For Na+ there were no significant differences among genotypes in the EC5 treatment, although in the EC10 treatment genotype DALZ 1309 had the greatest concentration (15.68 mg g−1) while DALZ 1701 had the least (7.92 mg g−1). For Cl in the EC5 treatment, the genotype Palisades had the greatest concentration (15.97 mg g−1) while Lazer and DALZ 1701 had the least (9.72 and 9.68 mg g−1, respectively). In the EC10 treatment, the genotypes DALZ 1309 and Innovation had the greatest concentrations (23.97 and 25.18 mg g−1), while Lazer, DALZ 1701, and Diamond, had the least (14.00, 13.25, 15.38 mg g−1, respectively). Regarding the significant interaction, this can be explained by most genotypes showing substantial increases of Na+ and Cl between the EC5 and EC10 treatments, while certain genotypes showed no differences, such as Lazer in the EC5 treatment and Palisades in the EC10 treatment.

      Table 3.  Sodium (Na+) and chloride (Cl) content in the tissue of the seven turfgrass genotypes that were treated with control or saline solutions (EC5 or EC10: electrical conductivity at 5 or 10 dS m−1) for a total of eight weeks.

      GenotypeControlEC5EC10
      Na+
      Lazer1.70 ± 0.04BCb8.76 ± 0.90Aa9.90 ± 0.42CDa
      DALZ 13092.07 ± 0.24ABCc8.89 ± 0.13Ab15.68 ± 1.43Aa
      DALZ 17011.39 ± 0.06Cc6.58 ± 0.41Ab7.92 ± 0.07Da
      DALZ 17132.79 ± 0.19Ac9.33 ± 0.28Ab12.86 ± 0.72ABCa
      Diamond1.81 ± 0.02BCc7.58 ± 0.33Ab10.78 ± 0.36CDa
      Innovation1.77 ± 0.15BCc9.46 ± 1.02Ab14.92 ± 0.56ABa
      Palisades2.45 ± 0.37ABc8.84 ± 1.12Ab12.19 ± 0.30BCa
      Cl
      Lazer7.55 ± 0.30Ac9.72 ± 0.56Db14.00 ± 0.28Ba
      DALZ 13095.65 ± 0.23Ab12.20 ± 0.02BCDb23.97 ± 2.84Aa
      DALZ 17015.67 ± 0.22Ac9.68 ± 0.68Db13.25 ± 0.77Ba
      DALZ 17136.55 ± 0.10Ac13.75 ± 0.31ABCb19.15 ± 0.40ABa
      Diamond6.43 ± 0.10Ac11.52 ± 0.10CDb15.38 ± 0.15Ba
      Innovation5.72 ± 0.41Ac14.88 ± 0.75ABb25.18 ± 1.45Aa
      Palisades7.27 ± 0.99Ab15.97 ± 1.36Aa18.87 ± 0.96ABa
      Means and standard errors are presented (n = 8). Different letters indicate significant differences Tukey's HSD test; uppercase among genotypes and lowercase among treatments.
    • Throughout the study, although the salt treatments remained fixed, salt accumulation occurred in the substrate as indicated by the increase in leachate EC of the salt treatments. Salt accumulation in the substrate depends on many factors, including salinity of the irrigation water, irrigation frequency, leaching fraction, and substrate type. The leaching fraction, simplified, is the percent of irrigation that drains out of the substrate[25]. Higher leaching fractions can flush ions, including Na+ and Cl, away from the root zone and out of the substrate. In this study, a leaching fraction of approximately 35% was applied to slow down the accumulation of salts in the substrate without wasting too much irrigation. Despite this, leachate EC increased throughout the study, most notably in the EC10 treatment. This imposed additional osmotic and/or ionic stress on the grasses beyond the fixed treatment salinities of 5.0 and 10 dS m−1. Nevertheless, this is representative of irrigation regiments in arid landscaping, where low irrigation volumes and leaching fractions are commonly applied[26].

      Biomass was reduced by the high salt treatment more notably in the root tissue compared to the shoot tissue. In fact, shoot tissue increased marginally relative to the control in most genotypes when treated with salt, which is indicative of salt tolerance and halophytes[24]. However, root tissue sensitivity to salt stress is rather unique, since typically shoot tissue is more sensitive[10]. Chavarria et al.[1] observed both increases and reductions in root mass among eight turfgrass genotypes when treated with 15 and 30 dS m−1 salinity, when compared to the control. Additionally, in the present study the cluster analysis based on shoot and root tissue biomass identified the genotypes DALZ 1309, DALZ 1713, and Diamond as the most salt tolerant, which corresponds with the relative root tissue DW in the EC5 treatment. Therefore, our results indicate that root tissue biomass in turfgrasses might be a greater indication of salinity tolerance than shoot tissue biomass. This could be because grasses have relatively small leaf surface area and large root/shoot ratios compared to other plants[27].

      Salt tolerance for ornamental crops not only depends on growth under saline conditions, but also visual appearance, as salinity can impose ionic stress to plants which can lead to leaf burn[12]. For turfgrasses, this is especially true due to its primary use for aesthetics and environmental benefits in residential, recreational, or commercial landscaping[28]. For turfgrass managers, greenness can be a more important trait than shoot yield. Marcumm & Pessarakli[29] reported GLA ranges of 7% to 84% in eight Distichlis spicata turfgrass genotypes treated with up to 1.0 mol L−1 (58.5 g L−1) NaCl for one week. In the present study, our results indicate relatively high GLA and hence, excellent visual quality in most genotypes even when treated with high (EC10) salinity for eight weeks. In contrast, the significant reductions in GLA in the genotype DALZ 1309 in the EC10 treatment indicates less salt tolerance and susceptibility to ionic stress.

      Continued growth under saline conditions is another desirable trait in turfgrasses, indicating long-term establishment. However, mowing is a necessary management practice for turfgrass and has been shown to affect the salinity tolerance of turfgrass varieties. For example, clipping yield of creeping bentgrass (Agrostis palustris) was reduced the most under low mowing height (6.4 mm) compared to high mowing height (25.4 mm), when treated with salinity ranging from 5 to 15 dS m−1[30]. Similarly, our results showed that high salinity (EC10) reduced clipping DW compared to the control at a mowing height of 2.0 cm. However, at the end of the study, clipping DW tended to decline in the EC10 treatment in the genotypes DALZ 1701, Innovation, and Palisades, indicating less tolerance to salinity at this mowing height, whereas the remaining genotypes showed marginal gains in clipping DW, most notably Diamond, indicating greater tolerance to salinity at the respective mowing height.

      A key mechanism of salt tolerance in plants is the ability to exclude Na+ and Cl from the leaf tissue by various means, such as sequestration in the cell vacuole or excretion through specific glands in the leaf[6]. Therefore, concentration of salts in the shoot tissue can be an indication of salt tolerance and/or mechanisms of salinity tolerance in a specific plant. In the present study, Na+ and Cl concentrations in the shoot tissue increased in all genotypes when treated with high salt (EC10), indicating salt accumulation in the shoot tissue. However, our results also showed genotypic variation in salt accumulation, indicating different mechanisms for dealing with the salts. For example, genotype DALZ 1701 had the lowest concentration of Na+ which correlated with its excellent visual and growth parameter, and a potential explanation for this is it could more affectively exclude Na+ or excrete it from the leaves, which could contribute to its excellent visual quality (high GLA) and overall good salt tolerance in the present study. Chavarria et al.[1] reported Na+ concentrations in the shoot tissue ranging from 7.2 to 22.4 mg g−1 of eight warm-season turfgrasses when treated with 15 dS m−1, which were comparable values to those reported here considering the higher salt treatment. However, they also reported that salt excretion correlated with salt tolerant genotypes. Nevertheless, the ability to maintain high Na+ and Cl concentrations in the shoot tissue while maintaining good visual quality and growth, indicates tolerance to osmotic and ionic salinity stress, which our results demonstrated.

    • Our results primarily indicated genotypic variation present within zoysiagrasses for the improvement of salt tolerance. The genotypes Zoysia matrella 'Diamond', Z. japonica 'Palisades', three Z. matrella x Z. japonica hybrids (DALZ 1701, DALZ 1713, and 'Innovation'), and two Z. minima × Z. matrella hybrids (DALZ 1309 and 'Lazer') showed variation in potential for use in landscaping with saline irrigation in arid regions for the purpose of conserving freshwater resources and maintaining aesthetic and environmental benefits of green groundcover. Based on the growth, visual quality (GLA), and physiological results of this study, the genotypes Diamond and DALZ 1713 exhibited superior salt tolerance across multiple growth and physiological traits evaluated while DALZ 1701 expressed potential for improved salinity tolerance from its ability to exclude salt and maintain high visual quality.

      • Funding for this project is provided by USDA NIFA to Project No. 2017-68007-26318, through the Agriculture and Food Research Initiative, Water for Agricultural Challenge Area, and hatch project TEX07726.

      • The authors declare that they have no conflict of interest.

      • Supplemental Fig. S1 Photos of the seven genotypes that were treated for 8 weeks: control (left), EC5 (middle), and EC10 (right). The genotype arrangement in each treatment was randomized. Photos were taken right after trimming.
      • Copyright: © 2022 by the author(s). Published by Maximum Academic Press, Fayetteville, GA. This article is an open access article distributed under Creative Commons Attribution License (CC BY 4.0), visit https://creativecommons.org/licenses/by/4.0/.
    Figure (4)  Table (3) References (30)
  • About this article
    Cite this article
    Hooks T, Masabni J, Ganjegunte G, Sun L, Chandra A, et al. 2022. Salt tolerance of seven genotypes of zoysiagrass (Zoysia spp.). Technology in Horticulture 2:8 doi: 10.48130/TIH-2022-0008
    Hooks T, Masabni J, Ganjegunte G, Sun L, Chandra A, et al. 2022. Salt tolerance of seven genotypes of zoysiagrass (Zoysia spp.). Technology in Horticulture 2:8 doi: 10.48130/TIH-2022-0008

Catalog

  • About this article

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return