—年华北饥荒的时空分布及其与气候灾害收成的关系

中国人民大学清史研究所,北京 100872

Institute of Qing History, Renmin University of China, Beijing 100872, China

作者简介 About authors

基于历史文献中的饥荒记载,提取1736—1911年176年华北境内饥荒事件2 950县次,根据饥荒严重程度等级,以加权平均法重建逐年饥荒指数序列,将其与同期气候变化、水旱灾害和粮食收成情况进行时空对比,结论如下:温度、降水与歉收、饥荒呈显著负相关,旱灾与歉收、饥荒的相关关系较水灾更为显著;回归分析显示歉收和饥荒对于旱灾的敏感性随清朝统治的衰落急剧上升,反映有效的社会应对有利于减少歉收和饥荒发生;核密度分析显示饥荒多发区主要分布在汾河谷地、河南中南部以及冀鲁豫交界处,与旱灾、歉收的空间对应关系较好,但河南中南部饥荒多发还可能与当地特殊的社会经济脆弱性有关。研究有助于深化对历史时期气候变化社会经济影响机制的科学认识,未来应进一步探讨历史饥荒的发生机制,并完善气候变化影响在社会系统中的传递路径。

关键词:历史气候变化;水旱灾害;粮食收成;饥荒;时空分布;华北

Based on records on famine kept in historical documents, 2950 famine events counted with “county time” in North China during 1736-1911 were collected, and annual famine index series in North China was reconstructed using weighted average method according to the severity of famine. Spatio-temporal distribution of famine and its relationship with climate, disaster and harvest was analyzed. The research shows that variation of temperature and precipitation has significant negative correlation with poor harvest and famine. The drought has a more significant correlation with harvest and famine than flood. In 1736-1795 when was at the peak of the Qing dynasty, the sensitivity of poor harvest and famine to drought was significantly lower than in 1796-1911 when was in decline, because effective social response to disaster could reduce the risk of poor harvest and famine. Using the kernel density estimation method, three main high risk areas of famine were indentified: Fen River Valley in south Shanxi Province, south-central He'nan Province, and the junction of Hebei, Shandong and He'nan Province. In general, there is significant corresponding relationship between the spatial distribution of drought, poor harvest and famine, however, abnormally high probability of famine in south-central He'nan Province might be related to local special socio-economic vulnerability. This study would enhance our understanding of the nature of interactions between climate change and socioeconomic system. Further research would focus on the mechanism by which famine originated in the past, and the pathway through which climatic impacts were delivered in human society.

Keywords:Historical climate change;Flood/drought;Grain harvest;Famine;Spatio-temporal distribution;North China.

本文引用格式

从原始记录可以看到,清代华北地区的饥荒发生带有强烈的季节性,多发于冬春,特别是灾害发生的次年春季青黄不接之时;而方志中记载的饥荒是以自然年为单位的,这就导致文字记录中饥荒的发生常常与灾害存在1年的时间差,也导致前一年冬季与次年春季的同一次饥荒事件可能在记录中被分为2条。由此对饥荒记录的发生时间进行适当修正——凡原始记录中明确记载为春季(麦收之前)发生的饥荒,其时间均向前修正1年,如前一年冬季亦有饥荒记录,则合并处理。最终获得饥荒记录2 950条,每1条记录对应该县级政区当年发生饥荒事件1次,记为1县次。

表1基于语义差异的饥荒等级划分标准

Table 1  Quantification of famine grade based on semantic differential method

注:历史资料中常见的“饿殍载道(道殣相望)”等词语需做具体分析;如果从字面理解,“饿殍”或“殣”都指饿死者,但如果联系上下文来看,这类描述在很多情况下都是对饥荒景象的夸张性描述,更多反映的是人口迁徙而非真正的大量死亡现象;因此如果没有确切的大量人口死亡记录作为参照,对其按第2级处理

对不同严重等级的饥荒事件频次分别赋以不同的权重,再进行累加,最终得到逐年华北饥荒指数,这一指数综合考虑了饥荒事件在2个维度上的严重程度(社会后果和波及范围)。公式如下:

