Ecological Security Assessment Based on Remote Sensing and Landscape Ecology Model

The various ecological processes of human beings are not only restricted by the landscape pattern on the regional scale but also affect the local and regional landscape together with global climate change. To date, most of the research on ecological security is based on the pressure-state-response (PSR) model, while there were a few studies based on the landscape ecology model approach. In addition, there has been little literature focus on the dynamic change process of ecological security, especially the simulation and prediction of the future development trend of ecological security. The purpose of this research is to establish a landscape ecological security evaluation method based on grid division, be aimed at breaking the inherent drawbacks of the administrative region as a unit mode approach, anticipated to better reflect the landscape ecological security status of the study area. A complex framework was constructed by integrating random forest algorithm, Fishnet model, landscape ecology model, and CA-Markov model. Multitemporal remote sensing data were selected as a data source, and land use maps of the study area were obtained through the random forest machine learning algorithm firstly. And then, the study area is divided into 307 grids of 2 km × 2 km using the Fishnet model. Next, the landscape disturbance index, landscape vulnerability index, and landscape loss index are used on the grid scale to establish a landscape ecological security evaluation model. Finally, ecological security assessment of Zhengzhou city was carried out, and the distribution map of the landscape ecological status in 1986, 1996, 2006, 2016, and predicted for 2026 was obtained. The results of the study showed that, as time goes by, the areas with high ecological safety gradually decrease. It is predicted that by 2026, the ecological security level of Zhengzhou will be dominated by lower ecological security areas. The research results can provide basic information and decision support for government agencies and land use planners to ensure responsible and sustainable development of the urban environment and ecology.


Introduction
Landscape ecology evaluates the ecological security status of a certain area through the landscape index and landscape pattern, which can provide hierarchical and integrated information data for the multiscale study of ecosystem functions [1]. Landscape pattern is the key content of landscape ecology research [2]. Landscape pattern and its corresponding changes are a comprehensive reflection of the ecological environment in a certain area due to the interaction of vari-ous factors such as natural and human factors. At the meantime, landscape pattern affects the regional ecological process [3] as well for the type, shape, size, quantity, and spatial combination of landscape patches are the interaction of various disturbing factors. Methods of landscape pattern index, which can reflect the composition and characteristics of the landscape in the spatial structure, can be used to quantitatively describe the change of landscape pattern and spatial heterogeneity [4]. Particularly, changes in landscape pattern will affect the ecological environment to a certain extent, thereby affecting ecological security. Researches on regional landscape pattern can effectively reveal ecological status and its corresponding spatial changes. Constructing an ecological security evaluation and analytics model by quantitative evaluation of regional landscape ecological index has certain advantages in the research of regional ecosystem security [5,6].
In the past decades, more attentions have been given to ecological security assessment, and a variety of methods have been developed to evaluate ecological security, such as ecological footprint [7,8], PSR models [9][10][11], and landscape ecology models [12]; for example, Sun et al. [13] used the PSR model combined methods of AHP and FCE to assess coastal wetland ecosystem health. Yang et al. [8] further improved the ecological footprint method to evaluate ecological security, and Peng et al. [14] took Liaoyuan City as a research case to discuss the possible ecological risks of landscapes at the county scale and the resulting ecological losses. Although different ecological safety assessment methods have different focuses, most of them are based on the PSR model. While the PSR model is applicable to regional ecological evaluations with small spatial scales with a few spatial variabilities and a few influencing factors, it is not applicable to natural disasters where human activities exceed the carrying capacity of the natural environment, and it is impossible to determine ecological security risks and uncertain threats. In addition, the logical connection between pressure and state is not inevitable; it is affected by a variety of complex factors. And the PSR evaluation system overemphasizes the impact of human activities on the ecological environment, and the evaluation results are somewhat one-sided [15]. Compared with the PSR model, the landscape ecology model starts from the ecosystem itself and is based on the structure, function, landscape patch dynamics and ecological succession of the landscape, the integrity and stability of the ecosystem, anti-interference ability, and resilience; the establishment of an evaluation index system can better reflect the temporal and spatial changes of landscape patterns. The landscape ecology model is not only suitable for larger spatial scale landscape design and land use planning but also suitable for smaller scale ecological evaluation [16].
In this research, we propose an improved landscape ecology ecological security evaluation model which can break through the long-term use of administrative regions as the research unit for data collection and analysis calculation [17,18] and can better reflect the landscape ecological security status of the research area [19]. For this purpose, multitemporal remote sensing images are used as a major data source, and a random forest algorithm machine learning algorithm is used to interpret basic information of ecological security evaluation. Next, the landscape disturbance index, vulnerability index, and loss index are integrated to establish the landscape safety index. Then, multitemporal ecological security distribution maps are obtained. Finally, with the help of the CA-Markov model, ecological security trend in future is simulated and predicted. We test the method in Zhengzhou urban area, the core city of Chinese urban concentration. The innovation of this method is as follows: (1) all the indicators used in the research are readily calculated from remote sensing satellite images; (2) the research obtained spatial distribution status of ecological security as well as the dynamic changes of ecological security; (3) simulate and predict the future development trend of ecological security; (4) break the solid-state model of data collection, analysis calculation based on administrative regions as research units in previous research.
The results will be significant to guide regional ecological planning, protection, and development and will be helpful to optimize the allocation of resources, to control the further deterioration of the ecology, and to realize regional sustainable development services [20]. The rest part of the paper is organized as follows: a brief introduction about research area, preliminary work, and preprocessing of datasets is given in Section 2. The proposed research framework based on remote sensing (RS) is drawn and applied to the research area in Section 3. The results are shown in Section 4, and analysis and discussion are given in Section 5. Finally, the conclusions are drawn in Section 6.

