1 Introduction

The first case of coronavirus disease 2019 (COVID-19) was reported in December 2019, and COVID-19 was declared a pandemic by the World Health Organization on 13 March 2020. Among its many socioeconomic impacts, the pandemic had a negative economic impact at the regional and global levels. Infectious diseases such as COVID-19 differ from natural hazards and disasters as they cause economic impact rather than physical damage (Ceylan et al. 2020). Severe acute respiratory syndrome coronavirus 1 (SARS-COV) dramatically and quickly, albeit briefly, decreased consumer spending (Siu and Wong 2004; Yang and Chen 2009), while the Middle East respiratory syndrome coronavirus (MERS-COV) caused a decline in both tourism spending (Joo et al. 2019) and retail sales (Jung and Sung 2017).

Previous studies on the economic impact of infectious diseases resulting in disasters—including SARS-COV, MERS-COV, and COVID-19—have mainly focused on the declining gross domestic product (GDP) of a country or the declining employment rate and sales in specific industries (Yang and Chen 2009; Ceylan et al. 2020; Shafi et al. 2020; Fairlie and Fossen 2021). They suggested policy implications that can mitigate the negative effect based on the time series analysis of the economic impact of the diseases. However, these studies were unable to identify forces that can absorb and sustain an unexpected shock such as the COVID-19 crisis. Additionally, the existing studies related to the economic impact of COVID-19 do not reflect the differences in industrial structure at the regional level, as they were conducted at the community or village level or focused on economic impact at the national level.

This study focused on the concept of resistance as the framework of the economic resilience. The concept of resilience, recently utilized in the field of disaster studies, provides a theoretical framework for this research question. In particular, according to the subconcepts constituting the concept of resilience, resistance refers to the extent of withstanding an initial shock (Martin et al. 2016), which will allow us to resolve the limitations of previous studies.

Micro-, small-, and medium-sized enterprises (MSMEs), which are known to have suffered the greatest economic impact from COVID-19 in South Korea in 2020, are the objects of this study. Although South Korea had its first outbreak of COVID-19 before it became a global pandemic, the country’s projected change in GDP was between − 1.2% and − 2.5% in 2020, the lowest among the Organisation for Economic Co-operation and Development (OECD) countries (OECD 2020a). Specifically, the MSMEs of South Korea faced significant sales and employment setbacks due to the COVID-19 shock (OECD 2020b). This study used local administrative units—the basis for MSME economic activity—as the units of analysis to consider the composition and differences of regional economic structures.

This study estimated the regional resistance of MSMEs against COVID-19 shocks and identified the regional economic structure that influences resistance. The regional resistance of MSMEs was estimated by comparing their sales revenues recorded in early 2020, during the first COVID-19 outbreaks (weeks 6–9), to the same weeks of 2019. Robust regression and spatial regression analyses were then performed to identify the regional economic structure that affects resistance. The regional economic structure variables were the agglomeration of manufacturing and service industries and regional industrial diversity.

2 Literature Review

The economic impact of COVID-19 has been frequently described in terms of decreasing GDP. Its negative effect on global GDP is projected to approximately range between 2.3 and 4.8% (ADB 2020). The United States and major European countries are projected to experience a 12% drop in GDP, while Japan and South Korea are expected to experience a 6.5% and 4.1% decrease in their GDP, respectively (IMF 2020). The negative economic impact of COVID-19 has been generally found in MSMEs (Bartik et al. 2020; Shafi et al. 2020; Fairlie and Fossen 2021). Compared to large enterprises, MSMEs tend to be financially vulnerable, as their size and resources are insufficient to withstand shocks (Asgary et al. 2020; Bartik et al. 2020). Safety measures such as lockdowns and social distancing are necessary to combat infectious diseases, but they curtail individual consumption (Martin et al. 2020; Hong and Choi 2021). Economic losses and impacts are reflected unevenly across regions (Napierała et al. 2020).

