#### Abstract

The current liquefaction evaluation methods mainly focus on the success rate for liquefied sites so that the evaluation result tends to be conservative at different seismic intensities. Therefore, a new formula about liquefaction evaluation by introducing logistic regression theory is proposed to solve the deficiencies of the current evaluation method, which is based on 225 sets of shear wave velocity data reported by Andrus. The reliability of the new formula is verified based on 336 sets of *V*_{s} data collected from the Kayen database. The performance of the new formula on liquefaction evaluation is compared with existing liquefaction evaluation methods including the Andrus method and the Chinese code method. Compared with the Andrus method and Chinese code method, the success rates of liquefaction evaluation given by the new formula under different seismic intensities are more balanced between liquefied site and nonliquefied site. The new formula at 50% probability of liquefaction is more adaptable for a wide range of seismic intensities, ground water table, and sand buried depth. In addition, the new formula at different probabilistic levels of liquefaction can be adopted based on the importance of the engineering site in risk analysis.

#### 1. Introduction

Soil liquefaction induced by the earthquake may often occur on some sites. The damage induced by soil liquefaction was not fully recognized until the occurrence of the 1964 Niigata and the 1964 Alaska earthquake [1], which destroyed structures and infrastructures and caused loss of life. It is necessary to develop a proper liquefaction evaluation method in order to prevent or minimize such damage.

At present, numerous researchers have tried to evaluate the soil liquefaction by standard penetration test (SPT), cone penetration test (CPT), Becker hammer test (BHT), and shear wave velocity (*V*_{s}) test. Among these methods, the shear wave velocity (*V*_{s}) test has wider applicability for liquefaction evaluation than other methods [2]. The shear wave velocity test can also be conducted in gravelly soils where SPT and CPT are unreliable [3–6] since they cannot be conducted in soils containing gravels, cobbles, or boulders. In contrast to SPT and CPT, *V*_{s}-based method has the advantage that no corrections are required based on soil characteristics [7–11]. In addition, shear wave velocity tests measure the shear modulus (stiffness) of the soil at low strains, which can also be obtained in the laboratory, so that comparison can be made between laboratory test and field test [5, 6]. Thanks to these advantages, *V*_{s}-based methods have great prospects for liquefaction evaluation [4, 12].

Dobry et al. [13] confirmed that the shear wave velocity is relevant to liquefaction resistance. In addition, both shear wave velocity values and blow counts are similarly influenced by many factors such as void ratio, state of stress, and stress history. Accordingly, a simple procedure incorporating a cyclic stress ratio (CSR)-*V*_{s} triggering resistance curve was proposed by Seed et al. [14]. Andrus and Stokoe [15] proposed guidelines for using the *V*_{s}-based liquefaction evaluation procedure, which is subsequently updated by Andrus et al. [16]. Yamazaki et al. [17] presented a new simplified liquefaction prediction and assessment method considering the influence of the waveforms and durations of earthquakes. EI-Sekelly et al. pointed out that the existing *V*_{s}-based liquefaction charts could favorably evaluate the liquefaction of uncompacted clean and silty sand recent fills but might be too conservative for heavily preshaken deposits in highly seismic areas. Chen et al. [18] confirmed that the cyclic resistance and the corrected shear wave velocity are uniquely related for a range of sandy soils and presented a new *V*_{s}-based approach to assess the liquefaction triggering of saturated sandy soils. Akbari-Paydar and Ahmadi [19] indicated that the commonly used CRR (cyclic resistance ratio)-*V*_{s} correlation curves might be conservative for silty and clayey sands; thus, it can only be used as an initial estimation of liquefaction resistance by performing a series of cyclic triaxial and bender element tests on reconstituted specimens of clean sand with different contents.

