Geothermal Characterization of the St. Lawrence Lowlands Sedimentary Basin,Québec,Canada

Defining temperature at depth to identify geothermal resources relies on the evaluation of the Earth heat flow based on equilibrium temperature measurements as well as thermal conductivity and heat generation rate assessment. Such high-quality geothermal data can be sparse over the region of interest. This is the case of the St. Lawrence Lowlands sedimentary basin covering 20,000 km² to the south of Québec, Canada, and enclosing only three wells up to a depth of 500 m with equilibrium heat flow measurements. However, more than 250 oil and gas exploration wells have been drilled in this area, providing for this study (parce que c'est 93 sinon) 81 locations with bottom-hole temperature up to a depth of 4300 m, however, not at equilibrium. Analyzing these data with respect to the deep geothermal resource potential of this sedimentary basin requires evaluating the thermal conductivity and heat generation rate of its geological units to properly extrapolate temperature downward. This was done by compiling literature and recent thermal conductivity measurements in outcrop and core samples as well as new heat generation rate estimates from spectral gamma ray logs to establish a first thermal assessment of geological units deep down into the basin. The mean thermal conductivity of the thermal units varies from 2.5 to 6.3 W/m·K, with peak values in the basal sandstones, while the heat generation rate varies from 1.6 to 0.3 µW/m³, decreasing from the upper caprocks toward the base of the sequence. After correcting the bottom-hole temperatures for drilling disturbance with the Harrison correction and subsequently for paleoclimate variations, results indicate a mean geothermal gradient of 23.1 °C/km, varying from 14 to 40 °C/km. Evaluating the basin thermal state from oil and gas data is a significant challenge facilitated by an understanding of its thermal properties.     

Numerical Model of Two-Phase Flow in Dissolvable Porous Media and Simulation of Reservoir Acidizing

Mathematical model of porous media dissolution coupled with two-phase flow is proposed. The model is based on the conception of dissolvable porous medium with deformable mass-variable porous skeleton. Model can be used for simulation of coupled chemo- and hydrogeomechanical processes which are difficult to examine experimentally. Acidizing of calcite oil reservoir is used as an example of the process. Water solution of hydrochloric acid and oil are two fluid phases of the model with several components. Dissolvable porous media is treated as deformable mass-variable solid phase. Change in mass of the solid phase is caused by hydrochloric acid dissolving the calcite part of the solid phase. Dissolution is supposed to be congruent; kinetics is governed by the Nernst law. Software for numerical solution of the model is developed. It uses AmgCL parallel library for high-performance computing in order to deal with large algebraic systems on the each time step of calculations. The library uses algebraic multigrid methods for preconditioning and parallel iterative solvers. NVidia CUDA framework is used as a backend to perform GPGPU calculations, because it proved to be faster than OpenCL framework on this problem. Numerical experiments on the basis of data set from real reservoirs are conducted with the developed software. Good correlation between field and calculated data is achieved. Numerical experiments for different configurations of heterogeneous layer are performed. Acidizing of layers with highly permeable conduit and with random distribution of permeability is modeled.     

Potential promoted productivity and spatial patterns of medium- and low-yield cropland land in China

With a continuously increasing population and better food consumption levels, improving the efficiency of arable land use and increasing its productivity have become fundamental strategies to meet the growing food security needs in China. A spatial distribution map of medium- and low-yield cropland is necessary to implement plans for cropland improvement. In this study, we developed a new method to identify high-, medium-, and low-yield cropland from Moderate Resolution Imaging Spectroradiometer (MODIS) data at a spatial resolution of 500 m. The method could be used to reflect the regional heterogeneity of cropland productivity because the classification standard was based on the regionalization of cropping systems in China. The results showed that the proportion of high-, medium-, and low-yield cropland in China was 21%, 39%, and 40%, respectively. About 75% of the low-yield cropland was located in hilly and mountainous areas, and about 53% of the high-yield cropland was located in plain areas. The five provinces with the largest area of high-yield cropland were all located in the Huang-Huai-Hai region, and the area amounted to 42% of the national high-yield cropland area. Meanwhile, the proportion of high-yield cropland was lower than 15% in Heilongjiang, Sichuan, and Inner Mongolia, which had the largest area allocated to cropland in China. If all the medium-yield cropland could be improved to the productive level of high-yield cropland and the low-yield cropland could be improved to the level of medium-yield cropland, the total productivity of the land would increase 19% and 24%, respectively.     

Data-Driven Predictive Modeling of Mineral Prospectivity Using Random Forests: A Case Study in Catanduanes Island (Philippines)

The Random Forests (RF) algorithm is a machine learning method that has recently been demonstrated as a viable technique for data-driven predictive modeling of mineral prospectivity, and thus, it is instructive to further examine its usefulness in this particular field. A case study was carried out using data from Catanduanes Island (Philippines) to investigate further (a) if RF modeling can be used for data-driven modeling of mineral prospectivity in areas with few (i.e., <20) mineral occurrences and (b) if RF modeling can handle predictor variables with missing values. We found that RF modeling outperforms evidential belief (EB) modeling of prospectivity for hydrothermal Au–Cu deposits in Catanduanes Island, where 17 hydrothermal Au–Cu prospects are known to exist. Moreover, just like EB modeling, RF modeling allows analysis of the spatial relationships between known prospects and individual layers of predictor data. Furthermore, RF modeling can handle missing values in predictor data through an RF-based imputation technique whereas in EB modeling, missing values are simply represented by maximum uncertainty. Therefore, the RF algorithm is a potentially useful method for data-driven predictive modeling of mineral prospectivity in regions with few (i.e., <20) occurrences of mineral deposits of the type sought. However, further testing of the method in other regions with few mineral occurrences is warranted to fully determine its usefulness in data-driven predictive modeling of mineral prospectivity.     

Environmental Geochemistry and Acid Mine Drainage Evaluation of an Abandoned Coal Waste Pile at the Alborz-Sharghi Coal Washing Plant,NE Iran

In this paper, an abandoned waste coal pile, which is resulted from Alborz-Sharghi coal washing plant, NE of Iran was mineralogically and geochemically characterized to evaluate pyrite oxidation, acid mine drainage (AMD) generation, and trace element mobility. After digging ten trenches and vertical sampling, a quantitative method including the atomic absorption test, and the quality-based methods including optical study were carried out for determination of pyrite fractions in the waste pile. The geochemical results revealed that the fraction of remaining pyrite increased with depth, indicating that pyrite oxidation is limited to the shallower depths of the pile which were confirmed by variations of sulfate, pH, EC, and carbonate with depth of the pile. To evaluate the trend of trace elements and mineralogical constituents of the waste particles, the samples were analyzed by using XRD, ICP-MS, and ICP-OES methods. The results showed the secondary and neutralizing minerals comprising gypsum have been formed below the oxidation zone. Besides, positive values of net neutralization potential indicated that AMD generation has not taken in the waste pile. In addition, variations of trace elements with depth reveal that Pb and Zn exhibited increasing trends from pile surface toward the bottom sampling trenches while another of them such as Cu and Ni had decreasing trends with increasing depth of the waste pile.     

Bitumen Prices and Structural Changes in North American Crude Oil Markets

In an earlier report, changes in bitumen prices at Hardesty, Alberta, Canada, were modeled as the responses to changes in monthly prices of Hardesty light/medium crude oil for the period 2000–2006 with a simple error correction econometric model. This note re-examines that price relationship for the period 2009–2014. Over the period 2006–2014, there was also rapid growth in North American light oil production from low-permeability carbonate, sandstone, and shale reservoirs. During that period, Canadian raw bitumen production grew by more than 12% per year and there was significant geographical diversification in its markets. Results of the statistical analysis showed that the change in the dynamic relationships between bitumen prices and Hardesty light oil prices probably reflected, in part, the maturation of bitumen markets and closer integration with North American light oil markets. The analysis also examines the dynamic relationships between bitumen prices and West Texas Intermediate and Brent international benchmark crude oil prices. Ideally, if bitumen prices are found to be closely related to a widely traded benchmark crude oil, the benchmark crude oil price forecasts could be used as a basis for predicting bitumen prices. However, neither of international benchmark crude oils tested had high explanatory power.     

Capacity of soil loss control in the Loess Plateau based on soil erosion control degree

The capacity of soil and water conservation measures, defined as the maximum quantity of suitable soil and water conservation measures contained in a region, were determined for the Loess Plateau based on zones suitable for establishing terraced fields, forestland and grassland with the support of geographic information system (GIS) software. The minimum possible soil erosion modulus and actual soil erosion modulus in 2010 were calculated using the revised universal soil loss equation (RUSLE), and the ratio of the minimum possible soil erosion modulus under the capacity of soil and water conservation measures to the actual soil erosion modulus was defined as the soil erosion control degree. The control potential of soil erosion and water loss in the Loess Plateau was studied using this concept. Results showed that the actual soil erosion modulus was 3355 t?km^-2?a^-1, the minimum possible soil erosion modulus was 1921 t?km^-2?a^-1, and the soil erosion control degree was 0.57 (medium level) in the Loess Plateau in 2010. In terms of zoning, the control degree was relatively high in the river valley-plain area, soil-rocky mountainous area, and windy-sandy area, but relatively low in the soil-rocky hilly-forested area, hilly-gully area and plateau-gully area. The rate of erosion areas with a soil erosion modulus of less than 1000 t?km^-2?a^-1 increased from 50.48% to 57.71%, forest and grass coverage rose from 56.74% to 69.15%, rate of terraced fields increased from 4.36% to 19.03%, and per capita grain available rose from 418 kg?a^-1 to 459 kg?a^-1 under the capacity of soil and water conservation measures compared with actual conditions. These research results are of some guiding significance for soil and water loss control in the Loess Plateau.     

Degrading river network due to urbanization in Yangtze River Delta

Evolution of river systems under the background of human activities has been a heated topic among geographers and hydrologists. Spatial and temporal variations of river systems during the 1960s–2010s in the Yangtze River Delta (YRD) were investigated based on streams derived from the topographic maps in the 1960s, 1980s and 2010s. A list of indices, drainage density (Dd), water surface ratio (WSR), ratio of area to length of main streams (R), evolution coefficient of tributaries (K) and box dimension (D), were classified into three types (quantitative, structural, and complex indices) and used to quantify the variations of stream structure. Results showed that: (1) quantitative indices (Dd, WSR) presented decreasing trend in the past 50 years, and Dd in Wuchengxiyu, Hangjiahu and Yindongnan have decreased most, about 20%. Structurally, the Qinhuai River basin was characterized by significant upward R, and K value in Hangjiahu went down dramatically by 46.8% during the 1960s–2010s. Decreasing tendency in D was found dominating across the YRD, and decreasing magnitude in Wuchengxiyu and Hangjiahu peaks for 7.8% and 6.5%, respectively in the YRD. (2) Urbanization affected the spatial pattern of river system, and areas with high level of urbanization exhibited least Dd (2.18 km/km²), WSR (6.52%), K (2.64) and D (1.42), compared to moderate and low levels of urbanization. (3) Urbanization also affected the evolution of stream system. In the past 50 years, areas with high level of urbanization showed compelling decreasing tendency in quantitative (27.2% and 19.3%) and complex indices (4.9%) and trend of enlarging of main rivers (4.5% and 7.9% in periods of the 1960s–1980s and the 1980s–2010s). In the recent 30 years, areas with low level of urbanization were detected with significant downward trend in Dd and K. (4) Expanding of urban land, construction of hydraulic engineering and irrigation and water conservancy activities were the main means which degraded the river system in the YRD.     

Construction and progress of Chinese terrestrial ecosystem carbon,nitrogen and water fluxes coordinated observation

Eddy Covariance technique (EC) achieves the direct measurement on ecosystem carbon, nitrogen and water fluxes, and it provides scientific data for accurately assessing ecosystem functions in mitigating global climate change. This paper briefly reviewed the construction and development of Chinese terrestrial ecosystem flux observation and research network (ChinaFLUX), and systematically introduced the design principle and technology of the terrestrial ecosystem carbon, nitrogen and water fluxes coordinated observation system of ChinaFLUX. In addition, this paper summarized the main progress of ChinaFLUX in the ecosystem carbon, nitrogen and water exchange and environmental controlling mechanisms, the spatial pattern of carbon, nitrogen and water fluxes and biogeographical mechanisms, and the regional terrestrial ecosystem carbon budget assessment. Finally, the prospects and emphases of the terrestrial ecosystem carbon, nitrogen and water fluxes coordinated observation of ChinaFLUX are put forward to provide theoretical references for the development of flux observation and research in China.     

Spatiotemporal distribution and historical evolution of polders in the Dongting Lake area,China

Dongting Lake is the largest lake in the middle reaches of the Yangtze River in China. For centuries, people inhabiting the Dongting Lake area have been reclaiming land and constructing dams for flood resistance, agricultural production, and rural settlement, forming geographical entities known as polders. In this study, the regional spatial distributions of polders in the Dongting Lake area in 1949, 1998, and 2013 were obtained using historical maps and modern remotely sensed data, revealing changes since the establishment of the People’s Republic of China. Nanxian County was then selected to demonstrate polder changes at the county level, because it has undergone the most dramatic changes in the area. Different polder change models for the Datonghu, Yule, and Renhe polders were analyzed for eight periods: 1644 (the early Qing Dynasty), 1911 (the late Qing Dynasty), 1930 (the Republic of China), 1949 (the People’s Republic of China), 1963, 1970, 1998, and 2013. Three resulting polder evolution models are: 1) reclaiming polders from lakes, 2) integrating polders by stream merging, and 3) abandoning polders for flood release. The polder evolution models demonstrate the wisdom of local people in using land resources according to the specific regional conditions. Throughout their long-term historical evolution, the spatial distribution of polders in the Dongting Lake area tended to be homogeneous, and the degree of human disturbance tended to be stable. However, a shift occurred, from pure polder area growth or removal to more comprehensive management and protection of the regional environment.