Introduction

The onslaught of the coronavirus disease 2019 (COVID-19) remains severe as of June 20211. Many countries have implemented lockdowns to prevent the spread of infection. While these lockdowns have been effective in reducing the spread of infection, many negative psychological effects have been reported, including a high prevalence of psychiatric symptoms, such as depression and anxiety2,3,4,5,6,7.

Previous studies reporting the negative effects during lockdown have been limited to cross-sectional surveys3,4,6 and follow-up studies under single lockdown situations5,7. There have been no longitudinal data demonstrating the effects of repeated lockdowns. Therefore, the question of whether multiple lockdowns have a negative influence on symptoms or whether the influence can be reduced as people become accustomed to lockdown circumstances remains unanswered8. Considering that lockdowns have a significant impact on people's lives and economic activities9, it can be assumed that repeated lockdowns have different effects on people's mental health and interpersonal interaction styles than single lockdowns. Since repeated lockdown policies have already been implemented in many parts of the world and may continue to be implemented in the future due to new outbreaks, it is essential to elucidate the influence of repeated lockdowns on mental health and interpersonal interactions.

Japan is one country that has implemented repeated lockdowns. The Japanese government declared a state of emergency on 7 April 2020 and 8 January 2021 and implemented a "mild lockdown”6, an unenforceable request for self-restraint (Fig. 1). During the first mild lockdown, a significant increase in psychological distress and depression was reported despite the lack of legal enforcement of self-restraint6,10,11, while the influence of the second mild lockdown on those participants is unclear.

Figure 1
figure 1

Status of COVID-19 infection in Japan from March 2020 to April 2021 and the study timeline.

In this lockdown context in Japan, longitudinal studies of mental health characteristics are critical to understanding the influence of repeated lockdowns. Such analyses will provide important fundamental data for improving lockdown approaches that still permit maximum infection prevention and maintain mental health.

The influence of lockdown is also particularly prevalent among certain groups, such as women, the young, the socially disadvantaged in terms of income, and those with a history of mental illness6,12,13. Because the effects of repeated lockdowns may be specific to these groups, a detailed examination of these associations would provide information regarding what groups require particular attention.

In this study, we examined the influence of repeated lockdowns on mental health and social isolation by conducting a large-scale panel survey of people in urban areas under two lockdowns due to declared states of emergency in Japan.

Results

Participant characteristics

The mean age of those who participated in both surveys (N = 7893, response rate 69.6%, Table 1) was 49.6 years (SD = 13.7), with the 30–49 years age group being the most highly represented (41.8%). Of the respondents, 65.6% were married, and 68.2% were employed. Compared to those who participated in the first wave only (N = 3440), those who participated in both were older and had a larger proportion of men and married people (all p-values < 0.001; Supplementary Table S1). In addition, those who participated in both experienced greater loneliness and smaller social networks, while fewer people exceeded the cut-off points for depression (Patient Health Questionnaire-9 [PHQ9] ≥ 10) and psychological stress (Kessler Psychological Distress Scale [K6] ≥ 5) and reported suicidal ideation in the first survey (all p-values < 0.001; Supplementary Table S1).

Table 1 Demographic characteristics of the participants.

Differences in living conditions and self-restraint between the first and second waves

Compared to the first wave, the frequency of experiencing household financial difficulties, living difficulties, and academic or occupational difficulties in the second wave decreased significantly (all p-values < 0.001; Supplementary Table S2). In addition, compared to the first wave, consumption activity index improved (wave 1: 82.7; wave 2: 94.4; Supplementary Fig. S1) and public transportation use tended to increase in the second wave (wave 1: − 43 to − 42%; wave 2: − 29 to − 25%; Supplementary Fig. S2).

Mental health under repeated mild lockdowns in the whole sample

Depression, suicidal ideation, psychological distress, and physical symptoms were significantly lower in the second than the first wave (all p-values < 0.001; Table 2). Alternatively, in the second wave, loneliness increased significantly (b = 0.24, 95% confidence interval [CI] 0.15 to 0.33, p < 0.001) and social network decreased significantly (b =  − 0.43, 95% CI − 0.54 to − 0.33, p < 0.001). The percentage of patients who exceeded the cut-off score for Generalised Anxiety Disorder scale-7 (GAD), measured only in the second wave, was 8.5% (Table 2).