A few studies on the impact of the COVID-19 crisis on MSMEs have focused on the effect of early lockdowns and social distancing measures. Studies generally examined their effect on the industries that have experienced revenue losses. Fairlie and Fossen (2021) found that shutdowns of California in the second quarter of 2020 resulted in an average loss of 17% in sales, with industries such as accommodation being heavily affected by mandatory lockdowns, recording 91% losses of sales. In particular, they found a relationship between COVID-19 patients and sales losses in essential business types and person-to-person contact industries. Motoyama (2020) analyzed cities with industrial structures that may be vulnerable to COVID-19. The study set 2019 figures as baseline values for simulation to study revenue changes and the correlation between population density, land area, and income. Lu et al. (2021) examined the impact of COVID-19 on small- and medium-sized enterprises and their survival. The study found that in 2020, the level of effects differed with respect to industry, and the measures undertaken by the government to control the pandemic are the most important factors for their survival. Prior studies on COVID-19 and its economic impact on MSMEs primarily conducted exploratory analyses on the changes observed in sales. In particular, they focused on the early shocks of COVID-19 and compared them to their corresponding 2019 figures. These studies explained the relationship between the various aspects of the pandemic and MSMEs’ vulnerable characteristics.

Resistance, the main object of this study, refers to the intensity of the reaction to a shock (Martin 2012; Kim and Lim 2016; Martin et al. 2016). The COVID-19 pandemic has acted as an external shock to MSMEs. Economic actors react differently to unexpected external shocks based on their resistance (Martin et al. 2016). Regional resistance, in particular, is affected by the regional economic structure, resources, capabilities, competencies, and local institutions. Studies on regional resistance against recessions have focused on the 2008 global financial crisis with an emphasis on industrial structures. Recent discussions have focused on the diversity in industrial structures (Di Caro 2017; Sedita et al. 2017), the composition of industries (Martin and Gardiner 2019; Giannakis and Bruggeman 2020), and agglomeration economies (Fratesi and Perucca 2018).

3 Materials and Methods

We chose the basic levels of South Korea’s local governments—si (city), gun (county), and gu (district)Footnote 1—as the units of analysis. Yeongyang-gun, Sinan-gun, Ulleung-gun, and Gunwi-gun were excluded from the sample because of missing data points. The dependent variable, MSME resistance, was obtained using Eq. 1, following the method introduced by Martin and Gardiner (2019). We compared the average MSME sales revenue recorded in weeks 6 to 9 of 2020 during the COVID-19 outbreak to that of weeks 6 to 9 of 2019.

$${\mathrm{RESISTANCE}}_{r}^{t,t-1}=\frac{\Delta {Y}_{r}-\Delta E({Y}_{r})}{|\Delta E\left({Y}_{r}\right)|},$$
(1)

where \(\Delta E\left({Y}_{r}\right)\) is the expected change in sales (\(Y\)) in region \(r\) during the economic recession of 2020 (\(t\)) compared to sales of 2019 (\(t-1\)), given by \(\Delta E\left({Y}_{r}\right)=\left(\frac{{Y}_{N}^{t}-{Y}_{N}^{t-1}}{{Y}_{N}^{t-1}}\right)\times {Y}_{r}^{t-1}\).

The MSMEs included in this study are defined by Article 2 of the Framework Act on Small and Medium Enterprises, which includes mining, manufacturing, construction, and transportation enterprises with fewer than 10 full-time workers and enterprises in other industries with fewer than 5 full-time workers. The model was constructed using the time periods when MSME sales revenue dropped significantly due to COVID-19 shocks. According to Kim (2021), South Korea experienced three waves during the COVID-19 outbreak. Figure 1 shows the percent changes of MSME average weekly sales revenue in 2020 compared to 2019. The drastic decrease in the four weeks between week 6 and week 9 can be attributed to the COVID-19 shocks observed after its first outbreak. Hence, we considered the resistance against the shocks observed during these four weeks as the dependent variable.