Study Area and Datasets
2.1. The General Situation of the Study Area. The study site is located in Zhengzhou city (Figure 1), the capital of Henan Province in the central part of the P.R. China, with a total area of 7507 km 2 as well as a population of 9878000 inhabitants. Zhengzhou is one of the National Central Cities in China and serves as the political, economic, technological, and educational center of the province, as well as a major transportation hub in China. The annual average temperature of the city is 14.5°C, and the general terrain trend is tilt from southwest to northeast. The study area is undergoing the accelerating of Chinese agglomeration, economic development, and urban expansion.
Zhengzhou city core area is located in the southern part of the North China Plain, the lower reaches of the Yellow River, and the north of central Henan Province, between 112°42′-114°14′ east longitude and 34°16 ′ -34°58 ′ north latitude. The city is 135-143 km long from east to west. The north-south width is 70-78 km, and the total area of the city is about 7446 km 2 . It has jurisdiction over 6 districts and 5 cities, 2 national-level development zones, and one national-level export processing zone. Zhengzhou city belongs to the warm temperate continental monsoon climate with four distinct seasons: spring, dry and drought; summers are hot, concentrated rainfall; autumn is cool, time is short; long winter, dry and cold, rain and snow rarely. The annual average temperature is 15.6°C, July is the hottest, and the monthly average temperature is 27.3°C. January is the coldest month, with an average monthly temperature of -0.2°C, and most of the rainfall is concentrated in July, August, and September. This paper selects the urban area of Zhengzhou city, including Huiji District, Jinshui District, Zhongyuan District, Guancheng District, and Erqi District as the study area ( Figure 1).  [21][22][23]. The prediction result of the random forest algorithm is obtained by solving the mode of the decision results of multiple decision trees that make up the forest. Each time, about 2/3 of the samples from the original training dataset were extracted with bootstrap sampling technology for decision tree model training in the corresponding epoch. So some samples in the original training set may be used to perform the decision tree at the same time. The remaining about 1/3 of the samples are drawn for the out-of-bag error evaluation of the random forest model, and then, the prediction error rate of the model is calculated [24]. The basic idea of random forest classification [21,25] can be described as three steps: (1) extracting k samples from the original training set using the bootstrap method, and ensuring the percentage of each class in subsets is the same as the original training set; (2) building k decision tree models for k sample subsets to obtain k classification results; (3) voting all classification results to obtain a final classified map. The two most important parameters in these procedures are number of trees used to construct RF and number of features selected for each individual tree. The second parameter is much more important than the first one by literature review [26][27][28]. In this research, the number of features is setting as log 2d, where d is all features available in this research. The evaluation of land use classification accuracy is a test of the classification results. There are two commonly used methods. One is to compare the classification results with the actual land use classification data; the other is to use high-resolution remote sensing images to select a certain proportion of training samples to test the accuracy of the classification data. At present, the most commonly used method of classification accuracy evaluation is the confusion matrix method, which is an error matrix, which is expressed by the number of correct classifications and the number of misclassifications in the experimental samples. In the confusion matrix, the overall classification accuracy (OA) and Kappa coefficient are mainly used to verify the classification results and evaluate the accuracy [23,29].

