Urban land cover has major impacts on a city’s ecosystem services and the inherent quality of its urban residential environment. The spatio-temporal distribution of impervious surface area and green areas in Chinese cities has exhibited a significantly marked difference in comparison with USA cities. This study focused on monitoring and comparing the spatio-temporal dynamics, land cover patterns and characteristics of functional regions in six Chinese (n=3) and USA (n=3) cities. The study data were collated from Landsat TM/MSS imagery during the period 1978-2010. Results indicate that Chinese cities have developed compactly over the past three decades, while development has been notably dispersed among USA cities. Mean vegetation coverage in USA cities is approximately 2.2 times that found amongst Chinese urban agglomerations. Land use types within Chinese cities are significantly more complex, with a higher density of impervious surface area. Conversely, the central business district (CBD) and residential areas within USA cities were comprised of a lower proportion of impervious surface area and a higher proportion of green land. Results may be used to contribute to future urban planning and administration efforts in both China and the USA.

The quantity and spatial pattern of farmland has changed in China, which has led to a major change in the production potential under the influence of the national project of ecological environmental protection and rapid economic growth during 1990-2010. In this study, the production potential in China was calculated based on meteorological, terrain elevation, soil and land-use data from 1990, 2000 and 2010 using the Global Agro-ecological Zones model. Then, changes in the production potential in response to farmland changes from 1990 to 2010 were subsequently analyzed. The main conclusions were the following. First, the total production potential was 1.055 billion tons in China in 2010. Moreover, the average production potential was 7614 kg/ha and showed tremendous heterogeneity in spatial pattern. Total production in eastern China was high, whereas that in northwestern China was low. The regions with high per unit production potential were mainly distributed over southern China and the middle and lower reaches of the Yangtze River. Second, the obvious spatiotemporal heterogeneity in farmland changes from 1990 to 2010 had a significant influence on the production potential in China. The total production potential decreased in southern China and increased in northern China. Furthermore, the center of growth of the production potential moved gradually from northeastern China to northwestern China. The net decrease in the production potential was 2.97 million tons, which occupied 0.29% of the national total actual production in 2010. Third, obvious differences in the production potential in response to farmland changes from 1990 to 2000 and from 2000 to 2010 were detected. The net increase in the production potential during the first decade was 10.11 million tons and mainly distributed in the Northeast China Plain and the arid and semi-arid regions of northern China. The net decrease in the production potential during the next decade was 13.08 million tons and primarily distributed in the middle and lower reaches of the Yangtze River region and the Huang-Huai-Hai Plain. In general, the reason for the increase in the production potential during the past two decades might be due to the reclamation of grasslands, woodlands and unused land, and the reason for the decrease in the production potential might be urbanization that occupied the farmland and Green for Grain Project, which returned farmland to forests and grasslands.

Surface albedo is a primary causative variable associated with the process of surface energy exchange. Numerous studies have examined diurnal variation of surface albedo at a regional scale; however, few studies have analyzed the intra-annual variations of surface albedo in concurrence with different land cover types. In this study, we amalgamated surface albedo product data (MCD43) from 2001 to 2008, land-use data (in 2000 and 2008) and land cover data (in 2000); quantitative analyses of surface albedo variation pertaining to diverse land cover types and the effect of the presence/absence of ground snow were undertaken. Results indicate that intra-annual surface albedo values exhibit flat Gaussian or triangular distributions depending upon land cover types. During snow-free periods, satellite observed surface albedo associated with the non-growing season was lower than that associated with the growing season. Satellite observed surface albedo during the presence of ground snow period was 2-4 times higher than that observed during snow-free periods. Surface albedo reference values in typical land cover types have been calculated; notably, grassland, cropland and built-up land were associated with higher surface albedo reference values than barren while ground snow was present. Irrespective of land cover types, the lowest surface albedo reference values were associated with forested areas. Proposed reference values may prove extremely useful in diverse research areas, including ecological modeling, land surface process modeling and radiation energy balance applications.

The aim of this study is to establish several important factors representing land use intensification in cultivated land (denoted by CII), using a multi-dimensional approach to achieve realistic and practical cultivated land use policies in China. For this reason, the theoretical framework was first built to explain the changes of land use intensification in the cultivated land, and then the variables and index were further developed for the purpose of characterizing the dynamic trends and driving forces of the land use intensification in the cultivated land at the provincial level. The study results indicate that the extent of CII significantly increased during the period of 1996 to 2008, due to the extensive use of fertilizers, machinery and pesticide, increased labor and capital input, and intensified land use. Moreover, the principal component regression results show that the productivity of cultivated land, economic benefits of cultivated land, labor productivity, and land use conversion are the main factors affecting the village development. The first three factors play a positive role, while the last one has a negative effect on the land use intensification in the cultivated land. According to these results, the main policies for sustainable intensification in cultivated land are proposed. First, the sustainable pathways for intensification should be adopted to reduce the unsustainable uses of chemical fertilizer, agricultural chemicals, etc. Second, the conditions for agricultural production should be further improved to increase the cultivated land productivity. Third, it is very necessary and helpful for improving labor productivity and land use efficiency from the viewpoint of accelerated the cultivated land circulation. The last step is to positively affect the production activities of peasants by means of reforming the subsidy standards.

Based on the daily precipitation from a 0.5°×0.5° gridded dataset and meteorological stations during 1961-2011 released by National Meteorological Information Center, the reliability of this gridded precipitation dataset in South China was evaluated. Five precipitation indices recommended by the World Meteorological Organization (WMO) were selected to investigate the changes in precipitation extremes of South China. The results indicated that the bias between gridded data interpolated to given stations and the corresponding observed data is limited, and the proportion of the number of stations with bias between -10% and 0 is 50.64%. The correlation coefficients between gridded data and observed data are generally above 0.80 in most parts. The average of precipitation indices shows a significant spatial difference with drier northwest section and wetter southeast section. The trend magnitudes of the maximum 5-day precipitation (RX5day), very wet day precipitation (R95), very heavy precipitation days (R20mm) and simple daily intensity index (SDII) are 0.17 mm·a^-1, 1.14 mm·a^-1, 0.02 d·a^-1 and 0.01 mm·d^-1·a^-1, respectively, while consecutive wet days (CWD) decrease by -0.05 d·a^-1 during 1961-2011. There is spatial disparity in trend magnitudes of precipitation indices, and approximate 60.85%, 75.32% and 75.74% of the grid boxes show increasing trends for RX5day, SDII and R95, respectively. There are high correlations between precipitation indices and total precipitation, which is statistically significant at the 0.01 level.

