1 Introduction

The pandemic that started in 2020 has led to the first world-wide economic downturn in recent times triggered by a deadly virus. As the outbreak of COVID-19 commenced, governments were confronted with the dilemma of how to balance the economic and health costs of a surge in COVID-19 cases versus the costs of preventive health measures to stop the spread of the virus (Layard et al., 2020). Two major complementary mechanisms for how COVID-19 may have affected short-run labour market outcomes are investigated in this paper.

First, there is a direct economic effect arising from the population’s health concerns due to COVID-19, which may differ by location, as people living in so-called COVID-19 hotspot areas may be more aware of the presence and detrimental consequences of the virus. As a result, they would take voluntary preventive measures of social distancing more seriously, which is likely to have a negative impact on social activity and the labour market. Second, the indirect economic effect through (regional) enforced lockdown and social distancing regulations by the government in response to the virus would lead to an immediate loss of the economic activities that are no longer allowed, and a negative impact on labour market outcomes of workers who can be dismissed easily.

This paper assesses the relevance of the above two mechanisms for how the COVID-19 outbreak has affected the labour market in the Netherlands. We use unique Dutch administrative weekly panel microdata covering the period until the end of March 2020, drawing a random sample of 50 per cent of all Dutch employees (about 4.2 million individuals). Taking a regional perspective, we examine whether the economic slowdown as measured by individual labour market outcomes has been stronger in COVID-19 hotspot areas where a substantial proportion of the population was affected by the virus. Such a geographical examination is important as different labour market changes across regions suggest that it may be worthwhile to consider imposing preventative measures, such as a lockdown, at the regional (rather than national) level, isolating hotspots from low-risk areas.Footnote 1 This is illustrated by the Australian case, where the economic recovery after the first lockdown in March–May 2020 was considerable in all states except for Victoria which experienced a second wave of COVID-19 and local ‘hard’ lockdown from July until October 2020 (Kalb, Guillou, and Meekes, 2020).

The labour-market effects during the COVID-19 outbreak have been documented for many countries by a rapidly expanding literature. The general picture that emerges from these studies is that there was a steep decline of the labour market outcomes in the first months of the outbreak. Furthermore, there was an unequal impact across workers– employees with relatively weak employment protection experienced worse labour market outcomes. We report the outcomes for employment.Footnote 2 For the US, the employment-population ratio fell by about 8 percentage points from February to April, where low-wage workers and minority workers were hit hardest (Coibion et al., 2020; Hershbein & Holzer, 2021). In Japan, in May 2020, regular employment declined by 1 per cent; for workers in non-standard jobs the decline was 4–5 per cent (Kikuchi et al., 2020). In Greece, employment had declined by about 12 per cent by June (Betcherman et al., 2020). For the Netherlands, Von Gaudecker, Holler, Janys, Siflinger, and Zimpelmann (2020), who used monthly survey data on 2918 salaried workers until late March 2020 from the Longitudinal Internet studies for the Social Sciences (LISS) panel, found a reduction in total working hours of 11 per cent or 3 h. In Australia, employment declined by about 6.7 per cent from March to May (Borland & Charlton, 2020). In Mexico, the formal job market declined by 5.4 per cent from March to November (Hoehn-Velasco, Silverio-Murillo, and Balmori de la Miyar, 2021).

Evidence on within-country regional differences in COVID-19 impacts on labour market outcomes is more limited.Footnote 3 For South Korea, Aum, Lee, and Shin (2020) investigate regional differences in the COVID-19 impact on the labour market. They use a differences-in-differences estimator to compare the COVID-19 impact for the local area Shincheonji –which experienced a significant COVID-19 outbreak from February 18 onwards without imposing a lockdown– with other areas that did not experience a significant number of COVID-19 cases. Aum et al. (2020) find that a one per thousand increase in infections caused a 2 to 3 per cent drop in local employment in the initial months after the outbreak (up to May 2020). As South Korea has not imposed any lockdowns, they argue that the estimates imply that at most half of the 5 to 6 per cent decrease in employment in the US and UK can be attributed to lockdowns. For the US, Cho, Lee, and Winters (2020) show that the negative impact of the COVID-19 crisis on employment is larger for metropolitan areas, arguing that the local COVID-19 infection rate explains half of the heterogeneity in employment changes across US metropolitan areas.

2 Regional Differences in COVID-19 Cases and Preventive Measures

