Engaging the next generation of editorial talent through a hands-on fellowship model.

Peer-review and subject-matter editing is the backbone of scientific publishing. However, early-career researchers (ECRs) are given few opportunities to participate in the editorial process beyond reviewing articles. Thus, a disconnect exists: science needs high-quality editorial talent to conduct, oversee and improve the publishing process, yet we dedicate few resources to building editorial talent nor giving ECRs formal opportunities to influence publishing from within. ECRs can contribute to the publishing landscape in unique ways given their insight into new and rapidly developing publishing trends (e.g. open science). Here, we describe a two-way fellowship model that gives ECRs a "seat" at the editorial table of a field-leading journal. We describe both the necessary framework and benefits that can stem from editorial fellowships for ECRs, editors, journals, societies, and the ​broader scientific community.

The role of phosphorus and nitrogen on chlorophyll a: Evidence from hundreds of lakes.

The effect of nutrients on phytoplankton biomass in lakes continues to be a subject of debate by aquatic scientists. However, determining whether or not chlorophyll a (CHL) is limited by phosphorus (P) and/or nitrogen (N) is rarely considered using a probabilistic method in studies of hundreds of lakes across broad spatial extents. Several studies have applied a unified CHL-nutrient relationship to determine nutrient limitation, but pose a risk of ecological fallacy because they neglect spatial heterogeneity in ecological contexts. To examine whether or not CHL is limited by P, N, or both nutrients in hundreds of lakes and across diverse ecological settings, a probabilistic machine learning method, Bayesian Network, was applied. Spatial heterogeneity in ecological context was accommodated by the probabilistic nature of the results. We analyzed data from 1382 lakes in 17 US states to evaluate the cause-effect relationships between CHL and nutrients. Observations of CHL, total phosphorus (TP), and total nitrogen (TN) were discretized into three trophic states (oligo-mesotrophic, eutrophic, and hypereutrophic) to train the model. We found that although both nutrients were related to CHL trophic state, TP was more related to CHL than TN, especially under oligo-mesotrophic and eutrophic CHL conditions. However, when the CHL trophic state was hypereutrophic, both TP and TN were important. These results provide additional evidence that P-limitation is more likely under oligo-mesotrophic or eutrophic CHL conditions and that co-limitation of P and N occurs under hypereutrophic CHL conditions. We also found a decreasing pattern of the TN/TP ratio with increasing CHL concentrations, which might be a key driver for the role change of nutrients. Previous work performed at smaller scales support our findings, indicating potential for extension of our findings to other regions. Our findings enhance the understanding of nutrient limitation at macroscales and revealed that the current debate on the limiting nutrient might be caused by failure to consider CHL trophic state. Our findings also provide prior information for the site-specific eutrophication management of unsampled or data-limited lakes.

Ecological prediction at macroscales using big data: Does sampling design matter?

Although ecosystems respond to global change at regional to continental scales (i.e., macroscales), model predictions of ecosystem responses often rely on data from targeted monitoring of a small proportion of sampled ecosystems within a particular geographic area. In this study, we examined how the sampling strategy used to collect data for such models influences predictive performance. We subsampled a large and spatially extensive data set to investigate how macroscale sampling strategy affects prediction of ecosystem characteristics in 6,784 lakes across a 1.8-million-km2 area. We estimated model predictive performance for different subsets of the data set to mimic three common sampling strategies for collecting observations of ecosystem characteristics: random sampling design, stratified random sampling design, and targeted sampling. We found that sampling strategy influenced model predictive performance such that (1) stratified random sampling designs did not improve predictive performance compared to simple random sampling designs and (2) although one of the scenarios that mimicked targeted (non-random) sampling had the poorest performing predictive models, the other targeted sampling scenarios resulted in models with similar predictive performance to that of the random sampling scenarios. Our results suggest that although potential biases in data sets from some forms of targeted sampling may limit predictive performance, compiling existing spatially extensive data sets can result in models with good predictive performance that may inform a wide range of science questions and policy goals related to global change.

Identifying and characterizing extrapolation in multivariate response data.

