Abstract
The strategic value of the digital economy in driving the transformation of the tourism industry has become increasingly prominent. This study investigates the impact of China’s National Big Data Comprehensive Pilot Zones (NBDCPZs) on tourism economic development. Using panel data from 284 cities (2009–2019), the study treats the policy as a quasi-natural experiment and applies multi-period DID, PSM, and spatial models. Results show that NBDCPZs significantly promote tourism growth, though effects decline over time and vary spatially—stronger in eastern regions and moderate in large cities. Industrial structure upgrading is identified as the main transmission channel, while financial support amplifies the policy effect, and human capital shows no significant moderation. The findings reveal how digital policy drives tourism transformation and provide practical insights for strengthening data—tourism integration, optimizing industrial structures, and fostering high-quality regional development.
Keywords
Introduction
As a strategic pillar of China’s national economy, the tourism industry is experiencing a fundamental shift from a traditional resource-dependent model to a data-driven paradigm. The convening of the 2024 National Tourism Development Conference—China’s first central-level meeting devoted exclusively to tourism—marked a turning point in the nation’s development agenda. The conference emphasized the need to “promote high-quality tourism development and accelerate the building of a strong tourism nation,” signifying that advancing high-quality growth in the tourism economy has become both a national priority and a key component of China’s modernization process. High-quality tourism development entails improving efficiency, optimizing industrial and consumption structures, and enhancing environmental sustainability (Liu & Han, 2020). Scholars have approached this issue from multiple perspectives, including the theoretical connotations of tourism economic development (Li & Wei, 2014), the spatial correlation structure and its effects (Wang et al., 2017), and the determinants influencing tourism growth (Zhou et al., 2021). At its core lies a fundamental question: how can data-driven innovation be harnessed to achieve sustainable tourism economic growth?
Amid this transformation, the digital economy has emerged as a decisive force reshaping the production logic and growth trajectory of the tourism industry. By improving productivity, optimizing resource allocation, and promoting structural upgrading, digitalization enables tourism enterprises and destinations to achieve higher levels of operational efficiency and service innovation (Chen et al., 2025). In this context, big data policies play a crucial role by facilitating the market-oriented allocation of data elements, deepening industrial integration, and enhancing policy precision (Wei et al., 2022). Understanding how these policies influence tourism economic development and through which mechanisms they operate is therefore of significant theoretical and practical importance for achieving data-enabled improvements in tourism quality and efficiency. Since 2015, China has successively established several National Big Data Comprehensive Pilot Zones (NBDCPZs)—including the Guizhou pilot zone, two interregional zones (Beijing-Tianjin-Hebei and the Pearl River Delta), four demonstration zones (Shanghai, Henan, Chongqing, and Shenyang), and one infrastructure-oriented zone (Inner Mongolia). These zones serve as experimental platforms for promoting digital transformation, integrating data resources, and advancing cross-sector collaboration. From a broader industrial policy perspective, National Big Data Comprehensive Pilot Zones (NBDCPZs) represent a quintessential top-down digital governance tool within China’s evolving industrial strategy framework. Existing research indicates that large-scale digital and industrial policies—such as Made in China 2025—play a pivotal role in advancing sustainable development by reshaping production efficiency, innovation structures, and green transition pathways (Zhang et al., 2025). Emerging evidence further demonstrates that digitally oriented trade policies foster green growth by reducing transaction costs, enhancing cross-border service efficiency, and accelerating industrial upgrading—particularly in service-intensive sectors (Zhang & Choi, 2025). Compared to manufacturing-focused strategies, NBDCPZs prioritize data factor allocation, digital infrastructure integration, and platform-based governance—characteristics highly aligned with the service-dominant nature of the tourism economy. Therefore, analyzing digital trade pilot free trade zones within the broader framework of China’s digital industrial policy system helps clarify their structural role in promoting high-quality, sustainable tourism development.
Given the lack of standardized indicators to comprehensively measure big data development, the establishment of NBDCPZs provides an ideal quasi-natural experiment for evaluating the economic impact of data policies. Drawing on panel data from 284 prefecture-level cities in China between 2009 and 2021, this study investigates the causal effects and mechanisms through which NBDCPZs influence tourism economic development. It contributes to the literature in three main ways. First, while existing studies have examined the effects of big data policies on innovation, industrial upgrading, and macroeconomic growth few have explored their impact on tourism (Afolabi, 2023; Chen et al., 2022; Zhao et al., 2024). Second, prior studies have investigated the formation and development of the tourism economy (Wei & Hu, 2025), but most remain at a qualitative level and have paid limited attention to the causal relationship between big data policy and tourism economic performance by empirically testing the relationship between digital policy interventions and tourism performance, this study clarifies the role of data as a production factor in the tourism economy. Third, whereas prior research on NBDCPZs has focused primarily on firm-level or regional outcomes (Yue et al., 2024), this study highlights their effects on local tourism economies and spatial heterogeneity, refining the theoretical understanding of how big data initiatives can drive tourism transformation in the era of digital intelligence. This study analyzes the innovation-driven logic of data elements in shaping tourism economic development, examines the policy effects and transmission mechanisms of the pilot zones, and thereby contributes to improving the theoretical framework of coordinated development between digital technology and the tourism economy. The findings also offer practical implications for optimizing regional big data policy design and advancing the digital transformation of the tourism industry.
Literature Review
Theoretical Analysis and Research Hypotheses
As a key platform for advancing the market-oriented allocation of data elements and optimizing digital infrastructure, the National Big Data Comprehensive Pilot Zones (NBDCPZs) have injected new momentum into China’s tourism transformation. In the context of the digital economy, tourism production modes and industrial boundaries are being reshaped, profoundly influencing the future trajectory of tourism development (Wei, 2022). Empirical evidence indicates that smart tourism models significantly enhance sustainable tourism growth, improving efficiency and innovation through digital technologies (Phaosathianphan & Leelasantitham, 2021). From an international trade perspective, tourism is widely regarded as an intangible export whose competitiveness is highly dependent on information and communication technology (ICT). A growing body of trade literature indicates that ICT is fundamentally reshaping international trade patterns by reducing information asymmetry and facilitating cross-border service transactions (Park & Choi, 2025). Digital tools such as online booking platforms, intelligent recommendation systems, and data-driven marketing channels enable tourism destinations to expand market reach and enhance service matching efficiency. Further evidence indicates that ICT development significantly enhances export performance, particularly in service-oriented and information-intensive industries (Nguyen & Choi, 2025). By strengthening digital infrastructure, data integration, and platform governance, New Frontier Development Cooperation Plan Zones (NBDCPZs) have elevated tourism destinations’ capacity to participate in digital services trade. Therefore, this study proposes the following hypothesis:
However, regional heterogeneity shapes the effectiveness of big data policies. Differences in resource endowments, economic foundations, and policy environments have created marked spatial disparities in China’s digital development (Lü et al., 2023; Song et al., 2020). Regions with weaker initial conditions often experience stronger path dependence, limiting policy effectiveness (Jedwab et al., 2017). Eastern China—characterized by advanced digital infrastructure, higher marketization, and stronger policy adaptability—translates data resources into tourism value more efficiently than central and western regions, which face deficits in infrastructure, human capital, and financial input (Zhi et al., 2021). Consequently, the policy effects of NBDCPZs are expected to be strongest in eastern cities. Therefore, this study proposes the following hypothesis:
The NBDCPZs also promote tourism growth through industrial restructuring. Industrial structure optimization is a key determinant of high-quality regional development. Big data policies accelerate the transformation of traditional industries toward high value-added, knowledge-intensive sectors, fostering integration between tourism, culture, and technology (Guo et al., 2024). Empirical research indicates that tourism-related foreign direct investment and green technological innovation jointly promote sustainable tourism development by fostering technology spillover effects and enhancing environmental performance (Zhou & Choi, 2025). Furthermore, by deepening the integration of manufacturing and modern service industries, pilot regions have established digital service systems that support sustainable tourism upgrades, thereby strengthening their competitiveness in attracting high-end tourism projects (Chen & Zhang, 2023). The resulting “structural dividend” expands the tertiary sector and enhances environmental and allocative efficiency, improving both the quality and spillover effects of tourism growth (Wang et al., 2023).Therefore, this study proposes the following hypothesis:
Information infrastructure plays a pivotal role in transmitting benefits. Research indicates that digitalization reduces information barriers and expands digital service coverage, thereby mitigating socioeconomic inequality (Yin & Choi, 2023). Simultaneously, data-driven systems significantly enhance audience engagement by enabling personalized interactions, real-time feedback, and behavioral response mechanisms. Empirical research indicates that digital platforms substantially boost audience participation and behavioral responses through personalized interaction mechanisms (Yun et al., 2025). Innovative initiatives such as smart scenic area systems, crowd flow monitoring, and location-based services vividly demonstrate how informatization enhances tourism industry performance (Gannon & Tahri, 2020). Collectively, these mechanisms foster the emergence of “new productive forces,” characterized by the organic integration of technological innovation, institutional efficiency, and inclusive participation. Therefore, this study proposes the following hypothesis:
Human capital is often regarded as a key driver of innovation, sustainable trade, and green economic growth. Existing research indicates that human capital plays a crucial role in promoting sustainable trade and facilitating the absorption of digital technologies (Choi et al., 2023). However, the effectiveness of human capital depends on its alignment with industry-specific skill requirements. Within the tourism sector, digital transformation increasingly relies on interdisciplinary competencies integrating data analytics, platform operations, creative content production, and service innovation. Furthermore, cultural and creative innovation has been demonstrated to enhance emission reductions and environmental performance in tourism-related activities (Park et al., 2024). The potential mismatch between educational structures and digital tourism demands provides a theoretical explanation for the heterogeneous or weak moderating effects of human capital observed in policy impact assessments. Therefore, this study proposes the following hypothesis:
Finally, financial support amplifies the policy effect of NBDCPZs. Digital inclusive finance and financial technology alleviate financing constraints, promote investment in smart tourism, and facilitate digital transformation. Robust financial systems enhance policy transmission, industrial integration, and consumer responsiveness, jointly driving tourism upgrading. Empirical findings confirm that regions with higher financial support experience stronger NBDCPZ policy effects, with notable spatial spillovers (Guo & Ming, 2025; Xu et al., 2024).Therefore, this study proposes the following hypothesis:
Research Design
Research Area
This study selects 284 prefecture-level and provincial capital cities across China as the research sample, ensuring representativeness with respect to the locations of the National Big Data Comprehensive Pilot Zones (NBDCPZs). The selected cities cover a wide geographical range, including economically developed areas such as the eastern coastal and southern riverside regions, as well as the central, western, and northeastern parts of China. In terms of urban scale, the sample encompasses mega and large cities in addition to medium- and small-sized cities, thereby providing a comprehensive reflection of the characteristics of cities across different regions and administrative levels. The study period spans from 2009 to 2021, covering the establishment cycles of multiple batches of pilot zones. This time frame provides a complete temporal series for evaluating the dynamic effects of the policy implementation. The post-2019 period is deliberately excluded from the baseline analysis. The outbreak of COVID-19 constitutes a nationwide exogenous shock that led to an unprecedented collapse of tourism demand and supply, generating a structural break that is orthogonal to the implementation of the NBDCPZ policy. Including the pandemic years would contaminate the identification of long-term policy effects with short-term crisis-induced fluctuations. Restricting the sample to the pre-pandemic period is therefore necessary to ensure internal validity.
Research Data
Data were sourced from the China Urban Statistical Yearbook, provincial statistical yearbooks, and the EPS database. A tourism economic indicator system was constructed following principles of systematicity, scientific rigor, representativeness, and data availability, using sub-indicators proposed by Huang et al. (2024) suitable for the Chinese context (Table 1). Drawing on previous research (Gao & Li, 2024; Zhang & Guo, 2023), the core explanatory variable is a dummy for national-level big data pilot zones (0/1), with designated cities as the treatment group and others as the control group. To reduce bias in estimating tourism economic development, several control variables were included: per capita GDP (log), foreign investment ratio, government expenditure ratio, urban education spending, social consumption ratio, and scientific expenditure ratio (Table 2), capturing key economic, social, and technological factors influencing tourism development.
Index System for Tourism Economy.
Variable Definitions and Indicator Descriptions (Including Fixed Effects).
Data Processing Methods
Tourism economic development is measured using multiple sub-indicators reflecting tourism scale and urban development, normalized and weighted via the entropy method to ensure comparability (Liang et al., 2024). Entropy assigns higher weights to indicators with greater variation (Zhi-Hong et al., 2006). Ratio-type controls are log-transformed with a constant added to avoid undefined values. After merging into a panel dataset, missing values are interpolated, and multicollinearity is checked to ensure robust model estimation.
Research Methods
This study evaluates the impact of National Big Data Comprehensive Pilot Zones (NBDCPZs) on tourism economic development using a multi-level framework. Baseline regressions are conducted with mixed, fixed, and random effects models, selecting the optimal specification via the Hausman test. Time-lag models assess immediate, cumulative, and long-term policy effects. Propensity score matching addresses potential endogeneity, followed by regression on matched samples to ensure robustness. Spatial econometric models, including SLM and SDM, examine spatial spillovers, while heterogeneity and spatial quantile regressions capture variations across regions, city sizes, and tourism development levels. Mechanism analysis explores mediation pathways, and interaction terms test moderating effects of education and finance. All baseline, dynamic, spatial, heterogeneity, and mechanism analyses are consistently conducted using the pre-pandemic period (2009–2019) to ensure identification consistency and internal validity. This multi-method design comprehensively examines policy impacts from dynamic, spatial, distributional, and mechanism perspectives.
Benchmark Regression Model Specification
Benchmark regression is conducted using Pooled OLS, Fixed Effects (FE), and Random Effects (RE) models, with the Hausman test guiding model selection. The FE model, controlling for time-invariant city characteristics and common time effects, is preferred if the Hausman test is significant (Baltagi & Liu, 2016). Considering that the COVID-19 pandemic represents a nationwide structural shock that fundamentally altered tourism demand and supply conditions, the post-2019 period is excluded from the baseline analysis rather than being controlled for using a simple dummy variable. Robust standard errors are applied to address heteroskedasticity, and the regression equation is specified to estimate the average effect of big data pilot zones on tourism economic development.
Parallel Trend Test and Event-Study Design
The validity of the multi-period DID approach relies on the parallel trend assumption, which requires that treated and control cities exhibit similar tourism development trajectories prior to policy implementation. To formally test this assumption, an event-study specification is employed. Specifically, the treatment indicator is interacted with a series of lead and lag dummy variables constructed based on the relative time distance between the observation year and the policy implementation year. The year immediately preceding policy implementation is set as the reference period and omitted from the regression. If the coefficients of the pre-treatment periods (e.g., t − 3 and t − 2) are statistically insignificant and fluctuate around zero, the parallel trend assumption can be reasonably satisfied, thereby supporting the causal interpretation of subsequent dynamic treatment effects.
Time-Lag Effect Model Specification
A distributed lag model is used to assess the dynamic effects of National Big Data Comprehensive Pilot Zones (NBDCPZs) on tourism economic development. Lagged terms capture immediate (lag1), cumulative (sum of lags), and long-term effects. The ARDL framework distinguishes short-term impacts within the same period from long-term equilibrium effects, enabling a comprehensive analysis of the policy’s sustained influence across different time horizons (Cho et al., 2023). Within the ARDL framework, short-term and long-term effects are distinguished: short-term effects capture the impact of changes in explanatory variables on the dependent variable within the same period, whereas long-term effects reflect the equilibrium relationship over an extended horizon (Brauer, 2013). Accordingly, this study calculates the immediate effect (lag1), cumulative effect (lag2), and long-term effect (lag3) to reveal the sustained influence of NBDCPZ policies on tourism economic development.
Endogeneity Treatment: Propensity Score Matching (PSM)
Since the designation of National Big Data Comprehensive Pilot Zones (NBDCPZs) is non-random and may be correlated with pre-existing city-level economic and institutional characteristics, baseline DID estimates could be subject to selection bias. To mitigate bias arising from observable factors, this study applies propensity score matching (PSM), using a logit model to estimate each city’s probability of being designated as a pilot zone based on pre-treatment covariates. One-to-one nearest neighbor matching with a 0.2 caliper ensures covariate balance, assessed via standardized mean differences, t-tests, and bias reduction rates. After matching, two-way fixed effects panel regressions are conducted on the matched sample to re-estimate the policy effect. Importantly, PSM-DID is treated as a robustness check rather than the primary identification strategy, which relies on the parallel trend assumption validated through event-study analysis. Besides, the matched sample is subsequently used for all baseline estimations and robustness checks.
Placebo Test and Anticipatory Effect
Although PSM-DID mitigates selection bias from observable characteristics, unobserved time-varying factors may still bias the estimates. To further strengthen causal identification, this study conducts placebo tests by assigning fictitious policy implementation years that precede the actual policy adoption. If the estimated effects are driven by spurious correlations or pre-existing trends, significant placebo effects would be observed. Conversely, statistically insignificant placebo coefficients would provide supportive evidence that the baseline results are not driven by anticipation or reverse causality.
Spatial Econometric Model Specification
Spatial autocorrelation of tourism economic indicators is tested using Moran’s I. If significant, a spatial weight matrix is constructed, and Spatial Lag (SLM; Lam & Souza, 2020) and Spatial Durbin (SDM) models (Wang et al., 2021) are estimated. SLM captures spatial dependence of the dependent variable, while SDM accounts for both explanatory variables and their spatial lags.
Where W denotes the spatial weight matrix, represents the spatial lag coefficient for policy variables, ρ denotes the spatial lag coefficient, and θ signifies the spatial lag coefficient for independent variables. By comparing the estimation results of SLM and SDM, we can test whether policy effects are spatially robust and identify the channels of spatial spillovers. During estimation, the aforementioned conventional control variables are similarly controlled, and robust standard errors are employed to address issues such as spatial heteroskedasticity.
Heterogeneity Analysis
Heterogeneity tests examine how policy effects vary by region, city size, and tourism development level. Following Yu et al. (2023), Cities are grouped by geographic region and population (mega to small), and region- or size-specific regressions are estimated. Spatial quantile regression, incorporating a spatial weight matrix, assesses differential policy impacts across the conditional distribution of tourism economic development, capturing effects on cities with low, medium, and high tourism revenues.
In the formula, τ denotes the τth quantile. By conducting the aforementioned heterogeneity test, regional, scale, and level differences in policy effects can be revealed, providing reference for policy optimization.
Mechanism Testing Model Specification
To examine how big data pilot zones influence tourism economic development, mediation analysis is conducted using industrial structure (tertiary-to-secondary output ratio) and digital infrastructure (postal and telecom service volume) as mediators. Stepwise regressions assess whether policy effects operate through these channels, with the Sobel test evaluating significance. A significant policy–mediator and mediator–tourism relationship indicates partial or full mediation, while insignificance suggests the path is not primary. To test moderation, interaction terms between the policy dummy and educational or financial variables—such as higher education enrollment and the ratio of financial deposits and loans to GDP—are included. Significant interaction coefficients reveal how regional education and financial capacity strengthen or weaken the policy’s impact on tourism, capturing variations in policy effectiveness across different institutional and economic environments.
Empirical Findings and Analysis
Benchmark Processing Effects of National Big Data Comprehensive Pilot Zones
Baseline regressions using POLS, FE, and RE models (Table 3) show the Hausman test favors RE, with robust standard errors applied. Column (3) of Table 3 indicates NBDCPZs significantly boost tourism development (coef. = .0074, 10%). GDP per capita and government expenditure also positively affect tourism, supporting H1. Although the Hausman test suggests the random-effects estimator is efficient in the baseline setting, fixed-effects specifications are consistently adopted in subsequent analyses to ensure conservative identification and comparability across robustness, spatial, and heterogeneity exercises.
Results of the Benchmark Regression Analysis.
Note. Values in parentheses are z-statistics.
, **, and * indicates significance at the 1%, 5%, and 10% levels, respectively.
Parallel Trend Test and Event-Study Results
Figure 1 reports the event-study estimates examining the dynamic effects of the National Big Data Comprehensive Pilot Zones on tourism economic development. The coefficients for the pre-treatment periods (t − 3 and t − 2) are statistically insignificant and fluctuate closely around zero, indicating no systematic differences in tourism development trends between treated and control cities prior to policy implementation. This result supports the validity of the parallel trend assumption underlying the multi-period DID framework. The absence of significant pre-treatment coefficients indicates no detectable anticipation effects.