用于与饥荒指数序列进行对比(主要采用双变量相关分析,以Pearson相关系数衡量关系强弱;并辅以线性回归分析,以回归系数判断敏感性差异)的气候变化、水旱灾害、粮食收成等时间序列来自已经发表的研究成果,主要包括以下几条:

考虑到某年某县域内发生的饥荒事件在空间上可以被抽象为点要素(对应县治所在地),以核密度估计法对研究时段内饥荒在华北境内的空间分布特征进行刻画,对饥荒多发区进行识别,并与灾害和收成进行空间对比。核密度估计是一种常见的空间分析手段,特别有利于对离散点的空间分布趋势进行呈现,其公式如下:

图11736—1911年华北的气候变化、水旱灾害、粮食收成与饥荒指数序列

Fig.1Annual famine index and its relationship with climate change, flood/drought, harvest in North China from 1736 to 1911

尽管整个序列并未呈现出显著的线性上升趋势,但不难发现随着时间推移,清代华北社会遭受饥荒的严重程度在持续提升。根据清代不同皇帝的在位时间,可将序列划分为3个大致相等的时段:1736—1795年(乾隆)、1796—1850年(嘉庆、道光)和1851—1911年(咸丰—宣统),分别对应清代的兴盛、衰落与崩溃阶段。在这3个时段中,饥荒指数均值分别为3.35、3.76和4.30,逐次升高;每个时段各出现1个显著峰值(分别为1785年、1813年和1877年),3个峰值的饥荒指数值也是逐次升高。

表21736—1911年华北气候、灾害与收成、饥荒的相关性

Table 2  Correlation coefficients between climate, disaster, harvest and famine in North China from 1736 to 1911

注:**表示相关关系显著性水平为0.01,*表示相关关系显著性水平为0.05

表3不同时段水旱灾害与收成、饥荒的相关性

Table 3  Correlation coefficients between flood, drought, harvest and famine in different periods

图21736—1911年华北水旱灾害、歉收、饥荒核密度分布

Fig.2Kernel density distribution of flood, drought, poor harvest and famine in North China from 1736 to 1911

饥荒多发区在空间上的集聚受到多方面因素的综合作用。极端灾害(特别是旱灾)→粮食歉收→饥荒的传递机制可以很好地解释上述第二类多发区的形成。海河水系下游地区虽然由于水灾集中而存在较高的歉收概率,但发生饥荒的可能性偏低,除了与这里极端旱灾少发有关之外,靠近京师(北京)的有利地理位置,也使其能够获得较为充分的救济,从而降低发生饥荒的概率;同理,鲁西北地区发生饥荒的概率以东、西两侧更高,而中部稍低,可能与这里临近省会济南和京杭运河仓储重地德州有关。

以《中国三千年气象记录总集》为主要资料源,提取以县级政区为单元的历史饥荒信息,采用语义差异法对文本描述的饥荒严重程度进行分级,以加权平均方法重建1736—1911年华北饥荒指数序列,将其与同期气候变化、水旱灾害和粮食收成进行时空对比,得到主要结论如下:

(1)研究时段内饥荒指数并未呈现显著的线性趋势,但3个典型时段(1736—1795年、1796—1850年和1851—1911年)中饥荒严重程度持续提升。核密度分析显示饥荒多发区主要分布在汾河谷地、河南中南部以及冀鲁豫交界处。

(2)相关分析结果显示研究时段内温度、降水变化与秋粮歉收指数、饥荒指数均为显著负相关;秋粮歉收指数与饥荒指数之间的相关系数达到0.88,说明粮食收成在气候影响传递到饥荒的过程中是一个关键的中间环节;旱灾指数则与歉收和饥荒均为显著正相关。回归分析结果显示3个典型时段中,歉收和饥荒对于旱灾的敏感性随着时间推移在不断增强。这种差异可能来自不同时段社会响应能力的差别,有效应对有利于减少灾害背景下歉收和饥荒发生。