Data of flood, drought, hailstorms, and low temperature events in Xinjiang from 1949 to 2012 were analyzed with the diffusion method to assess the risk of the most common types of disasters in Xinjiang. It was proved that the frequency and intensity of meteorological disasters of the study area showed an increasing trend associated with global warming. Among the four types of disasters, surpass probability of drought was the largest, followed by hailstorm, low temperature and flood in turn. Moreover, the wavelet method analysis revealed that greater oscillations had occurred since 2000, which may be associated with the occurrence of extreme climatic changes. The spatial distribution of frequencies reveals that the northern slope of Tianshan Mountains is a multiple disaster area, the southern slope of Tianshan is the area where more floods and hailstorms occur, and the west of Turpan-Hami Basin is the area wind is prevalent. The relationships between disaster-affected areas and corresponding meteorological and socio-economic indexes were also analyzed. It indicated that there were significant positive correlations between the areas affected and the most meteorological and socio-economic indicators except the grain acreage.

This paper examines the changes in the time series of water discharge and sediment load of the Yellow River into the Bohai Sea. To determine the characteristics of abrupt changes and multi-scale periods of water discharge and sediment load, data from Lijin station were analyzed, and the resonance periods were then calculated. The Mann-Kendall test, order clustering, power-spectrum, and wavelet analysis were used to observe water discharge and sediment load into the sea over the last 62 years. The most significant abrupt change in water discharge into the sea occurred in 1985, and an abrupt change in sediment load happened in the same year. Significant decreases of 64.6% and 73.8% were observed in water discharge and sediment load, respectively, before 1985. More significant abrupt changes in water discharge and sediment load were observed in 1968 and 1996. The characteristics of water discharge and sediment load into the Bohai Sea show periodic oscillation at inter-annual and decadal scales. The main periods of water discharge are 9.14 years and 3.05 years, whereas the main periods of sediment load are 10.67 years, 4.27 years, and 2.78 years. The significant resonance periods between water discharge and sediment load are observed at the following temporal scales: 2.86 years, 4.44 years, and 13.33 years. Water discharge and sediment load started to decrease after 1970 and has decreased significantly since 1985 for several reasons. Firstly, the precipitation of the Yellow River drainage area has reduced since 1970. Secondly, large-scale human activities, such as the building of reservoirs and floodgates, have increased. Thirdly, water and soil conservation have taken effect since 1985.

The distribution of plant species and relationships between species and soil factors in the east central part of Gurbantunggut Desert was studied to provide more insight into the flora and determine differences in vegetation across various parts of the desert. Two-way Cluster Analysis showed that the vegetation in the area could be divided into three groups, the first group was dominated by the shrub species, Ephedra przewalskii and the grass species, Carex physodes mainly in areas of flat grounds and gentle slopes; the second group was dominated by C. physodes, Artemisia songorica and A. xerophytica mainly on the slope of sand dunes and the third group was dominated by the shrub species, Haloxylon persicum mainly on the top of sand dunes. There was no difference in plant density between Groups 1 and 2 but there was a significant decrease in Group 3. Soil water under vegetation Group 3 was much lower than that in the other two groups at all soil depths. The EC, organic matter, total P and soluble Na, Ca and Mg varied very similarly with soil water. Canonical correspondence analysis (CCA) satisfactorily assessed the species-soil relations in the area. The distribution of plant species was strongly correlated with the soil factors of water content, organic matter, EC and nutrients. The variations in species occurrence explained by the three CCA axes were about 70%, indicating that some explanatory site variables may exist outside our studied parameters. Soil texture is suggested to be included in future studies to improve the explanation of CCA.

The Forest Landscape Model (FLM) is an efficiency tool of quantified expression of forest ecosystem’s structure and function. This paper, on the basis of identifying FLM, according to the stage of development, summarizes the development characteristics of the model, which includes the theoretical foundation of mathematical model, FLM of stand-scale, primary development of spatial landscape model, rapid development of ecosystem process model as the priority, and developing period of structure and process driven by multi-factor. According to the characteristics of different FLMs, this paper classifies the existing FLM in terms of mechanism, property and application, and elaborates the identifications, advantages and disadvantages of different types of models. It summarizes and evaluates the main application fields of existing models from two aspects which are the changes of spatial pattern and ecological process. Eventually, this paper presents FLM’s challenges and directions of development in the future, including: (1) more prominent service on the practical strategy of forest management’s objectives; (2) construction of multi-modules and multi-plugin to satisfy landscape research demand in various conditions; (3) adoption of high resolution’s spatial-temporal data; (4) structural construction of multi-version module; (5) improving the spatial suitability of model application.

In this study, we have used four methods to investigate the start of the growing season (SGS) on the Tibetan Plateau (TP) from 1982 to 2012, using Normalized Difference Vegetation Index (NDVI) data obtained from Global Inventory Modeling and Mapping Studies (GIMSS, 1982-2006) and SPOT VEGETATION (SPOT-VGT, 1999-2012). SGS values estimated using the four methods show similar spatial patterns along latitudinal or altitudinal gradients, but with significant variations in the SGS dates. The largest discrepancies are mainly found in the regions with the highest or the lowest vegetation coverage. Between 1982 and 1998, the SGS values derived from the four methods all display an advancing trend, however, according to the more recent SPOT VGT data (1999-2012), there is no continuously advancing trend of SGS on the TP. Analysis of the correlation between the SGS values derived from GIMMS and SPOT between 1999 and 2006 demonstrates consistency in the tendency with regard both to the data sources and to the four analysis methods used. Compared with other methods, the greatest consistency between the in situ data and the SGS values retrieved is obtained with Method 3 (Threshold of NDVI ratio). To avoid error, in a vast region with diverse vegetation types and physical environments, it is critical to know the seasonal change characteristics of the different vegetation types, particularly in areas with sparse grassland or evergreen forest.

Hydrological processes were compared, with and without the influence of precipitation on discharge, to identify the differences between glacierized and non-glacierized catchments in the Urumqi River source region, on the northern slope of the eastern Tianshan Mountains, during the melting season (May-September) in 2011. The study was based on hydrological data observed at 10-min intervals, meteorological data observed at 15-min intervals, and glacier melting and snow observations from the Empty Cirque, Zongkong, and Urumqi Glacier No.1 gauging stations. The results indicated that the discharge differed markedly among the three gauging stations. The daily discharge was more than the nightly discharge at the Glacier No.1 gauging station, which contrasted with the patterns observed at the Zongkong and Empty Cirque gauging stations. There was a clear daily variation in the discharge at the three gauging stations, with differences in the magnitude and duration of the peak discharge. When precipitation was not considered, the time-lags between the maximum discharge and the highest temperature were 1-3 h, 10-16 h, and 5-11 h at the Glacier No.1, Empty Cirque, and Zongkong gauging stations, respectively. When precipitation was taken into consideration, the corresponding time-lags were 0-1 h, 13 h, and 6-7 h, respectively. Therefore, the duration from the generation of discharge to confluence was the shortest in the glacierized catchment and the longest in the catchment where was mainly covered by snow. It was also shown that the hydrological process from the generation of discharge to confluence shortened when precipitation was considered. The factors influencing changes in the discharge among the three gauging stations were different. For Glacier No.1 station, the discharge was mainly controlled by heat conditions in the glacierized region, and the discharge displayed an accelerated growth when the temperature exceeded 5°C in the melt season. It was found that the englacial and subglacial drainage channel of Glacier No.1 had become simpler during the past 20 years. Its weaker retardance and storage of glacier melting water resulted in rapid discharge confluence. It was also shown that the discharge curve and the time-lag between the maximum discharge and the highest temperature could be used to reveal the evolution of the drainage system and the process of glacier and snow melting at different levels of glacier coverage.

