The low sensitivities and relatively low R2 values are less than those found in sand or clay soils

While circle-packing nearly describes our model problem, there is one major caveat: no physical justification prevents the circles from overlapping one another. This ’soft boundary’ makes it possible to achieve 100% coverage of the domain by allowing overlap. If the only objective was to maximize the effective areal coverage of the field, then one could distribute sensors next to one another without discretion. However, the monetary cost of sensors, sensor operation, and sensor maintenance makes this approach unreasonable, which motivates our stated objective of maximizing field coverage with the fewest number of sensors possible.The outcome of this algorithm is that sensors are placed throughout the field such that sensors are placed in the largest gaps between sensors. This is done by incrementing the acceptable distance between sensors by a small amount and then making many attempts at placing a sensor before repeating the process. We score the fitness of each placement design by the ratio of the number of field pixels that are within the half-variogram range of a sensor to the total number of pixels in the field. In other words, what percentage of the field area is within the half-variogram range of one or more sensors? This process is repeated until it is impossible to place a sensor outside of the range of all other sensors in the design, or until the field is completely covered. The flowchart for this algorithm is shown in Figure 5.3A, and a schematic depicting the evolution of sensor placement in an arbitrary field shape is shown in Figure 5.3B.Soil is a complex environment. It is a three-phase medium containing organic matter, minerals, metals, ceramics, gases, water, and a host of biological life. This is all to say that many things could complicate potentiometric nitrate sensor readings. Although the nitrate sensor nodes are sensitive to nitrate and are largely insensitive to most other ions – soil is more complicated than aqueous solutions. Ideally, nitrate sensors should not be sensitive to soil properties ,plastic garden pot but calibration or direct measurement of the interfering property could help return accurate nitrate measurement values.

To characterize the performance of the nitrate sensor nodes in soil, the potential of the nodes was recorded over time in varying nitrate concentrations and moisture levels in three types of agricultural soils. After calibrating the nitrate sensor nodes in aqueous solutions, they were immediately cycled through several containers of soil saturated with aqueous solutions of varying nitrate concentrations. Finally, they were measured in several containers of dry soil with varying volumes of 10 mM nitrate solution. The nitrate sensor nodes were tested in sand, peat, and clay soils. Sand tests were performed with commercially available desert sand , consisting of only sand-sized soil particles and no initial nitrate concentration. Clay tests were performed with an agricultural clay soil utilized for perennial alfalfa and from Bouldin Island in the Sacramento–San Joaquin Delta, California . Peat tests were performed in Miracle-Gro Nature’s Care organic potting mix. To circumvent the possibility of hysteresis attributed to the exposure of one soil type before measuring the sensor node in another soil type, the sensors were only measured in a single kind of soil on a given day, and the stakes were gently but thoroughly cleaned with deionized water and dried between measurements. In the nitrate concentration experiment, aqueous solutions of 0.1 mM – 1 M NaNO3 were used to saturate an array of containers containing sand, clay, or peat soils, depending on which soil type the sensors were being tested in that day. After saturating the soil with a surplus of nitrate solution, the sensor nodes were placed in the 0.1 mM container by displacing the soil with a spoon, inserting the stake, and then gently redistributing the soil over the electrode surfaces. After the nitrate sensor node was inserted into the container, the program ran on the micro-controller, uploading the recorded potential over WiFi to an online spreadsheet once per minute. Once the potential of the sensor stabilized, it was then moved to the 1 mM container, and so on. Figure 5.11 shows the nitrate sensor nodes being measured in the sand. After the measurements were made, the nitrate concentration of the soils was determined by taking KCl extractions and then measuring the extractions with ionexchange chromatography. The nitrate sensor nodes’ responses to the calibrated nitrate concentrations are plotted in Figure 5.12.

