Agriculture is a major source of NOx pollution in California

See allHide authors and affiliations

Science Advances  31 Jan 2018:
Vol. 4, no. 1, eaao3477
DOI: 10.1126/sciadv.aao3477


Nitrogen oxides (NOx = NO + NO2) are a primary component of air pollution—a leading cause of premature death in humans and biodiversity declines worldwide. Although regulatory policies in California have successfully limited transportation sources of NOx pollution, several of the United States’ worst–air quality districts remain in rural regions of the state. Site-based findings suggest that NOx emissions from California’s agricultural soils could contribute to air quality issues; however, a statewide estimate is hitherto lacking. We show that agricultural soils are a dominant source of NOx pollution in California, with especially high soil NOx emissions from the state’s Central Valley region. We base our conclusion on two independent approaches: (i) a bottom-up spatial model of soil NOx emissions and (ii) top-down airborne observations of atmospheric NOx concentrations over the San Joaquin Valley. These approaches point to a large, overlooked NOx source from cropland soil, which is estimated to increase the NOx budget by 20 to 51%. These estimates are consistent with previous studies of point-scale measurements of NOx emissions from the soil. Our results highlight opportunities to limit NOx emissions from agriculture by investing in management practices that will bring co-benefits to the economy, ecosystems, and human health in rural areas of California.


Nitrogen oxide (NOx = NO + NO2) gases are among the most important components of air pollution, which, according to the World Health Organization, is responsible for one in eight premature deaths worldwide (1). These nitrogen (N) gases have been linked to upper respiratory disease, asthma, cancer, birth defects, cardiovascular disease, and sudden infant death syndrome (2, 3). Global studies have pointed to similarities in the magnitude of NOx emissions from fossil fuel combustion and soil, with the largest soil emissions from regions with heavy N fertilizer applications (47). Despite the significance of soil microbial NOx emissions at the global scale, policies have focused largely on limiting NOx from mobile and stationary fossil fuel sources (8, 9). Where agriculture is an important source of NOx, strategies to reduce nonpoint emissions will need to incorporate soil management approaches and policies that are fundamentally different from fossil fuel sources.

California is considered the world’s sixth largest economy in terms of gross national product and supports 12.2% of the U.S. food economy (10). The state has instituted policies to reduce NOx pollution from fossil fuel sources, resulting in NOx declining by 9% per year in Los Angeles, San Francisco, and Sacramento over the period of 2005 to 2008 (11). Recent findings have suggested that agriculture is one of the dominant sources of NOx in the United States, particularly in the midwest region, where fertilizer inputs are substantial (6, 12). In California, local field measurements have similarly ascribed high NOx emissions to agricultural soil (13). Matson et al. (14) provided some of the first evidence of substantial NOx production from agricultural soils in California’s Central Valley; however, a statewide assessment, which is needed to drive new policies for NOx pollution, is hitherto lacking. The California Air Resources Board (CARB) estimates that ~3.8% of the state’s NOx budget can be attributed to cropland soils, but these estimates are based on data limited to farms located within 200 km of Sacramento and miss many of the most heavily fertilized areas in the state (15). Moreover, CARB does not include these estimated emissions in their official statewide database for air quality modeling (16).

Here, we provide the first large-scale quantification of soil NOx emissions for California through two different approaches: integrative “bottom-up” spatial modeling and “top-down” airborne NOx measurements. This two-pronged approach allows us to independently examine the contribution of biogenic NOx emissions in California while comparing these estimates to local empirical data. Our overarching hypothesis is that biogenic emissions of NOx from agricultural areas are much higher than we used to believe and could be a major source of atmospheric NOx statewide. Alternatively, if agricultural sources are of minor significance, then we would expect to find uniformly low emissions throughout natural and agricultural ecosystems.