Data Sources and
The overall accuracy represents the probability of the number of samples that are correctly classified in the overall sample size and can reflect the overall accuracy of the classification results. As shown in where OA is the total classification accuracy, q is the number of classification categories, n is the total number of samples, and n ii is the number of discriminative samples of the ith category.
The Kappa coefficient is the accuracy coefficient of the overall classification. It is a kind of ratio, which represents the ratio of error reduction between classification and completely random classification, as shown in where q is the number of classes, n represents the total number of considered pixel, n ii is the diagonal element of the confusion matrix, n i+ represents the marginal sum of the rows in the confusion matrix, and n +i represents the marginal sum of the columns in the confusion matrix.

Fishnet
Model. The Fishnet model refers to the grid division of the research area based on GIS technology, where the grid scale is used as the evaluation unit to analyze the ecological security status of a selected area [30,31]. With the continuous advancement of science and technology and the deepening of ecological safety research, the advantages of this method have become more and more obvious: it can solve the problems of data processing accuracy in space. As the smallest analysis unit is selected as a grid, the method facilitates the effective collection of relevant data and breaks the drawbacks of using only an administrative division unit as a research unit. In addition, grid division technology of the Fishnet model can treat each grid as an independent unit by establishing a grid and then calculate each grid through the establishment of an evaluation system. In order to better reflect the landscape ecological security status of the study area, we used the grid tool in ArcGIS to divide the study area into 307 grids of 2 km × 2 km (shown in Figure 2) to break the previous inherent pattern. The division of the grid provides a more reasonable scale for the assessment of spatial change of landscape ecological security [31]. After the grid is identified, a grid layer is used to batch 3 Journal of Sensors crop the land use classification data in each period to obtain corresponding grids. Land use and land cover maps of the small units are imported into the landscape pattern analysis software fragstats4.2 in batches mode for analysis. Then, landscape metrics of each grid in the study area can be obtained through calculation, and the calculated landscape metrics are used to construct the landscape ecological security index. Finally, the grid data is converted into a grid point file (located at the center of each grid), and the calculated landscape ecological safety index is assigned to the grid point file. During this process, spatial resolution is setting to 30 meters with the kriging interpolation method.

Analysis of Landscape Pattern Changes.
To analyze the changes of land use land cover landscape pattern, the ecological significance of various landscape pattern metrics in this study area was calculated [20]. By literature review [32,33], landscape pattern can be divided into (1) patch level, which reflects the structural characteristics of a single patch in the landscape and is also the basis for calculating other landscape level indicators; (2) class level, which reflects the respective structural characteristics of different patch types in the landscape; (3) landscape level, which reflects the overall structural characteristics of the landscape. According to these landscape pattern metrics, the landscape of the study area was analyzed from a global perspective. In order to explore the structural characteristics and evolution rule of landscape pattern in various periods in Zhengzhou city, representative indicators (shown in Table 1) are selected at class level and landscape level with fully understanding the ecological significance of each landscape pattern index and the focus of the landscape structure reflected.

Analysis of Landscape Ecological Security
4.2.1. Landscape Disturbance Index. The landscape disturbance index describes the degree to which different landscape ecological environments are affected by natural and human activity factors. Its mathematical model can be described as where W i is the landscape disturbance index; a, b, and c are the landscape fragmentation index, landscape separation index, and landscape dominance index, respectively; a, b, and c represent the weight value of constructing fragmentation, separation, and dominance index, respectively. According to previous research results, and referring to the expert scoring method, the three weights are assigned as a = 0:5, b = 0:3, and c = 0:2.
Within this model, the mathematical expression of the landscape fragmentation index C i is shown as where N i is the number of landscape patches of type i and A i is the total area of landscape patches of type i.  Figure 2: Fishnet grid division.

Journal of Sensors
The mathematical expression of the landscape separation index S i is shown as where P i is the relative coverage of the i landscape and A is the total landscape area. The mathematical expression of the landscape dominance index is shown as where L i is the relative density of the landscape. Related studies believe that the relative coverage has a greater contribution to the dominance of the landscape than the relative density, so the weights d and e are assigned values of 0.6 and 0.4.