Fig. 1
figure 1

Percentage change of average weekly sales revenue in micro-, small-, and medium-sized enterprises (MSMEs) in South Korea in 2020 compared to 2019. Week 2: 6–12 January; week 3: 13–19 January; week 4: 20–26 January; … week 53: 28 December 2020–3 January 2021

The variables and data used to construct the model are listed in Table 1. For each time period, we used the estimated value of RESISTANCE from Eq. 1. The independent variables consist of the confirmed COVID-19 cases recorded during the first COVID-19 outbreak, the location quotients of manufacturing (LQ_MA), the location quotients of the service sector (LQ_SE), industrial diversity (DIV), population density (POPDEN), and the number of people employed per thousand persons (EMP). We included the LQ of manufacturing, wholesale, and retail industries as agglomeration economy variables, suggesting that these variables are the regional economy’s structural factors that affect shock resistance. Diversity is measured by the entropy index. As we are interested in the resistance of MSME sales revenue against the COVID-19 shock, we included the confirmed COVID-19 cases, population density, and the number of people employed per thousand persons as the control variables. The data from 2020 were used for the confirmed COVID-19 cases and to estimate resistance, while the other variables were constructed using the data from 2018. Equation 2 represents the model used in this study.

$${\mathrm{RESISTANCE}}_{r}={\beta }_{0}+{\beta }_{1}{\mathrm{COVID}19}_{r}+{\beta }_{2}{\mathrm{LQ}\_\mathrm{MA}}_{r}+{\beta }_{3}{\mathrm{LQ}\_\mathrm{SE}}_{r}+{\beta }_{4}{\mathrm{DIV}}_{r}+{\beta }_{5}{\mathrm{POPDEN}}_{r}+{\beta }_{6}{\mathrm{EMP}}_{r}+{\epsilon }_{r}.$$
(2)

where \(r\) is region and \(\epsilon \) is error term.

Table 1 Variables and data

The model of Eq. 2 means that the resistance of MSMEs is treated at the level of individual jurisdiction (spatial unit): si (city), gun (county), and gu (district). Since COVID-19 exhibits a spatial clustering phenomenon (Andrews et al. 2021; Liu et al. 2021), the resistance of MSMEs within the COVID-19 period—the dependent variable of this study—may also have spatial correlation. If the dependent variable has spatial correlation, the linear regression model cannot be BLUE (Best Liear Unbiased Estimator) and the spatial regression model must be considered (Anselin 2002). Therefore, this study used spatial regression models expressed in Eqs. 3 and 4, based on Eq. 2. Two local jurisdictions, Jeju Island and Ulleung Island, were excluded due to the spatial adjacency problem in the spatial regression model. The spatial weight matrix used Queen’s method, which sets the weighted value based on neighboring units sharing any point among contiguity-based weighting calculation methods (Anselin 2005).

$$y=\rho {W}_{y}+X\beta +\epsilon , \epsilon \sim N(0,{\sigma }^{2}I)$$
(3)
$$y=X\beta +\epsilon , \epsilon =\lambda {W}_{\epsilon }+\mu $$
(4)

where \(y\): \(\mathrm{n}\times 1\) matrix of dependent variable (\({RESISTANCE}_{r}\)); \(X\): \(\mathrm{n}\times \mathrm{k}\) matrix of independent variable (\({COVID19}_{r}, {LQ\_MA}_{r}\), \({LQ\_SE}_{r}\), \({DIV}_{r}, {POPDEN}_{r}\), \({EMP}_{r}\)); \({W}_{y}\): \(\mathrm{n}\times \mathrm{n}\) spatial weight matrix; \(\rho \): spatial lag parameter; \(\beta \): \(k\times 1\) matrix of regression parameter; \(\epsilon \): matrix of error term; \(\lambda \): spatial error parameter; \({W}_{\epsilon }\): \(\mathrm{n}\times \mathrm{n}\) spatial weight matrix; and \(\mu \): matrix of error term.