Our combined bottom-up and top-down estimates uniformly point to high NOx emissions from California’s agricultural soil, revealing a significant unrecognized source of N pollution statewide. Our bottom-up model reveals that 161,100 metric tons of NOx-N year−1 is emitted from California soil with croplands accounting for 79% of total emissions. When combined with data on existing mobile and stationary fossil fuel sources (16), our results indicate that fertilized croplands account for 20 to 32% of total NOx-N emissions from all sectors of the state, whereas natural soils account for 5 to 9% (Fig. 1). A meta-analysis of soil NOx emissions from the existing literature demonstrates quantitative coherence between our model-based estimates and empirical measurements from different areas of the state (Table 1). Mean NOx emissions from California cropland soils were 19.8 (±27.3 SD) kg of N ha−1 year−1 and ranged from 0 to 276 kg of N ha−1 year−1 (Fig. 2), with 1/4 quartile and 3/4 quartile values of 4.3 and 24.9 kg of N ha−1 year−1, respectively. NOx emissions were largest from agricultural soils where N fertilizer applications can reach >600 kg of N ha−1 year−1 (average N fertilizer rates for fertilized soils, 131.8 kg of N ha−1 year−1; Fig. 3). A spatial maximum hot spot of soil NOx emissions is identified for southern reaches of the state, where climate is relatively hot and arid (17). The model also predicts local maxima in the Sacramento Delta region, the Salinas Valley, and the San Joaquin Valley, with the latter being confirmed by aircraft measurements (see below).

Fig. 1 Contribution of soils to statewide NOx emissions.

Based on CARB emission estimates and IMAGE-modeled emission estimates for cropland and natural ecosystems (A) without vegetation scavenging (gross rates) and (B) with 50% of NOx emissions scavenged by vegetation (net rates).

Table 1 Modeled values and observed values collected from the literature of NO emissions in California.

LAT, latitude; LONG, longitude; SFREC, Sierra Foothill Research and Extension Center.

View this table:
Fig. 2 Estimates of NOx emissions from California soils (natural and cropland) generated by using stable isotopic modeling and IMAGE model.
Fig. 3 Nitrogen fertilizer inputs to California soils.

Fertilizer application rates are generated based on crop type, using crop-specific data provided by the DWR of California and USDA fertilizer consumption database for 1964 to 2006.

Modeled NOx emissions track N fertilizer applications

Our findings support the hypothesis that biogenic sources represent a significant fraction of NOx emissions in California, particularly in areas with high N fertilizer applications. Although we report gross soil emission estimates, NOx uptake by vegetation can cut atmospheric NOx emissions in half (4, 18). We thereby provide a more conservative estimate that attributes 25% of statewide NOx to the soil (Fig. 1B), which assumes that half of the soil NOx is lost to dry deposition within nearby vegetation canopies. Reducing uncertainty regarding the soil contribution of NOx to the statewide budget will require spatial and temporal assessments that can distinguish between sources.

Our findings for California are consistent with previous global-scale estimates given the tremendous agricultural productivity of the state: Yienger and Levy (4) used a model to demonstrate that soils account for 50% of the total NOx budget in remote agricultural regions of the Northern Hemisphere, Jaeglé et al. (6) found that soils were responsible for 22% of worldwide NOx emissions, Wang et al. (7) estimated that agricultural soils were responsible for 20 to 30% of global NOx sources, and Davidson et al. (5) estimated that soils were responsible for 10% of emissions in the Southeastern United States, with agriculture contributing to 75% of total soil NOx emissions. Furthermore, our soil NOx emission estimates are on par with those from agrarian areas in Europe (24 to 62%) (19) and during the crop-growing season in the midwestern United States (15 to 40% in June and July) (12).

In contrast to high mean efflux from agricultural soils (average, 19.8 kg of N ha−1 year−1), NOx emissions from natural ecosystems were much lower (average, 1.0 kg of N ha−1 year−1). This points to the importance of N inputs (in fertilizer) in accelerating NOx emissions from soil microbial communities (20). Our model-based estimates compared favorably with published literature values for California and beyond, which range from 0 to 25 kg of N ha−1 year−1 (Table 1). Likewise, in the southern and midwestern regions of the United States, empirically measured NOx emissions from cultivated soils are on the order of 3 to 14 kg of N ha−1 year−1 versus 0.1 to 3 kg of N ha−1 year−1 from nearby unfertilized grasslands (2125). Other studies have also demonstrated that cultivated land produces NOx emissions an order of magnitude larger than forest soils in the same biome (26).