The deterministic methods discussed above only give the result of liquefaction or nonliquefaction of the site. However, it is unreasonable to give a deterministic evaluation result due to the uncertainty of seismic characteristics and the nonlinearity of soil. In addition, these methods commonly pursue the success rate of liquefaction, leading to a low success rate of nonliquefaction. Using these deterministic methods to treat the foundation or design buildings is likely to increase the construction cost. A probabilistic liquefaction evaluation method can potentially account for the uncertainties of soil properties, model parameters, and seismic input characteristics, which is also in line with the current development on earthquake engineering [1, 2, 5]. Thus, numerous scholars pay more attention to the research on probabilistic approaches to liquefaction analysis. Liao et al. [20] recommended regression models for evaluating the probability of liquefaction through corrected/normalized SPT (*N*_{1})_{60} values (the corrected SPT resistance normalized to a hammer energy ratio of 60%, as per Seed et al. [21]) and comparisons are made with other methods of liquefaction analysis. Similar studies have been conducted by logistic regression analysis based on field performance data [22–24]. Shen et al. [5] developed an updated log-log model for liquefaction potential evaluation based on shear wave velocity and concluded that the probabilistic models are quite satisfactory when they are checked against the case histories. Rahmanian and Rezaie [25] presented that the severity of liquefaction occurrence in the studied area by the *V*_{s}-based method was lower than that by SPT-based method through the comparison of the safety factor and liquefaction potential index. Johari et al. [2] developed a soil liquefaction model based on shear wave velocity using the jointly distributed random variables (JDRV) method. They pointed out that the probability distribution of the liquefaction safety factor obtained by the JDRV method is very close to that predicted by the Monte Carlo simulation. Andrus et al. [26, 27] concluded that the wave velocity-based method is more conservative than the penetration-based method for liquefaction evaluation based on complied Holocene data including 43 Holocene-age sand layers in California, South Carolina, Canada, and Japan.

The existing probability methods usually use CSR and CRR [7, 27] to establish the liquefaction evaluation formula. However, this is not in line with the Code for Investigation of Geotechnical Engineering in China (GB 50021-2009) for liquefaction evaluation [28], which is referred to as Chinese code hereafter. Therefore, a new liquefaction evaluation model combined with logistic regression theory, which is based on 225 sets of *V*_{s} data from the Andrus database [16], is proposed to solve the deficiencies of the current liquefaction evaluation method. The reliability of the new model is verified by 329 sets of data from the Kayen database and the performance of the new model on liquefaction evaluation is compared with the existing liquefaction evaluation methods.

#### 2. Chinese Code Method and Andrus Method on Liquefaction Evaluation

##### 2.1. Data Collection for Liquefaction Evaluation

The liquefaction evaluation method established in this study is based on the authoritative database published by Andrus et al. [16], which is summarized in Table 1. Many scholars use the Andrus database to develop a new model for liquefaction evaluation and the authority of the database is widely recognized [2, 29]. The database was involved in 225 groups of shear wave velocity data, which were collected from 26 earthquakes and more than 70 sites. The fundamental parameters such as (moment magnitude), CSR (cyclic stress ratio), (overburden pressure-corrected shear wave velocity), (sandy-layer buried depth), (depth of ground water table), and (peak horizontal ground surface acceleration) are provided in the database. According to the Chinese seismic intensity table [30], these data are divided into seismic intensities VI, VII, VIII, and IX, as shown in Table 2. It should be noted that all sites that are distributed in seismic intensity VI are nonliquefied. At present, the existing code for liquefaction evaluation in China does not consider the liquefaction in seismic intensity VI sites. On the one hand, few data about shear wave velocity are available in seismic intensity VI sites located in China. On the other hand, the site distributed in intensity VI is basically not liquefied. However, there are 34 sites provided in the Andrus database in seismic intensity VI and they will be considered in establishing a new model for liquefaction evaluation in this study in order to determine the lower boundary of the liquefaction critical curve more reasonably.

Figure 1 presents the data frequency distributions with the buried depth of sandy layer for case histories based on 225 groups of sites [16]. It is shown that there exists a slight difference in the distribution of sandy-layer buried depth between liquefied site and nonliquefied site. At liquefied sites and nonliquefied sites, most of the sandy-layer buried depth ranges from 2 m to 10 m. Figure 2 depicts the data frequency distributions with the depth of ground water table. The depth of ground water table ranges from 0 m to 8 m. At liquefied sites, the depth of water table mainly distributes from 1 m to 4 m, while at nonliquefied sites mainly from 0 to 3 m.

**(a)**

**(b)**

**(a)**

**(b)**

##### 2.2. Chinese Code Method for Liquefaction Evaluation

In Chinese code [28], the critical shear wave velocity value is used as the index for the soil liquefaction evaluation. The liquefaction evaluation formula based on shear wave velocity has been widely used in China, which is written as follows:where is the critical shear wave velocity and is a reference value for shear wave velocity, which are 65 m/s, 95 m/s, and 130 m/s corresponding to seismic intensity VII, VIII, and IX, respectively.