The outbreak of the COVID-19 pandemic has had a common pattern across countries. First, there was a phase of denial by the authorities and the public, downplaying the severity of the outbreak. Second, a substantial outburst of cases occurred in a local region. Governments responded by introducing preventive measures for this local area only. Third, COVID-19 cases spread to other parts of the country. During the outbreak of the virus, most governments started by imposing a regional lockdown, before broadening it to the entire country. Of all policy measures, the compulsory societal lockdown was the most disruptive to the economy, enforcing social distancing, staying at home and working from home rules.

The virus outbreak in the Netherlands followed this pattern. On February 27 2020, the first person tested positive. In the first weeks of March, the southern province Noord-Brabant had about half of all detected infections in the Netherlands despite this province only accounting for 15 per cent of the Dutch population. At the same time, the northern provinces were almost free of infections. With regard to reported confirmed COVID-19 cases, hospitalizations and deaths, Noord-Brabant was leading in absolute terms per 100,000 residents (Fig. 1) as well as in relative terms as a proportion of total Dutch confirmed COVID-19 cases (Fig. 2). Consequently, the government’s preventive measures were at first directed at Noord-Brabant only. On March 6, people living in this province were advised to stay home, particularly if they had colds, coughs or a fever. On March 9, the Dutch Prime Minister suggested the population of Noord-Brabant should work from home. On March 10, large gatherings were banned in Noord-Brabant. On March 12, restrictions were imposed on the entire country, including social distancing, banning of gatherings over 100 persons, and a work-from-home directive. From March 15 onwards, all restaurants, schools, childcare and sport facilities were closed. On March 23, physical distancing requirements were tightened, imposing the 1.5-m distance measure and cancelling all gatherings including those with fewer than 100 people.

Fig. 1
figure 1

Cumulative number of cases per 100,000 residents by province, March 2020 See Rijksinstituut voor Volksgezondheid en Milieu (RIVM) (2020) for the COVID-19 cases data. See CBS (2020) for the population data

Fig. 2
figure 2

Proportion of cases in Noord-Brabant (NB), relative to the total number of cases in the Netherlands, March 2020 Notes: See RIVM (2020) for the COVID-19 cases data. See CBS (2020) for the population data

3 Statistical Identification

We estimate the causal impact of the COVID-19 outbreak on three individual labour market outcomes: employment, measured by a 0–1 indicator which equals one if the person was employed–at least partly–in a given week; (logarithm of) the number of hours worked; and (logarithm of) the hourly wage.Footnote 4 At the individual level, the information on working hours and gross wages is derived from monthly income statement data. If a calendar week sits across two calendar months, data from the first calendar month is used. Weekly variation in hours worked and hourly wage is driven by changes in employment only (from job to job or to unemployment).

Identification of both mechanisms is based on a specific regional pattern during the virus outbreak in the Netherlands. In the first weeks of March, the southern province Noord-Brabant had about half of all detected infections in the Netherlands despite this province only accounting for 15 per cent of the Dutch population. At the same time, the northern provinces were almost free of infections. With regard to reported confirmed COVID-19 cases, hospitalizations and deaths, Noord-Brabant was leading in numbers per 100,000 residents as well as in terms of the proportion of total Dutch confirmed COVID-19 cases.

For a panel of individual employees, we apply a difference-in-difference specification at the weekly level by interacting the 0–1 indicator for 2020 (which is set to zero for 2019 observations) with 0–1 indicators for each of the first thirteen calendar weeks of the year. The interaction terms are used for a comparison of the outcome variable by calendar week relative to week 9–the week of the COVID-19 outbreak in the Netherlands. For each outcome variable, the specification is

$$ y_{it} = \alpha_{i} + \mathop \sum \limits_{\begin{subarray}{l} \tau = 1 \\ \tau \ne 9 \end{subarray} }^{13} \beta_{\tau } DW_{\tau } + \mathop \sum \limits_{\begin{subarray}{l} \tau = 1 \\ \tau \ne 9 \end{subarray} }^{13} \gamma_{\tau } DY_{c} \times DW_{\tau } + \delta DY_{c} + \eta \prime X_{ic} + \varepsilon_{it} $$
(1)
$$ \{ i \in 1, \ldots ,N;\,\,t \in 1, \ldots ,13\;\,\,\text {for} \, c = 2019;\;\,\,t \in 14, \ldots ,26\;\,\,\text {for} \, c = 2020\} $$

