Based on the above analysis results, we theoretically analyze the evolution mechanism of regional population aging in natural watersheds caused by population migration, and further build a series of empirical models in case of the Yellow River Basin.
5.1. The Theoretical Framework of Migration’s Effect on Population Aging
Population aging is a dynamic gradual process of demographic age structure. Based on demographic transition [
9,
10], the direct reason for population aging is the increase of aging population or the decrease of non-aged population, namely: aging population grow faster than non-aged population. Population aging can be divided into various types (
Figure 5). In general, the cause of relative population aging is either the slower growth of children at the bottom of demographic age pyramid (decreased fertility rate), the loss of labor force at the middle of demographic age pyramid, or the acceleration of the growth of aging population at the top of demographic age pyramid (decrease in mortality), namely: the absolute scale of aging population has not changed much, while whose relative proportion has increased a lot. However, the cause of absolute population aging is decrease in aging population, including both local population aging and agglomeration of adventive aging population, which is an increase in the absolute proportion of aging population at the top of demographic age pyramid. At present, relative population aging and absolute population aging have generally been variously undergone.
With the current slowdown in natural growth, mechanical growth has become the key to regional demographic transition. Migrants can be divided into two categories and three subcategories. Among them: labor migration is mainly one-way flow from developing areas to developed areas, mainly affected by economic gap (especially interregional wage gap). However, elderly migration is characterized by two-way flow, which is mainly affected by welfare conditions (especially social security), without considering elderly migration [
7] from developing high welfare areas to developed low welfare areas, while which exists.
There are many reasons for differentiation in population aging. Among them, population migration, family planning policy and mortality are the direct causes, while socio-economic development and population inertia (or population base) are the fundamental factors, and social security such as household registration system, social welfare are auxiliary. Economic development, that is, economic base has largely advanced or delayed population aging. However, human history, that is, human’s long development process makes population accumulation have a great inertia. We roughly divide the factors influencing population aging into four categories, namely: proximity effect, social security, economic development, and population inertia. The influence mechanism of these factors on population aging varies. Among them, the impact of economic development on population aging can be summarized as Productive Aging [
48] or Active Aging [
49], while the degree of geographical equalization of social security has an important impact on population aging, especially successful aging [
50]. Therefore, the general evolution mechanism of population aging is roughly as follows.
We assume that there are four major regions on a homogeneous plain, namely developed high-welfare region (A), developing high-welfare region (B), developed low-welfare region (C), and developing low-welfare region (D). Population migration is free among the four regions, only affected by wages or benefits, especially labor migration and elderly migration. The migration direction of labor force is only guided by wages, while which in developed regions is relatively high. The direction of elderly migration is divided into two categories. Elderly migration that belongs to successful aging is concentrated in high-welfare regions, while eliminated this part of the elderly can only stay in low-welfare regions. This also reflects the inequality of the elderly. In the long-wave period, labor force eventually become the elderly and is the main reason for change of regional population aging. The migration direction of labor force is mainly affected by their original emigration area, that is, the part of the successful aging labor force immigrates in high-welfare regions, while whose surplus is stranded in low-welfare regions. Through the interaction of the above factors, labor migration and migratory elderly result in the formation of an elderly gathering place in A at the primary and intermediate stages of population aging with regional economic development and social security improvement. They also result in the formation of an elderly gathering place in B at the intermediate and advanced stages of population aging when degree of geographical equalization of social security continues to improve, while regional economic development no longer significantly affects social security (
Figure 6). Among four regions, labor force and the elderly continuously gather in A, while B only gathers the elderly and outputs labor force. This leads to population aging in B higher than that in A, while the scale of the elderly in B may not be as large as that in A. D not only losses many labor force, but also losses many elderly, so the degree of its population aging depends on the loss of the two. However, C attracts many labor force and losses part of the elderly, so its population aging is the lowest in the entire plain. Thus, the polarization effect of population aging in the entire plain is formed. Meanwhile, a relatively stable regional spatial pattern of population aging is formed. This is a qualitative quadri-area model for regional population aging.
5.2. The Effect Mechanism of Migration on Population Aging in the Yellow River Basin
Based on the above theoretical mechanism analysis and data related to the Yellow River Basin during 2000–2020 by SPSS25.0, we further verify the effect mechanism of migration on population aging in natural watershed. We use regression models to quantitatively determine the linkage effects between population aging and various dominant factors.
RDE and
RDI are not suitable for regression analysis because of only 9 samples and 10 samples during 2000–2020 (
Table 6), while the remaining four types of migration areas are suitable for regression analysis with enough samples.
Different types of population migration have various effects on regional population aging. For these four types, we use
AP or
PA as dependent variables, and the remaining 6 indicators as independent variables. According to
Table 2 and
Table 6, our empirical model is divided into four categories as follows:
where
Areaij ∈ [
SDE,
SSE,
SSI,
SDI],
j ∈ [2000, 2005, 2010, 2015, 2020],
i varies depending on each type,
i ∈ [1, …, 89]
SDE, [1, …, 99]
SSE, [1, …, 67]
SSI, [1, …, 91]
SDI.
First, we perform variable standardization and correlation analysis, and then perform stepwise regression analysis to obtain the regression equation. We find that
AP or
PA has a significant positive correlation with
VGDP,
Pu and
Pd, and a low negative correlation with
Ph (
Table 7). Except for
SSI,
AP or
PA in other types has a low negative correlation with
Mr.
According to the specific regression results (
Table 8), we find that
PA or
AP is mainly affected by
VGDP,
Pu and
Ph, but only
VGDP is positively affected in
SDE. Among them, for every percentage point increase in
VGDP,
PA will increase by 0.441% in
SDE with other factors unchanged, while
AP will do by 0.535%. However, it is
Ph that has the greatest impact on
PA or
AP, although it has a negative impact.
PA will decrease by 0.674% in
SDE with other factors unchanged, while
AP will do by 1.555% with every one percentage point increase in
Ph.
PA is mainly affected by
VGDP and
Pu, while
AP is done by
Pd and
Pu, both of which have no significant correlation with
Ms and Ph in
SDI.
PA will decrease by 0.417% in
SDI, while
AP will do by 0.569% with every one percentage point increase in
Pu with other factors unchanged. Meanwhile, for every percentage point increase in
Pd,
AP will increase by 0.481% in
SDI.
PA is mainly affected by
VGDP, while
AP is done by
Ms and
Pd, both of which have no significant correlation with
Pu and
Ph in
SSE. For every percentage point increase in
VGDP,
PA will increase by 0.218% in
SSE with other factors unchanged. Under the same conditions as other factors,
AP will increase by 0.440% in
SSE with every one percentage point increase in
Pd, while which will do by 0.412% with every one percentage point increase in
Ms.
PA only has a significant positive correlation with
VGDP and
Pd in
SSI, while
AP does especially with
Ms and
Mr. Under the same conditions as other factors,
AP will decrease by 1.124% in
SSI with every one percentage point increase in
Mr, while which will increase by 1.107% with every one percentage point increase in
Ms. In addition,
Ph also significantly affects
AP in
SSI.
In general, population aging in these four types of areas is not only affected by migration factors, but also obviously interfered by socio-economic factors. When considering the influence of a single factor, AP can be regarded as a power function of Mr or Ms. From this, we can judge the shape of AP’s function curve of Mr or Ms, and then judge the trend of linkage between them.
Based on the first-order partial derivative and the second-order partial derivative of these power functions about
Mr or
Ms, we can judge the extent of the impact of population migration on regional population aging, especially the scale of aging population. However, due to the little difference among their exponents, these power function curves are difficult to distinguish in common drawing software. Therefore, we make 2D function image drawing with Origin 2019. In the process, we only regard the core explanatory variable
Mr or
Ms as independent variables, while the remaining control variables are temporarily ignored. After getting the original image, we select the “signal processing” in “analysis” menu bar of Origin 2019 toolbar to perform “smoothing”, in order to process problems such as irregular and noisy data. According to the characteristics of these power function, we find that only the power function graph of
AP in
SSI about
Ms is convex, and that of other types of
AP about
Mr or
Ms is monotonically decreasing, while the inclination degree of whose graph differs (
Figure 7). These function curves have no inflection points because of
. Among them,
AP in
SDE and
SSE is a monotonically decreasing concave function about
Mr, and a monotonically increasing convex function about
Ms;
AP in
SDI is a monotonically decreasing concave function of
Mr;
AP in
SSI is a monotonically decreasing concave function about
Mr, and a monotonically increasing concave function about
Ms (
Figure 7).
In short, as a direct factor, population migration obviously differs in the impact of regional population aging, while the impact of indirect factors such as economic development is more consistent. This also explains why the existing research models mostly mix these factors together instead of separately expressing them. Based on the existing research paradigm, the complex relationship between these direct and indirect factors limits the precise analysis of evolutionary mechanism of regional population aging. This is also a common problem encountered in quantitative research in social sciences because of the complexity of social phenomena. This has always been a difficulty in model construction because of some common problems such as the correlation between independent variables. In construction of regression model, how to deal with the comparative analysis of the impact of direct and indirect factors on regional population aging has become an important aspect of our future efforts.