Lithospheric structure beneath the northeastern Qinghai-Xizang Plateau is of vital significance for studying the geodynamic processes of crustal thickening and expansion of the Qinghai-Xizang Plateau. We conducted a j...Lithospheric structure beneath the northeastern Qinghai-Xizang Plateau is of vital significance for studying the geodynamic processes of crustal thickening and expansion of the Qinghai-Xizang Plateau. We conducted a joint inversion of receiver functions and surface wave dispersions with P-wave velocity constraints using data from the Chin Array Ⅱ temporary stations deployed across the Qinghai-Xizang Plateau. Prior to joint inversion, we applied the H-κ-c method(Li JT et al., 2019) to the receiver function data in order to correct for the back-azimuthal variations in the arrival times of Ps phases and crustal multiples caused by crustal anisotropy and dipping interfaces. High-resolution images of vS, crustal thickness, and vP/vSstructures in the Qinghai-Xizang Plateau were simultaneously derived from the joint inversion. The seismic images reveal that crustal thickness decreases outward from the Qinghai-Xizang Plateau. The stable interiors of the Ordos and Alxa blocks exhibited higher velocities and lower crustal vP/vSratios. While, lower velocities and higher vP/vSratios were observed beneath the Qilian Orogen and Songpan-Ganzi terrane(SPGZ), which are geologically active and mechanically weak, especially in the mid-lower crust.Delamination or thermal erosion of the lithosphere triggered by hot asthenospheric flow contributes to the observed uppermost mantle low-velocity zones(LVZs) in the SPGZ. The crustal thickness, vS, and vP/vSratios suggest that whole lithospheric shortening is a plausible mechanism for crustal thickening in the Qinghai-Xizang Plateau, supporting the idea of coupled lithospheric-scale deformation in this region.展开更多
The margin of the Tibetan Plateau of Southwest China is one of the most seismically active regions of China and is the location of the China Seismic Experimental Site(CSES).Many studies have developed seismic velocity...The margin of the Tibetan Plateau of Southwest China is one of the most seismically active regions of China and is the location of the China Seismic Experimental Site(CSES).Many studies have developed seismic velocity models of Southwest China,but few have compared and evaluated these models which is important for further model improvement.Thus,we compared six published seismic shear-wave velocity models of Southwest China on absolute velocity and velocity perturbation patterns.The models are derived from different types of data(e.g.,surface waves from ambient noise and earthquakes,body-wave travel times,receiver functions) and inversion methods.We interpolated the models into a uniform horizontal grid(0.5° × 0.5°) and vertically sampled them at 5,10,20,30,40,and 60 km depths.We found significant differences between the six models.Then,we selected three of them that showed greater consistency for further comparison.Our further comparisons revealed systematic biases between models in absolute velocity that may be related to different data types.The perturbation pattern of the model is especially divergent in the shallow part,but more consistent in the deep part.We conducted synthetic and inversion tests to explore possible causes and our results imply that systematic differences between the data,differences in methods,and other factors may directly affect the model.Therefore,the Southwest China velocity model still has considerable room for improvement,and the impact of inconsistency between different data types on the model needs further research.Finally,we proposed a new reference shear-wave velocity model of Southwest China(SwCM-S1.0) based on the three selected models with high consistency.We believe that this model is a better representation of more robust features of the models that are based on different data sets.展开更多
The Chinese mainland is subject to complicated plate interactions that give rise to its complex structure and tectonics. While several seismic velocity models have been developed for the Chinese mainland, apparent dis...The Chinese mainland is subject to complicated plate interactions that give rise to its complex structure and tectonics. While several seismic velocity models have been developed for the Chinese mainland, apparent discrepancies exist and, so far, little effort has been made to evaluate their reliability and consistency. Such evaluations are important not only for the application and interpretation of model results but also for future model improvement. To address this problem, here we compare five published shear-wave velocity models with a focus on model consistency. The five models were derived from different datasets and methods (i.e., body waves, surface waves from earthquakes, surface waves from noise interferometry, and full waves) and interpolated into uniform horizontal grids (0.5° × 0.5°) with vertical sampling points at 5 km, 10 km, and then 20 km intervals to a depth of 160 km below the surface, from which we constructed an averaged model (AM) as a common reference for comparative study. We compare both the absolute velocity values and perturbation patterns of these models. Our comparisons show that the models have large (> 4%) differences in absolute values, and these differences are independent of data coverage and model resolution. The perturbation patterns of the models also show large differences, although some of the models show a high degree of consistency within certain depth ranges. The observed inconsistencies may reflect limited model resolution but, more importantly, systematic differences in the datasets and methods employed. Thus, despite several seismic models being published for this region, there is significant room for improvement. In particular, the inconsistencies in both data and methodologies need to be resolved in future research. Finally, we constructed a merged model (ChinaM-S1.0) that incorporates the more robust features of the five published models. As the existing models are constrained by different datasets and methods, the merged model serves as a new type of reference model that incorporates the common features from the joint datasets and methods for the shear-wave velocity structure of the Chinese mainland lithosphere.展开更多
Using the southern limit of snowfall recorded in Chinese documents, chronologies of tree-ring width, and tree-ring stable oxygen isotope(δ^(18)O), the annual temperature anomaly in southern China during 1850–200...Using the southern limit of snowfall recorded in Chinese documents, chronologies of tree-ring width, and tree-ring stable oxygen isotope(δ^(18)O), the annual temperature anomaly in southern China during 1850–2009 is reconstructed using the method of signal decomposition and synthesis. The results show that the linear trend was 0.47℃(100 yr)^(-1)over 1871–2009,and the two most rapid warming intervals occurred in 1877–1938 and 1968–2007, at rates of 0.125℃(10 yr)^(-1)and 0.258℃(10 yr)^(-1), respectively. The decadal variation shows that the temperature in the moderate warm interval of the 1910s–1930s was notably lower than that of the 1980s–2000s, which suggests that the warming since the 1980s was unprecedented for the past 160 years, though a warming hiatus existed in the 2000s. Additionally, there was a rapid cooling starting from the 1860s,followed by a cold interval until the early 1890s, with the coldest years in 1892 and 1893. A slight temperature decline was also found from the 1940s to the late 1960s. This study provides an independent case to validate the global warming for the past 160 years and its hiatus recently, because the proxy data are not affected by urbanization.展开更多
Based on observations made during recent decades, reconstructed precipitation for the period A.D. 1736-2000, dry-wet index data for A.D. 500-2000, and a 1000-yr control simulation using the Community Earth System Mode...Based on observations made during recent decades, reconstructed precipitation for the period A.D. 1736-2000, dry-wet index data for A.D. 500-2000, and a 1000-yr control simulation using the Community Earth System Model with fixed pre-industrial external forcing, the decadal variability of summer precipitation over eastem China is stud- ied. Power spectrum analysis shows that the dominant cycles for the decadal variation of summer precipitation are: 22-24 and quasi-70 yr over the North China Plain; 32-36, 44-48, and quasi-70 yr in the Jiang-Huai area; and 32-36 and 4448 yr in the Jiang-Nan area. Bandpass decomposition from observation, reconstruction, and simulation re- veals that the variability of summer precipitation over the North China Plain, Jiang-Huai area, and Jiang-Nan area, at scales of 20-35, 35-50, and 50-80 yr, is not consistent across the entire millennium. We also find that the warm (cold) phase of the Pacific Decadal Oscillation generally corresponds to dry (wet) conditions over the North China Plain, but wet (dry) conditions in the Jiang-Nan area, from A.D. 1800, when the PDO became strengthened. However, such a correspondence does not exist throughout the entire last millennium. Data-model comparison sug- gests that these decadal oscillations and their temporal evolution over eastem China, including the decadal shifts in the spatial pattem of the precipitation anomaly observed in the late 1970s, early 1990s, and early 2000s, might result from internal variability of the climate system.展开更多
Jinping Underground laboratory for Nuclear Astrophysics(JUNA) will take the advantage of the ultra-low background of CJPL lab and high current accelerator based on an ECR source and a highly sensitive detector to dire...Jinping Underground laboratory for Nuclear Astrophysics(JUNA) will take the advantage of the ultra-low background of CJPL lab and high current accelerator based on an ECR source and a highly sensitive detector to directly study for the first time a number of crucial reactions occurring at their relevant stellar energies during the evolution of hydrostatic stars. In its first phase, JUNA aims at the direct measurements of^(25)Mg(p,γ)^(26)Al,^(19)F(p,α)^(16)O,^(13)C(α,n)^(16)O and ^(12)C(α,γ)^(16)O reactions. The experimental setup,which includes an accelerator system with high stability and high intensity, a detector system, and a shielding material with low background, will be established during the above research. The current progress of JUNA will be given.展开更多
The temperature variability over multidecadal and longer timescales(e.g., the cold epochs in the late 15 th, 17 th, and early 19 th centuries) is significant and dominant in the millennium-long, large-scale reconstruc...The temperature variability over multidecadal and longer timescales(e.g., the cold epochs in the late 15 th, 17 th, and early 19 th centuries) is significant and dominant in the millennium-long, large-scale reconstructions and model simulations;however, their temporal patterns in the reconstructed and simulated temperature series are not well understood and require a detailed assessment and comparison. Here, we compare the reconstructed and simulated temperature series for the Northern Hemisphere(NH) at multidecadal and longer-term timescales(>30 years) by evaluating their covariance, climate sensitivity and amplitude of temperature changes. We found that covariances between different reconstructions or between reconstructions and simulations are generally high for the whole period of 850–1999 CE, due to their similar long-term temporal patterns. However,covariances between different reconstructions or between reconstructions and simulations steadily decline as time series extends further back in time, becoming particularly small during Medieval times. This is related to the large uncetainties in the reconstructions caused by the decreased number of proxy records and sample duplication during the pre-instrumental periods.Reconstructions based solely on tree-ring data show higher skill than multiproxy reconstructions in capturing the amplitude of volcanic cooling simulated by models. Meanwhile, climate models have a shorter recovery(i.e., lag) in response to the cooling caused by volcanic eruptions and solar activity minima, implying the lack of some important feedback mechanisms between external forcing and internal climate processes in climate models. Amplitudes of temperature variations in the latest published tree-ring reconstructions are comparable to those of the multiproxy reconstructions. We found that the temperature difference between the Medieval Climate Anomaly(950–1250 CE) and the Little Ice Age(1450–1850 CE) is generally larger in proxybased reconstructions than in model simulations, but the reason is unclear.展开更多
Using modern census and environmental factor data,this study first identified the environmental factors that significantly affect the population distribution through Geodetector analysis and then constructed a populat...Using modern census and environmental factor data,this study first identified the environmental factors that significantly affect the population distribution through Geodetector analysis and then constructed a population spatial distribution model based on the random forest regression algorithm.Finally,with this model and historical population data that were examined and corrected by historians,gridded population distributions with a spatial resolution of 10 km by 10 km in the traditional cultivated region of China(TCRC,hereafter)were reconstructed for six time slices from 1776 to 1953.Using the reconstruction dataset,the spatiotemporal characteristics of the population distribution were depicted.The results showed that(1)the environmental factors that significantly affected the population density differences among counties in the TCRC mainly consisted of elevation,slope,relief amplitude,distances to the nearest prefectural and provincial capitals,distance to the nearest river and the climatology moisture index.(2)Using the census data of 1934 counties in the TCRC in 2000 and the abovementioned environmental factor data,a random forest regression algorithm-based population spatial distribution model was constructed.Its determination coefficient(R^(2))is 0.81.In 88.4%of the counties(districts),the relative errors of the model predictions were less than 50%.(3)From 1776 to 1953,the total population in the study area showed an uptrend.Prior to 1851,the population increased mainly in the Yangtze River Delta.During this period,the number of grid cells in which the population densities were greater than 500 persons per km^(2) increased from 292 to 683.From 1851 to 1953,the population increased extensively across the study area.In the North China Plain and the Pearl River Delta,the number of grid cells in which the population densities were greater than 500 persons per km^(2) increased from 36 to 88 and from 4 to 35,respectively.The spatial clustering pattern of the population distribution varied temporally.The potential reasons included the shifts in economic development hot spots,traditional beliefs,wars,famine,and immigration policies.(4)Between our reconstructions and the HYDE dataset,there are large differences in the data sources,selected environmental factors and modeling methods.As a consequence,in comparison to our reconstructions,there were fewer populations in the eastern area and more populations in the western area from 1776 to 1851 and more populations in urban areas and fewer populations in rural areas after 1851 in the HYDE dataset.展开更多
基金supported by the Natural Science Basic Research Program of Shaanxi(No.2023-JC-QN-0306)the Special Fund of the Institute of Geophysics,China Earthquake Administration(No.DQJB21B32)the National Natural Science Foundation of China(No.42174069).
文摘Lithospheric structure beneath the northeastern Qinghai-Xizang Plateau is of vital significance for studying the geodynamic processes of crustal thickening and expansion of the Qinghai-Xizang Plateau. We conducted a joint inversion of receiver functions and surface wave dispersions with P-wave velocity constraints using data from the Chin Array Ⅱ temporary stations deployed across the Qinghai-Xizang Plateau. Prior to joint inversion, we applied the H-κ-c method(Li JT et al., 2019) to the receiver function data in order to correct for the back-azimuthal variations in the arrival times of Ps phases and crustal multiples caused by crustal anisotropy and dipping interfaces. High-resolution images of vS, crustal thickness, and vP/vSstructures in the Qinghai-Xizang Plateau were simultaneously derived from the joint inversion. The seismic images reveal that crustal thickness decreases outward from the Qinghai-Xizang Plateau. The stable interiors of the Ordos and Alxa blocks exhibited higher velocities and lower crustal vP/vSratios. While, lower velocities and higher vP/vSratios were observed beneath the Qilian Orogen and Songpan-Ganzi terrane(SPGZ), which are geologically active and mechanically weak, especially in the mid-lower crust.Delamination or thermal erosion of the lithosphere triggered by hot asthenospheric flow contributes to the observed uppermost mantle low-velocity zones(LVZs) in the SPGZ. The crustal thickness, vS, and vP/vSratios suggest that whole lithospheric shortening is a plausible mechanism for crustal thickening in the Qinghai-Xizang Plateau, supporting the idea of coupled lithospheric-scale deformation in this region.
基金supported by Open Fund from Sino Probe Laboratory (No. Sinoprobe Lab 202201)the National Natural Science Foundation of China (No. U1939204)the Special Fund of the Institute of Geophysics, China Earthquake Administration (No. DQJB21B32)
文摘The margin of the Tibetan Plateau of Southwest China is one of the most seismically active regions of China and is the location of the China Seismic Experimental Site(CSES).Many studies have developed seismic velocity models of Southwest China,but few have compared and evaluated these models which is important for further model improvement.Thus,we compared six published seismic shear-wave velocity models of Southwest China on absolute velocity and velocity perturbation patterns.The models are derived from different types of data(e.g.,surface waves from ambient noise and earthquakes,body-wave travel times,receiver functions) and inversion methods.We interpolated the models into a uniform horizontal grid(0.5° × 0.5°) and vertically sampled them at 5,10,20,30,40,and 60 km depths.We found significant differences between the six models.Then,we selected three of them that showed greater consistency for further comparison.Our further comparisons revealed systematic biases between models in absolute velocity that may be related to different data types.The perturbation pattern of the model is especially divergent in the shallow part,but more consistent in the deep part.We conducted synthetic and inversion tests to explore possible causes and our results imply that systematic differences between the data,differences in methods,and other factors may directly affect the model.Therefore,the Southwest China velocity model still has considerable room for improvement,and the impact of inconsistency between different data types on the model needs further research.Finally,we proposed a new reference shear-wave velocity model of Southwest China(SwCM-S1.0) based on the three selected models with high consistency.We believe that this model is a better representation of more robust features of the models that are based on different data sets.
基金supportedby the Special Fund of the Instituteof Geophysics,China Earthquake Administration(Grant No.DQJB21B32)the National Natural Science Foundation of China(No.U1939204).
文摘The Chinese mainland is subject to complicated plate interactions that give rise to its complex structure and tectonics. While several seismic velocity models have been developed for the Chinese mainland, apparent discrepancies exist and, so far, little effort has been made to evaluate their reliability and consistency. Such evaluations are important not only for the application and interpretation of model results but also for future model improvement. To address this problem, here we compare five published shear-wave velocity models with a focus on model consistency. The five models were derived from different datasets and methods (i.e., body waves, surface waves from earthquakes, surface waves from noise interferometry, and full waves) and interpolated into uniform horizontal grids (0.5° × 0.5°) with vertical sampling points at 5 km, 10 km, and then 20 km intervals to a depth of 160 km below the surface, from which we constructed an averaged model (AM) as a common reference for comparative study. We compare both the absolute velocity values and perturbation patterns of these models. Our comparisons show that the models have large (> 4%) differences in absolute values, and these differences are independent of data coverage and model resolution. The perturbation patterns of the models also show large differences, although some of the models show a high degree of consistency within certain depth ranges. The observed inconsistencies may reflect limited model resolution but, more importantly, systematic differences in the datasets and methods employed. Thus, despite several seismic models being published for this region, there is significant room for improvement. In particular, the inconsistencies in both data and methodologies need to be resolved in future research. Finally, we constructed a merged model (ChinaM-S1.0) that incorporates the more robust features of the five published models. As the existing models are constrained by different datasets and methods, the merged model serves as a new type of reference model that incorporates the common features from the joint datasets and methods for the shear-wave velocity structure of the Chinese mainland lithosphere.
基金supported by the Strategic Priority Research Program of the Chinese Academy of Sciences(Grant No.XDA05090104)the National Natural Science Foundation of China(Grant Nos.41430528 and 41671201)
文摘Using the southern limit of snowfall recorded in Chinese documents, chronologies of tree-ring width, and tree-ring stable oxygen isotope(δ^(18)O), the annual temperature anomaly in southern China during 1850–2009 is reconstructed using the method of signal decomposition and synthesis. The results show that the linear trend was 0.47℃(100 yr)^(-1)over 1871–2009,and the two most rapid warming intervals occurred in 1877–1938 and 1968–2007, at rates of 0.125℃(10 yr)^(-1)and 0.258℃(10 yr)^(-1), respectively. The decadal variation shows that the temperature in the moderate warm interval of the 1910s–1930s was notably lower than that of the 1980s–2000s, which suggests that the warming since the 1980s was unprecedented for the past 160 years, though a warming hiatus existed in the 2000s. Additionally, there was a rapid cooling starting from the 1860s,followed by a cold interval until the early 1890s, with the coldest years in 1892 and 1893. A slight temperature decline was also found from the 1940s to the late 1960s. This study provides an independent case to validate the global warming for the past 160 years and its hiatus recently, because the proxy data are not affected by urbanization.
基金Supported by the National Natural Science Foundation of China(41430528 and 41471171)Institute of Geographic Sciences and Natural Resources Research,Chinese Academy of Sciences(TSYJS04,2014RC101,and 2015RC101)
文摘Based on observations made during recent decades, reconstructed precipitation for the period A.D. 1736-2000, dry-wet index data for A.D. 500-2000, and a 1000-yr control simulation using the Community Earth System Model with fixed pre-industrial external forcing, the decadal variability of summer precipitation over eastem China is stud- ied. Power spectrum analysis shows that the dominant cycles for the decadal variation of summer precipitation are: 22-24 and quasi-70 yr over the North China Plain; 32-36, 44-48, and quasi-70 yr in the Jiang-Huai area; and 32-36 and 4448 yr in the Jiang-Nan area. Bandpass decomposition from observation, reconstruction, and simulation re- veals that the variability of summer precipitation over the North China Plain, Jiang-Huai area, and Jiang-Nan area, at scales of 20-35, 35-50, and 50-80 yr, is not consistent across the entire millennium. We also find that the warm (cold) phase of the Pacific Decadal Oscillation generally corresponds to dry (wet) conditions over the North China Plain, but wet (dry) conditions in the Jiang-Nan area, from A.D. 1800, when the PDO became strengthened. However, such a correspondence does not exist throughout the entire last millennium. Data-model comparison sug- gests that these decadal oscillations and their temporal evolution over eastem China, including the decadal shifts in the spatial pattem of the precipitation anomaly observed in the late 1970s, early 1990s, and early 2000s, might result from internal variability of the climate system.
基金supported by the National Natural Science Foundation of China(Grant Nos.11490560 and 11321064)the National Basic Research Program of China(Grant No.2013CB834406)
文摘Jinping Underground laboratory for Nuclear Astrophysics(JUNA) will take the advantage of the ultra-low background of CJPL lab and high current accelerator based on an ECR source and a highly sensitive detector to directly study for the first time a number of crucial reactions occurring at their relevant stellar energies during the evolution of hydrostatic stars. In its first phase, JUNA aims at the direct measurements of^(25)Mg(p,γ)^(26)Al,^(19)F(p,α)^(16)O,^(13)C(α,n)^(16)O and ^(12)C(α,γ)^(16)O reactions. The experimental setup,which includes an accelerator system with high stability and high intensity, a detector system, and a shielding material with low background, will be established during the above research. The current progress of JUNA will be given.
基金supported by the National Key R&D Program of China (Grant No. 2017YFA0603302)supported by the National Natural Science Foundation of China (Grant Nos. 41888101 & 41602192 & 41977383)+4 种基金the Belmont Forum and JPI-Climate, Collaborative Research Action “INTEGRATE” (Grant No. 41661144008)support by the Youth Innovation Promotion Association Foundation of the Chinese Academy of Sciences (Grant No. 2018471)supported by the National Natural Science Foundation (Grant No. 41901095)supported by the National Natural Science Foundation (Grant No. 41877440)supported by Opening Fund of Key Laboratory of Desert and Desertification, Chinese Academy of Sciences (Grant No. KLDD-2019-04)。
文摘The temperature variability over multidecadal and longer timescales(e.g., the cold epochs in the late 15 th, 17 th, and early 19 th centuries) is significant and dominant in the millennium-long, large-scale reconstructions and model simulations;however, their temporal patterns in the reconstructed and simulated temperature series are not well understood and require a detailed assessment and comparison. Here, we compare the reconstructed and simulated temperature series for the Northern Hemisphere(NH) at multidecadal and longer-term timescales(>30 years) by evaluating their covariance, climate sensitivity and amplitude of temperature changes. We found that covariances between different reconstructions or between reconstructions and simulations are generally high for the whole period of 850–1999 CE, due to their similar long-term temporal patterns. However,covariances between different reconstructions or between reconstructions and simulations steadily decline as time series extends further back in time, becoming particularly small during Medieval times. This is related to the large uncetainties in the reconstructions caused by the decreased number of proxy records and sample duplication during the pre-instrumental periods.Reconstructions based solely on tree-ring data show higher skill than multiproxy reconstructions in capturing the amplitude of volcanic cooling simulated by models. Meanwhile, climate models have a shorter recovery(i.e., lag) in response to the cooling caused by volcanic eruptions and solar activity minima, implying the lack of some important feedback mechanisms between external forcing and internal climate processes in climate models. Amplitudes of temperature variations in the latest published tree-ring reconstructions are comparable to those of the multiproxy reconstructions. We found that the temperature difference between the Medieval Climate Anomaly(950–1250 CE) and the Little Ice Age(1450–1850 CE) is generally larger in proxybased reconstructions than in model simulations, but the reason is unclear.
基金This work was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences(Grant No.XDA19040101)the National Key Research and Development Program of China(Grant No.2017YFA0603301).
文摘Using modern census and environmental factor data,this study first identified the environmental factors that significantly affect the population distribution through Geodetector analysis and then constructed a population spatial distribution model based on the random forest regression algorithm.Finally,with this model and historical population data that were examined and corrected by historians,gridded population distributions with a spatial resolution of 10 km by 10 km in the traditional cultivated region of China(TCRC,hereafter)were reconstructed for six time slices from 1776 to 1953.Using the reconstruction dataset,the spatiotemporal characteristics of the population distribution were depicted.The results showed that(1)the environmental factors that significantly affected the population density differences among counties in the TCRC mainly consisted of elevation,slope,relief amplitude,distances to the nearest prefectural and provincial capitals,distance to the nearest river and the climatology moisture index.(2)Using the census data of 1934 counties in the TCRC in 2000 and the abovementioned environmental factor data,a random forest regression algorithm-based population spatial distribution model was constructed.Its determination coefficient(R^(2))is 0.81.In 88.4%of the counties(districts),the relative errors of the model predictions were less than 50%.(3)From 1776 to 1953,the total population in the study area showed an uptrend.Prior to 1851,the population increased mainly in the Yangtze River Delta.During this period,the number of grid cells in which the population densities were greater than 500 persons per km^(2) increased from 292 to 683.From 1851 to 1953,the population increased extensively across the study area.In the North China Plain and the Pearl River Delta,the number of grid cells in which the population densities were greater than 500 persons per km^(2) increased from 36 to 88 and from 4 to 35,respectively.The spatial clustering pattern of the population distribution varied temporally.The potential reasons included the shifts in economic development hot spots,traditional beliefs,wars,famine,and immigration policies.(4)Between our reconstructions and the HYDE dataset,there are large differences in the data sources,selected environmental factors and modeling methods.As a consequence,in comparison to our reconstructions,there were fewer populations in the eastern area and more populations in the western area from 1776 to 1851 and more populations in urban areas and fewer populations in rural areas after 1851 in the HYDE dataset.