Incorporating the benefits of vegetative filter strips into risk assessment and risk management of pesticides.

The pesticide registration process in North America, including the USA and Canada, involves conducting a risk assessment based on relatively conservative modeling to predict pesticide concentrations in receiving waterbodies. The modeling framework does not consider some commonly adopted best management practices that can reduce the amount of pesticide that may reach a waterbody, such as vegetative filter strips (VFS). Currently, VFS are being used by growers as an effective way to reduce off-site movement of pesticides, and they are being required or recommended on pesticide labels as a mitigation measure. Given the regulatory need, a pair of multistakeholder workshops were held in Raleigh, North Carolina, to discuss how to incorporate VFS into pesticide risk assessment and risk management procedures within the North American regulatory framework. Because the risk assessment process depends heavily on modeling, one key question was how to quantitatively incorporate VFS into the existing modeling approach. Key outcomes from the workshops include the following: VFS have proven effective in reducing pesticide runoff to surface waterbodies when properly located, designed, implemented, and maintained; Vegetative Filter Strip Modeling System (VFSMOD), a science-based and widely validated mechanistic model, is suitable for further vetting as a quantitative simulation approach to pesticide mitigation with VFS in current regulatory settings; and VFSMOD parametrization rules need to be developed for the North American aquatic exposure assessment. Integr Environ Assess Manag 2024;20:454-464. © 2023 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Quantifying the contribution of external loadings and internal hydrodynamic processes to the water quality of Lake Okeechobee.

The water quality of a waterbody is determined by internal hydrodynamic processes as well as external loadings. Understanding the interaction between the external loading and internal process of a waterbody is essential for efficient water management and water quality improvement. Studies and efforts have focused on water and nutrient loading from drainage watersheds, but the contribution of the waterbody's internal process to water quality is often ignored and not well documented. This study investigated how the water quality of Lake Okeechobee is controlled by external and internal factors using statistical and numerical modeling approaches. Water quantity and quality observed at the outlets of the Lake Okeechobee drainage basins and 19 monitoring sites located within the lake were statistically analyzed using multilinear regression. A three-dimensional numerical model, namely Environmental Fluid Dynamics Code (EFDC), was calibrated to the observations to mathematically represent the lake's internal hydrodynamic process. The multilinear regression found that the water quality was the most sensitive to air temperature, the total phosphorus (TP) concentration of inflow entering the lake from the Kissimmee River basins, and the amount of outflow discharged from the lake among external factors. However, the regression models and their explanatory power were substantially varied by the monitoring stations. The model parameter sensitivity analysis of the calibrated EFDC model showed that model parameters related to the lake's internal algal processes including algal growth, predation, and basal metabolism rates had greater impacts on algal biomass than other model parameters controlling nutrient-related processes such as nutrient half-saturation and hydrolysis rates. The EFDC input data sensitivity analysis found that wind (speed) is the major driving force for the internal hydrodynamic processes; its impact on algal biomass was greater than those of the external loadings. In addition, the algal biomass was found to have an inverse relationship with wind-induced horizontal currents. The results demonstrate the dynamic contribution of the internal and external drivers to the water quality of Lake Okeechobee, suggesting the need to consider both internal hydrodynamic and external loading processes for efficient water quality improvement of the lake.

Dynamic prediction of effective runoff sediment particle size for improved assessment of erosion mitigation efficiency with vegetative filter strips.

The most widely implemented mitigation measure to reduce transfer of surface runoff pesticides and other pollutants to surface water bodies are vegetative filter strips (VFS). The most commonly used dynamic model for quantifying the reduction by VFS of surface runoff, eroded sediment, pesticides and other pollutants is VFSMOD, which simulates reduction of total inflow (∆Q) and of incoming eroded sediment load (∆E) mechanistically during the rainfall-runoff event. These variables are subsequently used to calculate the reduction of pesticide load by the VFS (∆P). Since errors in ∆Q and ∆E propagate into ∆P, for strongly-sorbing compounds an accurate prediction of ∆E is crucial for a reliable prediction of ∆P. The most important incoming sediment characteristic for ∆E is the median particle diameter (d50). Current d50 estimation methods are simplistic, yielding fixed d50 based on soil properties and ignoring specific event characteristics and dynamics. We derive an improved dynamic d50 parameterization equation for use in regulatory VFS scenarios based on an extensive dataset of 93 d50 values and 17 candidate explanatory variables compiled from heterogeneous data sources and methods. The dataset was analysed first using machine learning techniques (Random Forest, Gradient Boosting) and Global Sensitivity Analysis (GSA) as a dimension reduction technique and to identify potential interactions between explanatory variables. Using the knowledge gained, a parsimonious multiple regression equation with 6 predictors was developed and thoroughly tested. Since three of the predictors are event-specific (eroded sediment yield, rainfall intensity and peak runoff rate), predicted d50 vary dynamically across event magnitudes and intensities. Incorporation of the improved d50 parameterization equation in higher-tier pesticide assessment tools with VFSMOD provides more realistic quantitative mitigation in regulatory US-EPA and EU FOCUS pesticide risk assessment frameworks. The equation is also readily applicable to other erosion management problems.

