Urban green space and happiness in developed countries

Urban green space is thought to contribute to citizen happiness by promoting physical and mental health. Nevertheless, how urban green space and happiness are related across many countries with different socioeconomic conditions has not been explored. By measuring the urban green space score (UGS) from high-resolution satellite imagery of 90 global cities covering 179,168 km2 and 230 million people in 60 developed countries, we find that the amount of urban green space and GDP are correlated with a nation’s happiness level. More specifically, urban green space and GDP are each individually associated with happiness. Yet, only urban green space is related to happiness in the 30 wealthiest countries, whereas GDP alone can explain happiness in the subsequent 30 countries in terms of wealth. We further show that the relationship between urban green space and happiness is mediated by social support and that GDP moderates this relationship. These findings corroborate the importance of maintaining urban green space as a place for social cohesion to support people’s happiness. Supplementary Information The online version contains supplementary material available at 10.1140/epjds/s13688-021-00278-7.


Introduction
The advantages of urban green space for public health and urban planning have been of great interest in recent years. Green spaces such as parks, gardens, street trees, riversides, and even private backyards facilitate physical activity, social events, mental relaxation, and relief from stress and heat, leading to direct and indirect benefits for mental and physical health and well-being [1][2][3]. Thus, worldwide policy changes and efforts have been made to provide more urban green space to create sustainable and comfortable living environments [4].
Urban green space and happiness are known to have an implicit positive correlation. Although this association is still unclear, there are at least three main reasons that relate green space with health, well-being, and happiness [5]: an innate human emotional affiliation with nature [6], reduced environmental "bads" [7,8], and facilitating behaviors beneficial to physical and mental health. Notably, five pathways through which greenery might have beneficial effects have been reported: relieving stress, stimulating physical activity, facilitating social interactions, generating aesthetic enjoyment, and facilitating a sense of shelter from and adjustment to environmental stressors [3,9,10]. Studies have suggested that the same pathways exist in numerous countries [11]. While various socioeconomic factors are related to happiness [12], social interaction facilitation has been confirmed with strong evidence. Studies [13,14] have shown that open green space promotes social cohesion by providing places for social contact; people can naturally encounter neighbors in local green spaces while walking dogs, gardening, and having outdoor parties, which enhances community engagement. Moreover, larger green areas such as parks can hold public events and activities, enabling social mixing between communities.
The amount of urban green space can be captured mainly by three kinds of measurements: qualitative ratings of observers [9,15], national land-use and land-cover databases [5,16,17], and geographic information system (GIS) techniques. Among these measurements, GIS techniques are the most recently developed method. One example is utilizing the normalized difference vegetation index (NDVI), a vegetation index computed from Landsat series satellite images (30 m resolution) [10,11,18]. Researchers such as Tsai et al. [19] introduced multiple landscape metrics based on GIS and showed a strong association between green space and mental health in U.S. metropolitan areas. These studies assume that the distance from an individual's residence to the nearest green space has associations with health data [3,20]. The green space level was then measured as the fraction of areas with NDVI values above a certain threshold (e.g., 0.2 to 0.4 for sparse vegetation and 0.6 for highly dense vegetation) [21]. However, this method raises the question of how to set an appropriate NDVI threshold for global cities.
Despite the rich literature on the mental benefits of green space, research is limited in terms of global-scale comparative research [2]. First, the analytical settings are based on a limited number of Western countries [10]; most of these studies have been conducted in the United States [18,19] and Europe [3,11]. Moreover, only a few are based on multi-country settings that enable comparative analysis [22]. As a result, it is unclear whether the association between green space and mental health is robust in developing countries or only in developed countries. The main limitation arises because there is no global medical dataset providing reliable and standardized mental health surveys from different countries. Moreover, no studies have established which green space measurement is appropriate for analysis across countries. Various methods of measuring green space -questionnaires, qualitative interviews, satellite images, Google Street View images, and even smartphone technology [23] -still rely on individual-level measurements (e.g., calculating the greenery level around residential buildings) and hence are not scalable to the global level. This paper presents a new way to analyze the effects of green space on happiness at the planetary scale, incorporate socio-economic contexts of different countries, and achieve robust results. First, we measure the amount of urban green space from high-resolution satellite images for different countries by developing a globally comparable green space metric. Our metric based on the total NDVI of built-up areas enables this comparison as it does not require an arbitrary threshold that varies for different regions. It also overcomes the limitations of official statistics based on national land-use land cover data that tend to have different criteria by countries and often include only official parks and open space. Our analysis on high-resolution (10 m) Sentinel-2 satellite images provides more accurate information on urban green space than previous studies on Landsat series images (i.e., a resolution of 30 m) [3,10,11,18].
Next, this study uses selected happiness scores from the World Happiness Report [24], which provides reliable and standardized data on multiple countries' happiness and allows comparisons among nations. As happiness is a criterion of how people feel about the course of their life, it is interconnected with various socioeconomic variables. From the perspective that economic studies distinguish between emotional well-being and life satisfaction [25], we focus on the impact of green space on life satisfaction. Specifically, we study this relationship in developed countries with the highest Human Development Index (HDI), where green environments in cities are considered more important for life satisfaction.
Using these datasets from satellite imagery, we explore the relationship between urban green space and happiness globally. Additionally, we identify conditional indirect effects by national wealth and social support by employing a moderated mediation regression model on socioeconomic indicators.

