Onsite wastewater treatment systems (OWTSs), or septic tank systems, are commonly used throughout the United States and are generally effective at remediating wastewater. However, malfunctioning OWTSs can introduce excess nutrients (i.e., nitrogen and phosphorous) and pathogens (i.e., E. coli) into the environment. There is increasing evidence that OWTSs can be a significant, and potentially underestimated, nonpoint source (NPS) of pollution. Thus, the objectives of this research were to (1) develop a model to assess the pollution potential from OWTSs using GIS-based multi-criteria decision analyses (MCDA) and (2) evaluate the relationship between the pollution potential from OWTSs and water pollutants. This study was completed in the Choccolocco Creek watershed, Alabama. The main tributary in this watershed, the Choccolocco Creek, is an impaired waterbody due to elevated E. coli concentrations. An MCDA was developed to model the pollution potential from OWTSs using environmental and OWTS variables. Similarly, an OWTS site unsuitability analysis, that only included environmental variables, was used to predict where OWTS may poorly perform, if OWTS data are not accessible in other areas. Water samples were taken along Choccolocco Creek to measure nitrogen, phosphorous, and E. coli concentrations. Pollutant concentrations were correlated to modeled pollution potential from OWTSs and OWTS site unsuitability, to compare how the exclusion of OWTS data changes the results. Additionally, land cover distribution was correlated to pollutant concentrations to account for other potential NPSs of water pollution. All water pollutants were significantly, positively correlated to OWTS count. Additionally, E. coli and nitrogen concentrations were significantly, positively correlated to pollution potential from OWTSs. This suggests that OWTSs may contribute to water pollution within the watershed. Furthermore, the location of areas most probable to have OWTS pollution varied between models, highlighting the importance of accounting for OWTSs as a NPS of water pollution. The methods presented could be adapted for other watersheds and used to guide best watershed management practices.
The current contribution presents wintertime climatology from 2012 to 2020 of mixed-phase clouds and their radiative effects when coupled to the sea ice states. Measurements from the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) at the North Slope Alaska (NSA) site in Utqiagvik, Alaska are being analyzed.Classification of cloud hydrometeors in the liquid, ice or mixed-phase states was primary determined by the Cloudnet algorithm, developed by the Finish Meteorological Institute, and applied to a set of ground-based remote sensing instruments from NSA . To evaluate the influence by sea ice, which plays an important role on the Arctic surface-atmosphere interaction, the statistics are separated into cases when clouds are coupled or decoupled to specific sea ice conditions, like presence of leads or polynyas in the vicinity of NSA .We found that clouds coupled to sea ice with presence of leads have shown distinguished features like the increase of total liquid content, lower cloud base heights and less ice content when compared to decoupled cases. Nevertheless, these results rely on Cloudnet accurately detecting cloud droplets within mixed-phase clouds. Arctic cloud radiative effects (CRE) have already been studied from short expeditions like the SHEBA campaign (Shupe et al., 2004) and middle-term ground observations in Barrow (Shupe et al., 2015) and Ny-Ålesund, Svalbard (Ebell et al., 2020). We extend similar CRE studies for 8 years during wintertime, where longwave up- and down-welling flux measurements from NSA are used to estimate surface net fuxes and other cloud radiative features for cases when clouds are coupled or decoupled to sea ice conditions and their sensitivity to different gradients of air-surface temperature when leads or polynyas are present.
United States Federal Emergency Management Agency provides model-output localized flood grids that are useful in characterizing flood hazards for properties located in the Special Flood Hazard Area (SFHA ─ areas expected to experience a 1% or greater annual chance of flooding). But these flood grids are often unavailable or fail to include return periods for particular applications, such as understanding flood risk of properties during the 70-year useful building life cycle. Furthermore, due to the unavailability of higher-return-period flood grids, the flood risk of properties located outside the SFHA cannot be quantified. Here, we present a method to estimate the flood hazard for U.S. properties that are located both inside and outside the SFHA. The flood hazard is characterized by the Gumbel extreme value distribution to project flood elevations to extreme flood events for which an entire area is assumed to be submerged. Spatial interpolation techniques impute elevation values in the extreme flood elevation surfaces and therefore can estimate the flood hazard for areas outside the SFHA. The proposed method can improve the assessment of flood risk for properties located in both inside and outside the SFHA and therefore, the decision-making process regarding flood insurance purchases, mitigation strategies, and long-term planning for enhanced resilience to one of the world’s most ubiquitous natural hazards.
ORCiD: https://orcid.org/0000-0003-4699-3733 Rootstocks are gaining importance in viticulture as a strategy to combat abiotic challenges, as well as enhancing scion physiology and attributes. Therefore, understanding how the rootstock affects photosynthesis is insightful for genetic improvement of either genotype in the grafted grapevines. Photosynthetic parameters such as maximum rate of carboxylation of RuBP (Vcmax) and the maximum rate of electron transport driving RuBP regeneration (Jmax) has been identified as ideal targets for breeding and genetic studies. However, techniques used to measure these photosynthetic parameters are time consuming and subjective to leaf level which is complex to implement at field scale. Hyperspectral remote sensing uses the optical properties of the entire vine to predict photosynthetic capacity at canopy level. In this study, estimates of Vcmax and Jmax were assessed using different machine learning models: PLS (Partial least Squares), LR (Least Angle Regression), LASSO (Least Absolute Shrinkage and Selection Operator), PCR (Principle Component Regression) based on leaf reflectance metrics obtained with hyperspectral wavelength ranging from 400 to 1000nm. Prediction models were developed for six different rootstock genotypes with common scion Marquette considering three different sampling dates carried out in Brookings, South Dakota in 2021. Preliminary results indicate that each rootstock has distinctly different Vcmax and Jmax profiles across the season. From the model assessment, PLS was found to have robust prediction of Vcmax with R2 of 0.53 and for Jmax with R2 of 0.63. Multiple year trials will be used to validate precise and rapid quantification of photosynthesis using hyperspectral remote sensing.
Root exudation refers to the processes by which plants release compounds called root exudates into the soil. These exudates are primarily carbon-containing compounds that interact with microbial communities in the rhizosphere. Microbial consumption of exudates reduces the concentration of the exudated compounds in the soil, causing the plant to exude more of those compounds. Currently, there is limited understanding of the interaction between plant-root exudation mechanisms and the surrounding microbial communities. Among the Sorghum Association Panel (SAP), an established and genetically characterized sorghum diversity panel, we observed a spectrum of root colors (tan, yellow, red, purple-brown, black) identical to the range of observed sorghum seed colors. Previous studies examining differentially expressed metabolites between colorful seeds showed that flavonoids and anthocyanins were higher in dark seeds than white seeds. Root color is genotype-dependent and consistent over time. We hypothesized that the observed color diversity of sorghum roots was due to differential metabolite profiles in the root exudates across genotypes. We designed an experiment to collect exudates from 15 genotypes (n=60). After three weeks of growth, sorghum roots were washed and submerged in ultrapure water for 24 hours. The hydroponic solution was filtered and incubated with methanol. The whole root system was also ground after exudation. The root exudate solutions and the ground-up roots underwent either HILIC and RPLC analysis to separate and detect polar and hydrophobic metabolites. Through metabolite profiling of root exudates, we aim to identify sorghum genotypes that more efficiently allocate carbon below ground via their root systems.
Hemp (Cannabis sativa L. <0.3% THC) is a versatile crop cultivated for grain, fiber, and cannabinoids used for health and wellness purposes. Following the passage of the 2018 Farm Bill, which removed hemp from the controlled substances list, there is potential for expanded hemp acreage and a concomitant need to breed cultivars with suitable agronomic performance in US growing regions. Understanding the relationships between morphological, biochemical, and spectral traits with respect to yield will allow for high-throughput phenotyping and advanced breeding efforts. For five years, the Cornell Hemp research team has evaluated high-cannabinoid hemp cultivars in replicated field trials, with populations originating from both commercial sources and the Cornell Hemp Breeding Program. These field trials evaluated plants for phenotypic traits, including plant height, morphology, flowering time, cannabinoid concentration, and total biomass yield. We have also determined key morphometric measurements that are correlated with end-of-season biomass yield. Also, both floral and foliar cannabinoid samples were correlated with end-of-season whole plant biomass cannabinoid concentration. To improve phenotyping efficiency, an unmanned aerial vehicle (UAV)-based multispectral system was deployed to characterize morphological and biochemical traits over time. These datasets are being used to develop high-throughput phenotyping methods to predict biomass yield, and in the future cannabinoid concentration and flowering time.
A dry beans (Phaseolus vulgaris L.) cultivar must fit the environment in which it will be grown. Therefore, days to maturity (DM) is the most important physiological component affecting yield and grain quality outcomes. Additionally, dry bean stand count (SC) at early growth stages estimation provides useful information for agronomic decision-making and can measure root rot loss due to damping-off. The visual inspection to determine the accurate maturity date and the final number of emerged plants is labor-intensive, time-demanding, and tedious. Therefore, there is an increasing demand for alternative approaches to estimating DM and SC in a high-throughput phenotyping mode (HTP). In this study, we developed a Deep Learning (DL) HTP pipeline to capture the sequential behavior of time series data for estimating DM and to identify target plants in the early growth stage for SC estimation using field dry bean data obtained from aerial RGB images at the plot-level. A state-of-the-art hybrid model combining Convolutional Neural Networks (CNN) and Long Short-Term Memory (LSTM) was used to extract DM features and capture the sequential behavior of time series data. Faster R-CNN object detection method was deployed to SC. The DL model to estimate DM was tested on five different environments across years, and SC was done comparing different ground sample resolutions in two trials. Results suggest the effectiveness of the CNN-LSTM and Faster R-CNN models employed compared to traditional methods. Furthermore, this study highlighted the technical parameters that can influence the DL model results in the breeding program decision-making.
Using Unoccupied Aerial Systems (UASs) to collect image data for field phenotyping enables the coverage of large spatial scales, but the extraction of phenotypes at the level of an individual plant remains a challenge. Phenotypes are generally extracted at the plot level after manual delineation of the plots in an orthomosaic image. However, this approach is impractical when individual plant phenotypes are of interest, such as in mutant screenings, heterozygous population trials, and other cases where individual plants in a plot or a landscape may be very diverse. For mutant studies in particular, UASs have had limited utility since one typically wants to identify individual plants with outlier phenotypes. Current mutant population screenings require scientists to walk through the field on a regular basis to manually identify and record plants of interest. If UAS software could be designed to identify these plants, much larger studies could be conducted with lower labor costs. Here, we use deep Convolutional Neural Network (CNN) models to detect individual maize plants in raw UAS images followed by the extraction of individual phenotypes. The predictions from the CNN models are used to identify plants in the orthomosaic image to derive information on plot location. For the extraction of phenotypes, plant pixels are segmented from the background and traits including color-based indices and canopy coverage are calculated. Individual plant height is derived from the corresponding Digital Surface Model (DSM). This project will contribute to a more efficient approach for the phenotyping of individual plants using UAS images.
ORCiD: [0000-0001-6507-9985 of Jinyoung Y. Barnaby], [0000-0001-9082-6583 of Scott E. Warnke] Precise assessment of large mapping populations, comprising a few thousand plants including replications (a prerequisite step for breeding) is time-consuming and labor-intensive. Furthermore, phenotyping results tend to be variable and subjective depending on who is doing the scoring. One way to overcome these limitations is by collecting more data in the form of digital images, and precisely evaluating phenotypic variation in stress severity as well as temporal progression of stress symptoms within the population through machine learning methods. 230,400 images representing temporal progression of drought stress symptoms of interspecific turfgrass hybrid mapping population were processed using Python OpenCV and NumPy packages for noise removal, edge-preserving smoothing, color space conversion, contrast enhancement, and identification mapping. Then machine learning-based algorithms and models were developed not only to quantify stress severity but also to monitor temporal progression rate of stress symptoms. Hierarchical clustering was then performed to assess genotypic variation in stress progression. Such machine learning-based high-throughput digital phenotyping platforms can significantly increase the success of quantitative trait locus mapping and candidate gene identification to develop potential molecular markers that will assist in a faster characterization of germplasm to ultimately breed for stress resilient cultivars.
A high-throughput image analysis pipeline was developed to facilitate root phenotyping by reducing time-consuming labeling while maintaining phenotyping accuracy. This pipeline leverages a deep learning-based tool named SLEAP (SLEAP Estimates Animal Poses) which is designed to automate the detection of distinct morphological landmarks. By training SLEAP to detect the root branch points, tips, and midline of each root imaged in a gel cylinder, we were able to robustly and efficiently recover the root system geometry. We trained models to identify these landmarks on primary, lateral, and seminal roots across a range of crop plants, including soybean, rice, canola, and pennycress. We find that our SLEAP models are robust across genotypes and experiments, enabling automated root system quantification at the rate of hundreds of plants per hour. Using predictions of root landmark locations, we developed Python-based pipelines to extract phenotypic traits, including tip depths, root lengths, convex hulls, root angles, measures of curviness, and lateral root distribution (available at https://github.com/talmolab/sleap-roots). In order to extract meaningful patterns from this high-dimensional description of plant phenotypes, we use machine learning-based methods for dimensionality reduction and manifold embedding, allowing us to capture the statistical structure of root phenotypes present in our screens. In future work, we will use these quantitative phenotypic traits as a predictor for root system traits that enhance carbon sequestration capabilities in genome-wide association studies.
Recent advancements in proximal remote sensing have increased the spatial and temporal resolution of data collection, as well as the availability of these technologies for applications to precision agriculture. These sensors have allowed the collection of new and large quantities of data, which have been used to successfully determine phenotypes and parametrize crop growth models. So far, these data streams have been mostly used separately, though they contain unique structural, spatial, and spectral information. Thus, this research aims to integrate these disparate data sources to improve estimations of agronomically important crop traits. In this study, we examine two high-throughput and relatively inexpensive remote platforms: unoccupied ground vehicles (UGV) and unoccupied aerial vehicles (UAV). Data were collected on maize hybrids from the Genomes to Fields initiative over 5 years, from 2018 to 2022, in Aurora, NY. We used ground rovers to collect lidar scans, which were converted to point clouds, to construct the three-dimensional sub-canopy architecture of maize plants. Multispectral sensors, covering red, green, blue, red-edge, and near infrared (NIR) were deployed on a UAV platform to characterize maize canopies. Machine learning methods, including autoencoders, will be used to extract latent phenotypes from the lidar point clouds and multispectral images. Ultimately, these will be used to predict manually measured traits, such as yield, in order to compare the prediction accuracies of models using these measurements separately and jointly.
Resource allocation drives the above-ground distribution of mass in grass plants across discrete developmental units called phytomers. Although the number of phytomers varies in genetically-identical grasses, there frequently isn't an associated variance in some summary phenotypes. To understand what may be driving this, we tracked the growth of 30 S. italica plants from genotypes B100 and A10.1. We experimentally observed that plants from the genotype B100 had between 20 and 22 phytomers, while plants from the genotype A10.1 had between 7 and 9 phytomers. B100 plants with more phytomers (e.g., 22) did not grow taller or have more total leaf length, despite having more leaves than plants with fewer phytomers (e.g., 20). A10.1 plants with more phytomers (e.g., 9) did grow taller and had more total leaf length than those with fewer phytomers (e.g., 7). We developed a dynamical model to determine if these patterns are emergent from the underlying growth structure. The model is parameterized using the number of phytomers and related developmental time parameters: leaf emergence, stem and leaf elongation time, panicle emergence, and flowering time. The model uses the semi-sequential nature of phytomer growth as its structure. The model predicts that differences in timing of the shift to reproductive growth could explain the patterns observed. Experimental measurements suggest this shift is primarily due to tuning the developmental time parameter controlling the units contributing to the stem, rather than leaves.
Models of fault slip development generally consider interfacial strength to be frictional and deformation of the bounding medium to be elastic. The frictional strength is usually considered as sliding rate- and state-dependent. Their combination, elastic deformation due to differential slip and rate-state frictional strength, leads to nonlinear partial differential equations (PDEs) that govern the spatio-temporal evolution of slip. Here, we investigate how data on fault slip rate and stress can directly discover the complex system (of PDEs) that governs aseismic slip development. We first prepare (synthetic) data sets by numerically solving the forward problem of slip rate and fault stress evolution with models, such as a thin laterally deformable layer over a thick substrate. We now identify the variables, for example, slip rate or friction state variable, and use nonlinearity identification algorithms to discover the governing PDE of the chosen variable.In particular, we use sparse identification of nonlinear dynamics algorithm (SINDy, Brunton et al., 2016) where we solve a regression problem, Ax=y. Here, y is the time derivative of the variable of interest, for example, slip rate. A is a large matrix (library) with all possible candidate functions that may appear in the slip rate evolution PDE. The entries in x, to be solved for, are coefficients corresponding to each library function in matrix A. We update A according to the solutions x so that A's column space can span the dynamics we seek to find. To find the suitable column space for A, we encourage sparse solutions for x, suggesting that only a few columns in matrix A are dominant, leading to a parsimonious representation of the governing PDE.We show that the algorithm successfully recovers the terms of the PDE governing fault slip and could also find the frictional parameter, for example, a/b, where a and b, respectively, are the magnitudes that control direct and evolution effects. Moreover, the algorithm can also determine whether the associated state variable evolves as aging- or slip-law types or their combination. Further, with the data set prepared from distinct initial conditions, we show that the SINDy can also determine the problem parameter’s spatial distribution (heterogeneities) from fault slip rate and stress data.