Table 2 Changes in outcome variables over waves 1 and 2 with odds ratios, betas, and 95% confidence interval.

Mental health under repeated mild lockdowns with a focus on subgroups

The results of the interactions are mainly described, and details of all analyses are provided in Supplement A in the Supplementary Materials.

Estimated rate of depression

Both men and women displayed a reduction in depression rate in the second wave (Table 3). However, this reduction was significantly smaller among women than men (b = 0.18, 95% CI 0.52 to 0.32, p < 0.001). Women also had a significantly higher estimated rate of depression in both the first and second waves than men (wave 1: odds ratio [OR] = 1.30, 95% CI 1.15 to 1.46, p < 0.001; wave 2: OR = 1.56, 95% CI 1.36 to 1.78, p < 0.001).

Table 3 Changes in outcome variables by age group and gender.

The highest estimated rate of depression from the first to the second wave was observed among the 18–29 years age group, and only this group demonstrated no reduction in the rate of depression (OR = 0.96, 95% CI 0.73 to 1.27, p = 0.778; Fig. 2; Table 3).

Figure 2
figure 2

Estimated rate of depression for each age group in wave 1–2. CI confidence interval, OR odds ratio.

In addition, those with a history of treatment for mental illness (current, past, and both) and those with an income of < ¥2 million had particularly high estimated rates of depression (Table 4, Supplementary Table S3, Supplement A).

Table 4 Changes in outcome variables by income group.

Suicidal ideation

Although suicidal ideation was reduced in the second wave for both men and women (Table 3), the degree of reduction in suicidal ideation was significantly smaller among women than men (b = 0.27, 95% CI 0.14 to 0.39, p < 0.001; Fig. 3). Women also reported significantly more suicidal ideation than men in the second wave (OR = 1.29, 95% CI 1.15 to 1.46, p < 0.001), although this difference was not observed in the first wave (OR = 0.99, 95% CI 0.89 to 1.11, p = 0.907).

Figure 3
figure 3

Suicidal ideation by gender in wave 1–2. CI confidence interval, OR odds ratio.

Only the 18–29 years age group displayed no reduction in suicidal ideation from the first to the second wave (OR = 0.99, 95% CI 0.76 to 1.29, p = 0.947; Table 3).

In addition, those with a history of treatment for mental illness and those with an income of < ¥2 million had a particularly high rate of suicidal ideation (Table 4; Supplementary Table S3, Supplement A).

Psychological distress

All age groups displayed lower stress in the second wave (Table 3), but this reduction was significantly smaller among the 18–29 and 30–49 years age groups than the 65 years age group (18–29 years: b = 0.45, 95% CI 0.21 to 0.69, p < 0.001; 30–49 years: b = 0.32, 95% CI 0.14 to 0.49, p < 0.001). Among those aged ≥ 65 years, stress was significantly lower than in the other age groups in both waves (all p-values < 0.001).

Moreover, stress was higher among women, those with a history of treatment for mental illness, and those with an income of < ¥2 million (Tables 3 and 4; Supplementary Table S3, Supplement A).

Somatic symptoms

Although physical symptoms decreased significantly from the first to the second wave in all age groups (Table 3), only the 18–29 years age group demonstrated a significantly smaller decrease in physical symptoms than those aged ≥ 65 years (b = 0.88, 95% CI 0.45 to 1.31, p < 0.001). Additionally, in both waves, physical symptoms were reported significantly less frequently among those aged ≥ 65 years than other age groups (all p-values < 0.001).

Further, physical symptoms decreased significantly from the first to the second wave for all income groups (Table 4), but compared with the group with an income of ≥ ¥8 million, the group with an income of < ¥2 million in the second wave displayed particularly high physical symptoms (b = 1.88, 95% CI 1.20 to 2.56, p < 0.001).

