Evaluating Oilseed Biofuel Production Feasibility in California’s San Joaquin Valley Using Geophysical and Remote Sensing Techniques

Though more costly than petroleum-based fuels and a minor component of overall military fuel sources, biofuels are nonetheless strategically valuable to the military because of intentional reliance on multiple, reliable, secure fuel sources. Significant reduction in oilseed biofuel cost occurs when grown on marginally productive saline-sodic soils plentiful in California’s San Joaquin Valley (SJV). The objective is to evaluate the feasibility of oilseed production on marginal soils in the SJV to support a 115 ML yr−1 biofuel conversion facility. The feasibility evaluation involves: (1) development of an Ida Gold mustard oilseed yield model for marginal soils; (2) identification of marginally productive soils; (3) development of a spatial database of edaphic factors influencing oilseed yield and (4) performance of Monte Carlo simulations showing potential biofuel production on marginally productive SJV soils. The model indicates oilseed yield is related to boron, salinity, leaching fraction, and water content at field capacity. Monte Carlo simulations for the entire SJV fit a shifted gamma probability density function: Q = 68.986 + gamma (6.134,5.285), where Q is biofuel production in ML yr−1. The shifted gamma cumulative density function indicates a 0.15–0.17 probability of meeting the target biofuel-production level of 115 ML yr−1, making adequate biofuel production unlikely.


Introduction
The United States is the world's largest consumer of crude oil, resulting in two major problems: (1) low energy security and (2) high greenhouse gas (GHG) emissions. To increase national energy independence and decrease GHG emissions, the U.S. Congress enacted the Energy Independence and Security Act (EISA) in 2007 (110. P.L. 140). The EISA aims to increase the production of clean renewable fuels within the USA. Biofuels, such as biodiesel and renewable jet fuel from oilseed crops, fuel source. USDOE and USDA [6] recommend the need for greater research to identify viable biofuel feedstock production and management systems to support bio-refineries at commercially viable capacities, i.e., 115 ML yr −1 or more. However, it is unknown if sufficient oilseed production could support a biofuel conversion facility of sufficient capacity to be economically viable, i.e., cost less than $1 L −1 . Before considering an in-depth analysis of the economic viability and commercialization of biofuel, answers to basic questions are needed. The most fundamental question is, can sufficient biofuel feedstock be grown to support a biofuel conversion facility within agricultural regions of the USA, whether the Southwest, Midwest, or Southeast? All are potential agricultural areas proposed for biofuel production. More specifically, can Ida Gold mustard oilseed grow with sufficient yields on marginally productive salt-affected soils in the SJV to support a 115 ML yr −1 conversion facility?

Objective
It is the objective of this study (1) to formulate a crop yield model relating Ida Gold mustard oilseed yield to edaphic properties for the WSJV and (2) to use the crop yield model to predict the yield of Ida Gold mustard oilseed on salt-affected soils (i.e., soils with an EC e greater than 4 dS m −1 ) for evaluating the feasibility of oilseed production on marginal soils to support a 115 ML yr −1 biofuel conversion facility in the SJV.

Materials and Methods
The feasibility evaluation involves four steps: (1) development of an Ida Gold mustard oilseed yield model for marginal SJV soils using apparent soil electrical conductivity (EC a ) directed soil sampling; (2) identification of marginally productive salt-affected soils for oilseed production in the SJV; (3) development of a spatial database of edaphic factors influencing mustard yield for the SJV derived from satellite imagery and the SSURGO database and (4) applying the Ida Gold mustard oilseed yield model on marginally productive soil for the SJV and performing Monte Carlo simulations to show the range and probability of potential biofuel production in the region. Figure 1  commercially viable capacities, i.e., 115 ML yr −1 or more. However, it is unknown if sufficient oilseed production could support a biofuel conversion facility of sufficient capacity to be economically viable, i.e., cost less than $1 L −1 . Before considering an in-depth analysis of the economic viability and commercialization of biofuel, answers to basic questions are needed. The most fundamental question is, can sufficient biofuel feedstock be grown to support a biofuel conversion facility within agricultural regions of the USA, whether the Southwest, Midwest, or Southeast? All are potential agricultural areas proposed for biofuel production. More specifically, can Ida Gold mustard oilseed grow with sufficient yields on marginally productive salt-affected soils in the SJV to support a 115 ML yr −1 conversion facility?

Objective
It is the objective of this study (1) to formulate a crop yield model relating Ida Gold mustard oilseed yield to edaphic properties for the WSJV and (2) to use the crop yield model to predict the yield of Ida Gold mustard oilseed on salt-affected soils (i.e., soils with an ECe greater than 4 dS m −1 ) for evaluating the feasibility of oilseed production on marginal soils to support a 115 ML yr −1 biofuel conversion facility in the SJV.

Materials and Methods
The feasibility evaluation involves four steps: (1) development of an Ida Gold mustard oilseed yield model for marginal SJV soils using apparent soil electrical conductivity (ECa) directed soil sampling; (2) identification of marginally productive salt-affected soils for oilseed production in the SJV; (3) development of a spatial database of edaphic factors influencing mustard yield for the SJV derived from satellite imagery and the SSURGO database and (4) applying the Ida Gold mustard oilseed yield model on marginally productive soil for the SJV and performing Monte Carlo simulations to show the range and probability of potential biofuel production in the region.

Development of a Field-Based Ida Gold Mustard Oilseed Yield Model
A field experiment was conducted to identify the edaphic properties that influence oilseed yield of Ida Gold mustard. The approach of Corwin et al. [7] was followed. The approach uses geospatial electromagnetic induction measurements of EC a to direct soil sampling for determining the soil properties influencing crop yield. The approach is based on the concept that if a correlation exists between crop yield and EC a , then EC a is measuring, either directly or indirectly, one or more soil properties that are influencing the crop yield. By conducting an EC a survey to direct soil sampling, crop yield/soil sampling sites can be identified that provide a range of edaphic properties and their influence on yield.
The development of a field-based Ida Gold mustard oilseed yield model involves the following steps: (1) selection of an appropriate field that has a full range of edaphic properties that are suspected of influencing Ida Gold mustard oilseed yield; (2) conducting an intensive EC a survey of the field; (3) identification of sites within the field where crop yield and soil core samples are taken that reflect the range and variability of edaphic influences on oilseed yield; (4) analysis of chemical and physical properties of the soil cores thought to influence yield and (5) statistical analysis and crop yield model formulation.

Study site description
The study site was a 16.2-ha field (latitude-longitude coordinates: 37º02 02.97 N, 120º47 31.56 ) located west of Los Banos in Merced County, California on the WSJV (Figure 2). The site provided a range of soil properties thought to influence the yield of Ida Gold mustard oilseed. In particular, the field was characterized by a broad range of salinity and boron values. Salinity and boron are properties common to marginally productive soils in the SJV that are known to significantly influence Ida Gold mustard oilseed yield. The soil at the study site is a Britto clay loam. The soil taxonomic class is fine, smectitic, thermic Typic Natraqualfs. The Britto series consists of deep, very poorly drained soils with high concentrations of salt and alkali in the lower horizons. The soil ranges from moderately saline to strongly saline (8-16 dS m -1 ). The texture in the top 55 cm is a clay loam and 0.55-1.55 m is a sandy clay loam. The parent material is alluvium derived from sedimentary rock. The mean annual precipitation is 25 cm.