We chose ordinary least squares (OLS), robust regression, and spatial regression as our methods of analysis. Multiple linear regression analysis generally follows OLS, but in the presence of outliers, the normality assumption may not be valid (Rousseeuw and Leroy 1987; Verardi and Croux 2009). The Shapiro-Wilk W test for normal data returned a W value of 0.988 (p = 0.067), which resolved concerns related to normality, and both the Breusch-Pagan and White tests returned values of 2.18 (p = 0.1398) and 23.24 (p = 0.6722), respectively, ruling out heteroskedasticity. The mean variance inflation factor (VIF) value was 1.47; hence, we did not detect multicollinearity. However, the diagnostic plot (Fig. 2) reveals the existence of vertical outliers and bad leverage. Resolving the outlier problem requires M-estimator, S-estimator, and MM-estimator, which maintain the breakdown point while improving the efficiency by performing a robust regression. Moran’s I value was used to confirm spatial correlation as a precondition, using the spatial regression model. For model specification of the spatial regression model, Lagrange multiplier (LM) value, robust LM-lag test, and robust LM-error test were used (Yang 2010).

Fig. 2
figure 2

Diagnostic plot of standardized robust residuals versus robust Mahalanobis distance

4 Findings

The summary statistics of the variables, including MSME resistance to the COVID-19 shock, and the spatial distribution of MSME regional COVID-19 resistance are presented in Table 2 and Fig. 3, respectively. Gwacheon-si, Gyeonggi Province had the lowest resistance, followed by Hamyang-gun, Gyeongsangnam Province; Wanju-gun, Jeollabuk Province; Seocho-gu, Seoul Metropolitan Area; Jung-gu, Seoul Metropolitan Area; and Suseong-gu, Daegu Metropolitan Area. The areas with a high degree of resistance were Jangheung-gun, Jeollanam Province; Sunchang-gun, Jeollabuk Province; and Boseong-gun, Jeollanam Province. The spatial distribution of resistance shows some clusters. High levels of COVID-19 were found in Daegu Metropolitan City’s Nam-gu, Dalseo-gu, Dong-gu, and Suseong-gu, which were the epicenters of the first outbreak.

Table 2 Summary statistics
Fig. 3
figure 3

Spatial distribution of micro-, small-, and medium-sized enterprises’ (MSMEs) 2020 regional resistance to early COVID-19 in South Korea. Si = city; Gun = county; Gu = district

The results of the OLS, robust regression, and spatial regression analyses are presented in Table 3. The Hausman test finds the Robust-S estimation to be better than the OLS model (\({\chi }^{2}\) = 35.704, p = 0.000). However, between Robust-S and Robust-MM estimations, the Hausman test fails to provide a definitive answer (\({\chi }^{2}\)= 10.613, p = 0.1011). The Moran’s I value of resistance of MSMEs was 0.33 (p = 0.01), which was suitable for the application of the spatial regression model. The value of LM was 21.82 (p = 0.01) in the spatial lag model and 15.21 (p = 0.01) in the spatial error model. However, the robust LM test value of the spatial error model was not statistically significant, and the robust LM test value of the spatial lag model was 7.95 (p = 0.01). Therefore, the spatial lag model was applied as the spatial regression model.

Table 3 Ordinary least squares (OLS), robust, and spatial analysis results

In all models, MSME regional resistance against the COVID-19 shock falls as the number of confirmed COVID-19 cases increases. Although there is some difference from the coefficient value of Robust-S and Robust-MM, it was still statistically significant, as shown in the spatial lag model considering the spatial effect. The confirmation of new cases and the revelation of their movement patterns, often covering long distances, instigated a fear of infection and social contact, resulting in individuals resolving to pursue home isolation, which reduced their consumption and caused an early shock. Hence, MSME resistance against shocks is lower in areas with a high number of confirmed COVID-19 cases.