The net accumulation record of ice core is one of the most reliable indicators for reconstructing precipitation changes in high mountains. A 20.12 m ice core was drilled in 2006 from the accumulation zone of Laohugou Glacier No.12 in the northeastern Tibetan Plateau, China. We obtained the precipitation from the ice core net accumulation during 1960-2006, and found out the relationship between Laohugou ice core record and other data from surrounding sites of the northeastern Tibetan Plateau. Results showed that during 1960-2006, the precipitation in the high mountains showed firstly an increasing trend, while during 1980 to 2006 it showed an obvious decreasing trend. Reconstructed precipitation change in the Laohugou glacier basin was consistent with the measured data from the nearby weather stations in the lower mountain of Subei, and the correlation coefficient was 0.619 (P<0.001). However, the precipitation in the high mountain was about 3 times more than that of the lower mountain. The precipitation in Laohugou Glacier No.12 of the western Qilian Mountains corresponded well to the net accumulation of Dunde ice core during the same period, tree-ring reconstructed precipitation, the measured data of multiple meteorological stations in the northeastern Tibetan Plateau, and also the changes of adjacent PDSI drought index. Precipitation changes of the Laohugou glacier basin and other sites of the northeastern Tibetan Plateau had significantly positive correlation with ENSO, which implied that the regional alpine precipitation change was very likely to be influenced by ENSO.

The objective of this study is to quantitatively evaluate Tropical Rainfall Measuring Mission (TRMM) data with rain gauge data and further to use this TRMM data to drive a Distributed Time-Variant Gain Model (DTVGM) to perform hydrological simulations in the semi-humid Weihe River catchment in China. Before the simulations, a comparison with a 10-year (2001-2010) daily rain gauge data set reveals that, at daily time step, TRMM rainfall data are better at capturing rain occurrence and mean values than rainfall extremes. On a monthly time scale, good linear relationships between TRMM and rain gauge rainfall data are found, with determination coefficients R² varying between 0.78 and 0.89 for the individual stations. Subsequent simulation results of seven years (2001-2007) of data on daily hydrological processes confirm that the DTVGM when calibrated by rain gauge data performs better than when calibrated by TRMM data, but the performance of the simulation driven by TRMM data is better than that driven by gauge data on a monthly time scale. The results thus suggest that TRMM rainfall data are more suitable for monthly streamflow simulation in the study area, and that, when the effects of recalibration and the results for water balance components are also taken into account, the TRMM 3B42-V7 product has the potential to perform well in similar basins.

Understanding the topographic context preceding the development of erosive landforms is of major relevance in geomorphic research, as topography is an important factor on both water and mass movement-related erosion, and knowledge of the original surface is a condition for quantifying the volume of eroded material. Although any reconstruction implies assuming that the resulting surface reflects the original topography, past works have been dominated by linear interpolation methods, incapable of generating curved surfaces in areas with no data or values outside the range of variation of inputs. In spite of these limitations, impossibility of validation has led to the assumption of surface representativity never being challenged. In this paper, a validation-based method is applied in order to define the optimal interpolation technique for reconstructing pre-erosion topography in a given study area. In spite of the absence of the original surface, different techniques can be nonetheless evaluated by quantifying their capacity to reproduce known topography in unincised locations within the same geomorphic contexts of existing erosive landforms. A linear method (Triangulated Irregular Network, TIN) and 23 parameterizations of three distinct Spline interpolation techniques were compared using 50 test areas in a context of research on large gully dynamics in the South of Portugal. Results show that almost all Spline methods produced smaller errors than the TIN, and that the latter produced a mean absolute error 61.4% higher than the best Spline method, clearly establishing both the better adjustment of Splines to the geomorphic context considered and the limitations of linear approaches. The proposed method can easily be applied to different interpolation techniques and topographic contexts, enabling better calculations of eroded volumes and denudation rates as well as the investigation of controls by antecedent topographic form over erosive processes.

Urban expansion models are useful tools to understand urbanization process and have been given much attention. However, urban expansion is a complicated socio-economic phenomenon that is affected by complex and volatile factors involving in great uncertainties. Therefore, the accurate simulation of the urban expansion process remains challenging. In this paper, we make an attempt to solve such uncertainty through a reversal process and view urban expansion as a process wherein the urban landscape overcomes resistance from other landscapes. We developed an innovative approach derived from the minimum cumulative resistance (MCR) model that involved the introduction of a relative resistance factor for different source levels and the consideration of rigid constraints on urban expansion caused by ecological barriers. Using this approach, the urban expansion ecological resistance (UEER) model was created to describe ecological resistance surfaces suitable for simulating urban expansion and used to simulate urban expansion in Guangzhou. The study results demonstrate that the ecological resistance surface generated by the UEER model comprehensively reflects ecological resistance to urban expansion and indicates the spatial trends in urban expansion. The simulation results from the UEER-based model were more realistic and more accurately reflected ecological protection requirements than the conventional MCR-based model. These findings can enhance urban expansion simulation methods.

Groundwater is the most appropriate and widely used source of drinking water, which is increasingly threatened by pollution from industrial and agricultural activities. To check the severity of the problem, 156 groundwater samples were collected from various depths (60-110 ft) of 52 different localities in Faisalabad city, the third largest metropolis in Pakistan, and analyzed for the metals (Zn, Cu, Cd, Ni, Pb, Mn and Fe) concentration in 2009. Quantification was done by using Flame Atomic Absorption Spectrophotometer technique and the results were compared with WHO standards for drinking water quality. Results showed that the levels of Cu, Mn and Fe were below the WHO standards while the concentrations of Zn, Cd, Ni and Pb were above the recommended levels of safe drinking water. Correlation analysis among the occurrence of these heavy metals revealed a highly significant and positive correlation of Mn with Zn and Fe. A significant and positive correlation of Cd was also found with Cu and groundwater depth showing that there is strong association between Cu-Cd pair and that the Cd concentration varies with depth of groundwater in the study area. Regional patterns of heavy metals occurrence were mapped using Geographical Information System (GIS) for the identification and demarcation of risk areas. The concentration maps may be used by policymakers of the city to mitigate groundwater pollution.

