An assessment of the carrying capacity of groundwater resources in North China Plain region–Analysis of potential for development

Over-exploitation of groundwater in North China Plain (NCP) has resulted in a series of eco-environment problems. Sustainable use of groundwater resources in NCP, in particular management of groundwater resource carrying capacity (GRCC), faces an unprecedented challenge. Here we define GRCC, and a new assessment method is tentatively proposed and applied to evaluate GRCC based on the whole NCP, city administrative units and county administrative units. Our study divided the NCP into three zones, i.e. non-overexploited non-overloaded zone (NNZ), overexploited but non-overloaded zone (ONZ), and overexploited overloaded zone (OOZ). Results confirmed 27.6% of counties belonged to NNZ. However, 58.9% of counties and NCP as a whole belonged to ONZ, and 13.5% of counties belonged to OOZ. Spatially, NNZs were mainly distributed in Beijing, parts of eastern coastal cities and Henan Province. OOZs were mostly distributed in middle-eastern part of Cangzhou, parts of Dezhou, Tianjin and Binzhou, and the remaining areas belonged to ONZs. We suggest two approaches for enhancing GRCC, i) increasing the amount of available groundwater and ii) improving the water use efficiency. An increase of 11.0 billion cubic meters to the available groundwater levels combined with water use efficiency improvements up to 479 CNY per cubic meter of the world mean, the gross domestic product (GDP) sustained by groundwater in the NCP could reach 11.1 trillion CNY and maintain a 20 years of GDP development assuming the current rate of growth.     

Quaternary Aquifer of the North China Plain—assessing and achieving groundwater resource sustainability

The Quaternary Aquifer of the North China Plain is one of the worlds largest aquifer systems and supports an enormous exploitation of groundwater, which has reaped large socio-economic benefits in terms of grain production, farming employment and rural poverty alleviation, together with urban and industrial water-supply provision. Both population and economic activity have grown markedly in the past 25 years. Much of this has been heavily dependent upon groundwater resource development, which has encountered increasing difficulties in recent years primarily as a result of aquifer depletion and related phenomena. This paper focuses upon the hydrogeologic and socio-economic diagnosis of these groundwater resource issues, and identifies strategies to improve groundwater resource sustainability.

---

Resumen El acuífero cuaternario de la Llanura Septentrional de China es uno de los mayores sistemas acuíferos del mundo y soporta una enorme explotación de su agua subterránea, las cuales han originado grandes beneficios socioeconómicos en términos de producción de grano, empleo en agricultura y mitigación de la pobreza rural, además de proveer agua para abastecimiento urbano e industria. Tanto la población como la actividad económica han crecido mucho en los últimos 25 años con una gran dependencia de las aguas subterráneas, que ha encontrado dificultades recientes por la explotación intensiva del acuífero y fenómenos relacionados. Este artículo se centra en la diagnosis hidrogeológica y socioeconómica de los problemas relacionados con las aguas subterráneas e identifica estrategias para mejorar la sustentabilidad de este recurco.

---

Résumé Laquifère Quaternaire de la Plaine du Nord de la Chine est lun des plus grands systèmes aquifères du monde; il permet une exploitation énorme deau souterraine, qui a permis des très importants bénéfices socio-économiques en terme de production de céréales, demplois ruraux et de réduction de la pauvreté rurale, en même temps que lapprovisionnement en eau potable et pour lindustrie. La population comme lactivité économique ont remarquablement augmenté au cours de ces 25 dernières années. Elles ont été sous la forte dépendance du développement de la ressource en eau souterraine, qui a rencontré des difficultés croissantes ces dernières années, du fait du rabattement de laquifère et des phénomènes associés. Cet article est consacré aux diagnostiques hydrogéologique et socio-économique des retombées de cette ressource en eau souterraine; il identifie les stratégies pour améliorer la pérennité des ressources en eau souterraine.

     

Helium isotopic fingerprints of the heavy land subsidence left in groundwater of the Saga Plain,near Beppu–Shimabara graben,Kyushu, Japan

