An Investigation of Links between Environmentally Responsible Behaviors and Built and Natural Features of Landscape in Central New Jersey

: Responses from Central New Jersey residents (n = 290) to nine survey items about pro-environmental behaviors were analyzed for their connection to survey respondents’ local built and natural landscape characteristics. These nine survey items were of interest due to their clustering in three groups in earlier dimension reduction analysis. Nine logistic regression models (one for each of these items) were built using a suite of built and natural features of landscape at the municipal level—including land use, population density, and access to parks. The logistic regression models using the landscape variables failed to provide effective explanations of engagement in pro-environmental behavior. In total, only ten landscape feature variables—out of a possible thirty-eight -were used in any of the most effective models. The logistic regression model of proper appliance disposal performed best and could correctly classify responses of whether respondents had disposed of an appliance correctly 84% of the time. The rest of the models were of little use. This suggests that there may not be a connection between these built and natural landscape features and these behaviors, or that the study did not look at location at a level that was granular enough to detect any patterns.


Introduction
There are a wide range of threats facing humanity that endanger our continued ability for communities to thrive-most of which can be conclusively linked to anthropogenic climate change [1,2]. In the face of this, it is necessary to come up with effective solutions for lowering consumption of resources, and increasing conservation and preservation behaviors. There are many possible ways of framing these types of behavior-as pro-environmental behaviors [3][4][5][6], environmentally responsible behaviors [7,8], and environmentally significant behaviors [9,10]. While there are some distinctions among these, fundamentally, they are talking about some set of actions and behaviors undertaken by humans that benefit some conception of the environment in some way (or, at the very least, reduce harm to the environment); some of these definitions include that behaviors be explicitly intentional in this regard, others are more outcome focused, regardless of intentionality. This research was undertaken with the express purpose of measuring the extent to which nine of these behaviors could be linked to built and natural features of landscape. Knowledge of links between features of landscape and behaviors are useful in generating tools that policy-makers and managers can use to manage landscapes to encourage specific behaviors and to achieve certain behaviors as outcomes. respondents as engaging or not engaging in a behavior at a high rate would be desirable, and one that only provided correct classification at a lower rate would not be particularly useful. In this research, we seek to build a series of logistics binary regression models for nine ERBs that would allow such a tool to be built.