In Robust-S and Robust-MM models, agglomeration and diversity, that is the factors that affect resistance during a recession, affect resistance differently. Agglomeration had statistically insignificant effects on manufacturing and was found to be positively related to resistance, while agglomeration in wholesale and retail had a negative relationship with resistance. The areas with a high concentration of wholesale and retail industries directly related to consumption showed low resistance. Wholesale and retail industries were directly related to consumption. The low resistance in areas where these industries are concentrated could explain the possibility of a complex relationship between the initial lockdown and social distancing due to COVID-19 and the decrease in consumption brought about by individuals’ psychological behavior. However, in the spatial lag model that considers the spatial effect on resistance of MSMEs, none of the agglomeration variables were statistically significant. This means that if the situation where COVID-19 affects the resistance of MSMEs in adjacent spaces through epidemiological transmission is considered, the regional industrial agglomeration cannot specify the relationship with the resistance of MSMEs in the region. The Robust-S and Robust-MM models assume that the phenomenon of disaster is independent per a spatial unit. This means that when the spatial dependences of disaster were not considered, the agglomeration of the industry related to the object of resistance can affect it. However, disasters such as COVID-19 that cluster with surrounding areas and show their socioeconomic impacts mean that the relationship with the surrounding areas is more important than the industrial agglomeration characteristics of the area.

Similar to previous studies (Di Caro 2017; Sedita et al. 2017), diversity was found to be positively related to resistance in all models. High levels of diversity in regional industrial structures mean that other substitutable industries can absorb some of the effects. It appears that MSME resistance to COVID-19 shocks is affected by the diversity in regional industrial structures. The control variables population density, and the number of people employed per thousand individuals, are negatively related in all models. Population density indicates an area’s level of urbanization, and the number of people employed per thousand persons indicates the size of the area’s economy. The negative relationship between both the variables suggests that urbanized areas and those with large economies exhibit low resistance. For the Robust-S and Robust-MM models, the coefficient of the COVID-19 variable denotes a decrease in MSME resistance with an increase of one confirmed COVID-19 case. The ratio of the estimated coefficient value to the standard deviation of the resistance indicator is 5% to 10%, which implies that the initial shock of COVID-19 had a large impact.

5 Discussion

This study helps us understand the differences in the regional economic impact of disasters, for example, the unexpected and sudden impacts of infectious diseases such as COVID-19. Governments in many countries have used social distancing and lockdown policies (Castex et al. 2021; Glogowsky et al. 2021; Woskie et al. 2021) to contain the spread of the disease. The initial impact of the unexpected shocks and unprecedented policies appears to have been large, and this study focuses on the resistance response to this shock. The resistance of MSMEs associated with the low- and middle-income classes has been negatively affected due to COVID-19, such that many governments are implementing financial subsidies for MSMEs. People with low socioeconomic status are less resilient and more likely to fall back into poverty (Ur Rahman et al. 2021). Low- and middle-income class people working in MSMEs are also likely to be less resilient and to have a worse socioeconomic status than in their pre-COVID-19 situation. In order to respond to the shock experienced by low socioeconomic status groups of MSMEs, policies that can increase income in the short term are needed (Shafi et al. 2021). Subsidies for employment and business are also required to minimize vulnerabilities (Ur Rahman et al. 2021). As well as for MSMEs to respond more sustainably, it is necessary to ensure diversity in the local industrial structure, as indicated by the results of this study.