The Saga Plain is near Beppu–Shimabara graben, a region of potential active volcanism. In the graben, mantle He, which has a high ³He/⁴He ratio of 1.1 × 10⁻⁵, escapes easily from the underlying subduction zone. In groundwater of the Saga Plain, except in the Shiroishi district, this ratio gradually increased as the dissolved He content increased, to a maximum of 5 × 10⁻⁶. In central Shiroishi, however, the ratio reached a minimum of 8.7 × 10⁻⁷ with increasing dissolved He content, suggesting that groundwater in central Shiroishi has selectively accumulated radiogenic He, which has a very low ratio of 1 × 10⁻⁸, rather than reflecting the regional He, which is rich in mantle He. This can be explained if groundwater in Shiroishi has become mixed with fossil pore water drawn from impermeable marine clay aquitard layers. The withdrawal of pore water has also caused severe land subsidence in central Shiroishi.     

Spatiotemporal relation of RADAR-derived land subsidence with groundwater and seismicity in Semnan—Iran

Natural Hazards - Due to excessive harvesting of underground water resources in many important aquifers inside Iran, ground subsidence is occurring at different speeds. In present study, InSAR...     

Urban land use/land cover changes in the Tema metropolitan area (1990–2010)

Land use and land cover changes are local and place specific, occurring incrementally in ways that often escape our attention. This study sought to detect changes in land cover in the Tema Metropolis of Ghana from 1990 to 2010. Multispectral Landsat Thematic Mapper data sets of 1990, 2000 and 2007 were acquired, pre-processed and enhanced. Unsupervised classification of the images was performed and six land cover classes (water, wetlands, closed vegetation, open vegetation, cropped lands, and built-up) were derived. The post-classification change detection technique was performed to derive the changes in land cover and their corresponding change matrices. Between 1990 and 2010, built-up areas expanded steadily to become the most prevalent land cover type in the metropolis, reducing vegetation cover dramatically. High population growth with its attendant rise in the demand for housing, and increasing commercial activities, were found to have influenced land cover changes over the period.     

Dual impact on the groundwater aquifer in the Kazan Plain (Ankara, Turkey): sand–gravel mining and over-abstraction

This study explains the impact of sand–gravel mining and over-abstraction, and the response of the groundwater system in the Kazan Plain, Turkey. The plain used to be known for its fresh groundwater potential, valuable agricultural lands and natural beauty until the 1980s. According to the estimation in 1975, there was 15.5 million m³ annual useable freshwater in the Quaternary sand–gravel aquifer. Groundwater level ranged between 0 and 5 m from the surface of the plain in the 1970s. Because of the large and deep excavations by the sand–gravel pits during the past 25–30 years, the aquifer has thinned and removed entirely in some places. In addition, over-abstraction has accelerated the decline of the groundwater level, particularly in the middle and the upper part of the plain in recent years. As of 2009, about 12% of the total volume of the aquifer area removed by the pits and groundwater table has been reduced to between 5 and 20 m. The decline of the water table reaches 15–20 m in the regions where over-abstraction has taken place. To reduce the hazards to the groundwater system, the sand–gravel pits have to be banned immediately, a reclamation project applied and abstraction must be reduced.     

Review: Water–rock interactions and related eco-environmental effects in typical land subsidence zones of China

Land subsidence is common in some regions of China. Various eco-environmental problems have arisen due to changes in water–rock interactions in these subsided areas, for which a comprehensive understanding of the hydrogeological setting is needed. This paper presents the general status of land subsidence in three typical subsided areas of China through the compilation of relevant data, and reviews some typical changes in the water–rock interactions in subsided areas along with related eco-environmental issues. It is found that the subsidence development and distribution are controlled by the groundwater-withdrawal intensity externally, and by the thickness and compressibility of unconsolidated sediments internally. The physical changes and related effects of water–rock interactions in subsided areas include: (1) the decreased ground elevation that caused floods, waterlogged farmland, etc.; (2) the differential subsidence that caused ground fissures; and (3) the change of seepage field that caused substantial reduction of the water resource. Chemically, the changes and related effects of water–rock interactions include: (1) the change to the chemical environment or processes due to the hydrogeologic structure alteration, which caused groundwater pollution; and (2) hydrologic mixing (seawater intrusion, artificial recharge; exchange with adjacent aquifers or aquitards), which degraded the groundwater quality. Further research on the subsided areas in China is suggested to reveal the mechanisms regarding biological and gaseous (meteorological) changes from the perspective of interacting systems among water, rocks, biological agents and gases.     

Characteristics and Classification of the Geothermal Gradient in the Beijing–Tianjin–Hebei Plain,China