Modeling the relationship between air temperature and grapefruit quality traits.

BACKGROUND: Consumers of grapefruit require consistent fruit quality with a good physical appearance and taste. The air temperature during the growing season affects both the external (external color index (ECI)) and internal (titratable acidity (TA) and total soluble solids ratio (TSS/TA)) fruit quality of grapefruit. The objective of this study was to develop computer models that encompass the relationship between preharvest air temperature and fruit quality to predict fruit quality of grapefruit at harvest. RESULTS: There was a logarithmic relationship between the number of days with a daily minimum air temperature ≤13 °C and ECI, with a greater number of days resulting in higher ECI. In addition, there was a second-order polynomial relationship between the number of hours ≥21 °C and both TA and TSS/TA, with a greater number of hours resulting in lower TA and higher TSS/TA. Model performance for predicting the ECI, TA, and TSS/TA during 2004-05 and 2005-06 growing seasons was good, with Nash and Sutcliffe coefficient of efficiency (NSE) values for each season of 0.835 and 0.917 respectively for ECI, 0.896 and 0.965 respectively for TA and 0.898 and 0.966 respectively for TSS/TA. Applying the model to statistical survey data covering 13 growing seasons demonstrated that the TSS/TA model was robust. CONCLUSION: Statistical models were developed that predicted the development of grapefruit ECI, TA, and TSS/TA. The TSS/TA model was confirmed after application to long-term statistical survey data covering 13 growing seasons. © 2022 Society of Chemical Industry.

Experimental evidence that rill-bed morphology is governed by emergent nonlinear spatial dynamics.

Past experimental work found that rill erosion occurs mainly during rill formation in response to feedback between rill-flow hydraulics and rill-bed roughness, and that this feedback mechanism shapes rill beds into a succession of step-pool units that self-regulates sediment transport capacity of established rills. The search for clear regularities in the spatial distribution of step-pool units has been stymied by experimental rill-bed profiles exhibiting irregular fluctuating patterns of qualitative behavior. We hypothesized that the succession of step-pool units is governed by nonlinear-deterministic dynamics, which would explain observed irregular fluctuations. We tested this hypothesis with nonlinear time series analysis to reverse-engineer (reconstruct) state-space dynamics from fifteen experimental rill-bed profiles analyzed in previous work. Our results support this hypothesis for rill-bed profiles generated both in a controlled lab (flume) setting and in an in-situ hillside setting. The results provide experimental evidence that rill morphology is shaped endogenously by internal nonlinear hydrologic and soil processes rather than stochastically forced; and set a benchmark guiding specification and testing of new theoretical framings of rill-bed roughness in soil-erosion modeling. Finally, we applied echo state neural network machine learning to simulate reconstructed rill-bed dynamics so that morphological development could be forecasted out-of-sample.

Wetland hydropattern and vegetation greenness predict avian populations in Palo Verde, Costa Rica.