Faced with limitations in data availability, funding, and time constraints, ecologists are often tasked with making predictions beyond the range of their data. In ecological studies, it is not always obvious when and where extrapolation occurs because of the multivariate nature of the data. Previous work on identifying extrapolation has focused on univariate response data, but these methods are not directly applicable to multivariate response data, which are common in ecological investigations. In this paper, we extend previous work that identified extrapolation by applying the predictive variance from the univariate setting to the multivariate case. We propose using the trace or determinant of the predictive variance matrix to obtain a scalar value measure that, when paired with a selected cutoff value, allows for delineation between prediction and extrapolation. We illustrate our approach through an analysis of jointly modeled lake nutrients and indicators of algal biomass and water clarity in over 7000 inland lakes from across the Northeast and Mid-west US. In addition, we outline novel exploratory approaches for identifying regions of covariate space where extrapolation is more likely to occur using classification and regression trees. The use of our Multivariate Predictive Variance (MVPV) measures and multiple cutoff values when exploring the validity of predictions made from multivariate statistical models can help guide ecological inferences.

Team climate mediates the effect of diversity on environmental science team satisfaction and data sharing.

Scientific research-especially high-impact research-is increasingly being performed in teams that are interdisciplinary and demographically diverse. Nevertheless, very little research has investigated how the climate on these diverse science teams affects data sharing or the experiences of their members. To address these gaps, we conducted a quantitative study of 266 scientists from 105 NSF-funded interdisciplinary environmental science teams. We examined how team climate mediates the associations between team diversity and three outcomes: satisfaction with the team, satisfaction with authorship practices, and perceptions of the frequency of data sharing. Using path analyses, we found that individuals from underrepresented groups perceived team climate more negatively, which was associated with lower satisfaction with the team and more negative perceptions of authorship practices and data sharing on the team. However, individuals on teams with more demographic diversity reported a more positive climate than those on teams with less demographic diversity. These results highlight the importance of team climate, the value of diverse teams for team climate, and barriers to the full inclusion and support of individuals from underrepresented groups in interdisciplinary science teams.

Do lakes feel the burn? Ecological consequences of increasing exposure of lakes to fire in the continental United States.

Wildfires are becoming larger and more frequent across much of the United States due to anthropogenic climate change. No studies, however, have assessed fire prevalence in lake watersheds at broad spatial and temporal scales, and thus it is unknown whether wildfires threaten lakes and reservoirs (hereafter, lakes) of the United States. We show that fire activity has increased in lake watersheds across the continental United States from 1984 to 2015, particularly since 2005. Lakes have experienced the greatest fire activity in the western United States, Southern Great Plains, and Florida. Despite over 30 years of increasing fire exposure, fire effects on fresh waters have not been well studied; previous research has generally focused on streams, and most of the limited lake-fire research has been conducted in boreal landscapes. We therefore propose a conceptual model of how fire may influence the physical, chemical, and biological properties of lake ecosystems by synthesizing the best available science from terrestrial, aquatic, fire, and landscape ecology. This model also highlights emerging research priorities and provides a starting point to help land and lake managers anticipate potential effects of fire on ecosystem services provided by fresh waters and their watersheds.

Spatial and temporal variation of ecosystem properties at macroscales.