China has witnessed unprecedented urbanization over the past decades. The rapid expansion of urban population has been dominantly contributed by the floating population from rural areas, of which the spatiotemporal patterns, driving forces, and multidimensional effects are scrutinized and evaluated in this study by using the latest national censuses conducted in 2000 and 2010. Analysis based on the county-level data comes to conclusions as follows. The spatial pattern of floating population has remained stable over the first decade of the new century. The top 1% cities with the largest floating population received 45.5% of all migrants in China. As the rapid development of mega-city regions, the coastal concentration areas of floating population tended to geographically united as a whole, whereas the spatial distribution of migrants within each region varied significantly. The migrant concentration area in the Yangtze River Delta was the largest and its expansion was also the most salient. However, the floating population has growingly moved into provincial capitals and other big cities in the inland regions and its gravity center has moved northward for around 110 km during the study period. The spatial pattern of floating population has been formed jointly by the state and market forces in transitional China and the impacts of state forces have been surpassed by those of market forces in the country as a whole. The attractiveness of coastal cities and counties to the floating population comes mainly from the nonagricultural employment opportunities and public services, reflecting that long-distance and long-term migrants have moved coastward not only to gain employment but also to enjoy city life. By contrast, in the central and western regions, places with a higher economic development level and at a higher administrative level are more attractive to floating populations, demonstrating that the state remains to play an important role in allocating economic resources and promoting regional development in inland China. As the main body of new urban residents, the floating population has contributed substantially to the elevation of the urbanization levels of migrant-sending and -receiving places, by 20.0% and 49.5% respectively. Compared with extensively investigated interprovincial migrants, intra-provincial migrants have higher intention and ability to permanently live in cities and thus might become the main force of China’s urbanization in the coming decades. The internal migration has also reshaped China’s urban system in terms of its hierarchical organization and spatial structure.

In this study, we analyzed the spatiotemporal variation of cold surges in Inner Mongolia between 1960 and 2012 and their possible driving factors using daily minimum temperature data from 121 meteorological stations in Inner Mongolia and the surrounding areas. These data were analyzed utilizing a piecewise regression model, a Sen+Mann- Kendall model, and a correlation analysis. Results demonstrated that (1) the frequency of single-station cold surges decreased in Inner Mongolia during the study period, with a linear tendency of -0.5 times/10a (-2.4 to 1.2 times/10a). Prior to 1991, a significant decreasing trend of -1.1 times/10a (-3.3 to 2.5 times/10a) was detected, while an increasing trend of 0.45 times/10a (-4.4 to 4.2 times/10a) was found after 1991. On a seasonal scale, the trend in spring cold surges was consistent with annual values, and the most obvious change in cold surges occurred during spring. Monthly cold surge frequency displayed a bimodal structure, and November witnessed the highest incidence of cold surge. (2) Spatially, the high incidence of cold surge is mainly observed in the northern and central parts of Inner Mongolia, with a higher occurrence observed in the northern than in the central part. Inter-decadal characteristic also revealed that high frequency and low frequency regions presented decreasing and increasing trends, respectively, between 1960 and 1990. High frequency regions expanded after the 1990s, and regions exhibiting high cold surge frequency were mainly distributed in Tulihe, Xiao’ergou, and Xi Ujimqin Banner. (3) On an annual scale, the cold surge was dominated by AO, NAO, CA, APVII, and CQ. However, seasonal differences in the driving forces of cold surges were detected. Winter cold surges were significantly correlated with AO, NAO, SHI, CA, TPI, APVII, CW, and IZ, indicating they were caused by multiple factors. Autumn cold surges were mainly affected by CA and IM, while spring cold surges were significantly correlated with CA and APVII.

Seasonal water-level fluctuations (WLF) play a dominate role in lacustrine ecosystems. River-lake interaction is a direct factor in changes of seasonal lake WLF, especially for those lakes naturally connected to upstream and downstream rivers. During the past decade, the modification of WLF in the Poyang Lake (the largest freshwater lake in China) has caused intensified flood and irrigation crises, reduced water availability, compromised water quality and extensive degradation of the lake ecosystem. There has been a conjecture as to whether the modification was caused by its interactions with Yangtze River. In this study, we investigated the variations of seasonal WLF in China’s Poyang Lake by comparing the water levels during the four distinct seasons (the dry season, the rising season, the flood season, and the retreating season) before and after 2003 when the Three Gorge Dam operated. The Water Surface Slope (WSS) was used as a representative parameter to measure the changes in river-lake interaction and its impacts on seasonal WLF. The results showed that the magnitude of seasonal WLF has changed considerably since 2003; the seasonal WLF of the Poyang Lake have been significantly altered by the fact that the water levels both rise and retreat earlier in the season and lowered water levels in general. The fluctuations of river-lake interactions, in particular the changes during the retreating season, are mainly responsible for these variations in magnitude of seasonal WLF. This study demonstrates that WSS is a representative parameter to denote river-lake interactions, and the results indicate that more emphasis should be placed on the decrease of the Poyang Lake caused by the lowered water levels of the Yangtze River, especially in the retreating season.

Detecting variation trend in dry-wet conditions can provide information for developing strategic measures to mitigate the impacts of global warming, particularly in dry regions. Taking the hilly region of northern Shaanxi on the Loess Plateau as a case area, this study analyzed the trend of aridity variation during 1981-2012, and explored the effect of vegetation restoration promoted by the Grain-for-Green (GFG) program implemented in 1999. The results indicated that the aridity in the region was non-significantly increased by 0.88% per year during 1981-2012, showing a drying trend. This drying trend and amplitude were changed by the influence of vegetation restoration promoted by the GFG program, based on two findings. The first one was that the aridity variation tended to increase during 1981-1999 while it turned to decrease during 2000-2012, with the regional mean relative change rate changed from 2.45% to -1.06%. This distinction was more remarkable in the loess gully region, where the vegetation was improved more obviously. The second one was that the mean vegetation coverage as indicated by EVI increased by 0.90% to 4.32% per year at county level, while the aridity decreased by 0.14% to 2.32% per year during 2000-2012. The regression analysis using the mean county data indicated that the change rate of aridity was negatively related to that of EVI with the coefficient of determination (R²) of 0.56, illustrating that around half of the aridity decline was explained by the EVI change. The mechanism of this effect was complicated, but it was found that the wind speed decline induced by the vegetation improvement could be an important contributor. It is concluded that the region became drier during 1981-2012, but the eco-restoration reduced the drying speed. However, this conclusion is involved in uncertainties, and further study based on experiments is needed to confirm the effect of the GFG-promoted vegetation restoration.

