Effects of school climate and teacher self-efficacy on job satisfaction of mostly STEM teachers: a structural multigroup invariance approach

Background: Identification and retention of effective teachers in STEM education play cardinal roles in teacher recruitment exercises worldwide. Studies on factors that characterize effective teachers have therefore gained popularity in recent times. Teacher self-efficacy, job satisfaction and school climate are among other factors that have attracted global attention. Thus, proper understanding of the relations between these factors is equally important. The purpose of this study is to validate and cross-validate a model of direct/indirect effects of school climate and teacher self-efficacy on job satisfaction. Results: The data used for the current study are extracted from a publicly available data of Teaching and Learning International Survey (TALIS) 2018 survey. Structural equation modeling approach was used in the analyses coupled with robust maximum likelihood to ensure accurate estimations in the models. The results of the validated models show a strong direct impact of school climate on job satisfaction, a direct impact of teacher self-efficacy on job satisfaction and a mediating effect of teacher self-efficacy between school climate and job satisfaction. This model exhibits structural invariance in factor loadings, intercepts and regression weights across two independent samples from a population of 3951 lower secondary school teachers in Norway. Conclusion: The findings of this study do provide empirical evidence for the relations between teacher selfefficacy, job satisfaction and school climate among Norwegian lower secondary school teachers. The crossvalidation of these relations was also established using an independent sample to enhance generalization of the findings. Two methodological observations concerning recoding of some items as well as an addition of item cross-loading in the measurement model of the job satisfaction scales are raised and addressed. It is therefore recommended that researchers who will be using TALIS 2018 data should take note of these observations.