Landscape Vulnerability Index.
Different landscape plays different roles in the ecosystem and has different resistance to external pressure and sensitivity to external changes at the meantime. In research, this difference is usually associated with the vulnerability of the landscape, combined with previous research results [41]; on the basis of the actual situation of Zhengzhou city, values were assigned to corresponding land use types. For instance, unused land has the strongest vulnerability, and unused land is 5; the fragility of water is less than unused land, and water is set as 4; construction land is relatively stable which is set as -1; and the vulnerability index of cultivated land and vegetation are 3 and 2, respectively. Due to the different dimensions of each index, the normalization of the assignment of land use type is used as the vulnerability index of this type of landscape. The calculation method is shown as where L i represents the score assigned to the i type of land use and V i represents the landscape vulnerability of the i type of land use.

Landscape Loss Index.
The two variables of the landscape disturbance index and vulnerability index are superimposed and integrated to reflect the degree of loss of natural attributes of a certain landscape type after external disturbance. The calculation method is shown as where the calculation method and meaning of W i and V i are landscape disturbance index and landscape vulnerability index, respectively.

Landscape Safety
Index. The landscape safety index is expressed using two variables, the landscape loss index and the landscape area. The calculation method is shown in where LES e represents the landscape ecological security value of the eth sampling area, e is the number of grids, A e is the area of the eth grid, A ei is the area of the ith landscape in the eth sampling area, and m is the number of land use types.

CA-Markov
Model. The Markov model is a long-term prediction method based on situation and probability of a certain event at a certain moment, using the Markov chain, to predict its changes at various moments (or periods) in the future [42]. The CA model [43] has a capacity for spatial and temporal change simulation which can be defined in equation (10). The key factor of the Markov model is determining the transition probability of an event. The main components of the cellular automata model (CA) are units, states, rules, and neighborhoods. Each unit has one of a finite number of states [5,42,44]. The state of the unit at a certain moment depends on the current and previous moments. All unit states can be updated at the same time according to the transition rules. The principle is that the state of a cell at the next moment is a function of its neighborhood state where S ðt+1Þ and S ðtÞ are the set of cell states at t + 1 and t, respectively; f is the transition rule function of the state of the space cell; N is the cell field. t and t + 1 represent start and end times. The Markov model [45] has been widely used for trend simulation at various scenarios and can be expressed as where S ðt+1Þ and S ðtÞ are the status at time t + 1 and t, respectively.  Figure 3: LULC classification results of Zhengzhou city from 1986 to 2016.  [46]. And when combined with the CA model, the CA-Markov model can overcome the limitation that the Markov model failed to catch the spatial distribution in the future and can be used to simulate spatiotemporal changes. In this research, LULC maps obtained with RF in 1986, 1996, 2006, and 2016 were selected as the study years to calculate the transition area matrix using the Markov model. And a standard 5 by 5 contiguity filter was selected for the CA model which means the condition of future pixel is not only decided by information from previous state but also considered by corresponding surrounding pixels.

LULC Classification Maps with Random Forest.
In this research, training samples were selected by an expert. And the separability of samples was also checked to verify which is greater than 1.8 to meet the basic requirements for classification. The selected four RS images of Zhengzhou city were input to the RF classifier model to obtain classification maps of corresponding date over Zhengzhou city (Figure 3). Verification samples were selected from Google Earth to verify the accuracy of the classification results. And the total classification accuracy (OA) and Kappa coefficient were assessed with the help of selected ground truth. The total classification accuracy of the four images is greater than 80% (more details can be found in Table 2) which meet the demand of further process such as landscape pattern analysis.     Table 4). Results of the PLAND index showed that the dominated land cover types are farmland and built-up area in Zhengzhou city. In the period of 1986, the rank of PLAND in descending order is vegetation, farmland, built-up area, unused land, and water body. While in the period of 2016, the maximum PLAND index is builtup area and followed by farmland. The cause of these changes can be ascribed to two aspects. One is that, as the core area of national food production, governments of Henan province at all levels and the Department of Land and Resources have comprehensively strengthened the goal of arable land protection, which has further promoted the conservation and intensive use of land. The other one is that developing the economy, various provinces and cities in Henan have accelerated reclamation and at the same time high-standard farmland construction to ensure the amount of arable land; on the other hand, the rapid economic development in the urban area of Zhengzhou has resulted in continuous urban construction land expand. The PLAND index of vegetation has changed significantly, from 32.231 in 1986 to 7.852 in 2016, indicating that this type of land use type accounts for a large reduction in the overall landscape area.