We present the first quantitative estimation of monsoon precipitation during the late glacial-Holocene in the sandy land of northern China, based on organic carbon isotopic composition data from a loess-sand sequence at margin of the Mu Us sandy land. We use the relationship between monsoon precipitation and the carbon isotopic composition of modern soils as an analogue, with a minor modification, to reconstruct precipitation back to c. 47 ka ago. The preliminary results indicate that annual monsoon precipitation was high after 8 ka, with an average of 435 mm; and it decreased during 18 and 8 ka with a mean value of 194 mm. The precipitation value of 47-18 ka varied between the two. We compare the reconstructed precipitation with other records and paleoclimatic modeling results, showing that our record agrees with reconstructions of the monsoon precipitation from other sources, even capturing short climatic events such as the Younger Dryas. We suggest that solar irradiance, high-latitude temperature/ice volume and local evaporation have together modified moistures in the sandy land.

Runoff calculation is one of the key components in the hydrological modeling. For a certain spatial scale, runoff is a very complex nonlinear process. Currently, the runoff yield model in different hydrological models is not unique. The Chinese LCM model and the American SCS model describe runoff at the macroscopic scale, taking into account the relationship between total actual retention and total rainfall and having a certain similarity. In this study, by comparing the two runoff yield models using theoretical analyses and numerical simulations, we have found that: (1) the SCS model is a simple linear representation of the LCM model, and the LCM model reflects more significantly the nonlinearity of catchment runoff. (2) There are strict mathematical relationships between parameters (R, r) of the LCM model and between parameters (S) of the SCS model, respectively. Parameters (R, r) of the LCM can be determined using the research results of the SCS model parameters. (3) LCM model parameters (R, r) can be easily obtained by field experiments, while SCS parameters (S) are difficult to measure. Therefore, parameters (R, r) of the LCM model also can provide the foundation for the SCS model. (4) The SCS model has a linear relationship between the reciprocal of total actual retention and the reciprocal of total rainfall during runoff period. The one-order terms of a Taylor series expansion of the LCM model describe the same relationship, which is worth further study.

According to the highway data and some socioeconomic data of 1990, 1994, 2000, 2005 and 2009 of county units in the Pearl River Delta, this paper measured urban integrated power of different counties in different years by factor analysis, and estimated each county’s potential in each year by means of expanded potential model. Based on that, the spatio-temporal association patterns and evolution of county potential were analyzed using spatio-temporal autocorrelation methods, and the validity of spatio-temporal association patterns was verified by comparing with spatial association patterns and cross-correlation function. The main results are shown as follows: (1) The global spatio-temporal association of county potential showed a positive effect during the study period. But this positive effect was not strong, and it had been slowly strengthened during 1994-2005 and decayed during 2005-2009. The local spatio-temporal association characteristics of most counties’ potential kept relatively stable and focused on a positive autocorrelation, however, there were obvious transformations in some counties among four types of local spatio-temporal association (i.e., HH, LL, HL and LH). (2) The distribution difference and its change of local spatio-temporal association types of county potential were obvious. Spatio-temporal HH type units were located in the central zone and Shenzhen-Dongguan region of the eastern zone, but the central spatio-temporal HH area shrunk to the Guangzhou-Foshan core metropolitan region only after 2000; the spatio-temporal LL area in the western zone kept relatively stable with a surface-shaped continuous distribution pattern, new LL type units emerged in the south-central zone since 2005, the eastern LL area expanded during 1994-2000, but then gradually shrunk and scattered at the eastern edge in 2009; the spatio-temporal HL and LH areas varied significantly. (3) The local spatio-temporal association patterns of county potential among the three zones presented significant disparity, and obvious difference between the eastern and central zones tended to decrease, whereas that between the western zone and the central and eastern zones further expanded. (4) Spatio-temporal autocorrelation methods can efficiently mine the spatio-temporal association patterns of county potential, and can better reveal the complicated spatio-temporal interaction between counties than ESDA methods.

A network perspective has increasingly become an organizational paradigm for understanding regional spatial structures. Based on a critical overview of existing empirical models for estimating intercity networks based on firm linkages, this study extends the recently proposed regional corporate city model algorithm by proposing a new method for approximating urban networks based on the locational strategies of firms. The new method considers both regional and hierarchical network features and avoids the information loss associated with the conversion from two-mode firm-city networks to one-mode city-city networks. In addition, networks estimated by using the method proposed herein are suitable when employing social network analysis. Finally, this method is empirically validated by examining intercity firm networks formed by advanced producer services firms in China’s two largest metropolitan areas, namely the Yangtze River Delta and Pearl River Delta. The presented empirical analysis suggests two main findings. First, in contrast to conventional methods (e.g., the interlocking city network model), our new method produces regional and hierarchical urban networks that more closely resemble reality. Second, the new method allows us to use social network analysis to assess betweenness and closeness centralities. However, regardless of the model applied, the validity of any method that measures urban networks depends on the soundness of its underlying assumptions about how network actors (firms, in our case) interact.

Using a sample of 14 prefecture-level cities in Liaoning Province, this study first explored the spatial hierarchy and structural characteristics of energy efficiency from the following three viewpoints: energy technical efficiency based on data envelopment analysis, energy consumption per unit of GDP, and energy utilization efficiency combining the previous two indexes. After measuring and analyzing the advancement, rationality, and concentration of the industrial structure in each city, we made some generalizations about the coupling features of the energy efficiency and industrial structure in Liaoning, using the coupling degree rating model. Some of our conclusions are as follows: (1) The 14 cities differ significantly in their energy efficiency, with Shenyang, Dalian, Anshan, and Jinzhou enjoying the highest energy efficiency. Northwestern Liaoning and other heavy-industrial cities such as Fushun and Benxi belong to low-efficiency and high-consumption areas. (2) In areas with higher efficiency, the spatial patterns of the energy technical efficiency, energy consumption per unit of GDP, and energy utilization efficiency are, respectively, "π"-,"Ⅱ"- and "H"- shaped. Geographically, the energy utilization efficiency shows different trends from east to west and from north to south. Factors such as the binuclear structure of economic development have a major effect on this spatial pattern of energy efficiency. (3) Southeastern Liaoning enjoys a highly advanced industrial structure. Areas with a highly rational industrial structure form an “H” shape, with Shenyang and Dalian at the two poles. The urban agglomerations in middle and southern Liaoning have a highly concentrated industrial structure. (4) Overall, the coupling between energy efficiency and industrial structure is low in Liaoning, except for Shenyang and Dalian at both ends, where the coupling between an advanced industrial structure and energy efficiency is higher than in other cities.

