key: cord-0841950-06c00on8 authors: Li, Haizheng; Ma, Mingyu; Liu, Qinyi title: How the COVID-19 pandemic affects job sentiments of rural teachers date: 2022-02-01 journal: China Econ Rev DOI: 10.1016/j.chieco.2022.101759 sha: 00712fa575256cf8273221e5ca4594017267d2eb doc_id: 841950 cord_uid: 06c00on8 This study investigates how the COVID-19 pandemic has affected teachers' job-specific stresses and their enthusiasm for the teaching occupation. We use unique data from China that cover the periods before and after the start of the pandemic and apply various estimation methods. We find that, among rural young teachers, the pandemic has caused higher teaching stress and career development stress and has reduced passion toward the teaching occupation. We investigate the working channels of the pandemic, including changes in job-related activities and social network. After controlling for possible working channels, the COVID-19 pandemic still shows a strong direct impact on job sentiments. to the COVID-19 pandemic. For example, Qiu et al. (2020b) studies the impact on the general population; Kang et al. (2020) focuses on vulnerable groups such as health professionals; Aucejo et al. (2020) studies the influence on student experiences. However, little is known about how the pandemic influences specific job-related stresses and job perceptions. Because the COVID-19 pandemic may affect occupations differently, it is important to investigate its psychological impacts on certain highly affected occupations, especially those where work morale has a deep impact on society. 2 Teaching is considered one of the most stressful professions and requires intensive interactive instructions (see Collie et al., 2012; Avdiu and Nayyar, 2020) . Teachers with high levels of stress show a reduced sense of job satisfaction and a tendency to exit the teaching profession. More importantly, teachers' work attitudes affect education quality and student performance (Harris and Adams, 2007) . Therefore, the influences of the pandemic on teachers have potential consequences for the educational outcomes of the future generations. In this study, we investigate the impact of the COVID-19 pandemic on teachers' job sentiments. We focus on how the pandemic influences teachers' job-specific stresses, including teaching stress and career development stress, and their job satisfaction as measured by enthusiasm for the teaching occupation. Due to the COVID-19 pandemic, lots of sudden changes have occurred for teachers, such as prolonged school closures and distance teaching (Orlov et al., 2021) . A recent survey of teachers in the United States during August 2020 found that approximately 32% of the respondents reported low morale, and 47% considered making a major job-related change, with 17% saying they would completely change their career away from teaching. 3 Our data come from large-scale annual surveys from approximately 7,500 elementary school and middle school teachers in rural China. 4 The surveys were conducted both before and after the breakout of the pandemic and thus allow us to estimate the total changes in job sentiments attributed to COVID-19. We apply the cross-section estimation as well as the difference-in-differences (DD) method to identify the impact of the pandemic. We further test the robustness of our results with various samples and estimation techniques. For most studies about the psychological effect of the COVID-19 pandemic, the data were collected after the pandemic started, and they are mostly cross-sectional (e.g., Wang, et al., 2020a; Tan et al., 2020) . For example, Zhang and Ma (2020) , using one province's data from China, finds around 70% of their sample reported no increase in stress from work between January 28 and February 5, 2020. Some longitudinal studies exist, but they generally cover different time points during the pandemic, not before the pandemic. For example, Wang et al. (2020b) found stable stress levels of the general population and no significant temporal changes despite the sharp increase in the number of COVID-19 cases from January 30 th to March 11, 2020. Such data can help identify the change of the psychological effect at different times during the pandemic but not the overall before-after effect of the pandemic. In comparison with the literature, our data have some unique features: 1) it is a longitudinal survey conducted before the pandemic and during the pandemic; 2) the survey provides detailed assessments of the respondents' job sentiments; 3) the survey focuses on young rural teachers in China, and thus the samples are homogenous in representing a particular population. In addition, as a routine annual survey of a large national training program for rural teachers, the purpose of the surveys is not directly related to the pandemic, so the responses are less likely to be induced toward a particular direction by the survey questions. We propose a theoretical framework on how the pandemic affects job stress, and measure job stress and enthusiasm levels with both categories and scales. The results show that the pandemic significantly increases teaching stress and career development stress. Moreover, the pandemic reduces passion toward teaching. In addition, local pandemic severity has statistically significant effects on teachers' career development stress and job enthusiasm, but the magnitudes are very small and economically insignificant. This result indicates that studies based on cross-sectional data during the pandemic may only reflect a smaller portion of the total COVID-19 pandemic effect. We further investigate the working channels of the pandemic, including behavioral changes anxiety, stress and even panic among the public. In general, the level of perceived risk increases with three factors: how dreaded, uncontrollable, and fatal the risk is, how unfamiliar and unknown it is, and the level of personal and social exposure to the risk (Wong, 2008) . The COVID-19 pandemic is at a remarkably high level for all the above characteristics and thus directly affects an individual's stress. Moreover, such stress can cause an overall social amplification of risks when the information is transmitted between individuals via interpersonal networks. For example, Holmes et al. (2020) finds that daily COVID-19 related media exposure and conflicting information in media were associated with acute stress and depressive symptoms in the United States. Therefore, the overall impact of the pandemic may not be just related to local exposure. We develop a conceptual framework of the dynamics of individual stress based on Cowan et al. (2011) . We model the change in stress over time as: where total stress level is , and the first term ( ) represents the effect of the COVID-19 pandemic . Because the overall pandemic impact is expected to come from the overall pessimistic and depressed atmosphere as well as from local exposure, we specify the impact of the COVID-19 pandemic into two components: the overall effect and the effect due to local COVID-19 severity , We define ( , ) as the rate at which stress levels fall due to self-coping. The ability of coping varies across seasons, for example due to weather or holidays, and the seasonal factor is represented by . represents the individual traits and experiences that affect the individual's ability to adjust their stress level. An individual can also reduce stress by sharing it with colleagues, friends and family members. ( ) captures these buffering effects, where represents the relative strength of one's social network, and ∈ [0,1], where 0 means the social relationship is broken, and 1 represents a strong social relationship. In addition, stress spills over between persons. Suppose an individual interacts with n colleagues at work and m friends and family members; their stress level would be influenced by the stress level of those social contacts. In particular, ∑ Additionally, we model the effect of the pandemic on job passion/enthusiasm in a similar framework to job stress in Equation (1). In psychology, passion/enthusiasm is defined as a strong inclination toward an activity (such as work) that one loves and that is self-defining (Vallerand and Houlfort, 2003) . Apart from individual characteristics, social recognition also represents an important determinant of passion (Tóth-Király et al., 2019) . A person's passion towards work can be determined by their own "enhancing" capacity as well as by "recognitions" from their social network (just like self-coping and buffering effects for stress). Moreover, the spillover effect also exists for job passion as for job stress. The first COVID-19 case was found in Wuhan, Hubei province in China. January 19 is the 5 We could disaggregate the buffering effects for different members in the network, and then aggregate to get the total buffering effect. To simplify the argument, we adopt the current model structure. 6 For simplicity, we use the linear function form in our empirical specifications. For example, in equation (4), we specify ( , ) = + . In empirical model (5), the stress coping rate ( , ) is represented in the linear form of S + X . first day that COVID-19 cases were reported outside of Wuhan (Qiu et al., 2020a) . On January 23, the Chinese government imposed a lockdown measure on Wuhan. By January 30, almost all provinces implemented the Level 1 Response to Public Health Emergency, the highest response level. 7 Figure 1 shows the trend of the cumulative, existing, and new COVID-19 cases. The spike of cumulative cases occurred around mid-February and then flattened out afterwards. The implementation of strict quarantine measures in China has kept many people in isolation. All schools in China postponed the start of the spring semester. The government encouraged schools to provide online instruction to students. Since February 17, a national online learning platform has been operated by the Ministry of Education to provide educational materials for students at primary and secondary institutions. 8 About 85 percent of students and teachers use mobile devices, such as smart phones and tablets, for online education (Huang et al., 2020) . As early as the end of March, with COVID-19 under control, primary schools started to reopen in some less affected provinces. Wuhan city reopened from a citywide lockdown on April 8. 9 Other provinces began reopening primary schools between April and early June. Different policies for teaching during the pandemic brought drastic changes in teachers' work patterns. The pandemic has forced teachers to switch to online teaching and has led to many abrupt changes at work and in life. Online teaching posed new challenges for teachers, as most of them were unfamiliar with the online teaching tools. After schools reopened, teaching did not go back to normal due to the new requirements for social distancing and the new hybrid format. Moreover, teachers' administrative workload increased substantially due to the need to prevent the COVID-19 from spreading in school and amongst students. Further, rural teachers usually work in less-developed, remote areas, and online teaching is new especially for them. Therefore, they are more vulnerable to direct exposure of the pandemic's influence. Our data are collected from the routine survey of a large scale online annual training program for young teachers in rural China, the Young Teacher Empowerment Program (hereafter as "YTEP"). The YTEP was initiated in 2017 through the sponsorship of non-profit organizations, universities, and corporations in China. 10 It aims to help young rural teachers better fit into rural environments and improve their teaching skills and work morale. The YTEP is a year-long training program starting in September and ending in the following June, and it provides online training courses via a broadcast platform. Participants watch the program videos online, either live or recorded, via a computer or cell phone. YTEP participants are selected by the local government and their schools. There are two types of teachers: permanent teachers (regular teachers) and special-term teachers. Special-term teachers work via a national program in which college graduates are recruited to teach in rural areas for three years to improve education quality. 11 After three years of service, teachers who pass the assessment can become permanent teachers, or they can choose other jobs. 12 Teachers work in various types of rural schools, including 1) rural schoolhouses, usually located in remote rural areas, which have the smallest school size and only offer primary school level education, 2) village schools, larger than schoolhouses, and 3) rural district schools, the largest among all three school types, which may offer both primary and middle school education. Our data come from the annual surveys of the YTEP participants as a routine evaluation of the program. One survey was conducted for the Class of 2018-19 (starting in September 2018 and ending in June 2019, hereafter as "Class-2019"). The other survey was for the Class of 2019-20 (hereafter as "Class-2020"). All surveys were administered online via the administrative team of the YTEP. 13 The survey of Class-2019 was conducted at the end of the program in June 2019. For Class-2020, the administrative team conducted two surveys; Wave-1 was added in the middle of the program in January 2020, and Wave-2 was conducted regularly at the end of the program in June 2020. They sent the survey links to participants during live-class times, and the survey was usually live for about one week. During the period that the survey was live, the administration team sent 2-3 reminders to training participants to fill out the survey. As a result, the respondents in each survey are not selected in any non-random way. A total of 7,502 rural teachers participated in all three surveys. The Wave-1 survey of Class-2019 started on January 2, 2020, and was completed by January 20, before COVID-19 became public in China. 14 The public in China was not yet informed about the new coronavirus during the Wave-1 survey. For example, according to Fang et al. (2020) , in the epidemic center Wuhan city, on January 18, 2020, more than 10,000 families gathered for the annual Lunar New Year banquet. In June 2020, the pandemic was mostly under control in China but was spreading rapidly all over the world. Therefore, those who participated in Wave-2 in June were influenced by the pandemic for approximately half a year. Details of the three surveys are shown in Appendix Table A1 . The initial purpose of the surveys is not related to the pandemic but is instead just a program evaluation. There is no question related to the pandemic in the surveys. Therefore, the surveys could have the advantage of receiving more accurate responses because there are no hints toward the pandemic. Given the nature of the homogenous sample, i.e., young rural teachers, the data provide a unique opportunity to study the changes of job stress and job enthusiasm due to the pandemic. In order to represent more accurately the same population, we restrict our samples to permanent teachers and special-term teachers. Temporary teachers are excluded from the sample as their job attitudes can be very different. We also keep teachers aged 35 or below and with no more than 5 years' teaching experience to focus on relatively inexperienced young teachers, because, as age and work experience increase, people are much more capable in handling stresses. 15 The final sample size used in the analysis is 5,767 after eliminating those with incomplete information. Because one objective of the YTEP is to help develop better morale amongst rural teachers, one part of the surveys specifically assesses job-specific stresses and attitudes toward the teaching occupation. The related survey questions are listed in Table 1 . In the literature, different sources have been cited as causes of teacher stress, e.g., stress related to workload or related to students' behavior and discipline (Klassen and Chiu, 2010) . Teachers may have different concerns and stress at different stages of career development (Holmes, 2005) . We classify related survey questions into three aspects pertaining to job attitudes, each aspect consisting of two questions. They are: 1) Teaching stress, 2) Career development stress, and 3) Job passion/enthusiasm for the teaching occupation. The measure of teaching stress focuses on stress 14 There were 6 participants that responded after January 20, and they were deleted from the sample. 15 Approximately 3.5% of the sample are outside this range. from students, such as helping them graduate and maintain discipline. The measure of career development stress concerns career advancements, including receiving promotions and awards. These stresses represent different aspects of job-related pressure. Job passion/enthusiasm represents a teacher's job attitude and is measured by questions like "I will not feel tired of being a teacher." In general, passion is associated with determination, motivation, and a high degree of self-control. When work is highly valued, it will be internalized in the person's identity in an autonomous fashion, leading to a harmonious passion . Research shows that passion for work is positively related to work satisfaction (Carbonneau et al., 2008) and will lower turnover intentions . For the above survey questions, the choices are divided into different categories, such as "no," "unsure," and "yes," as listed in Appendix Table A2 . We first classify those choices into two categories to indicate the stress status for a particular measure: high (as 1) and low (as 0). For example, the survey questions related to teaching stress are "I feel great pressure to help students graduate and enroll in the next level of education" and "I feel great pressure to maintain students' discipline." A respondent is considered to have a high value of stress if he/she chooses "completely agree," "agree," or "mostly agree" among all seven possible choices listed in the survey. If the high value of stress occurs for either of the two questions, we define an individual as having high teaching stress with a value of 1, and 0 otherwise. The detailed classifications for various choice sets are presented in Appendix Table A2 . The advantage of this approach is that the categories were directly obtained from the survey responses. However, such a measure can only capture large changes across measured categories. Therefore, we assign a scale for each category of survey responses so that we can measure the continuous changes of job sentiments. Based on the different choice sets for the related questions, we construct scale measures by assigning a value of 0-10 for each choice. Appendix Table A2 shows the assigned values in more detail. For example, we can assign a value of 0 for "no", a value of 10 for "yes", and a value of 5 for "unsure". To ensure the scaling was assigned as representatively as possible, all 10 survey team members assigned a scale value to each choice independently in three rounds at different times during a period of one month. Their average values are then used in our analysis. To assess the degree of agreement among the raters, we test the interrater reliability based on the Fleiss' kappa statistic. The Fleiss' kappa statistic measures the degree of agreement attained J o u r n a l P r e -p r o o f in excess of what would be predicted by chance (Fleiss, 1971) . 16 Our overall Fleiss' kappa test statistic among all raters is 0.51, indicating a moderate agreement among raters. Given the interrater reliability, we study the job sentiments based on both the category measures and scale measures. Table 1 shows the distribution of job stress and enthusiasm in the three surveys. Because each measure of job attitudes is based on two questions in the survey, we take the average of their values. It seems that a lower percentage of individuals indicates high job stress in Wave-2 compared to Wave-1. However, comparing with Class-2019, Class-2020 displays higher proportions of job stress during the pandemic. Similar results are found based on the scale measures, i.e., job stress levels are higher in June 2020 compared to those in June 2019, and job enthusiasm is also much lower. Table 2 shows the summary statistics of individual and job characteristics across the three surveys. As expected, the statistics between Wave-1 and Wave-2 of Class-2020 are almost identical. The samples from Class-2019 and Class-2020 surveys are comparable, except some differences in married and ethnic groups, as well as in the share of special-term teachers. The largest difference occurs for regional distribution, with 68% of the Class-2019 sample from non-western regions compared to approximately 25% for Class-2020. For further analysis, to better interpret the regression results, we standardize the 0-10 scale assigned to each job sentiment measure. Following Autor (2003) , we standardize the scales based on the pre-pandemic sample, i.e., the sample of Wave-1 (January 2020), so that the scales are comparable before and during the pandemic. Therefore, all measures from different surveys are conditional on the same distribution and the corresponding z-values reflect the changes in percentile in the distribution of job sentiments. Based on Equations (1) and (2) Journal Pre-proof COVID-19 due to local severity. However, the limitation is that it can only capture the marginal effect of the local severity on job sentiments, not the total effect of the pandemic. Local severities vary substantially, as measured by the number of cumulative COVID-19 cases at the provincial level. As of June 5, 2020, when the Wave-2 survey was conducted, the number of cumulative cases confirmed was the lowest in Qinghai province (18 cases) and the highest in Hubei province (68,135 cases). 17 In addition to local cumulative cases as a measure of pandemic severity, other potential measures include the number of new cases or existing cases. However, as seen in Figure 1 , because of the dramatic restrictions adopted by the Chinese government to control the COVID-19 pandemic, these numbers became almost zero for almost all provinces when the Wave-2 survey was conducted. Additionally, the numbers of cases are too small compared to the population of the province, and thus their ratio is not used as a measure. A potential concern is that the cumulative cases may also capture provincial effects that affect teachers' job sentiments. For rural elementary and middle school teachers, their job sentiments may be related to the social status of the teaching occupation in the local area and to the competitiveness of the job in the local labor market. Therefore, their job sentiments are likely to be affected by different stages of economic development. Thus, we add a regional dummy in the model to distinguish the relatively less developed western region from other regions to control for such effects. The empirical model using only the cross-section sample to estimate the effect of the local pandemic severity is specified as follows: where represents measures of job stress or enthusiasm for individual i, is the number of cumulative confirmed COVID-19 cases in the province where individual i lives to capture local pandemic severity of the pandemic, and includes a set of individual characteristics and job characteristics, and is individual-specific error term. We use both measures for job stress and enthusiasm: (i) categorial measures (high or low) and (ii) standardized scale measures. The result is reported in Table 3 . For the stress measures, the local severity has no significant impact on teaching stress but increases career development stress based on the scale measure. However, an increase of local COVID-19 cumulative cases 17 We also ran the regressions by excluding Hubei province, which includes Wuhan, due to its extremely large number of the COVID-19 cases, and the results are similar. J o u r n a l P r e -p r o o f has a statistically significant effect on lowering job enthusiasm based on both categorical and scale measures, but the marginal effect of the local severity is very small. We use the increment of 100 cases because the average cumulative case count as of June 2020 is 345, as shown in Table 2 . The result shows that the impact of local pandemic severity on career development stress is 0.001 standard deviations for an increase of 100 cumulative cases at the provincial level. The effect is miniscule and economically insignificant. To estimate both overall effect P T and local effect P L , we use data before and during the pandemic by applying the before-after estimator (Smith and Todd, 2005) . In our data, the Wave-1 sample was collected in January 2020 before the pandemic became public. During the period between the Wave-1 and Wave-2 surveys, the pandemic unfolded. However, between Wave-1 and Wave-2, some other effects overlapped with COVID, including 1) the YTEP training effect and 2) the seasonal pattern of stresses. Regarding the YTEP training effect, the first half of the program ends in January, and the entire program ends on June. Thus, the program may alter participants' job attitudes in the second half of training given its design. For the seasonal effect, people's feelings are affected by cold (January) or hot (June) weather. According to Nelson and Martin II (2010) , seasonal changes in the quality and quantity of stress are common, and the stress response is stronger during the winter. Cooke et al. (2007) compares the seasonal (winter and summer) effectiveness of aromatherapy massage on the stress and anxiety levels of emergency nurses and finds pre-massage anxiety was significantly higher in winter than summer. Conversely, January represents the beginning of a new year, which may bring more hope and inspiration and thus may contribute to higher job morale. The training effect and seasonal impact on job sentiments result in a time-specific intercept common across individuals which causes the before-after estimation to break down (Smith and Todd, 2005) . One way to avoid the time-specific intercept is to use samples from Class-2019 and the Wave-2 of Class-2020 surveys. Because both surveys were conducted at the end of the YTEP training program and in summertime, there will be no difference in seasonal and training effects between those two samples, assuming that seasonal patterns do not vary across years and the YTEP training effect is similar in different years. Therefore, we use Class-2019 as counterfactual in the same period and apply the before-after estimation to the following models by pooling Class-2019 and Class-2020 Wave-2 samples together: J o u r n a l P r e -p r o o f where is a dummy variable that equals 1 if surveyed in 2020, and measures the overall effect of the pandemic. Other variables have the same definitions as above in Model (3). The results are reported in Table 4 . The estimated pandemic effects are consistent using both categorical and scale measures. The pandemic increases both teaching stress and career development stress and reduces passion for the teaching occupation. All estimated overall effects are statistically significant. In contrast, the effects of the local severity are either statistically insignificant or economically insignificant. The estimated total effects are much larger than the marginal effects of local severity. Given that the pandemic had been spreading already, the effect of local confirmed cases is relatively smaller compared to the overall pandemic effect. Additionally, for most people, the major impact on attitudes comes via various preventative measures rather than experiences or awareness of the direct pandemic illness on individuals around them. This result also shows that cross-sectional data during the pandemic can only capture a very small part of the total impact of COVID-19. We also report the estimates for control variables in Table 4 . These variables may influence an individual's stress coping abilities as well as their social network. The estimated parameters generally have the expected sign and significance. In particular, teachers with more experience feel less teaching stress, probably because they are better trained for teaching. Female teachers are more likely to feel higher teaching stress. Compared to schoolhouse teachers, those who work in larger schools have higher teaching stress, possibly due to higher instructional expectations at their schools. Compared to the relatively less developed western region, rural teachers in other regions feel less teaching stress. However, those with more years of teaching experience have higher career development stress. Special-term teachers have less career development stress, likely because they are less concerned about their future career as a teacher. As for job passion, being married and having children are both positively associated with the teacher's job enthusiasm. Those who work in larger schools have lower enthusiasm towards the occupation compared to teachers at small schoolhouses, possibly due to the fact that a teaching job is highly respected in a remote village. Special-term teachers indicate higher job enthusiasm. One issue with pooling only the samples of Class-2019 and Wave-2 of Class-2020 in the J o u r n a l P r e -p r o o f above estimation is that the Wave-1 sample is not used. This sample provides additional information about job sentiments and their variations. Therefore, applying the DD estimation allows us to use all data to improve the estimation efficiency. We define the treatment group as survey participants who took the survey in both January and June of 2020 and who experienced the COVID-19 pandemic ("treatment") for five months. Ideally, we need a control group who took the January survey but did not go through the pandemic. However, such data cannot exist because everyone experienced the pandemic one way or another. Therefore, we consider one potential control group, the Class-2019 who participated in the survey conducted in June 2019. This group was not under any pandemic influences but participated in similar YTEP training and did the survey at the end of the training in the summer. Additionally, both the 2019 and 2020 surveys asked similar questions and were implemented in the same way. Given the homogenous nature of the samples, i.e., young rural teachers, it is reasonable to assume that the training effect and the seasonal influence are similar for the 2019 and 2020 cohorts. Therefore, we can difference out the seasonal and training effects and identify the net effect of the pandemic. Our DD framework is illustrated in Figure 2 . The survey respondents in January 2020 are divided into two groups, a treatment group and a control group. The treatment group went through the pandemic "treatment" and was surveyed again in June. The control group did not go through the treatment but were "as if" a counterfactual represented by Class-2019. As shown in Figure 2 , a total of 3,076 samples took the survey in January 2020 (Wave-1), and 2,155 of them took the survey again in June (Wave-2). There are two different ways to specify the control group before the treatment for participants in Wave-1. One way is to include only those who did not participate in Wave-2 (921 observations); the other way is to include all participants in Wave-1 (total 3,076 observations). The first option is closer to the standard DD estimation. However, the advantage of the second option is that the sample size is larger. Therefore, we adopt the second approach and use all observations of Wave-1 as control group for January. 18 Based on Duflo (2001), Bertrand et al. (2004) , Hansen (2007) , and Imbens and Wooldridge (2009) , we specify the DD empirical model as below: where represents measures of job stress or enthusiasm for individual i in group g surveyed 18 As a robustness check, we did the estimation using the first option for the control group which includes only those that did not participate in Wave-2 survey. The results are generally consistent. Due to the nonlinearity of the probit DD model, Ai and Norton (2003) and Puhani (2012) show that the true treatment effect is not the coefficient of the interaction term, although has the same sign. We estimate the treatment effect of the probit model in the DD estimation following Puhani (2012) . 20 However, the standard error using the Delta Method is computationally complicated, and thus we obtain the standard error of the treatment effect using bootstrapping. 21 Based on the probit estimation, the pandemic increases an individual's probability of having high teaching stress by 8.3 percentage points (treatment effect), which is statistically significant. The pandemic also has a statistically significant effect on job enthusiasm. In particular, it reduces a young rural teacher's job enthusiasm by approximately 0.486 standard deviations in the scale model, which is a decrease from the 50 th percentile to the 31 st percentile. This result is consistent with Alves et al. (2020) , which finds that the pandemic has reduced the perception of well-being in the profession. The coefficient of the treatment group dummy is statistically insignificant for all models. It supports our assumption that there is no systematic difference between the treatment and control group. The key assumption in the DD approach is that participants in all three surveys are representative samples from the same population. Given the design of the YTEP training program is for young rural teachers and the way the participants are selected, the assumption is likely to hold. One issue is that, as seen in Table 2 , the much higher proportion of non-western participants in 2019 could imply this cohort has better access to resources and social networks that reduce stress. To test this assumption, we estimate the above model using subsamples from different regions. The estimated COVID effects are comparable for the west and non-west samples, except that the pandemic's effect on teaching stress becomes statistically insignificant for the western sample. 22 Therefore, the potential regional difference does not alter the estimated results due to the pandemic. In this section, we further test the robustness of the findings reported above. In particular, given that for the control group, those before and after the treatment are not the same individuals, we adopt a matching process to select individuals who are similar. We match the sample of Class-2019 with the Wave-1 of Class-2020 using multiple matching (Stuart, 2010) . Because individuals in the samples are similar, we can exactly match multiple variables, including demographic characteristics and job characteristics that may influence job stress and enthusiasm. The matching process involves: 1) matching exactly demographic and job characteristics; 2) doing nearest neighbor matching with age and years of teaching experience. 22 Since the proportions of special-term and permanent teachers between Class-2019 and Class-2020 differ by more than 10 percentage points, we have also estimated the effects using separated samples. The results are very robust. Importantly, to focus on potential regional differences in job sentiments, we construct the control groups to have similar sample proportions between western and non-western regions as that of the treatment groups. Because the proportion of participants from the west is 72% in the treatment group but only 32% in Class-2019, we retain all the western samples in Class-2019 and match those with western samples in Wave-1 of Class-2020; we then match non-western participants. There are more non-western observations than needed in Class-2019 because of matching the regional proportion. We then randomly select a sample from the matched non-western observations. This process results in matched control group samples between Class-2019 and Class-2020 Wave-1, with the western vs. non-western ratio of 62:38, which is closer to that of 72:28 for the treatment group. 23 Moreover, we estimate a differenced DD model to show the control for observed heterogeneity between "matched" samples. The Model (5) is restructured as below: where and ′ represent measures of job sentiments after and before the treatment. This model shows that any imperfect matching in observed heterogeneity has been accounted for in the model due to before and after the pandemic. More specifically, Model (6) shows that differences in job sentiment before and after treatment are caused by: 1) time trend γ; 2) the treatment effect δ; 3) observed heterogeneity; and 4) unobserved heterogeneity. The treatment effect of the pandemic is consistently identified if the unobserved heterogeneity satisfies that ( − ′ ) = 0. The estimations based on Model (5) and Model (6) are asymptotically identical but may differ in finite samples (e.g., due to changes in degrees of freedom). In Table 6 , we report the results based on the matched sample using both the regular DD estimation (Model 5) and the differenced DD estimation (Model 6). The estimated effects in both models are very close, and moreover they are consistent with those in Table 5 . For example, in the differenced DD estimation, the pandemic increases an individual's teaching stress and career development stress levels by 0.116 standard deviations and 0.173 standard deviations, respectively, and the effects are significant at the 5% level. The pandemic also shows a significantly negative effect on job enthusiasm, approximately 0.46 standard deviations. Therefore, the results are robust to different samples used and to different estimations. 23 To balance the western vs. non-western sample ratio and sample size, we cap the difference in the regional sample ratios to be within 10 percentage points. Another potentially more difficult issue with our DD estimation is that, if the seasonal and/or training effects differ between 2019 and 2020, the estimated pandemic effect may still capture some of those influences. In this case, the expected value of the error terms in Model (6) is not 0, i.e., ( − ′ ) ≠ 0. However, our assumption that the seasonal pattern is constant across years is commonly used in techniques for seasonal adjustment. Additionally, the YTEP training effect is likely to be constant across years due to its design. We further investigate those assumptions to get additional information about the pandemic effect. One possible way to investigate the effect of unobserved differences in the seasonal pattern and training effects is to find another measure that could reflect such differences, so that we can remove them in the estimation. Based on the data availability, a potential candidate is the stress of social interaction in rural areas. The YTEP program is to help teachers who work and live in rural areas improve their social interactions, and thus the stress of social interactions may be influenced by the YTEP training. Stress of social interaction may also fluctuate seasonally. A particular question in the survey related to social interaction stress is, "I feel pressure to interact with others in rural areas." Following the same scaling method as for job sentiments, we define a standardized scale measure for social interaction stress (details are shown in Appendix Table A2 ). The DD model for social interaction stress can be written as: where measures social interaction stress. If the difference in the social interaction stress between Class-2019 and Wave-2 of Class-2020 captures the unobserved differences between those two groups, it could help difference out the bias in estimating the effects of the pandemic. This approach works similarly to the difference-in-difference-in-differences (hereafter "DDD") estimation (Wooldridge, 2010) . We specify our model as: One particular issue in the above model is whether the stress of social interaction is affected by the pandemic. As shown in the second term of Model (8), if the pandemic has no effect on social interaction stress, then the estimate of − represents the total pandemic effect on job sentiments; otherwise, it is a relative pandemic's effect on job sentiments netting out the pandemic's impact on social interaction stress. The pandemic causes various restrictions on social interactions. It might increase people's anxiety about communicating with others, or it might reduce stress by reducing the need to interact with each other. In other words, the hypothetical direction of the potential impact is mixed. We estimate the DD model for the social interaction stress and the results are reported in Column 1 of Table 7 . It appears that the COVID-19 pandemic has no significant effect on social interaction stress for rural young teachers. The results of the "DDD" estimation for job sentiments are shown in Columns 2-4 of Table 7 . Compared to the DD results in Table 6 , the pandemic's effects on teaching stress, career development stress and job enthusiasm are quite robust, with the former two increasing by 0.110 and 0.167 standard deviations, and the latter decreasing by approximately 0.466 standard deviations; both are highly significant and are similar in magnitude to the previously reported DD estimates. Therefore, after differencing out the potential unobserved heterogeneity between the two samples, the estimated pandemic effects on job stress and job passion remain quite robust. Given the above results that the pandemic has significant impact on rural teachers' job sentiments, in this section, we investigate some working channels underlying the effect. The YTEP surveys provide some information on rural teachers' work activities, which can help identify working channels of the pandemic on their job stress and passion. The COVID-19 has changed work patterns and workload because teachers need to learn new teaching formats, conduct instruction via various online platforms, manage students' learning J o u r n a l P r e -p r o o f online, etc. For example, Yang (2020) finds that around 63% of teachers find using online education platforms difficult. The reasons include factors such as instability of internet connections and online platforms, unfamiliarity with relevant technology, difficulty in controlling the progress of the course, and limited interaction with students. Moreover, after schools reopened, to be prepared for a pandemic resurgence, teachers were required to get ready to switch between face-to-face and online teaching modes at any time. In many places, schools offer hybrid classes with both in-person and online formats. In addition, teachers are also responsible for additional administrative work, such as epidemic prevention for the school and students and regular data reporting, etc. 24 Therefore, even if hours worked may not increase, the changes due to the pandemic create additional work responsibilities for teachers. As shown in Table 8 , rural young teachers taught a relatively larger number of classes across time, for example, 14.92 in June 2019 vs. 18.83 in June 2020. They also gave students many more weekly homework assignments in Class-2020 than in Class-2019. Besides the above channel of workloads, we also look at social network at work. Based on the theoretical framework, the social network can have both buffering effects and spillover effects on an individual's stress. We use the information reported by the teachers about how many close colleagues they have in the school. Table 8 shows that the proportion of teachers with some close colleagues is smaller for Class-2020. Job-related training programs may enhance an individual's stress coping ability. Data in Table 8 shows that the proportion of samples participating in other job training programs also declines across time. Based on the discussion about potential working channels on job sentiments, we first estimate the same DD model on these channel variables to assess how they are affected by the pandemic. The results show that the COVID-19 pandemic has a significant effect on increasing the number of weekly classes taught by rural young teachers, as well as on increasing the homework assigned, and both effects are statistically significant. 25 One potential explanation for the increasing teaching load is online teaching. Teachers can record instruction videos and share the same materials to students in more classes. Because of the difficulty in controlling online teaching quality, teachers gave homework assignments more frequently. Rural teacher's social network could be one working channel of the pandemic on teachers' 24 Source: https://www.sohu.com/a/378302780_260616. 25 The results based on the matched samples are similar. J o u r n a l P r e -p r o o f job sentiments. In particular, the pandemic displays a statistically significant effect in reducing close colleagues in school. It also significantly reduced the rural young teachers' participation in other training programs outside the YTEP. The supporting results show that all these potential channel variables are affected by the pandemic in the DD estimation. In order to investigate these working channels, we add them in the DD model for the pandemic. The results in Table 9 show that teaching load and homework generally result in higher job stresses and lower job passion. The estimated effects of average number of classes taught per week are much stronger and are statistically significant in affecting all three job sentiments, based on both categorical and level measures. The amount of homework assignments has a significant effect on teaching stress but mixed effects on career development stress and passion for the teaching job, as perhaps assigning and grading homework represents a smaller part of work compared to teaching classes. In addition, having close colleagues shows a statistically significant effect of reducing both teaching stress and career development stress, for both the category-and scale-based models. This result indicates that the buffering effects for stress from the colleague network exceed that of the spillover effects, and thus help reduce the stress. On the other hand, other training programs do not show a statistically significant effect on job sentiments. Comparing Table 9 to Table 5 Besides work activities, another possible working channel of the pandemic is that a higher teaching stress during the pandemic may transfer to higher career development stress. Therefore, we add the teaching stress to the model. The results reported in Table 9 indicate that teaching stress is strongly related to career development stress. More specifically, it reduces the direct effect of the pandemic for the career development stress from 0.157 standard deviations to 0.091 J o u r n a l P r e -p r o o f standard deviations. If the teaching stress measure is influenced by other factors such as difference in seasonality and YTEP training between the two years, this model will also help offset such influences on career development stress and mitigates the potential bias. As discussed in the conceptual framework, it is also possible that teachers' teaching stress and career development stress affect their passion toward the teaching occupation. Stressful feelings at work might make job activities less enjoyable and thus reduce an individual's work satisfaction and job passion. We further estimate the model of job passion by including job stress measures. We find that both teaching stress and career development stress negatively affect job enthusiasm and are statistically significant (Column 8, Table 9 ). Nevertheless, the joint effect of both stress measures is -0.041 (by adding -0.030 and -0.011), which is much larger than the decline in the direct effect of the pandemic (i.e., from -0.486 in Column 6 of Table 5 to -0.473 in Column 8 of Table 9 ). All the results show that, even after taking out the impact due to behavioral changes or other potential effects, we still find a statistically significant direct effect of COVID-19 on job sentiments. It appears that a part of the change in job sentiments is not directly related to behavioral changes on the job or changes at the local level but is likely caused by subtle anxieties and fears about the uncertainties and risks all over the world due to the COVID-19 pandemic. The COVID-19 pandemic has greatly changed the way people work and live. In this study, we use unique survey data to investigate the pandemic's effects on rural teachers' job-specific stresses and their enthusiasm for the teaching profession. We propose a theoretical framework for the dynamics of an individual's stress formation, construct empirical models, and apply various methods in the estimation. Our results show that the COVID-19 pandemic raises a rural young teacher's probability of having high teaching stress and career development stresses by 8-9 and 5-7 percentage points, respectively. With the standardized scales, the pandemic increases the teaching stress by 0.11-0.13 standard deviations and the career development stress by 0.15 to 0.18 standard deviations. In addition, the pandemic reduces a teacher's enthusiasm towards the teaching occupation. The local severity of COVID-19 affects job sentiments, but its relative magnitudes are much smaller than the overall pandemic effect and the impact is economically insignificant. We investigate the channels through which COVID-19 affects teachers' job attitudes. Changes in work-related activities can help explain some of the pandemic effect. However, even after controlling for those channels, the COVID-19 pandemic still displays a strong direct influence on teachers' job sentiments. This result sheds light on the ways that the pandemic fosters overall anxiety and a pessimistic social atmosphere, thereby exerting a direct impact on job sentiments. Because teachers' morale affects education outcomes for the next generations, their mental health and psychological changes during the COVID-19 pandemic need to be addressed. Governments at various levels can provide teachers with better tools and training for online teaching along with services to reduce the psychological effects during the pandemic. It is also important to reduce the burden of frequent switches between instructional formats. Furthermore, it is helpful to provide timely and scientifically based information to calm the pandemic-related anxiety for the population at large. Appendix Notes: To ensure that the scaling was assigned as objective as possible, each member of the research team assigned scale values in three rounds at different time independently. We then average scale values across members and rounds for use in the estimation. The Fleiss' Kappa statistic on interrater reliability for the choices is 0.51, which indicates moderate agreement among raters. Appendix Figure A1 The Process to Matching Samples J o u r n a l P r e -p r o o f Journal Pre-proof Step 1:  Define Dataset A: Class-2019; Dataset B: Class-2020 Wave-1. We select observations from B to match A.  Make matching according to the observed individual characteristics (including gender, education, ethnicity, marriage, having children) and job characteristics (including teacher type, school type), but exclude years of teaching experience  In Dataset A, a total of 1,490 observations are matched. Step 2:  Among the observations matched in step 1, most have multiple matches.  We further match these observations based on ages. To narrow down the number of matches, we kept pairs with the smallest age gap.  After matching, the age differences in the matching samples have a mean of 0.043 years and standard error of 0.649.  We then match years of teaching experience in a similar fashion, with a mean difference in teaching experience between the matched observations of 0.040 years and standard error of 0.691. Step 3:  After steps 1 and 2, some observations still match more than once.  In this case, we randomly select a match pair based on ID.  If one observation in Dataset B matches to more than one observation in A, then one of the matched in A will be removed and will be rematched in later steps. In the matched dataset A', we randomly keep 290 observations from the non-west region and all observations from the west region to match the regional distribution in Dataset B. Step 4:  For those not matched in the first round described above, we follow the same procedure to match them.  The process is completed after 11 rounds. Table 4 . 5. Robust standard errors in parentheses and * p < 0.1, ** p < 0.05, *** p < 0.01 Notes: 1. The treatment effects in probit models are reported and they are adjusted for nonlinearity in the DD estimation, and their standard errors are calculated using bootstrapping. The marginal effects of other variables in the probit model are reported and are calculated using the average marginal effects. 2. The results in Columns 4-6 are based on the standardized scale with Wave-1 (January 2020) sample as the base. J o u r n a l P r e -p r o o f 3. Robust standard errors in parentheses and * p < 0.1, ** p < 0.05, *** p < 0.01. Notes: 1. The results are based on the standardized scale with Wave-1 (January 2020) sample as the base. 2. Other control variables included are the same as those listed in Table 4 . 3. Robust standard errors in parentheses and * p < 0.1, ** p < 0.05, *** p < 0.01. Notes: 1. The results are based on the standardized scale with Wave-1 (January 2020) sample as the base. 2. Other control variables included are the same as those listed in Table 4 . 3. Robust standard errors in parentheses and * p < 0.1, ** p < 0.05, *** p < 0.01. J o u r n a l P r e -p r o o f Notes: 1. The variable "Close colleagues" is defined as 1 if the respondents report having some or many close colleagues in the Class-2019 survey; and as 1 if the respondents reported having three or more close colleagues in the Class-2020 surveys. 2. Standard deviations in parentheses. Inequality in the impact of the coronavirus shock: Evidence from real time surveys Interaction terms in logit and probit models Teachers' well-being in times of Covid-19 pandemic: factors that explain professional well-being Identification and Inference in Nonlinear Difference-in-differences Models The impact of COVID-19 on student experiences and expectations: Evidence from a survey The skill content of recent technological change: An empirical exploration When face-to-face interactions become an occupational hazard: jobs in the time of COVID-19 How much should we trust differences-in-differences estimates? Stress and well-being at work: A century of empirical trends reflecting theoretical and societal influences COVID-19, lockdowns and well-being: Evidence from Google Trends The role of passion for teaching in intrapersonal and interpersonal outcomes School climate and social-emotional learning: Predicting teacher stress, job satisfaction, and teaching efficacy The effect of aromatherapy massage with music on the stress and anxiety levels of emergency nurses: comparison between summer and winter Productivity effects of innovation, stress and social relations Schooling and labor market consequences of school construction in Indonesia: Evidence from an unusual policy experiment Human mobility restrictions and the spread of the Novel Coronavirus (2019-nCoV) in China Measuring nominal scale agreement among many raters The anxiety-buffering function of close relationships: evidence that relationship commitment acts as a terror management mechanism The Oxford handbook of work engagement, motivation, and self-determination theory The longitudinal effects of work-family conflict and positive spillover on depressive symptoms among dual-earner couples Generalized least squares inference in panel and multilevel models with serial correlation and fixed effects Understanding the level and causes of teacher turnover: A comparison with other professions Teacher well-being: Looking after yourself and your career in the classroom Multidisciplinary research priorities for the COVID-19 pandemic: a call for action for mental health science On passion and heavy work investment: personal and organizational outcomes Handbook on facilitating flexible learning during educational disruption: The Chinese experience in maintaining undisrupted learning in COVID-19 Outbreak Recent developments in the econometrics of program evaluation Impact on mental health and perceptions of psychological care among medical and nursing staff in Wuhan during the 2019 novel coronavirus disease outbreak: a cross-sectional study Effects on teachers' self-efficacy and job satisfaction: Teacher gender, years of experience, and job stress Seasonal changes in stress responses Learning during the COVID-19 pandemic: it is not who you teach, but how you teach The treatment effect, the cross difference, and the interaction term in nonlinear "difference-in-differences" models Impacts of social and economic factors on the transmission of coronavirus disease 2019 (COVID-19) in China A nationwide survey of psychological distress among Chinese people in the COVID-19 epidemic: implications and policy recommendations Does matching overcome LaLonde's critique of nonexperimental estimators? Matching methods for causal inference: a review and a look forward Is returning to work during the COVID-19 pandemic stressful? A study on immediate mental health status and psychoneuroimmunity prevention measures of Chinese workforce Longitudinal trajectories of passion and their individual and social determinants: A latent growth modeling approach Passion at work: Toward a new conceptualization The Oxford handbook of work engagement, motivation, and self-determination theory Immediate psychological responses and associated factors during the initial stage of the 2019 coronavirus disease (COVID-19) epidemic among the general population in China A longitudinal study on the mental health of general population during the COVID-19 epidemic in China Stress and strain crossover Has SARS infected the property market? Evidence from Hong Kong Econometric analysis of cross section and panel data Teachers' perceptions of large-scale online teaching as an epidemic prevention and control strategy in China Impact of the COVID-19 pandemic on mental health and quality of life among local residents in Liaoning Province, China: A cross-sectional study Notes:1. The treatment effects in probit models are reported and they are adjusted for nonlinearity in the DD estimation, and their standard errors are calculated using bootstrapping. The marginal effects of other variables in the probit model are reported and are calculated using the average marginal effects. 2. The results in Columns 4-8 are based on the standardized scale with Wave-1 (January 2020) sample as the base. 3. Other control variables included are the same as those listed in Table 4 . 4. Robust standard errors in parentheses and * p < 0.1, ** p < 0.05, *** p < 0.01. 5. Notes: "Cumulative cases" represents the cumulative number of confirmed cases in China since the outbreak. "Existing cases" is the current number of confirmed cases. "New cases" is calculated as the change in the number of cumulative cases compared to the previous day. Highlights  Using unique survey data from China and applying various estimation methods.  Among rural young teachers, the pandemic has caused higher teaching stress and career development stress and has reduced passion towards the teaching occupation.  The working channels of the pandemic include job-related activities and social network.  After controlling for possible working channels, the COVID-19 pandemic still shows a strong direct impact on job sentiments.