Questionnaire
In the summer of 2016, 290 respondents were surveyed on a range of items including demographic information, use of outdoor recreation spaces, and environmentally responsible decision-making behavior. We purchased the sample from Qualtrics, who collected these data via the internet (on computers and mobile devices). Qualtrics is a U.S. based company that bills themselves on their website as performing "experience management". Among other things, they conduct market research. The authors' institution has a subscription to Qualtrics, and at the time of this study, both authors were very familiar with use of this platform. Additionally, other similar services were contacted for price quotes and Qualtrics was the platform where the highest number of responses could be gathered for the amount of funds budgeted for this project.
Criteria for inclusion were residence within the study area and being an adult aged eighteen or older. Respondents were also required to respond affirmatively to a consent form based on IRB anonymous consent protocol. All research was carried out in accordance with Rutgers IRB [#E15-678]. Respondents were identified for their location at the level of the municipality. This survey instrument can be found in its entirety in the Supplementary, as can descriptive statistics for survey items.
The final instrument used was developed through several rounds of piloting and testing based on previous work, the scope and aims of this study, and input from pilot respondents. Questions on demographic information were based on the U.S. Census [42]. These items included age, race, ethnicity, income, household makeup, and pets. The items chosen for inclusion on Environmentally Responsible Behaviors (ERBs), were based on six distinct categories: Direct Actions, Community and Civic Behaviors, Political Engagement, Consumer-consumable goods, Consumer non-consumable goods, and Identity-surveying respondents on thirty-eight pro-environmental behavioral items. These items were based on the background discussed above [6,9,15,43,52].
Earlier work with this dataset uncovered three latent sets of ERBs within the survey items. We used generalized canonical correlation analysis to identify three distinctive sets of behavior, each made of three individual behaviors. Behavioral items were grouped into sets together when respondents answered similarly across the behavioral items within a group; they are not behaviors that are similar in terms of the type of behavior or function that they entail. These groups of behaviors could be graphically represented as occupying similar positions in our dimensionality reduction analysis, indicating that respondents answered survey items similarly. Furthermore, these were all highly weighted within that dimensionality reduction indicating that these were areas wherein there were differences among responses that followed some set pattern. Some responses were dropped from analysis due to non-responses. Therefore, our research postulated that there were underlying factors not readily apparent that could be used to explain each group of behaviors. A more detailed description of this work is given in the Supplementary section.
A graphical representation of these clusters is given below in Figure 1. The term CCV refers to the Canonical Correlation Variables. The locations within the space indicate how similar or dissimilar responses were to the survey items. Two very close CCVs would indicate a high level of similarity across two responses. The further a CCV is from the origin, the greater the amount of difference represented throughout the study population. All three clusters included behaviors that were closer to each other on the first canonical correlation variables than to any other behaviors, and had a total distance of fewer than 0.1 units. A list of all behaviors surveyed, and more detailed results of this analysis can be found in the Supplementary section. The specific wording of these items is given below in Table 1 These behaviors were expected to fall into six a priori groups based on previous work with this population, and similar studies within the field. [15,16,43,44]. These a priori groups were direct actions, community and civic engagement, political engagement, identity, and consumer behavior. Consumer behavior was split into consumer behavior surrounding consumable products such as food and non-consumable products such as appliances. For example, a person who engaged in proenvironmental consumer behavior in their purchasing of consumable products would be expected to report caring about the origin and production methods of their groceries as well as those of hygiene products; a person who engaged in pro-environmental identity behaviors would report that they think of themselves as an environmentalist and also thought about the environmental impacts of their actions. We then hypothesized that these behaviors were linked to differences in land cover and other features of the landscape and built environment, and that respondents who undertook the identified sets of activities lived in physical environments that were alike in a quantifiable way. Table 1. ERBs used in analysis.

Cluster 1
Number of different public officials a respondent indicated they had contacted with a concern in the last year I consider myself to be an environmentalist. [5 levels, "Strongly Agree" to "Strongly Disagree"] Number of environmental organizations an individual indicated belonging to. These behaviors were expected to fall into six a priori groups based on previous work with this population, and similar studies within the field. [15,16,43,44]. These a priori groups were direct actions, community and civic engagement, political engagement, identity, and consumer behavior. Consumer behavior was split into consumer behavior surrounding consumable products such as food and non-consumable products such as appliances. For example, a person who engaged in pro-environmental consumer behavior in their purchasing of consumable products would be expected to report caring about the origin and production methods of their groceries as well as those of hygiene products; a person who engaged in pro-environmental identity behaviors would report that they think of themselves as an environmentalist and also thought about the environmental impacts of their actions. We then hypothesized that these behaviors were linked to differences in land cover and other features of the landscape and built environment, and that respondents who undertook the identified sets of activities lived in physical environments that were alike in a quantifiable way.

Cluster 1
Number of different public officials a respondent indicated they had contacted with a concern in the last year I consider myself to be an environmentalist. [5 levels, "Strongly Agree" to "Strongly Disagree"] Number of environmental organizations an individual indicated belonging to. Think about the last appliance you bought-TV, microwave, refrigerator, etc. When you bought this appliance, how important was energy efficiency to you in considering your options? [5 levels, "Extremely Important" to "Not important at all"]

Cluster 3
Would you be willing to pay more for your food if you knew it was grown sustainably? [Options: No, up to 10%, 25%, 50%, 75%, or double] Think about the last time you disposed of an appliance, television, or computer. Did you dispose of it properly? [5 Levels, "Definitely Yes" to "Definitely No"] When you do dishes by hand, or dry dishes, how often do you use reusable/cloth dish towels/dish rag? [5 levels, "Always/Almost Always" to "Never/Almost Never"]

Landscape Information and Data
The landcover data for the study was downloaded from the New Jersey DEP Bureau of GIS, using the 2012 Land Use/Land Cover dataset [41]. These data are classified based on a modified version of the Anderson system. The data for roadway locations in this study came from the New Jersey Department of Transportation [53], incorporated with roadways maps from the DEP [41].
The locations for parks was based on maps generated in the Grant F. Walton Center for Remote Sensing and Spatial Analysis (CRSSA), which maintains maps of this type for New Jersey [54]. The PAD-US database [55] was used to check some locations, however, due to enhanced levels of specificity, ground truthing, and local knowledge of the CRSSA cartographers of this region, conflicts generally were defaulted to the map from CRSSA.
From these data sources, we were able to determine the landcover characteristics of each municipality within the study area. This included the landcover as given in the state landcover map [41], as well as these features collapsed into more general categories-all water-based cover, all forest cover, all residential cover, and similar categories along these lines. All of these measures were available at a resolution of at least 30 m; some were available at the level of 10 m. Additionally, we quantified the road network density in each municipality using the state road data [41,53]. These measures included the length of roadways per unit land, the number of intersections per unit land, and the number of intersections per unit roadway. These measures were based on accessibility measures pioneered by Forsyth et al. in their studies of accessibility in the Twin Cities area [23,[56][57][58][59][60][61].
Because this study had a focus on the influence of recreational opportunities, we also calculated a number of distance and accessibility metrics for parks. These included parks as a percentage of land cover, as well as the number and proportion of parcels within given distances of parks. Furthermore, we used the Network Analyst toolbox to calculate the number and proportion of parcels within a given distance of parks along roadways.
For each of the land characteristics, values were calculated for within municipal boundaries, as well as for a 9.6 km radius around the municipalities' geographic centers. The longest distance that can be traveled from the center of any municipality to its edge was slightly over 9.5 km, so 9.6 km was used such that this would provide a standardized area to account for any effects that may be dependent on the size of the municipality. Landcover metrics were calculated for the total area, as well as using the land area only, as some of the municipalities in the study area are coastal and made up significantly of water bodies.
A total of over 350 landscape and land use measures were available based on the data set used; 38 explanatory variables were included in this analysis that were hypothesized to have explanatory power. The features of landscape included: population, population density, total road length, number of intersections, intersection per meter roadway, number of parcels, mean parcel size, mean parcel-park distance, area of parks, percentage of parks. The following were calculated to provide measures of park accessibility the percentage of parcels within the following distances to a park both direct line and along the roadway network: 0.5 km, 1 km, 1.5 km, and 2 km. The following were calculated both at the level of total area including water, and by percentage of land area only: agricultural land, barren land, commercial/industrial land, forestland, other vegetated land, recreational land, residential land, the number of intersections. For use in the models, the land-cover measures in the previous sentence were all based on a 9.6 km circle around the municipality, as described above. Three measures of the percentage of water were also used, using only interior water, interior and coastal water, and interior water and coastal water only out 0.1 km, as some of the municipalities were over 80% coastal water if state maps were used.
The survey instrument collected respondents' locations to the municipality. We did not know where a respondent lived within their municipality.

Analysis
All geospatial analyses were run using ArcMap version 10; all statistical analyses were run using SAS version 9.4.
Preliminary analysis indicated that it would not make sense to treat these clusters of behaviors as units. When we treated them this way-each respondent receiving a score from the three survey items within a cluster-the scores all tended to clump in the middle. Therefore, these did not act as viable treatment clusters, as it would be difficult to discern any differences among groups exposed to differing landscapes. Thus, this left us to model each behavior within a cluster on its own.
Initially, we attempted various forms of multiple linear regression to model these data. Using stepwise model selection, we attempted to model each of the nine ERBs from the set of landscape-based variables. The resulting models from this work were uninformative and/or of little to no practical use. For example, none of the nine ERBs could be modeled in such a way that produced a linear regression model with an R-squared value greater than 0.5, and most produced models with an R-squared value of 0.2 or greater-and wherein the model p-value was greater than 0.1. This occurred even when using 25 or more landscape features as explanatory variables. Additionally, general and generalized models were attempted as well. However, these also produced models that used many (20 or more) explanatory variables, without adequate explanatory power. Furthermore, given that the research goals of this study were to build tools that could potentially be used in planning and outreach, the generalized non-linear multiple regression models we constructed appeared to be over-parameterized. Additionally, they were difficult to interpret in a way that was meaningful and could be connected back to the landscape and/or people in question.
At this point in our analysis, we realized that building these linear regression models had not produced effective or meaningful results. As the responses were already non-continuous, we took the step of re-binning them from individual response levels into a binary response. This would allow us to build binary logistic regression models instead of multiple linear regression models.
Since the models built above did not perform well, we recoded the responses to all ERBs as so that they would be binary responses. This allowed us to model the responses using binary logistic regression models. Using this method of modeling eliminated some of the differentiation-that is to say, differences between people who reported that the definitely engaged in a behavior compared to those who were pretty sure that they engaged in a behavior. Due to the sample size, in the end, this might produce models that were usable and had some level of statistical meaning, even if they were somewhat less differentiated.
The coding for each ERB was as follows: For number of organizations, number of people written, and car types, a positive response in any number or amount was coded as "1" and a response of none was coded as "0". For Self-ID as an environmentalist, at either of the two levels of agreement were coded "1", neutral were left out, and in either of two levels of disagreement were coded as "0". For appliance efficiency, the top two levels of importance were coded as "1" and the three lowest levels were coded as "0". For environmentally friendly groceries and using a reusable dish towel, those indicating they did this "Always/Almost Always", or "Most of the Time" were coded as "1" whereas those who reported "half the time", "Rarely" or "Never/Almost Never" were coded as "0". Persons indicating that they would pay 50% or greater more for sustainable food were coded "1", all else were coded as "0" (10% and none). In proper appliance disposal, those saying "Definitely Yes" and "Probably Yes" were coded as "1" whereas those reporting "Definitely Not", "Probably Not" or "Might or might not" were coded as "0".
We recognize that this recoding process grouped together respondents who had different responses, and that this may have, on some items, led to an over or under representation of respondents who engaged in or did not engage in a behavior. When interpreting the results of these models, they should not be read as absolutes, and should be read as an indicator of the proportion of respondents engaging in the behavior, but possibly not to the fullest extent of the behavior.
The set of landscape variables used was reduced from all available information-approximately 350 variables including a very wide range of landcover classifications-down a smaller set of variables of interest. The landscape characteristics within this subset were based on features that were related to use of parks-including percentage of parks, distance to parks, proportion of residents living near parks, and park accessibility. Additionally, we included canopy cover, as this has been found in previous studies to be related to environmental literacy [16,44], and related land cover variables-percentage forests, and proximity to water features. Housing density was included as well, since a focus of this investigation was how differences in landscape may drive these behaviors across suburban and urban landscapes, and this measure may serve as an indicator of the extent of urbanization in a given location within the study area. These logistic regression models used stepwise model selection to find the best-fitting model with at each level of explanatory variables from one to twenty-five.
In summation, the logistic regression models we built for this study used the nine ERBs identified in the clusters as the dependent variables. The independent variables were landscape variables, as described above. Using the Stepwise function in SAS, each of the dependent variables could be modeled from any subset of these landscape variables. The stepwise function selected the independent variables that produced the logistic regression model with the highest correct classification rate. A full list of landscape variables and their description is available in the Supplementary information section.

Logistic Regression Models
A model that can accurately classify responses as 0 or 1 from the explanatory variables is the desired outcome in this type of logistic regression using a binary response variable. There is no universally agreed upon threshold for a logistic regression model that performs adequately, however, the lower end of acceptable accurate classification seen in most studies is 80% [62][63][64][65]. As producing a planning tool was an objective for this research, having a model that could correctly classify respondents as engaging ("1") or not engaging ("0") in a given behavior at a high rate was seen as an important outcome.
We modeled the responses using logistic regression. To do this, we recoded most of the responses, as described below, and used only those landscape variables that had appeared in the linear regression models. All of these were coded as events or non-events, in order to use binary logistic regression using the logit link function, as specified in the documentation for SAS 9.4: where α signifies the intercept parameter and β = (β 1. , . . . , β k ) signifies the vector of s slope parameters [66].
With logistic regression, one way of measuring the usefulness of a model is how often it is able to classify the responses correctly. A model with a 50% classification rate correctly classifies only half of the data [67]. Using the LOGISTIC function in SAS, a model can be specified with a 50% classification rate using an intercept [66]. While there are no universally accepted thresholds for what constitutes a "good" level of classification, for a model to be useful or provide some insight, generally 75-80% is seen as a lower threshold in that this increases the minimum classification by half again [62,65].

Description of Study Area
Two-hundred and ninety responses for this survey were gathered, due to budgetary constraints; this represented 0.035% of the population of the study area. Not all respondents answered all questions, making the number used for analysis fewer than 290 in some cases. This made it difficult to detect effects, especially where the number of responses at a given level were quite low.
Additionally, there were 25 municipalities in the study area. Of these, respondents came from 23 municipalities. Due to the way responses were gathered and priced, it was not feasible to guarantee collection of responses for all municipalities and still gather a large enough pool of responses. The number of responses ranged from one to thirty-one. Table 2 below gives the range of some of the most general landscape characteristics. The study area is a Middlesex County, New Jersey. The total area is approximately 835 square kilometers, with less than five percent of that area consisting of water [41]. The earliest European settlement of the area dates to the 1680s, and is currently the second most populous county in New Jersey with a population of approximately 840,000 [42]. The Raritan River is a major waterway, and the landscape is primarily flat, with the highest elevation at 90 m above sea level [41]. Middlesex County is part of the New York City metropolitan area.
Middlesex County contains a mixture of municipalities, ranging from small and highly urbanized areas such as Perth Amboy, to large and sparsely populated areas such as Cranbury township [41,42,54]. We had respondents from 23 of the 25 municipalities within Middlesex County, with none from Jamesburg or Helmetta. Population in these municipalities ranges from approximately 100,000 residents in both Edison and Woodbridge to 3800 in Cranbury (or 2200 in Helmetta, where there were no responses) [42]. The smallest municipality with responses is Dunellen, at 2.7 square kilometers, up to Monroe, at 108.8 square kilometers. Despite being entirely urban [42], Middlesex County is approximately 9.5% agricultural land, with major activities including crops of hay, corn, soybeans, fruits and vegetables, as well as greenhouse, floriculture, nursery, and livestock operations [68].

Logistic Regression Models
As we can see from Table 3, three of the nine ERBs had models where the correct classification rate was <0.6, indicating correct classification less than 60% of the time. Of the remaining six, five had correct classification rates between 60% and 70%. Classification is defined as the ability of the model to correctly identify responses. In this case, it would mean the ability of the model to, based on landscape variables, identify whether the ERB response was negative or positive ("0" or "1"). This means that these models could only accurately identify a "0" or "1" from the data 60-70% of the time, which is not typically accepted as a level at which a model is effective [62,65]. Furthermore, if the goal is to use this for management or some type of implementation, these models would not provide adequate guidance to decision-makers. Only Appliance Disposal has a correct classification rate above 0.8 (or 80%), at 0.84-meaning that, using this model, a response of "0" or "1" can be correctly identified 84% of the time. This model used road length, parcels within 1km of a park, and the area of lakes as the predictor landscape variables. Road length is positive, indicating that as total length of roads increased, respondents were more likely to have properly disposed of their last appliance. Percentage of area in lakes and parcels within 1km of a park were negative, indicating that as these increased, respondents were more likely to have disposed of their last appliance improperly.
In total, only ten landscape variables were used in these models: Forest Cover, Intersections/land area, Mean parcel to park distance, Parcels within 1 km of park %, Percent of Area in Streams, Percentage of land in Lakes, Population Density, Recreation land area, Total amount of water, Total Road Length.
In summation, the logistic regression model for appliance disposal has a correct identification rate of 84%, as discussed below. The remaining 8 ERBs we modeled had correct classification rates of less than 68%.

Appliance Disposal Model
The logistic regression model for respondents saying that they disposed of their last appliance properly used three explanatory variables and can be described using Table 4 below. Stating "definitely" or "probably" that a respondent properly disposed of an appliance was coded as an event "1", whereas responding that they "definitely" or "probably" did not dispose of an appliance properly was coded as a "0". Based on this, as the length of roadway within the municipality increased, respondents were more likely to properly dispose of an appliance. As the percentage of parcels within 1km of a park and the percentage of lakes as the landcover increased, proper appliance disposal decreased.

Implications
Based on the research we performed, it does not appear that features of landscape are a large influence on this set of ERBs within this population. It would seem that responses to these ERB items can be modeled neither individually nor in aggregate using this set of landscape characteristics as explanatory variables. Therefore, this study provides some limited insight into how planners, managers, stewards or policy makers might attempt to manage the landscape as to promote these behaviors.
It did seem that there is some consistency in which landscape features provided some explanatory power for logistic regression. Only ten landscape variables were used in these models. This logistic regression models built for these behaviors from landscape characteristics were weak. Given this, these metrics or similar metrics may serve as a potential starting point for future research in this direction. The fact that there was some consistency among the landscape variables may indicate that these would serve as better predictors at a different scale or level of granularity.
The landscape variables included in these models do coincide to some extent with the work done on the role of landscape in influencing behavior. Forest cover, percent of area in streams, total amount of water, and percentage of land in lakes are all included in these ten landscape features. All three of these support the findings to this point that seeing certain colors, particularly green and blue, might influence individuals to feel more connected to their landscape and be in a better mood [4,38,39,69].
In addition to this, five of the variables among these ten are related to availability and accessibility of parks and recreation opportunities-percentage of recreational land, intersections per unit land area, mean parcel to park distance, percentage of parcels within one kilometer of a park, and total road length. All five of these align with earlier research in this field that indicates that spending time outside can influence behaviors-including health related behaviors, functions of identity, sense of place, and individual and community well-being [16,44,[70][71][72][73][74][75] Furthermore, this highlights the importance of accessibility to parks and recreation in promoting these behaviors [37,[76][77][78][79][80][81][82]. These finding may provide some confirmatory evidence to framings by earlier authors citing landscape as an explanatory factor for ERBs [1][2][3][4][5].
Perhaps, the most obvious reason for our findings is that there is no relationship at all between landscape and ERBs-as least when it comes to this population. Previous work in the field has linked individuals' built environments to a range of behaviors-include health and wellness, physical activity, property values, employment, and wealth [6,[22][23][24][25]27,37,44,50,58,59,76,[83][84][85][86]. However, in the face of the literature to this point, and the findings of this study, it seems more likely that there are issues with the study size, methodology, and landscape granularity.
It is possible, however, that the specific behaviors studied-writing to officials, identity, organizational membership, grocery purchasing, dish towel use, appliance purchasing and disposal, and car interest-are not specifically linked to landscape characteristics at such a granular level. For example, a survey of the wider region with landscape features calculated more broadly might show a relationship. This would indicate that the relevant level at which residents experience or use the landscape contains a wider area than used in this study. Any differences in ERB engagement due to landscape might only be detectable on a larger scale and when the differences in the landscape are quite great-comparing densely-populated, highly-developed New Jersey to White Pine County Nevada which is about the same size but has only 10,000 people [42].

Discussion of Appliance Disposal Model
The binary logistic regression model for Appliance Disposal seemed to be the only one for which landscape attributes were able to provide useful classification. For the other 8 ERBs, the level of correct classification was under 70%. This is not particularly useful to policy makers or stewards.
The landscape characteristics used for modeling appliance disposal do not seem like inherently logical explanatory variables for this behavior. Proximity to a park may be related to overall greater engagement in ERBs. Potentially, road length serves as some proxy for the ability of residents to surreptitiously dispose of appliances-persons living in areas with fewer roads may be more able to discreetly dispose of waste improperly for example. That is to say, this may be an indicator of an area being more rural overall. The connection to the amount of lakes might possibly indicate that lakes are a site of improper appliance disposal, and that having more available makes it easier to dispose of appliances improperly. To test this, further work would need to be done interviewing persons improperly disposing of appliances.
Since proper appliance disposal presented the only potentially useful model, we analyzed this a bit further. Using the logistic regression scoring, we ranked each municipality. The scores tanged from 5.71 to 23.19 (least to most likely to dispose of an appliance properly). The two towns ranked highest (23.19 and 23.09 respectively) both have mobile phone applications where residents can schedule pickup of appliances, find information on roadside pickup days, and weekend hours for appliance drop-off at municipal centers. Comparatively, these towns seem to make it easier for residents to dispose of large appliances.
The four lowest ranked towns (all <7, 5.71, 6.34, 6.69, 6.98), all use a county-based collaborative program for their waste-management. None has a mobile application. One has a monthly bulk goods pickup day, but the rest are by appointment only. Additionally, the websites for three of these have not been updated in over a year, and there appears to be a dearth of information for residents. Two of these municipalities mention a drop-off center for bulk goods, but it is not located within the municipality, and only one lists weekend hours. Overall, it would seem that these municipal appliance disposal regulations line up with some features of landscape. It is unclear, however, what the direct link might be between the landscape features used in this model, and how they specifically relate to municipal waste policies. Future work in linking specific municipal waste policies to features of landscape might be instructive, as would studying the origins of municipal waste policies.

Limitations and Future Directions
Future work attempting to study connections between features of the landscape and engagement in ERBs could be designed somewhat differently in order to produce more coherent results, and allow for more in-depth analysis. Primarily, a greater number of responses stratified at the municipal level, and tagging respondents more specifically to their location would produce more robust results.
There may be some issues in detecting effects based on the wide range of sizes, response rates, and population within this study area. A stratified sample with a response rate of one per thousand residents, with a minimum of 30 respondents/municipality would aide in fixing many of the issues related to the small sample size.
The issues with sample size would be less of a problem if respondents were tagged more specifically to their location in the landscape. Due to concerns of anonymity and privacy, as well as technological limitations, respondents were only identified to the level of the municipality. If these methods were used to increase the specificity of respondent locations, spatial modeling methods-such as kriging, point-pattern analysis, and spatial auto-correlation models-could be attempted in order to explain differences in respondent behaviors.
Future work of ours will take both of these considerations into account when designing the study. We would recommend that researchers, stewards and planners attempting to do similar studies also implement these procedures, particularly tagging respondents within landscape-focused studies more specifically with their location.
The landscapes features found to be most useful in these models provide some future direction for future studies on this topic. It may be instructive to take a deeper look at the specific distribution of recreational land in relation to respondents' residences as well as the type and distribution of water within the landscape, for example.
The limited findings of this work may provide some instructive quality in their limitations. Policy makers may attempt some intervention, or make a claim as towards the impact from some set of regulation or policy change. In order to quantify, or evaluate the effectiveness of these changes, having good baseline data is important, as is having a suitable ground for comparison. Perhaps the information from this study could be used as a starting point or baseline for proposed interventions within the study region. Going forward, interventions could be measured against these data. If nothing else, this could also serve as a call for greater data collection to be used for evaluation of present and proposed future policies.

Conclusions
This study did not conclusively link eight environmentally responsible behaviors to features of the landscape. One behavior-proper appliance disposal-could be accurately classified 84% of the time using landscape characteristics as explanatory factors, the other eight could not be explained this way, nor could any of them be adequately modeled using linear regression based on landscape characteristics. Overall, our work did seem to suggest accessibility to parks, water features, and overall network accessibility may play some minor role in influencing behavior, which is consistent with literature to this point.
The findings of this study could be used to steer further research into ERB encouragement through land management. Additionally, the landscape features found to be associated with ERBs may be a useful starting place for the further analysis of the role of landscape characteristics in influencing individuals' environmental behaviors.
Future research in this direction could perhaps use more regular behavior tracking surveys-perhaps daily and weekly applications-in order to paint a clearer picture of who engages in what particular behaviors, where, and at what time. Additionally, planners and policy-makers instituting a policy change may want to track the behavior for a period of time before the change and after the change in order to get a clearer picture of how the policy has impacted ERBs. More in-depth interviews coupled with tracking would also help in uncovering potential motivating factors for ERBs.
Alternatively, these results could be interpreted to mean that any people engaging in ERBs (and this set of nine ERBs in particular) can, and do, live everywhere within the study region, and are not confined to landscapes with a particular set of characteristics.
Supplementary Materials: The following are available online at http://www.mdpi.com/2413-8851/2/4/114/s1, Complete Survey used in Research, list of all behaviors studied in research, a more in-depth explanation of the canonical correlation analysis don in conjunction with this work, and a full list of landscape features considered in analysis.
Author Contributions: Both authors worked jointly on development of the survey instrument, piloting and deployment. D.C. carried out the statistical analysis, obtained spatial data and cleaned data. R.J. performed validation of analysis. Both authors worked jointly on the manuscript.