This paper investigated several stages in the formation of the geopolitical influence of oil and gas, including the basis of its gestation, the means of transformation, and the formation and exercise of power. Based on this theoretical framework, a system for assessing the geopolitical influence of oil and gas was developed. This system is comprised of 13 indicators, each with its own method of measurement. Then 21 representative oil and gas importing, exporting, and transit countries were selected as assessment subjects. A quantitative assessment of the geopolitical influence of oil and gas in the selected countries was carried out using the proposed indicators, and factor analysis was used to obtain the main factors of these indicators and the composite score of each country. The empirical results showed that the 13 indicators could be summarized into five main factors, in the order of contribution rate, specifically comprehensive national strength, energy, transportation, risk, and geopolitics, each with its own variance contribution rate. Results of the assessment indicated that the selected countries could be classified into five categories in terms of oil and gas geopolitical influence: strong, relatively strong, moderate, relatively weak, and weak.

China is a disaster prone country, and a comprehensive understanding of change of disasters is very important for China’s agricultural development. In this study, statistical techniques and geographic information system tools are employed to quantify the main agriculture disasters changes and effects on grain production in China during the period of 1990-2011. The results show that China’s grain production was severely affected by disasters including drought, flood, hail, frost and typhoon. The annual area covered by these disasters reached up to 48.7×10⁶ ha during the study period, which accounted for 44.8% of the total sown area, and about 55.1% of the per unit area grain yield change was caused by disasters. In addition, all of the disasters showed high variability, different changing trends, and spatial distribution. Drought, flood, and hail showed significantly decreasing trends, while frost and typhoon showed increasing trends. Drought and flood showed gradual changes and were distributed across the country, and disasters became more diversified from north to south. Drought was the dominated disaster type in northern China, while flood was the most important disaster type in the southern part. Hail was mainly observed in central and northern China, and frost was mainly distributed in southern China. Typhoon was greatly limited to the southeast coast. Furthermore, the resilience of grain production of each province was quite different, especially in several major grain producing areas, such as Shandong, Liaoning, Jilin and Jiangsu, where grain production was seriously affected by disasters. One reason for the difference of resilience of grain production was that grain production was marginalized in developed provinces when the economy underwent rapid development. For China’s agricultural development and grain security, we suggest that governments should place more emphasis on grain production, and invest more money in disaster prevention and mitigation, especially in the major grain producing provinces.

The sediment discharge from the Yangtze River Basin has a stepwise decreasing trend in recent years. The impounding of the Three Gorges Reservoir exacerbated this decreasing trend and affected the change of the suspended sediment concentration (SSC) in the Yangtze River Estuary through the transmission effect. The SSC data of the Yangtze River Estuary during 1959-2012 showed that: (1) The SSC in the South Branch of the Yangtze River in the estuary and in the off-shore sea area displayed decreasing trends and decreased less towards the sea. At the same time, the difference in decreasing magnitude between SSC and sediment discharge became bigger towards the sea. (2) For the North Branch the preferential flow did not change much but the SSC tended to decrease, which was mainly caused by the decrease of SSC in the South Branch and China East Sea. (3) Due to the decreased runoff and the relatively strengthened tide, the peak area of the SSC in the bar shoal section in 2003-2012 moved inward for about 1/6 longitude unit compared with that in 1984-2002, and the inward-moving distance was in the order of flood season > annual average > dry season. (4) In the inlet of the South Passage, the SSC decreased mainly because the increase caused by resuspension and shore-groove exchange was less than the decrease caused by the sharp SSC decrease in the basin and the sea areas. The reverse was true in the middle section, where the SSC showed an increasing trend. (5) In the inlet of the North Passage, under the combined influence of decreased flow split and sediment split ratios, the decreased SSC in the basin and the sea area and decreased amount of resuspension, the SSC displayed a decreasing trend. In the middle section, because the increased amount caused by sediment going over the dyke was markedly more than the decreased amount caused by external environments, the SSC tended to increase. Holistically, the sharp decrease in sediment discharge caused synchronized SSC decreases in the Yangtze River Estuary. But there were still areas, where the SSC displayed increasing trends, indicating synchronicity and difference in the response of SSC to the sharp decrease in sediment discharge from the basin.

Evaluating the effectiveness of protecting nature assets is a vital component of responsive, pro-active management of protected areas. It is important to protect reserves of geological (including fossils) assets as a means to also protect national and natural heritage. At the end of 2013, 2669 protected areas had been established in China and 123 of them are reserves for nature assets. This paper builds an evaluation framework for protection effectiveness of these assets. The current elements and characteristics for effective protection are analyzed, along with an analysis of existing problems so as to construct a scientific approach to protect these assets. The influencing elements and characteristics for effective protection of nature asset reserves are included in an index evaluation system for effective protection in such areas, which contains four parts-target layer, evaluation criteria, evaluation indexes, and evaluation parameters-based on related documents and files. For the target layer, it includes nature asset evaluation and management evaluation. In the end, it is discussed how to build a comprehensive evaluation model and achieve an effective quantitative evaluation.

Depended on the analysis of ground snow situation, soil moisture loss speed and soil structure after planting crops of Mu Us Sandy Land remedied with feldspathic sandstone in the fallow period, it is concluded that feldspathic sandstone mixed with sand improved the sand stabilization in the governance of Mu Us Sandy Land in the fallow period. The sandy land remedied with feldspathic sandstone had big snow coverage, 25%-75% higher than normal sand; soil moisture losses slowed down, and moisture content rose by over 3 times; soil structure had been improved, and water stable aggregate content increased by 6.52%-18.04%; survival rate of protection forest increased to 85%; and ground flatness is less than 1%. The above conditions weakened sand rising conditions of Mu Us Sandy Land in the fallow period and formed two protective layers of snow cover and soil frozen layer under cold weather so as to prevent against wind erosion.

To control soil erosion and restore the degraded environment in the Loess Plateau, a large number of measures related to soil and water conservation have been employed that have profoundly affected catchment properties. This study constructed three indicators to characterize changes to the catchment slope, proposed both a method for a regression analysis of adjacent images and a sequence model, and applied multisource remotely sensed images and GIS spatial clustering analysis technologies to extract thematic information and comprehensively analyze the catchment change characteristics. The results indicate that the catchment slope properties changed significantly. At catchment scale, the average values of ARC, DVC and ART were 6.43%, 25.57% and 4.30%, respectively. There were six clustering types of catchment slope property changes. The maximum and minimum of the average similarities of the clustering types were 0.992 and 0.935. Each slope control measures had a distinct effect on catchment slope; the dominating factor of each clustering type was identified as: Type 1: D-VC, Type 2: D-VCLU, Type 3: D-LUVC, Type 4: D-TAVC, Type 5: D-TAC and Type 6: D-MFC. Type 5 and Type 1 covered the largest areas, respectively occupying 37.28% and 31.01%. Catchment slope property changes also had distinct types that depended on their geomorphological conditions. These findings provide a useful basis from which to further study catchment slope hydrological and soil erosion processes.