Urban green space and happiness in countries
We examine the global relationship between urban green space and happiness in 60 developed countries ranked by the Human Development Index. Note that the term happiness represents how people feel about the course of their life [24]. Using the Sentinel-2 satellite imagery dataset, we define each country's urban green space score (UGS) as the logarithmic total vegetation index per capita in the most populated cities (i.e., those that include at least 10% of the national population). Among the various vegetation indices available such as the NDVI, soil-adjusted vegetation index (SAVI) and two-band enhanced vegetation index (EVI2), the NDVI [26] is used based on the robustness of the results for different tested indices. The happiness score and the gross domestic product based on purchasing power parity (GDP (PPP)) per capita of each country are from the World Happiness Report [24] and the International Monetary Fund (IMF) estimation [27], respectively (see the Methods section and Additional file 1 for details). Figure 1(a) shows an overall view of urban green space and the happiness of countries worldwide. This map highlights regional differences in the green space distribution due to climate; countries near the Equator in tropical climates have relatively high UGSs, while countries located in the 20-30 latitude range have exceptionally low UGSs due to the dry climate. The UGS increases with latitude in higher-latitude regions. On the other hand, Northern and Western European and North American countries display relatively high happiness scores. Western Asian countries also show relatively high happiness with a low UGS, indicating that the relationship between happiness and green space is not trivial. Figure 1(b)-(d) shows the distributions of happiness, UGSs, and log-GDP, and they all show unimodal distributions, which is appropriate for linear regression analyses. Note that the probability distributions of GDP per capita converge to a normal distribution after logarithmic scaling. Our comparison of several green space measures shows that the logarithmic NDVI per capita is most suitable for the following analysis in terms of its distribution and explanatory power. Therefore, we choose the logarithmic NDVI per capita as the primary green space indicator in this research (see Additional file 1). We also use the logarithmic GDP per capita (PPP) (hereinafter referred to as the log-GDP) as a measure of the wealth of the country, as noted in the Happiness Report [24].
As per-country wealth is an important indicator of citizens' quality of life, wealth (i.e., log-GDP) should be considered when analyzing urban green space and happiness. Our  Table 1 The regression analysis for the UGS, happiness and log-GDP. The values denote the regression coefficients and the confidence intervals of each independent variable with its significance (i.e., * * * p < 0.01; * * p < 0.05; * p < 0.1). The regression model (1-3), model (4-6) and model (7)(8)(9) are examined for the data of all countries, the lower 30 countries and the top 30 countries ranked by GDP, respectively regression analysis finds that the UGS, together with log-GDP, is related to happiness. We make new observations from Table 1. Although the UGS is not substantially correlated with happiness in the simple linear regression (i.e., model (2)), the multilinear model with log-GDP (i.e., model (3)) has a substantial increase in prediction ability compared to the simple regression analysis on log-GDP (i.e., model (1)). Therefore, urban green space adds explanatory power to the correlation between wealth and happiness across countries. The regression analyses with other green space-variant measures further confirm the robustness of this result, showing a substantial increase in the adjusted R-squared value when the UGS is included in the regression. Specifically, the UGS based on the logarithmic NDVI per capita shows the best regression performance. Moreover, this relationship between the UGS and happiness is robust for the control variables, including life expectancy, health expenditure, unemployment, gender inequality, and education (see Additional file 1 for the results for the different measures).

