Causal Analysis of Ecological Impairment in Land Ecosystem on a Regional Scale: Applied to a Mining City Daye, China

We adopted a weight of evidence approach to establish a causal analysis of an impaired land ecosystem on a regional scale; namely, Daye, a traditional mining city in China. Working processes, including problem statements, a list of candidate causes, and a conceptual model were developed to represent a causal hypothesis for describing land degradation. Causal criteria were applied to integrate multiple lines of evidence. Then, various pieces of evidence were scored to either strengthen or weaken our causal assumptions. Results showed that habitat alteration, heavy metal accumulation, organic pollutants, water eutrophication, and nutrient runoff were the probable causes of land ecosystem impairment in Daye. Meanwhile, noxious gas, toxicants, altered underground runoff, atmospheric deposition, and acid rain were identified as possible causes. The most unlikely causes were altered hydrology, altered earth surface runoff, and soil erosion. Soil salinization, soluble inorganic salts, biological species invasion, and pathogens were deferred as delayed causes due to lack of adequate information. The causal analysis approach was applied to identify the primary causes of land degradation and implement accurate protective measures in an impaired land ecosystem.


Introduction
Land ecosystems are highly important life-support systems worldwide, and they play an essential role in ecological performance and socioeconomic development [1,2]. However, the unsustainable mining activities of human beings, including industrial land sprawl, the overpumping of underground water, the discharge of industrial wastewater and exhaust, mineral processing and smelting, large areas of solid waste piles and transportation, and deforestation, have led to extensive ecological degradation [3,4]. Examples of ecological degradation include soil erosion, desertification, soil pollution, the loss of natural habitats, and the sharp decline of arable land, which lead to land productivity decline, biodiversity reduction, human health and food safety problems, population overload, and ecological refugees [5][6][7].
For the mining city, land degradation has a long incubation period and ultimately leads to the serious unsustainable development of human beings [1,8]. Thus, a method that can accurately identify the causes of such a phenomenon should be developed. A causal analysis of an impaired land ecosystem is essential to determine state and degradation trends and ecological mechanisms, identify the causes of such a phenomenon, establish a mutual relation among land usage and high-quality risk assessment, and adopt efficient solutions for land ecological crises.
Establishing causal relationships between environmental stressors and observed effects on regional land ecosystems is difficult. This problem is due to various human disturbances and the intrinsic complexity of land ecosystems. Moreover, no widely accepted approach for determining such relationships is available. Most diagnoses of land ecological problems focus on various single issues, such as barren fields, drought, desertification, saline-alkali land, contaminated sites, wetland function degradation, and landscape diversity loss [9][10][11]. These mono-factor analyses are significant but are unable to represent the mechanisms of land ecosystem degradation on a region scale. In addition, scale control is an important factor in the causal analysis of regional land ecosystems. Land degradation is only reflected on a large scale, which represents the entire ecological function and health. For example, some ecological factors may be tolerable when they occur individually. However, when several factors simultaneously occur in a particular region, they will produce integrative ecological effects, such as the transfer of heavy metal contamination, the expansion of contaminated water, and extensive river blanking, which may overwhelm local land ecosystems [12][13][14]. Our understanding of ecological linkages and fate and transport pathways broaden as these problems expand [15,16].
The following methods and technologies have been used to conduct causal analysis research in land ecosystems on a larger scale. Laboratory experiments and field studies were combined to diagnose potential hazards [17,18]. Statistical analyses combined with a geographic information system (GIS) was used frequently for highlighting the distribution of the environmental degradation and risk sources, and endpoints [19]. A conceptual model was developed to present the complex relationships of ecological degradation [20,21]. Characteristic causal analysis was designed following the scientific assumptions of the discipline itself [22,23]. Multi-metric indexes were designed based on local key traces to identify the environmental stressors such as chemical, physical, biological, and socialeconomic [24][25][26] via the structural analysis of the whole system to simulate the causal mechanisms [27,28]. Lastly, various risk-based methodologies displayed the process of the function loss in the terrestrial ecosystem [29].
Current diagnosis methods overemphasize quantitative approaches, such as field sampling, experimental data analysis, statistical analysis, and model building [30]. However, such methods are costly depending on the required data and their quality, thereby limiting evidence selection. Consequently, ecologists experience difficulty in establishing risk hypotheses. In addition, causal analysis is accompanied with the development of ecological risk assessment. Numerous studies in this field have mostly focused on the causes of toxic and harmful substances. However, regional ecological risk assessments for large-scale areas are currently increasing, thereby indicating that entire ecosystem problems should be given considerable attention. Therefore, causal analysis methods must also keep up with changes in large-scale ecosystems. Moreover, qualitative approaches, including social survey, experience cognition, and expert assessment, are essential considering that land ecosystems are complex syntheses of nature and the social economy.
The weight of evidence (WOE) approach interprets lines of evidence for identifying complicated interactions, thereby providing assistance to assessors in verifying multiple hypotheses in complex ecosystems [31][32][33]. Numerous studies have detailed the principle of this method, whereas different causal analysis studies have focused on its developments [34][35][36].
The WOE method is extremely useful in guiding causal analysis procedures, particularly the manner by which various pieces of evidence are combined to identify stressors, such as chemical, physical, and biological ones, which cause the impairment of a complex ecosystem [35,37]. This approach embodies the combination of qualitative and quantitative methods and compensates for the disadvantages of single methods [33,38]. The WOE approach not only introduces quantitative methods for experiments, modeling, and statistical analyses, but also includes social investigations, such as interviews, questionnaires, and feedback channel construction. Simultaneously, the application of this framework is easy to enhance and update in different causal analysis research [39,40]. In addition, such application is highly important for follow-up work, including environmental investigations, database formation, risk communication, ecological modeling, simulation calculation, uncertainty problems, whether to promote preferred management actions, and support for legal opinions on ecological protection [16,41,42].
In this study, the WOE method was used in the causal assessment of an impaired land ecosystem in a mining area, i.e., Daye, a traditional mining city in China. Multiple lines of evidence were evaluated on the basis of the causal criteria for identifying probable causes. Our work highlights complicated ecological mechanisms, such as physical, chemical, and biological mechanisms, that cause ecological and human health damages, on a large-scale area. It can even help establish follow-up research, including effective ecological risk assessments, the identification of risk sources, decision making in eco-protection, and ecological restoration.  Figure 1). The area is rich in mineral deposits, with well-developed mining and metallurgy industries. The ecosystem in and around the city is wide-ranging, with lakes, rivers, forests, mines, arable lands, gardens, and urban and rural residential areas. The entire study area is 1566.3 km 2 .