Development of a Field-Based Ida Gold Mustard Oilseed Yield Model
A field experiment was conducted to identify the edaphic properties that influence oilseed yield of Ida Gold mustard. The approach of Corwin et al. [7] was followed. The approach uses geospatial electromagnetic induction measurements of ECa to direct soil sampling for determining the soil properties influencing crop yield. The approach is based on the concept that if a correlation exists between crop yield and ECa, then ECa is measuring, either directly or indirectly, one or more soil properties that are influencing the crop yield. By conducting an ECa survey to direct soil sampling, crop yield/soil sampling sites can be identified that provide a range of edaphic properties and their influence on yield.
The development of a field-based Ida Gold mustard oilseed yield model involves the following steps: (1) selection of an appropriate field that has a full range of edaphic properties that are suspected of influencing Ida Gold mustard oilseed yield; (2) conducting an intensive ECa survey of the field; (3) identification of sites within the field where crop yield and soil core samples are taken that reflect the range and variability of edaphic influences on oilseed yield; (4) analysis of chemical and physical properties of the soil cores thought to influence yield and (5) statistical analysis and crop yield model formulation.

Study site description
The study site was a 16.2-ha field (latitude-longitude coordinates: 37º02'02.97" N, 120º47'31.56") located west of Los Banos in Merced County, California on the WSJV (Figure 2). The site provided a range of soil properties thought to influence the yield of Ida Gold mustard oilseed. In particular, the field was characterized by a broad range of salinity and boron values. Salinity and boron are properties common to marginally productive soils in the SJV that are known to significantly influence Ida Gold mustard oilseed yield. The soil at the study site is a Britto clay loam. The soil taxonomic class is fine, smectitic, thermic Typic Natraqualfs. The Britto series consists of deep, very poorly drained soils with high concentrations of salt and alkali in the lower horizons. The soil ranges from moderately saline to strongly saline (8-16 dS m -1 ). The texture in the top 55 cm is a clay loam and 0.55-1.55 m is a sandy clay loam. The parent material is alluvium derived from sedimentary rock. The mean annual precipitation is 25 cm.

Preliminary and Intensive Apparent Soil Electrical Conductivity Surveys
Preliminary and intensive EC a surveys were conducted on 14 January and 28 January 2014, following a pre-plant irrigation to bring the water content in the root zone to field capacity. The methods and materials used in the EC a surveys followed the protocols and guidelines outlined in Corwin and Lesch [8][9][10]. An EM38 Dual Dipole electrical conductivity meter (Geonics Ltd., Mississaugua, Ontario, Canada. Product identification is provided for the benefit of the reader and does not imply endorsement by USDA.) connected to a GPS and mounted on a non-metallic sled was used in the EC a surveys. Geospatial EC a measurements in the vertical (EM v ) and horizontal (EM h ) coil configurations were taken simultaneously every 3-5 m. Each EC a measurement was geo-referenced using GPS. The GPS receiver accuracy had sub-meter accuracy. The preliminary EC a survey determined whether the study site provided a sufficiently wide range in soil properties influencing Ida Gold mustard oilseed yield to meet the objective of formulating a crop yield model. The preliminary survey consisted of making six east-west traverses. From the geospatial EC a data set, six locations were selected to take core samples (0-1.5 m depth increment), which were immediately analyzed for pH e , saturation percentage (SP), B, and EC e .
Following the preliminary EC a survey, two separate intensive EC a surveys were conducted. One EC a survey for the entire 16.2 ha and another survey confined to the southeastern corner of the field. Figure 3 shows maps of a composite of the EC a survey data for EM v and EM h . The reason for taking the two intensive EC a surveys was because the preliminary cursory EC a survey indicated the greatest variability in salinity over the range that would affect Ida Gold mustard oilseed yield was in the southeast corner of the field; consequently, a separate survey and crop yield/soil sampling were performed for the southeast corner. The intention was to provide a range of soil properties, particularly with regard to salinity and boron, which would influence Ida Gold mustard oilseed yield to varying degrees thereby providing the data to formulate a more robust statistical model of crop yield.

Soil and Ida Gold Mustard Oilseed Yield Sampling Design
Once the two intensive EC a surveys were completed, ESAP software version 2.10R [11][12][13] was used to identify sites where crop yield and soil core samples were taken based on the spatial variation in the EC a survey data. The ESAP software package uses a model-based sampling strategy (i.e., response surface sampling design) to identify sample site locations. The ESAP software identifies sites that characterize the range and variation in the geospatial EC a measurements, reflecting the observed spatial variability in EC a , while minimizing any clustering of the sample sites by maximizing the spatial uniformity of the sampling design across the study area. A detailed discussion of the application of the response surface sampling design using EC a survey data is in Lesch et al. [12].
The sample design for the 16.2-ha field consisted of 20 sample site locations for the entire field and 20 sample site locations for the southeast corner as selected by ESAP ( Figure 3). Soil cores were taken at the 40 sites with a Giddings rig at six depth increments: 0-0.15, 0.15-0.30, 0.30-0.60, 0.60-0.90, 0.90-1.20, and 1.20-1.50 m. Duplicate soil samples were taken at eight sample site locations within 1 m of the original core to establish local-scale variability as explained in Corwin and Scudiero [14]. All soil samples were bagged in zip-lock bags and stored in an ice chest until refrigerated. A total of 288 soil samples were taken (6 depths at each site, 40 sites, and 8 duplicate sites). The depth to the water table was recorded as <1.5 m or >1.5 m.
At each of the 40 sample-site locations biomass and Ida Gold mustard oilseed yield were determined by hand within a 1 m 2 area where each soil sample location was the centroid of the 1 m 2 plant sample area. The biomass and oilseed yield were collected on 28-29 May 2014. Six of the 40 sample site locations had no Ida Gold mustard oilseed yield. All subsequent referral to yield is with respect to oilseed yield.

Soil Chemical and Physical Analyses
In the field, a subsample (100-300 g) of each soil core sample was taken for soil moisture determination. The subsamples were weighed in the field to minimize error due to moisture evaporation. The subsamples were subsequently oven-dried at 110ºC for 24 h and weighed again to determine θg. Saturation pastes were prepared for all 288 soil samples and saturation paste extracts were obtained following the procedure of Rhoades [15]. The saturation extracts were analyzed for the following properties: ECe, SP, pHe, 5 major anions (Cl − , HCO3 − , PO4 3− , NO3 − , SO4 2− ), 4 major cations (Na + , K + , Ca 2+ , Mg 2+ ), B, and sodium adsorption ratio (SAR). The chemical analysis procedures followed were those found in Sparks [16]. Leaching fraction (LF), defined as the ratio of the quantity of water draining past the root zone to that infiltrated into the soil's surface, was estimated using two techniques: (1) the ratio of the EMh ECa divided by EMv ECa and (2) the ratio of Cl concentration in the irrigation water and Cl concentration at 1.2-1.5 m. The LF reflects the excess water applied to translocate salts from the root zone. Each property selected for analysis had the potential to influence Ida Gold mustard oilseed yield.

Soil Chemical and Physical Analyses
In the field, a subsample (100-300 g) of each soil core sample was taken for soil moisture determination. The subsamples were weighed in the field to minimize error due to moisture evaporation. The subsamples were subsequently oven-dried at 110ºC for 24 h and weighed again to determine θ g . Saturation pastes were prepared for all 288 soil samples and saturation paste extracts were obtained following the procedure of Rhoades [15]. The saturation extracts were analyzed for the following properties: EC e , SP, pH e , 5 major anions (Cl − , HCO 3 − , PO 4 3− , NO 3 − , SO 4 2− ), 4 major cations (Na + , K + , Ca 2+ , Mg 2+ ), B, and sodium adsorption ratio (SAR). The chemical analysis procedures followed were those found in Sparks [16]. Leaching fraction (LF), defined as the ratio of the quantity of water draining past the root zone to that infiltrated into the soil's surface, was estimated using two techniques: (1) the ratio of the EM h EC a divided by EM v EC a and (2) the ratio of Cl concentration in the irrigation water and Cl concentration at 1.2-1.5 m. The LF reflects the excess water applied to translocate salts from the root zone. Each property selected for analysis had the potential to influence Ida Gold mustard oilseed yield.

Statistical Analysis and Ida Gold Oilseed Yield Model Formulations
Simple correlations were determined between yield and the edaphic properties of θ g , EC e , SP, pH e , Cl − , HCO 3 − , PO 4 3− , NO 3 − , SO 4 2− , Na + , K + , Ca 2+ , Mg 2+ , B, SAR, and LF. Correlations between the edaphic properties and EC a and between oilseed yield and EC a were determined. Scatter plots of yield vs. individual soil-related properties were also obtained. The correlation analyses and scatter plots served as a basis for the development of the yield model for Ida Gold mustard oilseed. The correlations were useful to determine what properties were likely significant in influencing oilseed yield, while the scatter plots helped to determine the general form of the oilseed yield model. The oilseed yield model was developed using (spatial) multiple regression techniques [17]. The edaphic properties were the regressor or independent variables, and yield was the response or dependent variable. Backward variable selection was used to screen out the clearly nonsignificant edaphic properties with t-score values below 1.8. This predictor screening helped to filter out any multicollinearity in the regressor variables. Statistical data analyses were performed on the individual depth increment data (i.e., 0-0.15, 0.15-0.3, 0.3-0.6, 0.6-0.9, 0.9-1.2, and 1.2-1.5 m) and composite depth increment data (i.e., 0-0.15, 0-0.3, 0-0.6, 0-0.9, 0-1.2, and 0-1.5 m). Multiple regression modeling was performed on the composite depth increment data, and the increment characterized by the best goodness-of-fit was retained for further analyses.
In some instances, sufficient input data for the Ida Gold mustard oilseed yield model did not exist or fell outside the range of data that were used to develop the oilseed yield model. In those instances, the two-piece linear salt tolerance model of Maas and Hoffman [18] was used. For soil salinities exceeding the threshold salinity level, relative yield is estimated by the following equation: where Y r is the relative crop yield, a is the salinity threshold (dS m −1 ), b is the slope expressed in yield decrement percentage per dS m −1 , and EC e is the mean electrical conductivity of the saturation extract for the root zone (dS m −1 ). Maas and Hoffman [18] proposed that crop salt tolerance was represented by two linear lines, one a tolerance plateau with a slope of zero and the other shown in Equation (1) as a salinity dependent line whose slope was the yield reduction per unit increase in salinity. The point where both lines intersect is the salinity threshold, which represents the maximum soil salinity that does not reduce yield. The parameters a and b were determined from a compilation of salt tolerance data available for Ida Gold mustard oilseed for the SJV, including the data collected within this study and the work of Maas [19] and Grieve et al. [20]. If salinity was not limiting yield, then the B level determined the oilseed yield when B data were available. The three-piece trace element tolerance model presented in Page et al. [21] and first suggested by Burton et al. [22] was used. The salt and B tolerance data to develop these models were obtained from the 40 soil core and oilseed yield sample sites identified from EC a -directed soil sampling and from 10 supplemental sites. The locations of the supplemental sites were from a transect covering a range of yields. Salt tolerance data were from those sites varying in oilseed yield where all soil properties were optimal except salinity, which varied over a wide range, while B tolerance data were from those sites where all soil properties were optimal except B, which varied over a wide range. If no salinity and B data were available, then all properties were considered optimal for yield.

Identification of Salt-Affected Soils for Oilseed Production in the SJV
To assess the potential production of Ida Gold oilseed in the SJV, the yield model was applied over all salt-affected soils in the valley. Bohn et al. [23] defines salt-affected soils as soils with a root-zone salinity (i.e., EC e ) above 4 dS m −1 . Above 4 dS m −1 very sensitive, sensitive, and moderately salt-sensitive crops will show yield decrements. From 4-8 dS m −1 the yields of many crops are restricted, from 8-16 dS m −1 only salt tolerant crops yield satisfactorily, and above 16 dS m −1 only a few very tolerant crops produce satisfactory yields [24].
One means of identifying salt-affected soils in the SJV is the use of SSURGO. However, the accuracy and reliability of soil salinity in SSURGO is dubious because salinity is a spatially and temporally variable soil property influenced by crop and irrigation management strategies. Recent NRCS reports (e.g., [25]), provided salinity estimations only for non-irrigated farmland because the influence of irrigation on soil salinity cannot be accounted for using traditional soil survey protocols. Therefore, an evaluation of SSURGO's accuracy with respect to soil salinity is needed. To evaluate the accuracy of SSURGO, salinity assessment surveys were performed on 22 agricultural fields (total area: 542 ha) scattered throughout the WSJV. In 2013, intensive EC a surveys conducted at the 22 fields, collected 41,779 EC a readings at an average density of 175 readings ha −1 . Simultaneous EM h and EM v measurements of EC a were taken. Across the 22 fields, 267 soil-sampling locations were identified using ESAP. Soil cores were taken to a depth of 1.2 m, representing the root-zone depth. Details of the EC a survey and soil sampling are in Scudiero et al. [26].
Soil samples were analyzed for salinity (EC e ; dS m −1 ), gravimetric water content (θ g ; g g −1 ), pH e , and saturation percentage (SP) using procedures presented in Methods of Soil Analysis Part 3 [16] and Part 4 [27]. The SSURGO salinity data for the 22 fields were compared to both root-zone EC e hard and EC a soft data. The hard data (i.e., laboratory measurements of EC e of the 0-1.2 m soil samples) were averaged for each field. The soft data (i.e., geospatial measurements of EC a ) were calibrated to the EC e data using spatial linear regression models [17] with an overall R 2 = 0.93 by Scudiero et al. [26].
A map of EC e using the soft data and EC e -EC a calibration prepared for each field established spatial patterns for comparison to SSURGO map units. If the SSURGO database proved unreliable, then the regional-scale salinity assessment approach developed by Lobell et al. [28], which combines EC a -directed soil sampling with satellite imagery, would be used. This work has already been completed and published by Scudiero et al. [29] for the WSJV.

Development of a Spatial Database of Edaphic Properties Influencing Ida Gold Mustard Yield for the SJV
The most extensive spatial database of edaphic properties in the USA is the SSURGO Database. Collection of the information in SSURGO occurred by walking over the land and observing variation in the soil and vegetation to delineate map units. Characterization of soil properties within map units occurred by the collection of numerous soil samples and their analysis in the laboratory. Occasionally observation trenches are dug to characterize the horizonation. Soil maps outline areas called map units, which describe soils and other components that have unique properties, interpretations, and productivity. Collection of the information occurred at scales ranging from 1:12,000 to 1:63,360. The soils maps are intended for natural resource planning and management. Examples of information available from the database include available water capacity, texture, pH, electrical conductivity, and frequency of flooding; yields for cropland, woodland, rangeland, and pastureland; and limitations affecting recreational development, building site development, and other engineering uses. SSURGO map data can be viewed in the Web Soil Survey or downloaded in ESRI ® Shapefile format. Attribute data can be downloaded in text format. For the marginal soils of the SJV, information for water capacity, texture, pH, and frequency of flooding; yields for cropland, woodland, rangeland, and pastureland; and limitations affecting recreational development, building development, and other engineering uses was obtained through SSURGO.
To supplement the SSURGO data an extensive spatial database of quantitative soils information (i.e., EC e , pH e , saturation percentage, B, available water content, LF) for the SJV exists. The supplemental data set, collected over a period of 25 years by Corwin and colleagues, is a compilation of data that appeared in publications by Bourgault et al. [30], Corwin [31], Corwin and Lesch [8][9][10]32], Corwin et al. [7,[33][34][35][36], Lesch and Corwin [17], Lesch et al. [37][38][39], Loague et al. [40], Rhoades et al. [41], Sanden et al. [42], and Scudiero et al. [26,29,43]. The supplemental data set consisted of edaphic property data from 83 fields within the SJV ranging in spatial extent from 0.4 to 65 ha with from 6 to 72 sample sites within a field. Soil samples were collected at 0.3-m increments to a minimum depth of 1.2 m and occasionally to 1.5 and 1.8 m. The supplemental data were used to determine frequency distributions (i.e., histograms), averages, ranges, and standard deviations for those properties found to influence Ida Gold mustard oilseed yield in the SJV. From this statistical information probability density functions (PDFs) were developed for LF, θ g , and B for the composite depth increment of 0-1.2 m, which was found to be the root zone depth for Ida Gold mustard oilseed. In general, SSURGO provides water content and B ranges associated with soil type. The PDFs, defined within the ranges of water content and B provided by SSURGO, were used as input for Monte Carlo simulations with the oilseed yield model. In instances where ranges of water content or B were not given in SSURGO, then B was assumed optimal and water content was estimated from a pedotransfer function using SSURGO texture data.
The use of degraded soil is crucial to driving down the cost of biofuel production in the SJV. Subsequently, the reliability of spatial salinity data for the SJV was of paramount importance for identifying salt-affected soils and for model input data. There were concerns regarding the reliability of the salinity ranges provided in SSURGO for the root zone due to anthropogenic influences (e.g., leaching of salts due to irrigation); consequently, as discussed in detail in Section 2.2 an evaluation of the reliability of SSURGO root-zone salinity was conducted by comparison to salinity ground-truth measurements of 22 fields in the WSJV presented by Scudiero et al. [26]. If SSURGO salinity in the root zone proved unreliable, then salinity predictions from the regional-scale soil salinity model of Scudiero et al. [29] were used. Scudiero et al. [29] found that as salinity increased, the prediction error increased, and quantified the error within salinity categories (i.e., 0-2, 2-4, 4-8, 8-16, >16 dS m −1 ). To incorporate this uncertainty into the Monte Carlo simulations, PDFs were established for each salinity category. The PDFs were defined by the average residuals and standard deviation of the residuals between observed salinities from the ground-truth salinity measurements of Scudiero et al. [26] and predicted salinities from the regional-scale soil salinity model of Scudiero et al. [29].
Leaching fraction is a difficult edaphic property to obtain and is not found in SSURGO. The supplemental data set of Corwin and colleagues was used to determine the frequency distribution, average, range, and standard deviation of LF for the SJV. Leaching fraction was determined from a ratio of EM h EC a to EM v EC a and from a ratio of Cl in irrigation water to Cl concentration below the root zone). Only those LFs where EM and Cl ratios agreed to within 5% were used.

