Greenness and job-related chronic stress in young adults: a prospective cohort study in Germany.

OBJECTIVES: We aimed to prospectively study the association between normalised difference vegetation index (NDVI) as a measure of greenness around homes and occupational stress.
SETTING: A population-based cohort in Munich and Dresden cities was followed from age 16-18 years to age 20-23 years (n=1632). PARTICIPANTS: At baseline, all participants attended high-school while at follow-up some had started working and others studying at university. At baseline and in each follow-up, we assigned NDVI based on participants' residential geocoded addresses and categorised it by quartiles. OUTCOME MEASURES: School-related, university-related or job-related self-reported chronic stress was assessed at the two follow-ups by the Trier Scale for Assessment of Chronic Stress using work discontent and work overload as outcomes. We modelled the association employing ordinal generalised estimating equations model accounting for changes in sociodemographics, non-job-related stress, job history and environmental covariates. Stratified analysis by each city was performed.
RESULTS: NVDI at baseline was higher for participants from Dresden (median=0.36; IQR 0.31-0.41) than Munich (0.31; 0.26-0.34). At follow-up, it decreased only for participants in Dresden (0.34; 0.30-0.40). Higher greenness (quartile 4 vs quartile 1) was associated with less work discontent (OR 0.89; 95% CI 0.80 to 0.99) and less work overload (OR 0.87; 95% CI 0.78 to 0.96). In stratified analyses, results were more consistent for Munich than for Dresden.
CONCLUSIONS: Our results suggest that residential green spaces, using the vegetation index as a proxy for exposure, are inversely associated with two types of job-related chronic stress in German young adults transitioning from school to university or working life. © Article author(s) (or their employer(s) unless otherwise stated in the text of the article) 2018. All rights reserved. No commercial use is permitted unless otherwise expressly granted.

This paper is the first study investigating the association of greenness and job-related chronic stress using a prospective study in young adults.

We used a validated instrument to measure the job-related stress as well as stress dimensions outside the workplace.

Our results remained robust when controlling for a wide range of confounding factors like sociodemographic variables, type of employment and non-job-related chronic stress.

We used satellite images from 2009 for the normalised difference vegetation index calculations, which were not temporally aligned with the survey used.

We did not include participants’ commuting times, nor green measurements at work places, making it impossible to control for these covariates.

Some contextual data like ethnic background and social cohesion were not taken into account in our analysis.

Introduction

Some projections estimate that around 66% of the worldwide population is expected to live in urban areas by 2060; in Germany, this proportion is expected to reach 92% by then.1 Moreover, the global number of large cities (ie, 5–10 million inhabitants) will be 65 by the year of 2030, housing approximately 400 million people worldwide.2 Large cities present advantages for economic growth and industrialisation. Also, they promote livability and sustainability and are bringing social and health benefits.3 However, big cities have adverse impacts on health and well-being of their inhabitants.4

Mental health is especially affected in big cities, in part because of reduction of social interaction, excessive commuting and through increments in psychological stressors.5 6 Psychological distress has been identified as an essential public health and economic problem in large cities because it leads to work absence, lower productivity or early retirement.7 Susceptible populations like students and young adults in transition from school to job/university life reported increments in job-related chronic stress.8 In line with this, evidence has shown significant mental health benefits of residing in cities with more natural environments through restorative effects on psychological health.9–17

Mechanisms explaining the restorative effect of a green environment are attributed mainly to physical activity, social contact and deliberately seeking environments to recover from demanding situations and tasks.10 12 18–21 Restorative niches were associated with emotional well-being, as well as increments of project management capabilities in adolescents transitioning from junior to secondary school.22 Additionally, few studies have studied associations between the access to green environments and job satisfaction.23–25 Researchers have shown the restorative effect of greenness mainly at schools or workplaces.19 26–28 These studies showed that being exposed to a green environment at the workplace/school is beneficial for workers or students. Despite these potential benefits, most of the office workers and students do not go outdoors during the job/study day, mainly because of a perception of having many tasks to do and a job/studying culture that does not include outdoor behaving.27 29