We expect similar results in measuring the sensor nodes in saturated sand as we would in aqueous solutions. This is because sand has a low cation exchange capacity, so few other ionic species are present. Furthermore, uncharged solids in the soil are unlikely to interfere with the potentiometric measurement.The nitrate sensor nodes in the sand showed strong linear relationships between nitrate concentration and recorded potential, with about half of the 17 sensor nodes included in this trial showing R2 values above 0.99. Eight sensor nodes with R2 > 0.99 were used for the remaining analysis. The average sensitivity for these sensors was -42 ± 8 mV/dec. Figure 5.12A shows the linear relationship between nitrate concentration and output potential for five of the sensors in the sand, with sensitivities ranging from -40 to -42 mV/dec. In clay, five of 14 sensors had R2 values > 0.9, and the sensitivity of the sensor nodes was -39 ± 8 mV/dec mV/dec, shown in Figure 5.12C. Clay has a much finer grain size and higher cation exchange capacity than sand. This likely explains the loss in sensitivity when measuring in clay compared to an aqueous solution or saturated sand. This suggests that soil texture alone is not a primary interfering factor in soil nitrate measurements, though other soil characteristics may be. If this is the case, calibration of the sensor in different soil types would be necessary to deploy such devices in agricultural applications. Figure 5.12B shows the output potential of three nitrate sensor nodes in peat soil at varying nitrate levels. The nitrate sensor nodes demonstrated sensitivities of -31 ± 8 mV/dec,draining pots with R2 values for each sensor’s best fit line of 0.8, 0.71, and 0.99.We suspect this is because the sensors had become damaged by the time of these measurements. The sensors were measured in sand and clay soils before the measurements in peat, so by the time of these measurements, the sensors had been inserted, removed, and rinsed from soil media many times, and it is possible that despite our best efforts, the sensing element may have become damaged. Compared to Figure 4.10, which plots ‘fresh’ gold electrode nitrate sensors in the same soil type measured with a Campbell Scientific data logger, the sensitivity was -47 mV/dec, R 2=0.95, and E0 variation 30 mV . The different results depending on the age of the sensor and the electronics used to measure it highlights the importance of improved stability for real-world use cases.