When the measured shear wave velocity of the soil layer is less than the corresponding , the site is judged as a liquefied site. Otherwise, the site is judged as a nonliquefied site.

##### 2.3. Andrus Method for Liquefaction Evaluation

Based on 225 sets of data, Andrus et al. [16] proposed a liquefaction evaluation method which is widely used all over the world and referred to as the Andrus method hereafter. The Andrus method is expressed as follows:where CRR is soil liquefaction resistance; is in situ measured shear wave velocity in m/s, which can be obtained by several seismic tests [6]; is reference stress of 100 kPa; is overburden effective stress in kPa; MSF is a magnitude scaling factor related to ; and is limiting upper value of in m/s, which is dependent of fine content, *F*_{c}, in percent by mass.

can be estimated as follows:

Liquefaction is predicted to occur when CSR > CRR, and liquefaction is predicted not to occur when CSR < CRR. CSR is cyclic shear stress ratio which can be calculated by the simplified formula proposed by Seed and Idriss [31] and written as follows:where is the average equivalent uniform cyclic shear stress suffered by an earthquake in kPa; is overburden total stress in kPa; is the acceleration of gravity in ; and is a shear stress reduction coefficient to adjust for the flexibility of the soil profile, which can be calculated as follows [32]:

##### 2.4. Liquefaction Evaluation Comparison between Chinese Code and Andrus Method

Chinese code method evaluates liquefaction sites which locate in seismic intensity VII, VIII, and IX. Thus, the 34 sets of data in seismic intensity VI are removed when compared with the Andrus method. It should be noted that clay content is not provided in the Andrus database. Therefore, the value of clay content is taken as 3% when using the Chinese code method [29]. The fine content is assumed as FC ≤ 5% when it is not provided in the database [29].

Table 3 shows the success rates of liquefaction evaluation based on the Andrus database by Chinese code method and Andrus method. It should be noted that the total success rate is defined as the ratio of site number evaluated successfully to total site number. Taking the total success rate in seismic intensity VII as an example, as displayed in Table 2, there are 23 liquefied sites and 46 nonliquefied sites in seismic intensity VII. Of these, 21 liquefied sites and 19 nonliquefied sites are successfully evaluated by the Andrus method. Thus, the total success rate, , is obtained as follows:

In seismic intensity VII, the success rates for liquefied sites given by the Chinese code method and Andrus method are 56.5% and 91.3%, respectively. The success rates at nonliquefied sites reach 39.1% and 41.3%, respectively. That is to say, the Chinese code method tends to be dangerous in seismic intensity VII. In seismic intensity VIII and IX, for liquefied sites, the success rates of the two methods are both 100%; for nonliquefied sites, the success rates of both methods are lower than 35%. The success rate for liquefied sites is satisfactory, but these two methods stay conservative in seismic intensity VIII and IX. For liquefied sites, the total success rates using the Chinese code method and the Andrus method are 89.6% and 97.9%, respectively. However, the total success rates for nonliquefied sites are only 29.5% and 34.7%, respectively. These two methods with total success rates smaller than 70% tend to be conservative and reach an unacceptable level. Therefore, it is essential to improve the existing liquefaction evaluation methods.

#### 3. A New Liquefaction Evaluation Formula

##### 3.1. The Theoretical Basis

The existing liquefaction evaluation methods based on shear wave velocity have made great progress; however, they are not widely used in China. This mainly results from that CSR and CRR indexes are not used in existing Chinese codes. This study intends to introduce the parameters with higher liquefaction discrimination to establish a new liquefaction evaluation formula. The parameters affecting the liquefaction probability are analyzed herein.

According to Seed and Idriss [31], the minimum cyclic shear stress of soil to resist liquefaction, , can be expressed as follows:where is the dynamic stress amplitude in kPa, is the effective confining pressure in cyclic triaxial test in kPa, and is a stress correction factor related to the relative density of soil. The values of in equation (14) are summarized in Table 4 which are recommended by Seed and Idriss [31].

Furthermore, the cyclic stress ratio, , is expressed as follows:where is the relative density of sand and is cyclic stress ratio at a relative density of 50%.

Liquefaction potential, , can be expressed as follows [33]:

Umar et al. [34] pointed out that earthquake magnitude, , plays an important role in the assessment of liquefaction. To consider the influence of on and facilitate comparison of sites from different earthquakes, is obtained as follows:

Combining equations (9)–(12) yields the following:

For the sake of brevity, a new parameter, *k*, was defined as follows:

Furthermore, the liquefaction potential is obtained as follows:where is a coefficient related to ; MSF is relevant to ; and and are dependent on and . The sandy-layer buried depth and depth of ground water table can reflect the effective overburden stress. It is difficult to directly obtain the value of . The shear wave velocity or SPT counts can be alternative to reflect the relative density of the soil.

can be calculated fromwhere *e*, , and are the natural, maximum, and minimum void ratios, respectively.

Hardin and Drnevich [35] have shown that for sands, the shear modulus at low strain, , can be calculated fromwhere *n* is a parameter related to soil type and is usually taken as 0.5 and is the principal effective stress in kPa. The fitting parameters, *A* and , are optimized in order to obtain the best fit of the experimental points. *A* is dependent on the particle shape and gradation of soil; is dependent on void ratio. The optimization technique is a classical minimization approach [36].

Combining equations (16)–(18) yields

It should be noted that is equal to for the isotropic consolidation of the cyclic triaxial test. By combining equation (3) with equation (19), can be written as follows:

It is worth noting that can be obtained by substituting *A*, , , , and into equation (20). The value of can be obtained from Table 4. It can be seen from equation (20) that largely depends on the relative density of the soil so that overburden stress-corrected shear wave velocity () is used to represent since is easy to get. Therefore, liquefaction potential can be expressed by the five descriptive parameters including , *k*, , , and . These parameters can be regarded as independent variables influencing liquefaction occurrence. The result of liquefaction evaluation can be classified into liquefaction and nonliquefaction so that liquefaction evaluation is a two-category variable event. The logistic regression method can deal with complex and multivariable event, which is suitable for liquefaction evaluation by giving the results of probability.

##### 3.2. Establishment of New Liquefaction Evaluation Formula

Due to the uncertainty of earthquake occurrence and the fact that soil is a highly nonlinear material, there is no clear boundary between liquefaction and nonliquefaction, so the expression of the probability of liquefaction evaluation is more suitable for reflecting the liquefaction issue. Based on the derivation of the liquefaction potential, the five parameters including , *k*, , , and , are selected to establish the logistic regression-based probability model.

The logistic regression model proposed by Liao et al. [20] can be expressed as follows:where *X* is a vector of explanatory variable, which includes , *k*, , , and ; is the probability of liquefaction, which is a function of *X*; and , , , , , and are regression coefficients of the logistic model.

The vector of coefficients is estimated through maximization of the likelihood of function *L*(*X*; *θ*), which is given bywhere *m* is the number of total sites and *Y* is equal to 1 or 0 corresponding to liquefied sites and nonliquefied sites, respectively.

According to the principle of maximization of the likelihood function, the obtained is the best estimate of when *L*(*X*; *θ*) takes the maximum value. By solving the partial derivative of *L*(*X*; *θ*) to , a series of likelihood function equations can be established as follows:

Substituting 225 sets of *V*_{s} data from the Andrus database into equation (23) yields the vector of coefficient through the Newton–Raphson algorithm [20]. The probability of liquefaction based on *V*_{s} can be written as follows:

The critical shear wave velocity, , at a given probabilistic level can be obtained by equation (24) as follows:where is the given probabilistic level.

When the overburden pressure-corrected shear wave velocity, *V*_{s1}, is less than *V*_{scr} obtained by equation (25), the site is judged as being liquefied; otherwise, the site is judged as not liquefied.

#### 4. Validation of the New Formula

##### 4.1. Liquefaction Evaluation at Different Probabilistic Levels

In comparison with the Andrus method and the Chinese code method, the reliability of the new formula is verified based on the Andrus database. In order to evaluate the balance of the liquefaction evaluation between liquefied site and nonliquefied site, the index of DI is defined herein as the difference in success rate between liquefied site and nonliquefied site. Figure 3 lists the success rate of liquefaction evaluation given by the new formula at different probabilistic levels. As the probabilistic level increases from 10% to 90%, the total success rate for liquefaction evaluation first increases and then decreases; correspondingly, the success rate for liquefied sites decreases from 100% to 19.8%, and nonliquefied sites increase from 34.7% to 96.8%, respectively. It can be concluded from equation (25) that the critical shear wave velocity of the site becomes smaller with the increase of given probabilistic levels. Thus, as the probabilistic level increases, the success rate for liquefied sites calculated by the new formula decreases and the success rate for nonliquefied sites increases.