Emission estimate variability

Modeled soil NOx emissions produced ranges that were similar to those reported in the literature (Table 1). Both modeled and observed (in situ chamber measurements) emissions were typically less than 5 kg of N ha−1 year−1, with larger values in the Imperial Valley (located in Southern California between the Salton Sea and Mexico; Fig. 2). Whereas most observed values were similar to those of the model, NOx emissions showed marked variability, consistent with spatial heterogeneities of soil microbial processes. Biogenic NOx emissions can also vary temporally, with the largest emissions spiking when fertilizers are applied (12). The timing of fertilization in California varies regionally, considering the variety of crop species grown and the different management practices used. Our model validation was restricted to a handful of empirical studies (13, 14, 2730), which demonstrate the need for more ground measurements throughout California to better assess the local impact and spatial distribution of soil NOx emissions.

Surface emissions estimates from airborne NOx observations

We used airborne measurements of NOx concentrations to estimate regional emissions in the San Joaquin Valley. This allows us to further verify our model and determine how NOx emissions from the soil might affect regional compliance with ambient air quality standards, which are based on 1-hour and annual average concentration thresholds (31). Our top-down approach involved repeated airborne measurements of NOx made between Fresno and Visalia during the summer of 2016 (fig. S3) in conjunction with the California Baseline Ozone Transport Study (CABOTS) coordinated by CARB. Careful accounting of the height of the atmospheric boundary, coupled with direct measurements and some judicious estimates of all the terms in the NOx concentration budget equation, allowed us to estimate surface emissions (32) of NOx in the region of the flight experiment (see the Supplementary Materials). The average of six flight days (three at the end of July and three at the beginning of August) over a region of ~720,000 ha yielded a NOx emission estimate of 190 ± 130 metric tons day−1. According to the CARB California Emissions Projection Analysis Model (CEPAM) (33), which includes fossil fuel but not natural sources, the sum of average summertime NOx emissions over all three counties in the surrounding area (Fresno, Tulare, and Kings covering over 3.1 million ha) amounts to 100 metric tons day−1. Although the exact area and diurnal timing of the emissions from the CEPAM inventory cannot be precisely compared to the spatial and temporal footprint of our airborne sampling, the comparison between the CEPAM inventory and airborne sampling shows that soil emissions are likely a very important source of atmospheric NOx, especially in the agriculturally intensive San Joaquin Valley. In this case, the agricultural soil source would need to account for at least 47% of the total NOx emissions or a regional flux of 12.4 kg of N ha−1 year−1 (table S4). We consider this to be a conservative estimate because the county inventories make up a much larger area than the flight domain. Furthermore, the flights did not span either commuter rush hour, when NOx emissions are at their daily peak. We conclude that soils most likely contribute a majority of all NOx emissions to the atmosphere in the agriculturally intensive central San Joaquin Valley.

Comparing emissions estimates between methods

We compared surface emissions estimates for the San Joaquin Valley with soil model estimates for the same region (fig. S3 and table S4). Using year-round mean emissions for both natural and cropland soils, our model generated an annual flux of 24 kg of N ha−1 year−1 for the central San Joaquin Valley, between Fresno and Visalia, and as high as 36 kg of N ha−1 year−1 during the season of the airborne measurements (July to August), which yielded fluxes ranging from 14 to 39 kg of N ha−1 year−1. The correspondence between the bottom-up and top-down estimates builds robustness into our estimates for statewide NOx emissions and confirms our working conclusion for substantial NOx emissions from fertilized croplands in the Central Valley.

Our soil model estimates are slightly higher than, although comparable with, the few number of empirical measurements of NOx emissions from the San Joaquin Valley’s cropland soils (made between July and September of 1995) (14), which ranged from 0.1 to 14 kg of N ha−1 year−1 (Table 1). That we estimate higher soil NOx fluxes via the top-down and bottom-up approaches is consistent with more recent empirical measurements (13), suggesting that increases in N fertilizer use and population growth have likely accelerated soil NOx emissions since the last time empirical measurements were reported for the region, some 20 years ago (14).