Mathematical Geosciences - Study of the characteristics and classification of geothermal gradients can effectively guide the exploration and development of geothermal resources. In this paper, we...     

Assessing the ‘wicked problems’ associated with the quality of groundwater used in irrigation: a case study from the North China Plain

Studies on the quality of groundwater have moved beyond the physical realm of contamination and purification, towards the economic concerns of choice and the management of the problem. With these approaches the complex biophysical processes are assessed from the users’ perspective and the policy outcomes that could be used to resolve the problems of groundwater contamination are evaluated. However, in a set of unrelated studies, it has been found that attempts by governments to resolve the problems of groundwater contamination in agriculture have a poor record of success. This could be because the problem is too extensive and diverse to handle or it could be a case of poor policy selection. Taking an example from the North China Plain to illustrate some of the major issues raised in this study, it is concluded that the problem itself is unresolvable on a large scale. In other words, groundwater contamination can be defined as a ‘wicked problem’, i.e. unresolvable by applying pure science, closely linked with social issues, and for which there are no optimal solutions. In this situation, the best solution is possibly to encourage farmers to live with and handle the problem as they best see fit.     

Monitoring of land subsidence and ground fissures in Xian,China 2005–2006: mapped by SAR interferometry

The City of Xian, China, has been experiencing significant land subsidence and ground fissure activities since 1960s, which have brought various severe geohazards including damages to buildings, bridges and other facilities. Monitoring of land subsidence and ground fissure activities can provide useful information for assessing the extent of, and mitigating such geohazards. In order to achieve robust Synthetic Aperture Radar Interferometry (InSAR) results, six interferometric pairs of Envisat ASAR data covering 2005–2006 are collected to analyze the InSAR processing errors firstly, such as temporal and spatial decorrelation error, external DEM error, atmospheric error and unwrapping error. Then the annual subsidence rate during 2005–2006 is calculated by weighted averaging two pairs of D-InSAR results with similar time spanning. Lastly, GPS measurements are applied to calibrate the InSAR results and centimeter precision is achieved. As for the ground fissure monitoring, five InSAR cross-sections are designed to demonstrate the relative subsidence difference across ground fissures. In conclusion, the final InSAR subsidence map during 2005–2006 shows four large subsidence zones in Xian hi-tech zones in western, eastern and southern suburbs of Xian City, among which two subsidence cones are newly detected and two ground fissures are deduced to be extended westward in Yuhuazhai subsidence cone. This study shows that the land subsidence and ground fissures are highly correlated spatially and temporally and both are correlated with hi-tech zone construction in Xian during the year of 2005–2006.     

Quantification of the water balance and hydrogeological processes of groundwater–lake interactions in the Pampa Plain, Argentina

This paper gives an account of the assessment and quantification of the water balance and the hydrogeological processes related to lake–groundwater interaction in the Pampa Plain by using hydrogeochemical, isotopic and flow numerical modeling techniques. La Salada is a permanent shallow lake, with an area of 5.8 km², located on the SE of Buenos Aires Province. A total of 29 lake water samples and 15 groundwater samples were collected for both hydrochemical analysis and environmental stable isotope determination. Water table depths were measured in wells closed to the lake. Groundwater samples appear grouped on the Local Meteoric Water Line, suggesting a well-mixed system and that rainfall is the main recharge source to the aquifer. Water evaporation process within La Salada is also corroborated by its isotopic composition. The model that best adjusts to La Salada Lake hydrochemical processes includes evaporation from groundwater, calcite precipitation with CO₂ release and cationic exchange. The annual water balance terms for the lake basin indicates for each hydrological component the following values: 1.16 E⁰⁸ m³ rainfall, 8.15 E⁰⁷ m³ evapotranspiration, 1.90 E⁰⁶ m³ runoff, 1.55 E⁰⁷ m³ groundwater recharge, 6.01 E⁰⁶ m³ groundwater discharge to the lake, 9.54 E⁰⁶ m³ groundwater discharge to the river, 5.00 E⁰⁵ m³ urban extraction and 4.90 E⁰⁶ m³ lake evaporation. Integrated analysis of hydrochemical and isotopic information helped to calibrate the groundwater flow model, to validate the conceptual model and to quantitatively assess the basin water balance.     

