도시화가 출산에 미치는 영향에 관한 연구

Abstract

This study aims to empirically examine the effects of urbanization on fertility and to investigate the mediating roles of human development, gender inequality, and income inequality in this relationship. While a substantial body of prior research has confirmed a negative association between urbanization and fertility, the structural pathways through which urbanization influences fertility have not been sufficiently elucidated. Drawing on relevant theoretical frameworks—including r/K selection theory, the gender equity theory of fertility, and Easterlin's hypothesis—this study constructs a multiple mediation analytical framework and conducts quantitative analysis. For the analysis, a panel dataset covering 105 countries worldwide over a 34-year period from 1990 to 2023 was assembled. To address missing data, multiple imputation using the Amelia algorithm was applied. The primary analytical method employed was pooled OLS panel regression, which simultaneously utilizes between-country structural variation and within-country temporal variation to estimate long-term structural relationships among variables. To verify the robustness of the findings, a multiple mediation structural equation model （SEM） assuming the same causal structure was employed as a robustness check and supplementary analysis. Additionally, two-way （country- and year-） fixed effects （FE） panel regression results are presented in an appendix to facilitate comparison according to whether fixed effects are applied. In the main analysis—pooled OLS panel regression following Amelia multiple imputation—the key independent variable, the urbanization rate, was found to exert a statistically significant positive （+） effect on the Human Development Index （HDI） and a statistically significant negative （−） effect on the Gender Inequality Index （GII）, while its effect on the disposable income Gini coefficient was negative （−） in direction but did not reach statistical significance, remaining only at a marginal level. In the integrated equation incorporating the urbanization rate and all three mediating variables as predictors of the fertility rate, the direct effect of the urbanization rate was positive （+） in direction but statistically non-significant, while among the mediating variables, only the HDI showed a strong and statistically significant negative （−） effect on the fertility rate. This indicates that the fertility-reducing effect of urbanization operates primarily through the indirect pathway of improving human development levels rather than through the direct influence of urbanization itself. In the pathway-specific indirect effects analysis, only the indirect pathway via the HDI （urbanization → human development → fertility） was strongly significant in the negative （−） direction at the 0.1% significance level, while the indirect pathways via the GII and the disposable income Gini coefficient were not statistically significant. The total indirect effect, aggregating all three mediation pathways, was significant in the negative （−） direction at the 0.1% significance level, and the total effect of the urbanization rate on the fertility rate was likewise confirmed to be significant in the negative （−） direction at the 5% significance level. These results provide empirical support for the conclusion that, among the three competing hypotheses, the mediation pathway of CH1 （the ecological living environment hypothesis）—positing that urbanization raises human development levels, which in turn reduces fertility—possesses the strongest explanatory power in the main analysis. Among the control variables, both the working-age population ratio and the male unemployment rate were found to exert statistically significant negative （−） effects on the fertility rate. In the robustness check and supplementary analysis—a multiple mediation SEM following Amelia multiple imputation—the estimated results for the a-paths （urbanization → mediating variables） were largely consistent with those of the main pooled OLS analysis in terms of both direction and statistical significance. The direct effect of the urbanization rate was significant in the positive （+） direction, which is consistent in direction with the pooled OLS results, though differences in significance levels were observed due to the underestimation of standard errors attributable to the non-application of cluster-robust standard errors in the SEM. The key findings—including the significant negative （−） indirect pathway via the HDI, the non-significant pathways via gender inequality and income inequality, and the significant negative （−） total indirect effect and total effect—were consistently reconfirmed, thereby supporting the robustness of the main analysis results. This study holds academic significance in that it analyzes the relationship between urbanization and fertility from a multiple mediation perspective, and explicitly demonstrates how the choice of analytical methodology affects the interpretation of research findings by systematically comparing results from long-term cross-sectional pathway analysis incorporating between-country structural variation （pooled OLS）, simultaneous equation system-based robustness verification （SEM）, and within-country temporal variation-focused analysis （FE）. The study also contributes theoretically by empirically establishing that, among the three competing hypotheses, the human development pathway （CH1） constitutes the most robust mediating pathway explaining the structural relationship between urbanization and fertility. These findings suggest that, in designing policies to address the problem of low fertility, it is necessary to recognize the limitations of short-term policy instruments and to pursue long-term structural policy approaches responsive to the enhancement of human development levels.|본 연구는 도시화가 출산에 미치는 영향을 분석하고, 그 과정에서 인적개발 수준, 성불평등 수준, 소득불평등 수준이 매개변수로서 어떠한 역할을 하는지를 실증적으로 규명하는 것을 목적으로 한다. 기존의 연구들 에서는 도시화와 출산 간에 음（-）의 관계가 있는 것으로 다수 확인된 바 있으나, 도시화가 어떠한 구조적 경로를 통해 출산에 영향을 미치는지에 대해서는 충분한 실증적 규명이 이루어지지 못하였다. 이에 본 연구는 r/K 선택 이론, 성평등화 출산이론, Easterlin의 출산이론 등 관련 이론을 토대로 다중매개 분석틀을 설계하고 정량분석을 실시하였다. 분석을 위해 전 세계 105개국을 대상으로 1990년부터 2023년까지 34 년 간의 패널데이터를 구축하였으며, 결측치 문제를 보완하기 위해 Amelia 알고리즘을 활용한 다중대치를 적용하였다. 주 분석방법으로는 국가 간 구조적 수준 차이와 국가 내 시간변화를 동시에 활용하여 변수 간 장기적·구조적 관계를 추정하는 pooled OLS 패널회귀분석을 사용하였고, 분석 결과의 강건성을 검증하기 위해 동일한 인과구조를 가정한 다중매개 구조방정식모형을 강건성 확인 및 보조분석으로 활용하였다. 아울러 국가· 연도의 이원 고정효과 패널회귀분석 결과도 부록으로 제시하여 고정효과 적용 여부에 따른 분석결과를 비교할 수 있도록 하였다. 주 분석인 ‘Amelia 다중대치 적용 후 pooled OLS 패널회귀’에서는 핵심 독립변수인 도시화율이 인적개발지수（HDI）에 양（+）의 유의미한 영향을, 성불평등지수에 음（-）의 유의미한 영향을 각각 미치는 것으로 나타났으나, 가처분소득 기준 지니계수에 대해서는 음（-）의 방향이나 경계적으로 유의하지 않았다. 도시화율과 세 매개변수를 출산율에 통합적으로 투입한 방정식에서, 도시화율의 직접효과는 양（+）의 방향이었으나 통계적으로 유의하지 않았으며, 매개변수 중에서는 인적개발지수만이 출산율에 음 （-）의 강한 유의성을 보였다. 이는 도시화의 출산율 저하 효과가 도시화 자체의 직접적 영향보다는 인적개발 수준 향상을 통한 간접효과를 통해 주로 작동함을 의미한다. 경로별 간접효과 분석에서는 인적개발지수를 통한 간접경로（도시화→인적개발→출산）만이 음（-）의 방향으로 0.1% 유의 수준에서 강하게 유의하였고, 성불평등지수 및 가처분소득 지니계수를 통 한 간접경로는 유의하지 않았다. 세 매개경로를 합산한 총 간접효과는 음 （-）의 방향으로 0.1% 유의수준에서 유의하였으며, 도시화율의 출산율에 대한 총효과 역시 음（-）의 방향으로 5% 유의수준에서 통계적 유의성이 확인되었다. 이는 세 경쟁가설 중 ‘도시화 → 인적개발 수준 향상 → 출산 저하’를 핵심 내용으로 하는 CH1（생태학 기반 삶의 환경 가설）의 매개경로가 주 분석에서 가장 강한 설명력을 갖는 경로임을 실증적으로 지지하 는 결과이다. 한편 통제변수인 생산연령인구비율과 남성의 실업률은 모두 출산율에 음（-）의 유의한 영향을 미치는 것으로 나타났다. 강건성 확인 및 보조분석으로 실시한 'Amelia 다중대치 적용 후 다중매 개 구조방정식모형'에서는, ‘도시화 → 매개변수’의 a경로의 추정결과가 주 분석인 pooled OLS와 방향성 및 유의성이 다수 일치하였다. 도시화율의 직접효과는 양（+）의 방향으로 유의하게 나타나, pooled OLS와 방향은 일 치하였으나 군집표준오차 미적용에 따른 표준오차 과소추정으로 인해 유의성 수준에는 차이가 있었다. 인적개발지수를 통한 간접경로의 음（-）의 유의성, 성불평등 및 소득불평등 경로의 비유의성, 총간접효과 및 총효과의 음（-）의 유의성 등 핵심 결과가 주 분석과 일관되게 재확인되어, 주 분석 결과의 강건성이 지지되었다. 본 연구는 도시화와 출산율 간의 관계를 다중매개 관점에서 분석하고, 국가 간 구조적 수준 차이를 포함한 장기적·횡단면적 경로 분석（pooled OLS）과 동시방정식 체계 기반의 강건성 확인（구조방정식모형） 및 동일국 가 내 시간변화를 중심으로 한 분석（이원 고정효과 패널회귀）의 결과를 체계적으로 비교함으로써, 분석방법론의 선택이 연구결과의 해석에 미치는 영향을 명시적으로 제시하였다는 점에서 학술적 의의를 갖는다. 또한 세 경쟁가설 중 인적개발 경로（CH1）가 도시화와 출산과의 구조적 관계를 설명하는 가장 강건한 매개경로임을 실증적으로 규명하였다는 점에서도 이론적 기여를 갖는다. 이러한 결과는 저출산 문제에 대응하기 위한 정책 설계에 있어 단기적 정책수단의 한계를 인식하고, 인적개발 수준 상향에 대응한 장기적·구조적 정책 접근의 필요성을 시사한다. 【주요어】 출산율, 도시화율, 인적개발지수（HDI）, 성불평등지수（GII）, 가처분소득 지니계수, 생산연령인구, 남성의 실업률, pooled OLS 패널회귀, 다중매개 구조방정식모델（SEM）, 이원 고정효과（FE） 패널회귀, 매개효과, Amelia 다중대치, 결측치