Controls on soil NOx emissions

Previous studies have demonstrated that NOx emissions are controlled by water-filled pore space (influenced by precipitation, irrigation, and soil texture) (34), N availability (20), and temperature (13), all of which were fundamental parameters in our model-based estimate of soil NOx emissions in California (that is, the model relies on functions related to soil organic carbon, soil texture, drainage, temperature, and precipitation). We performed a model sensitivity analysis by evaluating the response of cropland denitrification rates to model input parameters at the ±10% level and examined the effect size of this parameter variation on NOx emissions. We found that soil NOx emissions were least sensitive to changes in soil carbon and were much more responsive to changes in soil texture, soil drainage, and climate (fig. S2).

Nitrogen input rates and climate were primary determinants of soil NOx emissions in our model. The largest chamber-based measurements of soil NOx emissions come from the Imperial Valley in Southern California (Table 1) (12), which was accurately predicted by our model, implying that our model is capable of detecting hot spot emissions. High emissions in the Imperial Valley are likely explained by three factors. First, a biogenic source in these soils suggests a kinetic response to high temperatures that occur in this region. Second, arid soils not only produce more NOx relative to N2O and N2 but also allow for the build-up of inorganic N via nitrification; N that will then be released in large quantities when soils are irrigated and microbial denitrification is triggered. Third, high fertilizer inputs that increase N availability in the soil may help soils to develop a healthy community of nitrifying bacteria, providing a positive feedback to N availability and subsequent loss.

Implications for California

The CARB emission inventory provides an assessment of air pollution magnitudes and sources in California. Sources are inventoried based on four main categories: mobile, stationary, area-wide, and natural. In the current CARB NOx inventory, mobile emissions are thought to predominate (83%), whereas soil emissions are currently considered negligible (16). Here, we show that agricultural soils contribute a substantial amount of NOx to the atmosphere. We can expect to see the significance of biogenic NOx emissions increase as N fertilizer inputs increase to keep pace with food demands (35) and automotive NOx controls continue to attenuate mobile fossil fuel sources. Our findings suggest the need to reconsider the role of soil NOx sources and provide a pathway to constrain these diffuse pathways into CARB inventory analyses. Recent climate changes in California have caused pronounced heat waves and drought, factors which could exacerbate biogenic NOx emissions, leading to increased air pollution and N deposition rates in natural ecosystems (8). Considering the limited number of field-based NOx measurements and the difficulty involved with partitioning soil versus fossil fuel sources through satellite imagery, a more robust field sampling strategy of soil NOx emissions throughout the state could aid in efforts to understand agricultural impacts on air pollution in the Central Valley.

Several existing approaches could be used to reduce soil NOx emissions from fertilized croplands. There are many strategies to improve fertilizer efficiencies, which would minimize the unwanted risks of N fertilizer spillovers into the environment and benefit farmers by reducing fertilizer costs. Where mineral fertilizers are used exclusively, for example, applying different forms of fertilizer (for example, slow-release fertilizers) (36) or lowering N applications and using precision agriculture to target developmental stages (37), have been shown to cut N fertilizer losses from cropland soil. Where organic amendments are applied, separating the application timing of mineral N and organic fertilizer has been shown to reduce N emissions (38). Precision fertilization, as opposed to broadcasting, can also increase N uptake and minimize losses (39). Cover crops that consume residual N, which can subsequently be incorporated into the soil, are another option for reducing N fertilizer application rates (40). A complementary institutional strategy would be to incentivize plant production for human versus livestock consumption because livestock manure and the N used to grow livestock feed are major sources of N pollution in the air and water (3). Another strategy would be to promote the reduction of NOx to an environmentally benign gas such as dinitrogen (N2), which can be achieved by installing riparian zones to collect fertilizer runoff or introducing nitrification inhibitors to stem denitrification rates (3, 41). The ratio in which harmful (NOx and N2O) and inert (N2) gases are emitted from soils depends heavily on N availability, soil moisture, and temperature; thus, irrigation strategies are another important step to reduce N losses from agriculture.

