An Exploration of Human Well-Being Bundles as Identifiers of Ecosystem Service Use Patterns

We take a social-ecological systems perspective to investigate the linkages between ecosystem services and human well-being in South Africa. A recent paper identified different types of social-ecological systems in the country, based on distinct bundles of ecosystem service use. These system types were found to represent increasingly weak direct feedbacks between nature and people, from rural "green-loop" communities to urban "red-loop" societies. Here we construct human well-being bundles and explore whether the well-being bundles can be used to identify the same social-ecological system types that were identified using bundles of ecosystem service use. Based on national census data, we found three distinct well-being bundle types that are mainly characterized by differences in income, unemployment and property ownership. The distribution of these well-being bundles approximates the distribution of ecosystem service use bundles to a substantial degree: High levels of income and education generally coincided with areas characterised by low levels of direct ecosystem service use (or red-loop systems), while the majority of low well-being areas coincided with medium and high levels of direct ecosystem service use (or transition and green-loop systems). However, our results indicate that transformations from green-loop to red-loop systems do not always entail an immediate improvement in well-being, which we suggest may be due to a time lag between changes in the different system components. Using human well-being bundles as an indicator of social-ecological dynamics may be useful in other contexts since it is based on socio-economic data commonly collected by governments, and provides important insights into the connections between ecosystem services and human well-being at policy-relevant sub-national scales.

section_content: Human well-being is dependent on ecosystems for provisioning and regulating services like food, clean drinking water, and protection from hazards such as floods, as well as cultural services such as spiritual enrichment and recreation [1, 2]. Global environmental change threatens the supply of essential ecosystem services, and it is imperative to understand how such changes may influence human well-being into the future [3] [4] [5]. Yet despite extensive research in recent years, many gaps in our knowledge of the linkages between ecosystem services and human well-being remain [6] [7] [8]. In part, this is due to a lack of studies that empirically measure both ecosystem services and human well-being, and that consider the socialecological systems within which these interactions take place [7, 9]. 

An understanding of the type of social-ecologicalsystem that underpins the interactions between ecosystems and society can yield important insights into resource use patterns and dynamics. For example, in a recent review Cumming et al. [10] identified different types of social-ecologicalsystems based on the strength of the connection between people and their local environment. In rural "green-loop" systems, where people depend on local ecosystems and the services they provide, feedbacks between ecosystems and societies are direct and clear. However, as populations grow, production roles become more specialized, market economies develop, and fewer people come into direct contact with local ecosystems. These societies transition to a "red-loop" system, in which resources are no longer extracted mainly from local ecosystems, but instead obtained from far-away places and impacts are outsourced. The result is a weakening of feedbacks between ecosystem services and society, so that people's reliance on ecosystem services becomes obscured. 

Applying this system typology, Hamann et al. [11] mapped social-ecological systems in South Africa, based on bundles of ecosystem service use. The bundles consisted of six locallysourced provisioning services that cover people's basic needs: animal production, crop production, freshwater from a natural source (river or spring), natural building materials, wood for cooking, and wood for heating. The level of direct use of these six services among households was determined from census data for each municipality, and subjected to a cluster analysis to identify three distinct types of ecosystem service use bundles across South Africa. The three bundle types represented an overall high, medium and low level of direct ecosystem service use among households, which was taken to correspond, broadly, to "green-loop", "transition", and "red-loop" systems, as defined by Cumming et al. [10]. 

Beyond differences in the level of use of or dependence on ecosystem services, social-ecological systems are characterized by other fundamental differences in structure and dynamics. Cumming et al. 's [10] system typology indicates that the green-to red-loop transition brings with it major shifts in the way society operates: As populations grow and become more urbanized and disconnected from their local ecosystems (while becoming increasingly dependent on natural resources sourced from far-away ecosystems), economies scale up, the work force becomes specialized and wide-spread infrastructure development takes place. Hospitals and universities are built, and technology advances. In these red-loop dynamics, traditional wellbeing measures like income and life expectancy tend to improve. In contrast, green-loop dynamics keep the community closely connected to their local ecosystems, but a lack of specialization and technological advancement means that these communities often do not engage in broader market-based economies, so that well-being measures such as income and health may not improve. We therefore expect that different social-ecologicalsystems will not only express different ecosystem service use patterns, but also different human well-being patterns. 