The probability of liquefaction for saturated ground can be divided into five categories, which is shown in Table 5. It indicates the greater the *P*_{L} of the site, the more likely the site will liquefy. Thus, a lower probabilistic level can be adopted in order to enhance the success rate for liquefied sites. Important projects related to national security and buildings which may have serious secondary disasters induced by earthquake require special fortification. It can be seen from Figure 3 that the new formula at 10% probabilistic level gives a satisfying success rate for liquefied sites. Thus, the formula with a probabilistic level equal to 0–10% can be proposed for important engineering facilities, such as nuclear power plants and the three Gorges Dam in China. When the liquefaction probability level is between 10% and 30%, the new formula can satisfactorily predict the liquefaction with a success rate greater than 90% and nonliquefaction with 35%, respectively. Lifeline buildings whose functions cannot be interrupted during an earthquake and should be recovered as soon as possible, such as hospitals and schools, can adopt the new formula with the probabilistic level ranging from 10% to 30% for seismic design. The new formula with a probabilistic level ranging from 30% to 60% can be used for civil buildings that will not cause large losses after an earthquake. The probabilistic level can be adjusted according to the importance of the projects.

##### 4.2. Reevaluation Based on Andrus Database

The proposed formula at a given 50% probabilistic level is compared with the Andrus method and Chinese code method. The 50% probabilistic level is chosen because the critical line of the liquefaction evaluation represents the same probability of liquefaction and nonliquefaction. When the probabilistic level is set to be 50%, the difference in the success rate between liquefied and nonliquefied sites given by the new formula is only 0.2%. The success rates of liquefied sites and nonliquefied sites both approximate 80%, which fulfill the requirement of the critical line of liquefaction evaluation. It is proved that the new formula at a probabilistic level of 50% is theoretically safe and reliable.

Table 6 presents the success rates of liquefaction evaluation by the new formula. In seismic intensity VII, the success rates for liquefied sites and nonliquefied sites are 78.3% and 89.1%, respectively. In seismic intensity VIII, the success rates for liquefied sites and nonliquefied sites are 65.7% and 80%, respectively. Liquefaction evaluation gives a satisfactory success rate of 94.7% at liquefied sites, while 50% at nonliquefied sites in seismic intensity IX. The results are slightly conservative due to the lack of actual data of nonliquefied sites in seismic intensity IX. In conclusion, the new formula overcomes the conservative disadvantage of the Chinese code method or the Andrus method for liquefaction evaluation.

As mentioned in Section 2.4, the Andrus method and Chinese code method tend to be conservative, which are not optimal. For the Andrus method and Chinese code method, the higher success rate of liquefaction evaluation comes at the cost of the low success rate of nonliquefied sites. Obviously, it is not advisable or beneficial to cost control in engineering design. Therefore, the balance of the success rate of the new formula needs to be achieved between the liquefied sites and nonliquefied sites. It is undoubted that the new formula at = 50% is more competitive than the Andrus method and Chinese code method because it can fulfill the balance requirement.

#### 5. Comparison of Liquefaction Evaluation with Different Methods

##### 5.1. Liquefaction Data from Kayen Database

Andrus et al. [16] have compiled a significant number of liquefied case histories and recommended deterministic boundary curves with different fine contents for liquefaction evaluation. A total of 415 groups of liquefaction databases based on shear wave velocity published by Kayen et al. [9] provide an opportunity to compare the new formula proposed in this study with other methods. It should be noted that 75 groups of data are repetitive between the Andrus database and Kayen database. In addition, 4 sites in the Kayen database are around the boundary line between liquefaction and nonliquefaction, which are removed from the database. Thus, the remaining 336 groups of data summarized in Table 7 are used to compare the performance of liquefaction evaluation between the new formula and other methods.