Hydrogeochemical characteristics of groundwater and pore-water and the paleoenvironmental evolution in the past 3.10 Ma in the Xiong’an New Area,North China

The groundwater level has been continuously decreasing due to climate change and long-time overexploitation in the Xiong’an New Area, North China, which caused the enhanced mixing of groundwater in different aquifers and significant changes in regional groundwater chemistry characteristics. In this study, groundwater and sediment pore-water in drilling cores obtained from a 600 m borehole were investigated to evaluate hydrogeochemical processes in shallow and deep aquifers and paleo-environmental evolution in the past ca. 3.10 Ma. Results showed that there was no obvious change overall in chemical composition along the direction of groundwater runoff, but different hydrochemical processes occurred in shallow and deep groundwater in the vertical direction. Shallow groundwater (< 150 m) in the Xiong’an New Area was characterized by high salinity (TDS > 1000 mg/L) and high concentrations of Mn and Fe, while deep groundwater had better water quality with lower salinity. The high TDS values mostly occurred in aquifers with depth < 70 m and >500 m below land surface. Water isotopes showed that aquifer pore-water mostly originated from meteoric water under the influence of evaporation, and aquitard pore-water belonged to Paleo meteoric water. In addition, the evolution of the paleoclimate since 3.10 Ma BP was reconstructed, and four climate periods were determined by the δ¹⁸O profiles of pore-water and sporopollen records from sediments at different depths. It can be inferred that the Quaternary Pleistocene (0.78–2.58 Ma BP) was dominated by the cold and dry climate of the glacial period, with three interglacial intervals of warm and humid climate. What’s more, this study demonstrates the possibilities of the applications of pore-water on the hydrogeochemical study and further supports the finding that pore-water could retain the feature of paleo-sedimentary water.© 2021 China Geology Editorial Office.     

Arsenic pollution of groundwater and its release mechanism in the central part of the Hetao Plain, Inner Mongolia, China

Arsenic pollution of groundwater has become a serious health problem in several regions of the world. Wuyuan County located in the central part of the Hetao Plain, Inner Mongolia, is one of the most severely arsenic polluted areas. Peoples of the county have used groundwater of shallow （12-25 m） tube-well as drinking water. The groundwater is Na-HCO3 and Na-Cl type on the piper plots and aligns with steady Ca ： Mg ratio on cation diagram. The arsenic content in the water is more than 0.2 mg/L and the highest value reaches up to 0.52 mg/L. Drilling data show that Holocene sediments are divided into three zones on the basis of their color appearances; brown, grey and dark gray zones. The strata from 2m to 9-14m under the surface are mainly composed of clay and silt which are impermeable layers. The grey and dark grey zones composed of sand underlie the impermeable bed and the condition of their zones is reductive. At the depth of 30 m near the bottom of drilling, ORP of the groundwater is -300- -400 mV. The dark grey sand changes to yellow in color under the atmosphere. Under reducing environment, the color returns to its original dark grey. The characteristic color is due to amorphous or low-crystalline Fe sulfide coating sand grains. The sediments are of fluvial origin by the Yellow River. The arsenic content in the sediments is rich in fine-grained part and its values are more than 20 mg/kg in clay layer. The highest is 38 mg/kg. X-ray diffraction and macroscopic observation show that the clay layer is made of minerals such as chlorite and illite and contains high quantities of organic matter.     

Geological controls and effects of floodplain asymmetry on river–groundwater interactions in the southeastern Coastal Plain, USA

Channel sediment and alluvial aquifer hydraulic properties exert a major control on river–groundwater interactions. Channels and floodplains are often asymmetrical, resulting in differences in sediment hydraulic properties across the river. Floodplain asymmetry is common along Coastal Plain rivers in South Carolina and North Carolina, USA. The Tar River, North Carolina, has an asymmetrical valley. The study objective was to characterize the effects of floodplain asymmetry and geological controls on river–groundwater interactions. Floodplain and river channel sediments adjacent to the river were characterized with split spoon cores and hand auger samples along a 22-km reach. Hydrogeology was characterized with 38 piezometers and water level recorders in and adjacent to the river. Ground penetrating radar was used to define the shallow stratigraphy. Channel sediments were significantly different between the north and south sides of the river. Hydraulic conductivity and groundwater inputs were greater on the side of the river (north) that contained more permeable fluvial deposits. Groundwater chemistry (δ¹⁸O, specific conductance) data also suggested greater exchange between surface water and groundwater on the north side of the river channel. A conceptual hydrogeological model illustrates that groundwater movement and contaminant transport to the river differs across the channel due to asymmetrical geology.     