Urban green space is effective in rich countries
Our results show that happiness is correlated with urban green space and the GDP of a country. However, is this green space-happiness relation uniform across all countries? Previous studies on the marginal effect of income on happiness suggest that happiness may have a nonlinear relationship with GDP, presumably showing saturation after a specific GDP -a concept known as the Easterlin paradox [28]. This paradox tells us that increases in happiness through GDP reach a saturation point, yet what factors improve happiness beyond the saturation point is unknown.
To test the Easterlin paradox, we repeated the analysis over clusters of countries grouped by GDP. Figure 2(a) shows a high correlation between GDP and happiness in the 30 lower-GDP countries (i.e., ρ = 0.40). In contrast, the correlation is no longer evident in the 30 higher-GDP countries (i.e., ρ = -0.04). These results suggest that economic prosperity (as measured by GDP) is crucial for people's happiness but fails to explain happiness in rich countries. The GDP appears to reach a happiness-correlation threshold around the 30th wealthiest country, which corresponds to a GDP of 38,518 dollars. Previous research on the Easterlin paradox has shown that GDP per capita can increase happiness until it reaches a certain threshold but cannot further increase happiness above that threshold. We observe a similar pattern for wealth and happiness across countries. On the other hand, happiness in the 30 wealthiest countries is well explained by urban green space. As shown in Fig. 2(b), urban green space is positively correlated with happiness in the richest countries (i.e., ρ = 0.66, p < 0.01), but this correlation is not significant in the 30 lower-GDP countries (i.e., ρ = 0.19, p = 0.32). Thus, urban green space is a factor that further increases the happiness of a country after its GDP reaches a certain level.
The regressions for each of the 60 countries ranked by GDP in Table 1 confirm the individual effects of urban green space and GDP on happiness. GDP is the only substantial factor explaining happiness in the 30 lower-GDP countries (models 4-6). In contrast, for the 30 higher-GDP countries, happiness is associated only with the UGS (7-9). These findings suggest that GDP is critical for happiness until it reaches a certain GDP threshold (i.e., the Easterlin paradox), after which urban green space better explains happiness.
The correlation between the UGS and happiness also corroborates the effect of the UGS in rich countries. The correlation in Fig. 2(d) decreases as more countries in decreasing order of GDP are added. The correlation is substantial (i.e., ρ is approximately 0.8) among the countries excluding the top 30. Figure 2(c) summarizes the effects of urban green space and GDP, which intersect with each other around the 30th wealthiest country. For the top 30 countries, urban green space has positive coefficients, but the GDP effect is not significant. These relationships are reversed for less affluent countries.
In summary, economic support seems to be associated with happiness until the essential requirements and living standards are met. However, economic growth is not persistently linked to increasing happiness. After a certain level, urban green space appears to be related to happiness via other social factors.

Urban green space for social cohesion
Our findings highlight urban green space as an indicator that might be correlated with social factors explaining happiness beyond achieving economic success. The question then arises: which social factors connect urban green space with happiness? To identify this connection, we first examine the correlation between the UGS and the socioeconomic variables reported in the World Happiness Report: GDP per capita, social support, life expectancy, freedom, generosity, and perceptions of corruption. Of these six variables, only "social support" has a significantly positive correlation (ρ = 0.43, p < 0.01) with the UGS, as shown in Fig. 3(a), implying that social support could mediate the relationship Solid black arrows denote a statistically significant relationship between a pair of variables with the regression coefficient and the p-value (i.e., * * * p < 0.01). The gray dashed arrow represents a non-significant relationship. Note that the coefficients are calculated with z-scores of the variables to compare the effect size directly between urban green space and happiness. This relationship is consistent with several existing studies suggesting that urban green space is a place of social cohesion [13,14]. Additionally, the relationship is robust for other control variables, including health expenditure, unemployment, gender inequality, and education. On the other hand, as indicated by life expectancy, physical health does not display a significant relationship with green space (ρ = 0.19, p = 0.15), contradicting common sense. The regression analysis on happiness with the UGS and six socioeconomic variables also captures the interchangeability of urban green space and social support (see Additional file 1 for details). Therefore, social support is the major component for connecting urban green space and happiness, and we construct a more elaborate model only with social support.
Here, we employ a moderated mediation model [29] to characterize the complicated relationships among urban green space, social support, GDP, and happiness. In moderated mediation models, the moderator describes a variable's conditional effect through the interaction term, and the mediator describes a variable's indirect effect connecting the other two variables. Accordingly, moderated mediation models determine the pathway of directed interactions between multiple variables.
First, we examine the mediation effect of urban green space and social support as independent variables. The mediation regression model shows that social support mediates the relationship between urban green space and happiness such that (1) urban green space improves social support and (2) social support is positively correlated with happiness. The mediation effect is significant (3) only when GDP is considered in the model. Consequently, our moderated mediation model combining these three effects, shown in Fig. 3(b), presents the pathway by which green space is related to happiness through social support, given that GDP moderates the effect of social support. If we describe this relationship in equations, where H, M, S, and G represent happiness, log-GDP, social support, and the NDVI per capita, respectively, and the β values denote the regression models' coefficients (see Additional file 1 for details). Our moderated mediation model can be used to estimate the amount of urban green space that corresponds to a certain amount of increase in happiness according to In Equation (3), the required ratio of urban green space in a country decreases as its log-GDP increases. The required increase in urban green space per capita can be estimated for each country based on its current GDP value. For example, in the United States, an increase in the NDVI per capita of urban green space of 36.1908