where y is the outcome variable; the subscripts i, t and c refer to individual, week and year, respectively. \(\alpha\) is an individual fixed effect; \(\tau\) represents the calendar week number; DY and DW are 0–1 indicator variables for year and calendar week. \(\varepsilon\) is an idiosyncratic error term.

The vector X contains 50 variables which are time constant within a year but may vary between the two years, and which is included to reduce the impact of any variation in observables by calendar year. X includes dummy variables for age (6 categories), job characteristics (type of contract (2), type of job (4), full-time/part-time status (2)).Footnote 5 These variables are all measured in calendar week 9, preventing any endogeneity issues resulting from changes in covariates because of COVID-19. Additionally, X contains dummy variables for firm characteristics (size (3), economic sector (20) and a dummy variable for missing firm data although less than 1 per cent of observations fall in this category) and for household characteristics (married (1) and home location (11 provincial regions)), which are all measured on 31 December of the previous year. The results provided in Fig. 3 are robust to excluding \(X\) and are available upon request.

Fig. 3
figure 3

COVID-19 Difference-in-Difference (DD) effects on employment, log hours worked and log hourly wages (Eq. (1)) Notes: Parameter estimates of the double interaction terms between year and calendar week. Each graph represents a single regression for a different outcome variable. Reference year is 2019 and reference calendar week is 9. The 95% confidence intervals are computed based on standard errors clustered by individual. The total number of estimated parameters equals 75

Besides estimating baseline Eq. (1), we estimate a corresponding heterogeneous difference-in-difference equation to investigate the heterogeneity of the COVID-19 impacts. This model complements (1) by also including triple and double interactions between year, calendar week, and all variables in vector XFootnote 6:

$$ \begin{aligned} y_{it} = \alpha_{i} + & \mathop \sum \limits_{\begin{subarray}{l} \tau = 1 \\ \tau \ne 9 \end{subarray} }^{13} \left[ {\beta_{\tau } DW_{\tau } + \gamma_{\tau } DY_{c} \times DW_{\tau } + (\kappa_{\tau } \prime X_{ic} )DW_{\tau } + (\lambda_{\tau } \prime X_{ic} ) \times DY_{c} \times DW_{\tau } } \right] \\ & + (\mu \prime X_{ic} ) \times DY_{c} + \delta DY_{c} + \eta \prime X_{ic} + \varepsilon_{it} \\ \end{aligned} $$
(2)

where \(\kappa_{\tau }\)and \(\lambda_{\tau }\) and \(\mu\) are additional parameters to be estimated, with vector \(\lambda_{\tau }\) including the key parameters of interest.

4 Data

We use administrative data from Statistics Netherlands. For computational reasons, we take a 50 per cent random sample of Dutch employees. We select two cohorts of employees who are followed from January until March. Specifically, we include employees who were employed in calendar week 9 of 2019 and we select employees who were employed in calendar week 9 of 2020. We follow individuals from January 1 until March 31, calendar week 13, of each calendar year. Table 1 reports the individual summary statistics for calendar week 9 of 2019 and 2020, respectively. Although there are some small differences in background characteristics, on the whole, the employee population in week 9 is quite similar in 2019 and 2020.

Table 1 Individual summary statistics by year (proportions unless otherwise noted)

5 Results

We report on three sets of novel results. First, the estimated week effects of Eq. (1) show that employment decreased slightly in weeks 10 to 12 and more substantially, by about 2 percentage points, in week 13 (Fig. 3). The evidence shows a slightly higher employment rate (0.1 to 0.2 percentage point) in weeks 1 to 8 in 2020 (relative to 2019). For the number of paid working hours, a comparable development is observed with a 1.5 per cent decrease in hours in week 13. We observe a small negative effect from COVID-19 on hourly wages of about 0.3 per cent.Footnote 7 The reported COVID-19 effects on working hours and hourly wages are conditional on employment. The parallel trends restriction holds for the models on employment and hourly wages, but not for the model of working hours due to the fact that 2020 is a leap year. Including zeros for the unemployed, which limits the impact of selection into employment, Fig. 4 shows a reduction of 1.75 h in monthly working hours and a reduction of 0.25 euro in hourly wages. Relative to monthly mean working hours of 127 and a mean hourly wage of 19.65 in week 9 of 2020, this represents a decrease of 1.4 per cent and 1.3 per cent, respectively. The very small positive results for weeks 1 to 8 indicate that the labour market in 2020 was very similar to the labour market in 2019 before COVID-19 arrived; if anything, employment had been more stable for those employed in week 9 in 2020 than in 2019.