So far, there are few studies relating greenness around the home environment to job-related chronic stress, especially in longitudinal settings.16 30 31 No previous research has analysed the association between greenness at home and job-related stress development in young people in the transition from school to university or working life. Our aim was therefore to investigate the association between normalised difference vegetation index (NDVI) as a measure of greenness around homes and two different job/study-related chronic stress measurements: work discontent and work overload.32 33 Data were drawn from a population-based cohort in two major German cities, specifically Munich and Dresden, surveying young adults in the transition from high school to university or working life.

Methods

We analysed data of the SOLAR (Studie in Ost- und Westdeutschland zu beruflichen Allergierisiken) I and SOLAR II studies.34 These studies were aimed at studying the course of respiratory diseases and atopy in symptomatic and non-symptomatic children and young adults. Also, occupational risk factors were assessed to investigate associations among occupational factors, stress and the course of respiratory diseases. SOLAR II is the second follow-up of the German phase II of the International Study on Asthma and Allergies in Childhood (ISAAC II), a multicentre study planned to assess the prevalence of asthma and allergies with participants from elementary school students.35 The German branch of the ISAAC II was carried out in 1996/1997 in Munich and Dresden. In these cities, it included 7498 participants at age 10 years (fourth grade, age range 9–11 years). From August to January 2003, 3053 of these participants answered the questionnaire of the first follow-up of the study (SOLAR I; age range 16–18 years). From August 2007 to November 2008, 2051 of the participants agreed to participate in SOLAR II (age range 19–24 years).

In the present analysis, we considered questionnaire data from SOLAR I and II, as only these two surveys contained data on chronic stress (figure 1). Besides chronic stress items, the written questionnaires contained validated items on sociodemographics, type of job, occupational diseases and physical activity. More details on the study methods are given by Heinrich et al.34

Figure 1

Flow chart of follow-up from International Study on Asthma and Allergies in Childhood (ISAAC) phase II to SOLAR II including normalised difference vegetation index (NDVI) measurements, job-related and non-job-related stress variables and total participants in each follow-up and the current analysis. Excluded participants: n=318 ever worked previously to SOLAR I, n=21 no educational status information, n=24 more than two Trier Scale for Assessment of Chronic Stress (TICS) items missing and n=56 residence out of the study areas.

In the current analyses, we excluded participants who had ever worked before SOLAR I (n=318), those without information on educational status (n=21) and those with more than two items missing in one of the Trier Inventory for Chronic Stress (TICS)32 33 scales (n=24). We geocoded the families’ or participants’ most recent addresses in SOLAR I and SOLAR II. We then calculated the environmental covariates and matched them to the individual questionnaire information. We excluded 56 participants because they moved out of the study areas (figure 1).

All participants or their legal guardians provided written informed consent.

Variable definition

Job-related stress

For this study, it was essential to use an instrument applicable to the school, university and working environments. TICS is a well-established instrument that includes scales of job-related chronic stress and stress outside the workplace dimensions. These scales were selected using the model of health33 and validated using confirmatory analysis in a representative sample.32 Hence, it is well suited to study the change of stress from school life to working or university life. We used two job-related TICS subscales as chronic stress outcomes: work discontent (eg, "Satisfaction from the work that I have to perform daily") and work overload (eg, "Too little time to execute my daily tasks"). Responses were assessed on a 5-point Likert scale from ‘never’ (0 points) to ‘very often’ (4 points). Total scores were formed as sum of all particular items in each subscale. Based on the recommendations of the scales’ authors,33 up to two missing items per subscale were accepted otherwise, the subscale was considered as missing. Each subscale was categorised based on the frequency of stressful situations, as ‘low’ (≤median), ‘average’ (above median to median+1 SD) and ‘high’ (≥1 SD from the median).36 We used the median and SD obtained in SOLAR I to categorise the outcomes in SOLAR I and SOLAR II.

Greenness of the home environment