This study also showed that the role of industrial diversity, which is known to have positive externalities in regional economic growth, plays a positive role in the regional economic resistance to the unexpected shock of COVID-19. The positive role of industrial diversity in regional economic growth was explained by the knowledge externalities of the interactions performed among workers with diverse industries, jobs, and skills in the region (Boschma 2005; Frenken et al. 2007; Boschma et al. 2014). Such industrial diversity plays a positive role not only in growth, but also in answer to shock. Diversity provides redundancy to the economic system through overlapping and complementary functions and plays a positive role in absorbing the impact of a shock (Martin and Sunley 2015). A more diverse industry structure showed better employment performance with respect to the impact of the Great Recession (Brown and Greenbaum 2017), and economic diversity was an effective part of the resilience related to regional recovery and growth, even during and after natural hazard-related disasters such as floods. The fact that regional industrial diversity plays a positive role in the resistance of MSMEs to the shock of COVID-19 is attributed to the existence of multiple networks between local MSMEs and other industries and redundancy, which perform a complementary function in mitigating the impact of shocks.

This study has two main implications. First, industrial policies that encourage diversity are required to reduce the regional impact of unexpected shocks on MSMEs. Although the confirmed cases directly related to COVID-19 shocks caused related changes and had a negative effect on resistance, the magnitude of effect was less than diversity. Hence, diversity in a region’s industrial structure is more important for resistance than the direct effects of a variable related to shocks. Studies on recessions have shown that a region’s industrial structure diversity has positive effects on resistance and recovery, which was also found to be true in the context of COVID-19. However, MSMEs are structurally vulnerable to unexpected shocks and various types of recession. Our results emphasize the need for mid- to long-term plans that encourage regional diversity, rather than the implementation of plans to support MSMEs directly, to improve their regional resistance.

Second, in the expectation that the future will entail various threats with highly uncertain effects, alternative policies are needed to establish a more resilient structure by reducing the impact of early shocks. The impact of early shocks is driven by panic created by misinformation and related behaviors, which have various ripple effects. During the early shock of COVID-19, there was a fear of infection, an increase in misinformation, a decrease in consumption, and consequently, there were decreases in MSME revenues. Our results show that MSME revenues did not change as dramatically during the later COVID-19 outbreaks observed in 2020 compared to the early shock. Because MSMEs have weak financial structures and resources, their survival may depend on the length of the early shock. Thus, this study suggests that for the shocks that entail uncertainty and unexpectedness, such as COVID-19, reducing the early shock duration and providing correct information are essential to improving resistance.

6 Conclusion

This study examined the South Korean MSMEs most affected by COVID-19 by estimating their regional resistance to early COVID-19 shocks and identifying the regional economic factors affecting regional resistance. The spatial distribution of MSME regional resistance was partly concentrated in the southwest and northeast areas of South Korea, and spatial correlations of MSME regional resistance were identified. The increase in the number of confirmed COVID-19 cases showed a negative relationship with regional resistance. Among the regional economic factors that affect MSME regional resistance, diversity was the most positively related, while the agglomerations of the wholesale and retail industry sectors, consisting of many MSMEs, were negatively related. The agglomeration of the manufacturing industry was positively related to regional resistance but was not statistically significant. Densely populated urban area and the number of confirmed COVID-19 cases were negatively related to MSME regional resistance.

This study has limitations. While it addresses the outliers by utilizing robust regression, this is a limitation as the model does not consider business fluctuations, the second and third waves of the COVID-19 outbreaks through time series data, and the relationships of MSMEs with large enterprises. To include business fluctuations, the data on MSME revenues recorded before 2019 are required. However, these data are limited, and were not included in the model. Government relief programs and safety measures have changed over time, with the implementation of varying degrees of social distancing and lockdown measures. These factors must be controlled to include shocks from the second and third pandemic waves in 2020. However, because it is difficult to compare the different relief policies and measures intended to control infectious diseases, they were excluded from this study. Finally, it is important to study the change in MSME revenues in relation to not only individual consumers, but also large enterprises. However, because of insufficient data to study the relationships between MSMEs and large enterprises, they were excluded from this study. Therefore, the results of this study should be interpreted such that there are relevant linkages between resistance, COVID-19, and the related economic variables. However, no causal relationship was determined.