To find the impact of moisture content on the sensor nodes, 10 mM nitrate solution was used to water an array of containers containing sand, clay, or peat soils, depending on which soil type the sensors were being tested in that day. The containers were watered to 0 – 50% volumetric water content in 10% increments. The sensor nodes were initially placed in 0% VWC soil while the program ran on the micro-controller, uploading the recorded potential over WiFi to an online spreadsheet once per minute. Once the potential of the sensor node stabilized, it was then moved to the 10% VWC container, and so on. Ideally, the sensors’ output signal should not depend on soil moisture content. However, potentiometric sensors require ionic contact between the two electrodes to function: ions must be able to move freely between the ISE and RE, proportional to the finite current associated with potentiometric measurements. If soil does not hold enough water to support the flow of ions, the sensor becomes an open circuit, and there would be no signal. The results are shown in Figure 5.13. The potential abruptly increases at low moisture content because the micro-controller is programmed to report high signal output at open-circuit inputs. However, above a certain moisture threshold, sand and clay soils follow the expected pattern: relatively constant potential with respect to moisture content. For sand, the threshold is between 10 and 20% volumetric water content VWC, while for clay, it is between 20 and 30%. This makes sense because clay’s matrix potential is higher than sand’s, meaning water is bound more tightly to the surfaces of solid particles in clay, lowering the likelihood of an ionic pathway forming between the ISE and the RE of the nitrate sensor. The relationship between water content and sensor signal output in peat soil is less clear. The minimum water threshold seems to be between 10-20%, but the output potential is not as stable between 30-50% VWC. This could be attributed to the possibility of damage buildup in the nitrate sensors, or it could indicate that different sensors have different minimum thresholds. It could also be explained by the variations of water retention in high-organic matter soil. Further investigation is warranted.As noted earlier, there is significant E0 variability from one sensor node to another. Up until now, the behavior of the nitrate sensor nodes has been plotted altogether. However, when we look at individual sensors, E0 is relatively consistent, indicating that the variability has to do mainly with the sensor nodes rather than the properties of the soil. Individual sensor nodes have relatively consistent E0 values across the different trials. Subsets of the same 20 sensors were used in the six plots shown in Figures 5.12 and 5.13. The sensors with low R2 in the sand also had low R2 in clay, while most sensors with high R2 in the sand also had high R2 in clay . Similarly, sensors with relatively high potential outputs had this characteristic across soil types and nitrate/moisture measurements. Figure 5.14 shows the response of two sensor nodes measured in sand and clay soils. Sensor A has a much lower potential than Sensor B in all cases. A common misconception surrounding ‘green materials’ is that if something will decompose or degrade into a naturally-occurring material , it is safe to deploy into the environment. This is an incomplete way of thinking because, by that definition, every material on the planet is ‘green.’ Stop reading for a moment and take a look around you. Everything you see would eventually break down into naturally-occurring materials on a long enough time scale devoid of human intervention. It was made of materials sourced from our planet, after all. The time it takes for something to degrade is also essential and is the basis of many different standards surrounding bio-degradation. Different organizations worldwide have developed standards that benchmark different degrees of how things degrade. Table 5.2 highlights some of the more widely known standards. Interestingly, these standards vary widely in the description, and only the compostability standards explicitly state the timescale at which those materials degrade. The standards in Table 5.2 also lack to describe of how these materials might degrade. To design devices that degrade in a controlled manner, one must first understand the different ways a material might degrade. There are three ways in which materials degrade. These are compositional or micro-structural changes, time-dependent deformation and associated damage accumulation, and environmental attack. Generally, materials degrade by some combination of several mechanisms. There are few opportunities for time-dependent deformation in soil, though several mechanisms catalyze compositional changes or attack material bonds. These primary mechanisms by which materials degrade in soil are discussed here. For most materials, microbial activity is the primary contributor to the degradation rate. Microbial and enzymatic digestion describes the degradation carried out by microorganisms and naturally-occurring enzymes. Microbes can degrade most – if not all – naturally occurring organic chemicals and convert them to inorganic end products to supply the microbes with nutrients and energy. As with other life forms, the material’s molecular bonds are broken to release energy, transformed into a less thermodynamically energetic material, and excreted from the microbe’s system. Intuitively, different materials degrade at different rates. Generally, more complex molecular structures degrade on longer time scales than simpler ones. Also, different microbes and enzymes are preferential to different molecular compounds.

The coordinate system used to define the parent model was NAD83 California Albers

The Sacramento River is a vital source of water for landowners and irrigation water suppliers in the Colusa Subbasin in addition to providing 80 percent of inflow to the Sacramento-San Joaquin Delta . The Glenn-Colusa Canal originates from the Sacramento River and flows through the Colusa Subbasin into Colusa County via the local canal system, providing irrigation water for the Glenn-Colusa Irrigation District, which is the largest irrigation district in the Sacramento Valley . The only principal aquifer in the Colusa Subbasin is contained in the freshwater-bearing sediments and stream channel deposits of the Modesto Formation, Tehama Formation, Riverbank Formation and Tuscan Formation, with the majority of the fresh groundwater in the project area contained within the Tehama Formation . Soils conducive to recharge, based on the SAGBI rating, range from excellent to very poor suitability, with the most suitable soils in the project area falling within the good to poor suitability range .Hydrologic conditions of interest in the project area include a natural cone of depression in groundwater elevation levels in the middle to southern half of the parent model domain . Based on recent groundwater monitoring data from the California Department of Water Resources , seasonal variations in groundwater elevation contours show a cone of depression forming in the center of the parent model domain during the Fall 2020 season due to stresses to the aquifer, such as excessive groundwater pumping for agriculture during the summer. Groundwater levels usually recover in the springtime, depending on how much precipitation occurred during the winter, but Spring 2021 contours still show lower groundwater levels in the same area with steeper groundwater head gradients,plastic garden pot as opposed to other areas in the parent model domain. Groundwater typically flows away from topographically high areas like the foothills of the Coast Ranges in the west, away from the Sacramento River in the east, and flows towards the topographically low areas where the cone of depression forms.