Parallel trend test results.
Following policy implementation, the estimated coefficients exhibit short-term fluctuations, with a modest negative effect observed in the contemporaneous period (t = 0), and heterogeneous responses in subsequent years. These patterns suggest that the policy impact does not materialize instantaneously and may involve adjustment costs and delayed responses.
Dynamic Effects of the Policy
A distributed lag model assesses the dynamic effects of NBDCPZs on tourism economic development. Results show the first-lag coefficient is .0060 (5%), indicating a significant 1.28% increase in the first year, while the second-lag coefficient is .0062 (5%), suggesting a persistent but smaller effect. Notably, the third-lag coefficient further increases to .0178 and is statistically significant at the 1% level, indicating a cumulative and strengthening policy effect over time. Among control variables, GDP per capita exhibits a strong and positive effect, together with government expenditure and retail sales, highlighting the importance of overall economic scale and fiscal capacity in supporting tourism development. Education expenditure shows a weakly negative effect at the 10% significance level, while foreign direct investment and science and technology expenditure remain statistically insignificant. The model’s R2 = .2421, indicating that the model’s R2 equals .2421, indicating that the distributed lag specification explains a substantial proportion of the variation in tourism economic development (Table 4).
Analysis Results of the Time-Effect Model.
Note. Values in parentheses are z-statistics.
, **, and * indicates significance at the 1%, 5%, and 10% levels, respectively.
Endogeneity Correction: PSM Estimation Results
Cities designated as NBDCPZs differ from non-designated cities in per capita GDP and financial development, potentially biasing estimates. Following the approaches of Zhi et al. (2021), and Shi et al. (2018), propensity score matching (PSM) is applied to mitigate selection bias. Table 5 shows baseline disparities are largely eliminated post-matching, with remaining standardized biases below 15%, validating the matched sample (Caliendo & Kopeinig, 2008). Regressions on matched samples (Table 6) show the NBDCPZ remains positive and statistically significant, though its magnitude slightly decreases from .0156 to .0133 after matching, confirming the robustness of the positive policy effect on tourism economic development. GDP per capita continues to exhibit a significant positive association, underscoring the persistent role of economic fundamentals in shaping urban tourism growth.
PSM Balance Test Results.
Results of PSM-Based Regression Analysis.
Note. Values in parentheses are z-statistics.
, **, and * indicates significance at the 1%, 5%, and 10% levels, respectively.
Placebo Test Results
To further assess whether the estimated effects are driven by spurious correlations or anticipation behavior, a placebo test is conducted by assigning a fictitious policy implementation year 2 years prior to the actual designation of National Big Data Comprehensive Pilot Zones. The regression results indicate that the placebo treatment coefficient is statistically insignificant and close to zero (Table 7), suggesting that the baseline estimates are unlikely to be driven by pre-existing trends or reverse causality.
Results of PSM-Based Regression Analysis.
Note. Values in parentheses are z-statistics.
, **, and * indicates significance at the 1%, 5%, and 10% levels, respectively.
Verification of Spatial Spillover Effects
Moran’s I statistics (.2077, significant at the 1% level) confirm the presence of significant spatial autocorrelation in tourism economic development, justifying the use of spatial econometric models. As shown in Table 8, the Spatial Lag Model (SLM, Column 8) indicates that the designation of National Big Data Comprehensive Pilot Zones (NBDCPZs) exerts a positive and statistically significant effect on local tourism economic development (coef. = .0113, 5%). Several control variables, including fiscal expenditure intensity and education expenditure, also exhibit marginal significance. The spatial lag term of the dependent variable (WY) is negative and significant, suggesting competitive or substitution effects among neighboring cities in tourism development.
Spatial Regression Model Analysis Results.
Note. Values in parentheses are z-statistics.
, **, and * indicates significance at the 1%, 5%, and 10% levels, respectively.
The Spatial Durbin Model (SDM, Column 9) further corroborates the positive policy effect of NBDCPZs, with a larger estimated coefficient (.0273, significant at the 10% level), indicating that accounting for spatially lagged explanatory variables amplifies the estimated local impact. In addition, the spatial lag of GDP per capita (WGDP) is significantly negative, implying that economic development in neighboring cities may crowd out local tourism growth. Although most spatially lagged covariates are statistically insignificant, the significant spatial autoregressive term (WY) highlights the existence of spatial dependence in tourism outcomes. Comparatively, the SDM achieves a higher explanatory power (R2 = .2985 vs. .2561 in the SLM), suggesting that explicitly modeling spatial spillover channels provides a better representation of regional interdependencies in tourism economic development.
Heterogeneity Analysis of Policy Effects
Grouped regressions reveal pronounced spatial heterogeneity in the effects of NBDCPZs (Table 9). The policy effect is significantly positive in Northern China (coef. = .0243, t = 4.0731) and remains robust in Southern China (coef. = .0277, t = 2.3324), indicating that regions with stronger economic foundations and digital adoption capacity benefit more from the policy intervention. By contrast, the estimated effects in Central and Northwestern China are statistically insignificant, suggesting limited short-term tourism responses in these regions. Regarding control variables, GDP per capita exerts a consistently positive and significant effect in Central, Northern, and Southern regions, underscoring the fundamental role of economic development in shaping regional tourism performance. In Southern and Northwestern China, government budget expenditure shows a significant promoting effect, while education and science-related expenditures exhibit heterogeneous impacts across regions, reflecting regional differences in factor endowments and policy transmission pathways.
Analysis Results of Urban Samples Across Different Regions.
Note. Values in parentheses are z-statistics. “–” indicates that the corresponding variable is excluded to mitigate multicollinearity in regressions.
, **, and * indicates significance at the 1%, 5%, and 10% levels, respectively.
When grouped by urban resident population size (Table 10), the results reveal clear heterogeneity in policy effectiveness across city scales. In mega and extra-large cities (Columns 16 and 17), the estimated coefficients of the NBDCPZ policy are negative but statistically insignificant, suggesting that the tourism-enhancing effect of the policy is limited in highly developed metropolitan areas. By contrast, large cities (Column 18) exhibit a significantly positive policy effect (coefficient = .0171, t = 2.3447), indicating that cities at this scale are more responsive to big data policy interventions.
Analysis Results for Samples at Different City Levels.
Note. Values in parentheses are z-statistics. “–” indicates that the corresponding variable is excluded to mitigate multicollinearity in regressions.
, **, and * indicate significance at the 1%, 5%, and 10% levels, respectively.
Regarding control variables, GDP per capita is significantly positive in mega and extra-large cities, reflecting the dominant role of overall economic development in shaping tourism outcomes at higher urban tiers. Fiscal expenditure intensity remains positively significant in extra-large and large cities, while science and technology expenditure shows heterogeneous effects across city groups, highlighting scale-dependent differences in policy transmission mechanisms.
Quantile regression results reveal pronounced nonlinear effects of NBDCPZs across the distribution of tourism economic development. The policy effect is insignificant at the lower quantiles of .1 and .3 and at the upper quantile of 0.9, while it becomes positive and statistically significant at the median and higher middle quantiles, reaching its maximum at the .7 quantile with a coefficient of .0079 and a t-value of 2.7613. This pattern indicates that NBDCPZs exert the strongest promotional effect on cities with relatively high levels of tourism development, rather than on cities at the lower or extreme upper ends of the distribution. The spatial lag term WY remains highly significant across all quantiles and increases monotonically, indicating that spatial spillover effects strengthen as tourism development intensifies. In addition, science and technology expenditure exhibits evident threshold characteristics, remaining negative and statistically significant at the middle quantiles of .5 and .7, but turning strongly positive at the .9 quantile, which highlights the differentiated role of technological investment across stages of tourism development (Table 11).
Results of Quantile Regression Analysis.
Note. Values in parentheses are z-statistics.
, **, and * indicates significance at the 1%, 5%, and 10% levels, respectively.
Testing the Transmission Mechanisms of Policy Effects
Mediation analysis shows NBDCPZs partially affect tourism via industrial structure, significantly increasing the tertiary-to-secondary sector ratio (coef. = .1230, t = 2.8577) with the Sobel test confirming mediation (Z = 2.4091). When included jointly, both policy and industrial structure remain significant. In contrast, postal and telecommunications volume shows no mediating effect, with insignificant coefficients (coef. = .0287) and Sobel Z = .6820, indicating this channel does not transmit policy impacts (Table 12).
Results of Mediating Effect Tests.
Note. Values in parentheses are z-statistics.
, **, and * indicates significance at the 1%, 5%, and 10% levels, respectively.
Moderation analysis reveals heterogeneous conditioning effects of human capital and financial development on NBDCPZ outcomes. As shown in Table 13, the interaction between NBDCPZs and higher education enrollment is positive and statistically significant (coef. = .3495, t = 2.1922), indicating that stronger human capital endowment enhances the tourism effects of the policy. In contrast, although financial development exhibits a positive direct association with tourism growth, its interaction with NBDCPZs is statistically insignificant (coef. = −.2472, t = −.8771), suggesting that financial depth does not significantly amplify policy impacts. Overall, the results indicate that human capital plays a more effective moderating role than financial development in shaping the tourism effects of NBDCPZs.
Results of Moderation Effect Tests.
Note. Values in parentheses are z-statistics.
, **, and * indicates significance at the 1%, 5%, and 10% levels, respectively.
Further Robustness Checks and Accounting for the COVID-19 Shock
Although the baseline specification already controls for common time shocks through year fixed effects, the COVID-19 pandemic represents an unprecedented external disruption that may have altered both tourism demand and the effectiveness of digital governance policies. To ensure that the estimated impact of the National Big Data Comprehensive Pilot Zone (NBDCPZ) policy is not confounded by pandemic-related shocks, we conduct an additional robustness check by explicitly modeling the differential policy response during the COVID-19 period.
Specifically, we introduce an interaction term between the NBDCPZ policy indicator and a COVID-19 dummy variable that equals one for the years 2020 to 2021 and zero otherwise. The empirical specification is given by:
where
It is important to note that the COVID-19 dummy is not included as a standalone regressor, as its average effect across cities is fully absorbed by year fixed effects. Instead, the interaction term captures whether the tourism impact of the NBDCPZ policy differs between pilot and non-pilot cities during the pandemic period. This specification follows the logic of a difference-in-differences interaction embedded within a fixed-effects framework and allows us to assess the resilience or vulnerability of data-driven tourism governance under extreme external shocks. This approach does not alter the baseline identification strategy but serves as a robustness check to verify that the estimated policy effect is not driven by pandemic-induced structural breaks.
The results are reported in Table 14. After accounting for the COVID-19 period, the coefficient of the National Big Data Comprehensive Pilot Zone (NBDCPZ) policy remains positive and statistically significant (coef. = .0023, p < .05), indicating that the baseline tourism-promoting effect persists in pilot cities. The interaction term between NBDCPZ and the COVID-19 dummy is significantly negative, suggesting that the policy effect was temporarily attenuated during the pandemic years. This finding reflects the strong external shock imposed by COVID-19 on tourism activities rather than a reversal of the underlying policy effect. Overall, the results confirm that the main conclusions are robust to the inclusion of pandemic-related disturbances. Accordingly, the baseline analysis focuses on the pre-pandemic period to ensure internal validity, while the COVID interaction specification is implemented solely as an auxiliary robustness check rather than an alternative baseline.
Results of COVID-19 Robustness Test.
Note. Values in parentheses are z-statistics.
, **, and * indicates significance at the 1%, 5%, and 10% levels, respectively.
Alternative Interpretations of Scale-Related Effects
The baseline heterogeneity analysis reveals that the tourism-promoting effect of the NBDCPZ policy is insignificant or even attenuated in super-large cities. While this pattern is consistent with the notion of diminishing marginal returns and coordination costs in highly developed metropolitan areas, it may also reflect a nonlinear relationship between city scale and the effectiveness of digital policy interventions. To explore this possibility and to rule out alternative interpretations based on discrete city-size classifications, we conduct an additional robustness check by introducing a nonlinear scale term into the baseline regression. Specifically, we augment the main specification by including both the logarithm of city population and its squared term, allowing the marginal effect of city size to vary along the scale distribution:
In this framework, a statistically significant negative coefficient on the squared term (
The estimation results indicate that neither the linear population term nor its squared term is statistically significant, although their coefficients display opposite signs. Specifically, the coefficient on the logarithm of population is negative (β = −2.059, p = .216), while the squared term is positive (β = 2.074, p = .214), suggesting no empirically identifiable nonlinear scale relationship within the sample period. Although the direction of the estimates is consistent with a potential concave association, the lack of statistical significance implies that no clear population threshold or turning point can be established. Importantly, core economic fundamentals such as GDP remain statistically significant (β = .048, p < .01), indicating that the insignificance of the scale terms is unlikely to stem from model misspecification. Taken together, these findings suggest that the attenuated policy effects observed in super-large cities are better interpreted as reflecting structural maturity and absorptive capacity constraints rather than a mechanically nonlinear population effect. Therefore, the absence of statistical significance indicates that scale-related attenuation is more plausibly driven by institutional saturation rather than a mechanical population-size nonlinearity.
Conclusions and Implications
Conclusions
Based on panel data from 284 prefecture-level cities in China during 2009 to 2019, this study employs multi-period difference-in-differences (DID), propensity score matching (PSM), and spatial Durbin models (SDM) to systematically examine the effects and mechanisms of the National Big Data Comprehensive Pilot Zone (NBDCPZ) policy on tourism economic development. The empirical results demonstrate that the NBDCPZ policy exerts a significant and robust positive effect on urban tourism development, which remains stable after a series of endogeneity corrections, spatial dependence tests, and robustness checks. Additional robustness analyses accounting for the COVID-19 shock further confirm that the baseline policy effect is not driven by pandemic-induced disruptions.
First, the dynamic analysis reveals that the tourism effects of NBDCPZs exhibit a lagged and cumulative pattern over time. The policy impact becomes statistically significant after implementation and continues to strengthen in subsequent periods, indicating that the effectiveness of big data–driven governance in tourism relies on gradual processes of institutional adaptation, data integration, and industrial adjustment rather than instantaneous responses. This finding highlights the importance of sustained policy support in unlocking the long-term dividends of digital transformation. The spatial econometric results suggest that spatial dependence in tourism development mainly operates through market competition rather than policy diffusion. While the significant spatial autoregressive term indicates that tourism outcomes in one city are affected by neighboring cities, the insignificant spatial lag of the NBDCPZ variable implies that policy designation itself does not generate cross-city spillovers. Combined with the negative spatial effect of neighboring GDP, these findings indicate that adjacent cities primarily interact through competitive substitution in tourism markets, rather than through cooperative or emulative policy channels.
Second, the policy effect displays pronounced multidimensional heterogeneity. From a regional perspective, the impact is significantly positive in Northern and Southern China, while remaining statistically insignificant in Central and Northwestern regions. These differences reflect disparities in economic foundations, digital infrastructure conditions, and the capacity to absorb data-driven policy instruments. Regarding city size, the results show that large cities benefit significantly from the NBDCPZ policy, whereas the estimated effects in mega and extra-large cities are insignificant. This pattern reflects the fact that tourism development and data infrastructure in super-large cities are already highly mature and stable, such that the designation of an additional national data pilot zone does not constitute a binding or transformative institutional shock. In this context, the policy effect is diluted not because digitalization becomes counterproductive, but because baseline capacities and tourism demand are already saturated. In this sense, the attenuated policy effect in super-large cities reflects saturation and institutional redundancy rather than a reversal of the digitalization–tourism relationship. Furthermore, distributional evidence from quantile regressions indicates that the policy effect is most pronounced at medium-to-high levels of tourism development, while remaining weak or insignificant at the lower and upper extremes of the distribution, underscoring the nonlinear nature of digital policy impacts.
Third, mechanism analysis identifies industrial structure optimization as the core transmission channel through which the NBDCPZ policy promotes tourism economic development. The policy significantly increases the ratio of tertiary to secondary industry output, and mediation tests confirm that industrial upgrading plays a partial but substantive role in transmitting policy effects. It should be emphasized that the identified mediating role of industrial structure optimization reflects a macro-level structural transmission channel, rather than a direct observation of firm- or consumer-level behavioral mechanisms. Micro-level data would be required to further unpack the behavioral foundations of this structural effect. By contrast, the mediating role of information infrastructure (measured by postal and telecommunication business volume) is insignificant, suggesting that basic connectivity expansion cannot capture the deeper data-driven innovation (Van Nuenen & Scarles, 2021), precision marketing (Li et al., 2022), and tourist analytics emphasized by the policy (Zhang & Cheng, 2024), reflecting a stage-of-development distinction. Basic connectivity indicators capture infrastructure conditions that have largely converged across cities during the sample period, whereas the NBDCPZ policy operates on higher-order data governance and analytics capabilities that remain unevenly distributed and are directly relevant for decision-making optimization rather than mere information transmission within the tourism sector.
Finally, moderation analysis reveals heterogeneous conditioning effects across enabling factors. Human capital endowment significantly amplifies the tourism-promoting effect of the NBDCPZ policy, indicating that regions with stronger higher education capacity are better positioned to translate data resources into tourism productivity gains. In contrast, although financial development exhibits a positive direct association with tourism growth, its interaction with the NBDCPZ policy is statistically insignificant, suggesting that financial depth does not systematically strengthen policy effectiveness once data-driven mechanisms are in place. This pattern implies that, at the city level, financial resources may primarily support general tourism expansion rather than the specific data integration, analytics, and governance capacities targeted by the NBDCPZ policy. When core data infrastructures and institutional arrangements are already determined by policy designation, additional financial deepening alone is insufficient to further amplify marginal policy effects.
Contributions and Policy Implications
The marginal contributions of this study are twofold. First, through systematic empirical analysis, it establishes a direct nexus between the National Big Data Comprehensive Pilot Zone (NBDCPZ) policy and tourism economic development, identifying industrial structure optimization and human capital conditions as the core mediating mechanism. This finding advances the theoretical understanding of how data elements reshape the tourism value chain through sectoral upgrading and capability enhancement mechanisms. Second, by integrating multi-period DID, spatial Durbin, and quantile regression models within a multidimensional analytical framework, the study captures the dynamic, spatial, and distributional characteristics of policy effects, offering a methodological paradigm for evaluating digital economy policies. Nevertheless, as an initial exploration of the NBDCPZ’s impact on the tourism economy, this research has limitations that merit further refinement. The depth of micro-level mechanism analysis remains limited: although the mediating effect of industrial upgrading is verified, the micro-dynamics of data-driven value chain restructuring require finer-grained investigation. Future research should enhance analytical precision through micro-level data collection. Moreover, the data dimension is constrained by reliance on traditional statistical indicators, lacking real-time measures of digital consumption behavior. Subsequent studies should incorporate multi-source dynamic data to construct a tourism digital vitality index and capture the time-varying impacts of digital technology on market recovery, thereby providing stronger empirical support for policy innovation.
Based on the above findings, this study proposes a three-dimensional framework for optimizing the effectiveness of the National Big Data Comprehensive Pilot Zone (NBDCPZ) policy in promoting tourism economic development. First, to address the lagged and cumulative nature of policy effects, a multidimensional monitoring system should be established to track core indicators—such as tourism economic elasticity and industrial transformation rates—in real time, thereby identifying inflection points in policy effectiveness. For pilot zones with implementation cycles exceeding 2 years, third-party evaluations should be conducted, and a policy resource elasticity allocation model should be developed to ensure the precise alignment of policy instruments with regional development stages. Second, regions with stronger economic foundations and digital absorption capacity should emphasize technological integration, innovation, and market-oriented data allocation, including the construction of cross-regional cultural and tourism collaboration platforms. Regions with relatively weaker short-term policy responses should prioritize the deployment of digital infrastructure hubs and talent cultivation systems, supported by interregional transfer mechanisms to strengthen foundational capacities. Regarding city size, large cities should foster smart tourism industrial clusters, whereas mega and extra-large cities should mitigate scale-related coordination inefficiencies through computational resource redistribution and spatial reconfiguration to enhance technology diffusion and resource sharing. Finally, to deepen the core transmission pathway of industrial structure optimization, the establishment of market-oriented platforms for tourism data elements should be accelerated, promoting the full-chain digital transformation of cultural and tourism enterprises and expanding the modern service sector. In addition, the complementary role of human capital should be strengthened by improving the alignment between higher education training and the practical needs of data-driven tourism development, thereby enhancing cities’ capacity to translate digital policies into sustained tourism growth.
In summary, this study empirically demonstrates the promotional role of national big data comprehensive pilot zones in tourism economies and clarifies the heterogeneous conditions under which such effects are most likely to materialize. It provides theoretical underpinnings and policy targets for empowering high-quality tourism development through national digital policies. Future research should deepen micro-level mechanism analysis and multi-source data integration to foster the vigorous development of new tourism productive forces within China’s modernization process.
Footnotes
Author Contributions
Conceptualization: XL, YQ, and ZH-X; Methodology: XL and YQ; Software: XA; Formal analysis: XL, YQ, and ZH-X; Resources: YQ and ZH-X; Writing – original draft preparation: XL and YQ; Writing – review and editing: YQ and ZH-X; Visualization: XA. All authors have read and agreed to the published version of the manuscript.
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research was funded by the Fujian Provincial Social Science Fund, grant number FJ2024C088.
Declaration of Conflicting Interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Data Availability Statement
The datasets used and/or analyzed in the current study are available from the corresponding author upon reasonable request.*