We geocoded participants’ addresses in SOLAR I and SOLAR II. Based on the geocoding, the exposure to natural and green areas around each home address was estimated using NDVI, a satellite image-based vegetation index. is the ratio of differences between the near-infrared region (NIR) and red reflectance (RR) to the sum of these two measures, and it ranges between −1.0 (water) and 1.0 (dense green vegetation), with zero indicating sand, snow or rocks.37 Choosing Landsat 5 Thematic Mapper (http://earthexplorer.usgs.gov/) cloud-free satellite images, an average NDVI was obtained using a 30 m by 30 m resolution in a 500 m radius around home addresses. The 500 m radius was used because it is considered to reflect better greenness directly accessible outside each home.38 For the Munich area (including the city of Munich and an adjacent region of Upper Bavaria), we used two images from 31 August in 2009 and merged them to cover the entire study area. In Dresden (including the Saxony area), four images were combined to cover the study area. As no data for a single day were available, we used images from 24 and 31 August in 2009. Finally, the exposure was categorised using survey and city-specific quartiles to account for distributional differences between city and study period.39

Covariates

Sociodemographics

As sociodemographic variables, we included: sex (male vs female), highest educational status reported at SOLAR I and SOLAR II (elementary education vs secondary education, advanced technical or higher education) and if the participant had children by SOLAR II (no vs yes). To account for physical activity, we used the WHO recommendations40 considering the age-specific cut-off point for SOLAR I (≥4 hours/week) and SOLAR II (≥2 hours/week).

Non-job-related chronic stress

TICS also measures the stress outside the workplace/university environment using the following four subscales:

social overload (eg, "I quarrel with others because I do not behave the way others expect me to");

lack of social recognition (eg, "The experience that other people have no trust in me");

chronic worrying (eg, "Times when I am not able to suppress my worries");

stressful memories (eg, "Recurrent memories of failures").

We measured these non-job-related chronic stress scales in SOLAR I only (figure 1) because they are considered to be persistent over time8 and used them as covariates to control for self-perceived stressful circumstances outside the work or study environment. Categorisation of the non-job-related chronic stress scales was done using the same method as for job-related stress.

Current status and job type

At SOLAR II, we categorised the occupational status as follows: employee (reference), university student, vocational trainee (in a dual training programme), unemployed, self-employed and other (ie, on maternity leave or being work disabled).8 34 In addition, participants reported any jobs they had ever held up until SOLAR I and between SOLAR I and SOLAR II. Two trained persons coded these jobs (including regular employment, student job or internship) according to the International Standard Classification of Occupations (ISCO-88).41 Following the ISCO classification, we assigned each working participant to one of five job groups: clerks, professionals and technicians, healthcare professionals (ie, with direct patient contact), plant machine operators and elementary occupations (ie, routine tasks using mainly handheld tools and some involving physical effort).

Environmental covariates

We used the European Environmental Agency databases42 to link participants’ residential addresses to the following spatial covariates: distance to sports facilities,43 distance to nearest urban green space44 (ie, either a garden, a park, a cemetery or a plant nursery), distance to nearest lake or river11 and the percentage of tree coverage45 (ie, the percentage of flat ground covered by woody vegetation over a 5 m height). Additionally, we included the proportion of soil sealing derived in 2000 m buffer for the year 2009 as a measure of urbanisation degree. Finally, we used the NDVI measurement in a 500 m buffer at ISAAC II as a continuous variable to control for long-term greenness exposure.42

Statistical analysis

We employed generalised estimating equations (GEE) models for ordinal outcomes with an exchangeability correlation structure.46 Applying this strategy, we assessed the relationship between greenness and job-related stress, meanwhile controlling for the outcomes temporal correlation.46 47 GEE methods produce population-averaged estimates, that is, they describe changes in the population mean based on changes in covariates.47 48 In the present analysis, GEE estimates tell us how much the studied outcomes (work overload or work discontent) would change on average in the population for each NDVI quartile increase controlling for covariates. Positive values thus mean an increase in chronic stress over time.

We included all mentioned covariates in the final models. Given the possible interaction between city and greenness,39 we stratified the analyses by city (Munich and Dresden). Furthermore, we estimated the final models using only participants who never changed their home addresses during the whole study period (n=675). Additionally, we restricted the analyses to only students (n=845) and only workers, that is, participants that reported being employees or vocational trainees (n=670). Finally, we performed a mediation analysis using the approach suggested by Schluchter49 to test the hypothesis that physical activity could mediate the association between greenness and job-related stress.

Using complete-case analyses in the presence of missing data might bias the results.50 Therefore, we constructed five imputed datasets using multiple imputations by chained equations.51 Using Rubin’s rules,52 we obtained the combined adjusted ORs and CIs. We presented the comparison between complete cases and multiply-imputed estimates in the online supplementary material.

Statistical analyses were performed in R V.3.3.153 and the geographical calculations using ArcGIS 10.0 (Environmental Systems Research Institute, 2012) joint with the Geospatial Modelling Environment software (Spatial Ecology).

Patient and public involvement

Patients and public were not involved in the design of the study.

Results

The number of participants in both cities was similar (869 in Munich and 763 in Dresden). Surrounding average greenness in a 500 m buffer for SOLAR I was higher in Dresden (median=0.36; IQR 0.31–0.41) than Munich (0.31; IQR 0.26–0.34). In SOLAR II, NVDI decreased for participants in Dresden (0.34; IQR 0.30–0.40) while it remained the same for participants living in Munich (0.31; IQR 0.26–0.34) (see online supplementary table S1 and figure S1).

Comparing the covariates across the greenness quartiles, we found small difference between quartiles (table 1): participants with children tended to live in places with greenness values in the lower quartiles of the greenness distribution, and those with higher education and students in the upper quartiles. Likewise, the percentage of physically active people was higher with increasing level of greenness around the participants’ houses. Based on the p values, we did not find statistically significant differences in Munich between the distance to sports facilities or distance to the nearest urban green space and NDVI quartiles; however, all other environmental covariates showed statistically significant differences with respect to NDVI quartiles. In Dresden, all environmental covariates were associated with NDVI (table 1).

Table 1

Sociodemographic characteristics, type of job, psychological scores, environmental variables at baseline by city-specific NDVI quartile and city

	Total	Munich, n=869				P values	Dresden, n=763				P values
		Centre-specific NDVI quartile					Centre-specific NDVI quartile
		Q1	Q2	Q3	Q4		Q1	Q2	Q3	Q4
Individual characteristics, n (%)
Sex, woman	961 (58.9)	127 (59.1)	122 (56.7)	32 (61.7)	119 (55.3)	0.56	111 (60.3)	110 (59.8)	111 (60.7)	112 (60.9)	0.99
Having children* (NA=20), yes	51 (3.1)	5 (2.4)	6 (2.8)	2 (<1)	3 (1.4)	0.26	12 (6.8)	6 (3.4)	7 (4.0)	6 (3.4)	0.17
Education* (NA=6)
Elementary	8 (<1)	1 (<1)	2 (<1)	2 (<1)	1 (<1)	<0.01	0	1 (<1)	0	1 (<1)	0.08
Secondary	419 (25.8)	73 (34.3)	55 (25.9)	53 (25.2)	47 (22.4)		54 (30.7)	41 (23.4)	48 (27.3)	33 (18.9)
Advanced technical	194 (11.9)	31 (14.6)	32 (15.1)	32 (15.2)	21 (10.0)		17 (9.7)	14 (8.0)	18 (10.2)	21 (12.0)
Higher	1005 (61.8)	108 (50.7)	123 (58.0)	123 (58.6)	141 (67.1)		105 (59.7)	119 (68.0)	110 (62.5)	120 (68.6)
Physical activity* (NA=33), >2 hours/week	511 (32.0)	55 (26.1)	76 (35.7)	74 (35.1)	70 (32.9)	0.69	48 (27.4)	57 (31.7)	50 (27.8)	64 (35.6)	0.05
Current status* (NA=5)
Employed	365 (22.4)	52 (24.5)	44 (20.8)	42 (20.0)	43 (20.4)	0.35	49 (27.8)	35 (20)	43 (24.4)	38 (21.6)	0.10
Students	845 (51.9)	97 (45.8)	110 (51.9)	118 (56.2)	124 (58.8)		74 (42.0)	97 (55.4)	82 (46.6)	96 (54.5)
Apprentice	305 (18.7)	47 (22.2)	49 (23.1)	41 (19.5)	35 (16.6)		32 (18.2)	23 (13.1)	39 (22.2)	31 (17.6)
Unemployed	59 (3.6)	8 (3.8)	4 (1.9)	5 (2.4)	2 (<1)		10 (5.7)	12 (6.9)	7 (4.0)	8 (4.5)
Other	53 (3.3)	8 (3.8)	5 (2.4)	4 (1.9)	7 (3.3)		11 (6.2)	8 (4.6)	5 (2.8)	3 (1.7)
Type of job, n (%)
Clerks	159 (9.7)	19 (8.8)	16 (7.4)	21 (9.8)	25 (11.6)	0.24	18 (9.8)	20 (10.9)	14 (7.7)	21 (11.4)	0.64
Professionals and technicians	318 (19.5)	41 (19.1)	48 (22.3)	39 (18.2)	40 (18.6)	0.74	41 (22.3)	32 (17.4)	31 (16.9)	37 (20.1)	0.92
Health professionals	202 (12.4)	18 (8.4)	22 (10.2)	27 (12.6)	31 (14.4)	0.24	23 (12.5)	21 (11.4)	30 (16.4)	26 (14.1)	0.80
Plant machine operators	248 (15.2)	31 (14.4)	36 (16.7)	34 (15.9)	22 (10.2)	0.38	34 (18.5)	29 (15.8)	21 (11.5)	37 (20.1)	0.30
Elementary occupations	90 (5.5)	10 (4.7)	10 (4.7)	8 (3.7)	5 (2.3)	0.65	14 (7.6)	14 (7.6)	10 (5.5)	17 (9.2)	0.85
Psychological variables, n (%)
Lack of social recognition (NA=14)
Low	739 (45.7)	91 (42.9)	102 (47.4)	112 (53.1)	98 (46.0)	0.09	85 (46.7)	73 (40.1)	82 (44.8)	76 (41.5)	0.35
Average	656 (40.5)	84 (39.6)	85 (39.5)	78 (37.0)	90 (42.3)		69 (37.9)	84 (46.2)	74 (40.4)	79 (43.2)
High	223 (13.8)	37 (17.5)	28 (13.0)	21 (10.0)	25 (11.7)		28 (15.4)	25 (13.7)	27 (14.8)	28 (15.3)
Social overload (NA=14)
Low	638 (39.4)	76 (35.8)	73 (34.4)	89 (42.2)	76 (35.7)	0.25	78 (42.4)	80 (43.5)	74 (40.7)	73 (39.9)	0.49
Average	697 (43.1)	88 (41.5)	99 (46.7)	79 (37.4)	93 (43.7)		76 (41.3)	77 (41.8)	79 (43.4)	90 (49.2)
High	283 (17.5)	48 (22.6)	40 (18.9)	43 (20.4)	44 (20.7)		30 (16.3)	27 (14.7)	29 (15.9)	20 (10.9)
Chronic worrying (NA=16)
Low	822 (50.9)	109 (50.9)	106 (50.2)	99 (46.9)	104 (48.8)	0.29	100 (54.9)	97 (53.0)	92 (50.5)	92 (50.3)	0.40
Average	461 (28.5)	70 (32.7)	55 (26.1)	50 (23.7)	61 (28.6)		48 (26.4)	58 (31.7)	55 (30.2)	56 (30.6)
High	333 (20.6)	35 (16.4)	50 (23.7)	62 (29.4)	48 (22.5)		34 (18.7)	28 (15.3)	35 (19.2)	35 (19.1)
Stressful memories (NA=12)
Low	852 (52.6)	111 (52.1)	112 (52.1)	101 (47.6)	99 (46.3)	0.81	106 (58.2)	105 (57.7)	101 (55.8)	93 (50.8)	0.30
Average	423 (26.1)	56 (26.3)	51 (23.7)	56 (26.4)	65 (30.4)		42 (23.1)	44 (24.2)	43 (23.8)	55 (30.1)
High	345 (21.3)	46 (21.6)	52 (24.2)	55 (25.9)	50 (23.4)		34 (18.7)	33 (18.1)	37 (20.4)	35 (19.1)
Environmental variables, mean (SD)*
Distance to sports facilities (km)	0.34 (0.26)	0.37 (0.25)	0.36 (0.21)	0.33 (0.19)	0.35 (0.26)	0.31	0.45 (0.34)	0.24 (0.20)	0.26 (0.23)	0.38 (0.33)	<0.01
Distance to nearest urban green space (km)	0.28 (0.27)	0.25 (0.18)	0.25 (0.19)	0.27 (0.39)	0.23 (0.22)	0.15	0.36 (0.27)	0.35 (0.21)	0.25 (0.22)	0.34 (0.40)	<0.01
Distance to nearest lake or river (km)	1.43 (0.91)	1.26 (0.52)	1.58 (0.90)	1.69 (1.02)	1.82 (1.15)	<0.01	0.82 (0.60)	1.28 (0.75)	1.32 (0.78)	1.55 (0.99)	<0.01
Proportion of tree coverage in 500 m	2.5 (4.1)	0.75 (0.61)	1.09 (1.24)	1.48 (1.86)	4.54 (5.90)	<0.01	1.60 (1.02)	1.50 (1.12)	2.23 (1.96)	7.30 (7.45)	<0.01
Proportion of soil sealing	49.1 (18.6)	67.7 (10.1)	53.7 (11.7)	46.7 (11.8)	40.2 (14.2)	<0.01	60.8 (15.0)	53.1 (15.1)	43.5 (14.7)	24.5 (14.5)	<0.01
Greenness at ISAAC II (500 m buffer)	0.31 (0.08)	0.30 (0.06)	0.32 (0.06)	0.35 (0.06)	0.37 (0.08)	<0.01	0.22 (0.05)	0.28 (0.05)	0.32 (0.04)	0.35 (0.07)	<0.01

Total sample, n=1632.

*Calculated in the second follow-up (SOLAR II).

ISAAC, International Study on Asthma and Allergies in Childhood; NDVI, normalised difference vegetation index; NA, missing data; Q, quartile.

Prevalence of high levels of work discontent and work overload decreased by increasing level of greenness in a buffer of 500 m around the home, especially for SOLAR II (table 2). Results were confirmed when we took changes over time into account as shown by ORs <1 for work discontent (adjusted OR comparing the fourth to the first quartile of greenness: OR 0.89; 95% CI 0.80 to 0.99) and work overload (OR 0.87; 95% CI 0.78 to 0.96) (figure 2A).

Table 2

Prevalence of job-related chronic stress outcomes by city-specific NDVI quartile in SOLAR I and SOLAR II

	Munich, n=869				Dresden, n=763
	Centre-specific NDVI quartile				Centre-specific NDVI quartile
	Q1	Q2	Q3	Q4	Q1	Q2	Q3	Q4
SOLAR I
Work discontent (NA = 14)
Low	103 (48.4)	108 (50.9)	108 (50.7)	102 (48.1)	111 (61.0)	95 (52.2)	100 (55.2)	103 (56.3)
Average	56 (26.3)	60 (28.3)	62 (29.1)	69 (32.5)	46 (25.3)	57 (31.3)	55 (30.4)	51 (27.9)
High	54 (25.4)	44 (20.8)	43 (20.2)	41 (19.3)	25 (13.7)	30 (16.5)	26 (14.4)	29 (15.8)
Work overload (NA=13)
Low	116 (55.8)	114 (53.5)	125 (58.7)	118 (55.4)	104 (57.1)	97 (52.7)	102 (56.0)	92 (50.5)
Average	56 (26.9)	63 (29.6)	53 (24.9)	55 (25.8)	50 (27.5)	54 (29.3)	54 (29.7)	60 (33.0)
High	36 (17.3)	36 (16.9)	35 (16.4)	40 (18.8)	28 (15.4)	33 (17.9)	26 (14.3)	30 (16.5)
SOLAR II
Work discontent (NA=14)
Low	126 (59.4)	135 (64.0)	152 (71.7)	141 (66.8)	108 (62.4)	126 (72.4)	125 (71.4)	127 (72.2)
Average	47 (22.2)	43 (20.4)	34 (16.0)	46 (21.8)	40 (23.1)	27 (15.5)	33 (18.9)	30 (17.0)
High	39 (18.4)	33 (15.6)	26 (12.3)	24 (11.4)	25 (14.5)	21 (12.1)	17 (9.7)	19 (10.8)
Work overload (NA=13)
Low	99 (46.7)	110 (52.1)	114 (53.8)	105 (50.0)	90 (51.4)	100 (56.8)	101 (58.0)	98 (56.0)
Average	59 (27.8)	50 (23.7)	56 (26.4)	64 (30.5)	48 (27.4)	48 (27.3)	38 (21.8)	47 (26.9)
High	54 (25.5)	51 (24.2)	42 (19.8)	41 (19.5)	37 (21.1)	28 (15.9)	35 (20.1)	30 (17.1)

Total sample, n=1632.

NA, missing data; NDVI, normalised difference vegetation index; Q, quartile.

Figure 2

City-specific associations between greenness in a 500 m buffer and work discontent and work overload. Ordinal generalised estimating equation models adjusted for sex, having children, physical activity, education, current status, type of job, psychological variables and environmental variables using complete cases in selected subpopulations: (A) complete cases for the combined population (n=1430), Munich (n=779) and Dresden (n=651). (B) Complete cases that never moved in the combined population (n=629), Munich (n=443) and Dresden (n=186). NDVI, normalised difference vegetation index; Q, quartile.

Restricting the study population to participants who did not move between ISAAC II and SOLAR II did not affect the results (figure 2B). Stratifying by city, results were more consistent for Munich than for Dresden. In students (figure 3A), the OR for work overload was <1 for all quartiles compared with the first quartile of NDVI. For workers (figure 3B), no statistically significant difference for work overload and work discontent was seen by NDVI quartiles. Stratifying results for city, results for students followed the same pattern when restricting the analyses for Munich while for work discontent they did not reach the level of statistical significance for students from Dresden. For workers, we did not see any statistically significant associations for participants from Munich while work overload decreased with increasing NDVI quartile in Dresden (figure 3B).

Figure 3

City-specific associations between greenness in a 500 m buffer and work discontent and work overload. Ordinal generalised estimating equation models adjusted for sex, having children, physical activity, education, current status, type of job, psychological variables and environmental variables using complete cases in selected subpopulations: (A) complete cases students only for the combined population (n=845), Munich (n=463) and Dresden (n=382). (B) Complete cases workers only in the combined population (n=670), Munich (n=359) and Dresden (n=311). NDVI,  normalised difference vegetation index; Q, quartile.

In the mediation analysis, we found a mediated proportion of physical activity on stress of 1.46% (95% CI −10.2% to 5.6%) for work discontent, and a 0.13% (95% CI −0.6% to 3%) for work overload. Unadjusted results were similar to the adjusted ones (see online supplementary table S2). Likewise, complete-case analyses were not different from imputed results (see online supplementary table S3).

Discussion

In our cohort of young adults in the transition from school to working/university life, we found that more greenness around the place of living corresponds to lower levels of job-related stress. For work discontent, our results suggest a linear, inverse dose-response pattern across quartiles while for work overload we observed a ‘J-shape’ association. Our results remained robust when controlling for a wide range of confounding factors and with only small differences when stratifying for various subpopulations. These findings thus contribute to an improved understanding of the mental health benefits of green environments, especially in young adults transitioning from school to the university or working life. Our results expand and corroborate the finding of previous studies on the benefits of being exposed to greenness on mental health among different populations.17 29 31 54 55

Greenness exposure brings benefits on mental health. It could operate through increasing opportunities for physical activity and enhanced access to recreational and sports facilities.56–58 Therefore, we controlled for physical activity and distance to sports facilities which did not change the results. Some researchers suggested that physical activity may act as a mediator between greenness and mental health.6 17 However, this was not confirmed in our mediation analysis. Therefore, we hypothesise that the described associations between greenness and stress might be attributed to mechanisms beyond physical activity, which is in line with other studies.18 54

In our study, higher exposure to NDVI levels was inversely related to reporting stress at school, university or workplace, indicating that a favourable environment is an essential resource for recreation and recovery, setting off or balancing out stress related to school, university or workplace. The association remained stable when adjusting for type of job/being a student/being employed so that such potential differences do not explain the observed associations. Differences in the estimated associations between Munich and Dresden may result from different distributions of other factors, like socioeconomic status among both cities, which could affect the susceptibility to greenness exposure and job-related chronic stress.17 39 59

The ‘J-shape’ observed for work overload in our study was also seen in an Australian life-course study using GHQ-12 as a mental health outcome.31 We believe that participants in our cohort who were living in the highest greenness quartile have the most extended commuting times between their homes and place of study or work. Commuting plus work/studying may result in increased mental burden measured by work overload. Unfortunately, we did not ask participants for their commuting times. Future studies on this subject should include commuting times as a potential confounder. As suggested by other authors,18 29 30 58 60 we included environmental variables such as tree coverage, distance to sports facilities, distance to rivers, etc as potential confounders. However, removing those from the adjusted models did not change the results (data not shown).

Among the strengths of our study is its prospective design, which allows to differentiate the temporal sequence between the exposure to greenness and job-related chronic stress. Our questionnaire focused on asthma and allergies so that participants were not aware of the hypotheses of the current paper. Therefore, we expect differential misclassification of the outcomes to be limited. Standardised, objective and blinded assessments of greenness at each follow-up point, from good satellite resolution images are further strengths of our study. Using NDVI as a proxy for exposure to green spaces accounts for all kinds of vegetation, and it can promote greener landscape views and leisure activities, which has shown stress recovering benefits in workers.13 28 56

Among the limitations of our study are the use of satellite images for the NDVI calculations, which were not temporally aligned with the SOLAR I surveys. Nevertheless, we expected spatial greenness to remain stable over several years.58 We could not measure exposure to greenness at work because for ethical issues we could not assess the address of the workplace. This made it also impossible to calculate the commuting greenspace of participants. The shortcomings of our study are common, and almost no study in the field has obtained information other than the residential address for greenspace assessment.60 61 A distinctive feature of our study compared with similar ones is that job-related stress was assessed using a validated scale, instead of using general stress and mental health. Residential greenspaces could improve coping mechanisms to job-related chronic stress in different ways, through viewing, physical activity and social interaction.12 18 Inclusion of additional sources of greenness deserves more attention in future studies on greenness and job-related chronic stress. Our data draw on self-reports of job-related chronic stress, which may be biased by personality or other reporting bias. Moreover, we did not consider whether the participants answered during or after the exam period. Nevertheless, in SOLAR I, participants received the questionnaires between August 2002 and January 2003, which is a period without a high academic load in Germany. In SOLAR II, we sent the questionnaires between August 2007 and November 2008, this means, over a more than a 1 year period. Therefore, SOLAR II also included periods of higher academic load for students. We expected random fluctuations in stress measurement in general, but we do not believe that our results were influenced by this non-differential misclassification because our instrument measures chronic stress rather than study load or workload.33 Job-related chronic stress may be influenced by several personal, social and cultural characteristics. Accordingly, we adjusted all estimates for a set of predefined variables such as sex, education, having children and non-job-related chronic stress measurements. In our study, we did not collect contextual data such as ethnic background, cultural conceptions and social cohesion due to legal and logistical reasons. Those suffering from job-related stress at follow-up might have been less likely to participate. However, we do not assume that this depended on our exposure, so potential selection bias should be limited.

Conclusion

Our findings suggest that residential green spaces, measured using the vegetation index, is inversely associated with job-related chronic stress in German young adults transitioning from school to university or working life. Further studies on this topic could further contribute to improve urban planning or to develop recommendations for health promotion through favourable living and working spaces.