Carbon Sequestration and Nitrous Oxide Emission Mitigation Potential

8 downloads 0 Views 1MB Size Report
Shredding. 1. 1. 1.0. (1) Crop cycles per year refers to the number of times per year a particular farming operation is performed. A value less than one indicates ...
Carbon Sequestration and Nitrous Oxide Emission Mitigation Potential of California's Croplands: A County Scale Assessment Using DNDC and GIS Databases. William Salas1, Marc Los Huertos2 and Changsheng Li3 1Applied

GeoSolutions, LLC, 10 Newmarket Road, Durham, NH, 03824 [email protected] 2Center for Agroecology and Sustainable Food Systems, University of California Santa Cruz 3 Complex Systems Research Center, University of New Hampshire, Durham, NH 03824

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Outline • Overview of Model: DeNitrificationDeComposition or DNDC – Model structure and modeling approach – Validation – Scaling up from site to regional scale

• Ongoing Carbon Scoping Project: County scale assessment of carbon sequestration potential and trace gas emission from California croplands…

• Linking Science and Policy:Role of DSS Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Figure 2

The DNDC Model Ecological drivers

Climate

Vegetation

Soil

Anthropogenic activity

Plant growth Daily water demand

Daily potential ET

Annual average temperature

Daily biomass accumulation ( LAI )

Water uptake by roots

LAI-regulated albedo

Evaporation

Stalks

Soil temperature profile

Soil moisture profile

Oxygen diffusion

Soil Eh profile

Resistant litter

Labile microbes

Resistant microbes

Labile humads

Resistent humads

NH4

Daily N uptake by roots Root respiration

Labile litter

Grain

Water stress

Water flow between layers

Transpiration

Very labile litter

CO2

N demand

Roots

Oxygen consumption

DOC Passive humus

Soil climate

Effect of temperature and moisture on decomposition

Soil environmental variables

Temperature

NO2 -

Nitrate denitrifier

NO

Moisture

NO3

-

Nitrite denitrifier

N2O

pH

DOC

Decomposition

Substrates (NH4+, NO3- and DOC)

Eh

Nitrifiers

NH4+

NH3

NO3-

Clay-NH4+

Soil Eh

CH4 production

Aerenchyma

CH4 oxidation

CH4

DOC N2

Denitrification

CH4 transport

DOC

N2O denitrifier

N2O

NO

NH3

Nitrification

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

Fermentation

©Applied GeoSolutions, LLC

How DNDC Links Management to C and N Dynamics: Elemental cycling

Biochemical & geochemical reactions

Environmental factors

Ecological drivers

Mechanical movement Dissolution / crystallization

Gravity Radiation

Combination / decomposition Transport and transformation of chemical elements

Climate

Temperature

Soil

Moisture Oxidation / reduction Adsorption / desorption Complexation / decomplexation

Vegetation

Eh

Anthropogenic activity

pH Substrate concentration gradient

Assimilation / dissimilation Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

DNDC Links Ecological Drivers to Crop Yield/Trace Gas Emission

INPUT

PROCESSES

Climate - Temperature - Precipitation

Soil properties - Texture - Organic matter - Bulk density - pH

Management - Crop type - Tillage - Fertilization - Manure use - Irrigation - Grazing

OUTPUT

Used by microbes

Production of CO2, CH4, N2O, NO, N2, and NH3

DNDC containing fundamental biochemical & geochemical processes

Dynamics of soil water, NH4, NO3, and DOC

Competition

Used by plants

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

Growth of crop biomass

©Applied GeoSolutions, LLC

Input Parameters Required by DNDC 1. Climate:

- Daily max and min temperature; - Daily precipitation; - Atmospheric N deposition;

2. Soil:

- Bulk density; - Texture (clay fraction); - Total organic C content; - pH;

3. Management::

- Crop: type, cultivars, and rotation; - Tillage: timing and depth; - Fertilization: timing, type, amount, and depth; - Manure amendment: timing, type, and amount; - Irrigation: timing, amount, and pH; - Weeding: timing; - Grazing: livestock type, intensity, and season.

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

DNDC Predicts 1. 2. 3. 4.