Feasibility of Biofuel Production for the SJV: Application of the Yield Model and Monte Carlo Simulations
Monte Carlo simulations with the Ida Gold mustard oilseed yield model were performed with 10,000 iterations to provide a range of potential yields (kg ha −1 ) and probability of those yields, which are easily converted to L ha −1 of biofuel. This was done first for the WSJV and then for the entire SJV. The mean, median, standard deviation, skewness, and kurtosis of the Monte Carlo simulation distribution were calculated to characterize quantitatively the PDF and subsequently derive the cumulative density function (CDF) of biofuel production. Once the CDF is known, the probability (and thereby feasibility) of oilseed production in the SJV to support sufficiently a conversion facility becomes evident.
Sufficient input data was not always available for each field in the WSJV or SJV or sometimes the input data was outside the range of data used to develop the crop yield model for Ida Gold mustard oilseed. For these instances, an alternative model was needed (see Figure 1). When a complete set of input data was available at a field location, then the full crop yield model was used. For instances where input data was available but was outside the range of data from which the crop yield model was formulated, then either the EC e (i.e., Equation (1)) or B tolerance model was used, whichever was more limiting at the site. For instances where insufficient input data existed and only EC e and B input data were available, then either the EC e (i.e., Equation (1)) or B tolerance model was used, whichever was more limiting at the site. If insufficient input data existed and only EC e input data was available, then the EC e tolerance model (Equation (1)) was used.