In addition, women and those with a history of treatment for mental illness showed higher physical symptoms (Table 3; Supplementary Table S3, Supplement A).

Loneliness

Among men, loneliness did not vary between waves, (b = 0.09, 95% CI − 0.03 to 0.21, p = 0.162), whereas for women, loneliness increased significantly (b = 0.42, 95% CI 0.07 to 0.55, p < 0.001) (Table 3). Loneliness was significantly higher among men than women during both waves (wave 1: b = 0.76, 95% CI 0.51 to 1.01, p < 0.001; wave 2: b = 0.43, 95% CI 0.17 to 0.69, p = 0.001).

In addition, loneliness was higher among those aged < 65 years, those with a history of treatment for mental illness, and those with an income of < ¥2 million (Tables 3 and 4; Supplementary Table S3, Supplement A).

Social network

From the first wave to the second wave, both men and women significantly decreased their social networks (Table 3), but the degree of decrease was significantly greater for women than for men (b = −0.43, 95% CI − 0.54 to − 0.33, p < 0.001). Additionally, in both waves, men had significantly lower social networks than women (wave 1: b = −1.36, 95% CI − 1.62 to − 1.09, p < 0.001; wave 2: b = −1.14, 95% CI − 1.40 to − 0.87, p < 0.001).

Additionally, those aged 30–64 years, those with a history of treatment for mental illness, and those with an income of < ¥2 million had lower social networks (Tables 3 and 4; Supplementary Table S3, Supplement A).

Exhaustive interaction of factors associated with depression

The final convergence results of the extracted cluster structures are shown in Fig. 4 and Supplementary Table S4.

Figure 4
figure 4

Visualization of an exhaustive interaction structure of factors associated with depression. The rows represent the data of the participants. The columns represent depression, suicidal ideation, and their associated factors. The black lines show the clusters of participants and factors, respectively. The colours in the cells indicate the magnitude of the z-values or the demographic attributes (see Box at the bottom of the figure). The appended (a), (c), (e), and (g) indicate clusters with significantly higher or lower loneliness and social networks, and (b), (d), (f), and (h) with significantly higher or lower depression and suicidal ideation. M, million; SI, suicidal ideation; SocNetw, social network; Treatment, history of treatment for mental illness; W, wave of the survey.

Cluster 9 (CL9, n = 385) and CL13 (n = 136), with significantly higher loneliness and the lowest social network (Fig. 4a), had the highest estimated rates of depression in both waves (CL9: 47.8–54.8%; CL13: 35.3–42.6%), and both CLs had the highest number of people reporting suicidal ideation (CL9: 48.6–50.6%; CL13: 40.4–47.8%) (Fig. 4b). Similarly, CL3 (n = 718), the next highest in loneliness (Fig. 4c), also showed higher estimated rates of depression (28.6–37.6%) and suicidal ideation (30.4–36.8%) in both waves (Fig. 4d).

Alternatively, CL11 (n = 281), which had the lowest loneliness and also the most preserved social network (Fig. 4e), showed the lowest estimated rates of depression (0.4–0.7%) and suicidal ideation (both 0.4%) among all CLs in both waves (Fig. 4f). Similarly, CL7 (n = 637) and CL5 (n = 794), who were the next least lonely and maintained social networks (Fig. 4g), had lower estimated rates of depression (CL7: 1.9–3.3%; CL5: 2.0–2.3%) and suicidal ideation (CL7: 2.8–3.9%; CL5: 2.0–2.9%) (Fig. 4h).

Discussion

The results of this study indicate that the overall psychological and physical symptoms decreased during the second mild lockdown compared to the first. During the first mild lockdown due to the declaration of the first-ever state of emergency in Japan, there were major economic and social disruption impacts14, resulting in a significant increase in depression and stress. Alternatively, difficulties in home, work, and school life were reduced in the second wave compared to the first. Further, in the second wave, both consumption activity and mobility tended to increase. Therefore, it is likely that these symptoms were reduced to some extent during the second wave, due to adaptation to various environments and changes caused by the COVID-19 pandemic.