Crop growth and yield; Soil organic C and N pools; N leaching; Gas emissions: CO2, CH4, N2O, NO, N2 and NH3

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

DNDC Validation • Approach: – Site scale validation – Independent researchers – Sites have covered wide range of agro-ecosystems, soils and climate conditions – California: need for field data…

• Regional Applications: scaling site to region – Uncertainty analyses based on variability of input parameters: Most Sensitive Factor (MSF), Monte Carlo, and Latin Hypercube Sampling (LHS)

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Testing DNDC at Site Scale

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Comparison on CO2 emissions from a silty loam soil in a tilled and fertilized winter wheat field in Columbia, Missouri 40 Tillage

Winter wheat field CO2 emission rate, kg C/ha/day

30

20

10

0 30

60

90

120

150

180

210

240

270

300

330

360

Julian day Measured CO2

Simulated total CO2

Simulated root respiration

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

150-Year Simulation for Soil C Dynamics in A Winter Wheat Field with Different Cropping Practices at Rothamsted Station, UK 90,000

80,000

Soil C, kg C/ha

70,000

Fallow year

Winter wheat field

60,000

50,000

40,000

30,000

20,000 1840

1850

1860

1870

1880

1890

1900

1910

1920

1930

1940

1950

1960

1970

1980

1990

Year Simulated Simulated Simulated Field Field Field control fertilized manured control fertilized manured Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Validated SOC dynamics across range of crops and climate conditions…

Comparison of observed (points) and modeled (lines) long-term SOC dynamics at the agricultural sites in the U.S., the U.K, Canada, Germany, Australia, and China. Field datasets are from USDA 1992, Odell et al. 1984, Jenkinson 1991, Grant (unpublished), Li et al. 1997, Eduardo (unpublished), and Li et al. 2003. Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Observed and DNDC-Modeled N2O Fluxes from Agricultural Soils in the U.S., Canada, the U.K., Germany, New Zealand, China, Japan, and Costa Rica

Modeled N2O flux, kg N/ha/year

1000

100

0.4

0.34 0. 0.43

0.41 0.

0.4

0.37

10

0.032

0.035 0.037

0.011

0.015

0.035 0.029 0.033

0.035

1 0.1

1

0.01 0.015

0.01 0.02 0.025 0.006

0.019

0.032 0.032 0.05

0.032

0.031

0.05 0.029

0.

0.028

0.029

0.

10

100

0.025

0.019

0.1

Observed N2O flux, kg N/ha/year

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

1000

DNDC: Modeling CH4 emissions from rice paddy

CH4 emission

Ebullition

Plant-mediated transport

CH4 oxidation

Soil CH4

Eh

Aerenchyma development

CH4 production

CO2

DOC

soilOxygen moisture

Flooding duration

Rhizodeposition

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

Decomposition

Root respirartion

©Applied GeoSolutions, LLC

CH4 fluxes from a paddy rice (cultivar Mars) plot (Plot 1) at Beaumont, Texas, 1994 Field data from Ron Sass, Rice University 20

C H4 f lu x , k g C /h a /d a y

15 10 5 0 -5 -10 -15 Julian day Modeled CH4 product ion

Modeled CH4 oxidat ion

Modeled CH4 f lux by plant

Modeled CH4 f lux by ebulit ion

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

Modeled CH4 f lux

Flied CH4 f lux

©Applied GeoSolutions, LLC

Comparison of Measured and Modeled CH4 Emissions from 5 Sites In the U.S., China and Italy

-1

-1

simulated mean CH4-Emissionen [g C ha d ]

6

5

Texas (Plot 1)

4

Texas (Plot 2)

3 Vercelli Texas (Plot 3) 2

1 Wuxian 0 0

1

2

3

4

5 -1

6

-1

measured mean CH4-Emissionen [g C ha d ]

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

N2O and NO from Forests: Comparisons between observed and modeled fluxes from 28 forest stands in Europe and the U.S. (from dissertation of Florian Stange, Fraunhofer Institute for Atmospheric Environmental Studies, GarmischPatenkerchin, Germany, 2000)

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

DNDC GHG, C and NO3 Validation Country Canada

Validation N2O fluxes from cropland

The U.K.