Fig. 4
figure 4

COVID-19 Difference-in-Difference (DD) effects on hours worked and hourly wages including zeros for the unemployed (Eq. (1)) Notes: Parameter estimates of the double interaction terms between year and calendar week. Each graph represents a single regression for a different outcome variable. Each outcome variable is in levels and zeros are used for unemployed individuals. Reference year is 2019 and reference calendar week is 9. The 95% confidence intervals are computed based on standard errors clustered by individual. The total number of estimated parameters equals 75

The effects of COVID-19 on working hours are consistent with, but somewhat smaller than, those reported by Von Gaudecker et al. (2020) as they find a reduction in total working hours of 11 per cent or 3 h. This difference in results could be explained by changes to respondents’ actual working hours whereas in our study employees’ paid working hours as per income statement could remain the same. Overall, the evidence suggests that employment and working hours are the relevant margins of labour adjustment rather than hourly wages in the first response to the COVID-19 shock.

Second, separate estimation of Eq. (1) for each of the provinces indicate small regional differences in the changes in the outcome variables (Fig. 5). Importantly, these regional differences do not seem to be strongly related to the Dutch COVID-19 hotspot provincial areas of March 2020.

Fig. 5
figure 5

COVID-19 Difference-in-Difference effect stratified by province (Eq. (1)) Notes: Each graph represents a different outcome variable and each line represents a single regression for a different province. Several provinces are left out from Fig. 5 to ensure clear graphs. The figure for all provinces is available from the authors upon request

Third, Table 2 presents the estimated parameters on the triple difference interactions of Eq. (2) revealing which categories of employees had the strongest decline in the three outcome variables in week 13 of 2020. Consistent with the results provided in Fig. 5, the evidence in Table 2 does not suggest a region-specific impact of COVID-19 on the outcome variables when comparing COVID-19 hotspots such as Noord-Brabant and Limburg to other Dutch provinces, including some of the northern provinces such as Groningen and Friesland which had very few cases but experienced a larger negative impact on employment than Noord-Brabant. Other characteristics of employees are shown to be more relevant. Individuals who (i) are aged below 20 years, (ii) have a non-permanent contract; and (iii) are in a flexible or payrolling job, were most negatively affected by the economic effects of the COVID-19 shock. Overall, for a country with a relatively small area size like the Netherlands, the results suggest that the employee’s job characteristics are more important than the regional location of residency for the effects of COVID-19 on individual labour market outcomes.

Table 2 The role of observed individual characteristics in the effects of COVID-19 (Eq. (2))

6 Conclusion

Our study examined the impact of the COVID-19 pandemic in the Netherlands in the initial month after the outbreak, analysing the role of coronavirus hotspot areas in the impact on local labour markets. There are two major outcomes. First, we find that employment decreased by 2 percentage points in March 2020, which is a relatively low reduction compared to studies from other countries. However, the stronger negative effects in other research were estimated using information covering a longer period of time after the COVID-19 outbreak commenced.

Second, the results in this paper indicate limited regional differences in the negative impacts on the Dutch labour market during the outbreak of COVID-19. It appears that higher virus case numbers did not reinforce the decline of the labour market beyond the impacts from the government-enforced lockdown. As a result, the northern Dutch provinces, which experienced a limited number of COVID-19 cases, suffered a similar (or even worse) decline in labour market conditions, as compared with the provinces that were severely affected by the virus.

We examined the net outcome of employment, but did not disentangle the supply and demand side impacts in the labour market. Neither did we exploit any sectoral specialization in economic activities across regions. These factors could potentially contribute to a (partial) explanation of the within-country differences in COVID-19 impacts on regional labour markets, but investigation of these is outside the scope of this paper.

The outcome of an absence of regional differences in the first month of the outbreak suggests that policy makers should be cautious when implementing preventive measures nationwide as the economic costs can be substantial. Thus, where feasible, preventive health measures should be at the regional level, isolating hotspots from low-risk areas. This would allow relatively unaffected parts of the country to continue economic activities as much as possible, ultimately benefitting the nation as a whole.