Predicted
Particularly, by prediction in 2026, there will be almost no changes in the area of water and unused land; the PLAND value of vegetation has increased, the area occupied by builtup area has further increased, and the PLAND index of farmland has decreased significantly. This result indicates regional urbanization and industrialization. During the process, the farmland was severely disturbed by human activities, and the farmland was squeezed by large-scale development activities by humans.
For the LPI index, the value of water has decreased greatly from 2.8 in 1986 to 0.48 in 2016, which indicated that the largest patch area of the water area has decreased,  Journal of Sensors gradually shrinking, and the distribution is gradually fragmented; the LPI index of vegetation was 0.029 by 2016; this showed that the maximum patch area of vegetation first increases and then decreases, and its fragmentation degree first relaxes and then increases year by year. In addition, the LPI index value of built-up area increased greatly during the period from 1986 to 2016, which to a certain extent reflects that the urban construction of Zhengzhou city is more concentrated and the area of urban contiguous areas has increased significantly. The LPI index of farmland decreased from 13.03 in 1986 to 3.29 in 1996 and then increased to 7.26 in 2016, indicating that the largest patch area of farmland is fluctuating. According to the prediction results, the increase in the LPI value of water and vegetation indicates that the degree of fragmentation has decreased; the maximum patch index of unused land and farmland has decreased, and the distribution of farmland has become more concentrated.
From the perspective of the LSI index, the value of water is generally lower than that of other land types, mainly because its spatial distribution is relatively scattered, and the degree of aggregation is lower than that of other land types. During the study period, the PAFRAC index of various land use types in the urban area of Zhengzhou showed that vegetation was always the landscape with the largest PAFRAC index, indicating that the shape of vegetation in such patches was the most complicated. The PAFRAC index of built-up land has declined in fluctuations, indicating that the shape of this landscape type tends to be simple, and the complexity is reduced. The PAFRAC of water, farmland, and unused land have changed, showing an increasing trend, indicating that the shapes of these three landscape types are   Journal of Sensors becoming more complicated. It is estimated by 2026, the LSI index of land use types other than built-up land and vegetation will decrease, indicating that the degree of aggregation is reduced; the PAFRAC of unused land and farmland has increased, indicating that the complexity of patch shape has increased, and its shape has become more and more irregular.
In terms of the IJI index, vegetation and unused land showed a downward trend, indicating that the adjacent landscape types of these two land types decreased year by year. The IJI index value of water decreased from 83.53 in 1986 to 67.01 in 1996 and increased to 78.85 in 2016 and then became relatively stable. By 2016, the IJI index of built-up land increased significantly compared to 1986, mainly due to the mutual transformation of land types around built-up land and the increase in the types of adjacent landscapes. It is estimated that by 2026, compared with 2016, the spread of land use types other than built-up land and the parallel index will change slightly, indicating that the adjacent patch types of land types other than built-up land have decreased.

Changes of Landscape Index Landscape Level.
From the analysis of the landscape pattern in the landscape level, the calculation results of landscape metrics in different periods of Zhengzhou city are shown in Table 5.
The CONTAG index of Zhengzhou city increased from 31.02 in 1986 to 38.21 in 2016. It can be seen that the concentration of different landscape types in the study area has increased. During the study period, the IJI index of Zhengzhou city showed a decreasing trend, from 81.84 in 1986 to 74.97 in 2006, which indicated that the number of adjacent patches of a landscape type has decreased. However, the probability of all patch types adjacent to each other has stabilized after 2006. In terms of the Shannon diversity index (SHDI), whose value decreased in fluctuation, it indicated that the landscape types in the study area are gradually decreasing and the landscape diversity has decreased. In terms of the SHEI index, the uniformity of the distribution of various landscape types showed a decreasing-increasingdecreasing trend, which decreased from 0.94 in 1986 to 0.89 in 1996 and then increased to 0.95 in 2006. And SHEI increased to 0.7747 in 2016, which indicated that there were no particularly obvious advantages in the landscape in 1986, 1996, and 2006, while the advantages of landscape were obvious in 2016. According to the predicted results in 2026, SHDI and SHEI of the study area will be increased, which indicated that the intensity of human disturbance to the landscape has increased. And the predicted result also showed that the intensity of human disturbance to the landscape will increase, the diversity of the landscape will increase, and the control effect of the landscape dominance category on the overall landscape will be weakened, different land use types are more concentrated, and the types of adjacent features were increased.