Results and Discussion
The preliminary EC a survey revealed that the greatest variation in salinity over the range that would result in oilseed yield decrements was in the southeast corner of the study site; consequently, two intensive EC a surveys were conducted with separate soil sampling designs for each survey. One EC a survey and associated soil sampling covered the full field and the other focused on the SE corner. Figure 3 shows maps of the horizontal coil configuration (EM h ) and vertical coil configuration (EM v ) EC a surveys. The combined soil sampling designs (i.e., full field and southeast corner) provided a full range of soil properties and oilseed yields to build a robust Ida Gold mustard oilseed yield model. The only property potentially influencing Ida Gold mustard oilseed yield that did not vary significantly at the study site was texture. In general, the fine-textured soils (mainly loams and clay loams), which predominate the WSJV, do not vary to a significant extent because the soil is a consequence of lacustrine deposits and of fine-grained alluvium material originating from the Coastal Ranges [44,45].
Geospatial EC a measurements, both EM h and EM v , were higher on the west side of the 16.2-ha field than on the east side (see Figure 3). The lowest EC a measurements were in the southeast corner, where EC a ranged from 0.15 to 0.76 dS m −1 . Over this range of EC a the yield was found to vary the greatest (see Figure 4) and therefore provided the most useful data for oilseed yield model formulation.
The range of EC a over the entire 16.2 ha was 0.15 to 3.97 dS m −1 . Oilseed yield over the EC a range of 0.76 to 3.97 dS m −1 tended to be low. Yield significantly correlated to both EM h and EM v EC a , with correlation coefficients of 0.68 and 0.51, respectively. The higher correlation of EM h EC a to oilseed yield suggests that the root zone for Ida Gold mustard was around 1 m since the EM h measurement penetrates to a depth of approximately 1 m, while EM v penetrates to approximately 1.5 m. The fact that EC a and oilseed yield correlated indicates that EC a must be measuring a soil property or properties that influence oilseed yield; therefore, the response surface sampling design will successfully map the property or properties [9]. oilseed yield model. The only property potentially influencing Ida Gold mustard oilseed yield that did not vary significantly at the study site was texture. In general, the fine-textured soils (mainly loams and clay loams), which predominate the WSJV, do not vary to a significant extent because the soil is a consequence of lacustrine deposits and of fine-grained alluvium material originating from the Coastal Ranges [44,45]. Geospatial ECa measurements, both EMh and EMv, were higher on the west side of the 16.2-ha field than on the east side (see Figure 3). The lowest ECa measurements were in the southeast corner, where ECa ranged from 0.15 to 0.76 dS m −1 . Over this range of ECa the yield was found to vary the greatest (see Figure 4) and therefore provided the most useful data for oilseed yield model formulation. The range of ECa over the entire 16.2 ha was 0.15 to 3.97 dS m −1 . Oilseed yield over the ECa range of 0.76 to 3.97 dS m −1 tended to be low. Yield significantly correlated to both EMh and EMv ECa, with correlation coefficients of 0.68 and 0.51, respectively. The higher correlation of EMh ECa to oilseed yield suggests that the root zone for Ida Gold mustard was around 1 m since the EMh measurement penetrates to a depth of approximately 1 m, while EMv penetrates to approximately 1.5 m. The fact that ECa and oilseed yield correlated indicates that ECa must be measuring a soil property or properties that influence oilseed yield; therefore, the response surface sampling design will successfully map the property or properties [9]. The measured edaphic properties that were felt to potentially influence Ida Gold mustard oilseed yield included: ECe, θg, SP, pHe, trace elements (B, Se, As, Mo), major cations (Na + , K + , Ca 2+ , Mg 2+ ), major anions (Cl − , HCO3 − , PO4 3− , NO3 − , SO4 2− ), SAR, micro-elevation, depth to groundwater, and groundwater EC. Table 1 is a summary by depth of these edaphic properties, except for micro-elevation, depth to groundwater, and groundwater EC. The measured edaphic properties that were felt to potentially influence Ida Gold mustard oilseed yield included: EC e , θ g , SP, pH e , trace elements (B, Se, As, Mo), major cations (Na + , K + , Ca 2+ , Mg 2+ ), major anions (Cl − , HCO 3 − , PO 4 3− , NO 3 − , SO 4 2− ), SAR, micro-elevation, depth to groundwater, and groundwater EC. Table 1 is a summary by depth of these edaphic properties, except for micro-elevation, depth to groundwater, and groundwater EC.  = gravimetric water content, SP = saturation percentage, EC e = electrical conductivity of the saturation extract, SAR = sodium adsorption ratio, LF = leaching fraction. ‡ Significant. Skewness is significant if skewness divided by standard error of skewness > 2. Kurtosis is significant if kurtosis divided by stand error of kurtosis > 2. § Leaching fraction was determined by dividing the Clconcentration of the irrigation water by the Cl − concentration of the saturation extract at the 1.2-1.5 m depth increment at each site where the oilseed yield was greater than zero. £ Sites where oilseed yield was greater than zero. These sites were used for Ida Gold mustard oilseed yield model development. Table 1 reveals patterns in the field-wide soil profile. Field-wide average soil salinity (EC e ) increases with depth up to the 0.3-0.6, levels off at the 0.3-0.6 and 0.6-0.9 m depth increments, and then decreases with depth. Saturation percentage (SP) is reasonably constant over depth ranging from means of 48.13% at 1.2-1.5 m to 54.85% at 0.3-0.6 m, indicating uniform texture through the soil profile. Gravimetric water content (θ g ) at field capacity increases with depth from 0.10 kg kg −1 at 0-0.15 m to 0.24 kg kg −1 at 1.2-1.5 m. Boron levels tend to be lower below 0.3-0.6 m. pH increases with depth from 7.29 at 0-0.15 m to 7.98 at 1.2-1.5 m. SAR increases with depth from 6.53 at 0-0.15 m to 19.04 at 0.3-0.6 m, then decreases to 10.72 at 1.2-1.5 m. Table 2a presents mean and range statistics, standard deviation, standard error, coefficient of variation, skewness, and kurtosis for the composite 0-1.5 m depth. The highest coefficients of variation (CVs) are for EC e , various cations and anions (e.g., Na + , Ca 2+ , SO 4 2− , PO 4 − , and Cl − ), and SAR, while the lowest CVs are for pH e , SP, and θ g . All edaphic properties in Table 2a