N2O fluxes from cropland and pasture

Germany

NO and N2O from forests

Denmark Austria Italy

NO and N2O from forests NO and N2O from forests N2O from cropland

Australia

CO2 and N2O from cropland and grassland N2O fluxes from cropland and pasture

New Zealand

USA

NO3 leaching from cropland

USA

CO2 and N2O from cropland and grassland

China

CO2, CH4and N2O from cropland and grassland

Japan Thailand Costa Rica

CO2 and N2O from cropland CH4 from cropland N2O from cropland

Application National N2O inventory for agricultural lands Developed UK-DNDC for national N2O inventory for agricultural lands National NO and N2O inventory for cropland and forests

National N2O inventory for agricultural lands

Reference Smith et al. 2002; Grant et al. 2002 Brown et al. 2001

Butterbach-Bahl et al. 2001; Stange et al. 2000 Stange et al. 2000 Stange et al. 2000 Mulligan 2002 (unpublished) Wang et al. 1997

Developed NZ-DNDC for national N2O inventory for agricultural lands

National C sequestration and N2O inventory for cropland National C sequestration, CH4 and N2O inventory for cropland

Regional N2O inventory for agricultural land

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

Andrew 2002 (unpublished)

Farahbakhshazad et al. 2003 (unpublished) Li et al. 1996 and 2002a

Xiu et al. 1999; Li et al. 2002b; Cai et al. 2002 Cai et al. 2002 Cai et al. 2002 Plant 1998 and 2000

©Applied GeoSolutions, LLC

Regional Applications of DNDC • Estimating uncertainties is a critical Issue in scaling site to region – Uncertainty analyses based on variability of input parameters: Most Sensitive Factor (MSF), Monte Carlo, and Latin Hypercube Sampling (LHS) – Outputs provided in ranges. (e.g. 3.5 to 5.4 kg N2O/ha)

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Scaling Up from Site to Regions Field & lab experiments

Statistical data collection

Remote sensing data acquisition

Model development: predicting biochemical & geochemical processes at site scale

GIS database construction: providing climate, soil, vegetation, and management data at regional scale

Soil fertility determined by soil organic matter storage

Crop yield

Modeling with DNDC

Remote sensing analysis: improving crop acreage data & providing phenology data

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

Emissions of CO2, CH4, N2O, NO, N2, and NH3

Leaching of nitrate

©Applied GeoSolutions, LLC

Sensitivity to Inputs: MSF 1200

C H 4 f lu x , k g C /h a /y r

1000

800

600

400

200

Cultivar

Draining times

Straw amandment

Soil texture

SOC

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

Temperature change

M in im u m

M a x im u m

M in im u m

M a x im u m

1 2

-2 -1

4%

2%

1%

C la y

S n a d y c la y lo a m

S a n d y lo a m Loam

1000 2000

0

1 2

0

C-4 8

V -7 7 G E-T

M a rs

0

Texture+SOC Extreme scenario

©Applied GeoSolutions, LLC

Uncertainty Analyses: Monte Carlo vs MSF Sensitivity Analysis: Fresno Cotton Monte Carlo Analysis: 1000 DNDC runs 300 MSF Range: 2272 to 2655 kg C/ha)

250 75% in MSF Range

Frequency

200 150 100 50 0 2100

2200

Demonstration Results Only:

2300

2400

2500

2600

2700

2800

Annual C Sequestration kg C/ha

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Frequency of N 2O Emissions: 1997 Cotton in Fresno 250

Frequency: 1000 Monte Carlo samples

MSF Range: 3.2 to 3.8 (kg N/ha) 200

150

100

50

0 2.8

2.9

3

Demonstration Results Only:

3.1

3.2

3.3

3.4

3.5

3.6

3.7

3.8

3.9

4

N2O Emission Rate (kg N/ha)

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Frequency of change in SOC for all pasture fields in Modoc, CA, 1997 600

500

Frequency, 1/2000

400

300

200

100

0 1700

1800

1900

2000

2100

2200

2300

2400

dSOC, kg C/ha/yr

Demonstration Results Only: Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Frequency of N2O fluxes from all pasture fields in Modoc, CA, 1997 600

500