Introduction
Identification and retention of efficacious teachers have long been playing cardinal roles in teacher recruitment exercises worldwide. This occurrence can be linked to the pronounced role of teacher self-efficacy in improving quality of learning outcomes (e.g. Caprara, Barbaranelli, Steca, & Malone, 2006). Globally, increased attention has been witnessed in recent times especially in the aspect of factors that enhance or inhibit self-efficacy of primary and secondary school teachers. For example, multiple results from the Organization for Economic Co-operation and Development (OECD) and Teaching and Learning International Survey (TALIS) 2013 and 2018 studies have shown that teacher perceived disciplinary (which is an aspect of school climate factor) has a negative correlation with teacher self-efficacy (OECD, 2014(OECD, , 2019a. Further, teachers who had high scores on their job satisfaction scales are identified with high selfefficacy as well. Explorations of the relationships between these constructs-teacher self-efficacy, school climate and job satisfaction-have many implications for teachers' identification, retention and recruitment and by extension quality of learning outcomes. Several attempts have been made to establish these relationships in the literature. However, the results are inconsistence or rather inconclusive. For instance, Caprara et al. (2006) in a study involving 2184 junior high school teachers reported a significant impact of teacher selfefficacy on their job satisfaction. Similar finding of a strong predictive power of teacher job satisfaction by teacher self-efficacy has also been reported in (Skaalvik & Skaalvik, 2014;You, Kim, & Lim, 2016). On the other hand, Blackburn and Robinson (2008) in their study involving 80 early career teachers found a significant positive correlation between 'teacher self-efficacy in student engagement' and job satisfaction, and a significant negative correlation between 'teacher self-efficacy in classroom management' and job satisfaction. Some aspects of the school climate such as teacher-student relations have also been reported to correlate positively with teacher self-efficacy (e.g. Hosford & O'Sullivan, 2016). More recently, in contrary to the earlier findings Shaukat, Vishnumolakala, and Al Bustami (2019) reported no significant relationship between teacher self-efficacy and job satisfaction.
Several other researchers have adopted structural equation modeling approach to unravel the causal relationships between these constructs. An empirical longitudinal study on 642 teachers in lower secondary schools revealed that school climate exerts a positive impact on job satisfaction. However, this impact is mediated by teacher self-efficacy (Malinen & Savolainen, 2016). This finding is consistent with a significant mediating effect of teacher self-efficacy between school climate and job satisfaction reported in (Aldridge & Fraser, 2016). The main issue with most of these structural equation modeling approaches to testing hypothesized relationships between these constructs is lack of cross-validation. It is a common practice among researchers to test a structural model and improve on it to achieve model fits by conducting post hoc analyses without any recourse to further validations in independent samples (e.g. Aldridge & Fraser, 2016). Meanwhile, the post hoc analysis procedure is purely exploratory since it involves adding or freeing constraints as suggested by the data and therefore limits generalization. A way of circumventing this problem is to conduct a cross-validation of a well fitted model through multigroup invariance. Another solution to this problem could also involve conducting longitudinal studies (Byrne, 2012). This lack of cross-validation of the tested structural model on school climate, teacher self-efficacy and teacher job satisfaction coupled with inconclusive results on the relationships between these constructs have created knowledge gaps begging for more investigations. The current study is aimed at providing empirical evidence to fill these gaps within a national context and by extension to a global community. In specific terms, the current study seeks to provide empirical evidence for: (a) a structural model that explains the direct effect of school climate on job satisfaction, direct effect of teacher self-efficacy on job satisfaction and indirect effect of school climate on job satisfaction as mediated by the teacher self-efficacy; (b) cross-validation of this structural model using independent sample from the same population; (c) multigroup invariance testing of the estimated parameters and regression weights of the structural model. In the forthcoming section, a theoretical background coupled with conceptual propositions of this study are presented.

Conceptual framework
A theoretical model that could explain teacher job satisfaction and its associated teachers' personal traits, e.g. teacher self-efficacy and environmental factors, e.g. school climate, is the social ecological theory (Bronfenbrenner, 1986;Darling, 2007). This grand theory sees an actor (e.g. a teacher) as a person whose intrinsic feelings (e.g. job satisfaction, self-efficacy) are constantly being shaped by interaction with the environment (e.g. school climate). From this theoretical perspective, one can justify the rationality of establishing empirical relations between these constructs. However, in order to maintain consistency with the source of data used in the current study, subsequent discussions are based on the conceptual framework and technical report as released by the TALIS 2018 team (Ainley & Carstens, 2018;OECD, 2019b). TALIS is an international survey that is sponsored by the OECD to provide empirical evidence for the teaching and learning practices among teachers and school principals. It is cyclical survey that is conducted every 5 years and covers several aspects of teachers'/ principals' factors including instructional practices, job satisfaction, professional practices, self-efficacy, school climate, school leadership, teacher feedback and development. Thus, the current study relies on TALIS for its data source being the only large scale publicly available data on the research constructs under investigation.

Teacher job satisfaction
Teachers' intrinsic evaluations of their jobs have been found to be strong indicators of retention on teaching jobs and predictive of teachers' tendency to quit their jobs (e.g. McConnell, 2017;Skaalvik & Skaalvik, 2011).
Teacher job satisfaction according to Ainley and Carstens (2018) connotes "the sense of fulfilment and gratification that teachers experience through their work as a teacher"(p.43). It is a mixture of positive and negative evaluations teachers make about the teaching profession (Skaalvik & Skaalvik, 2015). It is a multidimensional construct such that a teacher may make positive evaluations of one or more aspects of it and yet maintains negative evaluations of the others. Van der Ploeg and Scholte (2003) in (Veldman, van Tartwijk, Brekelmans, & Wubbels, 2013) identified five dimensions of job satisfaction: "management support, autonomy, relations with colleagues, nature of work, and working conditions" (p.56). These dimensions are reduced to three in the construction of scales used in the current study. These three dimensions are (a) "teacher job satisfaction with work environment" (JSENV), (b) their "job satisfaction with profession" (JSPRO), and (c) their "job satisfaction with target class autonomy" (JSTCA), see OECD (2019b, p. 293).
The relations between teacher job satisfaction, and by extension its dimensions, with other constructs have been extensively studied in the literature. For example, Skaalvik and Skaalvik (2011) found that teachers with high job satisfaction exhibit low motivation to leave the profession (β = − .53, p < .05), those with high sense of belonging exhibit high job satisfaction (β = .28, p < .05) and teachers with high emotional exhaustion developed low job satisfaction (β = − .53, p < .05). Similar corroborative findings have also been reported elsewhere (e.g. Skaalvik & Skaalvik, 2015). A handful of empirical studies have also confirmed positive direct and interactive effects between teacher job satisfaction and school climate-related factors, e.g. teachers' collaboration in improving students' mathematics and reading growth as well learning outcomes (e.g. Banerjee, Stearns, Moller, & Mickelson, 2016). In a more recent study, it is found that job satisfaction positively correlates with teachers' perceived school leadership styles and that this correlation is mediated by teacher professional collaboration (see García Torres, 2019). These findings suggest that the implications of high and/or low teacher job satisfaction are not only limited to teachers' well-being-related factors but also extend to students' learning outcomes, school organization and the community at large.

Teacher self-efficacy
The origin of self-efficacy, i.e. internal convictions to perform tasks with an estimation of desired outcomes, has been linked to the works of Albert Bandura. It is postulated to have four sources of influence: mastery experience, vicarious experience, verbal/social persuasions and physiological or affective states (Bandura, 1997;Zakariya, Goodchild, Bjørkestøl, & Nilsen, 2019). It is an important affective factor that boosts job performance and a strong militating factor against attrition and absenteeism. Teacher self-efficacy according to Ainley and Carstens (2018) connotes 'the beliefs that teachers have of their ability to enact certain teaching behaviour that influences students' educational outcomes, such as achievement, interest, and motivation' (p. 51). The sources of self-efficacy as postulated by Bandura are recently investigated in a cross-national study among the Japanese and Finnish teachers. It was found that mastery experience has the strongest influence on the teachers' self-efficacy across the two countries (Yada et al., 2019). Therefore, from a theoretical perspective, teacher selfefficacy is a multifaceted affective construct. As such, three dimensions are conceptualized and operationalized to reflect teachers' self-efficacy by the TALIS 2018 team. These dimensions are: (a) "self-efficacy in classroom management" (TSECLS), (b) "self-efficacy in instruction" (TSEINS), and (c) "self-efficacy in student engagement" (TSEENG), see OECD (2019b, p. 275).
Several factors have been linked empirically to teacher self-efficacy in recent times. These include (a) factors that increase teacher self-efficacy, e.g. innovative teaching strategies and supportive school environment (Choi, Lee, & Kim, 2019;Fackler & Malmberg, 2016); (b) factors that decrease teacher self-efficacy, e.g. burnout and stress (Zhu et al., 2018); (c) factors that are positively affected by teacher self-efficacy, e.g. job satisfaction (Türkoğlu, Cansoy, & Parlar, 2017). Teachers' job satisfaction appears to be the most widely studied factor in relation to teacher self-efficacy. A handful of studies have established a strong positive relation between teacher self-efficacy and job satisfaction with former enhancing the latter (e.g. Edinger & Edinger, 2018;Skaalvik & Skaalvik, 2014). A team of researchers in an empirical study involving 500 high school teachers using cluster analysis found that low burnout increases self-efficacy which in turn improves job satisfaction (Molero Jurado, Perez-Fuentes, Atria, Oropesa Ruiz, & Gazquez Linares, 2019). These findings partly agree with an earlier study that confirmed a significant positive direct impact of teacher self-concept on self-efficacy and a significant negative direct impact of burnout on self-efficacy (Zhu et al., 2018). Thus, based on the aforementioned body of literature, it is hypothesized that: (H01) Teachers' self-efficacy in classroom management, their self-efficacy in instruction and their selfefficacy in student engagement have direct positive effects on teachers' job satisfaction.

School climate
School climate which is sometimes referred to as school culture has consistently been identified to contribute significantly in sustaining school improvement, teachers' well-being and students' learning outcomes. A conducive and inclusive school climate that is free of hostility and bullying enhances teachers' job satisfaction, boosts selfefficacy, improves teacher-student relationship, reduces attrition and burnout (Lee & Louis, 2019;Van Beurden, Van Veldhoven, Nijendijk, & Van De Voorde, 2017). School climate is equally beneficial not only to the teachers but also to students' behaviour, approaches to learning, total well-being and improved achievement (Zakariya et al., 2019). According to Ainley and Carstens (2018), school climate is a multidimensional construct that encapsulates "safety, relationships, engagement with teaching and learning, institutional environment, and school improvement activities" (pp. 37-38). Thus, TALIS 2018 team conceptualized and operationalized three dimensions of this construct which are (a) "teachers' perceived disciplinary climate" (TPDC), (b) "teacher-student relations" (TSR), and (c) "participation among stakeholders" (TSTAKE), see OECD (2019b, p. 329).
There have been mixed results on the impacts of school climate on other constructs such as teacher selfefficacy and job satisfaction. These mixed results stemmed from different facets of the school climate. For example, it was found in a qualitative longitudinal study involving four highly job satisfied teachers that teacherstudent relationship is strongly positively related to teacher job satisfaction (Veldman et al., 2013). More so, teachers with high sense of belonging to the school system and have more opportunity to participate among stakeholders are reported to have high self-efficacy which translates to high job satisfaction (e.g. Skaalvik & Skaalvik, 2011;You et al., 2016). On the other hand, teachers' perceived disciplinary climate has been reported to negatively impact self-efficacy which in turn translates to low job satisfaction (Aldridge & Fraser, 2016;Malinen & Savolainen, 2016). Building on this body of literature, the following hypotheses are formulated: (H02) There is a direct positive effect of teacherstudent relations on job satisfaction.
(H03) There is a direct negative effect of teacher perceived disciplinary climate on job satisfaction.
(H04) The direct effects of each dimension of the school climate on job satisfaction are mediated by teacher self-efficacy.
(H05) The pattern of factor loadings, other estimated parameters, and regression coefficients are likely to be invariance across an independent sample. Figure 1 explicitly demonstrates a graphical structural model of the hypothesized relations between the constructs of the current study. To the best of author's knowledge, no large-scale empirical study has ever been reported within Norway on the structural invariance of a model relating school climate, teacher self-efficacy and job satisfaction. The current study seeks to complement decades of research in understanding, identifying, recruiting, and retaining effective teachers.
The oval shapes in Fig. 1 depict the dimensions of each of the constructs: school climate, job satisfaction and teacher self-efficacy. Both JSTCA and TSTAKE are Fig. 1 A conceptual model of the relations between school climate (TPDC, teachers' perceived disciplinary climate and TSR, teacher-student relations ), job satisfaction (TJSENV, job satisfaction with work environment and TJSPRO, job satisfaction with profession), and teacher self-efficacy (TSECLS, self-efficacy in classroom management; TSEINS, self-efficacy in instruction; and TSEENG, self-efficacy in student engagement) Zakariya International Journal of STEM Education (2020) 7:10 Page 4 of 12 not included in Fig. 1 because their corresponding data are not released by the TALIS 2018 team for lower secondary school teachers. The single-headed arrows indicate the direction of the theorized causal effects from one dimension of a construct to another. For instance, the single-headed arrow between TSECLS and TJSENV depicts a direct theorized casual effect of the 'teacher self-efficacy in classroom management' on the 'teacher job satisfaction with work environment' as hypothesized in H01. In the forthcoming section, data preparation, measures, procedure of analysis and other methodological related issues are presented.

Data preparation
The data used for the current study came from the TALIS 2018 international survey. The survey took place towards the end of 2017/2018 academic year for teachers and principals in many countries across primary (ISCED 1), lower secondary (ISCED 2) and upper secondary schools (ISCED 3). Two hundred schools per country and 20 teachers in each school are sampled using probabilistic technique. In line with the aim of the current study, the ISCED 2 teacher data for Norway are extracted and teacher responses related to the variables of interest are prepared. ISCED 2 teacher data are used because no data are available for both ISCED 1 and ISCED 3 in Norway. These data are screened for outliers and 203 respondents out of a total of 4154 are found with missing values on all variables. These respondents are removed, and the effective data become 3951 including 2541 males (64.3%) and 1410 females (35.7%). The main source of the data used for the current study is publicly available on the OECD website (http://www.oecd.org/ education/).

Measures
Teacher job satisfaction is measured by TJSENV (four items) and TJSPRO (four items) subscales each measured on a four-point Likert scale ranges from 1 (strongly disagree) to 4 (strongly agree). Teacher selfefficacy is measured by TSECLS, TSEINS and TSEENG subscales each measured by four items on a four-point Likert scale ranges from 1 (not at all) to 4 (a lot). School climate is measured by TPDC (four items) and TSR (four items) subscales each measured on a four-point Likert scale ranges from 1 (strongly disagree) to 4 (strongly agree). All the subscales of the ISCED 2 teachers have high internal consistency with Omega reliability coefficients range from .706 to .901 and .642 on TSEENG subscale. They also exhibit measurement invariance up to metric level for Norwegian data (OECD, 2019b). Table 1 presents items measuring each subscale (including means and standard deviations for our sample, N = 3951), their wordings and response options.

Recoding and missing data
It appears that items TG53C, TG53D, TG53F and TG41B are not reverse recorded to reflect their negatively worded nature before TALIS team made the 2018 data public. This is observed when these items show negative covariances with other subscale items in their respective constructs during data screening process. Thus, each of these items is recoded to reflect these changes. The data are screened for missing values using Little's missing completely at random (MCAR) tests. These tests show no significant statistics for each components of the teacher self-efficacy (χ 2 = 302.98, df = 317, p = .750 at 0.1-0.5% missing), borderline significant statistics for job satisfaction (χ 2 = 195.79, df = 165, p = .051 at 1.6-2.0% missing), and significant statistics for school climate (χ 2 = 134.25, df = 97, p = .007 at 1.0-14.2% missing). Thus, for convenient use of these data in subsequence analysis, full information maximum likelihood (FIML) estimation with expectation maximum (EM) algorithm is utilized as opposed to multiple imputation to replace the missing data (Cham, Reshetnyak, Rosenfeld, & Breitbart, 2017).

Data analysis procedure
The extracted data are randomly split into two date sets using split at random command in SPSS version 23. The first data set (data 1) is made up of 1975 respondents (1300 males and 675 females) while the second data set (data 2) is made up of 1976 respondents (1241 males and 735 females). The data analyses proceed in two stages. Stage 1 involves using data 1 to calibrate measurement models for each of the constructs and to validate the structural model as well as to conduct post hoc analysis. Stage 2 involves using data 2 to cross-validate the established structural model in stage 1. Structural multigroup invariance techniques are adopted in stage 2 which involves comparing configural model models with more restrictive models. A combination of methods are used to access model fits in both stage 1 and stage 2 which includes: (a) significant parameter estimates; (b) goodness of fits (GOF) indices CFI, TLI close to or ≥ .95, SRMR ≤ .08 and RMSEA ≤ .06; (c) chi-squared difference tests and changes in GOF indices, e.g. ΔCFI > − .010, ΔRMSEA < .015, and ΔSRMR < .010. These criteria are chosen because of their satisfactory performance reported in the literature (e.g. Chen, 2007;Hu & Bentler, 1999). Considering the large sample size involved in the current study, chi-square statistic is not used to assess model fits because it has been deemed sensitive to sample size (Brown, 2015). Mplus version 8.3 software is TG49D If a student needs extra assistance, the school provides it. 1 (strongly disagree)-4 (strongly agree)

.64
Note: Item with "*" indicates reverse coded. Construct labels and item wordings are adapted from (OECD, 2019b) Zakariya International Journal of STEM Education (2020) 7:10 Page 6 of 12 used for these analyses and the ensuing results are presented in next section.

Stage 1: Measurement models
Three separate measurement models are analyzed using maximum likelihood with robust standard errors (MLR), one model for each of the variables: teacher job satisfaction, self-efficacy and school climate. MLR estimator was used as opposed to maximum likelihood and others because of its robustness to violations of multiple assumptions of normality distribution that are inherent in ordinal scales with less than 5-point categories (Suh, 2015). A good model fit (χ 2 (17) = 152.18, p < .001, RMSEA = .063 (90% C.I. = 0.054-0.073, CFit = .008), CFI = .97, TLI = .94, and SRMR = .033) was found for job satisfaction after accounting for an error covariance between TG53D and TG53F and a cross-loading of item TG53J on 'job satisfaction with profession' construct. The error cross-loading was suggested by the TALIS 20018 team in (OECD, 2019b) while the item crossloading was suggested by modification indices in the output of the analysis. The inclusion of this item crossloading gave a significant MLR chi-square difference test with Satorra-Bentler correction Δχ 2 ½1 ¼ 171:50; p < :001. From a conceptual point of view, item TG53J: All in all, I am satisfied with my job, seems to fit very well to measuring 'teacher job satisfaction with profession' rather than 'job satisfaction with work environment' as postulated by the TALIS team. Thus, from conceptual/ theoretical and statistical perspectives, inclusion of item cross-loading on 'job satisfaction with profession' construct in the measurement model appears to be duly justified.
The statistics of the analyzed measurement model for teacher self-efficacy variable also suggest a good model fit: χ 2 (48) = 330.76, p < .001, RMSEA = .055 (90% C.I. = 0.049-0.060, CFit = .082), CFI = .95, TLI = .94, and SRMR = .037. This model fit was achieved after accounting two error covariances between items TG34A with TG34B and items TG34D with TG34I as recommended by the TALIS 2018 team (OECD, 2019b). Further, a good model fit (χ 2 (19) = 65.69, p < .001, RMSEA = .035 (90% C.I. = .026-.045, CFit = .995), CFI = .99, TLI = .99, and SRMR = .026) was found for the school climate variable. However, there was a negative significant standardized correlation (r = − .18, p < .01) between 'teachers' perceived disciplinary climate' and 'teacher-student relations' constructs of the school climate. This has an implication for the impact of each of these constructs on other variables as it could be interpreted to mean those teachers with smooth student relations also have good perception of the disciplinary measures in their classrooms. All the estimated parameters, e.g. factor loadings, correlation coefficients, and factor variances, are significant. The validated measurement models were used in the structural model whose results are presented in the next section.

Structural model
Three structural models are evaluated in order to test hypothesis one (H01) to hypothesis four (H04). In model 1, the structural direct effects of each dimension of teacher self-efficacy and job satisfaction are validated. On the other hand, model 2 establishes the direct effects of each dimension of school climate on job satisfaction while model 3 validates how these effects are mediated by each component of the teacher self-efficacy. Models 1 and 2 are found to demonstrate good fits with the following global fits statistics for model 1; χ 2 (155) = 651.39, p < .001, RMSEA = .040 (90% C.I. = 0.037-0.043, CFit = 1.000), CFI = .96, TLI = .95, and SRMR = .032, and for model 2; χ 2 (96) = 409.67, p < .001, RMSEA = .041 (90% C.I. = 0.037-0.045, CFit = 1.000), CFI = .97, TLI = .96, and SRMR = .037. Two non-significant structural paths are found in model 1, and these are deleted before validating model 3. The global fit statistics for model 3 are  Table 2 while Fig. 2 displays only the standardized significant path weights and R 2 values of the predicted variables for model 3.
The results in Table 2 (model 1) show that teachers' self-efficacy in classroom management and teachers' selfefficacy in instruction have significant direct effects on their job satisfaction with work environment while the direct of effect of the 'teachers' self-efficacy in student engagement' on their satisfaction with work environment is not significant. However, there is significant direct effect of teachers' self-efficacy in student engagement on their job satisfaction with profession. Also, the direct effect of teachers' self-efficacy in instruction on their job satisfaction with profession is significant while the direct effect of their 'self-efficacy in classroom management' is not significant. Further, among the three constructs that measure teacher self-efficacy it can be deduced that teachers' selfefficacy in instruction has the greatest effect on the overall job satisfaction of the teachers. This finding could be interpreted to mean that teachers who often craft good questions for students, use variety of assessment strategies and help confused students by giving alternative explanation are most probable of having high teaching job satisfaction in Norway. These results provide detail empirical evidence that confirms the first hypothesis (H01) of the current study.
The results presented in Table 2 (model 2) also reveal at a global level that school climate has significant direct impact on teachers' job satisfaction. Teacher-student relations appear to exert greater effect on teachers' job satisfaction than their perceived disciplinary climate which in parts confirm hypothesis two (H02). This finding could be interpreted to mean that schools whose teachers have good rapports with their students and care for students' well-being, are interested in the opinions of the students and provide extra assistance for students are more probable to have teachers with high job satisfaction. It is also confirmed, as postulated in hypothesis three (H03), that there is a significant negative effect of 'teachers' perceived disciplinary climate' on their job satisfaction. However, this negative effect needs to be interpreted with cautions going by the way the survey questions are framed. The finding confirms that schools whose teachers waste a lot of time getting students ready for lessons by managing disruptive noises, students' interruptions, and unpleasant atmosphere are more probable to have teachers with low job satisfaction.
The results presented Fig. 2 show at a global level that school climate have significant direct effects on teachers' Fig. 2 A structural model of the effects of school climate (TPDC, teachers' perceived disciplinary climate and TSR, teacher-student relations) on job satisfaction (TJSENV, job satisfaction with work environment and TJSPRO, job satisfaction with profession) as mediated by the teacher selfefficacy (TSECLS, self-efficacy in classroom management; TSEINS, self-efficacy in instruction, and TSEENG, self-efficacy in student engagement). All parameter estimates are standardized and significant at p < .05 Zakariya International Journal of STEM Education (2020) 7:10 Page 8 of 12 self-efficacy in classroom management, instruction and student engagement. These can be deduced from the significant regression weights from each dimension of the school climate to the respective dimensions of teacher self-efficacy. It is found that 21% variance associated with 'teacher self-efficacy in classroom management', 11% variance associated with 'teacher self-efficacy in student engagement' and 7% variance associated with 'teacher self-efficacy in instruction' are explained by the school climate-related factors. These percentages of explained variance appear low, however, they are statistically significant which could be indicative of the relevance of other predictor variables that are not accounted for in the current study. The results in Fig. 2 also generally confirm hypothesis four (H04) that teacher self-efficacy mediates the impact of school climate on teacher job satisfaction. This can be deduced from the absence of considerable changes in the standardized regression weights except that the direct impact of 'teacher self-efficacy in instruction' on 'job satisfaction with work environment' observed in Table 2 is no longer significant.
To further substantiate the mediating effect of teacher self-efficacy between the school climate and the teacher job satisfaction, a mediation analysis was conducted. The results show that 'teachers' perceived disciplinary climate' has a significant total effect of − .128 (p < .001) on 'job satisfaction with profession', a significant total effect of − .106 (p < .001) on 'job satisfaction with work environment', and significant total indirect effects of − .034 (p < .001) and − .030 (p = .009) through teacher self-efficacy on 'job satisfaction with profession' and 'job satisfaction with work environment' respectively. It is also found that teacher-student relations have a significant total effect of .227 (p < .001) on 'job satisfaction with profession', a significant total effect of .405 (p < .001) on 'job satisfaction with work environment', and a significant total indirect effect of .020 (p = .004) through teacher self-efficacy on job satisfaction with profession. However, the mediating effect of teacher self-efficacy between teacher-student relations and 'job satisfaction with work environment' is not significant .010 (p = .104). In summary, the results demonstrate evidence of mediating effect of the overall teacher self-efficacy construct between school climate and teacher job satisfaction. Further, 20% variance of 'teacher job satisfaction with work environment' and 9% variance associated 'job satisfaction with profession' are explained by the predictor variables.

Stage 2: Cross-validation of the structural model
In this section, results of the analyses on cross-validation of the structural model 3 established in stage 1 are presented. These analyses are in two steps. The first step involved fitting model 3 simultaneously on data 1 (calibration group, N = 1975) and data 2 (validation group, N = 1976) such that all parameters are freely estimated in the two groups. The second step involves multigroup invariance in which all factor loadings, regression weights and intercepts are constrained equal in the two groups. The results from the first (configural model) and second (constraints model) steps analyses are presented in Table 3. Note. Pchange, change in goodness of indices across the configural and the constraints models.
Zakariya International Journal of STEM Education (2020) 7:10 Page 9 of 12 The results presented in Table 3 (configural model) reveal the extent to which the structural model 3 fits the data from both the calibration and the validation groups simultaneously with freely estimated parameters. There is a slight change between the chi-square values of the two groups (calibration = 1147.257, validation = 1111.656, change = 35.601). However, the GOF indices indicated a good model fit of the two groups. This finding suggests that if the model parameters, e.g. factor loadings and regression weights, are allowed to freely estimate in both calibration and validation groups the overall model exhibit well-fitting statistics.
A review of the results in Table 3 (constraints model) confirms the extents to which factor loadings, intercepts and regression weights are invariance across the calibration and the validation groups. The results suggest a good fit of the model and demonstrate equality of constrained parameters across the two groups. This conclusion is evident from the non-significant MLR chisquare difference test with Satorra-Bentler correction Δ χ 2 ½59 ¼ 47:179; p ¼ :866 coupled with slight or no changes in GOF indices that are within recommended values to establish measurement invariance (Chen, 2007). According to Meredith (1993) classification of multigroup invariance model, these results could be interpreted to be indicative of strong measurement and structural invariance which confirm hypothesis five (H05) of the current study.

Discussion
An understanding of the relations between school climate, teachers' self-efficacy and their job satisfaction has a vital role to play in identifying, recruiting and retention of teachers. This is important for improving teachers' well-being, quality of school management, teaching and learning outcomes. However, previous studies on these relations are either led to inconclusive findings or lack cross-validation of the established models. Thus, findings of the current study provide empirical evidence for a cross-validated structural model that explains the direct/ indirect effects of school climate and teacher selfefficacy on job satisfaction. This model is also found to be invariance in estimated parameters, e.g. factor loadings, regression weights and intercepts, across two independent samples from a population.
In specific terms, it was found that 'teachers' selfefficacy in classroom management', their self-efficacy in instruction and their self-efficacy in student engagement have direct effects on teaching job satisfaction. The effect of 'teachers' self-efficacy in instruction' on job satisfaction appeared to be the greatest among the three dimensions of the teacher self-efficacy. These findings are consistent in parts with the results reported in (Edinger & Edinger, 2018;Skaalvik & Skaalvik, 2014). However, the dimensions of teachers' self-efficacy used in those two studies are slightly different from that of the current study. Nonetheless, the findings of the current study corroborate the two previously reported and seem to contradict that of Shaukat et al. (2019) who reported no significant relationship between teacher selfefficacy and job satisfaction. This disparity between the findings of the current study and that of Shaukat et al. (2019) could be linked to either small sample size (N = 118) involved in their study or some methodological issues (e.g. no account for measurement errors in predictor variables) inherent in the use of analysis of variance (ANOVA).
Another important finding that is revealed in the current study is the direct effect of school climate on teachers' job satisfaction. Both the teacher-student relations and the 'teachers' perceived disciplinary climate' have strong relations with the two dimensions of job satisfaction. That is, schools with conducive environment that is free of disruptive noises and unruly student behaviour, whose teachers have good rapports with their students and care for students' well-being, are interested in the opinions of the students and provide extra assistance for students, are more probable to have teachers with high job satisfaction. Consistent with the findings of Veldman et al. (2013), it is confirmed in the current study that teacher-student relations exert greater effect on job satisfaction than the other dimension of school climate variable. These findings in parts also corroborate the report of strong relations between school climate/ culture and teacher job satisfaction reported in You et al. (2016).
Further, It is found that 'teachers perceived disciplinary climate' has negative effect on job satisfaction just as it was reported in Aldridge & Fraser (2016) and Malinen & Savolainen (2016). However, this negative effect needs to be interpreted with cautions going by the way the survey questions are framed as mentioned before. In fact, the negative effects are positive conceptually because the questionnaire requires the teachers to rate the level of their agreement from 1 (strongly disagree) to 4 (strongly agree) with negatively worded statements, see Table 1. Thus, these findings revealed that schools whose teachers waste a lot of time getting students ready for lessons by managing disruptive noises, students' interruptions and unpleasant atmosphere are more probable to have teachers with low job satisfaction with work environment and profession.
The global effect of school climate on teacher job satisfaction is mediated by teacher self-efficacy. However, empirical evidence is only provided for strong indirect effects of teacher-student relation on 'job satisfaction with profession' through 'self-efficacy in student engagement' and Zakariya International Journal of STEM Education (2020) 7:10 Page 10 of 12 strong indirect effect of 'teacher perceived disciplinary climate' on 'job satisfaction with work environment' through 'self-efficacy in classroom management'. These findings are consistent with the results of most studies on the relations between these constructs (e.g. Aldridge & Fraser, 2016;Malinen & Savolainen, 2016). Finally, a relatively unique finding of this study is the invariance of the model that explains the relations between school climate, teacher self-efficacy and job satisfaction. This model was found to be invariant in factor loadings, intercepts, and regression weights across two independent samples from the same population. It could be interpreted to be evidence of cross-validation that enhance the possibility of getting similar results when the study is replicated. Some limitations of this study are briefly discussed in the next section.

Limitations of the study
There are a number of limitations to the current study. The first limitation concerns the unavailability of the response data that capture 'teacher job satisfaction with target class autonomy' and 'participation among stakeholders' that measure dimensions of job satisfaction and school climate respectively. These data are not included in the public data provided by TALIS 2018 team. Perhaps, the availability of these data would have improved the model results especially the explained variance of the job satisfaction and the teacher self-efficacy that appear low. Another limitation concerns non-feasibility of cross validating the established model of this study across other educational levels, e.g. primary and upper secondary school levels. This would have increased the generalization of the results to a wider range of teachers. However, this is not possible because no national data are available for these levels. Further, it is acknowledged that a threat to the external validity of the findings of the current study could stem from the restriction of the study sample to a Scandinavian country. Perhaps, more refine findings could have been discovered if crosscountry analyses are investigated.

Conclusions
The findings reported in this article have provided empirical evidence for the structural relations between school climate, teacher self-efficacy and job satisfaction. Proper understanding of this relationship is very crucial in identification, recruitment and sustainability of effective teachers. The school climate factors are found to have direct effect on teacher job satisfaction. However, this direct effect is mediated by some dimensions of teacher selfefficacy. A strength of these findings lies in the large sample involved coupled with the robust structural equation modeling approach used unlike mere correlational studies in the literature (e.g. Hosford & O'Sullivan, 2016;Shaukat et al., 2019). Another important contribution of this study to the literature is the cross-validation of the established structural model using independent sample from the same population. This type of structural model cross-validation technique is very important for generalization of results, and it is rarely reported in the literature (Byrne, 2012), perhaps for the obvious reason of lack of enough sample to conduct the analysis. This section is concluded by making two important observations on the general use of TALIS 2018 data for future studies. The first observation concerns the recoding of items TG53C, TG53D, TG53F and TG41B to reflect their negatively worded nature. TALIS 2018 team mentioned in their technical report that these items were reversed coded during the scale constructions and validations. However, what seems apparent is that responses on these items are not reversed coded in the 2018 survey data released to the public. This conclusion is reached because of the negative covariances found between these items and other items in their respective constructs during data screening process. The second observation concerns cross-loading of the item TG53J: All in all, I am satisfied with my job (that was hypothesized to measure 'job satisfaction with work environment') onto 'job satisfaction with profession' construct. This cross-loading is found to significantly improve parsimony and goodness of fits indices of the measurement model when included. It is therefore recommended for researchers who will be using TALIS 2018 data to take note of these observations.