Landscape Ecological Security Analysis.
Due to the characteristics of spatial heterogeneity in ecological security assessment, in this study, the Fishnet tool was used to refine the units of the ecological security assessment. By calculating the corresponding landscape index in the constructed ecological security assessment model, the landscape ecological security index value of each grid is obtained, and this value is interpolated to the center of the grid using the kriging interpolation method to obtain the ecological security distribution maps.
Furthermore, according to the natural discontinuity classification method [47], the ecological security of the study area is divided into five levels: low ecological safety area, lower ecological safety area, medium ecological safety area, higher ecological safety area, and high ecological safety area ( Figure 5).
From the characteristics of spatial distribution of LULC types, the high, higher, and medium security areas are mainly distributed in vegetation land cover area in the north and south, and the low and lower safety areas are distributed in the central built-up area and unused land areas.
By quantitative analysis ( Figure 6 and Table 6), in 1986, the ecological security of Zhengzhou city was dominated by high-safety, higher-safety, and medium-safety areas, and the low-security zones were less distributed, indicating that the overall ecological security level of Zhengzhou city during that period was relatively high. Compared ecological security of 1996 with 1986, the medium and higher ecological security area status of ecological security in 1996 increased, while lower and low ecological security areas are almost stable. And the location of the low ecological safety zone has hardly changed when compared with 1986, and it is roughly in the center of Zhengzhou city. Compared with the ecological of 2006 to 1996, the lower and medium ecological safety zones increased significantly in 2006, and the medium ecological safety zones were widely distributed, and the high ecological safety zones decreased. This was mainly because with the rapid development of urbanization, the urban construction land expanded greatly; by 2016, the area of urban construction land will reach to 40%, frequent human activities have intensified the degree of land use, and ecological security is in a low state. Compared ecological security of 2016 with The most obvious area is still located in the middle of Zhengzhou city, indicating that the urbanization of Zhengzhou city is developing rapidly, the population is growing rapidly, and human activities are frequent. Severe expansion has put a lot of pressure on the regional ecological environment, leading to a decrease in its ecological security. It is predicted that by 2026, the proportion of lower ecological safety areas will be the largest. Compared with 2016, the area of low ecological safety areas in this region has not increased much, the medium ecological safety areas have dropped by about 10%, and the percentage of higher ecological safety areas has declined, but the proportion of high ecological safety areas has increased slightly.

Conclusion
In this research, we have proposed an ecological security assessment framework based on random forest algorithm, Fishnet model, landscape ecology model, and CA-Markov model. With the support of multitemporal remote sensing and GIS technology, ecological security of research area in 1986,1996,2006,2016, and for predicted 2026 were calculated and simulated. The research results showed that due to the significant urbanization process and certain problems in the economic development model, the largest patch area of the water has been significantly reduced, from 2.78 to 0.48. The fragmentation of water is increasing year by year. However, the construction land is more concentrated, the area of urban contiguous areas has increased significantly, and the spread and juxtaposition index value of construction land has increased from 53.8871 to 83.5898. It is predicted that by 2026, the shape of unused land and farmland will become more irregular, and the area of construction land will further increase and the distribution will be more concentrated. The results of landscape ecological security evaluation showed that between 1986 and 1996, the dominated security level is high grade, while between 2006 and 2016, the area of medium ecological security zone was in major. It is predicted that the landscape ecological security situation of Zhengzhou city will continue to decline in 2026, but the percentage of high ecological security areas will increase slightly. The research results can provide basic information and decision support for government agencies and land use planners to ensure responsible and sustainable development of the urban environment and ecology. And the novelty of this paper is in the direction of successfully implementing the landscape ecology model for urban security assessment.

Data Availability
Researchers who are interested in this code can download it from hyperlink https://pan.baidu.com/s/ 1YeGDwd4DoOEkZpOh4vswKg, using password xdrg, using password i98q, or contact the corresponding author to obtain source data to conduct secondary analysis.