However, at both time points of the study, social network levels indicated social isolation (Lubben Social Network Scale [LSNS-6] < 12)15, which was even more pronounced during the second mild lockdown. Furthermore, loneliness scores increased in the second wave compared to the first wave. These considerations indicate that repeated mild lockdowns cause cumulative changes in interpersonal interactions and that subjective loneliness is maintained and exacerbated.

A striking result when focusing on subgroups was that younger people (aged 18–29 years) displayed the highest levels of depression and suicidal ideation, and unlike other groups, these indicators did not decrease during the second wave. Additionally, the younger age group demonstrated the smallest decrease in stress and physical symptoms in the second wave of the survey. Women also indicated higher levels of depression, suicidal ideation, and increased loneliness compared to men. These results are consistent with previous reports of deteriorating mental health among young people during the pandemic5,7,16, indicating that younger people are more vulnerable to the negative impacts of repeated mild lockdowns. Further, considering that the estimated rate of depression among younger people in this study (22.8–24%) was significantly greater than that of the general population in Japan in the past (7.9%)17, the present results may indicate that young people experience difficulties with their current situation (e.g., various obstacles related to school and work, financial difficulties due to inability to work), and concerns about the future (e.g., uncertain career paths and employment status) during lockdown situations. Given these considerations, immediate support for young people is highly desirable during repeated lockdown situations. For example, support is needed from educational institutions, companies, and the government by guaranteeing educational and employment opportunities, providing mental health care, and providing information to assuage their concerns.

In addition, those with a history of treatment for mental illness and those who were socially disadvantaged (especially those with incomes < ¥2 million) had higher levels of depression, suicidal ideation, stress, physical symptoms, and loneliness, and smaller social networks at both waves. These findings are consistent with that of previous studies18,19 and indicate that these populations are also more susceptible to negative effects of repeated lockdowns. The COVID-19 pandemic has also been reported to cause various economic crises, resulting in lower wages, cutbacks, and general job losses in various industries9, making life increasingly difficult. Therefore, during repeated lockdown situations, a detailed understanding of people's difficulties with these characteristics is important, and additional social support from social workers and local governments should be considered.

Furthermore, the comprehensive extraction of interactions among variables associated with depression demonstrated that the intensity of loneliness and the size of social networks were most closely related to mental health during both mild lockdown situations. Loneliness and social networks have been considered risk factors associated with mental and physical health20,21 and have been identified as important factors affecting mental health during the COVID-19 pandemic4,22. Considering these previous studies and our results indicating that repeated mild lockdowns cumulatively reduced social networks and maintained or exacerbated loneliness, we believe that the health status of people who display reduced social networks and high loneliness requires special attention in these situations.

Limitations

In this study, causal factors, such as which approach mitigates or exacerbates the effects of repeated mild lockdowns, have not yet been examined. Therefore, it is necessary to examine the effects of psychosocial variables, such as coping strategies and lifestyles, in the future.

Additionally, since this study used only an online survey method, the mental health of groups without online access has not yet been examined. Further, since mental health indicators were based on self-reports rather than clinical diagnoses, they may not necessarily correspond to objective assessments by mental health professionals. Therefore, it is necessary to widely use non-online methods in the future to accumulate findings that can be generalised to a larger population.

Finally, those who did not participate in the second wave of the survey had worse mental health than those who participated in both waves. Therefore, it should be noted that the present results using only those who participated in both waves may underestimate the impact of repeated lockdowns. Additionally, although the results of previous studies indicate that mental health worsened with the first mild lockdown6,23, we did not obtain pre-pandemic data for the participants in our study. Therefore, the exact transition in mental health needs to be judged comprehensively, taking into account the results of future studies.

Conclusion

