Domestic tourism demand in the North and the South of Europe in the Covid-19 summer of 2020

This study investigates empirically changes in domestic summer tourism demand following the Covid-19 pandemic in 305 regions across six European countries (Denmark, Finland, France, Italy, Spain and Sweden) based on official data. Five different groups of NUTS 3 regions are identified in accordance with a typology suggested by the OECD where density and connectivity are aspects of importance. Dynamic panel data estimations show that large metropolitan regions experience strong decreases in demand (approximately 30 per cent) both in July and August 2020. There are, however, clear differences between the Northern and Southern European countries. In the North, the remote regions encounter an increased demand that is partially offsetting losses in the large metropolitan regions. This pattern cannot be found in the South. The decline in domestic tourism flows to the major metropolitan areas is also more pronounced in the South of Europe, approximately 50 per cent per summer month compared with 20 per cent (July 2020) and stagnation (August 2020) in the North regions.


Introduction
During the first pandemic summer of 2020, domestic tourism flows develop unevenly across regions. Rural and mountain regions experience an increase in domestic overnight stays and arrivals, while cities attract fewer visitors (OECD 2020). Several studies examine the progress of domestic tourism flows to different regions during the pandemic (Altuntas and Gok;Arbulú, et al. 2021;Falk et al. 2022), 1 3 although none of them use a typology that takes into account both connectivity to large metropolitan or urban regions and population density. Since the Covid-19 pandemic affects the perception of crowds (Zenker and Kock 2020) and assemblies are linked to infection rates (Rader et al. 2020), connectivity to and density of regions are aspects of particular importance in analyses of visitor flows.
This empirical study investigates changes in domestic summer tourism demand following the Covid-19 pandemic in 305 regions across six European countries (Denmark, Finland, France, Italy, Spain and Sweden). These countries in the North and in the South of Europe are chosen because of their large domestic markets for summer tourism. They also have long expanses of coastal lines as well as deserted mountainous inland areas (except in Denmark). Five different groups of NUTS 3 regions are identified in accordance with a typology suggested by the OECD where density and connectivity are aspects of importance. This division leads to two metropolitan and three non-metropolitan regions (Fadic et al. 2019). The large metropolitan region includes a city with more than 1.5 million inhabitants (e.g. Barcelona, Bouches-du-Rhône, all departments of Paris, Madrid, Naples, København and surroundings, Rhône, Seville, Stockholm, Turin and Valencia). Domestic summer tourism is approximated by the number of overnight stays based on data from official sources for the months of July and August during the years 2014 to 2020. Dynamic panel data methods are applied to the tourism demand model in which the time effects are allowed to vary across regions.
Many studies explore the determinants of domestic tourism demand based on regional data (Massidda and Etzo 2012;de la Mata and Llano-Verduras 2012;Patuelli et al. 2013;Massidda and Piras 2015;Campaniello and Richiardi 2018), although most of them use relatively high aggregation levels of data (NUTS 2) and mainly refer to the pre-pandemic period of time. Based on city-level data, Plaza et al. (2015) show that the specific media attention in connection with the opening of the Guggenheim Museum attracts domestic visitors to Bilbao. Population density is a common measure in modelling domestic tourism demand (Massidda and Etzo 2012). This measure characterises rural and urban regions, but does not allow for spatial relationships between regions, for which specific econometric models are required (Romão and Saito 2017).
This study provides a quantification of domestic summer tourism demand across large parts of Europe in times of a worldwide crisis (pandemic) by use of a regional typology where both density and connectivity are built-in aspects of importance (Fadic et al. 2019). Theoretically, the study extends the relationships derived within the fields of medicine and psychology, between crowding, connectivity, population density and infection rates to include also regional tourism flows. The theoretical predictions are contrasted with the observed changes in the tourism flows during the pandemic.
Although the countries differ in size, they have many characteristics in common, making it possible to draw more general conclusions from the study. The Covid-19 pandemic stimulates studies on domestic tourism demand using regional data. This study also contributes to the growing literature in the field by using more detailed NUTS 3 regional data than is commonly employed in the literature (Canale et al.
The implication of the approach used in this study is that not only the opposite ends (urban versus rural) but also the range of regions in between are considered. Another contribution to the literature is the exploration of the possible benefits of using dynamic panel data models explicitly developed for a short time period and a large number of groups as is the case in this study. Dynamic panel data models also have the advantage of controlling for persistence of domestic tourism flows.
After this introductory section, the conceptual background is presented, followed by sections on the empirical approach, data and descriptive statistics, results and concluding remarks.

Conceptual background
This study focuses on the development of domestic summer tourism demand in the first year of the Covid-19 pandemic (2020). During the summer period, there are fewer pandemic related restrictions on mobility, in particular domestically (Hale 2021). 1 Literature provides several theories on possible changes in domestic travel flows during the Covid-19 pandemic. Zenker and Kock (2020), for instance, argue that the pandemic influences tourism flows in various ways, such as affecting the image of cities as tourist destinations negatively because they were the first to report high infection rates. Another aspect is that the threat of pathogens makes people more alert and adverse to crowds (Wang and Ackerman 2019;Zenker and Kock 2020). This could cause a change in travel behaviour towards more remote, less populated destinations (Zenker and Kock 2020).
According to Florida et al. (2021), large cities with international airports exhibit much higher infection rates than rural areas during the first phase of the pandemic (see also Nathan and Overman 2020). The higher overall incidence rate and more rapid speed of the virus in densely populated areas are confirmed by the medical literature (Rader et al. 2020). Sy, White and Brooke (2021) show that population density in the USA is significantly positively related to disease transmission, independent of transport connectivity and average income. This suggests that the possibility of effective spread is primarily driven by close contacts in densely populated areas, rather than by connectivity and mobility between different geographical zones. Other studies also demonstrate that urbanisation is the factor most strongly correlated with the number of Covid-19 cases and deaths (Viezzer and Biondi 2021).
In the consecutive phases of the Covid-19 pandemic, the urban-rural divide for infection rates is less pronounced, indicating that the type of density rather than the scale is of importance, compactness of jobs versus residential density, for instance (Florida et al. 2021;Florida and Mellander 2022). Proximity to or connectivity across regions is another possible dimension for the spread of the virus (Florida and Mellander 2022). These aspects are crucial for tourism flows, since it is easier to avoid crowding in remote or disconnected than in urban areas.
Several restrictions on large meetings and events, cancellation of conferences and trade fairs as well as increased use of digital meetings are imposed to curb the spread of the Covid-19 pandemic (Table 1; Hale et al. 2020;Sharifi and Khavarian-Garmsir 2020).
The regulatory measures lead to a decline in business customers, which are particularly common in larger cities and metropolitan regions. Besides the fear of the virus and the need to keep a distance, the partial closure of cultural attractions also makes visiting museums, galleries, restaurants and bars indoors less attractive, even if they are open, something that affects highly populated areas the most.
Early classifications of regions are either based on tourism intensity, continuous measures of population density or on a dichotomy of urban versus rural (Capone and Boix 2008;Santana-Jiménez and Hernández 2011;Marrocu and Paci 2013;Dijkstra and Poelman 2014;Laurin et al. 2020). Veneri and Ruiz (2016) distinguish between predominantly urban, intermediate and primarily rural regions on the basis of population density per local unit. The OECD is involved in several attempts to classify regions that go beyond the simple dichotomy of urban and rural ranging from semi-urbanised to remote peripheral areas (Dijkstra, Poelman and Veneri 2019). Remote peripheral regions, for example, are generally rural in character, but have two distinct characteristics-rurality and geographical remoteness (Brezzi, Dijkstra and Ruiz, 2011). This typology is based on the proximity to functional urban areas and distinguishes between: (i) predominantly urban regions (within a functional urban area), (ii) intermediate regions (close to a functional urban area) or (iii) predominantly rural regions (far from a functional urban area).
In the specific case of analysing domestic demand for summer tourism during a pandemic, a classification is required that takes into account not only urban versus rural, but something more detailed including the degrees of density and possibly also connectivity. Fadic et al. (2019) take the earlier typologies further by Table 1 Cancellation of public events and restrictions on gatherings. Source: Hale et al. 2020; Definition of gatherings: 0-no restrictions, 1-restrictions on very large gatherings (> 1000 persons), 2-restrictions on gatherings between 101 to 1000 persons, 3-restrictions on gatherings between 11-100 persons and 4-restrictions on gatherings of 10 persons or less. Definition of public events: 0-no measures, 1-recommend cancelling and, 2-require cancelling

Cancellation of public events
Restrictions on gatherings July 2020 August 2020 July 2020 August 2020 Domestic tourism demand in the North and the South of Europe… identifying five categories of regions including population density, urban accessibility as well as connectivity. Access is defined in terms of the time it takes to reach the nearest urban area, based on both geographical features and the state of the physical road infrastructure. Thus, the five-region typology allows a detailed examination of domestic tourism flows during different waves of the pandemic (Table 2). This classification is recently used to characterise both the degree of rurality and connectivity and estimations show that only 8 per cent of the OECD population live in remote regions where there is little interaction with nearby cities (Garcilazo, Moreno-Monroy and Martins 2021).
The course of events following the Covid-19 pandemic, with mobility and gathering restrictions as well as an initial strong spread of the virus in densely populated areas, as discussed in the literature, is expected to change the pattern of demand for domestic summer tourism in 2020. This change in demand may be apparent, but its extent and variations across regions can only be identified empirically when regional characteristics and other factors of importance are taken into account. Thus, the main hypotheses of the study can be formulated: H1: Demand for domestic summer tourism exhibits a new pattern in 2020 compared with earlier years.
H2: Demand for domestic summer tourism evolves differently across regions in 2020.
H3: Demand for domestic summer tourism in the North of Europe is affected dissimilarly by the Covid-19 pandemic than in the South.
The validity of the hypotheses is tested on regional data for the three Nordic countries (Denmark, Finland and Sweden) and the three South European countries (France, Italy and Spain) by use of a dynamic panel data model.

Empirical model
There are several methodological approaches for analysing changes in regional tourism flows. One popular approach is the shift-share analysis that is commonly used in the field of regional economics (Firgo and Fritz, 2017;Costantino et al. 2021). The advantage of this approach is that it does not require the formulation of an econometric model with its underlying assumptions. A clear disadvantage of shift-share analysis is that different control factors and stochastic factors cannot be taken into account. An alternative approach presented in the regional economics literature is the tourism gravity model, in which bilateral tourism flows depend on the income of the origin and destination regions, as well as on the actual distance between these two geographical points (De la Mata and Llano-Verduras 2012; Cafiso et al. 2018). Although many merits, this estimation method is seldom applicable for domestic travel flows due to lack of data. Yet another approach derives a regional tourism demand equation depending on prices, real income (GDP) and control factors such as population density (Massidda and Piras 2015). Given the research question and data availability, this latter approach is selected for the present study. Domestic tourism demand depends on economic (real income and prices) as well as non-economic factors.
The literature exhibits several examples where dynamic panel data models are used to estimate the determinants of domestic tourism flows at the regional level (Massidda and Etzo 2012;Canale et al. 2019). These kinds of models have the advantage of controlling for time-invariant regional factors such as natural amenities, which are important determinants of destination attractiveness (Patuelli et al. 2013;Naranpanawa et al. 2019). In dynamic panel data models, it is also possible to control for persistence in tourism demand . By taking the typology of regions into account (Fadic et al. 2019), demand for domestic summer tourism across the six countries is modelled as follows: where ln DNS it is the natural logarithm of the number of domestic overnight stays in the region i (i = 1, 2, …, 305) in the months of July or August of year t. Y ct denotes national GDP for the first two quarters in constant prices and CPI ct is the national consumer price index for the first two quarters in country c. Parameters 1 and 2 are the short-term income and price elasticities of tourism demand. Regional fixed effects are captured by i , which is part of the residual it ≡ i + u it . A set of time dummy variables d t are introduced to capture presumptive effects of the pandemic in 2020. The years 2018 and 2019 are treated as reference years.
The development in 2020 is modelled as the products of the 5 regional and the 2020 dummies, d MRL,20 ,d MRM,20 ,d NMRM,20 ,d NMRS,20 ,and d NMRR,20 reflecting the variations in kinds of regions in terms of their density, remoteness or connectivity. Parameters of these dummy variables inform if there is a distinct change in demand for domestic summer tourism in 2020, as compared with the average development in the base years 2018-2019.
Due to the correlation between the lagged lnDNS and the residual it , the ML estimator of the Dynamic Fixed Effects panel data model developed by Hsiao et al. (2002) is employed where all continuous variables are used in their first differences. In comparison with other estimators, such as the GMM, the ML estimator has a better finite sample property (Hayakawa and Pesaran 2015). At the same time, it ln DNS it = + ln DNS it−1 + 1 ln Y ct + 2 ln CPI ct + is designed to fit the datasets with a relatively large number of groups and a short time period, which coincides with the characteristics of the data in this study. This estimator is also robust to the cross-sectional heterogeneity (Hayakawa and Pesaran 2015). A common practice is that the beta coefficients can be interpreted as shortrun (semi-) elasticities. Dividing the beta coefficients by 1 minus the coefficient of the lagged dependent variable ( ∕( 1−γ) ) gives the long run elasticity.
An important aspect in any estimation is the presumptive appearance of endogeneity. The possible endogenous relationship between tourism demand, prices, and income is highlighted in the literature (Assaf et al. 2019;Pesämaa et al. 2021). Endogeneity may arise from reverse causality between tourism flows and economic activity. However, by using regional tourism flows and economic activity at the country level, this problem does not appear as tourism flows in a region are unlikely to affect aggregate GDP growth. What is more, the time dummy variables are exogenous, and the OECD classifications referring to population density, remoteness or connectivity remain unchanged at least in the short and medium term.
The specification does not take spatial dependence into account in the estimation method (Patuelli et al. 2013;Romão and Saito 2017;Pompili et al. 2019). This is valid in the period of Covid-19, given the restrictions or advice against travel even to major parts of the neighbouring countries. Another explanation behind this choice is the definition of regions used, that to some extent already includes spatial elements. Results will be presented for two groups of countries, Northern and Southern Europe, as beyond the features in common they also have dissimilarities relating to for instance size, attractiveness, number of regions and climate zone.

Data and descriptive statistics
This study employs regional panel data on the number of domestic overnight stays in Denmark, Finland, France, Italy, Spain and Sweden for a series of summer months (July and August 2013-2020). Data originate from the national statistical offices in each country and cover 305 regions and 2135 observations. The dependent variable appears in the logarithm of the number of domestic overnight stays for either July or August at the NUTS 3 level. Domestic overnight stays data for Finland, and Sweden are based on all accommodation establishments (hotels, holiday villages, youth hostels, camping sites and commercially arranged private cottages and apartments) while data for Denmark are more limited and only encompass information on hotels with minimum 40 beds. 2 French data on domestic overnight stays, excluding camping sites, originate from the Statistical Office of France (INSEE). 3 Spain, as well as Italy, hold registers on hotels and similar accommodation, holiday and other 1 3 short-stay accommodation, camping grounds, recreational vehicle parks and trailer parks by community and province. 4 Information on gross domestic product in constant prices refers to the first two quarters of each year, originating from OECDSTAT (Dataset: Quarterly National Accounts: Volume and price indices-GDP expenditure approach, series VOBARSA). 5 The consumer price information relates to the first six months of each year and is also brought from OECDSTAT. 6 Descriptive statistics show that the evolution of domestic overnight stays is uneven across the five regions in the first summer of the Covid-19 pandemic (Table 3A, B, C and Fig. 1). As opposed to the trend of earlier years, large metropolitan regions experience a strong decline in the number of domestic overnight stays in 2020, while remote regions suddenly reverse their development to an increase and thus are the only ones that are growing in July 2020. In August, the downturn is slightly less pronounced, but follows the same pattern as in July, with the exception of the NMRM region in the middle of the scale, which also experiences a slight upswing in demand. There are also apparent variations between the countries in the north and those in the south. Despite the downturn in domestic demand in metropolitan areas, the Northern countries exhibit a slight overall increase in domestic overnight stays both in July and in August compared with the year before. The opposite situation is found for the Southern countries, where demand is contracting in 2020, although the magnitude is somewhat smaller in August. Multivariate tests for equality of the mean values of the five groups show that the change in the number of domestic overnight stays in the summer of 2020 compared with the same period in 2019 differs significantly across regions (Fig. 1, p-value < 0.01).

Empirical results
Dynamic panel data estimations for the 305 regions reveal that there is a change in the pattern of domestic summer tourism demand in the pandemic summer of 2020 (Tables 4 and 5). The long-term increase in demand for summer tourism to large metropolitan areas (MRL) is broken both in the South and the North of Europe, while tourism flows to remote regions (NMRS and NMRR) increase. Given the high degree of persistence, the interpretation focuses on the short-term effects. Coefficients for metropolitan areas are approximately -0.34, suggesting that the number of domestic overnight stays is 30 per cent lower in the short term than in the reference period. Estimates of the dummy variables for the remotest regions (NMRR) are 0.20 in July and 0.12 in August. This implies that domestic tourism demand increases in the short-run by 18 and 12 per cent, respectively (calculated according to Halvorsen and Palmquist 1980). 5 https:// stats. oecd. org/ index. aspx? query id= 60702#. 6 https:// stats. oecd. org/ viewh tml. aspx? datas etcode= PRICES_ CPI& lang= en.

3
Domestic tourism demand in the North and the South of Europe…  Overall, the results indicate that not only population density but also the degree of connectivity to urban areas influences domestic tourism demand when a sudden shock like the Covid-19 pandemic occurs. This is a new finding in the literature. Thus, the first two hypotheses cannot be rejected.
Hidden behind the overall figures, there are more pronounced variations between regions in the South and the North of Europe. This means that the third hypothesis cannot be rejected either. In July, both metropolitan regions (MRL and MRM) in the South experience a decrease in demand for tourism, although the impact is strongest on the larger regions. As opposed to the regions in the North, there is also no direct significantly positive effect on the overnight stays in the three remote regions in the South. In the North, these three remote regions are all gaining in overnight stays, while the next largest metropolitan area is not directly affected and only the largest regions experience a decline.   The decrease in demand for the largest metropolitan regions in the South is equally strong in August, while there is no longer a significant direct effect in the North, as compared with the reference period. Instead, the demand for overnight stays surges in all regions in the group of Northern countries, but the effect is no Table 4 Demand for domestic overnight stays in July, Dynamic panel data estimations. Source: National statistical offices, OECD Stats and own calculations Asterisks ***, ** and * denote significance at the 1, 5 and 10 per cent levels. The standard errors are reported within the parentheses. The estimations are based on robust standard errors and conducted using the Stata command xtdpdqml (Kripfganz, 2016) that implements the QML method of Hsiao et al. (2002). The regions are defined as the following: with a very large city > 1.5 M inhabitants (MRL), with a large city > 250 K inhabitants (MRM), near a city > 250 K inhabitants (NMRM), near a city between 50 and 250 K inhabitants (NMRS) and remote region (NMRR 1 3 Domestic tourism demand in the North and the South of Europe… longer strongest in the utmost remotest areas, but in the third remotest regions (NMRM). Coefficients for real GDP and CPI show the expected signs, but the significances are weaker for the Southern countries. However, the income and price elasticities should be interpreted with caution as they only vary across time and countries. The coefficients for lagged domestic overnight stays range between 0.50 and 0.91 indicating a relatively high degree of persistence, which is also consistent with previous Table 5 Demand for domestic overnight stays in August, Dynamic panel data estimations. Source: National statistical offices, OECD Stats and own calculations Asterisks ***, **, and * indicate the significance at 1%, 5%, and 10%, respectively. Points estimates are reported for coefficients. The standard errors are reported in the parentheses. The significance is based on the t-test. The regions are defined as the following: with a very large city > 1.5 M inhabitants (MRL), with a large city > 250 K inhabitants (MRM), near a city > 250 K inhabitants (NMRM), near a city between 50 and 250 K inhabitants (NMRS) and remote region (NMRR). See Table 4 All 1 3 studies (Gil-Alana and Huijbens 2018; Falk et al. 2022). Besides this, there are no known similar studies to compare present results with. As a robustness check, the same specification is estimated with four classifications (three non-metropolitan and one metropolitan region). Unreported results show that these estimations are largely consistent with those presented. In addition, other dynamic panel data methods are used such as the system GMM estimator (Blundell and Bond 1998;Roodman 2009). These estimations give similar results and are available upon request. However, results for the Nordic sample should be interpreted with caution as the system GMM estimator requires a relatively large number of groups which is not fully satisfied for this area.

Conclusion and discussion
This study investigates changes in domestic summer tourism demand following the Covid-19 pandemic in 305 regions across six European countries (Denmark, Finland, France, Italy, Spain and Sweden) based on official data. These countries are chosen because of their large domestic markets for summer tourism. They also have long expanses of coastal lines as well as deserted mountainous inland areas (except Denmark). Literature in medicine and psychology on the relationship between crowding, connectivity, population density and infection rates is used to formulate the hypotheses on the shift in regional tourism flows. Five different groups of NUTS 3 regions are identified in accordance with a typology suggested by the OECD where density and connectivity are aspects of importance. Dynamic panel data estimations reveal that the pattern of domestic tourism demand during the summer of 2020 is uneven across regions and differs from the trend of earlier years. Large metropolitan regions face a marked decline of more than 30 per cent in both July and August, when real GDP (capturing the impact of the Covid-19 recession), prices and the level of overnight stays in the same summer months earlier years are controlled for. In contrast, remote areas experience a short-term increase of 18 per cent in July and 12 per cent in August.
The downturn in domestic tourism demand for metropolitan areas can partly be explained by the absence of events and attractions due to restrictions with the purpose to calm down the pandemic. Densely populated areas could also become frightening for visitors per se. In such situations, rural and remote regions offer more space and opportunities for individual outdoor activities, which might be far more appealing. However, only in the North, there are indications of a substitution in demand from metropolitan to remote areas. While all countries experience a pronounced downturn in demand for tourism to large metropolitan areas, only the countries in the North meet increased demand for remote areas. The magnitude of the decline is also larger for the regions in the South of Europe.
Several conclusions can be drawn from the study. First, regional typology indicators are stronger determinants of domestic tourism and travel demand during the Covid-19 pandemic than classical economic factors such as income and prices. Second, an in-depth distinction of different types of regions that integrate more aspects than population density is important. Third, the results suggest that more aggregated regional data such as the NUTS 2 level are not sufficient to explain variations in domestic travel flows.
Some limitations need to be mentioned. The study covers a short time period, only the first summer of the Covid-19 pandemic. Domestic tourism demand may vary over different stages of the pandemic. Future work should apply the model to other regions of the world. In addition, the remote region can be further refined by distinguishing between coastal and mountain regions.
Funding Open access funding provided by Södertörn University.

Declarations
Conflict of interest There are no financial interests or benefits that have resulted from the direct application of our research.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.