Data Source
A total of 225 valid soil samples were identified throughout the study area in 2016. Soil degradation was investigated on the basis of the following indicators: soil quality (organic matter content), soil acidification (pH value), soil erosion (surface soil and water loss amounts), and soil texture (soil layer thickness, gravel content, and soil coarsening).
Our research group selected Cu, Pb, Cd, and As to determine their contents in the 225 soil samples. Atomic absorption spectrophotometry was performed in the laboratory to determine the contents of Cu, Pb, and Cd. Meanwhile, the barium chromate spectrophotometric method was used to identify content of As. The Nemerow index was applied to calculate the single values of the heavy metals, i.e., Cu, Pb, Cd, and As, and the Nemerow composite index ( Table 1). The aforementioned soil component measurements and heavy metal content data based on the geographical coordinates (Global Positioning

Data Source
A total of 225 valid soil samples were identified throughout the study area in 2016. Soil degradation was investigated on the basis of the following indicators: soil quality (organic matter content), soil acidification (pH value), soil erosion (surface soil and water loss amounts), and soil texture (soil layer thickness, gravel content, and soil coarsening).
Our research group selected Cu, Pb, Cd, and As to determine their contents in the 225 soil samples. Atomic absorption spectrophotometry was performed in the laboratory to determine the contents of Cu, Pb, and Cd. Meanwhile, the barium chromate spectrophotometric method was used to identify content of As. The Nemerow index was applied to calculate the single values of the heavy metals, i.e., Cu, Pb, Cd, and As, and the Nemerow composite index (Table 1). The aforementioned soil component measurements and heavy metal content data based on the geographical coordinates (Global Positioning System data) of the soil samples and the remote sensing image and elevation data of the study area were processed into spatial data and integrated into a geographic information system (GIS). Our research group obtained the quarterly monitoring data of water quality in water bodies located in 37 monitoring points, including major rivers, lakes, and sensitive waters, from 2013 to 2018 on the basis of interviews and a survey conducted by the Daye Environmental Protection Bureau. The group selected chemical oxygen demand, NH3-N, and water eutrophication as criteria for water quality deterioration. A similar GIS spatial processing was performed on the basis of the geographical information of the sampling points.
The research team, together with the Daye City Land Resources Bureau, frequently organized investigations on the ecological environmental damage caused by mining in Daye (a traditional mining city) from 2013 to 2018. Investigation sites included Cu, Fe, coal, Au, and Ag mines and related smelting sites, such as ore dressing sites, quarries, tailing reservoirs, coal gangue dumps, and open metal mines. Problems included over 800 damages to the Earth's surface and vegetation, 349 geological hazards (e.g., collapse, mine goaf, excavation, landslide, and water depletion), 550 land damages caused by solid waste dumping, and 1294 land plaques covering an area of 6943.58 ha. The research group conducted field investigations, plotting, and scene photography to obtain the information of each plot in terms of actual harm.
The data application of land use category was based on the land use maps of Daye in 2012. A substantial amount of data were obtained from the Daye Statistical Yearbook and social surveys, such as data on acid rain, particularly acid rain intensity released by the Daye Environmental Protection Bureau, and natural disasters collected from the historical information of Daye (1949-2016).

Causal Analysis Framework
Our work started with a plan of transformation, particularly the long-term ecological restoration, of resource-exhausted cities in the current Chinese government [43][44][45][46]. An operational framework was developed to serve as a guideline for identifying the primary causes for an impaired land ecosystem on a regional scale ( Figure 2). We adopted the WOE approach to establish causation: (1) existing serious problems of impaired land ecosystems were introduced through group seminars; (2) our database based on field investigations and literature accumulation was proposed; (3) our list of candidate causes was obtained from various hypotheses; (4) a conceptual model was developed to exhibit the causal relationships of multiple stressors that caused impairment; (5) a WOE approach was applied to the analysis of various lines of evidence, and such pieces of evidence were scored to determine whether such causation was supported or weakened; and (6) discussions and uncertainty analysis were conducted [47,48].
Land 2021, 10, x FOR PEER REVIEW 6 of 26 Figure 2. The procedure for regional causal analysis of land degradation in Daye area.

Problem Formulation
Problem formulation is a critical step in determining the evaluation range and target of ecological risk assessments [49,50]. Without this step, a problem will be exceedingly broad and will invite debates about its occurrence. Our study was motivated by the observation of serious problems about land degradation on a regional scale. Problem formulation was supported by field investigations, subsequent laboratory examinations, sufficient literature review, historical data provided by related local government departments, and current popular topics in environment fields. Seminars were conducted to refine the identified problems and focus on major ones. The following 11 problems were determined in the present study: (1) soil quality deterioration, (2) soil contamination, (3) water quality deterioration, (4) water shortage, (5) usable land shortage, (6) air pollution, (7) biodiversity decrease, (8) public safety and health threats (e.g., living condition deterioration), (9) land productivity decline, (10) eco-resiliency reduction (landscape heterogeneity loss), and (11) landscape aesthetic dysfunction. The identified problems reflect the current condition of land degradation in Daye City.

List Candidate Cause
After the problem statement was formulated, the next step was to establish a list of candidate causes that influence the stability and health of land ecosystems. A series of possible hypotheses should be effectively proposed to explain the manner by which these candidate causes lead to various observed negative ecological effects while diverging from the earlier hypothesis [51,52]. Candidate causes were determined in accordance with the aforementioned problems by a panel of experts (i.e., professors and local officials)

Problem Formulation
Problem formulation is a critical step in determining the evaluation range and target of ecological risk assessments [49,50]. Without this step, a problem will be exceedingly broad and will invite debates about its occurrence. Our study was motivated by the observation of serious problems about land degradation on a regional scale. Problem formulation was supported by field investigations, subsequent laboratory examinations, sufficient literature review, historical data provided by related local government departments, and current popular topics in environment fields. Seminars were conducted to refine the identified problems and focus on major ones. The following 11 problems were determined in the present study: (1) soil quality deterioration, (2) soil contamination, (3) water quality deterioration, (4) water shortage, (5) usable land shortage, (6) air pollution, (7) biodiversity decrease, (8) public safety and health threats (e.g., living condition deterioration), (9) land productivity decline, (10) eco-resiliency reduction (landscape heterogeneity loss), and (11) landscape aesthetic dysfunction. The identified problems reflect the current condition of land degradation in Daye City.

Causality Establishment
After the initial list of candidate causes was completed by the study group, a hypothesis conceptual model that included the candidate causes was developed to represent the complete profile of land degradation in the study area. Meanwhile, various sources in Daye City were identified through familiarity with the local economic development patterns, particularly those concerning major pollutant-emitting enterprises [7,31,53]. Candidate causes from different sources that result in various hazardous effects depict complex interactions that can help understand potential mechanisms and promote the refinement of our list of candidate causes [31,54] (Figure 3).

Causality Establishment
After the initial list of candidate causes was completed by the study group, a hypothesis conceptual model that included the candidate causes was developed to represent the complete profile of land degradation in the study area. Meanwhile, various sources in Daye City were identified through familiarity with the local economic development patterns, particularly those concerning major pollutant-emitting enterprises [7,31,53]. Candidate causes from different sources that result in various hazardous effects depict complex interactions that can help understand potential mechanisms and promote the refinement of our list of candidate causes [31,54] (Figure 3).

Causal Criteria
The WOE approach combined different types of evidence to determine whether an apparent association of causes and effects is actually true [36,38,55] This approach has been used in various applications related to environmental impacts on human health [56,57]. Then, USEPA developed this approach and frequently applied it to the causal assessments of various impaired ecosystems [16,58]. The identification result of the pathways of effects is also critical in determining the possible causes of any observed changes. Our study adopted Suter and Cormier's criteria [38,59], terrestrial systems by [16], and the USEPA causal analysis decision information system [54] for weighting evidence. We applied these mature standards to land ecosystems on a regional scale to make complicated ecological processes more practical and easier to understand. Meanwhile, various types of evidence were used to assess environmental causes of serious land degradation.

Weighting and Scoring of Criteria
The WOE approach combined different types of evidence to determine whether an apparent association of causes and effects is actually true [36,38,55] This approach has been used in various applications related to environmental impacts on human health [56,57]. Then, USEPA developed this approach and frequently applied it to the causal assessments of various impaired ecosystems [16,58]. The identification result of the pathways of effects is also critical in determining the possible causes of any observed changes. Our study adopted Suter and Cormier's criteria [38,59], terrestrial systems by [16], and the USEPA causal analysis decision information system [54] for weighting evidence. We applied these mature standards to land ecosystems on a regional scale to make complicated ecological processes more practical and easier to understand. Meanwhile, various types of evidence were used to assess environmental causes of serious land degradation.
This principle of causal criteria was primarily obtained from Suter [33,38], Cormier [59], and the USEPA method [54] for identifying the performance of various pieces of evidence.
The following are brief descriptions include 7 "types of evidence that use data from the case", and 3 "types of evidence that use data from elsewhere".
Types of evidence that use data from the case: Spatial/temporal co-occurrence-ecological effects must be observed at the same time that causes occur. For example, in an aquaculture area, the boom season of breeding is also the time algae outbreaks occur. Similarly, industrial emissions upstream will lead to water quality deterioration downstream.
Evidence of temporal sequence-the cause must precede the ecological effect. Biological population density and ecological function are always stable. With urban sprawl, considerable habitats are converted into construction land, thereby resulting in the decline of biodiversity and the loss of ecological function.
Exposure-response relationships from the field-the intensity or frequency of hazardous effects at a site increases with increasing levels of exposure to the causal agent or decreases with decreasing levels. For example, Pb concentration decreases with distance from the expressway.
Causal pathway-the data of the causes are limited, and the causal pathway can provide supplementary or surrogate evidence. For example, high-intensity irrigation is one of the major causes of soil salinization, and the application of pesticides poses high risks to human health through the food chain.
Evidence of exposure mechanism-relevant exposure of candidate cause is founded, such as a large number of industrial exhaust emissions triggering fog and haze weather through atmospheric deposition, while causing respiratory illness of local residents.
Manipulation of exposure-a cause is manipulated to determine whether it can decrease or increase its ecological effect by eliminating or altering exposure. The limited discharge of industrial sewage has improved water quality in polluted water.
Symptoms-a set of symptoms may be used to diagnose a particular cause if they are unique to that cause. For example, mineralized water and large-scale sulfur precipitation are common around mining areas.
Types of evidence that use data from elsewhere: Exposure-response relationships from other filed studies-at the diagnosed area, the cause must be at levels sufficient to cause similar hazards effects in other filed studies. For example, a huge amount of plastic mulch was abandoned in the damage sites with the widespread use of farm greenhouses. This scale of agricultural plastic waste per hectare has led to a serious "white pollution" in other similar areas, which causes long-term potential hazards due to the inability of soil to transform plastics.
Exposure-response relationships from laboratory studies-the candidate cause must be at levels associated with related hazardous effects in laboratory studies. For example, sewage irrigation brings great harm to agricultural production and is not easily detected. Through the toxicological analysis of the lab work, the concentration of the toxicants in the wasted water is much higher than the safety standard in the water, which causes a great risk to food safety.
Analogous causes-agents similar to the causal agent at the impaired site should lead to similar effects at other sites. For example, Almond eucalyptus was introduced to China from Australia in the early stage. It grows fast and has high economic value, especially for paper making or as raw materials for construction. It was once the preferred economic tree for farmers. Later, it was found that the root system developed to absorb a lot of nutrients and water in the land; the surrounding plants will not survive because of the lack of nutrients. Moreover, it will lead to the desertification of soil in planting areas, which will cause the deterioration of the forest environment.
These applications were similar to those in previous studies using WOE.

Analysis of Evidence
We performed an analysis to evaluate data from the case and elsewhere. Then, such evidence was used to support or weaken the hypothesis that causes certain damages [38,42,54,59]. Our group presented no details of candidate causes due to constraints in time and working conditions. The research team will continue searching for evidence to test our causal hypothesis in the subsequent ecological restoration of land ecosystems. This approach was selected to obtain a relatively complete causal analysis of land degradation on the regional level.

•
Habitat alteration The remote sensing data of Daye City from 1995 to 2015 indicate that the urban sprawl has been persistently popular, and large areas of water, forest, grassland, wetland, tidal flat, and paddy field are transformed into construction land [61][62][63]. The data of five-phase remote sensing interpretations from 1995 to 2013 indicated that the analysis of watershed conversion showed that the intensity of watershed to farmland was 3.85%, and the number of lakes decreased from 309 in 1995 to 84 in 2013. For example, the areas of Bao'an and Daye Lakes decreased by 23.67% and 16.25%, respectively, compared with their areas in 1995.

•
Heavy metal accumulation Mining and processing industries, such as ore dressing and washing and metal smelting and forging, cause the discharge of a substantial amount of untreated industrial wastewater, thereby resulting in serious heavy metal pollution in the surrounding soil and water environments [64,65]. Based on the National Soil Environmental Quality Standard (GB15618-1995) and the background value of soil heavy metals (the values of Cu, As, Cd, and Pb were adopted as 30.7, 12.3, 0.17, 26.7), the content of heavy metals in regional soil was evaluated using Nemerow index [66,67]. Interestingly, no metal deposit was found in Chengui Town. However, this town is highly contaminated with Cu, Pb, and As because it is the area of concentration of the rough ore processing industry given the convenient advantage of the main traffic road in Daye City. The heavy metal pollution phenomenon is complex.
The environmental protection department conducted sampling and testing of the water and sediment quality of Daye Lake due to the river pollution problem caused by local Cu mining. The water quality data show that Cu concentration ranges from 34 mg/L near the wastewater discharge of Cu mines to 18 mg/L in several downstream estuaries. That is more than 15,500 times the normal water area. Meanwhile, Cu pollutants present a significant gradient, and the polluted area exhibits an expanding trend [68,69]. In addition, the pH value of acidic Cu wastewater ranges from 1.8 to 2.2, which is a mixed solution that contains high concentrations of ferric sulfate, Cu, and Al. This solution enriches heavy metal ions in the river and seriously deteriorates water quality.

•
Organic pollutants Soil sampling in agricultural areas shows high concentrations of polycyclic aromatic hydrocarbons (PAHs) and polychlorinated benzene residues, thereby indicating the widespread use of pesticides and herbicides [70,71]. In traditional forest areas, the concentration of arsenide, which is the main component of insecticides applied every quarter of the year, is high. The field sampling of chemical fertilizer compost indicates that soil is enriched with phosphorus and Cd. This area is neither a mining nor an industrial processing area, but it is speculated to be affected by the use of a toxic phosphate fertilizer that contains Cd [72]. Statistics on the composition of agricultural fertilizers in various towns for several years are as follows: 45% nitrogen, 64.5% organic phosphate, 55.5% potassium, and 12% compound fertilizer.
The deterioration of water bodies caused by the application of pesticides and fertilizers is primarily due to rain-wash and seepage. The group investigated many lakes and pits around the cultivated land to search for evidence of organic pollution because organic pollutants are difficult to dissolve in water. The subsequent experiments indicate the impairment of the water environment by organic toxicants, particularly the identified three types of organic poisons: PAHs, organic chlorobenzene, and dichlorodiphenyltrichloroethane. The group conducted several on-site comparisons of wastewater discharge from paper mills and ammonia plants by focusing on the hazards posed to aquatic organisms by organic poison emissions [73,74]. Nearly no aquatic plants were found in the sewage outlet waters, and no evident aquatic plant was visible until 1800 m away [75,76].

• Water eutrophication
Daye City has extensive water areas. Thus, the large-scale development of traditional aquaculture has contributed to water quality deterioration. The application of a large amount of organic feeds and long-term high-phosphorus fertilizers caused the high phosphorus enrichment of the lakes [77,78]. This phenomenon is the major cause of the high eutrophication of a water body [79]. Statistics show that up to 70% of organic phosphate fertilizers are applied in towns, and the high phosphorus content exceeds standards. Moreover, the amount of excreta produced by livestock farming without strict environmental protection treatment is 318,600 tons/year. These waste materials enter water areas through artificial dumping and rainstorms.

• Nutrient runoff
A high-intensity farming model is supported by the government to ensure grain yield, which results in an ineffective supplementation of land fertility, and eventually leads to the barren condition of cultivated lands [80,81]. The four-phase strength indexes of reclamation in Daye City are 3.2, 3.0, 2.8, and 3.2, which are considerably higher than those in other counties and cities in the province. For soil samples in the study area, nitrogen, phosphorus, and potassium contents were substantially reduced, but Cl content was several times higher than those in other areas, thereby indicating the leaching effect.

•
Noxious gas A substantial amount of harmful gas, such as CO, NO, H2S, and SO2, is produced by the coal-burning process of the metallurgical industry. These gases disperse in the atmosphere around mining areas [82,83]. The spontaneous combustion of coal gangue piles is widespread, and 27 cases are confirmed [84]. We were eager to verify the harmful components of these irritating gases. However, verification was not performed due to the limited experimental conditions. Moreover, investigations on chronic disease and lung cancer cases were prohibited by the local health authorities. However, during the extensive visits made by our team, the residents around the mining areas generally claimed that air quality had deteriorated. Factory exhaust gases frequently irritated the nasal mucosa, and the frequency of acute respiratory diseases increased in winter.
• Toxicants Substantial amounts of mineral processing fluid and slag directly and indirectly flow into the surrounding soil and lakes through the illegal emissions of private metallurgy processing firms without qualified environmental equipment [85,86]. Such sewage and slag contain acute toxicants, mostly heavy metals and organic poisons [87]. In a spot test of several waste minerals in 2011-2015, the Environmental Protection Agency found that discharged wastewater contained highly toxic cyanide, up to 16.5 mg per liter of water, which exceeded the national safety standard by over 33 times.

•
Altered underground runoff The high intensity pumping of groundwater can decrease regional groundwater level and underground runoff changes and cause large-scale ground subsidence, collapse, and cracking [88,89]. For example, the groundwater level of one banned coal mine dropped from 19 m to nearly 55 m (i.e., a decline of nearly 46 m), and their water supply system was basically scrapped. At present, 136 irrigation wells and 25 pumping stations are available.

•
Atmospheric deposition Metal oxide dust emitted by the metallurgy industry enters the soil in the form of dust reduction through gravity, thereby forming a point-like pollution centered on the emitting plant with a radius of 2-4 km. Soil sampling analysis performed around the metallurgy industry belt indicates that the Pb content of the soil surface ranges from 147.65 to 3490.40 mg/kg. This finding shows a considerable decrease in gradient. During an investigation in an open mining area, a substantial amount of coal dust and other industrial dusts are produced through various processes, such as coal mining, blasting, crushing, loading, and transportation. Various cases of acute respiratory diseases and irritations are frequent among locals (whether toxic gas causes acute poisoning remains uncertain).

•
Acid rain Statistics show that the acid rain area is gradually expanding. The pH value of acid rain in Daye City is 2.1, which is under the severe category [90][91][92]. Metallurgy, petrification, thermal power, and mining are highly concentrated, and the sulfur dioxide emissions are 81,300 tons/yea. Thus, sulfur dioxide exceeds the national atmospheric standards and worsens acid rain pollution. The pH of soil samples is generally 2.7-3.2, with a serious depletion trend [93]. Notably, acid rain contributes to heavy metals entering the human body through the food chain, thereby posing a considerable health risk [94]. In addition, sulfur deposition in the study area exceeds 11-24% of the critical load, which seriously pollutes the surrounding soil environment [95].

• Altered hydrology
Large-scale engineering construction changes original habitats [96]. Meanwhile, water conservation and hydropower projects also exert various negative effects, such as hydrological changes, groundwater recession, and surface water shrinkage [97]. For the actual situation in Daye, some water conservation and hydropower projects selected for risk source investigation seem to be changing the hydrological environment.

•
Altered earth surface runoff The water supply of some important inland lakes is disturbed by several activities, such as the development of greenhouses, vegetable gardens, flowerbeds, and aquaculture. Over-retention may cause the shrinkage of rivers and lakes [98]. However, the 1995-2015 Water Resources Bureau yearbook survey indicates that the major rivers exhibit no cutoff period, and the lake water level presents no remarkable decrease during the period, except during extreme drought.

Soil erosion
Soil erosion has resulted in a decrease in the area of cultivated land, soil fertility and crop yield, and a serious contradiction between man and land [99]. The degree of soil erosion was calculated by field soil samples include soil moisture and soil surface thickness. According to national industry standard sl190-96 soil erosion classification and classification standard, the majority of the study area was in mild soil erosion. Meanwhile, compared with other mining cities, there is also serious surface damage, but it is not soil erosion [100].

• Soil salinization
Soil salinization is a land quality degradation process that involves soil salinity and alkalinity caused by human-made improper water diversion irrigation [101]. The groundwater level rises above the critical depth due to high-intensity irrigation. During evaporation, salt enriches the surface layer of the soil, thereby forming saline soil [101,102]. In 2009-2018, irrigation measures, such as wells and pumping stations, increased annually and intensified groundwater extraction, which is one of the major soil salinization risks. Unreasonable irrigation had not alarmed local authorities until widespread degradation occurred in the study area. •

Soluble inorganic salts
The concentrations of Fe and sulfur in the soil test are considerably high and eventually exhibit enrichment. In the study area, the severe deposition and migration of sulfur led to the death of soil microbes, nutrient loss, and large-scale land impoverishment [103,104]. However, sulfur-containing pollutants remain unknown. Accumulation in several areas, such as open mines, concentrated plants, washing ponds, tailing ponds, and coal gangue piles, is common throughout the year. A large amount of soluble inorganic salts from waste ore slags flowing into the surrounding area are determined using long-term weathering and rainfall leaching approaches [105]. In this case, mineral salts infiltrate into the groundwater layer, thereby forming highly saline water, which poses a safety hazard to residents and their livestock [106]. During groundwater sampling, acid testing with high 2-methyl-4chloropropionate was not conducted due to constraints in the experiment schedule.

• Biological species invasion
The continuous development of the aquaculture industry, particularly the large-scale reproduction of freshwater products, such as scorpionfish, squid, and Mandarin fish, exerts a serious impact on native aquatic species [107]. Grass cultivation with a certain degree of sugar content for fish feeding is extensive, and the original vegetation along the riverside has already disappeared. These conditions are likely to change the food chain in the original aquatic habitat. Vegetation changes in the southern mountain forest area have also elicited attention. The large-scale planting of economic trees, such as camphor, bamboo and cypress trees, is implemented because the self-contracted business model has been promoted by the local government for many years.

• Pathogens
A large area of vegetation and crops has withered, and the yellow-leaf disease frequently occurs. This condition is attributed to water, soil, and crops, which are highly susceptible to pathogens [108]. Another factor is the huge amount of sulfur dioxide in exhaust gas from intensive smelting enterprises. The area also receives heavy acid rain. However, emission data, particularly those of SO 2 and H 2 S, were unavailable from the Department of Environment. Hence, further work should be performed during the latter stage.
Subsequent evidence of candidate causes will be identified depending on future investigations and experimental analyses in our study region.

Identification of the Candidate Causes
A WOE table was constructed to interpret our evidence and the candidate causes were presented ( Table 2).
The probable cause has convincing evidence to indicate the cause of obvious damages, such as habitat alteration, heavy metal accumulation, and water eutrophication. For example, the analysis of evidence of habitat alteration was presented. A large number of wetland habitats have been surrounded by lake reclamation under the continuous pressures of the agricultural economy. This is obviously key evidence of temporal sequence. Moreover, increased agricultural water withdrawals have led to long-term insufficient water supply, thereby resulting in a large shrinkage area of local lakes and wetlands. This is evidence of spatial/temporal co-occurrence of cause and effect. These pieces of evidence were convincing support for the cause of habitat alteration mostly destroying the ecological structure.
Next, let us analyze the cause of heavy metal accumulation. In a soil sampling test, the heavy metal elements Cd, Cu, Pb, and As seriously exceeded the values of local soil background by 6.15, 28.76, 7.04, and 9.53 times, respectively [13]. It is not only evidence of spatial/temporal co-occurrence, but also exposure-response relationships from the field. Additionally, our research team repeatedly sampled river bottom sludge at different locations. The result shows that Cu concentration in the river sediments gradually decreases from west to east, with evident accumulation in the river channel (related to the river flow direction). Obviously, this is a causal characteristic of the exposureresponse relationships from the field. This finding suggests a potential pollution hazard to downstream water quality.
Then, the evidence of spatial/temporal co-occurrence for water eutrophication was presented. Notably, algae outbreak frequently occurs in large areas and seriously affects the aquatic ecosystem of the study area. Similarly, the presence of a set of symptoms also increases confidence of the alternative cause. Therefore, this convincingly supports the cause-effect that water eutrophication causes serious water deterioration.
The possible cause has some convincing pieces of evidence that imply the possibility of damage, such as atmospheric deposition and altered underground runoff. For example, the different evidence for the cause of atmospheric deposition was described as follows. Soil sampling analysis along the polluted industrial and mining zones indicates that dust contains automobile exhaust emissions, industrial dust, petroleum hydrocarbon pollution, and pollutant leakage on both sides of the expressway. The test data of soil heavy metals are given in the previous section of analysis of evidence. It is obvious that the effect occurs at the same time and in the same place as the candidate causes. Moreover, a clear Pb pollution zone has formed on both sides of the expressway, and Pb concentration in the surrounding area decreases as distance from the road increases. Strong effect gradients of exposure for atmospheric deposition were observed at spatially correlated sites. We also actively use the types of evidence that use data from elsewhere. Just as observation data throughout the years show that the surfaces of surrounding large-scale plants and field crops are covered with dust, and the leaves of plants commonly turn yellow and fall to the ground. Simultaneous comparisons were made with similar biological effect studies related to other mining cities.  The symbol "/" indicates that no score assigned, and "+/-" indicates disagreement among workshop participants. When consistency is high, the letter R indicates exclusion, and the letter D indicates assertion.
Then, we saw the candidate cause of altered underground runoff. We failed to obtain access to determine the specific engineering situation of mining companies. It certainly somewhat weakens the causality of the case. However, our team found pieces of evidence of abandoned mining. The groundwater level of one banned coal mine dropped. This is evidence of analogous causes that are similar to the causal agent in other sites/mines for a traditional mining city. Of course, we also have to use evidence from the types of evidence that use data from elsewhere, just like the following two pieces of evidence. One is in traditional Fe mining areas, where numerous mines pump out water all year round, thereby resulting in a decline in groundwater level. The other is irrigation, which is also one of the major factors that exert pressure on groundwater consumption learned, from the empirical literature [109].
The unlikely cause only has weak and non-persuasive evidence, such as altered hydrology and altered earth surface runoff. For example, consider the following two pieces of evidence spatial/temporal co-occurrence. No dry or broken flow in relevant rivers, significant change in water potential during rainstorm and flood seasons, and evidence of considerable changes in the hydrological pattern were not found. Ecological effects do not occur in the same place at the same time as candidate causes. Likewise, a substantial amount of water is consumed by numerous metallurgy industries because the water source is directly introduced from the Yangtze River, and no shortage of industrial water occurs. As mentioned earlier, Daye is along the Yangtze River system, which comprises vertical and horizontal lakes and abundant water resources. Such evidence indisputably suggested that there is no causal relationship.
The delayed cause lacks sufficient information/data to provide valid evidence at present, such as soil salinization and biological species invasion. The field survey did not diagnose soil salinization across the region. However, we calculated the irrigation intensity of each town on the basis of the statistics of the agricultural irrigation areas from 2001 to 2018. The results are higher than those of neighboring cities. This is explicitly an analogy using types of evidence that use data from elsewhere. It is generally recognized that a large amount of blind irrigation has caused serious soil salinization in many different countries.
In subsequent experiments, soil salinity analysis should be continuously conducted in our sampling point.
Additionally, there is little experimental data to support a causal link for the cause of biological species invasion on the larger scale. However, we still cannot ignore the emergence of some symptoms evidence, as follows: statistics from the local agriculture bureau indicate a sharp decline in the population of original endemic fish species, such as pygmy icefish; the use of various fish feeds leads to changes in the water environment [110]; on the one hand, the large-scale planting of economic trees; on the other hand, the distribution of original plant species sharply decreases, thereby resulting in the loss of forest diversity and even the lack of food for animals [111,112]. Therefore, we still could not conclude that the phenomenon of biological species invasions did not occur in the Daye area. At a later stage, these possibilities of causal effects must be studied, and new causal analysis repeated where appropriate.

Uncertainty Analysis
Our group discussed the uncertainties of this causal analysis research in the following aspects. (1) In terms of time, data are not from one time dimension. Some survey data were obtained in summer, whereas others were collected during the rainy season, and thus, the influence of season is considerable. For example, the study area belongs to the drainage system of the Yangtze River, and the distance is only 39 km. Therefore, runoff data change substantially with the flood season in the Yangtze River. (2) Data collection was limited. For example, water quality in 37 sampling sites does not cover all the water in Daye City. (3) Moreover, sensitive data, such as epidemic data on local health, leakage accidents, serious poisoning events, and food safety checks, were not obtained because permission from relevant departments was denied. (4) Although we were strict with regard to standards, three universities collaborated to perform our experiments, which may lead to errors in the experimental analysis. In particular, the level of ecological knowledge among different research teams varies. These reasons lead to uncertainties in our causal analysis research, and many hypotheses have not been proven by strong evidence. In summary, the application of the WOE method to causal analysis may be as quantitative as possible, logical, fair, convenient, and transparent. This method can draw on a wide range of interdisciplinary expertise to overcome the characterization of complicated ecological mechanisms on a broad scale.

Discussion
With the emergent contradiction between the tight supply of living resources and rapid population growth, the exploitation of wetland resources has long been an indisputable fact. Over the years, the farming model of lake area reclamation has caused the water area to shrink sharply [113].
Mining, as an inherent model of regional economic development, has intensified long-term heavy metal accumulation in soil in terms of degree and scope in the study area. In the traditional Cu mining areas, the concentration of Cu in soil is highly cumulative, thereby causing serious threats to river ecology, residential areas, cultivated land, and orchards along riverbanks. This condition will cause growth inhibition and toxicity, which will result in land productivity decline. Heavy metals pose harm to human health through the food chain. Therefore, government officials should be highly vigilant about food safety in this region. In particular, the accumulation of heavy metals in surrounding areas, such as ore piles, tailings ponds, and coal gangue piles, was serious and exhibited an evident decreasing trend. This phenomenon is attributed to long-term rainwater leaching, which causes the accumulation of heavy metals in large areas of soil surrounding a solid waste pile.
Estuary water quality exceeded the standard of heavy metal pollution, and the sediment analysis showed that Cu exhibited serious accumulation [114,115]. River pollution caused by Cu mining should be mitigated because it poses serious risks to downstream drinking water and breeding. Environment departments should set up additional monitoring points for water pollution and constantly monitor Cu concentration.
Fertilizer, pesticide, and film applications are the major sources of nonpoint source pollution. Statistics on the amount of agricultural fertilizer, pesticide, and film applied in Daye in 2012-2018 indicate that application strength was considerably higher than those in surrounding towns. In recent years, the cultivation technology of plastic mulch ground covering has developed rapidly. Hence, a substantial number of films have been abandoned in the field due to poor management. Solving this problem within a short period is difficult, thereby resulting in "white pollution" of the soil [116,117]. Pesticide residues are not allowed in the spot check of agricultural products due to the limitation of the local administration. Similarly, other categories of data on human health hazards from pesticides or other toxicants, such as acute or chronic poisoning incidents in the study area, are unavailable. However, we still speculate that huge amounts of pesticide residues have entered crops, thereby causing considerable harm to human health through the food chain [118][119][120].
The change in phosphorus content results from the alteration of land use patterns. The serious deterioration of water quality in many lakes is closely related to the development of the aquaculture industry in Daye area. Large-scale aquaculture has resulted in the eutrophication of water quality, which affects its self-purification capacity. Low phosphorus content inhibits the growth of aquatic plants. Therefore, local regulators should strengthen the supervision of feed application and consider the use of biological fertilizers.
Daye is located along the Yangtze River system, which includes lakes and rivers, and is rich in water resources. However, extensive irrigation systems have led to an annual dramatic increase in water use for agriculture due to the continuous drainage of groundwater by the mining industry [121]. Although direct underground runoff data are difficult to obtain, our group still infers that a significant change is highly probable in underground runoff. Irrigation intensity is considerably higher than those of surrounding cities. Many wells have dried up due to the declining groundwater level, and irrigation systems have become useless.
Dust storms occur in different seasons and are difficult to prevent. The topsoil in the surrounding area exhibits high Pb concentration, which easily forms dust and causes air pollution under the action of wind [122,123]. Although the expressway provides convenience for the production, processing, and transportation of the mining industry, which promotes the development of the local industrial economy, it is also an important potential risk source. Therefore, local environmental protection departments should pay attention to the usual dust reduction process, particularly in areas polluted by heavy metals. Severe acid deposition is primarily caused by urban industrial and mining agglomerations, with high concentrations of sulfur dioxide and nitrogen oxides in the air [91,93]. Acid rain accelerates soil acidification and the loss of mineral nutrients, thereby ultimately rendering the soil barren and uncultivable. Moreover, water acidification exerts a cumulative effect, which dissolves toxic substances into water and causes the death of aquatic organisms [94].
In summary, the aggregation of the mining and metallurgical processes results in the long-term deterioration of water, soil, and atmospheric environments in the study area, thereby forming extensive land degradation zones with complex mechanisms. This phenomenon not only affects the growth of plants and animals and poses considerable threat to human health, but also causes irreparable losses to the ecosystem.

Conclusions
Our study provides a suggested guideline for investigating possible causal relationships between ecological stressors and their effects on land ecosystems. Multiple lines of evidence based on the application of the WOE approach are applied to the causal analysis to determine land degradation on a regional scale. Although not all environmental causes can be identified, and many assumptions still need to be tested further, this research is sufficiently useful in acquiring a considerable amount of information, understanding mechanistic processes, and forming various ecological hypotheses.
Simultaneously, this work provides an effective decision-making basis for the management and restoration of damaged land ecosystems. First, the causal analysis of land ecosystems has expanded the methods for regional land ecological investigation, which exhibit practical operational significance. Second, hypothesis design for the complex mechanisms of land degradation has also enhanced the understanding of the ecological process on the regional scale. Third, the advantage of combining qualitative and quantitative methods/analyses is that they not only capture key issues but also ensure the accuracy of quantitative risk research, including exposure response, evaluation, prediction, decision making, and uncertainty analysis, in the latter stage. Lastly, the design of this causal analysis approach has been widely recognized by different stakeholders of land resources as the basis of decision-making tools. The causal analysis established for an impaired land ecosystem has effectively identified major hazard causes and their action mechanisms on a regional scale. This finding contributes to the formation of a stable and healthy land use relationship. It is also crucial for ecological protection and the sustainable utilization of resources.
In future research work, attention should be paid to other heavy metal elements involved in urban soil pollution, such as the investigation and analysis of heavy metal Cr. In addition, too high sulfur concentration in soil samples is also worthy of attention, although its migration is still uncertain. At present, most urban residents drink water from the Yangtze River, but most rural residents still drink groundwater. However, groundwater mineralization causes a huge health hazard, which has not been confirmed. Moreover, the attention to rural residential areas is not enough, especially ignoring the health hazard of humans and livestock during land degradation. In the following work, a health survey will be included in the survey on regional land resources and environment. Additionally, it is necessary to pay attention to providing a legal basis for ecological damage caused by specific harmful ways.
Exposure-response analysis in the risk research on regional land degradation is also a difficult link. In the following work, it is necessary to continue to improve the causal investigation to further reveal the complex degradation mechanisms (physics, chemistry, biology and humanity) at the regional level. On this basis, a complete exposure-response analysis can better display the process of large-scale land degradation.
Currently, the research on tracing the causes of large-scale land degradation, identifying the risk sources and revealing the ecological mechanisms are mostly based on qualitative approaches, which are highly subjective. This requires us to design excellent experiments for attribution analysis. However, experimental design is difficult for the complex land habitat system, and it is not clear how much confidence will be obtained from the poorly designed experiment. Therefore, the next experiments should rely on the causal probability of statistical correlation, and we should also pay attention to the development of quantitative causal evaluation tools.