Many wetlands around the world that occur at the base of watersheds are under threat from land-use change, hydrological alteration, nutrient pollution, and invasive species. A relevant measure of whether the ecological character of these ecosystems has changed is the species diversity of wetland-dependent waterbirds, especially those of conservation value. Here, we evaluate the potential mechanisms controlling variability over time and space in avian species diversity of the wetlands in the Palo Verde National Park, a Ramsar Site of international importance in Costa Rica. To do so, we assessed the relative importance of several key wetland condition metrics (i.e., surface water depth, wetland extent, and vegetation greenness), and temporal fluctuations in these metrics, in predicting the abundance of five waterbirds of high conservation value as well as overall waterbird diversity over a 9-yr period. Generalized additive models revealed that mean NDVI, an indicator of vegetation greenness, combined with a metric used to evaluate temporal fluctuations in the wetland extent best predicted four of the five waterbird species of high conservation value as well as overall waterbird species richness and diversity. Black-bellied Whistling-ducks, which account for over one-half of all waterbird individuals, and all waterbird species together were better predicted by including surface water depth along with wetland extent and its fluctuations. Our calibrated species distribution model confidently quantified monthly averages of the predicted total waterbird abundances in seven of the 10 sub-wetlands making up the Ramsar Site and confirmed that the biophysical diversity of this entire wetland system is important to supporting waterbird populations both as a seasonal refuge and more permanently. This work further suggests that optimizing the timing and location of ongoing efforts to reduce invasive vegetation cover may be key to avian conservation by increasing waterbird habitat.

Parsimonious Mechanistic Modeling of Bacterial Runoff into Irrigation Ponds To Inform Food Safety Management of Agricultural Water Quality.

Pond irrigation water comprises a major pathway of pathogenic bacteria to fresh produce. Current regulatory methods have been shown to be ineffective in assessing this risk when variability of bacterial concentrations is large. This paper proposes using mechanistic modeling of bacterial transport as a way to identify improved strategies for mitigating this risk pathway. If the mechanistic model is successfully tested against observed data, global sensitivity analysis (GSA) can identify important mechanisms to inform alternative, preventive bacterial control practices. Model development favored parsimony and prediction of peak bacterial concentration events. Data from two highly variable surface water irrigation ponds showed that the model performance was similar or superior to that of existing pathogen transport models, with a Nash-Sutcliffe efficiency of 0.48 and 0.18 for the two ponds. GSA quantified bacterial sourcing and hydrology as the most important processes driving pond bacterial contamination events. Model analysis has two main implications for improved regulatory methods: that peak concentration events are associated with runoff-producing rainfall events and that intercepting bacterial runoff transport may be the best option to prevent bacterial contamination of surface water irrigation ponds and thus fresh produce. This research suggests the need for temporal management strategies. IMPORTANCE Preventive management of agricultural waters requires understanding of the drivers of bacterial contamination events. We propose mechanistic modeling as a way forward to understand and predict such events and have developed and tested a parsimonious model for rain-driven surface runoff contributing to generic Escherichia coli contamination of irrigation ponds in Central Florida. While the model was able to predict the timing of peak events reasonably well, the highly variable magnitude of the peaks was less well predicted. This indicates the need to collect more data on the fecal contamination inputs of these ponds and the use of mechanistic modeling and global sensitivity analysis to identify the most important data needs.

Seasonal dynamics of terrestrially sourced nitrogen influenced Karenia brevis blooms off Florida's southern Gulf Coast.

Harmful algal blooms (HABs) threaten coastal ecological systems, public health, and local economies, but the complex physical, chemical, and biological processes that culminate in HABs vary by locale and are often poorly understood. Despite broad recognition that cultural eutrophication may exacerbate nearshore bloom events, the association is typically not linear and is often difficult to quantify. Off the Gulf Coast of Florida, Karenia brevis blooms initiate in the open waters of the Gulf of Mexico, and advection of cells supplies nearshore blooms. However, past work has struggled to describe the relationship between terrestrial nutrient runoff and bloom maintenance near the Gulf Coast. This study applied a novel nonlinear time series (NLTS) analytical framework to investigate whether nearshore bloom dynamics observed near Charlotte Harbor, FL were causally and systematically driven by terrestrially sourced inputs of nitrogen, phosphorus, and freshwater between 2012 and 2018. Singular spectrum analysis (SSA) isolated low-dimensional, deterministic signals in K. brevis log10-density dynamics and in the dynamics of nine of 10 candidate drivers. The predominantly seasonal K. brevis signal was strong, explaining 77.6% of the total variance in the observed time series. Causal tests with convergent cross-mapping provided evidence that nitrogen concentrations measured at the discharge point of the Caloosahatchee River systematically influenced K. brevis bloom dynamics. However, further causal testing failed to link these nitrogen dynamics to an upstream basin, possibly due to data limitations. The results support the hypothesis that anthropogenic nitrogen runoff facilitated the growth of K. brevis blooms near Charlotte Harbor and suggest that bloom events would be mitigated by nitrogen source and transport controls within the Caloosahatchee and/or Kissimmee River basins. More broadly, this work demonstrates that management-relevant causal inferences into the drivers of HABs may be drawn from available monitoring records.