This area is heavily influenced by agriculture, and high amounts of pumping for agricultural use may be the cause of groundwater level depletion and the subsequent formation of the cone of depression and areas of low groundwater elevation.In 2019 and 2020, TNC recruited several landowners into their fall season bird habitat and ag-MAR incentive program. A total of 8 field sites were enrolled in TNC’s incentive program, which are all located in Colusa County, California . The nearest major cities to the ag-MAR sites are Williams, Colusa, and Meridian. The field sites modeled in this study are located near the city of Meridian, along the western bank of the Sacramento River . The fields are part of the Davis Ranches property. This area was chosen for the child model domain because Field 2 was the only field site that participated in flooding during consecutive years . The crop type grown on these fields is rice, but the fields were fallow during the period of flooding during the fall. The SAGBI recharge rating for the area in and around Field 2 and Field 15 are good to moderately good suitability for deep percolation in the soils .The data collected and used for the two groundwater models were primarily extracted from the Sacramento Valley Groundwater-Surface Water Simulation Model , which is similar to the California Central Valley Groundwater-Surface Water Simulation Model . Both SVSim and C2VSim were developed by the California Department of Water Resources, but SVSim is more specific to the northern Sacramento Valley, while C2VSim encompasses the entire Central Valley. Other sources of data include the California DWR, the United States Geological Survey , and Davids Engineering, Inc., the consulting company hired by TNC to conduct the field methods and technical analysis of the recharge program at each field site during 2019 and 2020. The data extracted from these sources were used in the MODFLOW model development process for parameterizing the packages for the parent and child models. More information on the development and calibration of C2VSim and SVSim can be found in Brush et al., 2013, and California Department of Water Resources, 2022, respectively. Layer stratigraphy data was extracted from SVSim, as were the layer and aquifer parameters, such as specific storage , specific yield , and vertical and horizontal hydraulic conductivity .

The layer stratigraphy and parameter data were used to define the aquifer parameters of all 9 groundwater model layers in both the parent and child models in their discretization packages. The uppermost layers in both models are relatively thin, with an average thickness ranging from 10 to 50 m. The lower layers are the thickest, ranging from about 50 to 200 m thick. The total thickness of all layers combined in both models is approximately 600 m . Each field site was flooded with a depth of 4 inches of water maintained for 30 continuous days during the fall, and 60 days of deep percolation data was collected , with deep percolation rates calculated for each site . These calculations were made using a mass balance approach – in which the inflows to the system are equal to the outflows – along with groundwater level measurements that were taken either periodically or continuously throughout the flooding period at groundwater monitoring wells adjacent to the field sites . Recharge took place at Field 2 during two periods, 9/18/2019-11/16/2019 and 9/17/2020-11/16/2020. Recharge took place at Field 15 during one period, 9/18/2019-11/16/2019. The average recharge rates at the three field sites ranged from 0.011 to 0.023 m/day . Deep percolation rates for all other days, or model stress periods, in the remainder of the domain are rates averaged by month from SVSim. Davids Engineering also monitored groundwater levels in four wells adjacent to the field sites . We used the groundwater level data provided by Davids Engineering for our child model’s head observation package. The average depth to groundwater between all four of the monitoring wells ranged from 5.3 to 7.7 m of depth . The software used for the models constructed for this study is the MODFLOW-2005 Three-Dimensional Finite-Difference Ground-Water Model developed by the U.S. Geological Survey . MODFLOW-2005, hereafter referred to as MODFLOW, uses packages that simulate the effects on groundwater flow processes of wells, rivers, lakes, and other relevant aspects and boundary conditions. A list of the packages used in both models is included in Table 4. MODFLOW can also vertically simulate aquifer systems with different geologic layers that may be specified as confined or unconfined. MODFLOW’s mathematical solution for simulating groundwater flow through the center of each cell in a model follows the partial-differential groundwater flow equation .The larger of the two models developed for this study was a steady state model that encompassed most of the area of Colusa County . As is the case with steady state models, all model parameters remained constant over a specified amount of time,draining pots so all input data that varied over time were averaged to one value for input.