Before the emergence of modern modes of transport, the traditional road infrastructure was the major historical means of carrying out nationwide socio-economic exchange. However, the history of transport infrastructure has received little attention from researchers. Given this background, the work reported here examined the long-term development of transport networks in China. The national road network was selected for study and the 3500 years from 1600 BC to 1900 AD was chosen as the study period. Indicators were designed for the maturity level of road networks and an accessibility model was developed for the paths of the shortest distance. The evolution of the road network in China since the Shang Dynasty (1600 BC) was described and its major features were summarized to reveal long-term regularities. The maturity level of the road network and its accessibility was assessed and regions with good and poor networks were identified. The relationship between China’s natural, social, and economic systems and the road network were discussed. Our analysis shows that the road network in China has a number of long-term regularities. The continuously expanding road network follows a path of inland expansion especially towards the border areas. However, its coverage and accessibility are characterized by a core-peripheral configuration, which has close relationships with, not only the natural conditions, but also national defense and warfare. The centralization of national power, national land governance, postal transport, the transport of specialized cargos, and international trade are also related to the development of the road network. This research draws attention to the evolving regularities of transport networks.

The heating effect (or mass elevation effect, MEE) of the Tibetan Plateau (TP) is intense due to its massive body. Some studies have been undertaken on its role as the heat source in summer and its implications for Asian climate, but little has been known of the implications of its MEE for the distribution of mountain altitudinal belts (MABs). Using air temperature data observed and remotely sensed data, MAB/treeline data, and ASTER GDEM data, this paper compares the height of MABs and alpine treelines in the main TP and the surrounding mountains/lowland and explains the difference from the point of view of MEE. The results demonstrate: 1) at same elevation, air temperature and the length of growing season gradually increase from the eastern edge to the interior TP, e.g., at 4500 m (corresponding to the mean altitude of the TP), the monthly mean temperature is 3.58°C higher (April) to 6.63°C higher (June) in the interior plateau than in the Sichuan Basin; the 10°C isotherm for the warmest month goes upward from the edge to the interior of the plateau, at 4000 m in the Qilian Mts. and the eastern edges of the plateau, and up to 4600-5000 m in Lhasa and Zuogong; the warmth index at an altitude of 4500 m can be up to 15°C·month in the interior TP, but much lower at the eastern edges. 2) MABs and treeline follow a similar trend of rising inwards: dark-coniferous forest is 1000-1500 m higher and alpine steppe is about 700-900 m higher in the interior TP than at the eastern edges.

Based on the daily reanalysis data released by NCEP/NCAR and the daily precipitation of 753 Chinese stations from May to August during the period of 1960 to 2012, the statistical characteristics of the cold vortex in northeastern China were analyzed. In addition, the strength index, which described the characteristics of the vortex consistently and frequently, and the geographical distribution were given by continuous anomalies of circulation. Based on this index, the activity routines of the cold vortex, characteristics of atmospheric circulation, and their effects on precipitation in northeastern China were analyzed. The results show that: the activities of the cold vortex exhibit remarkable features of annual and interdecadal oscillation, and the vortex high frequency and its characteristics of atmospheric circulation are described more accurately by the strength index of the cold vortex, which shows a high correspondence with the vortex precipitation during early summer and midsummer in the northeast. In strong (weak) vortex years, the general circulation in the middle and high latitudes of Eurasia is to the advantage (disadvantage) of the formation, development and maintenance of the cold vortex, thus it is easy (difficult) to form the circulation which is beneficial to transmit vapor from south to north during the period of July to August. Blocking over the Ural Mountains prevails (does not prevail) in early summer, and blocking over the Sea of Okhotsk prevails (does not prevail) in midsummer. Areas where the subtropical high is too small (large) and moves toward the north too late (early) are better (worse) for the maintenance of the cold vortex in northeastern China.

The objective of this study was to investigate the concentration and spatial distribution patterns of 9 potentially toxic heavy metal elements (As, Cd, Co, Cr, Pb, Cu, Zn, Mn, and Ni) in road dust in the Bayan Obo Mining Region in Inner Mongolia, China. Contamination levels were evaluated using the geoaccumulation index and the enrichment factor. Human health risks for each heavy metal element were assessed using a human exposure model. Results showed that the dust contained significantly elevated heavy metal elements concentrations compared with the background soil. The spatial distribution pattern of all tested metals except for As coincided with the locations of industrial areas while the spatial distribution of As was associated with domestic sources. The contamination evaluation indicated that Cd, Pb, and Mn in road dust mainly originated from anthropogenic sources with a rating of “heavily polluted” to “extremely polluted,” whereas the remaining metals originated from both natural and anthropogenic sources with a level of “moderately polluted”. The non-cancer health risk assessment showed that ingestion was the primary exposure route for all metals in the road dust and that Mn, Cr, Pb, and As were the main contributors to non-cancer risks in both children and adults. Higher HI values were calculated for children (HI=1.89), indicating that children will likely experience higher health risks compared with adults (HI=0.23). The cancer risk assessment showed that Cr was the main contributor, with cancer risks which were 2-3 orders of magnitude higher than those for other metals. Taken in concert, the non-cancer risks posed by all studied heavy metal elements and the cancer risks posed by As, Co, Cr, Cd, and Ni to both children and adults in Bayan Obo Mining Region fell within the acceptable range.

Slope spectrum has been proved to be a significant methodology in revealing geomorphological features in the study of Chinese loess terrain. The determination of critical areas in deriving slope spectra is an indispensable task. Along with the increase in the size of the study area, the derived spectra are becoming more and more alike, such that their differences can be ignored in favor of a standard. Subsequently, the test size is defined as the Slope Spectrum Critical Area (SSCA). SSCA is not only the foundation of the slope spectrum calculation but also, to some extent, a reflection of geomorphological development of loess relief. High resolution DEMs are important in extracting the slope spectrum. A set of 48 DEMs with different landform areas of the Loess Plateau in northern Shaanxi province was selected for the experiment. The spatial distribution of SSCA is investigated with a geo-statistical analysis method, resulting in values ranging from 6.18 km² to 35.1 km². Primary experimental results show that the spatial distribution of SSCA is correlated with the spatial distribution of the soil erosion intensity, to a certain extent reflecting the terrain complexity. The critical area of the slope spectrum presents a spatial variation trend of weak-strong-weak from north to south. Four terrain parameters, gully density, slope skewness, terrain driving force (T_d) and slope of slope (SOS), were chosen as indicators. There exists a good exponential function relationship between SSCA and gully density, terrain driving force (T_d) and SOS and a logarithmic function relationship between SSCA and slope skewness. Slope skewness increases, and gully density, terrain driving force and SOS decrease with increasing SSCA. SSCA can be utilized as a discriminating factor to identify loess landforms, in that spatial distributions of SSCA and the evolution of loess landforms are correlative. Following the evolution of a loess landform from tableland to gully-hilly region, this also proves that SSCA can represent the development degree of local landforms. The critical stable regions of the Loess Plateau represent the degree of development of loess landforms. Its chief significance is that the perception of stable areas can be used to determine the minimal geographical unit.