Importance of genetic parameters and uncertainty of MANIHOT, a new mechanistic cassava simulation model.

We identified the most sensitive genotype-specific parameters (GSPs) and their contribution to the uncertainty of the MANIHOT simulation model. We applied a global sensitivity and uncertainty analysis (GSUA) of the GSPs to the simulation outputs for the cassava development, growth, and yield in contrasting environments. We compared enhanced Sampling for Uniformity, a qualitative screening method new to crop simulation modeling, and Sobol, a quantitative, variance-based method. About 80% of the GSPs contributed to most of the variation in maximum leaf area index (LAI), yield, and aboveground biomass at harvest. Relative importance of the GSPs varied between warm and cool temperatures but did not differ between rainfed and no water limitation conditions. Interactions between GSPs explained 20% of the variance in simulated outputs. Overall, the most important GSPs were individual node weight, radiation use efficiency, and maximum individual leaf area. Base temperature for leaf development was more important for cool compared to warm temperatures. Parameter uncertainty had a substantial impact on model predictions in MANIHOT simulations, with the uncertainty 2-5 times larger for warm compared to cool temperatures. Identification of important GSPs provides an objective way to determine the processes of a simulation model that are critical versus those that have little relevance.

Comparison of empirical and mechanistic equations for vegetative filter strip pesticide mitigation in long-term environmental exposure assessments.

Recent advances in mechanistic modeling of vegetated filter strips (VFS) have made it possible to incorporate VFS mitigation into environmental exposure assessments (EEAs). However, outside of fixed efficiency approaches, there are no widely adopted and standardized procedures for incorporating VFS quantitative mitigation into long-term, higher-tier EEAs. A source of hesitation involves the use of empirical equations for predicting pesticide trapping by the VFS. A recent study evaluated existing empirical equations and a mechanistic mass-balance approach using the most extensive field database available of VFS pesticide efficiency from single-event storms. That study concluded that an updated empirical equation (Sabbagh equation) and a mechanistic mass-balance approach performed reasonably well. The objective of this research was to study the effect of upscaling the VFS trapping equations from single events into long-term EEAs. The U.S. EPA Pesticide in Water Calculator (PWC) model linked with the Vegetative Filter Strip MODeling system (VFSMOD) long-term EEA modeling framework (30 yr) was updated to incorporate the alternative trapping equations and tested VFS mitigation results under contrasting agroecological settings with varying erosion/sediment transport conditions. Differences in both acute and chronic 90th percentile estimated environmental exposure concentrations (EECs) were relatively small when comparing predictions using the four pesticide trapping equations. A global sensitivity analysis (GSA) also indicated that selection of a specific trapping equation for predicting EECs was less important than other important input factors such as the VFS length and pesticide properties. However, in terms of the percent reductions in EECs, the choice of pesticide trapping equation was as important as the VFS length. This research builds upon the conclusion of previous single-event studies that the mechanistic mass-balance and refit Sabbagh empirical equation were both valid for EEAs. The mass balance approach represents a reasonable option for regulatory agencies that prefer mechanistic approaches.

A spatiotemporal natural-human database to evaluate road development impacts in an Amazon trinational frontier.

Road construction and paving bring socio-economic benefits to receiving regions but can also be drivers of deforestation and land cover change. Road infrastructure often increases migration and illegal economic activities, which in turn affect local hydrology, wildlife, vegetation structure and dynamics, and biodiversity. To evaluate the full breadth of impacts from a coupled natural-human systems perspective, information is needed over a sufficient timespan to include pre- and post-road paving conditions. In addition, the spatial scale should be appropriate to link local human activities and biophysical system components, while also allowing for upscaling to the regional scale. A database was developed for the tri-national frontier in the Southwestern Amazon, where the Inter-Oceanic Highway was constructed through an area of high biological value and cultural diversity. Extensive socio-economic surveys and botanical field work are combined with remote sensing and reanalysis data to provide a rich and unique database, suitable for coupled natural-human systems research.