These and many other strategies can help to reduce potentially harmful N losses from agriculture (3, 42, 43). Losses of N fertilizer are not only costly to farmers but can also create economic costs to the greater United States on the order of $210 billion dollars per year in health and environmental damages (43, 44). Reducing NOx emissions therefore offers a win-win situation for farmers, environmental health, and the economy.


This study builds on local point-scale measurements (14) to provide the first spatially explicit evidence of substantial NOx emissions from agricultural soils in California, a previously unrecognized source that is estimated to contribute 20 to 51% of the state’s total NOx budget. These soil NOx emissions are sourced to N fertilizer applications in Central Valley croplands. The effect of large soil NOx emissions on air quality and human health remain unclear, but the magnitude of the flux alone raises concern about its potential impact, particularly in rural California. Where biogenic sources affect air quality and health, the implementation of strategies to reduce these emissions will be imperative. A better understanding of the sources, distribution, and impact of biogenically produced NOx will improve our ability to mitigate emissions in the future.


To model the spatial distribution of soil NOx emissions, we used an N isotope model (17, 45) in natural areas and an Integrated Model for the Assessment of the Global Environment (IMAGE) (46) in cropland areas to estimate total N losses from soils based on the surplus of N in the environment [see the study of Wang et al. (7) for model details]. The N surplus was a function of N inputs (deposition, fixation, fertilizer, and irrigation) minus N outputs besides denitrification and leaching (crop harvest and ammonia volatilization). Manure and grazing were not included as inputs/outputs; instead, we considered them as recycling functions of internal N cycling. Surplus N was then partitioned between leaching and gaseous losses based on temperature, precipitation (irrigation), evapotranspiration, soil texture, soil drainage, soil organic carbon content, soil total available water content, and land cover (see the Supplementary Materials). Gaseous losses were partitioned between NO, N2O, and N2 based on water-filled pore space (fig. S1). The resulting NO flux in California was between 0 and 276 kg of N ha−1 year−1 with one outlier (525.7 kg of N ha−1 year−1) being removed.

All data sets were transformed to 4000-m × 4000-m spatial grids before model runs. We used the California Augmented Multisource Land Cover Map published by the Information Center for the Environment, University of California, Davis (UC Davis). Temperature and precipitation data (mean from 1971 to 2000) were from the PRISM (Parameter-elevation Regressions on Independent Slopes Model) Climate Group, Oregon State University, created on 4 June 2010 (47). Nitrogen deposition data were based on the Community Multiscale Air Quality model by Tonnesen et al. (48). Nitrogen fertilizer application rates were from the Department of Water Resources (DWR) of California and U.S. Department of Agriculture (USDA) fertilizer consumption database for 1964 to 2006. Nitrogen harvest was calculated by multiplying crop production by percent dry matter by percent N for each crop type at each year. We then used the mean value for each crop type from 1980 to 2007. The data were from the USDA National Agricultural Statistics Service (2012) and the UC Davis Agricultural and Resource Economics. Nitrogen irrigation data were from the DWR irrigation database for 2001, assuming that N content in irrigation water was 1 mg of N liter−1 water. We used a fixed N fixation rate for different crops: rice, 2.5 kg of N ha−1 year−1; beans, 40 kg of N ha−1 year−1; alfalfa, 400 kg of N ha−1 year−1; clover, 15 kg of N ha−1 year−1; and natural vegetation, 5 kg of N ha−1 year−1 (10). The ammonia volatilization rate was estimated to be 3.6% of N fertilizer, as reported in the California Department of Food and Agriculture Fertilizer Research and Education Program, and to be 15 kg of N ha−1 year−1 in manure applied areas and 1.6 kg of N ha−1 year−1 in urban areas (49). Soil properties data, including soil texture, soil drainage, SOC (soil organic carbon) content, and total available water content, were from the USDA Soil Survey. We collected data from the literature, which we then compared to our model generated fluxes. Data were collected using Google Scholar and some combination of the following search terms: “NOx,” “nitrogen oxides,” “nitric oxide,” “emissions,” “trace gas,” “California,” and “site name.” References were then followed from papers that were found to be useful.