Frequency, 1/2000

400

300

200

100

0 0.7

0.9

1.1

1.3

1.5

1.7

1.9

2.1

2.3

2.5

2.7

N2O flux, kg N/ha/yr

Demonstration Results Only: Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Estimating Carbon Sequestration Potential at County Scale: Why use GIS Process-based models like DNDC? • Spatial and temporal variability in – Climate (inter-annual variability) – Soils – Management Impacts

• Long-term impacts? Decade? Century? • C and N coupling: Effect on GWP!

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Carbon Sequestration and Trace Gas Emissions 60

118000

116000

50 114000

40

110000 Intensive tillage Notill

108000

106000

N2O flux, kg N/ha/yr

SOC, kg C/ha

112000

30

20

104000

10 102000

100000

0 1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

1

2

3

4

Year

5

6

7

8

9

10

11

12

13

14

15

16

17

Year

Figure 1: DNDC modeled changes in soil organic carbon content (SOC) and nitrous oxide (N2O) emissions from a corn-soybean rotation under two different management systems. Initial soil conditions were set to be identical with the same nominal climate conditions. It is clear that a constant or even site specific emission factor for N2O would not capture the temporal dynamic of emissions.

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

18

19

20

Evaluating Management Alternatives for Mitigation of GHG emissions GWPi : Σ (CO2i + N2Oi * 310 + CH4i * 21; where GWPi (kg CO2 equivalent/ha/yr) is the Global Warming Potential induced by scenario i; CO2i, N2Oi and CH4i are CO2 flux (kg C/ha/yr), N2O flux (kg N/ha/yr) and CH4 flux (kg C/ha/yr), respectively, induced by scenario i. Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

DNDC-modeled C sequestration, N2O emissions and their global warming potentials (GWP) for a corn-soybean rotation system with different tillage approaches in Adair County, Iowa from 1994-2014 C sequestration kg C/ha/yr

N2O flux

SOC-GWP

kg N/ha/yr

N2O-GWP

Net GWP

kg CO2 equivalent/ha/yr

Intensive tillage

125

11.5

-459

5615

5156

Notill

468

21.1

-1716

10301

8585

Critical need for models to assess long-term impacts of management decisions!

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

he

at Al -ric fa e lfa Fa ric llo e w C -ri on ce ve n C on tion se a rv l til at l io age n til la ge N C o on til tin l uo 1 u s m f id se loo 2 di a m n s id se on g d 3 m aso rai n id se n d as rai o n ns 20 N dr o 00 ai st ns r k 40 g s aw 00 tra am en kg wst C a dm ra w me ent -C nd am me en nt 0 dm r 50 esid en % ue t re 90 si inco % du rp re e in or a si c du or ted po e in co rate rp d or at ed Am m on iu m Am bic Ure m a rb a on o iu nat m e su lfa t N e itr at e

W

GWP, kg CO2 equivalent/ha/yr

Summary of Management Alternatives and GWP 90000

80000

70000

60000

50000

40000

30000

20000

10000

0

Alternative management

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003 ©Applied GeoSolutions, LLC

Ongoing Carbon Scoping Project: :

¾Objective: County scale assessment of carbon sequestration and trace gas emission from California croplands ¾Data… ¾Demonstrate model capabilities ¾Next steps…

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Input Data • • • • •

Climate Data Soils Crop Areas Management Practices Scenarios for Carbon Sequestration

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Climate Data Inputs: • Minimum Needs: Daily Tmax, Tmin, Precip • Sources: Station data: NCDC and CIMIS, gridded DAYMET (NCAR/U Mont) data • DAYMET: – Produces daily temp, precipitation, humidity and radiation based on station data – Performs interpolation based on “spatial convolution of a truncated Gaussian filter”

• Initial Analysis: Using 1997 and 1983 daily DAYMET data for station nearest County centroid. Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Temperature (deg C)

Daily Min/Max Temperature: Fresno County (DAYMET Data: 1997 and 1983) 45 40 35 30 25 20 15 10 5 0 -5 1

31

61

91

121

151

181

211

241

271

301

331

361

DOY 1997 Avg 1997: max 25.9 °C min 11.2 °C

1997 maxT (°C)

1997 minT (°C)

1983 maxT (°C)

1983 minT (°C)

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

Avg 1983: max 24.2 °C min 10.5 °C

©Applied GeoSolutions, LLC

Daily Precipitation: Fresno County (DAYMET 1997 and 1983)

Daily Precipitation (cm)

3.5 3 2.5 2 1.5 1 0.5 0 0

25

50

75

100 125 150 175 200 225 250 275 300 325 350 DOY

1997 Total: 27.5cm 1983 Total: 64.8cm

1997 Precip (cm)

1983 Precip (cm)

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Soils • NRCS STATSGO Soils data ¾DWR crop area mask ¾Derived area weighted statistics of range (min, max) in SOC, pH, Texture (%clay), and bulk density by county

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

5.0%

4.0%

3.0%

0.0%

0.0%

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003 SAN BENITO SAN SAN DIEGO SAN JOAQUIN SAN LUIS OBISPO SANTA SANTA CRUZ SHASTA SOLANO STANISLAUS SUTTER TEHAMA TRINITY TULARE TUOLUMNE YOLO YUBA

1.0% Content

2.0%

10.0%

AMADOR BUTTE COLUSA CONTRA COSTA DEL NORTE FRESNO GLENN HUMBOLDT IMPERIAL KERN KINGS LAKE LASSEN Clay LOS ANGELES MADERA MARIN MARIPOSA MERCED MODOC MONTEREY AMADOR PLACER BUTTE PLUMAS COLUSA RIVERSIDE CONTRA COSTA SACRAMENTO DEL NORTE SAN BENITO FRESNO SAN GLENN SAN DIEGO HUMBOLDT IMPERIAL SAN JOAQUIN KERN SAN LUIS OBISPO KINGS SANTA LAKE SANTA CRUZ LASSEN SHASTA LOS ANGELES SOLANO MADERA STANISLAUS MARIN SUTTER MARIPOSA TEHAMA MERCED MODOC TRINITY MONTEREY TULARE PLACER TUOLUMNE PLUMAS YOLO RIVERSIDE YUBA SACRAMENTO

SOC

County Soils: SOC

(Derived from NRCS STATSGO)

6.0%

County Soils: Clay Content (Derived from NRCS STATSGO)

60.0%

50.0%

40.0%

30.0%

20.0%

©Applied GeoSolutions, LLC

Bulk Density

7.5

6.5

5.5

1.0 AMADOR BUTTE COLUSA CONTRA COSTA DEL NORTE FRESNO GLENN HUMBOLDT IMPERIAL KERN KINGS LAKE LASSEN LOS ANGELES MADERA MARIN MARIPOSA MERCED MODOC MONTEREY PLACER PLUMAS RIVERSIDE SACRAMENTO SAN BENITO SAN SAN DIEGO SAN JOAQUIN SAN LUIS OBISPO SANTA SANTA CRUZ SHASTA SOLANO STANISLAUS SUTTER TEHAMA TRINITY TULARE TUOLUMNE YOLO YUBA

4.5 AMADOR BUTTE COLUSA CONTRA COSTA DEL NORTE FRESNO GLENN HUMBOLDT IMPERIAL KERN KINGS LAKE LASSEN LOS ANGELES MADERA MARIN MARIPOSA MERCED MODOC MONTEREY PLACER PLUMAS RIVERSIDE SACRAMENTO SAN BENITO SAN SAN DIEGO SAN JOAQUIN SAN LUIS OBISPO SANTA SANTA CRUZ SHASTA SOLANO STANISLAUS SUTTER TEHAMA TRINITY TULARE TUOLUMNE YOLO YUBA

Soil pH

County Soils: pH

(Derived from NRCS STATSGO)

9.5

8.5

County Soils: Bulk Density (Derived from NRCS STATSGO)

1.7

1.6

1.5

1.4

1.3

1.2

1.1

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003 ©Applied GeoSolutions, LLC

County Agricultural Data • Various Sources of California data: County Commissioners Reports, FRAP (Fire Resource & Assessment Program, CDF), NASS, DWR Æ All have pluses and minuses! • Used DWR mid-1990s data: – Sub-county spatial resolution – Based on Aerial Photos coupled with field surveys – Total crop area: 38,344km2 Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

county fips cotton sugarbeet corn sorghum beans,dry sunflower other field barley wheat oats other grain alfalfa other pastuci AMADOR 6005 0 23 196 0 0 0 191 81 81 81 81 117 914 BUTTE 6007 0 1139 917 0 1056 359 1374 0 5011 0 124 1321 3439 COLUSA 6011 2564 470 1940 234 2748 606 6534 84 11580 0 1989 3185 1633 CONTRA COS 6013 0 0 3995 0 164 0 1671 715 715 715 715 1622 2865 DEL NORTE 6015 0 0 0 0 0 0 22 0 0 0 0 0 3806 FRESNO 6019 178197 7272 9913 0 726 5 11147 8706 8706 8706 8706 35386 7172 GLENN 6021 333 376 7880 227 1653 2933 1682 489 6670 529 2881 5601 9959 HUMBOLDT 6023 0 0 106 0 37 0 154 0 0 0 957 192 16914 IMPERIAL 6025 123 27535 1988 0 50 0 653 1418 40825 152 1728 75197 18479 KERN 6029 97782 2300 13834 403 3510 0 9150 47908 0 0 0 43473 3730 KINGS 6031 108314 140 24753 1391 5 0 31972 7237 7237 7237 7237 17637 1619 LAKE 6033 0 0 0 0 0 0 291 6 25 300 499 149 2404 6037 0 0 0 0 0 0 0 0 0 0 0 0 0 LOS ANGELES LASSEN 6035 0 63 10 3 1296 0 0 24 1785 2253 1514 15252 27467 MADERA 6039 20393 159 6173 15 3515 0 733 11298 0 0 0 14637 4688 MARIN 6041 0 0 0 0 0 0 0 0 0 0 1912 0 645 MARIPOSA 6043 0 0 0 0 0 0 0 1 1 1 1 0 1327 MERCED 6047 37110 3486 24124 0 68 0 2365 3792 3792 3792 3792 32950 24151 MODOC 6049 0 1691 0 0 1895 0 201 8274 1619 376 2067 15480 39038 MONTEREY 6053 0 35 439 0 0 0 2154 3795 3795 3795 3795 937 1169 PLACER 6061 0 0 241 0 0 0 191 0 0 0 2467 6 7516 PLUMAS 6063 0 0 0 0 0 0 0 268 0 0 109 2586 13802 RIVERSIDE 6065 0 0 343 0 323 0 475 3879 10962 2146 1427 3511 3054 6067 0 1190 17179 390 0 0 5030 160 2651 0 5163 2988 12491 SACRAMENTO SAN BENITO 6069 0 85 1065 0 0 0 656 2355 2355 2355 2355 531 439 6071 0 0 549 0 0 0 481 114 0 1082 64 1076 1637 SAN BERNARD SAN DIEGO 6073 0 0 66 0 6868 648 4 38 1421 2107 0 29 2748 6077 0 2729 30212 247 20 0 8245 0 0 0 30445 24850 14200 SAN JOAQUIN 6079 0 0 53 0 1697 0 291 22788 2847 10186 16778 1711 2888 SAN LUIS OBIS SANTA BARBA 6083 0 0 413 0 0 0 34 0 2827 2928 0 1037 2565 SANTA CRUZ 6087 0 0 8 0 0 0 38 34 34 34 34 4 454 SHASTA 6089 0 227 4 5 2042 1633 171 258 23 185 874 2366 17016 SOLANO 6095 0 4107 8006 20 12675 1 3384 0 0 0 27628 9233 7586 STANISLAUS 6099 28 27 23957 0 3977 263 694 1406 1406 1406 1406 13989 25657 SUTTER 6101 268 1528 2830 99 173 88 10068 0 0 2863 5061 2040 2636 TEHAMA 6103 0 157 1027 0 5101 15 1107 247 2800 439 2848 2147 13006 TRINITY 6105 0 0 0 0 0 0 4 0 0 69 0 16 946 TULARE 6107 31613 1821 45749 1167 786 2855 2481 28769 0 0 0 39562 3311 TUOLUMNE 6109 0 0 0 0 0 0 0 0 0 0 0 230 230 YOLO 6113 2051 2052 17798 617 39 0 14367 0 0 0 30002 16865 0 YUBA 6115 0 0 379 0 0 0 246 0 0 154 399 5802 0

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

GIS Database

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

DWR Crop Areas (km^2) Total Crop Area: 38,344 km^2 6000 5000 4000 3000 2000 1000

ck

t ru

f ru

it

&

nu

tT

re

e co s gr - l tt on ap et es tu ce -w i n al f a e & lf a ta pa bl e st ur e m c or isc n gr r ai n ic e (o at s ba ) rle m y isc c it ru f ie ld w s - s he or a gh t um be et oa t be s s u an nf s lo s o we rg r hu m

0

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Management Data COTTON Farming Operations Land Preparation Land Preparation Seed Bed Preparation

Crop Cycles Per Year(1)

Passes Per Crop Cycle(2)

Fraction Acreage Per Cycle(3)

1 1

4 2

1.0 1.0

1

1

1.0

1

3

1.0

1

1

1.0

1

1

1.0

Passes During Month Jan Feb

Mar

Apr

May

Jun

Jul

Aug

Sep

Oct

Nov

Dec

Planting

Growing Season Operation

Harvesting

Postharvest Shredding

(1) Crop cycles per year refers to the number of times per year a particular farming operation is performed. A value less than one indicates an operation is performed less than once per year. Values greater than one indicate the operation is done more than once per year.

Land Prep Acre-Passes Land prep 1x4x1=4 Seed bed prep 1 x 2 x 1 = 2 ----------------Total acre passes 6

(2) Passes per crop cycle refers to the actual number of passes by a farm implement Necessary to accomplish a particular farming operation. (3) Fraction acreage per cycle refers to the fraction of the acreage covered by the particular farming operation. For example, in an orchard or a vineyard, operations usually only disturb the ground between the rows. In those cases only 50% of the acreage is actually affected by the operation. In contrast, a discing operation usually affects 100% of the acreage.

Source: CARB Fugitive Dust Study Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Management Data FIELD CORN Farming Operations Land Preparation Stubble Disc Finish Disc List & Fertilize Mulch Beds

Crop Cycles Per Year(1)

Passes Per Crop Cycle(2)

Fraction Acreage Per Cycle(3)

Passes During Month Jan Feb

1 1 1 1

1 1 1 1

1.0 1.0 1.0 1.0

1

1

1.0

1

2

1.0

1

1

1.0

Mar

Apr

May

Jun

Jul

Aug

Sep

Oct

Nov

Dec

Planting

Cultivation

Harvesting

Postharvest

(1) Crop cycles per year refers to the number of times per year a particular farming operation is performed. A value less than one indicates an operation is performed less than once per year. Values greater than one indicate the operation is done more than once per year.

Land Prep. Acre-Passes Stubble Disc 1 x 1 x 1 = 1 Finish Disc 1 x 1 x 1 = 1 List 1x1x1=1 Mulch 1x1x1=1 ----------------Total acre passes 4

(2) Passes per crop cycle refers to the actual number of passes by a farm implement Necessary to accomplish a particular farming operation. (3) Fraction acreage per cycle refers to the fraction of the acreage covered by the particular farming operation. For example, in an orchard or a vineyard, operations usually only disturb the ground between the rows. In those cases only 50% of the acreage is actually affected by the operation. In contrast, a discing operation usually affects 100% of the acreage.

Source: CARB Fugitive Dust Study Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Fertilizer Application Rates DWR Crop Class C D D12 F G I P P1 R T V

Description

Fertilizer Application Rates (kg N/ha)

Citrus & Subtropical Fruit Decid Fruit & Nuts (non-almond) Almonds Field Crops Grain and Hay Idle Pasture Grass (non N-fixing) Alfalfa Pasture Rice Truck Crops (vegetables) Vineyard

SJV 150 140 100 140 100 0 100 0 100 250 70

Sac 140 110 80 140 100 0 100 0 100 210 60

CCst 140 110 80 140 100 0 100 0 100 300 60

Imp 125 140 100 140 100 0 100 0 100 300 75

Initial Baseline Analysis did not differentiate across regions Source: Potter et al 2001. (analysis of NH3 emissions for CARB) Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Next Steps • Fertilizer: use different application rates across regions DWR Crop Class C D D12 F G I P P1 R T V

Description

Citrus & Subtropical Fruit Decid Fruit & Nuts (non-almond) Almonds Field Crops Grain and Hay Idle Pasture Grass (non N-fixing) Alfalfa Pasture Rice Truck Crops (vegetables) Vineyard

Fertilizer Application Rates (kg N/ha) SJV 150 140 100 140 100 0 100 0 100 250 70

Sac 140 110 80 140 100 0 100 0 100 210 60

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

CCst 140 110 80 140 100 0 100 0 100 300 60

Imp 125 140 100 140 100 0 100 0 100 300 75

©Applied GeoSolutions, LLC

Next Steps cont. • Soils: Use crop class specific soils data at the county scale. Merge DWR and STATSGO • Merced Soil Characteristics: CROP G R F P T D/C/V I

AREA(m^2) 154433722.4 26424198.3 721508949.3 574194178.3 163349689.3 537401774.7 NA

CLAY FRACTION A.W.Min A.W.Max 17.63561 25.79246 19.61089 28.58731 21.73627 28.47172 18.52837 25.73746 19.36045 25.1621 8.545795 14.71284 NA NA

BULK DENSITY A.W.Min A.W.Max 1.394909 1.495448 1.343949 1.443949 1.400843 1.503235 1.415494 1.517749 1.43223 1.532635 1.513845 1.61709 NA NA

ORGANIC MATTER A.W.Min A.W.Max 0.671436 1.804702 0.603981 1.252003 0.797725 2.348781 0.734353 2.109567 0.755769 1.796073 0.625974 1.565343 NA NA

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

pH A.W.Min 6.966689 7.760544 7.069387 6.981879 6.771935 6.189352 NA

A.W.Max 8.463362 9.318588 8.475595 8.473096 8.097954 7.754332 NA

©Applied GeoSolutions, LLC

Next Steps cont. • Validation analyses for California. Critical! Need to collect existing data. Source? – Long-term SOC changes. – N2O data – CH4 from rice

• Evaluate scenarios for C sequestration: cover crops, conservation tillage, notill, climate change, … • Run 20 and 40 year scenarios to examine C sequestration capacity and net GWP (N2O offsets)

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC

Linking Science to Policy: Role for Decision Support Systems GIS Data Server: Climate: NCDC & DAYMET Soils: NRCS STATSGO & SUSRGO

DNDC Biogeochemical Processor:

Ag Census NASS. DWR Land use, etc

Long-term soil fertility: SOC Emissions of CO2, CH4, N2O, N2, and NH3

Field Studies: Model refinement Validation Auditing/verification

Remote Sensing Product Generation Data: IKONOS, Landsat ETM, MODIS RS Products: Crop Type, Crop Phenology, including planting and harvesting dates and LAI, Management Data, including • Tillage practice (conventional, vs no-till) • Tillage timing • Use of cover crops • Irrigation • Residue management

Nitrate Leaching

Management Parameter Server: -system will contain defaults for all management parameters (e.g. fertilizer, manure, tillage, crop cycles, crop types, irrigation, etc. ) -users can then change management parameters and modify soil properties on a field by field basis if they have access to better soils data. Individual fields can be selected using GIS tools and aerial photos from aerial photo server (see figure for example).

Record Keeping System -tracks historical management -keeps inventories of GHG, SOC, NH3 and nitrate fluxes.

GIS Map Server - creates user defined maps. Sensitivity and Error Tracking System -

provides estimates sensitivity estimates tracks differences between field and Model estimates

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

Report Generator -produces summary reports with maps - incorporates regulatory guidelines regarding fert and manure applications

©Applied GeoSolutions, LLC

Thank You! Acknowledgments: Ongoing California Countyscale analysis is supported by joint Kearney Foundation California Energy Commission and California Department of Food and Agriculture project.

Presented at Soil Carbon Sequestration Workshop, Kearney Foundation Soil Science, September 2003

©Applied GeoSolutions, LLC