Ida Gold Mustard Oilseed Yield Models for Marginal SJV Soils
Exploratory statistical analyses revealed that Ida Gold mustard oilseed yield was most significantly correlated to individual edaphic properties for the top 1.2 m. Table 2b presents mean and range statistics, standard deviation, standard error, coefficient of variation, skewness, and kurtosis for the composite 0-1.2 m depth. Table 3 presents simple correlations between edaphic properties and both EC a and oilseed yield for the 0-1.2 m depth increment. The edaphic properties most significantly correlated to EC a include θ g , EC e , B, SAR, LF (determined by the ratio EM h EC a /EMv EC a ), and SP. The edaphic properties most significantly correlated to oilseed yield include θ g , EC e , B, LF, and SP. Corwin and Lesch [9] indicate that the depth increment or composite depth increment associated with the best-fitting yield model (i.e., highest R 2 ) reflects the root zone of the crop. The top 1.2 m resulted in the most statistically significant and best-fit Ida Gold mustard oilseed yield model (see Equation (3)). Consequently, the 0-1.2 m soil interval was taken to represent the root zone of Ida Gold mustard at the study site. Subsequently, all data presented and discussed are with respect to the 0-1.2 m composite depth or to individual depth increments that lie within the composite depth of 0-1.2 m.
Exploratory statistical analyses revealed that θ g , EC e , B, and LF were the most influential edaphic properties on oilseed yield. Scatter plots of these properties vs oilseed yield indicated quadratic relationships of EC e and B to oilseed yield and linear relationships of LF and θ g to yield. Based on the initial exploratory correlation and multiple linear regression analysis, the following regression model structure was proposed to describe edaphic property effects on oilseed yield: where Y is the Ida Gold mustard oilseed yield (kg ha −1 ); B is boron concentration (mg L −1 ); EC e is electrical conductivity of the saturation extract (dS m −1 ); LF is the leaching fraction; θ g is the gravimetric water content (kg kg −1 ); β 0 , β 1 , β 2 , . . . , β 6 are the regression model parameters, and ε is the random error component, initially assumed to be normally distributed and spatially independent. Ordinary least squares (OLS) regression techniques resulted in a fitted regression equation with a R 2 = 0.89 and adjusted R 2 = 0.78. Adjusting for spatial autocorrelation using the maximum-likelihood approach resulted in the following Ida Gold mustard oilseed yield model (Equation (3) Equation (3) represents the most parsimonious and robust model for marginally productive salt-affected soils of the WSJV. Any locations where no oilseed yield was obtained were not used in the model development. The LF and θ g parameters are highly significant at or near the 0.01 level, and the EC e (linear and quadratic) and B (linear and quadratic) parameter estimates are significant at or near the 0.05 level. The LF and θ g parameters are both positive, implying that the yield increased as either LF or θ g increased, which is physically sound since increased leaching reduces osmotic stress and increased water content increases the plant-available water reducing matric stress. The positive linear and negative quadratic EC e terms imply that the yield increased at low EC e up to a point of maximum yield with respect to EC e and decreased beyond the maximum. The point of maximum yield with respect to EC e was calculated by setting the first partial derivative of the fitted regression to zero with respect to EC e , which resulted in a value of 6.8 dS m −1 . Similarly, the positive linear and negative quadratic B terms imply that the yield increased under low B up to a point of maximum yield with respect to B and decreased beyond the maximum. The point of maximum yield with respect to B was 4 mg L −1 .
In those instances where sufficient input data for Equation (3) did not exist or fell outside the range of data that was used to develop the model, the two-piece linear salt tolerance model (i.e., Equation (1)) of Maas and Hoffman [18] was used to predict oilseed yield. Salt tolerance data at the study site established a salinity threshold of 8.3 dS m −1 , which is the term a in Equation (1), and a yield decrement slope of 17%, which is the term b in Equation (1) (Figure 5a). The salinity threshold of 8.3 dS m −1 corresponds reasonably well with the salinity of maximum oilseed yield of 6.8 dS m −1 in Equation (3). Equation (1) established the upper limit of the salinity range of salt-affected soils that would grow Ida Gold mustard oilseed. A 17% yield decrement for each 1 dS m −1 increase in root-zone soil salinity beyond 8.3 dS m −1 resulted in no oilseed yield above 14.3 dS m −1 . It is important to note that the salinity range of 4-14.3 dS m −1 was not necessarily economically viable. Once the feasibility of reaching the 115 ML yr −1 goal is established, then the maximum yield decrement that is economically viable could be determined. Viability would take into account other economically relevant factors such as (1) selling the seed meal remaining, after the oilseed has been pressed to extract the oil, as Se-enriched meal for livestock or as a herbicide used in organic agriculture and (2) incorporation of gasification to generate power to run the oil press. The two-piece linear salt tolerance model (i.e., Equation (1)) of Maas and Hoffman [18] was not the best model to fit the data as seen in Figure 5a. A quadratic model (Figure 5b) actually fits the data best: Y = 74.0 + 254.6 EC e − 18.8 EC e 2 (R 2 = 0.87) By taking the derivative of Equation (4) with respect to EC e and setting it equal to 0, EC e corresponds to the peak mustard oilseed yield, which is 6.8 dS m −1 . This is equal to the EC e producing the maximum yield in Equation (3). Similarly, a quadratic relationship between yield and EC e was also found by Corwin et al. [7] for cotton seed yield. The explanation for a quadratic relationship between EC e and yield of cotton seed and mustard oilseed is that when the plant is osmotically stressed it puts energy into the development of the reproductive part of the plant rather than vegetative tissue; consequently, as EC e increases up to the salinity threshold biomass yield of both crops remains constant or may decrease, while mustard oilseed and cotton seed yield steadily increase. After the salinity threshold is reached then oilseed and cotton seed yield and the biomass of both crops decrease as salinity increases. This quadratic relationship indicates that more breeding research is needed to obtain new cultivars of mustard that have their peaks of potential yield at different ranges of salinity to maximize production in the SJV. found by Corwin et al. [7] for cotton seed yield. The explanation for a quadratic relationship between ECe and yield of cotton seed and mustard oilseed is that when the plant is osmotically stressed it puts energy into the development of the reproductive part of the plant rather than vegetative tissue; consequently, as ECe increases up to the salinity threshold biomass yield of both crops remains constant or may decrease, while mustard oilseed and cotton seed yield steadily increase. After the salinity threshold is reached then oilseed and cotton seed yield and the biomass of both crops decrease as salinity increases. This quadratic relationship indicates that more breeding research is needed to obtain new cultivars of mustard that have their peaks of potential yield at different ranges of salinity to maximize production in the SJV.
By taking the derivative of Equation (5) with respect to B and setting it equal to 0, B corresponds to the peak oilseed yield, which is 6.5 mg L −1 .
The piece-wise linear models of Figures 5a and 6a for salinity and B, respectively, are traditional models found throughout the literature. The piece-wise linear models were used in the Monte Carlo simulations. However, the quadratic models provided slightly better fits to the data; consequently, a second set of Monte Carlo simulations was performed with the salinity and B quadratic models, Equations (4) and (5), respectively.  (Figure 6a). Below 4.2 mg L −1 of B and above 8 mg L −1 the yield of Ida Gold mustard oilseed drops, where B < 4.2 mg L −1 is a B deficiency and B > 8 mg L −1 is a B toxicity. Above 8 mg L −1 the drop was 28% for every 1 mg L −1 increase in B with no yield occurring above 11.6 mg L −1 of B. A quadratic model (Equation (5)) produced a comparable fit to the B tolerance data (Figure 6b): By taking the derivative of Equation (5) with respect to B and setting it equal to 0, B corresponds to the peak oilseed yield, which is 6.5 mg L −1 .
The piece-wise linear models of Figures 5a and 6a for salinity and B, respectively, are traditional models found throughout the literature. The piece-wise linear models were used in the Monte Carlo simulations. However, the quadratic models provided slightly better fits to the data; consequently, a second set of Monte Carlo simulations was performed with the salinity and B quadratic models, Equations (4) and (5), respectively.

Evaluation of SSURGO Soil Salinity Accuracy and Identification of Salt-Affected Soils
An evaluation of the accuracy of SSURGO spatial data for root-zone soil salinity revealed that only 5 out of 22 fields assessed the mean salinity or range in salinity accurately, suggesting that the more transient salinity levels and patterns in the root zone are not captured in the one-time measurements of NRCS soil surveys. However, SSURGO was able to assess 15 out of 22 fields accurately for salinity below the root zone, indicating that the salt levels below the root zone remained relatively unchanged and unaffected by anthropogenic influences. The failure of SSURGO to provide accurate root-zone soil salinity spatial data necessitated a reliance on the regional-scale salinity model developed by Scudiero et al. [29] as input data for the Ida Gold mustard oilseed yield model.
Because of the inaccuracy of salinity in SSURGO the Scudiero et al. [29] regional-scale salinity model was used to identify salt-affected soils (EC > 4 dS m −1 ) for the SJV. Figure 7 shows the extent of salt-affected soils for the WSJV as estimated by the Scudiero et al. [29] regional-scale salinity model. A comparison of total land cover greater than 4 dS m −1 for SSURGO and the Scudiero et al. [29] regional-scale salinity model indicates that SSURGO estimated 33% less salt-affected land for the WSJV. Furthermore, the distribution of salt-affected soils was concentrated in contiguous patterns along the eastern half of the WSJV for SSURGO, whereas the Scudiero et al. [29] salinity patterns were more diffusely spread as shown in Figure 7. Only fields identified as salt-affected were subsequently used in the Monte Carlo simulations. Any field where the field average root zone ECe was estimated to be <4 dS m −1 was disregarded and not included in the Monte Carlo simulations.

Evaluation of SSURGO Soil Salinity Accuracy and Identification of Salt-Affected Soils
An evaluation of the accuracy of SSURGO spatial data for root-zone soil salinity revealed that only 5 out of 22 fields assessed the mean salinity or range in salinity accurately, suggesting that the more transient salinity levels and patterns in the root zone are not captured in the one-time measurements of NRCS soil surveys. However, SSURGO was able to assess 15 out of 22 fields accurately for salinity below the root zone, indicating that the salt levels below the root zone remained relatively unchanged and unaffected by anthropogenic influences. The failure of SSURGO to provide accurate root-zone soil salinity spatial data necessitated a reliance on the regional-scale salinity model developed by Scudiero et al. [29] as input data for the Ida Gold mustard oilseed yield model.
Because of the inaccuracy of salinity in SSURGO the Scudiero et al. [29] regional-scale salinity model was used to identify salt-affected soils (EC > 4 dS m −1 ) for the SJV. Figure 7 shows the extent of salt-affected soils for the WSJV as estimated by the Scudiero et al. [29] regional-scale salinity model. A comparison of total land cover greater than 4 dS m −1 for SSURGO and the Scudiero et al. [29] regional-scale salinity model indicates that SSURGO estimated 33% less salt-affected land for the WSJV. Furthermore, the distribution of salt-affected soils was concentrated in contiguous patterns along the eastern half of the WSJV for SSURGO, whereas the Scudiero et al. [29] salinity patterns were more diffusely spread as shown in Figure 7. Only fields identified as salt-affected were subsequently used in the Monte Carlo simulations. Any field where the field average root zone EC e was estimated to be <4 dS m −1 was disregarded and not included in the Monte Carlo simulations.

Input Data for Ida Gold Mustard Oilseed Yield Model
The reliability and accuracy of the spatial data that serves as input into any model are just as critical as the model itself, exemplified by the old adage, garbage in garbage out. Sensitivity analysis of the crop yield model provides an indication of the input variables that need the greatest level of accuracy and therefore need particular scrutiny when building a spatial database of inputs for the crop yield model. Sensitivity analysis (Table 4) established the degree of influence that each edaphic property had on oilseed yield in Equation (3). The influence was determined by calculating how much the predicted oilseed yield changed when the value for each independent variable in Equation (3) was individually shifted by 1 standard deviation from its mean level or point of maximum yield with respect to the edaphic property. The means and standard deviations were obtained for the 0-1.2 m depth increment excluding any points where no yield occurred (Table 2b). A baseline yield was used as the point of reference to establish the percentage of change. A baseline value of 6.8 dS m −1 was used for salinity, rather than the mean ECe level of 10.0 dS m −1 from Table 2b because the value 6.8 dS m −1 represents the point of maximum yield with respect to the quadratic salinity

Input Data for Ida Gold Mustard Oilseed Yield Model
The reliability and accuracy of the spatial data that serves as input into any model are just as critical as the model itself, exemplified by the old adage, garbage in garbage out. Sensitivity analysis of the crop yield model provides an indication of the input variables that need the greatest level of accuracy and therefore need particular scrutiny when building a spatial database of inputs for the crop yield model. Sensitivity analysis (Table 4) established the degree of influence that each edaphic property had on oilseed yield in Equation (3). The influence was determined by calculating how much the predicted oilseed yield changed when the value for each independent variable in Equation (3) was individually shifted by 1 standard deviation from its mean level or point of maximum yield with respect to the edaphic property. The means and standard deviations were obtained for the 0-1.2 m depth increment excluding any points where no yield occurred (Table 2b). A baseline yield was used as the point of reference to establish the percentage of change. A baseline value of 6.8 dS m −1 was used for salinity, rather than the mean EC e level of 10.0 dS m −1 from Table 2b because the value 6.8 dS m −1 represents the point of maximum yield with respect to the quadratic salinity response pattern. For the same reason a B baseline of 4.0 mg L −1 was used in the sensitivity analysis. The calculated percentage yield change shown in Table 4 indicates that B is the most significant factor influencing Ida Gold mustard oilseed yield, followed by EC e , then LF, and finally θ g . In the case of EC e , oilseed yield model input values were obtained from the 30 x 30 m predictions from Scudiero et al. [29] and used to determine the average EC e for each field for the root zone (0-1.2 m) within the SJV. For the Monte Carlo simulations, PDFs were defined by the average residual and standard deviation of the residuals for each category of salinity (i.e., 0-1, 1-2, 2-3, 3-4, 4-5, 5-6, 6-7, 7-8, 8-9, 9-10, 10-11, 11-12, 12-13, 13-14, 14-15, 15-16, and >16 dS m −1 ) using the predicted and ground-truth EC e data from Scudiero et al. [5,29]. The input EC e of each pixel in the Monte Carlo simulations was determined from the predicted EC e and residual PDF. Table 5 is a summary of the average residuals, standard deviation of the residuals, and data count for each salinity category and for the entire data set. All PDFs were normally distributed. The PDFs for categories 0-1, 1-2, 2-3, and 3-4 d S m −1 were not used since only those fields with an average root zone EC e of greater than 4 dS m were used in the Monte Carlo simulations. There was not a substantial difference in the Monte Carlo simulation findings using the PDF for the entire data set in Table 5 as compared to the individual salinity categories; consequently, all subsequent Monte Carlo simulation discussions will be for the use of the PDF for the entire EC e residual data set. Table 5. Summary of the average residual, standard deviation of the residual, and data count for each category of predicted salinity (i.e., EC e ) from Scudiero et al. [29], indicating the uncertainty of the EC e prediction and used in defining the EC e residual PDFs. Definition: EC e = electrical conductivity of the saturation extract.
Oilseed yield must be converted to biofuel yield. The oil content of Ida Gold mustard oilseed was 26% and the extraction efficiency was 64%. Subsequently, for every 1000 kg ha −1 of oilseed produced then 166.4 kg ha −1 of 100% biofuel resulted, which represented 175.3 L ha −1 of biofuel. It is of importance to note that transesterified biofuel is generally blended with diesel (1:5 ratio of biofuel to diesel) and sold as 20% biofuel in California.

Monte Carlo Simulations Feasibility of Biofuel Production for the SJV
Monte Carlo simulations for the WSJV never produced a simulation that resulted in more than 71.9 ML yr −1 of biofuel produced, which is well below the goal of 115 ML yr −1 . Subsequently, Monte Carlo simulations were performed for the entire SJV. To do this, the regional-scale salinity model of Scudiero et al. [29] was applied to the entire SJV to identify salt-affected soils, which totaled 9.7 × 10 5 ha. Monte Carlo simulations resulted in the PDF and associated CDF shown in Figure 8a,b, respectively. The PDF and CDF are based on Monte Carlo simulations that incorporate the piece-wise models of salinity and B tolerance. The histogram of Monte Carlo simulations is best fit with the shifted gamma PDF in Equation (6) where Q is the biofuel production in ML yr −1 . A comparison of the means, medians, standard deviations, skewness, and kurtosis for the measured Monte Carlo simulation data and for the estimates from Equation (6) (Figure 8b) there is a 17% probability of meeting the minimum production level when all salt-affected soils in the SJV are utilized for oilseed production. When the Monte Carlo simulations include quadratic salt and B tolerance models (i.e., Equations (5) and (6), respectively), there is little difference, with a 15% probability. Several potential weaknesses in the approach need discussion as well as how the impacts of these weaknesses were mitigated in the Monte Carlo simulations. First, the oilseed yield model was developed from a single field with a fine-texture soil that did not vary significantly. Even though the study site may have been representative of many of the fine-textured soils of the WSJV, it was not representative of coarse-textured soils found in the east side of the SJV; consequently, the use of the oilseed yield model for the entire SJV is dubious. Second, the regional-scale salinity model of Scudiero et al. [29] was developed from a database that did not include tree crops or vineyards, which makes the input salinities for the oilseed yield model dubious for areas containing vineyards or tree crops. Scudiero et al. [5] showed that the regional-scale salinity model of Scudiero et al. [29] over estimates salinity levels for orchards and vineyards, which would render reduced oilseed yield estimates of Ida Gold mustard when planted between tree rows or would identify these lands as salt-affected when they are not. Third, applying the oilseed yield model outside the range of data that was used to develop it is problematic.
To rectify these problems, several precautions were taken. No areas where orchards or vineyards occurred were included in the Monte Carlo simulations. In fields where the texture was coarse, the oilseed yield model of Equation (3) was not applied. Rather, the two-piece linear salt tolerance model of Maas and Hoffman [18] presented in Equation (1) or the three-piece linear B tolerance model presented in Figure 6 were used, depending on which was most limiting to oilseed yield. Similarly, in cases where the input data for the oilseed yield model (Equation (3)) fell outside the range of data used to develop the model, then the two-piece linear salt tolerance model of Maas and Hoffman [18] or the three-piece linear B tolerance model were used instead. Even though the feasibility of meeting desired oilseed production levels is dubious, there are circumstances that could improve the likelihood of meeting the 115 ML lower limit of production. Water in terms of water content and water available for leaching are two impactful properties pertaining to oilseed yield in the model. Maintaining the root zone water content at high levels with adequate leaching could drive up yields sufficiently high to make oilseed a viable biofuel in the SJV. In particular, high frequency irrigation to maintain root zone water contents near field capacity would have significant impact. However, this is only feasible when the SJV has sufficient irrigation water supplies, which may be the exception rather than the rule as shown by the recent 6-year California drought (2011-2016).
Orchards, whether in salt-affected (ECe > 4 dS m −1 ) or non-salt-affected (ECe < 4 dS m −1 ) soils, could play a major role in enhancing the feasibility of biofuel production in the SJV. Planting Even though the feasibility of meeting desired oilseed production levels is dubious, there are circumstances that could improve the likelihood of meeting the 115 ML lower limit of production. Water in terms of water content and water available for leaching are two impactful properties pertaining to oilseed yield in the model. Maintaining the root zone water content at high levels with adequate leaching could drive up yields sufficiently high to make oilseed a viable biofuel in the SJV. In particular, high frequency irrigation to maintain root zone water contents near field capacity would have significant impact. However, this is only feasible when the SJV has sufficient irrigation water supplies, which may be the exception rather than the rule as shown by the recent 6-year California drought (2011-2016).
Orchards, whether in salt-affected (EC e > 4 dS m −1 ) or non-salt-affected (EC e ≤ 4 dS m −1 ) soils, could play a major role in enhancing the feasibility of biofuel production in the SJV. Planting mustard oilseed between rows in orchards, sometimes referred to as intercropping or alley cropping, provides tremendous acreage for crop growth on lands that are otherwise left fallow in many instances. Intercropping has several additional advantages, including weed control, wind and water erosion control, dust management, improved percolation, more effective use of land and resources, and additional crop revenue. There are roughly 130,000 ha of orchards in the SJV. By planting an oilseed crop between rows as a secondary crop, biofuel production could increase roughly 15-30%, which would increase the probability of meeting the 115 ML yr −1 target from 15-17% to 60-80% for the entire SJV. Even though intercropping with oilseed may make biofuel production feasible, it certainly does not necessarily make it economically viable. Furthermore, secondary crops from intercropping are a management challenge for producers filled with additional management problems.
Aside from meeting the 115 ML yr −1 target, economic viability is a concern. Growing oilseed on marginally productive soil is not the only means of lowering cost. Another resource that can lower biofuel production cost is the reuse of degraded water for irrigation. Corwin [31] showed in a long-term (12 years) study that the reuse of 3-5 dS m −1 drainage water was a viable means of reclaiming saline-sodic soil on the WSJV, supporting a salt-tolerant crop during the reclamation process, and providing financial return to the producer on otherwise non-productive farmland. In addition, drainage water is a particularly valuable alternative water resource during drought years in the SJV, provided drainage water is in sufficient supply.
Even though biofuel feasibility evaluation is a valuable application of sensor technology, the unique aspect of this study is that it presents an innovative regional-scale approach for modeling the interaction of edaphic properties on crop yield by characterizing the spatial variability of soil properties influencing yield using proximal and satellite sensors. In essence, each sampling location identified from an EC a survey serves as an independent crop yield study looking at the interaction of edaphic properties on yield, thereby rendering a more robust model with lower long-term labor requirements than plant salt (or B) tolerance studies provided in the past. Furthermore, past plant salt tolerance studies did not evaluate the interaction of edaphic properties on crop yield.
The same approach can be used to predict drought impacts on crop productivity within an agricultural region, to identify reclamation needs to optimize crop productivity, and to manage resources (e.g., irrigation water, land, crop selection) and edaphic properties (e.g., salinity, B, etc.) to optimize crop yield. Because of its regional-scale application and quantitative probability assessment, the approach provides land and water resource specialists with a tool to make regionaland landscape-scale resource decisions with a clear understanding of the likelihood of impact on agricultural productivity.

Conclusions
This study exemplifies the combined use of proximal and satellite sensors to answer a practical agricultural question of national significance with strategic military implications concerning the feasibility of biofuel production in California's San Joaquin Valley to help meet the aviation fuel needs of the U.S. military. A quantitative assessment of the probability of attaining a minimum of 115 ML yr −1 of biofuel from oilseed indicates that there is a 15-17% probability of meeting the minimum production level when all salt-affected soils in the SJV are utilized for oilseed production (Figure 8). The low probability of meeting the minimum production level and the fact that all salt-affected soils throughout the SJV, not just WSJV salt-affected soils, would need to be cropped with Ida Gold mustard to produce sufficient oilseed biofuel to support a 115 ML conversion plants makes the idea untenable. Furthermore, a significant portion of the 115 ML of biofuel would be produced on soils of moderate salinity (i.e., 4-6 dS m −1 ), moderate B levels (i.e., 4-8 mg L −1 ), high LF, and high available water content, which are soils that can either be easily reclaimed to produce higher cash crops or used to grow salt and B tolerant crops of higher cash return than oilseed, such as pistachio.