Discussion
This paper revealed a global relationship between urban green space and happiness in over 60 countries using high-resolution satellite imagery. Urban green space shows strong correlations with happiness in developed countries (i.e., countries with higher GDPs), which suggests that urban green space is a key to understanding happiness beyond economic success. Our moderated mediation model further elucidates this relationship as social support mediates the green-happiness relationship and GDP moderates social support and happiness. This sophisticated model can estimate the differential amount of additional green space corresponding to the same increment of happiness in each country. The current study newly defined the UGS concept, which can be used to calculate the amount of green space at any spatial scale accounting for population density. We compared several green space measures and proposed using the logarithmic NDVI per capita as a preferred measure of the UGS. This index was validated through experiments and it makes it possible to investigate green space at a global level, allowing us to perform crosssectional research on green space. Furthermore, the method of obtaining the UGS can be utilized to investigate any spatial areas such as blue space (i.e., aquatic environments such as lakes and shores) [30,31].
Our findings have multiple policy-level implications. First, public green space should be made accessible to urban dwellers to enhance social support. In doing so, one critical aspect is public safety. If public safety in urban parks is not guaranteed [32,33], its positive role in social support and happiness may diminish. The meaning of public safety may change; for example, ensuring biological safety will be a priority in keeping urban parks accessible during the COVID-19 pandemic [34]. The high indoor transmission rate of the virus [35] will increase awareness and the importance of open spaces such as urban parks. While some urban parks may be closed during lockdowns, some reports suggest that viewing them from home could also help reduce stress during the pandemic [36]. Second, urban planning for public green space is needed for both developed and developing countries. While our findings confirmed a strong correlation between urban green space and happiness in developed countries, the same positive relation holds for developing countries, albeit to a smaller degree. Furthermore, it is challenging or nearly impossible to secure land for green space after built-up areas are developed in cities. Therefore, urban planning for parks and green recovery (new greening in built-up areas) should be considered in developing economies where new cities and suburban areas are rapidly expanding [37,38].
In addition to the above, recent climate changes can create substantial volatility in sustaining urban green space. Extreme events such as wildfires, floods, droughts, and cold waves could endanger urban forests around the world [39]. On the other hand, global warming could also accelerate tree growth in cities more than in rural areas due to the urban heat island effect [40]. Ultimately, the environmental influence is bidirectional; urban green spaces affect local climates by reducing carbon dioxide levels [41] and providing a cooling effect inside the city that indirectly affects people's well-being. Thus, more attention must be paid to predicting climate changes and discovering their impact on public places since extreme changes could hamper the benefits of urban green space.
As an exciting future direction, satellite images of higher spatiotemporal resolutions can be used to compute urban green space scores. This paper focused on the correlation across countries fixed in time, given the short span of the Sentinel-2 dataset launched in 2015. A causal analysis could be performed with satellite imagery data for a longer span. Additionally, our dataset does not cover all countries in the world. Fortunately, our observations from the 30 lower-income countries anticipate the substantial effect of GDP in other developing countries excluded in our analysis. Although we analyzed the highestresolution public dataset of satellite imagery in this study, our method still has room for application to higher-resolution nonpublic datasets such as the household level (less than 10 m resolution) available in national-scale health datasets [42]. Since satellite imagery cannot account for green space inside buildings (such as green walls), future research could quantify the effect of these mini-scale green spaces using computer vision [43].
Finally, this study is not without limitations. Our country-scale analysis inevitably has a small size of data, while our findings are statistically significant. A dataset with different levels of coverage and data on happiness at the national level and urban green space in cities that contain at least 10% of the population may also include a potential bias due to within-country heterogeneity. Better accessibility to satellite imagery data at different scales and high-resolution global happiness data will enable a better estimation accounting for variations across different cities in a country.