During the two rounds of mild lockdown in Japan, depression, suicidal ideation, stress, and physical symptoms decreased overall in the second compared to the first wave, while loneliness increased and social networks decreased. These results suggest that repeated mild lockdowns exacerbates social isolation and has a cumulative negative effect on loneliness. In addition, depression and suicidal ideation were highest among younger individuals, and as a striking result, they did not decrease in the second wave only in the younger age group. Furthermore, women, those with a history of treatment for mental illness, and those from socially disadvantaged backgrounds displayed poorer mental and physical health in both waves. These results indicate that these populations are particularly vulnerable to the negative effects of repeated mild lockdowns. Additionally, loneliness and social networks, in particular, were closely related to mental health at both time points. Given the potential for social isolation and loneliness to be exacerbated by repeated mild lockdowns, these results suggest that during a pandemic situation, people who exhibit these characteristics may require special attention. In the future, it will be essential to conduct continuous research on the health status of potentially vulnerable populations using a combination of various measures to elucidate protective and risk factors and to develop support systems tailored to individuals’ difficulties.

Methods

Study design, participants, and data collection

This longitudinal study was conducted online between 11 and 12 May 2020 (wave 1) and between 24 February and 1 March 2021 (wave 2). Both online surveys were conducted just before the emergency declaration was phased out, and each was set up to assess the psychological impact of a mild lockdown of approximately 1 month on participants (Fig. 1). Participants were recruited by Macromill, Inc. (Tokyo, Japan), a global marketing research company. In the first wave of the survey, approximately 80,000 people registered with Macromill, Inc. were recruited via email (the target sample was 11,000 people). A total of 11,333 people participated in the initial recruitment (the first wave), and 7893 of them participated in the follow-up recruitment (the second wave). Participants were recruited only from the seven prefectures where the emergency declaration was first applied (Tokyo, Kanagawa, Osaka, Saitama, Chiba, Hyogo, and Fukuoka prefectures) to sensitively detect the impact of mild lockdown. Since these cities have large populations and a large number of reported cases, it was assumed that they would sensitively reflect the impact of mild lockdown. The exclusion criteria for participants were: aged < 18 years, high school students, and living outside the seven prefectures.

The amount of data collected in each prefecture in the first wave was determined according to the ratio of the number of people living in each province to the total population of the seven prefectures (e.g., in the case of Tokyo, 11,428,937 (the population of Tokyo)/46,548,456 (the population of the seven prefectures) = 24.6%), Tokyo (n = 2783, 24.6%), Kanagawa (n = 1863, 16.4%), Osaka (n = 1794, 15.8%), Saitama (n = 1484, 13.1%), Chiba Prefecture (n = 1263, 11.1%), Hyogo (n = 1119, 9.9%), and Fukuoka (n = 1027, 9.1%). The data from the first wave of the survey in this study are part of the dataset used in previous publications6,11, and details of this data can be found in these papers.

A link to an online platform to participate in the survey was distributed to those who wished to participate in the study, and the survey was conducted online. All participants completed the anonymous survey voluntarily and provided informed consent online. The online survey form was designed so that all items were required to be answered before proceeding to the next step; thus, there were no missing data in this study. The survey procedures were clearly explained, and participants could interrupt or terminate the survey at any time without providing an explanation. This study was approved by the Research Ethics Committee of the Graduate School of Social and Industrial Science and Technology, Tokushima University (approval number: 212) and has been performed in accordance with the ethical standards laid down in the 1964 Declaration of Helsinki and its later amendments.

Measures

Depressive symptoms were assessed using the Japanese version of the PHQ-9 24, which consists of nine items; a score of ≥ 10 was considered to indicate a high likelihood of depression24. To show the trend of the percentage of people with severe symptoms and to allow comparison with previous studies3,5, we used this as a cut-off point in this study.

Suicidal ideation was assessed using item 9 of the PHQ-9 and coded as a binary variable, with a response of “not at all” as “no suicidal ideation” and all other responses as “suicidal ideation”.

Psychological distress was measured using the Japanese version of the K625, which consists of six items. A threshold of five points, commonly used for screening mild to moderate mood/anxiety disorders, was adopted as the cut-off point for psychological distress26.

Somatic symptoms were measured using the Japanese version of the Somatic Symptom Scale-827, which consists of eight items.

Loneliness was measured using the 10-item Japanese version of the UCLA Loneliness Scale Version 328.

Social networks and the presence of social support were measured using the Japanese version of the abbreviated LSNS-615, which consists of six items. LSNS-6 scores < 12 indicate social isolation, and higher scores indicate a larger social network.