The purpose for making the large model steady state was to have all inputs to the system equal to the outputs, stabilizing the system as a whole. This provided a more general simulation of the regional hydrologic processes , which allowed for faster model runs and processing times. The results of this large steady state model provided information on flow magnitude and direction that was used to constrain the boundaries of the child model. The domain for the large model spans an area of approximately 160,000 square meters and includes 9 layers, modeled after the layers in SVSim.The parent model has 140 rows and 113 columns, with a cell size of 400 m × 400 m. The number of stress periods, or days, is 1 since the parent model is run under steady state conditions. The boundary conditions used to define the model domain include 2 general head boundaries , 2 constant head boundaries , and the river package to simulate the effects of rivers, streams, and canals that flow through the model domain . The smaller of the two models developed for this study was a transient model that was located within the area of the parent model domain, including the ag-MAR field sites . Since the child model is transient, all data input for each package varied over time. The starting date of the child model was defined to be 8/18/2019, one month before the flooding occured, and the end date of the model is 11/16/2020, the last day of the flooding period. Like the parent model, the coordinate system used to define the child model was NAD83 California Albers. The child model has 58 rows and 58 columns, and each cell in the model has a size of 100 m × 100 m. The number of stress periods, or days, included in the child model is 456. Groundwater fluxes and directions of flow were extracted from the results of the parent model in order to inform the boundary conditions of the child model . The well package was used to simulate those fluxes. We defined the flux from each well in all model cells on the child model’s boundaries to represent the general flow that occurred in the same area of the parent model. In more specific terms, cell-by-cell flow data was extracted for each cell in every layer from the parent model, and the cell-by-cell flows were used to inform the groundwater fluxes at the boundaries of the child model. We also used the River package to represent the small part of the Sacramento River that acts as the eastern boundary of the child model. However, the RIV package is only present in the first layer of the model, so we used the GHB package for the eight layers directly below the river cells. The recharge package was used to represent the 2019 and 2020 ag-MAR field sites, with deep percolation data from Davids Engineering.Three different scenarios were run using the child model to test different parameters and their influence on the hydrologic response seen in the model outputs. The three scenarios were designed to assess whether the timing, frequency, and amount of recharge significantly influences the hydrologic results of the child model. Scenario 1 was the baseline model run, using the original and unaltered deep percolation data from Davids Engineering. No alterations or calibrations to the model inputs were made for this scenario. For Scenario 2, we increased deep percolation rates from Davids Engineering by one order of magnitude to ascertain whether we could see a significant hydrologic response. Only the deep percolation data was altered for Scenario 2; all other inputs remained the same as Scenario 1. It should also be clarified that a tenfold increase in recharge rates is not realistic. However, we wanted to test this hypothetical extreme solely to see a significant response. For Scenario 3, we increased the duration of recharge to span ten years.

We chose to use the 2019 deep percolation data, and repeated that year ten times to see how increasing the duration of recharge would affect the child model results. The 2019 deep percolation data input for Scenario 3 was unaltered; only the duration of recharge and the time discretization of the child model was increased to span ten years. Scenario 3 was designed to represent ten consecutive years of implementing ag-MAR on these field sites.The results analyzed from the steady-state parent model include model performance plots and model outputs. We assessed the performance of the parent model with a one-to-one plot of the observed heads that were input into the model, with the simulated equivalent heads that the model calculated at each groundwater monitoring well in the parent model domain . The different colors on the plot represent the binned data of the residual heads, which were calculated as the difference between the observed and simulated equivalent heads. The different bins represent different groups of residual data ranges. The gray line represents the line of equal values, or one-to-one line, between the observed and simulated equivalent heads.The same binned residual head data were plotted on a map of the parent model domain, showing the location of the groundwater monitoring wells and the residual head value . Like the residual head bins plotted in the one-to-one plot in Figure 12, the bins for the residual map were calculated from the difference between the observed and simulated equivalent head values at each groundwater monitoring well on the map.