Collecting happiness and remote sensing data
We use the happiness scores from the World Happiness Report [24] and the NDVI scores from Sentinel-2 satellite imagery as remote sensing data to identify the relationship between happiness and green space. The 2018 World Happiness Report covered 156 countries. The report provides an annual survey of how happy citizens perceive themselves to be and ranks the countries by their happiness scores. The score is the average of the participants' responses on how happy they are on a scale from 0 to 10. While many socioeconomic indicators (e.g., unemployment and inequality) may affect happiness, not all of these factors are measured annually across 156 countries. Instead, the report describes happiness with six primary socioeconomic indicators: GDP per capita, social support, life expectancy, the freedom to make life choices, generosity, and perceptions of corruption. For example, the social support variable is based on binary responses (yes/no) to a Gallup World Poll question: "If you were in trouble, do you have relatives or friends you can count on to help you whenever you need them or not?" To quantify urban green space in global cities, we use the Sentinel-2 dataset, which provides the highest spatial resolution (10 m) among the publicly available satellite imagery datasets (e.g., 30 m resolution in Landsat series) [3,10,11,18]. With this high resolution, we can identify granular green space, including street vegetation and home gardens that could not be detected in other public datasets. When using satellite imagery to detect small-scale vegetation, it is critical to consider the season in which the images were obtained [3,10,11,18]. We used images from summer: June to September 2018 for the Northern Hemisphere and December 2017 to February 2018 for the Southern Hemisphere. Satellite images with less than 10% cloud cover were used; when such images could not be obtained for the study period, data from 2019 were used instead.
The NDVI is a well-known remote sensing indicator of green vegetation areas in satellite images [26]. It detects vegetation as the difference between near-infrared and red light in the value range from -1 to +1. High NDVI scores include urban green spaces such as official parks, backyards, street trees, mountains, riverbanks, golf courses, and urban farmlands. There are a few well-known variants of NDVI [23], such as the SAVI [44], which is corrected for soil brightness, and the enhanced vegetation index (EVI) [45], which is corrected for atmospheric effects. All NDVI, SAVI, and EVI2 scores can be calculated from the two spectral bands of Sentinel-2, red (band 4) and near-infrared (NIR, band 8), as follows: The robustness of the results based on the three green space measures was verified using the NDVI as the primary metric.

Measuring the amount of green space
The vegetation indices are measured in three steps, as illustrated in Fig. 4(a). The first step is to identify target cities containing at least 10% of each country's total population and represent the country's overall happiness. The second step is to extract only the built-up areas within the identified cities' administrative boundaries. As the boundaries of cities are historically and culturally constructed and are often arbitrary, the cities' size needs to be standardized; some cities include vast suburban areas (e.g., Istanbul) or natural areas (e.g., deserts in Dubai). Thus, referring to global land cover data from the EU's Copernicus Programme [46], we focus on urban built-up areas to quantify urban green space. Finally, the vegetation indices (NDVI, EVI2, and SAVI) are calculated for the extracted urban areas. The final step is to compute the amount of green space in each country, determined from the measured vegetation indices. Here, we define the amount of green space as the logarithm of the total NDVI of built-up areas in the target cities divided by the cities' total population, called the UGS, as a metric for urban green space. The UGS is calculated as follows: where NDVI(b) is the NDVI value of pixel b within built-up areas b(c) in city c and N c is the population of city c. In this calculation, we adjusted negative NDVI values to zero [23] to prevent errors caused by the accumulation of negative values in areas next to bodies of water (see Additional file 1 for the entire dataset).
Additional file 1. Supporting figures and tables for the main manuscript (PDF 679 kB)