The Japanese version of the GAD-729, consisting of seven items, was used to assess symptoms of generalised anxiety disorder. Scores > 10 are considered to indicate a moderate level of anxiety30 and were used as a cut-off point in this study. Since the GAD-7 was used only in the second wave, it was used to describe the characteristics of participants in the second wave.

These self-report measures are more extensively described in Supplement B in the Supplementary Materials.

To account for differences in living conditions between waves 1 and 2, we asked respondents to rate the frequency of experiencing “household financial deterioration”, “living difficulties”, and “academic or occupational difficulties” from the start of mild lockdown to the time of the survey on a scale of 1 (not at all) to 7 (extremely), based on Yamamoto et al.6. Similarly, to examine differences in social self-restraint, consumption activity (nominal consumption activity index), and going-out behavior (public transportation congestion) in waves 1 and 2 were calculated based on basic data from the Bank of Japan31 and Google Community Mobility Reports32, respectively.

Further, to determine the influence of mild lockdowns on groups that have been previously identified as vulnerable to the effects of lockdown6,12,13, demographic attributes, such as age, gender, and income, were collected, as well as whether or not individuals were currently being treated for mental illness and whether or not they had a history of treatment for a mental illness.

Statistical analysis

Analyses were performed on a dataset of 7893 participants from both waves of the survey. First, comparisons were made using t-tests and chi-square tests for the data in the first wave for those who participated in both surveys (N = 7893) and those who only participated in the first wave survey (N = 3440).

Next, life and academic/occupational difficulties in the first and second waves were compared using t-tests.

For each outcome variable, a generalised estimating model (GEE) was constructed to test for changes in the variable between waves for the whole sample. This approach is applicable to binary variables and variables other than normally distributed ones and has excellent analytical precision, even for longitudinal data with high within-subject correlation33,34. In the GEE approach, binomial logit modelling was used for the binary outcome variables (PHQ-9 cut-off score, suicidal ideation, and K6 cut-off score), and linear Gaussian identity modelling was used for the continuous outcome variables (physical symptoms, loneliness, and social network).

Additionally, GEE models were used to examine differences in changes in outcome variables in subgroups focused on gender (men, women), age (18–29, 30–49, 50–64, ≥ 65 years), income (< ¥2 million (approximately £15,000), ¥2–4 million, ¥4–6 million, ¥6–8 million, ≥ ¥8 million), and treatment for mental illness (currently in treatment, past treatment, past treatment and still in treatment, no treatment). For analysis of the whole sample, binary outcome variables were assessed using a binomial logit GEE model, and continuous outcome variables were assessed using a linear Gaussian identity GEE model. Based on Ballinger33, the Wald test was used to examine the significance of the model. Significant interactions between each subgroup and wave were mainly reported. When analysing the subgroups of income groups, we excluded 1781 participants who answered that they did not know.

To visualise the interaction structure of the variables most closely associated with depression, we used nonparametric Bayesian co-clustering35 for gender, age, income, history of treatment for mental illness, depression, suicidal ideation, loneliness, and social networks. Since this method is capable of exhaustively extracting various clusters represented by multidimensional factors6,35, it is expected to find the most influential factors for depression. A total of 15,000 iterations based on the Bayesian optimization principle were performed, and the marginal log-likelihood (Bayesian factor) was calculated as an objective criterion for good classification. After confirming that the marginal log-likelihood was completely stable, the model with the highest marginal log-likelihood was finally adopted. This procedure allowed us to select a classification result with high objectivity and reproducibility.

The significance level for all tests was set at α = 0.05, two-tailed. The analysis was performed using Matlab R2017a (Mathworks Inc.) and SPSS (version 22.0; SPSS Japan Inc., Tokyo, Japan).

Ethical standard

This study was approved by the Research Ethics Committee of the Graduate School of Social and Industrial Science and Technology, Tokushima University (approval number: 212) and has been performed in accordance with the ethical standards laid down in the 1964 Declaration of Helsinki and its later amendments.