This paper aims to test this proposition by comparing the distribution of human well-being bundles and ecosystem service use bundles in South Africa. This approach is significant on two levels: Firstly, taking a bundles approach may help to disentangle some of the complexities inherent in the relationship between ecosystem services and human well-being. Using ecosystem service and human well-being bundles to identify different kinds of social-ecologicalsystems in a landscape may tell us more about the area's ecosystem services, human well-being, and households, as well an annual household income, education and unemployment levels, and property ownership percentages can be found for South Africa's municipalities. The direct web address for the SuperWEB2 database is: [http://interactive2. statssa. gov. za/webapi/jsf/login. xhtml]. Login is not required -all data is accessible under the "Guest" login. The other database is called 'Nesstar'. Here, users can find agricultural data and mortality data. These data are a 10% sample of the census, but can be weighted to be representative of the total census at the municipal level. All data is downloadable, and again, no login is required ("guests" login is sufficient to access the datainstructions are clearly provided on the site). The direct web address for Nesstar is: [http:// interactive. statssa. gov. za:8282/webview/]. Finally, StatsSA can be contacted directly for any data that users are not managing to download from the online databases (for technical reasons, for example). 

linkages between the two than an analysis focusing on individual indicators. Secondly, data on a variety of human well-being and development indicators are often more widely available than data on ecosystem services and their use. If different social-ecologicalsystems and their characteristic resource use dynamics can be identified using human well-being indicators, instead of ecosystem services, then it could assist with ecosystem management in data poor regions. 

Our primary objective is therefore to ascertain to what extent distinct human well-being bundles correlate with ecosystem service use bundles, and whether human well-being bundles can be used as a proxy to identify the same social-ecologicalsystems that were identified with the ecosystem service use bundles defined by Hamann et al. [11]. To achieve this objective we construct a multidimensional well-being bundle based on publicly available indicators that together represent the Millennium Ecosystem Assessment's (MA) [1] five constituents of human well-being: security (e. g. personal safety, security from disasters); basic material for a good life (e. g. sufficient nutritious food, shelter); health (e. g. strength, access to clean air and water); good social relations (e. g. social cohesion, mutual respect, ability to help others); and freedom of choice and action. While human well-being is a composite measure that has been defined in many different ways [12, 13], we chose an indicator bundle that captures the MA's well-being constituents because these constituents cover a broad and inclusive definition of human well-being and can be represented by indicators that are commonly collected by governments. In addition, even though their conceptual connection to ecosystem services was outlined in the MA, few empirical studies have explored these links at above-local scales. We also investigate the relationship between the ecosystem service and human well-being bundles, as well as their individual indicators, and find significant connections at a sub-national scale. In addition, we evaluate the bundles approach by analysing similarities and differences between the human well-being bundles we identified and more traditional well-being measures like economic activity and a South African deprivation index.
section_name: Introduction

section_content: We mapped human well-being in South Africa using a bundle of well-being indicators and compared the distribution of well-being bundles to that of ecosystem service use bundles, as well as other widely-used well-being indicators. Our unit of analysis was the municipality, which is the smallest administrative unit controlled by a local government council in South Africa. There are a total of 234 municipalities in the country, with an average size of 5217 km 2 (range: 252-36 128 km 2 ) and an average of number of 61 753 households (range: 1784-1 434 856).
section_name: Materials and Methods

section_content: South Africa is a market-based economy, the 2 nd largest in Africa after Nigeria, and the dominance of cash-based exchanges reaches into the most rural areas, supported by a large-scale social grant system [14]. Aspirations for modern lifestyles and economic opportunities are evidenced by the high migration rates to urban areas, especially among the young [15, 16]. This context played an important role in selecting indicators that reflect current South African realities in terms of human well-being. 

The five constituents of human well-being identified by the MA [1]-i. e. basic material for a good life, health, security, good social relations, and freedom of choice and action-were used as a guide in the choice of indicators for a human well-being bundle (Table 1 ). All data were obtained from the most recent South African population census conducted in 2011 [17], and can be freely accessed via www. statssa. gov. za. The indicators were selected based on their ability to represent one or more of the MA's well-being constituents. For example, average household income was selected as an indicator of having a livelihood and being able to afford the materials to cover one's basic needs (such as food, shelter and clothes). However, income may also contribute to health by providing access to healthcare. Furthermore, having an income provides opportunities that support freedom of choice and action. Average age of death is another indicator of health, and was chosen because life expectancy at birth, a more common health indicator, is not available for South Africa at the municipal level. Property ownership can be considered an indicator of security, since owning property indicates secure tenure and assets. In South Africa, property ownership is not necessarily tied to income, since traditional tenure systems and government housing schemes allow property ownership even when households are poor. Employment is a well-being indicator for a number of different well-being constituents, including basic material for a good life, and security. In addition, the percentage of people who are unemployed or discouraged work-seekers can be considered an indicator for the status of social relations in the municipality, especially in the context of a predominantly cash-based economy like South Africa. The MA characterized "good social relations" by social cohesion, mutual respect and the ability to help others. Being unemployed or even discouraged in one's efforts to find a job negatively impacts social relations, since having a job not only contributes to income but also conveys status in the community, elicits respect by others, and fosters self-esteem [18, 19]. At the aggregate level (e. g. municipalities or nations), unemployment may lead to increased inequality and social tension, and has been linked to increased crime rates [20] [21] [22]. While all the indicators described so far also contribute to freedom of choice and action, education is an especially important indicator for this well-being constituent, as a good education leads to employment opportunities and allows people to realize their life goals and participate more fully in society. Higher education levels within a population are also strongly correlated with civic engagement and democracy [23], which enhance personal freedom of choice and action. A correlation analysis of the indicators making up the well-being bundle was performed and visualized using the corrplot function in R [24] (S1 Fig). 

Our indicator selection is subject to data availability at a national scale and we acknowledge that our choice represents one of many possible human well-being bundles, but we feel that the bundle chosen here reasonably captures the well-being dimensions identified by the MA relevant to the South African context. To assess the sensitivity of our results to indicator selection, we repeated the analysis with a different "basic material for a good life" indicator, namely a poverty line measure. Average income may mask a skewed distribution of income within a municipality. This alternative indicator represented the level of poverty, measured in terms of the percentage of households per municipality that lived above a certain poverty threshold. This threshold was defined as an annual household income of more than ZAR 76 400, based on national poverty line data and census data income bands [17, 25].
section_name: Selection of human well-being indicators for a well-being bundle

section_content: : Basic material for a good life; health; security; good social relations; freedom of choice and action. All data were sourced from the South African National Census 2011 [17].
section_name: Table 1. Selected indicators for a human well-being bundle capturing the five well-being constituents of the Millennium Ecosystem Assessment (MA) [1]

section_name: Indicator

section_content: After calculating the values for each human well-being indicator per municipality, the data were visualized and analysed in ArcGIS 10. 0 [26]. The spatial clustering of all services was determined using spatial autocorrelation (Global Moran's I statistic [27] ). To identify groups of municipalities that shared similar combinations of values for the five well-being indicators, the values for all indicators were first scaled between 0-1, and then a k-means cluster analysis was performed using the Hartigan-Wong algorithm [28] with 25 random starts and a maximum of 10000 iterations to find the cluster solution with the lowest within-cluster sum of squares. The analysis was performed in R using the kmeans function from the stats package [29]. K-means clustering (number of clusters = 3) was chosen to repeat the clustering procedure used by Hamann et al. [11]. Results of the cluster analysis are visualised in a multidimensional scaling ordination diagram based on the dissimilarities in the human well-being measures between all the municipalities [30]. Multidimensional scaling was performed and plotted in R using the cmdscale and plot functions from the stats and graphics packages [29]. Since the clustering method can have an effect on which clusters are found in the data, hierarchical clustering was also performed to compare results between the different clustering algorithms. In addition, cluster validation procedures using the clValid package in R [31] indicated support for two clusters within the data, which was likewise investigated for comparative purposes. Human well-being bundles were visualized using star plots in R and clusters were mapped in ArcGIS.
section_name: Mapping human well-being bundles

section_content: The spatial distribution of human well-being bundles was compared to that of ecosystem service use bundles as determined by Hamann et al. [11]. The overlap between the different human well-being and ecosystem service bundle types was calculated in terms of area and population. To analyse relationships between individual indicators and bundle types, a multinomial logit model using the mlogit package in R [32] was run on the human well-being bundle types as the dependent variable and individual ecosystem service use indicators as predictor variables, as well as vice versa for ecosystem service bundles and individual well-being indicators (standardized to z-scores).
section_name: Exploring the relationship between human well-being and ecosystem service use

section_content: To compare our human well-being bundle to more traditional measures of well-being like GDP and the Human Development Index (HDI), we calculated the Gross Value Added (GVA) for each municipality from the South African Geospatial Analysis Platform [33]. GDP figures are not publicly available at the municipal level, but GVA is related to GDP and is a measure of regional economic activity. As an example of a well-being index, we used the South African Index of Multiple Deprivation (SAIMD) which is a publicly available index of well-being at the municipal level that was derived for the national Department of Social Development from the nation-wide 2007 community survey [34, 35]. It is composed of a total of eleven indicators that reflect four dimensions of deprivation: income and material deprivation, employment deprivation, education deprivation, and living environment deprivation (the latter referring to poor living conditions). Even though the data used to compile the GVA figures and SAIMD ranks precede the census data by 2 to 4 years, we assume the patterns to not have changed significantly between data compilations. Differences in per capita GVA scores and SAIMD ranks between municipalities in different human well-being bundle types were visualized in R using boxplots. Data were tested for normality using the Shapiro-Wilk test [36]. Since the per capita GVA data were of unequal sample size, found to contain outliers and non-normally distributed, they were log-transformed before being subjected to Welch's t-tests on unique pairs [37]. Differences in deprivation ranks between pairs of human well-being bundle types were assessed for significance using Mann-Whitney-Wilcoxon tests [38].
section_name: Comparing human well-being bundles and other well-being indicators

section_content: Mapping the five human well-being indicators showed that elevated levels of well-being were mostly found in the western and north-eastern parts of the country (Fig 1 ). The exception was property ownership, where high percentages of households that own their dwelling were found along the south-east coast where incomes are typically low.
section_name: Results

section_content: A cluster analysis of the municipalities based on their values for the five human well-being indicators resulted in three distinct types of well-being bundles (Fig 2 ). The first bundle type exhibited high levels of income and education, but unemployment was high and property ownership was the lowest of all bundle types. The second bundle was characterised by medium income, the highest age of death and lowest unemployment rate. The third bundle type exhibited overall low levels of well-being compared to the other bundle types, with the highest unemployment, and the lowest household income, age of death, and education levels. However, the percentage of households in a municipality that own their dwelling and have paid it off fully was higher in this bundle than in the other two-due to historical and cultural tenure arrangements in South Africa (see discussion). In order to facilitate the discussion of results, these bundles are henceforth referred to as the "high income", "medium income" and "low income" bundle, respectively. We use income to name the bundles, as it is one of the indicators that varied most clearly between the different bundle types, but emphasize that the bundles are multidimensional and represent far more than differences in income. 

When the different well-being bundle types were mapped across South Africa (Fig 2 ), the resulting pattern of human well-being was significantly clustered in the landscape (Moran's I, p < 0. 01). In total, 60 of the 234 municipalities fell into the high income bundle, 63 municipalities fell into the medium income bundle, and 111 municipalities fell into the low income bundle. That corresponded to 15. 3%, 46. 2% and 38. 5% of the total land surface area, and 60. 7%, 7. 0% and 32. 3% of the total population, respectively. The high income bundle occurs mainly in densely populated regions around metropolitan centres, while the medium income bundle occurs in the sparsely populated and rural central and western regions of the country. The low income bundle occurs mostly in more densely populated rural and semi-rural areas in the east of the country. 

A classical multidimensional scaling analysis of the municipalities and their human wellbeing measures revealed a clear delineation of the three bundle types identified in the cluster analysis (S2 Fig). When the clustering algorithm is restricted to two clusters, the medium income bundle area splits roughly equally into either high or low income bundle areas (S3 Fig). When hierarchical clustering is applied, the three-cluster pattern looks similar to that obtained after k-means clustering (S4 Fig), though more areas are classified within the low income bundle. If income is replaced with a measure of poverty to test the sensitivity of the analysis to indicator selection, 11% of the municipalities were assigned to a different bundle type than in the original analysis.
section_name: Human well-being bundles

section_content: We found that the distribution pattern of human well-being bundles mirrored, to a substantial extent, the pattern in ecosystem service use bundles as mapped by Hamann et al. [11] (Fig 3) : 93% of high and medium direct ecosystem service use areas coincided with the low income bundle areas. 

The relative area comparison between the two maps shows that there was no overlap between high or medium ecosystem service use areas and high income bundle areas (Fig 4 ). Only a very small percentage (3. 5%) of the area characterized by the medium income bundle also exhibited high or medium levels of ecosystem service use. In low income bundle areas, low and medium levels of ecosystem service use were found in nearly equal proportions (41. 4% and 43. 5%, respectively), while 15. 2% of the area exhibited high levels of ecosystem service use. 

When considering the percentage of total land area and population in South Africa (Table 2 )-as opposed to the relative area of different ecosystem service use and well-being 1 for an explanation of the units), calculated across all municipalities that were clustered together within a bundle type. Petal length represents the average value (relative to the absolute maximum) of each indicator among the municipalities that share a similar bundle type. The map shows the distribution of the three different bundle types in South Africa. Areas of the same colour share similar wellbeing bundles at the municipal level and are significantly clustered in space (Moran's I, p < 0. 01). Provincial borders and major metropolitan centres are shown. 2 ).
section_name: Comparison of human well-being and ecosystem service bundles

section_content: In terms of determinants of bundle membership, we found that an increase in average household income and age of death significantly increases the log-odds of a municipality being classified as part of the medium or low ecosystem service use bundle (Table 3 ). The other well-being indicators do not play a significant role in determining the bundle of ecosystem service use in South African municipalities. The results also show that the human well-being bundle type is sensitive to changes in the level of animal and crop production. Similarly, an increase of the percentage of households that make use of natural sources like rivers or springs for freshwater or that use wood for heating significantly decreases the log-odds of a municipality being classified as part of the high income and education well-being bundle. Human well-being bundles that are referred to as "high", "med", and "low" correspond to the "high income and education", "medium income and low unemployment" and "low income and high property ownership" bundles, respectively. Significance at p < 0. 05 denoted by * and at p < 0. 01 by **.
section_name: Relationships between single indicators and bundle types

section_content: This study used a bundles approach to map human well-being at a sub-national scale, and investigates whether different human well-being bundles and ecosystem service use bundles show similar distributions across South Africa. We found three distinct human well-being bundles which reflected differing levels of multidimensional well-being across the country. These human well-being bundles were found to mirror the distribution of ecosystem service use bundles to a substantial degree, suggesting that well-being bundles could be used as proxies for ecosystem service use and underlying social-ecologicaldynamics in data poor areas. Our results show that well-being across South Africa is differentiated and nuanced, and may not be adequately reflected by single well-being indicators or composite indices that mask diversity between well-being constituents.
section_name: Discussion

section_content: The main distinction between the medium or high income bundle and low income bundle was property ownership. The high levels of property ownership in the otherwise low well-being bundle may be explained by the high proportion (26. 4%) of land under the rule of traditional authorities that makes up the low income bundle areas in South Africa. In these communal lands a property may be acquired from a local traditional leader or village chairman in charge of land allocations, usually for a nominal fee [39], which means that people with otherwise few assets have access to property without the usual risks incurred by lower-income homeowners in other countries (like low rates of appreciation and subprime mortgages) [40, 41]. Globally, land titles are argued to convey social and economic empowerment to their owners, particularly if they are poor, as the titles represent assets that can be used as collateral security and thereby promote upward mobility [42, 43]. However, it is debated whether this holds true in the global South, and especially in South Africa [44] [45] [46] [47]. The title deeds held in communal areas do not necessarily promote lending to the poor, since banks are generally not willing to expose themselves to the high risk of non-repayment on land parcels with relatively low market value. However, while these homes may not lead to upliftment out of poverty in the form of income or credit, they nevertheless provide security of shelter (especially in old age) and represent a tie to one's community that may contribute to well-being [48] [49] [50]. 

Our results support the fact that it is the poorer and often most vulnerable sectors of society who depend most on their immediate natural environment to meet their basic needs [1, 51, 52], but not in all cases. We find that high and medium ecosystem service use areas are found almost exclusively in low income bundle areas (Fig 3b ), which constitute 24. 1% of the population (Table 2 ), or 12. 5 million people. Conversely, all of the high income bundle areas, as well as the vast majority of the medium income bundle areas, overlapped with the low ecosystem service use areas (Fig 4 ). However, a large proportion of the low income bundle area also overlaps with low levels of ecosystem service use. This most likely reflects the fact that in South Africa, many people may cover their basic needs with the help of social grants, which are distributed to almost 30% of the population [14, 53]. Low levels of human well-being across the majority of indicators chosen here therefore do not necessarily imply high levels of ecosystem service use or dependence, but high levels of ecosystem service use among households almost exclusively implies relatively low overall well-being.
section_name: Well-being bundles and overlap with ecosystem service use bundles

section_content: Hamann et al. [11] used bundles of ecosystem service use to identify green-loop, transition and red-loop systems in South Africa. These systems are theorized to represent distinct areas characterized not only by differences in the level of ecosystem service use among households, but also by differences in the broader socio-economic dynamics that affect human well-being [10]. We investigated whether well-being bundles would be able to act as proxies for ecosystem service use bundles in identifying these different system types. As outlined above, the distribution of human well-being bundles was only partly concordant with the distribution of low, medium and high ecosystem service use bundles-representing red-loop, transition and green-loop systems, respectively. While generally high levels of well-being (high and medium income bundles) are able to identify 79% of red-loop areas, low levels of well-being (low income bundles) were not able to distinguish between green-loop, transition or red-loop systems. 

This lack of congruence in low income bundle areas could be due to multiple factors, including the particular well-being indicators used in this study. We also hypothesize that this discrepancy between ecosystem service use and low human well-being may be due to a difference in slow and fast-changing variables [54]. The level of direct use of locally-sourced ecosystem services may change within very short timeframes at the household level. When electrification, sanitation and infrastructural development of an area takes place, households may switch within months from using fire wood for cooking and fetching water from a stream to using electric ovens and the municipal water supply, i. e. the use of locally available natural resources diminishes rapidly. However, the indicators of well-being that were included in this study (especially education, average life span and unemployment) are slower to improve, which may explain the areas in which both human well-being and ecosystem service use are relatively low. Our results therefore suggest that transformations from the green-loop to the red-loop system do not always imply an immediate improvement in human well-being. Current well-being may reflect ecosystem service use levels from ten years ago, not present levels. It has been proposed that similar time lags contribute to the reverse situation, i. e. when use of ecosystem services is high and the environment is being degraded, yet human well-being does not diminish (termed the "environmentalist's paradox") [8]. 

Differences between slow and fast-changing variables, and specifically the temporal dynamics of human well-being compared to ecosystem service use patterns, raises important questions about attempts to investigate the linkages between ecosystem services and human well-being using data from a single point in time. Maps, especially, portray a snapshot in time, often with the implicit assumption that processes captured by the image are in some kind of equilibrium. If, as our results suggest, there is a delay in the human well-being variables when systems shift from green-loop to red-loop dynamics, then it is more appropriate to observe and measure ecosystem service use and human well-being patterns over an extended period of time. This has implications for research design in the field of ecosystem service science, and the importance of longitudinal or repeated cross-sectional data in untangling the links between ecosystem services and human well-being.
section_name: Identifying green-loop and red-loop systems with human well-being bundles

section_content: An increase in household income appears to play the dominant role in moving from high to low levels of ecosystem service use (Table 3 ). The only other significant predictor was average age of death. In our case, the overall level of ecosystem service use in a municipality therefore appears to relate most closely to well-being indicators representing the basic material for a good life, and health. 

At the same time, as income and overall well-being increases, animal and crop production at the household level decreases (Table 3 ). As households transition from lower to higher income levels this may compromise food security and resilience to socio-economic upheaval. In South Africa-like in most developing countries-areas where well-being is generally low are mainly rural communities in which some proportion of household income is made up of remittances from household members working in urban centres as labour migrants [55] [56] [57]. Increases in household income through remittances may be relatively small compared to incomes in wealthy urban centres, but they may be large enough for rural households to transition into low ecosystem service use dynamics, abandoning local small-scale agriculture in the process, particularly if the remaining rural household members are elderly, sickly, or too young to farm [58, 59]. If, however, the remittances cease due to unforeseen circumstances (e. g. the wage earner's death, illness, or job loss) the household is left with a much reduced income and no safety net in the form of small-scale subsistence farming [60].
section_name: Predictors of human well-being and ecosystem service use bundles

section_content: Using bundle types to map human well-being allows a comparison of different constituents of well-being across an area. Such an analysis of differences between well-being components is often lacking in well-being analyses based on widely-used one-dimensional indicators like income or composite indices like the Human Development Index (HDI) [61], which tend to reduce the complexity and mask nuances [62] [63] [64]. In our case, average per capita GVA and South Africa's deprivation index corresponded to the income component of our bundles (Fig 5 ). However, both the per capita GVA and deprivation index miss the covariation of wellbeing constituents across the landscape as identified by the bundles approach. For example, the medium income bundle areas exhibit lower levels of economic activity and have a lower annual household income than the high income bundle areas, but unemployment is lower and life spans are longer in these mostly rural areas. Overall well-being may therefore arguably be higher in the medium income bundle than in the high income bundle, even though the economic data would suggest otherwise. On the other hand, single indicators have the advantage of being transparent and easily communicated. Choosing well-being measures will therefore always be a subjective exercise and must be tailored to requirements [63].
section_name: Comparison of human well-being bundles and other well-being indicators

section_content: The bundles approach applied here contributes to a richer understanding of the pattern of multidimensional human well-being in South Africa, and the linkages between human well-being and ecosystem services. However, it is challenging to construct a well-being bundle that is made up of indicators which do not bias the analysis towards modern aspirations and definitions of human well-being. It could be argued that indicators like income, education and unemployment may be important measures of well-being in wealthy modern, often urban societies (red loop systems), but are less important at determining levels of well-being in green-loop systems, i. e. rural, close-knit communities where not all interactions are cash-based and other aspects of well-being may be more highly valued. As our results show, simply identifying areas in which traditional red-loop well-being indicators exhibit low values is not necessarily sufficient to detect green-loop systems. Our approach is reliant on existing, publicly available data that are not tailor-made to measure all aspects of human well-being, which may partly explain why our wellbeing bundle was not entirely effective in identifying green-loop systems. A modified human well-being bundle that includes indicators more suited to represent green-loop communities (indicators like community cohesion and identity, resource tenure, or access to nature, for example) might help differentiate the green-loop/red-loop social-ecologicalsystems more successfully than the bundle chosen here. However, the challenge of finding data at relevant (sub-national) scales for these more unorthodox well-being measures remains difficult to overcome [65]. 

Furthermore, we examined human well-being through objective measures only, without taking into account subjective well-being or self-reported happiness. Research on subjective well-being is extensive and growing [66] [67] [68] [69], but due to data limitations at the scale of our analysis we chose to focus on objective well-being in South Africa. However, some of the objective indicators chosen for this study (level of income, education and unemployment) have been shown to be highly correlated with subjective well-being, both in developing countries generally [70, 71], and within South Africa more specifically [72, 73]. The low income bundle areas in this study are therefore likely to be characterised not only by an overall low level of objective well-being, but also by low levels of subjective well-being. It is more difficult to deduce the level of subjective well-being in the other two bundle types, since they each present a mix of factors that have been shown to contribute to subjective well-being, i. e. the high income bundle has high levels of income and education but also high levels of unemployment, while the medium income bundle has low unemployment but also low education levels. In further research it would be pertinent to include subjective well-being indicators when considering the links between human well-being and ecosystem services, as there is some evidence for a positive relationship between ecosystem services and life satisfaction at the country level [74]. A case study at local scale where both objective and subjective measures of well-being are considered with respect to ecosystem service use would greatly further this line of research. 

In general, we acknowledge that the scale of our analysis (at the municipal level) excludes a finer analysis of more localized and disaggregated interactions between well-being and ecosystem services, which are the focus of much of the research done in this field (e. g. [75] [76] [77] [78] ). Especially in South Africa, where inequality is high, variability in the human well-being indicators within each of the municipalities is likely to be high as well. More fine-scale data at the sub-municipal level would help to identify particularly unequal and diverse municipalities in terms of well-being and resource use dynamics, but is unavailable at present for all indicators assessed here. Further research is also required to investigate what drives the differentiation into green-loop, transition and red-loop systems within low income bundle areas, and could include an analysis of variables such as population density, land tenure, cultural identity, or land degradation.
section_name: Limitations and further research needs

section_content: The human well-being bundles identified in this study represent a novel approach to mapping well-being, and illustrate that well-being across South Africa is differentiated and nuanced in a way that may be missed by single well-being indicators or composite indices that mask diversity between well-being constituents. Furthermore, this approach allowed a spatially explicit comparison of human well-being and ecosystem service use patterns at the sub-national scale, and we found that human well-being bundles could partly identify the green-loop, transition and red-loop systems that were previously identified using ecosystem service use bundles: Almost all medium and high income bundle areas coincided with low ecosystem service use (red-loop systems), and the majority of low income bundle areas coincided with medium and high ecosystem service use (transition and green-loop systems). We hypothesize the remaining discrepancies in the well-being and resource use patterns may be due to time lags between fast-changing ecosystem service use dynamics and slower-changing well-being variables. Such asynchronous dynamics are to be expected in social-ecologicalsystems [79], but are not yet sufficiently acknowledged and accounted for in ecosystem service studies that aim to establish a connection between ecosystem services and human well-being. We suggest that this approach of using human well-being bundles to explore green-loop, transition and red-loop systems and their ecosystem service use can be usefully translated to other parts of the world, since it is based on existing socio-economic data that most countries collect regularly and at sub-national levels. Countries may use this approach to assess potential links between human well-being and ecosystem service use, initially to identify areas of concern (e. g. where wellbeing is low and land degradation is high), and then to focus more detailed investigations in those areas.
section_name: Conclusions

section_name: Acknowledgments

section_content: We thank Statistics South Africa for their assistance, and are grateful to Gretchen Daily, Stuart Kininmonth, Kristoffer Hylander, Garry Peterson and two anonymous reviewers for their constructive comments. This work was supported by a grant from the Swedish International Development Cooperation Agency (Sida) and the Swedish Research Council (VR). RB is supported by a Branco Weiss Society in Science Fellowship and the DST/NRF South African Research Chairs Initiative (Grant No 93097 ). BR received support from the WISER project
section_name: Acknowledgments

section_content: This work was supported by a grant from the Swedish International Development Cooperation Agency (Sida) and the Swedish Research Council (VR). RB is supported by a Branco Weiss Society in Science Fellowship and the DST/NRF South African Research Chairs Initiative (Grant No 93097 ). BR received support from the WISER project ( NE/L001322/1 ) funded by the Ecosystem Services for Poverty Alleviation program (ESPA). ESPA receives its funding from the Department for International Development (DFID), the Economic and Social Research Council (ESRC) and the Natural Environment Research Council (NERC). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

section_content: The census data are third party data. They are not owned by the authors, but are owned and collected by StatsSA, the Statistics Agency of South Africa. StatsSA makes all data publicly accessible and there are no restrictions on the data. The authors did not receive any special access privileges. There are a number of different online databases that other people can use to access the data. These databases are maintained by StatsSA, and links to the databases are found on their website: [http://www. statssa. gov. za/]. One of these databases is 'SuperWEB2', where data on the natural resource use of

section_name: Supporting Information