Using a coupled dynamic factor - random forest analysis (DFRFA) to reveal drivers of spatiotemporal heterogeneity in the semi-arid regions of southern Africa.

Understanding the drivers of large-scale vegetation change is critical to managing landscapes and key to predicting how projected climate and land use changes will affect regional vegetation patterns. This study aimed to improve our understanding of the role, magnitude, and spatial distribution of the key environmental and socioeconomic factors driving vegetation change in a southern African savanna. This research was conducted across the Kwando, Okavango and Zambezi catchments of southern Africa (Angola, Namibia, Botswana and Zambia) and explored vegetation cover change across the region from 2001-2010. A novel coupled analysis was applied to model the dynamic biophysical factors then to determine the discrete / social drivers of spatial heterogeneity on vegetation. Previous research applied Dynamic Factor Analysis (DFA), a multivariate times series dimension reduction technique, to ten years of monthly remotely sensed vegetation data (MODIS-derived normalized difference vegetation index, NDVI), and a suite of time-series (monthly) environmental covariates: precipitation, mean, minimum and maximum air temperature, soil moisture, relative humidity, fire and potential evapotranspiration. This initial research was performed at a regional scale to develop meso-scale models explaining mean regional NDVI patterns. The regional DFA predictions were compared to the fine-scale MODIS time series using Kendall's Tau and Sen's Slope to identify pixels where the DFA model we had developed, under or over predicted NDVI. Once identified, a Random Forest (RF) analysis using a series of static social and physical variables was applied to explain these remaining areas of under- and over- prediction to fully explore the drivers of heterogeneity in this savanna system. The RF analysis revealed the importance of protected areas, elevation, soil type, locations of higher population, roads, and settlements, in explaining fine scale differences in vegetation biomass. While the previously applied DFA generated a model of environmental variables driving NDVI, the RF work developed here highlighted human influences dominating that signal. The combined DFRFA model approach explains almost 90% of the variance in NDVI across this landscape from 2001-2010. Our methodology presents a unique coupling of dynamic and static factor analyses, yielding novel insights into savanna heterogeneity, and providing a tool of great potential for researchers and managers alike.

Effect of vegetative filter strip pesticide residue degradation assumptions for environmental exposure assessments.

Understanding and simulating the fate and transport of pesticides from a field to adjacent receiving water bodies is critical for estimating long-term environmental exposure concentrations (EECs) in regulatory higher-tier environmental exposure assessments (EEA). The potential of field mitigation practices like vegetative filter strips (VFS) to reduce pesticide pollution is receiving increasing attention. Previous research has proposed a modeling framework that links the US Environmental Protection Agency's (US-EPA) PRZM/EXAMS higher-tier EEA with a process-based VFS model (VFSMOD). This framework was updated to consider degradation and carryover of pesticide residue trapped in the VFS. However, there is disagreement on pesticide degradation assumptions among different regional EEA regulations (i.e. US or European Union), and in particular on how temperature and soil moisture dynamics may affect EECs. This research updated the VFS modeling framework to consider four degradation assumptions and determine if VFS residues and/or EECs differed with each assumption. Two model pesticides (mobile-labile and immobile-persistent) were evaluated for three distinct agroecological scenarios (continental row-crop agriculture, wet maritime agriculture, and dry Mediterranean intensive horticulture) with receiving water bodies and VFS lengths from 0 to 9m. The degradation assumption was important in long-term assessments to predict VFS pesticide residues (statistically different at p<0.01). However, due to the relatively small contribution of residues on the total pesticide mass moving through the VFS, degradation assumptions had a negligible impact on EECs. This indicates that, while important differences exist between EU or US EEAs, the choice of pesticide degradation assumption is not a main source of these differences.

Using Cluster Analysis to Compartmentalize a Large Managed Wetland Based on Physical, Biological, and Climatic Geospatial Attributes.