Although spatial and temporal variation in ecological properties has been well-studied, crucial knowledge gaps remain for studies conducted at macroscales and for ecosystem properties related to material and energy. We test four propositions of spatial and temporal variation in ecosystem properties within a macroscale (1000 km's) extent. We fit Bayesian hierarchical models to thousands of observations from over two decades to quantify four components of variation - spatial (local and regional) and temporal (local and coherent); and to model their drivers. We found strong support for three propositions: (1) spatial variation at local and regional scales are large and roughly equal, (2) annual temporal variation is mostly local rather than coherent, and, (3) spatial variation exceeds temporal variation. Our findings imply that predicting ecosystem responses to environmental changes at macroscales requires consideration of the dominant spatial signals at both local and regional scales that may overwhelm temporal signals.

Geographic patterns of the climate sensitivity of lakes.

Climate change is a well-recognized threat to lake ecosystems and, although there likely exists geographic variation in the sensitivity of lakes to climate, broad-scale, long-term studies are needed to understand this variation. Further, the potential mediating role of local to regional ecological context on these responses is not well documented. In this study, we examined relationships between climate and water clarity in 365 lakes from 1981 to 2010 in two distinct regions in the northeastern and midwestern United States. We asked (1) How do climate-water-clarity relationships vary across watersheds and between two geographic regions? and (2) Do certain characteristics make some lakes more climate sensitive than others? We found strong differences in climate-water-clarity relationships both within and across the two regions. For example, in the northeastern region, water clarity was often negatively correlated with summer precipitation (median correlation = -0.32, n = 160 lakes), but was not correlated with summer average maximum temperature (median correlation = 0.09, n = 205 lakes). In the midwestern region, water clarity was not related to summer precipitation (median correlation = -0.04), but was often negatively correlated with summer average maximum temperature (median correlation = -0.18). There were few strong relationships between local and sub-regional ecological context and a lake's sensitivity to climate. For example, ecological context variables explained just 16-18% of variation in summer precipitation sensitivity, which was most related to total phosphorus, chlorophyll a, lake depth, and hydrology in both regions. Sensitivity to summer maximum temperature was even less predictable in both regions, with 4% or less of variation explained using all ecological context variables. Overall, we identified differences in the climate sensitivity of lakes across regions and found that local and sub-regional ecological context weakly influences the sensitivity of lakes to climate. Our findings suggest that local to regional drivers may combine to influence the sensitivity of lake ecosystems to climate change, and that sensitivities among lakes are highly variable within and across regions. This variability suggests that lakes are sensitive to different aspects of climate change (temperature vs. precipitation) and that responses of lakes to climate are heterogeneous and complex.

LAGOS-NE: a multi-scaled geospatial and temporal database of lake ecological context and water quality for thousands of US lakes.

Understanding the factors that affect water quality and the ecological services provided by freshwater ecosystems is an urgent global environmental issue. Predicting how water quality will respond to global changes not only requires water quality data, but also information about the ecological context of individual water bodies across broad spatial extents. Because lake water quality is usually sampled in limited geographic regions, often for limited time periods, assessing the environmental controls of water quality requires compilation of many data sets across broad regions and across time into an integrated database. LAGOS-NE accomplishes this goal for lakes in the northeastern-most 17 US states.LAGOS-NE contains data for 51 101 lakes and reservoirs larger than 4 ha in 17 lake-rich US states. The database includes 3 data modules for: lake location and physical characteristics for all lakes; ecological context (i.e., the land use, geologic, climatic, and hydrologic setting of lakes) for all lakes; and in situ measurements of lake water quality for a subset of the lakes from the past 3 decades for approximately 2600-12 000 lakes depending on the variable. The database contains approximately 150 000 measures of total phosphorus, 200 000 measures of chlorophyll, and 900 000 measures of Secchi depth. The water quality data were compiled from 87 lake water quality data sets from federal, state, tribal, and non-profit agencies, university researchers, and citizen scientists. This database is one of the largest and most comprehensive databases of its type because it includes both in situ measurements and ecological context data. Because ecological context can be used to study a variety of other questions about lakes, streams, and wetlands, this database can also be used as the foundation for other studies of freshwaters at broad spatial and ecological scales.

Unexpected stasis in a changing world: Lake nutrient and chlorophyll trends since 1990.

The United States (U.S.) has faced major environmental changes in recent decades, including agricultural intensification and urban expansion, as well as changes in atmospheric deposition and climate-all of which may influence eutrophication of freshwaters. However, it is unclear whether or how water quality in lakes across diverse ecological settings has responded to environmental change. We quantified water quality trends in 2913 lakes using nutrient and chlorophyll (Chl) observations from the Lake Multi-Scaled Geospatial and Temporal Database of the Northeast U.S. (LAGOS-NE), a collection of preexisting lake data mostly from state agencies. LAGOS-NE was used to quantify whether lake water quality has changed from 1990 to 2013, and whether lake-specific or regional geophysical factors were related to the observed changes. We modeled change through time using hierarchical linear models for total nitrogen (TN), total phosphorus (TP), stoichiometry (TN:TP), and Chl. Both the slopes (percent change per year) and intercepts (value in 1990) were allowed to vary by lake and region. Across all lakes, TN declined at a rate of 1.1% year-1 , while TP, TN:TP, and Chl did not change. A minority (7%-16%) of individual lakes had changing nutrients, stoichiometry, or Chl. Of those lakes that changed, we found differences in the geospatial variables that were most related to the observed change in the response variables. For example, TN and TN:TP trends were related to region-level drivers associated with atmospheric deposition of N; TP trends were related to both lake and region-level drivers associated with climate and land use; and Chl trends were found in regions with high air temperature at the beginning of the study period. We conclude that despite large environmental change and management efforts over recent decades, water quality of lakes in the Midwest and Northeast U.S. has not overwhelmingly degraded or improved.

Creating multithemed ecological regions for macroscale ecology: Testing a flexible, repeatable, and accessible clustering method.

Understanding broad-scale ecological patterns and processes often involves accounting for regional-scale heterogeneity. A common way to do so is to include ecological regions in sampling schemes and empirical models. However, most existing ecological regions were developed for specific purposes, using a limited set of geospatial features and irreproducible methods. Our study purpose was to: (1) describe a method that takes advantage of recent computational advances and increased availability of regional and global data sets to create customizable and reproducible ecological regions, (2) make this algorithm available for use and modification by others studying different ecosystems, variables of interest, study extents, and macroscale ecology research questions, and (3) demonstrate the power of this approach for the research question-How well do these regions capture regional-scale variation in lake water quality? To achieve our purpose we: (1) used a spatially constrained spectral clustering algorithm that balances geospatial homogeneity and region contiguity to create ecological regions using multiple terrestrial, climatic, and freshwater geospatial data for 17 northeastern U.S. states (~1,800,000 km2); (2) identified which of the 52 geospatial features were most influential in creating the resulting 100 regions; and (3) tested the ability of these ecological regions to capture regional variation in water nutrients and clarity for ~6,000 lakes. We found that: (1) a combination of terrestrial, climatic, and freshwater geospatial features influenced region creation, suggesting that the oft-ignored freshwater landscape provides novel information on landscape variability not captured by traditionally used climate and terrestrial metrics; and (2) the delineated regions captured macroscale heterogeneity in ecosystem properties not included in region delineation-approximately 40% of the variation in total phosphorus and water clarity among lakes was at the regional scale. Our results demonstrate the usefulness of this method for creating customizable and reproducible regions for research and management applications.

Data-Intensive Science and Research Integrity.

In this commentary, we consider questions related to research integrity in data-intensive science and argue that there is no need to create a distinct category of misconduct that applies to deception related to processing, analyzing, or interpreting data. The best way to promote integrity in data-intensive science is to maintain a firm commitment to epistemological and ethical values, such as honesty, openness, transparency, and objectivity, which apply to all types of research, and to promote education, policy development, and scholarly debate concerning appropriate uses of statistics.

Lake nutrient stoichiometry is less predictable than nutrient concentrations at regional and sub-continental scales.

Production in many ecosystems is co-limited by multiple elements. While a known suite of drivers associated with nutrient sources, nutrient transport, and internal processing controls concentrations of phosphorus (P) and nitrogen (N) in lakes, much less is known about whether the drivers of single nutrient concentrations can also explain spatial or temporal variation in lake N:P stoichiometry. Predicting stoichiometry might be more complex than predicting concentrations of individual elements because some drivers have similar relationships with N and P, leading to a weak relationship with their ratio. Further, the dominant controls on elemental concentrations likely vary across regions, resulting in context dependent relationships between drivers, lake nutrients and their ratios. Here, we examine whether known drivers of N and P concentrations can explain variation in N:P stoichiometry, and whether explaining variation in stoichiometry differs across regions. We examined drivers of N:P in ~2,700 lakes at a sub-continental scale and two large regions nested within the sub-continental study area that have contrasting ecological context, including differences in the dominant type of land cover (agriculture vs. forest). At the sub-continental scale, lake nutrient concentrations were correlated with nutrient loading and lake internal processing, but stoichiometry was only weakly correlated to drivers of lake nutrients. At the regional scale, drivers that explained variation in nutrients and stoichiometry differed between regions. In the Midwestern U.S. region, dominated by agricultural land use, lake depth and the percentage of row crop agriculture were strong predictors of stoichiometry because only phosphorus was related to lake depth and only nitrogen was related to the percentage of row crop agriculture. In contrast, all drivers were related to N and P in similar ways in the Northeastern U.S. region, leading to weak relationships between drivers and stoichiometry. Our results suggest ecological context mediates controls on lake nutrients and stoichiometry. Predicting stoichiometry was generally more difficult than predicting nutrient concentrations, but human activity may decouple N and P, leading to better prediction of N:P stoichiometry in regions with high anthropogenic activity.

Honorary Authorship Practices in Environmental Science Teams: Structural and Cultural Factors and Solutions.

Overinclusive authorship practices such as honorary or guest authorship have been widely reported, and they appear to be exacerbated by the rise of large interdisciplinary collaborations that make authorship decisions particularly complex. Although many studies have reported on the frequency of honorary authorship and potential solutions to it, few have probed how the underlying dynamics of large interdisciplinary teams contribute to the problem. This article reports on a qualitative study of the authorship standards and practices of six National Science Foundation-funded interdisciplinary environmental science teams. Using interviews of the lead principal investigator and an early-career member on each team, our study explores the nature of honorary authorship practices as well as some of the motivating factors that may contribute to these practices. These factors include both structural elements (policies and procedures) and cultural elements (values and norms) that cross organizational boundaries. Therefore, we provide recommendations that address the intersection of these factors and that can be applied at multiple organizational levels.

Spatial Variation in Nutrient and Water Color Effects on Lake Chlorophyll at Macroscales.

The nutrient-water color paradigm is a framework to characterize lake trophic status by relating lake primary productivity to both nutrients and water color, the colored component of dissolved organic carbon. Total phosphorus (TP), a limiting nutrient, and water color, a strong light attenuator, influence lake chlorophyll a concentrations (CHL). But, these relationships have been shown in previous studies to be highly variable, which may be related to differences in lake and catchment geomorphology, the forms of nutrients and carbon entering the system, and lake community composition. Because many of these factors vary across space it is likely that lake nutrient and water color relationships with CHL exhibit spatial autocorrelation, such that lakes near one another have similar relationships compared to lakes further away. Including this spatial dependency in models may improve CHL predictions and clarify how well the nutrient-water color paradigm applies to lakes distributed across diverse landscape settings. However, few studies have explicitly examined spatial heterogeneity in the effects of TP and water color together on lake CHL. In this study, we examined spatial variation in TP and water color relationships with CHL in over 800 north temperate lakes using spatially-varying coefficient models (SVC), a robust statistical method that applies a Bayesian framework to explore space-varying and scale-dependent relationships. We found that TP and water color relationships were spatially autocorrelated and that allowing for these relationships to vary by individual lakes over space improved the model fit and predictive performance as compared to models that did not vary over space. The magnitudes of TP effects on CHL differed across lakes such that a 1 µg/L increase in TP resulted in increased CHL ranging from 2-24 µg/L across lake locations. Water color was not related to CHL for the majority of lakes, but there were some locations where water color had a positive effect such that a unit increase in water color resulted in a 2 µg/L increase in CHL and other locations where it had a negative effect such that a unit increase in water color resulted in a 2 µg/L decrease in CHL. In addition, the spatial scales that captured variation in TP and water color effects were different for our study lakes. Variation in TP-CHL relationships was observed at intermediate distances (~20 km) compared to variation in water color-CHL relationships that was observed at regional distances (~200 km). These results demonstrate that there are lake-to-lake differences in the effects of TP and water color on lake CHL and that this variation is spatially structured. Quantifying spatial structure in these relationships furthers our understanding of the variability in these relationships at macroscales and would improve model prediction of chlorophyll a to better meet lake management goals.

Conceptions of Good Science in Our Data-Rich World.

Scientists have been debating for centuries the nature of proper scientific methods. Currently, criticisms being thrown at data-intensive science are reinvigorating these debates. However, many of these criticisms represent long-standing conflicts over the role of hypothesis testing in science and not just a dispute about the amount of data used. Here, we show that an iterative account of scientific methods developed by historians and philosophers of science can help make sense of data-intensive scientific practices and suggest more effective ways to evaluate this research. We use case studies of Darwin's research on evolution by natural selection and modern-day research on macrosystems ecology to illustrate this account of scientific methods and the innovative approaches to scientific evaluation that it encourages. We point out recent changes in the spheres of science funding, publishing, and education that reflect this richer account of scientific practice, and we propose additional reforms.

Effects of Land Use on Lake Nutrients: The Importance of Scale, Hydrologic Connectivity, and Region.

Catchment land uses, particularly agriculture and urban uses, have long been recognized as major drivers of nutrient concentrations in surface waters. However, few simple models have been developed that relate the amount of catchment land use to downstream freshwater nutrients. Nor are existing models applicable to large numbers of freshwaters across broad spatial extents such as regions or continents. This research aims to increase model performance by exploring three factors that affect the relationship between land use and downstream nutrients in freshwater: the spatial extent for measuring land use, hydrologic connectivity, and the regional differences in both the amount of nutrients and effects of land use on them. We quantified the effects of these three factors that relate land use to lake total phosphorus (TP) and total nitrogen (TN) in 346 north temperate lakes in 7 regions in Michigan, USA. We used a linear mixed modeling framework to examine the importance of spatial extent, lake hydrologic class, and region on models with individual lake nutrients as the response variable, and individual land use types as the predictor variables. Our modeling approach was chosen to avoid problems of multi-collinearity among predictor variables and a lack of independence of lakes within regions, both of which are common problems in broad-scale analyses of freshwaters. We found that all three factors influence land use-lake nutrient relationships. The strongest evidence was for the effect of lake hydrologic connectivity, followed by region, and finally, the spatial extent of land use measurements. Incorporating these three factors into relatively simple models of land use effects on lake nutrients should help to improve predictions and understanding of land use-lake nutrient interactions at broad scales.

Building a multi-scaled geospatial temporal ecology database from disparate data sources: fostering open science and data reuse.

Although there are considerable site-based data for individual or groups of ecosystems, these datasets are widely scattered, have different data formats and conventions, and often have limited accessibility. At the broader scale, national datasets exist for a large number of geospatial features of land, water, and air that are needed to fully understand variation among these ecosystems. However, such datasets originate from different sources and have different spatial and temporal resolutions. By taking an open-science perspective and by combining site-based ecosystem datasets and national geospatial datasets, science gains the ability to ask important research questions related to grand environmental challenges that operate at broad scales. Documentation of such complicated database integration efforts, through peer-reviewed papers, is recommended to foster reproducibility and future use of the integrated database. Here, we describe the major steps, challenges, and considerations in building an integrated database of lake ecosystems, called LAGOS (LAke multi-scaled GeOSpatial and temporal database), that was developed at the sub-continental study extent of 17 US states (1,800,000 km(2)). LAGOS includes two modules: LAGOSGEO, with geospatial data on every lake with surface area larger than 4 ha in the study extent (~50,000 lakes), including climate, atmospheric deposition, land use/cover, hydrology, geology, and topography measured across a range of spatial and temporal extents; and LAGOSLIMNO, with lake water quality data compiled from ~100 individual datasets for a subset of lakes in the study extent (~10,000 lakes). Procedures for the integration of datasets included: creating a flexible database design; authoring and integrating metadata; documenting data provenance; quantifying spatial measures of geographic data; quality-controlling integrated and derived data; and extensively documenting the database. Our procedures make a large, complex, and integrated database reproducible and extensible, allowing users to ask new research questions with the existing database or through the addition of new data. The largest challenge of this task was the heterogeneity of the data, formats, and metadata. Many steps of data integration need manual input from experts in diverse fields, requiring close collaboration.

It's Good to Share: Why Environmental Scientists' Ethics Are Out of Date.

Although there have been many recent calls for increased data sharing, the majority of environmental scientists do not make their individual data sets publicly available in online repositories. Current data-sharing conversations are focused on overcoming the technological challenges associated with data sharing and the lack of rewards and incentives for individuals to share data. We argue that the most important conversation has yet to take place: There has not been a strong ethical impetus for sharing data within the current culture, behaviors, and practices of environmental scientists. In this article, we describe a critical shift that is happening in both society and the environmental science community that makes data sharing not just good but ethically obligatory. This is a shift toward the ethical value of promoting inclusivity within and beyond science. An essential element of a truly inclusionary and democratic approach to science is to share data through publicly accessible data sets.

Using uncertainty and sensitivity analyses in socioecological agent-based models to improve their analytical performance and policy relevance.

Agent-based models (ABMs) have been widely used to study socioecological systems. They are useful for studying such systems because of their ability to incorporate micro-level behaviors among interacting agents, and to understand emergent phenomena due to these interactions. However, ABMs are inherently stochastic and require proper handling of uncertainty. We propose a simulation framework based on quantitative uncertainty and sensitivity analyses to build parsimonious ABMs that serve two purposes: exploration of the outcome space to simulate low-probability but high-consequence events that may have significant policy implications, and explanation of model behavior to describe the system with higher accuracy. The proposed framework is applied to the problem of modeling farmland conservation resulting in land use change. We employ output variance decomposition based on quasi-random sampling of the input space and perform three computational experiments. First, we perform uncertainty analysis to improve model legitimacy, where the distribution of results informs us about the expected value that can be validated against independent data, and provides information on the variance around this mean as well as the extreme results. In our last two computational experiments, we employ sensitivity analysis to produce two simpler versions of the ABM. First, input space is reduced only to inputs that produced the variance of the initial ABM, resulting in a model with output distribution similar to the initial model. Second, we refine the value of the most influential input, producing a model that maintains the mean of the output of initial ABM but with less spread. These simplifications can be used to 1) efficiently explore model outcomes, including outliers that may be important considerations in the design of robust policies, and 2) conduct explanatory analysis that exposes the smallest number of inputs influencing the steady state of the modeled system.