Figure 4 presents the data frequency distribution with the sandy-layer buried depth. Nearly 94% of sandy-layer buried depths at liquefied sites lie between 2 m and 10 m. Approximate 80% of sandy-layer buried depths range from 2 m to 8 m for nonliquefied sites. Figure 5 presents the data frequency distribution with ground water table. Statistically, most of the ground water table for both liquefied and nonliquefied sites lies in depths of 0 m to 4 m. Overall, the liquefied sites do not differ obviously from nonliquefied sites in the distribution of the sandy-layer buried depth and ground water table.

**(a)**

**(b)**

**(a)**

**(b)**

##### 5.2. Comparison with Different Methods

In the Kayen database, there are 4 groups of nonliquefied sites in seismic intensity VI and 3 groups of liquefied sites in seismic intensity X. All 7 groups of samples were successfully predicted by using the new formula. Chinese code method based on shear wave velocity only covers the areas of seismic intensity VII, VII, and IX. Thus, the samples in the seismic intensity VI and X are not discussed herein.

Table 8 compares the success rates of liquefaction evaluation based on the Kayen database in seismic intensity VII, VIII, and IX. In seismic intensity VII, the new formula gives a more dangerous evaluation result than the other two methods. The lower success rate at liquefied sites is mainly due to the lack of enough data in seismic intensity VII than other seismic intensities. However, in seismic intensity VIII and IX, the new formula can satisfactorily evaluate the liquefaction with success rates greater than 80% at liquefied sites. Moreover, in seismic intensity VIII or IX, the success rate obtained by the new formula at nonliquefied sites is far higher than that by the Chinese code method and Andrus method.

As displayed in Table 8, the total success rate at liquefied sites is 97.2% and 96.4%, corresponding to the Chinese code method and Andrus method, respectively. But the Chinese code method and Andrus method only give success rates of 23.5% and 15.4% for unliquefied sites, respectively. These two methods are obviously conservative. To be noted, the new formula proposed in this study gives the success rate of 79% and 67.9% for liquefied sites and nonliquefied sites, respectively. Therefore, considering the balance of the liquefaction evaluation between liquefied and nonliquefied sites, it is undoubtedly that the new formula proposed is more advantageous than the Chinese code method and Andrus method.

#### 6. Conclusions

This study discussed the influence of several parameters on the liquefaction potential. Based on the logistic regression theory, a new formula was established to evaluate the site liquefaction. The new formula gave more balanced evaluation results between liquefied and nonliquefied sites by comparing with the Andrus method and Chinese code method. The main conclusions can be drawn as follows:(1)The liquefaction potential is explained by the five descriptive variables including , *k*, , , and . These parameters are selected to build the new formula based on the logistic regression method.(2)The new formula can evaluate the liquefaction at different probabilistic levels. The result at the probabilistic level of 50% is relatively balanced than other probabilistic levels based on 225 sets of *V*_{s} data from the Andrus database. The success rates of liquefied sites and nonliquefied sites given by the new formula are both above 80% and the total success rate is 80.1%, which can be applied to a general project. For another type of project, the probabilistic level can be adjusted based on the importance.(3)The proposed formula at a = 50% is verified by 336 case histories collected from the Kayen database. Three methods discussed in this study almost have the same success rate for liquefaction evaluation at all sites. However, the evaluating result obtained by the Chinese code method and the Andrus method is obviously conservative. The proposed formula predicts a success rate of 67.9% for nonliquefied sites, which is far higher than that of Chinese code method and Andrus method.

In summary, the evaluating results gained by the Chinese code method and Andrus method tend to be conservative in different seismic intensities. The new formula overcomes the disadvantage and gives a more balanced evaluation result.

#### Data Availability

The data of the Andrus database used to build the new formula are on pages 113–121 of the literature [16] (https://nehrpsearch.nist.gov/static/files/NIST/PB99117897.pdf). And the data of the Kayen database collected from Kayen et al. [9] are available online in the ASCE Library (https://ascelibrary.org/doi/abs/10.1061/%28ASCE%29GT.1943-5606.0000743).

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

#### Acknowledgments

The authors would like to express their thanks to the financial support for this study from the Project of the National Natural Science Foundation of China (Grant nos. 51678300 and 51378257), Natural Science Foundation of the Jiangsu Higher Education Institutions of China (Grant nos. 19KJB560015 and 18KJA560002), and Middle-Aged & Young Science Leaders of Qinglan Project of Universities in Jiangsu Province, China (Grant nos. QL20200203 and QL20210210).