Hierarchical and partitional cluster analyses were used to compartmentalize Water Conservation Area 1, a managed wetland within the Arthur R. Marshall Loxahatchee National Wildlife Refuge in southeast Florida, USA, based on physical, biological, and climatic geospatial attributes. Single, complete, average, and Ward's linkages were tested during the hierarchical cluster analyses, with average linkage providing the best results. In general, the partitional method, partitioning around medoids, found clusters that were more evenly sized and more spatially aggregated than those resulting from the hierarchical analyses. However, hierarchical analysis appeared to be better suited to identify outlier regions that were significantly different from other areas. The clusters identified by geospatial attributes were similar to clusters developed for the interior marsh in a separate study using water quality attributes, suggesting that similar factors have influenced variations in both the set of physical, biological, and climatic attributes selected in this study and water quality parameters. However, geospatial data allowed further subdivision of several interior marsh clusters identified from the water quality data, potentially indicating zones with important differences in function. Identification of these zones can be useful to managers and modelers by informing the distribution of monitoring equipment and personnel as well as delineating regions that may respond similarly to future changes in management or climate.

Revealing Biotic and Abiotic Controls of Harmful Algal Blooms in a Shallow Subtropical Lake through Statistical Machine Learning.

Harmful algal blooms are a growing human and environmental health hazard globally. Eco-physiological diversity of the cyanobacteria genera that make up these blooms creates challenges for water managers tasked with controlling the intensity and frequency of blooms, particularly of harmful taxa (e.g., toxin producers, N2 fixers). Compounding these challenges is the ongoing debate over the efficacy of nutrient management strategies (phosphorus-only versus nitrogen and phosphorus), which increases decision-making uncertainty. To improve our understanding of how different cyanobacteria respond to nutrient levels and other biophysical factors, we analyzed a unique 17 year data set comprising monthly observations of cyanobacteria genera and zooplankton abundances, water quality, and flow in a bloom-impacted, subtropical, flow-through lake in Florida (United States). Using the Random Forests machine learning algorithm, an ensemble modeling approach, we characterized and quantified relationships among environmental conditions and five dominant cyanobacteria genera. Results highlighted nonlinear relationships and critical thresholds between cyanobacteria genera and environmental covariates, the potential for hydrology and temperature to limit the efficacy of cyanobacteria bloom management actions, and the importance of a dual nutrient management strategy for reducing bloom risk in the long term.

Evaluating the U.S. Food Safety Modernization Act Produce Safety Rule Standard for Microbial Quality of Agricultural Water for Growing Produce.

The U.S. Food and Drug Administration (FDA) has defined standards for the microbial quality of agricultural surface water used for irrigation. According to the FDA produce safety rule (PSR), a microbial water quality profile requires analysis of a minimum of 20 samples for Escherichia coli over 2 to 4 years. The geometric mean (GM) level of E. coli should not exceed 126 CFU/100 mL, and the statistical threshold value (STV) should not exceed 410 CFU/100 mL. The water quality profile should be updated by analysis of a minimum of five samples per year. We used an extensive set of data on levels of E. coli and other fecal indicator organisms, the presence or absence of Salmonella, and physicochemical parameters in six agricultural irrigation ponds in West Central Florida to evaluate the empirical and theoretical basis of this PSR. We found highly variable log-transformed E. coli levels, with standard deviations exceeding those assumed in the PSR by up to threefold. Lognormal distributions provided an acceptable fit to the data in most cases but may underestimate extreme levels. Replacing censored data with the detection limit of the microbial tests underestimated the true variability, leading to biased estimates of GM and STV. Maximum likelihood estimation using truncated lognormal distributions is recommended. Twenty samples are not sufficient to characterize the bacteriological quality of irrigation ponds, and a rolling data set of five samples per year used to update GM and STV values results in highly uncertain results and delays in detecting a shift in water quality. In these ponds, E. coli was an adequate predictor of the presence of Salmonella in 150-mL samples, and turbidity was a second significant variable. The variability in levels of E. coli in agricultural water was higher than that anticipated when the PSR was finalized, and more detailed information based on mechanistic modeling is necessary to develop targeted risk management strategies.

Brief history of agricultural systems modeling.