(3)饥荒多发区的空间分布特征与旱灾、歉收的对应关系较好,但河南中南部在灾害相对少发、歉收发生概率也不高的前提下成为饥荒极端多发区,可能与当地特殊的社会经济脆弱性有关。

未来进一步开展研究的主要方向有2个:

(1)本文主要是从气候、灾害、收成与饥荒的关系角度分析饥荒成因,但历史时期饥荒的成因机制十分复杂,受到多方面因素(特别是社会经济因素)的综合影响,未来可以围绕典型时段、案例及区域(如本文中河南中南部饥荒极端多发区),更为全面地探讨饥荒的成因。

(2)饥荒是气候影响人类社会经济系统的链条上关键的一环,可能成为人口迁徙、社会动乱等事件的先声,未来可以基于历史气候变化影响与响应研究框架,将更多的定量重建序列(如粮价、人口、动乱)加入进来,最终形成一个比较完善的区域案例。

A profile of international climate change science in the past one hundred years

国际气候变化科学百年研究态势分析

A review of climatic impacts on Chinese socio-economic development over the past 2000 years

过去2000年气候变化对中国经济与社会发展影响研究综述

Historical climate change impact-response processes under the framework of food security in China

粮食安全视角下中国历史气候变化影响与响应的过程与机理

hunger and mass migration

太阳、气候、饥荒与民族大迁移

The causality analysis of climate change and large-scale human crisis

Transmission of climate change impacts from temperature change to grain harvests, famines and peasant uprisings in the historical China

冷暖—丰歉—饥荒—农民起义: 基于粮食安全的历史气候变化影响在中国社会系统中的传递

Climate and famines: A historical reassessment

Numerical processing of documentary sources on natural disasters

自然灾害史料的信息化处理

Study on the statistical law of historical famine in Gansu, Ningxia and Qinghai

甘宁青历史饥荒统计规律研究

Study on the statistical law of historical famine in Shaanxi

陕西历史饥荒统计规律研究

Study on the relationship between environmental change and famine occurrence in the middle reaches of the Yellow River in ancient times

古代黄河中游饥荒与环境变化关系及其影响

AD)

中国西汉—清代饥荒序列的重建及特征分析

Transmission pathways of China's historical climate change impacts based on a food security framework

Variation of precipitation for the last 300 years over the middle and lower reaches of the Yellow River

黄河中下游地区过去300年降水变化

Reconstruction of poor/bumper autumn harvest index series based on historical chorography and its relationship with climate change in North China from 1736 to 1911

基于地方志的1736—1911年华北秋粮丰歉指数序列重建及其与气候变化的关系

A study on the droughts and their impacts from 1470 to 1979 in North China

近500年华北地区最严重的干旱及其影响

1876-1878 severe drought in North China: Facts, impacts and climatic background

1876—1878年华北大旱: 史实、影响及气候背景

Famine, migration and war: Comparison of climate change impacts and social responses in North China between the late Ming and late Qing dynasties

Cannibalism in northern China between 1470 and 1911

Methodology for quantifying social and economic series on the base of semantic differential to serve the research on impacts of climate change

历史气候变化影响研究中的社会经济等级序列重建方法探讨

Reconstruction and analysis on the series of winter-half-year temperature change during the Qing Dynasty in the northern China region

清代华北地区冬半年温度变化重建与分析

Spatial-temporal distribution of locust plague and its relationship with flood /drought in North China during the Qing Dynasty

清代华北蝗灾时空分布及其与水旱灾害的关系

Comparison of governmental relief food scheduling and social consequences during droughts of 1743-1744 AD and 1876-1878 AD in North China

华北1743—1744年与1876—1878年旱灾中的政府粮食调度与社会后果对比