During the past decade, great efforts have been made to boost the land use transformation in the Loess Plateau, especially for reducing soil erosion by vegetation restoration measures. The Grain-for-Green project (GFG) is the largest ecological rehabilitation program in China, which has a positive impact on the vegetation restoration and sustainable development for the ecologically fragile region of west China. Based on the Landsat TM/ETM images for three time periods (2000, 2005 and 2010), this study applied the GIS technology and a hill-slope analytical model to reveal the spatio-temporal evolutional patterns of returning slope farmland to grassland or woodland in Baota District, Yan’an city of Shaanxi province. Results showed that: (1) from 2000 to 2010, the area of farmland decreased by approximately 35,030 ha, which is the greatest decrease among all the land-use types, whereas grassland, woodland and construction land increased, of which grassland expanded rapidly by 26,380 ha. (2) The annual variation rate of land-use dynamics was 1.98% during the period 2000-2010, of which the rate was 1.05% for the 2000-2005 period and 2.92% for the 2005-2010 period, respectively. Over the past decade, returning farmland to woodland or pastures was the main source of increased grassland and woodland, and the reduction of farmland contributed to the increase in grassland and woodland by 97.39% and 85.28%, respectively. (3) As the terrain slope increases, farmland decreased and woodland and grassland increased significantly. Areas with a slope ranging from 15° to 25° and less than 15° were the focus of the GFG project, accounting for 85% of the total area of farmland reduction. Meanwhile, the reduction in farmland was significant and spatially correlated with the increase in woodland and grassland. (4) Between 2000 and 2010, the area of destruction of grass and trees in grasslands and woodlands for the reclamation of farmland was approximately 4596 ha. The area subject to the GFG policy was 4456 ha with a slope greater than 25° over the decade, but the area of farmland was still 10,357 ha in 2010. Our results indicate that there has still a great potential for returning the steep-slope farmlands to woodlands or grasslands in the Loess Plateau.

Detailed analysis of Land Use/Land Cover (LULC) using remote sensing data in complex irrigated basins provides complete profile for better water resource management and planning. Using remote sensing data, this study provides detailed land use maps of the Lower Chenab Canal irrigated region of Pakistan from 2005 to 2012 for LULC change detection. Major crop types are demarcated by identifying temporal profiles of NDVI using MODIS 250 m × 250 m spatial resolution data. Wheat and rice are found to be major crops in rabi and kharif seasons, respectively. Accuracy assessment of prepared maps is performed using three different techniques: error matrix approach, comparison with ancillary data and with previous study. Producer and user accuracies for each class are calculated along with kappa coefficients (K). The average overall accuracies for rabi and kharif are 82.83% and 78.21%, respectively. Producer and user accuracies for individual class range respectively between 72.5% to 77% and 70.1% to 84.3% for rabi and 76.6% to 90.2% and 72% to 84.7% for kharif. The K values range between 0.66 to 0.77 for rabi with average of 0.73, and from 0.69 to 0.74 with average of 0.71 for kharif. LULC change detection indicates that wheat and rice have less volatility of change in comparison with both rabi and kharif fodders. Transformation between cotton and rice is less common due to their completely different cropping conditions. Results of spatial and temporal LULC distributions and their seasonal variations provide useful insights for establishing realistic LULC scenarios for hydrological studies.

China has been experiencing an unprecedented urbanization process. In 2011, China’s urban population reached 691 million with an urbanization rate of 51.27%. Urbanization level is expected to increase to 70% in China in 2030, reflecting the projection that nearly 300 million people would migrate from rural areas to urban areas over this period. At the same time, the total fertility rate of China’s population is declining due to the combined effect of economic growth, environmental carrying capacity, and modern social consciousness. The Chinese government has loosened its “one-child policy” gradually by allowing childbearing couples to have the second child as long as either of them is from a one-child family. In such rapidly developing country, the natural growth and spatial migration will consistently reshape spatial pattern of population. An accurate prediction of the future spatial pattern of population and its evolution trend are critical to key policy-making processes and spatial planning in China including urbanization, land use development, ecological conservation and environmental protection. In this paper, a top-down method is developed to project the spatial distribution of China’s future population with considerations of both natural population growth at provincial level and the provincial migration from 2010 to 2050. Building on this, the spatial pattern and evolution trend of Chinese provincial population are analyzed. The results suggested that the overall spatial pattern of Chinese population will be unlikely changed in next four decades, with the east area having the highest population density and followed by central area, northeast and west area. Four provinces in the east, Shanghai, Beijing, Tianjin and Jiangsu, will remain the top in terms of population density in China, and Xinjiang, Qinghai and Tibet will continue to have the lowest density of population. We introduced an index system to classify the Chinese provinces into three categories in terms of provincial population densities: Fast Changing Populated Region (FCPR), Low Changing Populated Region (LCPR) and Inactive Populated Region (IPR). In the FCPR, China’s population is projected to continue to concentrate in net immigration leading type (NILT) area where receives nearly 99% of new accumulated floating population. Population densities of Shanghai, Beijing, Zhejiang will peak in 2030, while the population density in Guangdong will keep increasing until 2035. Net emigration leading type (NELT) area will account for 75% of emigration population, including Henan, Anhui, Chongqing and Hubei. Natural growth will play a dominant role in natural growth leading type area, such as Liaoning and Shandong, because there will be few emigration population. Due to the large amount of moving-out labors and gradually declining fertility rates, population density of the LCPR region exhibits a downward trend, except for Fujian and Hainan. The majority of the western provinces will be likely to remain relatively low population density, with an average value of no more than 100 persons per km².

This study has revealed spatial-temporal changes in Recreational Business Districts (RBDs) in Beijing and examined the relationship between the location of urban RBDs and traffic conditions, resident and tourist density, scenic spots, and land prices. A more reasonable classification of urban RBDs (LSC, CPS, and ULA) is also proposed. Quantitative methods such as Gini Coefficient, Spatial Interpolation, Kernel Density Estimation, and Geographical Detector were employed to collect and analyze the data from three types of urban RBDs in Beijing in 1990, 2000, and 2014, respectively, and the spatial-temporal patterns as well as the distribution characteristics of urban RBDs were analyzed using ArcGIS software. It was concluded that (1) both the number and scale of urban RBDs in Beijing have been expanding and the trend for all types of urban RBDs in Beijing to be spatially agglomerated is continuing; (2) the spatial-temporal evolution pattern of urban RBDs in Beijing is “single-core agglomeration-dual-core agglomeration-multi-core diffusion”; and (3) urban RBDs were always located in areas with low traffic density, tourist attractions, high resident and tourist population density, and relatively high land valuations; these factors also affect the scale size of RBDs.