Agricultural systems science generates knowledge that allows researchers to consider complex problems or take informed agricultural decisions. The rich history of this science exemplifies the diversity of systems and scales over which they operate and have been studied. Modeling, an essential tool in agricultural systems science, has been accomplished by scientists from a wide range of disciplines, who have contributed concepts and tools over more than six decades. As agricultural scientists now consider the "next generation" models, data, and knowledge products needed to meet the increasingly complex systems problems faced by society, it is important to take stock of this history and its lessons to ensure that we avoid re-invention and strive to consider all dimensions of associated challenges. To this end, we summarize here the history of agricultural systems modeling and identify lessons learned that can help guide the design and development of next generation of agricultural system tools and methods. A number of past events combined with overall technological progress in other fields have strongly contributed to the evolution of agricultural system modeling, including development of process-based bio-physical models of crops and livestock, statistical models based on historical observations, and economic optimization and simulation models at household and regional to global scales. Characteristics of agricultural systems models have varied widely depending on the systems involved, their scales, and the wide range of purposes that motivated their development and use by researchers in different disciplines. Recent trends in broader collaboration across institutions, across disciplines, and between the public and private sectors suggest that the stage is set for the major advances in agricultural systems science that are needed for the next generation of models, databases, knowledge products and decision support systems. The lessons from history should be considered to help avoid roadblocks and pitfalls as the community develops this next generation of agricultural systems models.

Towards a new generation of agricultural system data, models and knowledge products: Design and improvement.

This paper presents ideas for a new generation of agricultural system models that could meet the needs of a growing community of end-users exemplified by a set of Use Cases. We envision new data, models and knowledge products that could accelerate the innovation process that is needed to achieve the goal of achieving sustainable local, regional and global food security. We identify desirable features for models, and describe some of the potential advances that we envisage for model components and their integration. We propose an implementation strategy that would link a "pre-competitive" space for model development to a "competitive space" for knowledge product development and through private-public partnerships for new data infrastructure. Specific model improvements would be based on further testing and evaluation of existing models, the development and testing of modular model components and integration, and linkages of model integration platforms to new data management and visualization tools.

Toward a new generation of agricultural system data, models, and knowledge products: State of agricultural systems science.

We review the current state of agricultural systems science, focusing in particular on the capabilities and limitations of agricultural systems models. We discuss the state of models relative to five different Use Cases spanning field, farm, landscape, regional, and global spatial scales and engaging questions in past, current, and future time periods. Contributions from multiple disciplines have made major advances relevant to a wide range of agricultural system model applications at various spatial and temporal scales. Although current agricultural systems models have features that are needed for the Use Cases, we found that all of them have limitations and need to be improved. We identified common limitations across all Use Cases, namely 1) a scarcity of data for developing, evaluating, and applying agricultural system models and 2) inadequate knowledge systems that effectively communicate model results to society. We argue that these limitations are greater obstacles to progress than gaps in conceptual theory or available methods for using system models. New initiatives on open data show promise for addressing the data problem, but there also needs to be a cultural change among agricultural researchers to ensure that data for addressing the range of Use Cases are available for future model improvements and applications. We conclude that multiple platforms and multiple models are needed for model applications for different purposes. The Use Cases provide a useful framework for considering capabilities and limitations of existing models and data.

Accounting for the Impact of Management Scenarios on Typha Domingensis (Cattail) in an Everglades Wetland.

The coupled regional simulation model, and the transport and reaction simulation engine were recently adapted to simulate ecology, specifically Typha domingensis (Cattail) dynamics in the Everglades. While Cattail is a native Everglades species, it has become invasive over the years due to an altered habitat over the last few decades, taking over historically Cladium jamaicense (Sawgrass) areas. Two models of different levels of algorithmic complexity were developed in previous studies, and are used here to determine the impact of various management decisions on the average Cattail density within Water Conservation Area 2A in the Everglades. A Global Uncertainty and Sensitivity Analysis was conducted to test the importance of these management scenarios, as well as the effectiveness of using zonal statistics. Management scenarios included high, medium and low initial water depths, soil phosphorus concentrations, initial Cattail and Sawgrass densities, as well as annually alternating water depths and soil phosphorus concentrations, and a steadily decreasing soil phosphorus concentration. Analysis suggests that zonal statistics are good indicators of regional trends, and that high soil phosphorus concentration is a pre-requisite for expansive Cattail growth. It is a complex task to manage Cattail expansion in this region, requiring the close management and monitoring of water depth and soil phosphorus concentration, and possibly other factors not considered in the model complexities. However, this modeling framework with user-definable complexities and management scenarios, can be considered a useful tool in analyzing many more alternatives, which could be used to aid management decisions in the future.