THE END
0.第十六届水科学发展论坛中国科学院南京地理与湖泊研究所李相虎副研究员作报告 三峡大学彭涛副教授作《变化环境下汉江流域非平稳水文干旱演变特征及归因分析》主题报告。在报告中,彭涛副教授首先说明研究背景及研究数据;然后,从水文干旱指数、非平稳水文干旱指数、干旱识别、统计方法和归因分析五方面介绍了研究方法;在此基础上,从非平稳标准化径jvzquC41yy}60ƒw0gjv0ls1yu5jpot132791:>450nuo
1.年际SPEI变化特征文献摘录spei热点图文章详细分析了辽宁省1970-2020年不同时间尺度下的SPEI值,揭示了月、季、年尺度的干旱特征及其驱动因素。月尺度波动剧烈,反映了降水对干旱的影响;季尺度和年尺度则显示出季节性和长期趋势。文章还讨论了干旱的周期性变化以及未来气候趋势,强调了SPEI在干旱预警和农业发展中的作用。 jvzquC41dnuh0lxfp0tfv8|gkzooa=947;;348ftvkimg8igvcomu865:6=44@:
2.中国近代北方极端干湿事件的演变规律[陈文海, 柳艳香, 马柱国. 中国1951-1997年气候变化趋势的季节特征. 高原气象, 2002, 21(3): 251-257.] [10] Ma Zhuguo, Wei Helin, Fu Congbin. Relationship between regional soil moisture variation and climatic variability over East China. Acta Mete. Sinica, 2000, 58(3): 278-287. jvzquC41yy}/inti0eun0ls1EP522763:470zk7225=t29=
3.中国水文地质及水循环中国不同区域的水文地质特征差异显著。华北地区由于长期过量开采地下水,导致地下水位持续下降,形成了大范围的地下水漏斗区;长江流域则拥有丰富的地表水和地下水资源,但近年来受人类活动和气候变化影响,水资源分布不均和污染问题日益突出;西北地区干旱少雨,水资源匮乏,地下水成为当地生活和农业用水的主要来源;青藏高原地区jvzquC41yy}/5?5fqe4dp8ftvkimg8;98:83a:6443=55A;0jvsm
4.1955—2021年黄土高原地区相对湿度时空演变规律黄土高原是我国气候变化敏感区,研究该区相对湿度(RH)时空演变及其与地理因子、气象要素之间的关系,有利于全面了解黄土高原气候变化规律。基于黄土高原及其周边地区90个气象站点逐日RH观测资料,利用趋势分析、敏感性分析和空间插值等方法,分析了1960—2021年及生态工程实施前后RH变化特征。结果表明:(1) 时间上,黄土高原地区jvzquC41yy}/inttgu4dqv3ep1mismq1EP52276433>0l7nuup42295/82<10;5450717
5.江苏省土地利用功能变化及其空间格局特征上述一系列城镇发展、农业生产、生态环境问题,归根结底是不同土地利用方式直接或间接影响国土空间发展而形成的,其结果是导致提供产品和服务满足人类需求的能力发生变化,从而影响区域可持续发展目标的实现。研究江苏省土地利用功能及其变化特征,确定和衡量多元化土地利用提供产品和服务的能力及其变化规律,有助于准确把握土地jvzq<84uejupn7ktggqbq‚fp0eun1ko1kiyot{4424703;24;17769<:5:963>5923=/uqyon
6.中国森林生态系统净初级生产力时空分布及其对气候变化的响应研究综合国内外目前的NPP研究, 研究内容主要包括四个方面:1)在时间序列上的波动情况;2)空间分布特征, 包括经纬向变化规律, 区域之间以及气候带之间的变化与差异;3)驱动因子;4)估算方法。由于森林生态系统本身的复杂性、野外测定困难, 同时也受人类活动影响等原因, 目前国际上对森林生态系统NPP的时空格局分异及生物地理jvzquC41yy}/gltnqiodc7hp1uzyd8hj1jznn872425268xvzd813B572:6:4A3jvo
7.我国干旱气候年际和年代际变化特征年际和年代际变化特征 从干旱的成因来看,由于干旱的某些影响因素,如太阳黑子的变化、日月食的出现和厄尔尼诺现象等都具有一定的周期性,因此干旱的发生也具有一定的周期性。又因为干旱的发生受到多种因素的影响,因而干旱发生的周期不是单一周期,而是复杂的混合周期。重旱面积百分率为达到重旱等级以上的面积占除常年干旱区jvzquC41yy}/q|lgq0io1ytuv1>1f@h