Long-term spatio-temporal hydrochemical and 222Rn tracing to investigate groundwater flow and water–rock interaction in the Gran Sasso (central Italy) carbonate aquifer

In the Gran Sasso fissured carbonate aquifer (central Italy), a long-term (2001–2007) spatio-temporal hydrochemical and ²²²Rn tracing survey was performed with the goal to investigate groundwater flow and water–rock interaction. Analyses of the physico-chemical parameters, and comparisons of multichemical and characteristic ratios in space and time, and subsequent statistical analyses, permitted a characterisation of the hydrogeology. At the regional scale, groundwater flows from recharge areas to the springs located at the aquifer boundaries, with a gradual increase of mineralisation and temperature along its flowpaths. However, the parameters of each group of springs may significantly deviate from the regional trend owing to fast flows and to the geological setting of the discharge spring areas, as corroborated by statistical data. Along regional flowpaths, the effects of seasonal recharge and lowering of the water table clearly cause changes in ion concentrations over time. This conceptual model was validated by an analysis of the ²²²Rn content in groundwater. ²²²Rn content, for which temporal variability depends on seasonal fluctuations of the water table, local lithology and the fracture network at the spring discharge areas, was considered as a tracer of the final stages of groundwater flowpaths.     

Chemical and strontium isotopic characteristics of shallow groundwater in the Ordos Desert Plateau,North China: Implications for the dissolved Sr source and water–rock interactions

In this study, the chemical and Sr isotopic compositions of shallow groundwater and rainwater in the Ordos Desert Plateau, North China, and river water from the nearby Yellow River, are investigated to determine the dissolved Sr source and water–rock interactions, and quantify the relative Sr contribution from each end-member. Three groundwater systems have been identified, namely, GWS-1, GWS-2 and GWS-3 according to the watershed distribution in the Ordos Desert Plateau. Ca²⁺ and Mg²⁺ are the most dominant cations in GWS-1, while Na⁺ is dominant in GWS-3. In addition, there is more SO₄²⁻ and less Cl⁻ in GWS-1 than in GWS-3. The shallow groundwater in GWS-2 seems to be geochemically between that in GWS-1 and GWS-3. The ⁸⁷Sr/⁸⁶Sr ratios of the shallow groundwater are high in GWS-1 and GWS-2 and are low in GWS-3. By geochemically comparing the nearby Yellow River, local precipitation and deep groundwater, the shallow groundwater is recharged only by local precipitation. The ionic and isotopic ratios indicate that carbonate dissolution is an important process controlling the chemistry of the shallow groundwater. The intensity of the water–rock interactions varies among the three groundwater systems and even within each groundwater system. Three end-members controlling the groundwater chemistry are isotopically identified: (1) precipitation infiltration, (2) carbonate dissolution and (3) silicate weathering. The relative Sr contributions of the three end-members show that precipitation infiltration and carbonate dissolution are the primary sources of the shallow groundwater Sr in GWS-3 whereas only carbonate dissolution is responsible for the shallow groundwater Sr in GWS-1 and GWS-2. Silicate weathering seems insignificant towards the shallow groundwater's chemistry in the Ordos Desert Plateau. This study is helpful for understanding groundwater chemistry and managing water resources.     

Isotopes (δD and δ18O) in precipitation, groundwater and surface water in the Ordos Plateau, China: implications with respect to groundwater recharge and circulation

The characteristics of δD and δ¹⁸O in precipitation, groundwater and surface water have been used to understand the groundwater flow system in the Ordos Plateau, north-central China. The slope of the local meteoric water line (LMWL) is smaller than that of the global meteoric water line (GMWL), which signifies secondary evaporation during rainfall. The distribution of stable isotopes of precipitation is influenced by temperature and the amount of precipitation. The lake water is enriched isotopically due to evaporation and its isotopic composition is closely related to the source of recharge and location in the groundwater flow systems. River water is enriched isotopically, indicating that it suffers evaporation. The deep groundwater (more than 150?m) is depleted in heavy isotopes relative to the shallow groundwater (less than 150?m), suggesting that deep groundwater may have been recharged during the late Pleistocene and early Holocene, when the climate was wetter and colder than at present. All groundwater samples plot around the LMWL, implying groundwater is of meteoric origin. Shallow groundwater has undergone evaporation and the average evaporation loss is 53%. There are two recharge mechanisms: preferential flow, and the mixture of evaporated soil moisture and subsequent rain.     