The airborne NO measurements were made with an Eco Physics (model CLD 88) chemiluminescence instrument with a stabilized photomultiplier tube and reaction chamber temperatures and other operating parameters to ensure a steady calibration point and high reproducibility. A blue-light light-emitting diode photolytic converter (Model 42i BLC2-395 manufactured by Air Quality Design Inc.) was used to selectively convert NO2 to NO for alternating measurements of NOx (= NO + NO2). Further, a Teflon prereaction chamber was installed to run the chemiluminescence reaction to completion before the detection cell to keep track of any interferences or changes in the background signal, thereby increasing the confidence in the measurements and lowering the detection limit generally to less than 0.05 parts per billion by volume. The instrument was cycled through the three states of NO and NOx, and background measurements were done every 20 s. Calibrations were performed by O3 titration with a National Institute of Standards and Technology traceable NO standard (Scott-Marrin Inc.) certified to within 5%. Full calibrations were performed before and after the entire flight series, with zero and span checks run routinely before and after each flight. The aircraft used to conduct the experiment is operated by Scientific Aviation Inc. (

Emission estimates were made using a simple boundary layer budget equation for NOx (see the Supplementary Materials). This technique was outlined in the study of Lenschow et al. (50), can be generalized to any scalar (5154), and was recently used to estimate regional methane emissions in the San Joaquin Valley (32). The technique involves thoroughly probing the atmospheric boundary layer (ABL) over a particular region via aircraft, horizontally and vertically, to determine the time rate of change, horizontal advection, and vertical mixing for various scalars, as well as the boundary layer height and its growth. This technique permits the calculation of residual terms within the scalar budgets for the region of interest (32). See the Supplementary Materials and Trousdell et al. (32) for greater details of the budget method, error estimates, and the other aircraft measurements.


Supplementary material for this article is available at

Supplementary Methods

table S1. Crop classification and fertilizer rate data (mean for 1964 to 2006) collected from the DWR of California and USDA fertilizer consumption database.

table S2. ABL heights, zi, and budget terms for the six flights.

table S3. NOx budget table and the consequent total regional emissions for each flight.

table S4. Flight estimates of total NOx and soil NOx and model estimates of soil NOx for the flight area in fig. S3 [Coordinates box: (36º51′52.09″N, 120º43′19.65″W), (37º0′6.85″ N, 119º50′53.87″W), (35º57′49.03″N, 120º1′37.93″W), and (36º5′27.03″N, 118º58′2.91″W) compared with CARB inventory of total NOx].

fig. S1. Model of how nitrogen oxide (NO), nitrous oxide (N2O), and dinitrogen (N2) partitioning varies with water-filled pore space.

fig. S2. Sensitivity of NO emission from croplands to different input parameters: soil organic carbon (fsoc), soil texture (ftxt), soil drainage (fdrain), and climate (fclim).

fig. S3. Airborne NOx observation sampling area.

References (5558)

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: We would like to thank T. O’Green and D. Liptzin for assistance with data acquisition and S. Chen for sharing her Weather Research and Forecasting (WRF) modeling runs. Funding: This work was supported by the NSF Postdoctoral Research Fellowship in Biology (1611895), the Major State Basic Research Development Program of China (2014CB954400), the NSF CAREER Grant (1411368), and the Dave and Lucile Packard Foundation. I.F. was supported in part by the California Agricultural Experiment Station, Hatch project CA-D-LAW-2229-H. The airborne NOx measurements were supported by a grant sponsored by the U.S. Environmental Protection Agency through the Bay Area Air Quality Management District (grant #2016-129) for the CARB’s CABOTS project. Author contributions: M.A., E.B., and B.Z.H. conceived the research and designed the study. M.A. collected literature data on NOx observations recorded in California. E.B. and B.Z.H developed the model. E.B and C.W. generated model estimates of soil NOx emissions. I.F., J.T., and S.C. performed airborne measurements of NOx and generated estimates of soil NOx emissions for the Central Valley in California. M.A. led the manuscript preparation with input from all co-authors. Competing interests: The authors declare that they have no competing interests and have filed no patents related to this work. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors. Please contact the corresponding author for additional data not reported in the manuscript or Supplementary Materials.
View Abstract

Navigate This Article