Hydrochemical evolution of groundwater in a riparian zone affected by acid mine drainage (AMD), South China: the role of river–groundwater interactions and groundwater residence time

Investigations were undertaken in the riparian zone near Shangba village, an AMD area, in southern China to determine the effects that river–groundwater interactions and groundwater residence time have had on environmental quality and geochemical evolution of groundwater. Based on the Darcy’s law and ionic mass balances as well as the method of isotopic tracer, the results showed that there were active interactions between AMD-contaminated river water and groundwater in the riparian zone in the study area. River water was found to be the main source of groundwater recharge in the northwestern part of the study area, whereas groundwater was found to be discharging into the river in the southeastern part of the study area throughout the year. End-member mixing analysis quantified that the contributions of river water to groundwater decreased gradually from 35.9% to negligible levels along the flow path. The calculated mixing concentrations of major ions indicated that water–rock reactions were the most important influence on groundwater quality. The wide range of Ca²⁺?+?Mg²⁺ and HCO₃^? ratios and the change of groundwater type from Ca²⁺–SO₄^2? type to a chemical composition dominated by Ca²⁺–HCO₃^? type indicated a change of the major water–rock reaction process from the influence of H₂SO₄ (AMD) to that of CO₂ (soil respiration) along a groundwater flow path. Furthermore, the kinetics interpretation of SO₄^2? and HCO₃^? concentrations suggested that the overlapping time of their kinetics triggered the hydrochemical evolution and the change of major weathering agent. This process might take approximately 8 years and this kinetic time will be continued when a steady source of contamination enter the aquifer.     

Geology,geochemistry, and geochronology of the Wangjiazhuang porphyry–breccia Cu(–Mo) deposit in the Zouping volcanic basin,eastern North China Block

The Wangjiazhuang porphyry–breccia Cu(–Mo) deposit is located in the Zouping volcanic basin, western Shandong Province. Seven molybdenite samples yield a Re–Os weighted mean age of 127.8 ± 0.7 Ma (2σ), which is identical within error to the zircon weighted mean ²⁰⁶Pb/²³⁸U age of 128.3 ± 1.3 Ma (2σ) determined for quartz monzonite samples. The host rock is characterized by high concentrations of K₂O (4.26–4.53 wt.%), Na₂O (4.97–5.76 wt.%), LILEs and LREEs, and high Mg# (> 40), and low concentrations of HFSEs and HREEs, with K₂O/Na₂O ratios of 0.76–0.88. The quartz monzonite also has high Sr/Y (69.9–112.5) and (La/Yb)_N (22.0–30.0) ratios, similar to adakitic rocks worldwide. Relatively low initial ⁸⁷Sr/⁸⁶Sr ratios (0.70549–0.70556), high ε_Nd(t) values (2.58–3.06), high radiogenic Pb [(²⁰⁶Pb/²⁰⁴Pb)i = 18.3424–18.4606, (²⁰⁷Pb/²⁰⁴Pb)_i = 15.5692–15.5985, (²⁰⁸Pb/²⁰⁴Pb)_i = 38.1714–38.2734] and high zircon ε_Hf(t) values (− 2.1 to + 4.3) indicate that the magma was likely derived from the partial melting of subducted oceanic crust which then reacted with the peridotitic mantle wedge. Both the breccia and porphyry ores have a narrow range of δ³⁴S (− 4.8 to + 2.1‰) and Pb isotopic compositions (²⁰⁶Pb/²⁰⁴Pb = 18.295–18.402, ²⁰⁷Pb/²⁰⁴Pb = 15.551–15.573, and ²⁰⁸Pb/²⁰⁴Pb = 38.215–38.331), suggesting that the ore metals were extracted primarily from the quartz monzonite or similar source. Subduction of the Paleo-Pacific slab during the Early Cretaceous resulted in the formation of the Wangjiazhuang quartz monzonite and associated Cu(–Mo) deposit in western Shandong Province.