Système d'information stratégique et agriculture (serveur d'exploration)

Attention, ce site est en cours de développement !
Attention, site généré par des moyens informatiques à partir de corpus bruts.
Les informations ne sont donc pas validées.

Present-day oxidative subsidence of organic soils and mitigation in the Sacramento-San Joaquin Delta, California, USA

Identifieur interne : 000055 ( Pmc/Corpus ); précédent : 000054; suivant : 000056

Present-day oxidative subsidence of organic soils and mitigation in the Sacramento-San Joaquin Delta, California, USA

Auteurs : Steven J. Deverel ; Timothy Ingrum ; David Leighton

Source :

RBID : PMC:4944668

Abstract

Subsidence of organic soils in the Sacramento-San Joaquin Delta threatens sustainability of the California (USA) water supply system and agriculture. Land-surface elevation data were collected to assess present-day subsidence rates and evaluate rice as a land use for subsidence mitigation. To depict Delta-wide present-day rates of subsidence, the previously developed SUBCALC model was refined and calibrated using recent data for CO2 emissions and land-surface elevation changes measured at extensometers. Land-surface elevation change data were evaluated relative to indirect estimates of subsidence and accretion using carbon and nitrogen flux data for rice cultivation. Extensometer and leveling data demonstrate seasonal variations in land-surface elevations associated with groundwater-level fluctuations and inelastic subsidence rates of 0.5–0.8 cm yr–1. Calibration of the SUBCALC model indicated accuracy of ±0.10 cm yr–1 where depth to groundwater, soil organic matter content and temperature are known. Regional estimates of subsidence range from <0.3 to >1.8 cm yr–1. The primary uncertainty is the distribution of soil organic matter content which results in spatial averaging in the mapping of subsidence rates. Analysis of leveling and extensometer data in rice fields resulted in an estimated accretion rate of 0.02–0.8 cm yr–1. These values generally agreed with indirect estimates based on carbon fluxes and nitrogen mineralization, thus preliminarily demonstrating that rice will stop or greatly reduce subsidence. Areas below elevations of –2 m are candidate areas for implementation of mitigation measures such as rice because there is active subsidence occurring at rates greater than 0.4 cm yr–1.


Url:
DOI: 10.1007/s10040-016-1391-1
PubMed: 27471427
PubMed Central: 4944668

Links to Exploration step

PMC:4944668

Le document en format XML

<record>
<TEI>
<teiHeader>
<fileDesc>
<titleStmt>
<title xml:lang="en">Present-day oxidative subsidence of organic soils and mitigation in the Sacramento-San Joaquin Delta, California, USA</title>
<author>
<name sortKey="Deverel, Steven J" sort="Deverel, Steven J" uniqKey="Deverel S" first="Steven J." last="Deverel">Steven J. Deverel</name>
<affiliation>
<nlm:aff id="Aff1">HydroFocus, Inc., 2827 Spafford Street, Davis, CA 95618 USA</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Ingrum, Timothy" sort="Ingrum, Timothy" uniqKey="Ingrum T" first="Timothy" last="Ingrum">Timothy Ingrum</name>
<affiliation>
<nlm:aff id="Aff1">HydroFocus, Inc., 2827 Spafford Street, Davis, CA 95618 USA</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Leighton, David" sort="Leighton, David" uniqKey="Leighton D" first="David" last="Leighton">David Leighton</name>
<affiliation>
<nlm:aff id="Aff1">HydroFocus, Inc., 2827 Spafford Street, Davis, CA 95618 USA</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt>
<idno type="wicri:source">PMC</idno>
<idno type="pmid">27471427</idno>
<idno type="pmc">4944668</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4944668</idno>
<idno type="RBID">PMC:4944668</idno>
<idno type="doi">10.1007/s10040-016-1391-1</idno>
<date when="2016">2016</date>
<idno type="wicri:Area/Pmc/Corpus">000055</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000055</idno>
</publicationStmt>
<sourceDesc>
<biblStruct>
<analytic>
<title xml:lang="en" level="a" type="main">Present-day oxidative subsidence of organic soils and mitigation in the Sacramento-San Joaquin Delta, California, USA</title>
<author>
<name sortKey="Deverel, Steven J" sort="Deverel, Steven J" uniqKey="Deverel S" first="Steven J." last="Deverel">Steven J. Deverel</name>
<affiliation>
<nlm:aff id="Aff1">HydroFocus, Inc., 2827 Spafford Street, Davis, CA 95618 USA</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Ingrum, Timothy" sort="Ingrum, Timothy" uniqKey="Ingrum T" first="Timothy" last="Ingrum">Timothy Ingrum</name>
<affiliation>
<nlm:aff id="Aff1">HydroFocus, Inc., 2827 Spafford Street, Davis, CA 95618 USA</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Leighton, David" sort="Leighton, David" uniqKey="Leighton D" first="David" last="Leighton">David Leighton</name>
<affiliation>
<nlm:aff id="Aff1">HydroFocus, Inc., 2827 Spafford Street, Davis, CA 95618 USA</nlm:aff>
</affiliation>
</author>
</analytic>
<series>
<title level="j">Hydrogeology Journal</title>
<idno type="ISSN">1431-2174</idno>
<idno type="eISSN">1435-0157</idno>
<imprint>
<date when="2016">2016</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc>
<textClass></textClass>
</profileDesc>
</teiHeader>
<front>
<div type="abstract" xml:lang="en">
<p>Subsidence of organic soils in the Sacramento-San Joaquin Delta threatens sustainability of the California (USA) water supply system and agriculture. Land-surface elevation data were collected to assess present-day subsidence rates and evaluate rice as a land use for subsidence mitigation. To depict Delta-wide present-day rates of subsidence, the previously developed SUBCALC model was refined and calibrated using recent data for CO
<sub>2</sub>
emissions and land-surface elevation changes measured at extensometers. Land-surface elevation change data were evaluated relative to indirect estimates of subsidence and accretion using carbon and nitrogen flux data for rice cultivation. Extensometer and leveling data demonstrate seasonal variations in land-surface elevations associated with groundwater-level fluctuations and inelastic subsidence rates of 0.5–0.8 cm yr
<sup>–1</sup>
. Calibration of the SUBCALC model indicated accuracy of ±0.10 cm yr
<sup>–1</sup>
where depth to groundwater, soil organic matter content and temperature are known. Regional estimates of subsidence range from <0.3 to >1.8 cm yr
<sup>–1</sup>
. The primary uncertainty is the distribution of soil organic matter content which results in spatial averaging in the mapping of subsidence rates. Analysis of leveling and extensometer data in rice fields resulted in an estimated accretion rate of 0.02–0.8 cm yr
<sup>–1</sup>
. These values generally agreed with indirect estimates based on carbon fluxes and nitrogen mineralization, thus preliminarily demonstrating that rice will stop or greatly reduce subsidence. Areas below elevations of –2 m are candidate areas for implementation of mitigation measures such as rice because there is active subsidence occurring at rates greater than 0.4 cm yr
<sup>–1</sup>
.</p>
</div>
</front>
<back>
<div1 type="bibliography">
<listBibl>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Boumans, Rmj" uniqKey="Boumans R">RMJ Boumans</name>
</author>
<author>
<name sortKey="Day, Jwj" uniqKey="Day J">JWJ Day</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Browder, Ja" uniqKey="Browder J">JA Browder</name>
</author>
<author>
<name sortKey="Volk, Bg" uniqKey="Volk B">BG Volk</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Carlton, Ab" uniqKey="Carlton A">AB Carlton</name>
</author>
<author>
<name sortKey="Schultz, Hb" uniqKey="Schultz H">HB Schultz</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Cook, J" uniqKey="Cook J">J Cook</name>
</author>
<author>
<name sortKey="Coleman, D" uniqKey="Coleman D">D Coleman</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Couwenberg, J" uniqKey="Couwenberg J">J Couwenberg</name>
</author>
<author>
<name sortKey="Hooijer, A" uniqKey="Hooijer A">A Hooijer</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Davidson, Ea" uniqKey="Davidson E">EA Davidson</name>
</author>
<author>
<name sortKey="Sudeep, S" uniqKey="Sudeep S">S Sudeep</name>
</author>
<author>
<name sortKey="Caramori, Ss" uniqKey="Caramori S">SS Caramori</name>
</author>
<author>
<name sortKey="Savage, K" uniqKey="Savage K">K Savage</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Deverel, Sj" uniqKey="Deverel S">SJ Deverel</name>
</author>
<author>
<name sortKey="Rojstaczer, Sa" uniqKey="Rojstaczer S">SA Rojstaczer</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Deverel, Sj" uniqKey="Deverel S">SJ Deverel</name>
</author>
<author>
<name sortKey="Wang, B" uniqKey="Wang B">B Wang</name>
</author>
<author>
<name sortKey="Rojstaczer, Sa" uniqKey="Rojstaczer S">SA Rojstaczer</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Drexler, Jz" uniqKey="Drexler J">JZ Drexler</name>
</author>
<author>
<name sortKey="De Fontaine, Cs" uniqKey="De Fontaine C">CS de Fontaine</name>
</author>
<author>
<name sortKey="Brown, Ta" uniqKey="Brown T">TA Brown</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Drexler, Jz" uniqKey="Drexler J">JZ Drexler</name>
</author>
<author>
<name sortKey="De Fontaine, Cs" uniqKey="De Fontaine C">CS de Fontaine</name>
</author>
<author>
<name sortKey="Deverel, Sj" uniqKey="Deverel S">SJ Deverel</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Ewing, Jm" uniqKey="Ewing J">JM Ewing</name>
</author>
<author>
<name sortKey="Vepraskas, Mj" uniqKey="Vepraskas M">MJ Vepraskas</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Fleck, Ja" uniqKey="Fleck J">JA Fleck</name>
</author>
<author>
<name sortKey="Bossio, Da" uniqKey="Bossio D">DA Bossio</name>
</author>
<author>
<name sortKey="Fujii, R" uniqKey="Fujii R">R Fujii</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Florsheim, Jl" uniqKey="Florsheim J">JL Florsheim</name>
</author>
<author>
<name sortKey="Dettinger, Md" uniqKey="Dettinger M">MD Dettinger</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Gambolati, G" uniqKey="Gambolati G">G Gambolati</name>
</author>
<author>
<name sortKey="Putti, M" uniqKey="Putti M">M Putti</name>
</author>
<author>
<name sortKey="Teatini, P" uniqKey="Teatini P">P Teatini</name>
</author>
<author>
<name sortKey="Gasparetto Stori, G" uniqKey="Gasparetto Stori G">G Gasparetto Stori</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Hatala, Ja" uniqKey="Hatala J">JA Hatala</name>
</author>
<author>
<name sortKey="Detto, M" uniqKey="Detto M">M Detto</name>
</author>
<author>
<name sortKey="Sonnentag, O" uniqKey="Sonnentag O">O Sonnentag</name>
</author>
<author>
<name sortKey="Deverel, Sj" uniqKey="Deverel S">SJ Deverel</name>
</author>
<author>
<name sortKey="Verfaillie, J" uniqKey="Verfaillie J">J Verfaillie</name>
</author>
<author>
<name sortKey="Baldocchi, Dd" uniqKey="Baldocchi D">DD Baldocchi</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Hirano, T" uniqKey="Hirano T">T Hirano</name>
</author>
<author>
<name sortKey="Segah, H" uniqKey="Segah H">H Segah</name>
</author>
<author>
<name sortKey="Kusin, K" uniqKey="Kusin K">K Kusin</name>
</author>
<author>
<name sortKey="Limin, S" uniqKey="Limin S">S Limin</name>
</author>
<author>
<name sortKey="Takahashi, H" uniqKey="Takahashi H">H Takahashi</name>
</author>
<author>
<name sortKey="Osaki, M" uniqKey="Osaki M">M Osaki</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Hobbs, Nb" uniqKey="Hobbs N">NB Hobbs</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Hooijer, A" uniqKey="Hooijer A">A Hooijer</name>
</author>
<author>
<name sortKey="Page, S" uniqKey="Page S">S Page</name>
</author>
<author>
<name sortKey="Jauhiainen, J" uniqKey="Jauhiainen J">J Jauhiainen</name>
</author>
<author>
<name sortKey="Lee, Wa" uniqKey="Lee W">WA Lee</name>
</author>
<author>
<name sortKey="Lu, Xx" uniqKey="Lu X">XX Lu</name>
</author>
<author>
<name sortKey="Idris, A" uniqKey="Idris A">A Idris</name>
</author>
<author>
<name sortKey="Anshari, G" uniqKey="Anshari G">G Anshari</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Kirk, Er" uniqKey="Kirk E">ER Kirk</name>
</author>
<author>
<name sortKey="Van Kessel, C" uniqKey="Van Kessel C">C van Kessel</name>
</author>
<author>
<name sortKey="Horwath, Wr" uniqKey="Horwath W">WR Horwath</name>
</author>
<author>
<name sortKey="Linquist, Ba" uniqKey="Linquist B">BA Linquist</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Knox, Sh" uniqKey="Knox S">SH Knox</name>
</author>
<author>
<name sortKey="Sturtevant, C" uniqKey="Sturtevant C">C Sturtevant</name>
</author>
<author>
<name sortKey="Hatala Matthes, J" uniqKey="Hatala Matthes J">J Hatala Matthes</name>
</author>
<author>
<name sortKey="Koteen, L" uniqKey="Koteen L">L Koteen</name>
</author>
<author>
<name sortKey="Verfaillie, J" uniqKey="Verfaillie J">J Verfaillie</name>
</author>
<author>
<name sortKey="Baldocchi, D" uniqKey="Baldocchi D">D Baldocchi</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Prokopovitch, Np" uniqKey="Prokopovitch N">NP Prokopovitch</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Rojstaczer, Sa" uniqKey="Rojstaczer S">SA Rojstaczer</name>
</author>
<author>
<name sortKey="Deverel, Sj" uniqKey="Deverel S">SJ Deverel</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Rojstaczer, Sa" uniqKey="Rojstaczer S">SA Rojstaczer</name>
</author>
<author>
<name sortKey="Deverel, Sj" uniqKey="Deverel S">SJ Deverel</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Schipper, La" uniqKey="Schipper L">LA Schipper</name>
</author>
<author>
<name sortKey="Mcleod, M" uniqKey="Mcleod M">M McLeod</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Schorthorst, Cj" uniqKey="Schorthorst C">CJ Schorthorst</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Shlemon, Rj" uniqKey="Shlemon R">RJ Shlemon</name>
</author>
<author>
<name sortKey="Begg, El" uniqKey="Begg E">EL Begg</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Stephens, Jc" uniqKey="Stephens J">JC Stephens</name>
</author>
<author>
<name sortKey="Allen, Lh" uniqKey="Allen L">LH Allen</name>
</author>
<author>
<name sortKey="Chen, E" uniqKey="Chen E">E Chen</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Strack, M" uniqKey="Strack M">M Strack</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Syvitski, Jp" uniqKey="Syvitski J">JP Syvitski</name>
</author>
<author>
<name sortKey="Kettner, Aj" uniqKey="Kettner A">AJ Kettner</name>
</author>
<author>
<name sortKey="Overeem, I" uniqKey="Overeem I">I Overeem</name>
</author>
<author>
<name sortKey="Hutton, Ew" uniqKey="Hutton E">EW Hutton</name>
</author>
<author>
<name sortKey="Hannon, Mt" uniqKey="Hannon M">MT Hannon</name>
</author>
<author>
<name sortKey="Brakenridge, Gr" uniqKey="Brakenridge G">GR Brakenridge</name>
</author>
<author>
<name sortKey="Day, J" uniqKey="Day J">J Day</name>
</author>
<author>
<name sortKey="Vorosmarty, J" uniqKey="Vorosmarty J">J Vörösmarty</name>
</author>
<author>
<name sortKey="Saito, Y" uniqKey="Saito Y">Y Saito</name>
</author>
<author>
<name sortKey="Giosan, L" uniqKey="Giosan L">L Giosan</name>
</author>
<author>
<name sortKey="Nicholls, Rj" uniqKey="Nicholls R">RJ Nicholls</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Tate, Rl" uniqKey="Tate R">RL Tate</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Tate, Rl" uniqKey="Tate R">RL Tate</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Tate, Rl" uniqKey="Tate R">RL Tate</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Terzaghi, K" uniqKey="Terzaghi K">K Terzaghi</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Volk, Bg" uniqKey="Volk B">BG Volk</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Weir, Ww" uniqKey="Weir W">WW Weir</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article">
<pmc-dir>properties open_access</pmc-dir>
<front>
<journal-meta>
<journal-id journal-id-type="nlm-ta">Hydrogeol J</journal-id>
<journal-id journal-id-type="iso-abbrev">Hydrogeol J</journal-id>
<journal-title-group>
<journal-title>Hydrogeology Journal</journal-title>
</journal-title-group>
<issn pub-type="ppub">1431-2174</issn>
<issn pub-type="epub">1435-0157</issn>
<publisher>
<publisher-name>Springer Berlin Heidelberg</publisher-name>
<publisher-loc>Berlin/Heidelberg</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="pmid">27471427</article-id>
<article-id pub-id-type="pmc">4944668</article-id>
<article-id pub-id-type="publisher-id">1391</article-id>
<article-id pub-id-type="doi">10.1007/s10040-016-1391-1</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Paper</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Present-day oxidative subsidence of organic soils and mitigation in the Sacramento-San Joaquin Delta, California, USA</article-title>
<trans-title-group>
<trans-title xml:lang="fr">Affaissement actuel de sols organiques par oxydation et atténuation dans le Delta Sacramento-San Joaquin, Californie, Etats-Unis d’Amérique</trans-title>
</trans-title-group>
<trans-title-group>
<trans-title xml:lang="es">La subsidencia actual de los suelos orgánicos y la mitigación en el Delta de Sacramento-San Joaquin, California, EE.UU</trans-title>
</trans-title-group>
<trans-title-group>
<trans-title xml:lang="zh">美国加利佛尼亚州沙加缅度河—圣华金河三角洲有机土壤目前的氧化性下沉及减缓措施</trans-title>
</trans-title-group>
<trans-title-group>
<trans-title xml:lang="pt">Atual subsidência de solos orgânicos por oxidação e sua mitigação no Delta de Sacramento-San Joaquin, Califórnia, EUA</trans-title>
</trans-title-group>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Deverel</surname>
<given-names>Steven J.</given-names>
</name>
<address>
<email>sdeverel@hydrofocus.com</email>
</address>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Ingrum</surname>
<given-names>Timothy</given-names>
</name>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Leighton</surname>
<given-names>David</given-names>
</name>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<aff id="Aff1">HydroFocus, Inc., 2827 Spafford Street, Davis, CA 95618 USA</aff>
</contrib-group>
<pub-date pub-type="epub">
<day>28</day>
<month>3</month>
<year>2016</year>
</pub-date>
<pub-date pub-type="pmc-release">
<day>28</day>
<month>3</month>
<year>2016</year>
</pub-date>
<pub-date pub-type="ppub">
<year>2016</year>
</pub-date>
<volume>24</volume>
<fpage>569</fpage>
<lpage>586</lpage>
<history>
<date date-type="received">
<day>13</day>
<month>7</month>
<year>2015</year>
</date>
<date date-type="accepted">
<day>26</day>
<month>2</month>
<year>2016</year>
</date>
</history>
<permissions>
<copyright-statement>© The Author(s) 2016</copyright-statement>
<license license-type="OpenAccess">
<license-p>
<bold>Open Access</bold>
This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.</license-p>
</license>
</permissions>
<abstract id="Abs1">
<p>Subsidence of organic soils in the Sacramento-San Joaquin Delta threatens sustainability of the California (USA) water supply system and agriculture. Land-surface elevation data were collected to assess present-day subsidence rates and evaluate rice as a land use for subsidence mitigation. To depict Delta-wide present-day rates of subsidence, the previously developed SUBCALC model was refined and calibrated using recent data for CO
<sub>2</sub>
emissions and land-surface elevation changes measured at extensometers. Land-surface elevation change data were evaluated relative to indirect estimates of subsidence and accretion using carbon and nitrogen flux data for rice cultivation. Extensometer and leveling data demonstrate seasonal variations in land-surface elevations associated with groundwater-level fluctuations and inelastic subsidence rates of 0.5–0.8 cm yr
<sup>–1</sup>
. Calibration of the SUBCALC model indicated accuracy of ±0.10 cm yr
<sup>–1</sup>
where depth to groundwater, soil organic matter content and temperature are known. Regional estimates of subsidence range from <0.3 to >1.8 cm yr
<sup>–1</sup>
. The primary uncertainty is the distribution of soil organic matter content which results in spatial averaging in the mapping of subsidence rates. Analysis of leveling and extensometer data in rice fields resulted in an estimated accretion rate of 0.02–0.8 cm yr
<sup>–1</sup>
. These values generally agreed with indirect estimates based on carbon fluxes and nitrogen mineralization, thus preliminarily demonstrating that rice will stop or greatly reduce subsidence. Areas below elevations of –2 m are candidate areas for implementation of mitigation measures such as rice because there is active subsidence occurring at rates greater than 0.4 cm yr
<sup>–1</sup>
.</p>
</abstract>
<trans-abstract xml:lang="fr" id="Abs2">
<title>Résumé</title>
<p>L’affaissement des sols organiques dans le Delta de Sacramento-San Joaquin menace la durabilité du système d’approvisionnement en eau et de l’agriculture de la Californie (Etats-Unis d’Amérique). Des données d’altitude de la surface topographique ont été collectées dans le but d’estimer le taux d’affaissement actuel des sols et d’évaluer l’occupation des sols par des rizières en tant que moyen d’atténuation de l’affaissement des sols. Pour décrire les taux actuels d’affaissement des sols à l’échelle du Delta, le modèle SUBCALC développé précédemment a été affiné et calé grâce à l’utilisation des données récentes sur les émissions de CO
<sub>2</sub>
et des variations d’altitude de la surface topographique enregistrées par des extensomètres. Les données de variations d’altitude de la surface topographique ont été évaluées par rapport aux estimations indirectes de l’affaissement du sol et de son accroissement basées sur les données de flux de carbone et d’azote au niveau des rizières. Les données d’extensiométrie et de nivellement montrent que les variations saisonnières de l’altitude de la surface topographique associées aux fluctuations des niveaux piézométriques et des taux d’affaissement inélastique du sol compris entre 0.5 et 0.8 cm an
<sup>–1</sup>
. Le calage du modèle SUBCALC indique une précision de ±0.10 cm an
<sup>–1</sup>
là où la profondeur de la nappe, la teneur en matière organique du sol et la température sont connues. Les estimations régionales de l’affaissement du sol sont compris entre <0.3 à >1.8 cm an
<sup>–1</sup>
. L’incertitude principale est la distribution de la teneur en matière organique du sol qui conduit à un lissage spatial dans la cartographie des taux d’affaissement des sols. L’analyse des données de nivellement et d’extensiométrie dans les rizières conduit à un taux d’accroissement estimé entre 0.02 et 0.8 cm an
<sup>–1</sup>
. Ces valeurs sont généralement en accord avec les estimations indirectes basées sur les flux de carbone et de la minéralisation de l’azote qui indiquent a priori que le riz stoppera ou réduira de manière substantielle l’affaissement du sol. Les zones d’altitude inférieure à –2 m sont des zones candidates pour la mise en œuvre de mesures d’atténuation telles que la culture du riz car il s’y produit un affaissement du sol actif avec des taux supérieurs à 0.4 cm an
<sup>–1</sup>
.</p>
</trans-abstract>
<trans-abstract xml:lang="es" id="Abs3">
<title>Resumen</title>
<p>La subsidencia de los suelos orgánicos en el Delta de Sacramento-San Joaquín amenaza la sostenibilidad del sistema de suministro de agua y la agricultura en California (EE.UU.). Se recolectaron datos de elevación de la superficie del terreno para evaluar las tasas de subsidencia actuales y evaluar el cultivo de arroz como un uso de la tierra para la mitigación de la subsidencia. Se refinó y calibró el modelo SUBCALC desarrollado previamente para describir las tasas de la subsidencia en el ancho actual del Delta, usando los datos más recientes de las emisiones de CO
<sub>2</sub>
y los cambios en la elevación de la superficie del terreno medidos con extensómetros. Se evaluaron los datos de cambio de la elevación de la superficie del terreno con respecto a las estimaciones indirectas de la subsidencia y la acumulación de carbono a partir de datos y el flujo de nitrógeno para el cultivo de arroz. Los datos de extensómetros y de nivelación demuestran variaciones estacionales en las elevaciones de la superficie terrestre asociados a las fluctuaciones del nivel del agua subterránea y tasas de subsidencia inelásticas de 0.5 a 0.8 cm año
<sup>–1</sup>
. La calibración del modelo SUBCALC indicó una precisión de ±0.10 cm año
<sup>–1</sup>
, donde son conocidas la profundidad del agua subterránea, el contenido de materia orgánica del suelo y la temperatura. Las estimaciones regionales de subsidencia varía de <0.3 a >1.8 cm yr
<sup>–1</sup>
. La principal incertidumbre es la distribución del contenido de materia orgánica en el suelo que resulta del promedio espacial en el mapeo de las tasas de subsidencia. El análisis de los datos de nivelación y extensómetro en los campos de arroz resultó en una tasa de acreción estimada entre 0.02 y 0.8 cm yr
<sup>–1</sup>
. Estos valores en general concuerdan con las estimaciones indirectas basadas en los flujos de carbono y nitrógeno, lo que demuestra en forma preliminar que la mineralización del arroz va a detener o reducir en gran medida la subsidencia. Las áreas por debajo de elevaciones de –2 m son áreas examinadas para la implementación de medidas de mitigación, tales como el cultivo de arroz debido a que no hay lugar para un ritmo de subsidencia activa de más de 0.4 cm año
<sup>–1</sup>
.</p>
</trans-abstract>
<trans-abstract xml:lang="zh" id="Abs4">
<title>摘要</title>
<p>沙加缅度河—圣华金河三角洲有机土壤的下沉威胁着(美国)加利佛尼亚州供水系统和农业的可持续性。收集了地面高程资料以评价目前的下沉速率并且评估了土地种植水稻减缓下沉的作用。为了描述三角洲范围内目前的下沉速度,采用目前的CO
<sub>2</sub>
排放资料和伸长仪测量的地面高程变化对先前开发的SUBCALC模型进行了改进和校正。针对利用水稻种植碳和氮通量资料得到的间接下沉和吸积量估算结果,对地表高程变化资料进行了评估。伸长仪和水准测量资料显示了地表高程变化的季节性变化与地下水位波动及0.5 至 0.8 cm yr
<sup>–1</sup>
的非弹性下沉速度有关。在地下水深度、土壤有机物含量和温度已知的情况下,SUBCALC模型校正表明精确度为±0.10 cm yr
<sup>–1</sup>
。区域下沉的估算范围为 < 0.3 至 >1.8 cm yr
<sup>–1</sup>
。主要的不确定性为土壤有机物含量的分布状况,这种不确定性会导致绘制下沉速度时的空间平均值。水稻田中的水准测量和伸长仪资料分析得出的吸积率估算值为0.02 至 0.8 cm yr
<sup>–1</sup>
。这些值通常与根据碳通量和氮矿化作用得到的间接估算值一致,从而初步展示,水稻能阻止或减少下沉。高程低于–2 m的区域是实施减缓措施诸如种植水稻的候选区域,因为在这些区域将会有速度大于0.4 cm yr
<sup>–1</sup>
的下沉。</p>
</trans-abstract>
<trans-abstract xml:lang="pt" id="Abs5">
<title>Resumo</title>
<p>A subsidência dos organossolos no Delta dos rios Sacramento-San Joaquin põe em risco a sustentabilidade do sistema de abastecimento de água e da agricultura do estado da Califórnia (EUA). Foram coletados dados de elevação da superfície do solo para estimar as taxas atuais de subsidência e avaliar o uso de cultivo de arroz para mitigar a subsidência. Para simular as taxas atuais de subsidência da região do Delta, o modelo SUBCALC, desenvolvido anteriormente, foi refinado e calibrado utilizando dados recentes de emissões de CO
<sub>2</sub>
e mudanças na elevação do terreno obtidas com extensômetros. Os dados de mudança na elevação foram avaliados em relação a estimativas indiretas de subsidência e acreção, utilizando dados de fluxo de carbono e nitrogênio no cultivo de arroz. Os dados dos extensômetros e de nivelamento demonstraram variações sazonais de elevações do terreno, associados a flutuações do nível da água subterrânea e taxas de subsidência inelásticas que variam de 0.5 a 0.8 cm ano
<sup>–1</sup>
. A calibração do modelo SUBCALC indicou precisão de ±0.10 cm por ano
<sup>–1</sup>
, onde a profundidade das águas subterrâneas e o teor de matéria orgânica e temperatura do solo eram conhecidos. As estimativas regionais de subsidência variam de <0.3 a >1.8 cm ano
<sup>–1</sup>
. A principal incerteza nas estimativas de subsidência modeladas é a distribuição do conteúdo de matéria orgânica do solo. Porem, as taxas mapeadas de subsidência representam médias em áreas de solos mapeados com teor de materia organica semelhantes. Análises dos dados de nivelamento e dos extensômetros em campos de produção de arroz, resultaram em uma taxa de acreção estimada de 0.02 a 0.8 cm por ano
<sup>–1</sup>
. Estes valores geralmente concordam com estimativas indiretas baseadas em fluxos de carbono e a mineralização do nitrogênio, portanto, indicando que o cultivo de arroz ira parar ou reduzira significativamente o processo de subsidência. Áreas com elevação inferior a –2 m são candidatas para a implementação de medidas de mitigação como o cultivo de arroz, por haver subsidência ativa ocorrendo com taxas superiores a 0.4 cm por ano
<sup>–1</sup>
.</p>
</trans-abstract>
<kwd-group xml:lang="en">
<title>Keywords</title>
<kwd>Subsidence</kwd>
<kwd>Geohazard</kwd>
<kwd>USA</kwd>
<kwd>Land use</kwd>
<kwd>Soil processes</kwd>
</kwd-group>
<funding-group>
<award-group>
<funding-source>
<institution-wrap>
<institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100006196</institution-id>
<institution>Jet Propulsion Laboratory</institution>
</institution-wrap>
</funding-source>
</award-group>
<award-group>
<funding-source>
<institution>U.S. Department of Agriculture</institution>
</funding-source>
<award-id>2011-67003-30371</award-id>
</award-group>
<award-group>
<funding-source>
<institution>California Department of Water Resources</institution>
</funding-source>
<award-id>TW-08-03</award-id>
</award-group>
</funding-group>
<custom-meta-group>
<custom-meta>
<meta-name>issue-copyright-statement</meta-name>
<meta-value>© Springer-Verlag Berlin Heidelberg 2016</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body>
<sec id="Sec1" sec-type="introduction">
<title>Introduction</title>
<p>Worldwide, many deltas are sinking due to reduced aggradation and compaction and soil loss resulting from fluid withdrawal, drainage, and oxidation of organic matter (Syvitski et al.
<xref ref-type="bibr" rid="CR60">2009</xref>
). Drainage of organic soils has resulted in soil subsidence due to changes in physical conditions and enhanced rates of microbial decomposition (Hirano et al.
<xref ref-type="bibr" rid="CR32">2012</xref>
; Rojstaczer and Deverel
<xref ref-type="bibr" rid="CR49">1993</xref>
; Stephens et al.
<xref ref-type="bibr" rid="CR58">1984</xref>
). In total, 14–20 % of the world’s organic soils or peatlands are currently drained for agriculture or forestry (Strack
<xref ref-type="bibr" rid="CR59">2008</xref>
).</p>
<p>Organic soils or histosols are defined by the Food and Agriculture Organization as having 40 cm or more of organic materials with an organic carbon content of 12–18 % or more. The terms organic soil and peat are often used interchangeably and generally refer to a soil which formed under saturated wetland conditions and is acidic and rich in humus. Peat soils cover an estimated 400 million ha, equivalent to 3 % of the Earth’s land surface (Kaat and Joosten
<xref ref-type="bibr" rid="CR37">2009</xref>
), with most peatlands occurring in the northern hemisphere.</p>
<p>Farming and subsidence of peats has been studied in multiple locations. The oldest records are from the Netherlands where peat soils were drained starting between the 9th and 14th centuries. Schorthorst (
<xref ref-type="bibr" rid="CR53">1977</xref>
) documented subsidence rates varying from 0.17 to 0.7 cm yr
<sup>–1</sup>
since the 1800s in these peat soils. Stephens et al. (
<xref ref-type="bibr" rid="CR58">1984</xref>
) summarized worldwide subsidence rates which ranged from less than 0.5 to 10 cm yr
<sup>–1</sup>
in California, Louisiana, Michigan, New York, Indiana and Florida in the USA and the Netherlands, Republic of Ireland, Norway, England (UK), Israel and Russia. Subsidence in agricultural peatlands has also been studied and estimated in New Zealand (Schipper and McLeod
<xref ref-type="bibr" rid="CR52">2002</xref>
), Southeast Asia (e.g. Hooijer et al.
<xref ref-type="bibr" rid="CR34">2012</xref>
) and Italy (Zanello et al.
<xref ref-type="bibr" rid="CR72">2011</xref>
; Gambolati et al.
<xref ref-type="bibr" rid="CR29">2006</xref>
).</p>
<p>Reported causes of peat subsidence include (1) shrinkage due to dewatering, (2) consolidation due to loss of buoyant force and loading, (3) wind and water erosion, (4) oxidation of soil organic matter and (5) burning. Ewing and Vepraskas (
<xref ref-type="bibr" rid="CR25">2006</xref>
) differentiated between relatively larger rates of primary subsidence or shrinkage upon drainage and lower rates of secondary subsidence or oxidation. During the 1970s, Schorthorst (
<xref ref-type="bibr" rid="CR53">1977</xref>
) reported that compaction, shrinkage, and microbial oxidation caused 28, 20, and 52 % of subsidence in the Netherlands, respectively.</p>
<p>US Department of Agriculture and University of Florida researchers extensively studied subsidence in Florida Everglades peats (Stephens et al.
<xref ref-type="bibr" rid="CR58">1984</xref>
) and reported that oxidation accounted for 53 % of historical subsidence. Florida researchers also demonstrated the relation of subsidence and carbon dioxide production (Stephens and Stewart
<xref ref-type="bibr" rid="CR57">1976</xref>
), soil temperature and moisture (Knipling et al.
<xref ref-type="bibr" rid="CR40">1970</xref>
; Volk
<xref ref-type="bibr" rid="CR69">1973</xref>
) and microbial activity (Tate
<xref ref-type="bibr" rid="CR61">1979</xref>
,
<xref ref-type="bibr" rid="CR62">1980a</xref>
,
<xref ref-type="bibr" rid="CR63">1980b</xref>
). Stephens et al. (
<xref ref-type="bibr" rid="CR58">1984</xref>
) and Couwenberg and Hooijer (
<xref ref-type="bibr" rid="CR13">2013</xref>
) reported an inverse correlation of subsidence rates, and oxidative carbon loss, to depth to groundwater.</p>
<sec id="Sec2">
<title>Sacramento-San Joaquin Delta</title>
<p>Substantial understanding and quantification of subsidence has occurred since the early 1900s in the Sacramento-San Joaquin Delta, California, USA (hereafter Delta; Fig. 
<xref rid="Fig1" ref-type="fig">1</xref>
). Subsidence of organic and highly organic mineral soils (hereafter referred to as organic soils or peat) is a primary landscape-altering process that threatens Delta infrastructure and water supply for over 25 million Californians. Drainage and cultivation of Delta soils since 1850 resulted in subsidence on over 60 islands from 1 to over 8 m (Thompson
<xref ref-type="bibr" rid="CR65">1957</xref>
; Deverel and Leighton
<xref ref-type="bibr" rid="CR16">2010</xref>
; Fig. 
<xref rid="Fig2" ref-type="fig">2</xref>
). Key factors influencing subsidence include percent soil organic matter, depth of peat, year of initial drainage and management practices such as burning or growing crops that leave the soil exposed to wind erosion (Deverel and Leighton
<xref ref-type="bibr" rid="CR16">2010</xref>
).
<fig id="Fig1">
<label>Fig. 1</label>
<caption>
<p>Map showing locations of the Sacramento-San Joaquin Delta and subsidence and carbon fluxes measurements</p>
</caption>
<graphic xlink:href="10040_2016_1391_Fig1_HTML" id="MO1"></graphic>
</fig>
<fig id="Fig2">
<label>Fig. 2</label>
<caption>
<p>Map showing land surface elevations where organic soils predominated at the time of data collection in the Sacramento-San Joaquin Delta (based on LiDAR data from CDWR 2007)</p>
</caption>
<graphic xlink:href="10040_2016_1391_Fig2_HTML" id="MO2"></graphic>
</fig>
</p>
<p>The overall objectives of this paper are to: (1) summarize the state of the science for oxidative subsidence of Delta peat soils and report the analysis of recent land-surface elevation data, refinement and recalibration of the subsidence model SUBCALC, (2) estimate present-day subsidence rates and (3) preliminarily assess rice cultivation for stopping or reducing subsidence.</p>
<p>Delta peat soils formed from decaying wetland plants (Atwater
<xref ref-type="bibr" rid="CR1">1982</xref>
; Shlemon and Begg
<xref ref-type="bibr" rid="CR56">1975</xref>
; Drexler et al.
<xref ref-type="bibr" rid="CR23">2009a</xref>
). During the 6,000–7,000 years prior to the 1850s, about 5 billion m
<sup>3</sup>
of tidal marsh sediment accumulated in the Delta (Deverel and Leighton
<xref ref-type="bibr" rid="CR16">2010</xref>
; Mount and Twiss
<xref ref-type="bibr" rid="CR46">2005</xref>
). Since the mid-19th century, half of this volume disappeared (Deverel and Leighton
<xref ref-type="bibr" rid="CR16">2010</xref>
; Mount and Twiss
<xref ref-type="bibr" rid="CR46">2005</xref>
). Present-day soils reflect organic matter accumulation through millennia, spatially variable fluvial deposition and oxidation; thus, soil type and organic matter content vary substantially (Deverel and Leighton
<xref ref-type="bibr" rid="CR16">2010</xref>
). Highly organic mineral surface soils generally predominate in the western and northern Delta and true surface organic soils, or histosols predominate in the central, eastern and southern-central Delta. The lowest organic matter content soils, which subside at relatively low rates, generally predominate in areas drained prior to 1880 near the Sacramento River where there was greater fluvial deposition (Deverel and Leighton
<xref ref-type="bibr" rid="CR16">2010</xref>
). Higher subsidence rates are associated with more recently drained and higher-organic matter soils in the central Delta where there was less fluvial deposition (Deverel and Leighton
<xref ref-type="bibr" rid="CR16">2010</xref>
).</p>
<p>Using data presented in Atwater (
<xref ref-type="bibr" rid="CR1">1982</xref>
) and boring logs, Deverel and Leighton (
<xref ref-type="bibr" rid="CR16">2010</xref>
) and Deverel et al. (
<xref ref-type="bibr" rid="CR22">2015</xref>
) presented maps of peat thickness. The thickest peat resides in the western and northwestern Delta where thicknesses range to over 7 m on Sherman Island (Fig. 
<xref rid="Fig1" ref-type="fig">1</xref>
). Three to over 7 m of peat remains on Ryer, southern Grand, western Brannan, and Twitchell islands (Fig. 
<xref rid="Fig1" ref-type="fig">1</xref>
). For most of the central, eastern and southern-central Delta, less than 1–2 m of peat remains (see Fig. 
<xref rid="Fig4" ref-type="fig">4</xref>
in Deverel et al.
<xref ref-type="bibr" rid="CR22">2015</xref>
).</p>
<p>Deverel and Rojstaczer (
<xref ref-type="bibr" rid="CR17">1996</xref>
) reported that microbial oxidation of organic matter is the primary present-day cause of subsidence. Consistently, Rojstaczer and Deverel (
<xref ref-type="bibr" rid="CR50">1995</xref>
) and Deverel and Leighton (
<xref ref-type="bibr" rid="CR16">2010</xref>
) demonstrated that spatial variations in soil organic matter content ranged from 4 to 60 % and explained over 55 % of the variation in average subsidence rates from 1978 to 2006. Deverel and Leighton (
<xref ref-type="bibr" rid="CR16">2010</xref>
) assessed recent and historic causes of and factors affecting subsidence rates using elevation and soils data collected during 2006 and data reported in Rojstaczer et al. (
<xref ref-type="bibr" rid="CR51">1991</xref>
), Rojstaczer and Deverel (
<xref ref-type="bibr" rid="CR50">1995</xref>
), Deverel and Rojstaczer (
<xref ref-type="bibr" rid="CR17">1996</xref>
) and the University of California (Weir
<xref ref-type="bibr" rid="CR70">1950</xref>
). Deverel and Leighton (
<xref ref-type="bibr" rid="CR16">2010</xref>
) developed a computer model, SUBCALC, and geographic information system (GIS) to simulate Delta-wide subsidence. The SUBCALC model was calibrated using historic land-surface elevation data collected on three islands in the central Delta (Bacon and Mildred islands and Jones Tract) by Weir (
<xref ref-type="bibr" rid="CR70">1950</xref>
) and Sherman Island in the western Delta (Rojstaczer et al.
<xref ref-type="bibr" rid="CR51">1991</xref>
).</p>
<p>Subsidence rates have decreased with time associated with decreasing soil organic matter content and changing land-management practices (Deverel and Leighton
<xref ref-type="bibr" rid="CR16">2010</xref>
). Prior to the early 1960s, burning and wind erosion caused soil loss. Burning no longer occurs, and there is minimal wind erosion. Wind erosion was associated with high-velocity spring winds and bare asparagus fields (Schultz and Carlton
<xref ref-type="bibr" rid="CR54">1959</xref>
; Schultz et al.
<xref ref-type="bibr" rid="CR55">1963</xref>
; Carlton and Schultz
<xref ref-type="bibr" rid="CR10">1966</xref>
). Asparagus was widely planted in the Delta during the 1920s through the 1950s. Due to economic reasons, since the 1960s, asparagus cultivation decreased to a small area. Model results and their agreement with measurements demonstrate that oxidation accounts for the majority of recent subsidence and the remaining portion is due to consolidation (Deverel and Leighton
<xref ref-type="bibr" rid="CR16">2010</xref>
). The application of Michaelis–Menton kinetics was used to simulate oxidation of soil organic carbon. Consistently, Volk (
<xref ref-type="bibr" rid="CR69">1973</xref>
) reported that Michaelis–Menton kinetics appropriately described oxidation of peat soils in the Florida Everglades.</p>
<p>Little has been documented about consolidation, a secondary cause of organic soil subsidence. Water in organic soils is held in three phases: (1) intercellular, (2) inter-particle water in micropores, and (3) bound or absorbed. Consolidation expulses pore water and particles rearrange (Hobbs
<xref ref-type="bibr" rid="CR33">1986</xref>
). As farmers deepened drainage ditches to compensate for land-surface elevation loss due to oxidation, wind erosion and burning, organic soils consolidated due to dewatering resulting from increased drainage-ditch depth, which reduces pore pressure and buoyancy, thus transferring load to the soil skeleton. Drexler et al. (
<xref ref-type="bibr" rid="CR24">2009b</xref>
) presented evidence for consolidation below the upper oxidized layer on farmed subsided islands.</p>
<p>To estimate consolidation of subsurface deposits, Deverel and Leighton (
<xref ref-type="bibr" rid="CR16">2010</xref>
) assumed compaction proceeds similar to dewatering and irreversible consolidation of a vertical soil column as described by Terzaghi (
<xref ref-type="bibr" rid="CR64">1925</xref>
). They employed Terzaghi’s effective stress principle and extensometer data on Twitchell Island to estimate consolidation in SUBCALC using a linear model relating compaction to the change in hydraulic head. The estimated percentages of the different causes of subsidence came from model simulations of elevation change data from the 1920s to 2006.</p>
<p>Since the late 1990s, researchers have investigated large-scale and small-scale spatial and temporal trends of subsidence on Delta islands and levees using remotely sensed data. Attempts to use satellite radar data for investigating oxidative subsidence of organic soils suffered from rapid decorrelation in agricultural areas (Cohen et al.
<xref ref-type="bibr" rid="CR11">1998</xref>
; Brooks et al.
<xref ref-type="bibr" rid="CR4">2012</xref>
)—for example, Brooks et al. (
<xref ref-type="bibr" rid="CR4">2012</xref>
) showed limited ability to estimate subsidence on Delta island interiors (see Fig. 
<xref rid="Fig2" ref-type="fig">2</xref>
in their paper) presumably due to decorrelation. The Brooks et al. (
<xref ref-type="bibr" rid="CR4">2012</xref>
) estimates were primarily for Delta levees around the periphery of islands where decorrelation has less of an effect. The applicability of satellite InSar is primarily limited to more stable structures such as levees.</p>
<p>More recently, the relatively longer L-band (23.8 cm) wavelength of the UAVSAR, combined with regular acquisitions, high spatial resolution and data processing techniques developed for low coherence regions, show promise for application of radar interferometry to monitor subsidence in Delta agricultural lands on island interiors at sub-centimeter vertical resolution levels on organic soils (Jones et al.
<xref ref-type="bibr" rid="CR35">2011</xref>
,
<xref ref-type="bibr" rid="CR36">2012</xref>
). UAVSAR data have been used to estimate rates of elevation change on Sherman Island organic soils (Priyanka Sharma, Jet Propulsion Laboratory, unpublished data, 2015). Values ranged from 0 to 5 cm yr
<sup>–1</sup>
and the spatially averaged rate was 1.5 cm yr
<sup>–1</sup>
. The UAVSAR-estimated rates were generally consistent with ground-based estimates described here and Deverel and Leighton (
<xref ref-type="bibr" rid="CR16">2010</xref>
) and Rojstaczer and Deverel (
<xref ref-type="bibr" rid="CR50">1995</xref>
). Terrestrial light detection and ranging technology (Terrestrial LiDAR) has also been used in the Delta to quantify short-term local-scale levee deformation (e.g. Bawden et al.
<xref ref-type="bibr" rid="CR2">2014</xref>
).</p>
<sec id="Sec3">
<title>Effects of subsidence and mitigation</title>
<p>By reducing the landmass and resistance to hydraulic pressure from adjacent channels, subsidence has contributed to levee failure and island inundation. Weir (
<xref ref-type="bibr" rid="CR70">1950</xref>
) reported the results of land surveys and field observations on three central Delta islands from 1922 to 1948 where elevations declined 1.8–2.3 m, resulting in elevations ranging from 3 to 3.4 m below sea level. Weir (
<xref ref-type="bibr" rid="CR70">1950</xref>
) warned about increased seepage, levee instability and island flooding associated with subsidence. At the time of Weir’s publication, 50 islands had flooded primarily due to overtopping of levees. From 1930 to the early 1980s, over 50 Delta islands or tracts flooded—the majority due to levee foundation instability (Prokopovitch
<xref ref-type="bibr" rid="CR48">1985</xref>
)—whereas from 1900 to 2006, over 100 island levees failed (Gaddie et al.
<xref ref-type="bibr" rid="CR28">2006</xref>
; Florsheim and Dettinger
<xref ref-type="bibr" rid="CR27">2007</xref>
), causing local infrastructural damage which historically cost hundreds of millions of dollars (Prokopovitch
<xref ref-type="bibr" rid="CR48">1985</xref>
).</p>
<p>The flooding of the 4,860-ha Jones Tract due to a levee breach that occurred during the early morning of June 2004 illustrates recent consequences. About $90 million were expended by the State of California to close the levee breach and complete removal of 197 million m
<sup>3</sup>
of flood water by December 2004. Before island reclamation was complete, litigation began when BNSF Railway alleged that the State of California’s operation of the State Water Project resulted in channel scour that induced levee failure. Jones Tract landowners and others joined in the lawsuit. In its statement of decision and judgement, the State Superior Court ruled against the plaintiffs, dismissed the theory of State-Water-Project channel scour and stated that: “Subsidence resulting from the loss of peat soil over the years makes Delta islands more susceptible to flooding due to levee failure”.</p>
<p>Island flooding in the western Delta may also cause eastward movement of saline water into the Delta, and thus impede water exports—for example, levee failure and flooding on Brannan-Andrus Island (Fig. 
<xref rid="Fig1" ref-type="fig">1</xref>
) in June 1972 caused movement of salt water into the Delta (Cook and Coleman
<xref ref-type="bibr" rid="CR12">1973</xref>
), resulting in cessation of state and federal water exports. An additional 370 million m
<sup>3</sup>
of water was released from reservoirs to mitigate the salinity intrusion. The total cost of flooding was over $97 million in 2015 dollars.</p>
<p>Delta subsidence will continue until management practices are adopted that stop subsidence or the organic deposits disappear. Continuing subsidence can make farming more difficult and expensive. For example, as peat disappears, drainage ditches may be excavated into underlying mineral sediments which can be unstable. Also, greater seepage due to increased hydraulic gradients and flow under levees is resulting in larger marginally or non-farmable acreage (Deverel et al.
<xref ref-type="bibr" rid="CR22">2015</xref>
).</p>
<p>Hydraulic forces on and seepage through and under levees will increase with continuing subsidence and sea level rise (Deverel et al.
<xref ref-type="bibr" rid="CR19">2007a</xref>
,
<xref ref-type="bibr" rid="CR21">2014</xref>
,
<xref ref-type="bibr" rid="CR22">2015</xref>
). Increased seepage under levees onto islands will decrease levee stability. Deverel et al. (
<xref ref-type="bibr" rid="CR20">2007b</xref>
) predicted that seepage onto Twitchell Island (Fig. 
<xref rid="Fig1" ref-type="fig">1</xref>
) will increase by 22–34 % during the next several decades. Levee failure or instability from seepage occurs when hydraulic gradients are large enough to erode or move levee internal and/or foundational materials. The Delta Risk Management Strategy (CDWR
<xref ref-type="bibr" rid="CR8">2009</xref>
) estimated the combined probability of levee failure and island flooding from earthquake, high-water flooding, and sunny-day levee failure for most of the deeply subsided central and western Delta ranges from 53 % to over 84 % during the next 20 years. Economic costs of future Delta levee failures were estimated in billions of dollars.</p>
<p>Island drainage volumes and dissolved organic carbon loads will increase with continuing subsidence (Deverel et al.
<xref ref-type="bibr" rid="CR19">2007a</xref>
) and drainage costs will increase due to greater pumping lifts and volumes. Oxidation of organic soils generates dissolved organic carbon which forms, upon treatment for drinking, carcinogenic disinfection byproducts (Fleck et al.
<xref ref-type="bibr" rid="CR26">2004</xref>
; Deverel et al.
<xref ref-type="bibr" rid="CR20">2007b</xref>
). Winter and spring rains and irrigation mobilize dissolved organic carbon and trihalomethane precursor loads into island drainage water (Deverel et al.
<xref ref-type="bibr" rid="CR20">2007b</xref>
). Ongoing oxidative subsidence therefore perpetuates the annual cycle of generation of dissolved organic carbon in soil and mobilization to drainage water which is exported to Delta channels which deliver drinking water to 25 million Californians via the State Water Project (CDWR
<xref ref-type="bibr" rid="CR9">2015</xref>
). Similar processes apparently operate to generate methyl mercury, and data presented in Heim et al. (
<xref ref-type="bibr" rid="CR31">2009</xref>
) indicate that methyl mercury loads from Delta islands will increase with continuing subsidence.</p>
<p>Managed and permanently flooded wetlands will stop and reverse the effects of subsidence (Deverel et al.
<xref ref-type="bibr" rid="CR18">1998</xref>
,
<xref ref-type="bibr" rid="CR21">2014</xref>
; Miller et al.
<xref ref-type="bibr" rid="CR44">2000</xref>
,
<xref ref-type="bibr" rid="CR45">2008</xref>
). These wetlands accumulate carbon, and accretion rates are about 3 cm yr
<sup>–1</sup>
(Deverel et al.
<xref ref-type="bibr" rid="CR21">2014</xref>
); wetlands near levees will reduce seepage onto islands (Deverel et al.
<xref ref-type="bibr" rid="CR21">2014</xref>
).</p>
<p>Agriculture is the predominant land use in the Delta and is important to the local economy. However, current farming practices which require an aerated root zone cause subsidence by exposing organic soils to oxygen; therefore, agricultural practices that stop subsidence are highly desirable. The original impetus for investigating rice as a subsidence mitigation land use came from Miller et al. (
<xref ref-type="bibr" rid="CR44">2000</xref>
) who indicated that wetlands which were flooded from early spring through midsummer resulted in no net carbon loss. Rice growers use a similar water management practice, flooding rice fields during the warmest months when soil oxidation rates are highest. In the past, cool night temperatures precluded Delta rice cultivation; however, development of cold-tolerant rice varieties resulted in increased Delta rice production with yields generally comparable to other rice-growing areas in California. Rice has been successfully grown on about 3,000 acres on central and eastern Delta islands since the mid-1990s.</p>
<p>Heightened recent interest in subsidence mitigation prompted further investigation into rice production on state-owned Twitchell Island in 2009. Micrometeorological data presented by Hatala et al. (
<xref ref-type="bibr" rid="CR30">2012</xref>
) and soil nitrogen dynamics reported by Kirk et al. (
<xref ref-type="bibr" rid="CR39">2015</xref>
) suggest that rice cultivation will greatly reduce oxidative subsidence in Delta organic soils. Extensometer and leveling data collected in rice fields and an adjacent cornfield are reported here.</p>
</sec>
<sec id="Sec4">
<title>Need for dissemination of quality information</title>
<p>Misinformation about present-day Delta subsidence points to the need for collection and dissemination of high quality information. For example, in the Economic Sustainability Plan, the Delta Protection Commission (
<xref ref-type="bibr" rid="CR15">2012</xref>
) cited work which attempted to delineate areas of active subsidence based on comparison of 2007 LiDAR data (CDWR
<xref ref-type="bibr" rid="CR7">2007</xref>
) and USGS Quadrangle maps surveyed between 1974 and 1976 (California Central Valley Flood Control Association
<xref ref-type="bibr" rid="CR6">2011</xref>
). The accuracy of the 2007 LiDAR data is about ±0.15 m. The error due to estimating the elevations from the 1974 and 1976 quadrangle contour maps is about one-half of the contour interval (1.5 m) for the topographic maps or 0.76 m (J. Vukovitch, USGS, Denver, personal communication, 1996). During 1974 to 2007, subsidence rates ranged from about 1–3 cm per yr
<sup>–1</sup>
, resulting in 0.3–0.9 m of elevation change which is similar to the estimation error from the Quadrangle maps and LiDAR data.</p>
<p>Since publication of Deverel and Leighton (
<xref ref-type="bibr" rid="CR16">2010</xref>
), additional land-surface elevation-change data has been collected on selected Delta islands; also, greenhouse gas emissions have been measured. While there is a need to assess, analyze and synthesize these data and improve subsidence modeling capability, the overall approach taken here was to collect and analyze land-surface elevation data in rice, corn and pasture fields, refine and recalibrate SUBCALC and estimate present-day Delta subsidence rates.</p>
</sec>
</sec>
</sec>
<sec id="Sec5" sec-type="materials|methods">
<title>Data sources and methods</title>
<sec id="Sec6">
<title>Land-surface elevation data</title>
<sec id="Sec7">
<title>Extensometer data</title>
<p>To monitor small-scale variations in land surface elevation during 2009–2015, HydroFocus personnel installed and operated two extensometers on Twitchell Island, one in a rice field and one in the cornfield adjacent to the rice field (Fig. 
<xref rid="Fig1" ref-type="fig">1</xref>
). The cornfield was converted to a wetland in 2014. At the extensometer location on Sherman Island (Fig. 
<xref rid="Fig1" ref-type="fig">1</xref>
) described in Deverel and Rojstaczer (
<xref ref-type="bibr" rid="CR17">1996</xref>
), new instrumentation recorded elevation changes starting in 2011. At all locations, land-surface elevations were measured relative to the extensometer structure which was anchored below the peat.</p>
<p>In the rice field, a steel base support pipe was driven to refusal into the mineral layer underlying the peat soil. A modified sedimentation-erosion table (SET; Boumans and Day
<xref ref-type="bibr" rid="CR3">1993</xref>
) was inserted into a grooved stainless steel sleeve in the base support pipe which ensured instrument stability and replacement to the exact same position after movement to accommodate field operations. The SET arm extended horizontally about 1 m and was adjusted to level. A metal rod with a 5-inch-diameter (12.7 cm) metal disk that rested on the ground freely moved vertically in a sleeve on a metal plate at the end of the arm. HydroFocus personnel fastened a Macro Sensors GHSI 750 linear variable differential transformer (LVDT) to the rod above the plate. The piston arm rested on the plate so that the sensor body would move with the rod and the piston arm would remain stationary. A Campbell CR510 data logger recorded LVDT measurements every 15 min. In the cornfield, HydroFocus personnel constructed an extensometer similar to the one described in Deverel and Rojstaczer (
<xref ref-type="bibr" rid="CR17">1996</xref>
).</p>
</sec>
<sec id="Sec8">
<title>Leveling data</title>
<p>Annual leveling surveys were conducted at seven locations in the Twitchell Island rice fields during 2009–2013 (Fig. 
<xref rid="Fig1" ref-type="fig">1</xref>
) by California Department of Water Resources personnel. Spirit leveling surveys in multiple directions relative to fixed monuments anchored in the mineral material were conducted in the spring of each year after cultivation and before flooding. Land-surface elevations were determined using GPS in 2001 and 2012 at seven monitoring-well locations on Twitchell Island (Fig. 
<xref rid="Fig1" ref-type="fig">1</xref>
).</p>
</sec>
</sec>
<sec id="Sec9">
<title>Modeling of delta subsidence</title>
<p>Deverel and Leighton (
<xref ref-type="bibr" rid="CR16">2010</xref>
) developed the computer model, SUBCALC, to integrate available data and quantify and simulate subsidence rates and causes. SUBCALC simulates aerobic microbial oxidation of organic carbon, consolidation, wind erosion and burning. Present-day subsidence is the result of oxidation and consolidation. SUBCALC simulates microbial oxidation of soil organic carbon to carbon dioxide using Michaelis–Menton (M–M) enzyme kinetics in which the rate of soil organic-matter oxidation is limited by soil organic carbon content (Browder and Volk
<xref ref-type="bibr" rid="CR5">1978</xref>
):
<disp-formula id="Equ1">
<label>1</label>
<alternatives>
<tex-math id="M1">\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \frac{V}{V_{\max }}=\frac{\left[S\right]}{K_{\mathrm{m}}+\left[S\right]} $$\end{document}</tex-math>
<mml:math id="M2">
<mml:mfrac>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mtext>max</mml:mtext>
</mml:msub>
</mml:mfrac>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mfenced close="]" open="[">
<mml:mi>S</mml:mi>
</mml:mfenced>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi mathvariant="normal">m</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mfenced close="]" open="[">
<mml:mi>S</mml:mi>
</mml:mfenced>
</mml:mrow>
</mml:mfrac>
</mml:math>
<graphic xlink:href="10040_2016_1391_Article_Equ1.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
</p>
<p>Michaelis-Menton equation (Eq. 
<xref rid="Equ1" ref-type="">1</xref>
) parameters,
<italic>K</italic>
<sub>m</sub>
(the M–M constant) and
<italic>V</italic>
<sub>max</sub>
, (the maximum oxidation rate), and effects of temperature were originally estimated based on data reported in Deverel and Rojstaczer (
<xref ref-type="bibr" rid="CR17">1996</xref>
). The [
<italic>S</italic>
] term (substrate concentration term) is the soil organic carbon fraction. For each annual time step, the different contributions to subsidence (e.g. oxidation and consolidation) were estimated based on newly calculated mass of organic matter and bulk densities (Deverel and Leighton
<xref ref-type="bibr" rid="CR16">2010</xref>
).</p>
<p>To estimate the consolidation of subsurface deposits, it was assumed that compaction processes are similar to dewatering and irreversible consolidation of a vertical soil column as described by Terzaghi (
<xref ref-type="bibr" rid="CR64">1925</xref>
). The use of Terzaghi’s effective stress principle is generally restricted by assumptions of Newtonian behavior of the liquid phase. Water in organic soil does not strictly follow Newtonian mechanical principles, especially during large changes in stress; however, it was assumed that for a small increment of stress change, dewatering would generally follow Newtonian behavior. SUBCALC simulates this process using a linear equation relating compaction to the change in hydraulic head based on data from the Twitchell Island extensometer (Kerr et al.
<xref ref-type="bibr" rid="CR38">2003</xref>
). Effects of varying depth to groundwater were accounted for using the relation of subsidence rates to depth to groundwater described in Stephens et al. (
<xref ref-type="bibr" rid="CR58">1984</xref>
). Substantial detail is provided for the original model in Deverel and Leighton (
<xref ref-type="bibr" rid="CR16">2010</xref>
, see their Appendix B in their ‘Supplemental materials’ section).</p>
<sec id="Sec10">
<title>Recent model updates</title>
<p>In light of recently available data for land-surface elevation changes and greenhouse-gas emissions from drained Delta organic soils, SUBCALC was modified and re-calibrated for site specific data (Table
<xref rid="Tab1" ref-type="table">1</xref>
). Specifically, data in Table
<xref rid="Tab1" ref-type="table">1</xref>
and information presented in Davidson et al. (
<xref ref-type="bibr" rid="CR14">2012</xref>
) were used to estimate
<italic>V</italic>
<sub>max</sub>
for Eq. 
<xref rid="Equ2" ref-type="">2</xref>
as follows.
<table-wrap id="Tab1">
<label>Table 1</label>
<caption>
<p>Observed versus SUBCALC-simulated subsidence and carbon flux rates</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th rowspan="2">Site (source)</th>
<th colspan="2">Subsidence
<sup>a</sup>
<break></break>
(cm yr
<sup>–1</sup>
)</th>
<th colspan="2">Carbon Flux
<sup>b</sup>
<break></break>
(g C cm
<sup>–2</sup>
yr
<sup>–1</sup>
)</th>
<th rowspan="2">Soil organic matter content (%)</th>
<th rowspan="2">Average annual depth to groundwater
<break></break>
(cm)</th>
<th rowspan="2">Average annual soil temperature
<break></break>
(°C)</th>
<th rowspan="2">
<italic>V</italic>
<sub>max</sub>
<break></break>
(g C cm
<sup>–2</sup>
yr
<sup>–1</sup>
)</th>
<th rowspan="2">
<italic>K</italic>
<sub>m</sub>
</th>
</tr>
<tr>
<th>Observed</th>
<th>Simulated</th>
<th>Observed</th>
<th>Simulated</th>
</tr>
</thead>
<tbody>
<tr>
<td>Twitchell (this study and Knox et al.
<xref ref-type="bibr" rid="CR41">2015</xref>
)</td>
<td>0.830</td>
<td>0.802</td>
<td>0.057</td>
<td>0.0607</td>
<td char="." align="char">39.7</td>
<td>85</td>
<td char="." align="char">16.4</td>
<td char="." align="char">0.182</td>
<td char="." align="char">0.076</td>
</tr>
<tr>
<td>Sherman (this study)</td>
<td>0.515</td>
<td>0.330</td>
<td></td>
<td></td>
<td char="." align="char">16.9</td>
<td>55</td>
<td char="." align="char">16.0</td>
<td char="." align="char">0.174</td>
<td char="." align="char">
<italic>0.128</italic>
</td>
</tr>
<tr>
<td>Sherman (Hatala et al.
<xref ref-type="bibr" rid="CR30">2012</xref>
)</td>
<td></td>
<td></td>
<td>0.024</td>
<td>0.023</td>
<td char="." align="char">22.5</td>
<td>50</td>
<td char="." align="char">16.0</td>
<td char="." align="char">0.174</td>
<td char="." align="char">0.115</td>
</tr>
<tr>
<td>Staten (Pellerin et al.
<xref ref-type="bibr" rid="CR47">2013</xref>
)</td>
<td></td>
<td></td>
<td>0.066</td>
<td>0.056</td>
<td char="." align="char">34.0</td>
<td>100</td>
<td char="." align="char">15.0</td>
<td char="." align="char">0.157</td>
<td char="." align="char">0.089</td>
</tr>
<tr>
<td>Orwood (Deverel and Rojstaczer
<xref ref-type="bibr" rid="CR17">1996</xref>
)</td>
<td>0.800</td>
<td>0.858</td>
<td></td>
<td></td>
<td char="." align="char">24.4</td>
<td>139</td>
<td char="." align="char">14.9</td>
<td char="." align="char">0.155</td>
<td char="." align="char">
<italic>0.111</italic>
</td>
</tr>
<tr>
<td>Jersey (Deverel and Rojstaczer
<xref ref-type="bibr" rid="CR17">1996</xref>
)</td>
<td>0.680</td>
<td>0.648</td>
<td></td>
<td></td>
<td char="." align="char">20.0</td>
<td>108</td>
<td char="." align="char">15.4</td>
<td char="." align="char">0.164</td>
<td char="." align="char">
<italic>0.121</italic>
</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p>Site-specific
<italic>K</italic>
<sub>
<italic>m</italic>
</sub>
values were used to determine a regression equation for model input km as a function of soil organic matter content.
<italic>Non-italicized K</italic>
<sub>
<italic>m</italic>
</sub>
<italic>values</italic>
are from sites where carbon flux was measured, so
<italic>K</italic>
<sub>
<italic>m</italic>
</sub>
was calculated directly (see Eqs.
<xref rid="Equ1" ref-type="">1</xref>
and
<xref rid="Equ2" ref-type="">2</xref>
).
<italic>Italicized K</italic>
<sub>
<italic>m</italic>
</sub>
<italic>values</italic>
are from sites where carbon flux was not measured. These were calibrated to field subsidence data through manual trial-and-error SUBCALC simulations.</p>
<p>
<sup>a</sup>
For subsidence (observed and simulated), the root mean square error (RMSE) is 0.0991</p>
<p>
<sup>b</sup>
For carbon flux (observed and simulated), the RMSE value is 0.00605</p>
</table-wrap-foot>
</table-wrap>
<disp-formula id="Equ2">
<label>2</label>
<alternatives>
<tex-math id="M3">\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {V}_{\max_{\mathrm{x}}}={a}_{\mathrm{x}}\times {e}^{-\mathrm{E}{\mathrm{a}}_{\mathrm{x}}/RT} $$\end{document}</tex-math>
<mml:math id="M4">
<mml:msub>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mo>max</mml:mo>
<mml:mi mathvariant="normal">x</mml:mi>
</mml:msub>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi mathvariant="normal">x</mml:mi>
</mml:msub>
<mml:mo>×</mml:mo>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo></mml:mo>
<mml:mi mathvariant="normal">E</mml:mi>
<mml:msub>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mi mathvariant="normal">x</mml:mi>
</mml:msub>
<mml:mo stretchy="false">/</mml:mo>
<mml:mi>R</mml:mi>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msup>
</mml:math>
<graphic xlink:href="10040_2016_1391_Article_Equ2.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
</p>
<p>Where
<italic>a</italic>
<sub>x</sub>
is the pre-exponential factor, Ea
<sub>x</sub>
is the activation energy for the soil organic carbon oxidation reaction,
<italic>T</italic>
is soil temperature and
<italic>R</italic>
is the universal gas constant. Values for
<italic>a</italic>
<sub>x</sub>
and Ea
<sub>x</sub>
were initially extracted from Davidson et al. (
<xref ref-type="bibr" rid="CR14">2012</xref>
). Site recorded soil temperatures shown in Table
<xref rid="Tab1" ref-type="table">1</xref>
were used.</p>
<p>Using the fraction organic carbon values and calculated
<italic>V</italic>
<sub>max</sub>
values,
<italic>K</italic>
<sub>m</sub>
(Eq. 
<xref rid="Equ1" ref-type="">1</xref>
) was used as a calibration term to match carbon fluxes and subsidence rates in Table
<xref rid="Tab1" ref-type="table">1</xref>
. A linear relationship between calculated/calibrated
<italic>K</italic>
<sub>m</sub>
values and soil organic matter content arose. The regression equation was used to estimate
<italic>K</italic>
<sub>m</sub>
from soil organic matter content values for estimating Delta-wide subsidence rates.</p>
</sec>
</sec>
<sec id="Sec11">
<title>Estimation of current subsidence rates</title>
<p>The primary spatially variable inputs for SUBCALC are depth to groundwater, soil temperature and soil organic matter content. Michaelis–Menton inputs are calculated from soil temperature and soil organic matter content. The depth of the organic soil where oxidation is simulated to occur is determined by the depth to groundwater and the oxidation rate is governed by depth-to-groundwater/carbon loss-subsidence relations described in Stephens et al. (
<xref ref-type="bibr" rid="CR58">1984</xref>
) and Couwenberg and Hooijer (
<xref ref-type="bibr" rid="CR13">2013</xref>
). The depth to groundwater was estimated from soil surveys described in the following.</p>
<sec id="Sec12">
<title>Depth to groundwater</title>
<p>Depth to groundwater on Delta subsided islands is controlled primarily by networks of drainage ditches that feed to island drainage pumping stations that in turn continuously discharge drainage water to Delta channels. Drainage ditches collect water that seeps from adjacent channels and deep percolation of applied irrigation water. There are few depth to groundwater measurements in Delta organic soils and, in general, groundwater levels have been maintained at about 0.8–1.2 m below land surface as the result of drainage system operation. Based on the first author’s experience in working in the organic soils throughout the Delta since the early 1980s, depth to groundwater has not changed substantially over time in most places. Also, data presented in Deverel et al. (
<xref ref-type="bibr" rid="CR22">2015</xref>
) indicate lack of change in Delta groundwater levels since the late 1980s. To estimate depth to groundwater throughout the Delta for input to the SUBCALC model, information was obtained for each soil type from the soil surveys for Sacramento, San Joaquin, Solano, Yolo, and Contra Costa counties. The soil surveys have an average, or range, of depth to water value for each soil type. Depth to groundwater values were incorporated into a GIS file used to generate a map of estimated depth to groundwater (USDA Soil Survey Staff, Natural Resources Conservation Service
<xref ref-type="bibr" rid="CR67">2006</xref>
,
<xref ref-type="bibr" rid="CR68">2007</xref>
; Welch
<xref ref-type="bibr" rid="CR71">1977</xref>
; McElhiney
<xref ref-type="bibr" rid="CR43">1992</xref>
; Tugel
<xref ref-type="bibr" rid="CR66">1993</xref>
).</p>
</sec>
<sec id="Sec13">
<title>Soil organic matter content and bulk density</title>
<p>Soil organic matter content percentages provided in soil surveys were modified with the results for recently collected soil samples on Twitchell, Staten, Bacon and Sherman islands. Due to oxidation of soil organic matter since collection of data for the soil surveys, available data indicate that present-day soil organic matter content is likely equal to or lower than the mid-range values used to map soil organic matter content in Deverel and Leighton (
<xref ref-type="bibr" rid="CR16">2010</xref>
). Soil organic matter determinations on Twitchell and Staten islands during 2012 through 2014 and on Bacon and Sherman islands in 2006 were compared with values reported in the soil surveys and a regression relation was used to estimate present-day values (Figs. 
<xref rid="Fig3" ref-type="fig">3</xref>
and
<xref rid="Fig4" ref-type="fig">4</xref>
). Data presented in Drexler et al. (
<xref ref-type="bibr" rid="CR24">2009b</xref>
) for soil percent organic matter and bulk density (Fig. 
<xref rid="Fig5" ref-type="fig">5</xref>
) were used to develop a regression relation for estimating initial soil bulk density in SUBCALC.
<fig id="Fig3">
<label>Fig. 3</label>
<caption>
<p>Scatter plot showing relation of average measured soil organic matter content values to middle of range of reported values in the soil surveys</p>
</caption>
<graphic xlink:href="10040_2016_1391_Fig3_HTML" id="MO3"></graphic>
</fig>
<fig id="Fig4">
<label>Fig. 4</label>
<caption>
<p>Scatter plot showing relation of measured soil organic matter content values to middle of range of reported values in the soil surveys</p>
</caption>
<graphic xlink:href="10040_2016_1391_Fig4_HTML" id="MO4"></graphic>
</fig>
<fig id="Fig5">
<label>Fig. 5</label>
<caption>
<p>Scatter plot showing the relation of soil bulk density and soil organic matter content on agricultural islands reported in Drexler et al. (
<xref ref-type="bibr" rid="CR24">2009b</xref>
) for samples collected above and below 120 cm</p>
</caption>
<graphic xlink:href="10040_2016_1391_Fig5_HTML" id="MO5"></graphic>
</fig>
</p>
</sec>
<sec id="Sec14">
<title>Soil temperature</title>
<p>Spatially referenced monthly near-surface air temperature data were obtained from sources described by Maurer et al. (
<xref ref-type="bibr" rid="CR42">2002</xref>
). Data included daily minimum and maximum temperature observations at National Oceanic and Atmospheric Administration (NOAA) Cooperative Observer (co-op) stations. Co-op stations occur at a density of approximately one per 700 km
<sup>2</sup>
. For the Maurer dataset, these observations were spatially interpolated into a 1/8 ° (degrees latitude-longitude) square grid and averaged monthly.</p>
<p>For each mapped grid cell overlaying Delta organic soils, an annual average temperature for the period 2007–2010 was calculated from all monthly averages. This period was chosen because the 1950–1999 data exhibit an upward trend. It was assumed that the 2007–2010 data adequately represent present-day Delta temperatures. Each soil feature in the GIS shapefile from this study had a temperature assigned to it from its corresponding grid cell. Air and soil temperature data from Knox et al. (
<xref ref-type="bibr" rid="CR41">2015</xref>
) and Hatala et al. (
<xref ref-type="bibr" rid="CR30">2012</xref>
) were compared to adjust the air temperatures for SUBCALC model input. On average, soil temperature was 1.2 °C greater than air temperature. This difference was used to convert air to soil temperature for input into Arrhenius function calculations within SUBCALC (Eq. 
<xref rid="Equ2" ref-type="">2</xref>
).</p>
</sec>
<sec id="Sec15">
<title>Mapping of present-day subsidence</title>
<p>Spatially variable present-day subsidence rates were estimated using the recalibrated SUBCALC model and ArcGIS Spatial Analyst. It was assumed that oxidation and consolidation are the only present-day causes of subsidence. It was also assumed that there will be zero subsidence in rice-growing areas and permanently flooded wetlands, and the subsidence rate is zero where the soil organic matter content is less than or equal to 2 %.</p>
</sec>
</sec>
</sec>
<sec id="Sec16" sec-type="results">
<title>Results</title>
<sec id="Sec17">
<title>Recently measured and estimated subsidence rates</title>
<sec id="Sec18">
<title>Extensometer and leveling data</title>
<p>Twitchell Island cornfield extensometer and observation-well data (Fig. 
<xref rid="Fig6" ref-type="fig">6</xref>
) illustrate seasonal variations in land-surface elevations associated with groundwater-level fluctuations from 2009 through 2013. Groundwater levels rose during fall primarily due to decreased crop evapotranspiration and winter precipitation recharge. Groundwater levels decreased in the spring with diminishing rain and increasing evapotranspiration. Using land-surface elevation measurements at times of equal groundwater levels in October 2009 and December 2013, an average inelastic subsidence rate of 0.83 cm yr
<sup>–1</sup>
was estimated. The soil organic matter content was 39.7 % (Table
<xref rid="Tab1" ref-type="table">1</xref>
). For the seven locations where elevations were determined at observation wells in 2001 and 2012, subsidence rates ranged from 0.11 to 1.94 cm yr
<sup>–1</sup>
in agricultural fields where estimated soil organic matter content ranged from 6 to 20 %. Generally consistent with extensometer data, the average of all seven measurements during 2001–2012 was 0.7 cm yr
<sup>–1</sup>
.
<fig id="Fig6">
<label>Fig. 6</label>
<caption>
<p>Time-series plot of land-surface elevation changes and groundwater level, Twitchell Island extensometer in cornfield</p>
</caption>
<graphic xlink:href="10040_2016_1391_Fig6_HTML" id="MO6"></graphic>
</fig>
</p>
<p>The Sherman Island extensometer was located next to a drainage pumping station which maintained groundwater levels relatively constant with time since April 2011 (Fig. 
<xref rid="Fig7" ref-type="fig">7</xref>
). The inelastic subsidence rate during April 2011 to April 2015 was 0.52 cm yr
<sup>–1</sup>
. Soil organic matter content was 16.9 % in 2015 (Table
<xref rid="Tab1" ref-type="table">1</xref>
).
<fig id="Fig7">
<label>Fig. 7</label>
<caption>
<p>Time-series plot of land-surface elevation changes and groundwater level, Sherman Island</p>
</caption>
<graphic xlink:href="10040_2016_1391_Fig7_HTML" id="MO7"></graphic>
</fig>
</p>
</sec>
<sec id="Sec19">
<title>Twitchell rice leveling and extensometer data</title>
<p>The results of surveying in Twitchell Island rice fields indicated average elevation changes ranging from −1.7 to 2.1 cm yr
<sup>–1</sup>
from 2009–2013. The average elevation-change rate for all seven locations in the rice fields was 0.05 cm yr
<sup>–1</sup>
. Spatial variations in rates are primarily due to land disturbance resulting from agricultural activities which include disking for weed control and preparation for planting, harvest and subsequent incorporation of plant residues. During 2012–2015, extensometer data from the Twitchell Island experimental rice field exhibits elastic and inelastic land-surface elevations and changes associated with seasonal cycles in groundwater level changes (Fig. 
<xref rid="Fig8" ref-type="fig">8</xref>
). Flooded conditions were maintained in the rice fields during the growing season (late spring through late summer) and during the winter. Fields were drained before seeding in the spring and harvest in late summer/early fall. Six periods of flooding and draining are evident from 2012–2015 (Fig. 
<xref rid="Fig8" ref-type="fig">8</xref>
).
<fig id="Fig8">
<label>Fig. 8</label>
<caption>
<p>Time-series plot of land-surface elevation and depth to groundwater at Twitchell Island rice extensometer</p>
</caption>
<graphic xlink:href="10040_2016_1391_Fig8_HTML" id="MO8"></graphic>
</fig>
</p>
<p>The 2012 data illustrate key events in the rice cultivation cycle. Land surface elevation initially followed the decline in the groundwater levels as the field was dewatered prior to planting. After planting and through the growing season, groundwater levels were about 20 cm below land surface at the observation well. During that period, land-surface elevations steadily increased. At the end of the growing season, the field was dewatered for harvest in September and groundwater levels and land-surface elevations decreased. The SET was moved off the field for harvesting and cultivation. When the SET was returned later to the exact same location in the field on 6 December 2012, land surface had risen by 17 mm, relative to the previous measurement on 11 November 2012 due to increased groundwater levels resultant from flooding and land disturbance. The water table was maintained through the winter at roughly the same level as during the growing season. The field was dewatered for preparation for planting beginning in late February at which point land-surface elevations declined to about the same level measured in November 2012 as shown by the manual measurement in March 2012. Similar oscillations were observed during 2013–2015.</p>
<p>Net changes in land-surface elevations were estimated from the data shown in Fig. 
<xref rid="Fig8" ref-type="fig">8</xref>
by comparing annual land-surface elevation measurements at times when groundwater levels were equal. Land-surface elevations measured when the groundwater levels were shallowest (0.2 m) and deepest (1.15 m) indicated an overall average net accretion of about 0.8 cm yr
<sup>–1</sup>
(Fig. 
<xref rid="Fig9" ref-type="fig">9</xref>
). Specifically, the average accretion rates for the shallow (0.2 m) and deep groundwater depths (1.15 m) were 0.47 and 1.2 cm yr
<sup>–1</sup>
, respectively. The average accretion rate for all measurements was 0.84 cm yr
<sup>–1</sup>
.
<fig id="Fig9">
<label>Fig. 9</label>
<caption>
<p>Time-series plot of land-surface elevation for equal depth to groundwater</p>
</caption>
<graphic xlink:href="10040_2016_1391_Fig9_HTML" id="MO9"></graphic>
</fig>
</p>
</sec>
</sec>
<sec id="Sec20">
<title>Estimation of Delta-wide subsidence</title>
<sec id="Sec21">
<title>Distribution of soil organic matter</title>
<p>Figure
<xref rid="Fig3" ref-type="fig">3</xref>
indicates that relative to values reported in the soil surveys, average soil organic matter has decreased over time and that average values are likely about 72 % of the reported values where there is over 20 % soil organic matter content; however, Fig. 
<xref rid="Fig4" ref-type="fig">4</xref>
indicates that overall, the medians of recently measured values were slightly lower than the mid-range of soil-survey reported values as indicated by the slope of 0.96. The average soil organic-matter percentage values in Fig. 
<xref rid="Fig10" ref-type="fig">10</xref>
were based on the data from soil surveys.
<fig id="Fig10">
<label>Fig. 10</label>
<caption>
<p>Map showing distribution of soil organic matter percentage, Sacramento-San Joaquin Delta</p>
</caption>
<graphic xlink:href="10040_2016_1391_Fig10_HTML" id="MO10"></graphic>
</fig>
</p>
<p>Mapped soil organic matter content varied from less than 6 % to over 52 % throughout the Delta (Fig. 
<xref rid="Fig10" ref-type="fig">10</xref>
). The distribution of soil organic matter reflects geomorphologic and subsidence history. Highly organic mineral surface soils generally predominate in the western and northern Delta and organic soils are prevalent in the central, eastern and south-central Delta. The lowest organic matter content soils generally prevail in areas drained prior to 1880 near the Sacramento River. In contrast, central and eastern Delta islands, where higher organic-matter soils dominate, were reclaimed during the late 19th century or early 20th century. Prior to reclamation, islands near the Sacramento River (e.g. Sherman Island) were subject to greater fluvial deposition relative to the more quiescent environment in the central and eastern Delta.</p>
</sec>
<sec id="Sec22">
<title>Distribution of soil temperature</title>
<p>There is little spatial variation in near-surface annual air temperature. Average annual temperatures ranged from 16.25 to 17.25 °C for the entire Delta.</p>
</sec>
<sec id="Sec23">
<title>Distribution of depth to groundwater</title>
<p>There are generally small variations in depth to groundwater in the Delta due to the influence of drainage systems. Based on information in the soil surveys, average depth to groundwater ranged from 0.8 to 1.2 m for most of the organic-soil area. Exceptions included the western Delta (notably Sherman and Jersey islands where pasture is the predominant land use). In these areas, depth to groundwater levels varied from 0 to 0.8 m.</p>
</sec>
<sec id="Sec24">
<title>Estimated subsidence rates</title>
<p>The SUBCALC model was re-calibrated for data for land-surface elevation change and greenhouse gas emissions (Table
<xref rid="Tab1" ref-type="table">1</xref>
) collected throughout the Delta (Fig. 
<xref rid="Fig1" ref-type="fig">1</xref>
). Average depth to groundwater varied from 0.50 to 1.39 m and soil organic matter content varied from 14.9 to 39.7 %. There was a small average annual temperature variation among the sites (14.9–16.4 °C) and therefore small variability in
<italic>V</italic>
<sub>max</sub>
. The values of
<italic>K</italic>
<sub>m</sub>
varied inversely with soil organic matter content. The root-mean square error (RMSE) was calculated as a goodness of fit parameter. For subsidence where depth to groundwater, soil temperature and soil-organic matter content are known, RMSE values indicated that model predictions were accurate within ±0.10 cm yr
<sup>–1</sup>
, whereas carbon fluxes were accurately predicted within ±0.006 g C cm
<sup>–2</sup>
yr
<sup>–1</sup>
.</p>
<p>Estimated present-day subsidence rates, which varied from 0.28 to 1.8 cm yr
<sup>–1</sup>
based on inputs for depth to groundwater, average annual soil temperature and soil organic matter content are shown in Fig. 
<xref rid="Fig11" ref-type="fig">11</xref>
. The highest rates (over 0.9 cm yr
<sup>–1</sup>
) correspond to high organic-matter soils in the central Delta where estimated soil organic matter content was over 40 % (Figs. 
<xref rid="Fig10" ref-type="fig">10</xref>
and
<xref rid="Fig11" ref-type="fig">11</xref>
). Rates generally equal to or lower than 0.9 cm yr
<sup>–1</sup>
corresponded to the western Delta where soil organic matter content generally ranged from less than 6 to over 17 %. Estimated subsidence was also low or nil in the northern, eastern and southern Delta where organic matter contents were generally less than 15 %. At locations where rice cultivation and wetlands have been implemented (Sherman, Twitchell, Brack, Canal Ranch, Wright Elmwood), zero subsidence was assumed; active subsidence occurs where there is peat at or below elevation –2 m and the highest rates occur below –4 m (Fig. 
<xref rid="Fig11" ref-type="fig">11</xref>
).
<fig id="Fig11">
<label>Fig. 11</label>
<caption>
<p>Map showing distribution of present-day modeled subsidence rates</p>
</caption>
<graphic xlink:href="10040_2016_1391_Fig11_HTML" id="MO11"></graphic>
</fig>
</p>
</sec>
</sec>
</sec>
<sec id="Sec25" sec-type="discussion">
<title>Discussion</title>
<sec id="Sec26">
<title>Subsidence rates</title>
<p>Measured present-day subsidence rates vary substantially in the Delta and are primarily related to soil organic matter content and secondarily to water- and land-management practices, which determine depth to groundwater. Deverel and Leighton (
<xref ref-type="bibr" rid="CR16">2010</xref>
) measured land-surface elevations at 51 locations in 2006 where Weir and colleagues determined elevations in 1978 on Bacon Island (Fig. 
<xref rid="Fig1" ref-type="fig">1</xref>
). Subsidence rates ranged from 1.5 to 3.0 cm yr
<sup>–1</sup>
where soil organic matter content varied from 15 to over 60 %. On Sherman Island, Deverel and Leighton (
<xref ref-type="bibr" rid="CR16">2010</xref>
) resurveyed power-pole foundations in 2006 originally surveyed in 1988 by Rojstaczer et al. (
<xref ref-type="bibr" rid="CR51">1991</xref>
) and reported subsidence rates ranging from 0.6 to 2 cm yr
<sup>–1</sup>
where soil organic matter content varied from 3 to 20 %.</p>
<p>Figure
<xref rid="Fig11" ref-type="fig">11</xref>
shows generally lower rates for most of the Delta relative to those reported by Deverel and Leighton (
<xref ref-type="bibr" rid="CR16">2010</xref>
). Two key factors explain lower subsidence rates: lower soil organic matter content resultant from ongoing oxidation, and recalibration of SUBCALC with more recent data for subsidence and carbon emissions. For example on Bacon Island (Fig. 
<xref rid="Fig1" ref-type="fig">1</xref>
), Deverel and Leighton (
<xref ref-type="bibr" rid="CR16">2010</xref>
) reported soil organic matter contents over 60 %. Using the adjusted values described in the ‘
<xref rid="Sec5" ref-type="sec">Data sources and methods</xref>
’ section resulted in lower present-day average organic matter content and subsidence rates for the eastern portion of Bacon Island as ranging from 0.9 to 1.4 cm yr
<sup>–1</sup>
. The average subsidence rate was 1.3 cm yr
<sup>–1</sup>
. For comparison, Deverel and Leighton (
<xref ref-type="bibr" rid="CR16">2010</xref>
) reported an average rate of 2.2 cm yr
<sup>–1</sup>
for this area from 1978 to 2006. Within any soil type, soil organic matter content can vary substantially (Fig. 
<xref rid="Fig4" ref-type="fig">4</xref>
). The SUBCALC-based subsidence maps presented here therefore provide spatially averaged rates that are representative of areas of similar soil organic matter content at the level of tens to hundreds of hectares.</p>
<p>To better assess SUBCALC-simulated present-day subsidence rates for varying organic matter content, on Bacon Island, present-day rates were simulated at the locations where Deverel and Leighton (
<xref ref-type="bibr" rid="CR16">2010</xref>
) reported 1978–2006 subsidence rates ranging from 1.5 to 3.8 cm yr
<sup>–1</sup>
and the average rate was 2.2 cm yr
<sup>–1</sup>
, whereas soil organic matter content ranged from 14 to 60 % and the average soil organic matter content was 39 %. SUBCALC present-day rates ranged from 0.7 to 2.0 cm yr
<sup>–1</sup>
for the same range of organic matter content percentages, with the average rate being 1.2 cm yr
<sup>–1</sup>
.</p>
<p>A similar comparison was made between SUBCALC-simulated present-day subsidence rates and Sherman Island 1988–2006 subsidence rates, also reported in Deverel and Leighton (
<xref ref-type="bibr" rid="CR16">2010</xref>
). Among Sherman Island measurement points, subsidence rates ranged from 0.6 to 2.0 cm yr
<sup>–1</sup>
and the average rate was 1.2 cm yr
<sup>–1</sup>
, while soil organic matter content ranged from 3 to 20 % with an average of 10 %. SUBCALC present-day rates ranged from 0.3 to 0.6 cm yr
<sup>–1</sup>
for the same range of organic matter content percentages, and the average rate was 0.4 cm yr
<sup>–1</sup>
. These comparisons illustrated the likely variation within a central (Bacon) and western (Sherman) Delta island and indicates that current subsidence rates are lower than rates in previous decades.</p>
<p>Priyanka Sharma and colleagues (P. Sharma, Jet Propulsion Laboratory, California Institute of Technology, personal communication, 2015) used UAVSAR data to estimate subsidence rates from 2009 to 2014 over all of Sherman Island. Reported rates ranged from 0 to 5 cm yr
<sup>–1</sup>
and averaged 1.5 cm yr
<sup>–1</sup>
. By comparison, SUBCALC-simulated subsidence rates ranged from 0.3 to 1.1 cm yr
<sup>–1</sup>
and averaged 0.5 cm yr
<sup>–1</sup>
at Sherman Island. Within the Sharma study, the UAVSAR method was found to have over-predicted the subsidence rate (0.7 cm yr
<sup>–1</sup>
) at the location of an extensometer, which measured a rate of 0.4 cm yr
<sup>–1</sup>
over the same period. The reported uncertainty in the UAVSAR estimates is about ≤ 0.2 cm yr
<sup>–1</sup>
.</p>
</sec>
<sec id="Sec27">
<title>Rice as subsidence-mitigation land use</title>
<p>Direct (elevation-change measurements) and indirect estimates (micrometeorological and organic-matter mineralization) of subsidence and accretion are generally consistent in preliminarily evaluations indicating that rice stops or greatly reduces subsidence by providing above- and below-ground plant residue which is incorporated into the soil and by reducing the rate of peat oxidation under saturated conditions. Data derived from leveling surveys and the extensometer demonstrate substantial temporal and spatial variability in land-surface elevation changes and indicate a small net accretion rate. These data illustrate the difficulty in estimating short-term land-surface elevation changes in this and other systems where there is land disturbance and elastic shrinking/swelling changes associated with short-term and seasonal groundwater elevation and soil moisture content changes (e.g. Zanello et al.
<xref ref-type="bibr" rid="CR72">2011</xref>
). Similar difficulties were noted with data collected by Weir (
<xref ref-type="bibr" rid="CR70">1950</xref>
) in Deverel and Leighton (
<xref ref-type="bibr" rid="CR16">2010</xref>
). Longer-term, high-quality land-surface elevation measurements for rice will provide more definitive answers.</p>
<p>Reported subsidence estimates in the Twitchell Island rice fields from indirect methods indicate small rates of subsidence. Based on the eddy-covariance determination of carbon loss, Hatala et al. (
<xref ref-type="bibr" rid="CR30">2012</xref>
) estimated subsidence in the Twitchell Island rice field at 0.1 to 0.14 cm yr
<sup>–1</sup>
during 2 years. During a 1-year study, Kirk et al. (
<xref ref-type="bibr" rid="CR39">2015</xref>
) estimated soil organic matter-nitrogen mineralization rates at four locations in the Twitchell Island rice field and used these with soil carbon:nitrogen ratios to estimate subsidence rates ranging from 0.07 to 0.11 cm yr
<sup>–1</sup>
, in close agreement with results reported by Hatala et al. There is uncertainty in these indirect estimates.</p>
<p>In estimating subsidence rates using nitrogen mineralization from peat, Kirk et al. (
<xref ref-type="bibr" rid="CR39">2015</xref>
) determined the annual nitrogen budget and by accounting for fertilizer application and plant uptake, calculated the annual mineralized nitrogen as a source to plant nitrogen uptake. Nitrogen from groundwater was estimated in situ using groundwater and soil-water samples in mesocosms. Kirk et al. (
<xref ref-type="bibr" rid="CR39">2015</xref>
) assumed that groundwater nitrogen contributions resulted from mineralization during the year of investigation. Given low hydraulic conductivity and high porosity of the organic soils, groundwater nitrogen likely resulted from mineralization during previous years. Additionally, Kirk et al. (
<xref ref-type="bibr" rid="CR39">2015</xref>
) applied an f
<sub>min</sub>
factor of 0.67 based on Deverel and Leighton (
<xref ref-type="bibr" rid="CR16">2010</xref>
), a factor that is not applicable because any subsidence in rice fields results from oxidation of the organic soil. In Deverel and Leighton (
<xref ref-type="bibr" rid="CR16">2010</xref>
), the value of 0.67 was the estimated fraction of organic soil oxidation contributing to subsidence originally reported in Deverel and Rojstaczer (
<xref ref-type="bibr" rid="CR17">1996</xref>
). The remaining fraction was attributed to consolidation due to deepening of drainage ditches. Because the rice field drainage ditches have not been altered, the sole cause of subsidence would be oxidation of soil organic matter. Removing the groundwater contribution and f
<sub>min</sub>
factor from the Kirk et al. (
<xref ref-type="bibr" rid="CR39">2015</xref>
) calculations resulted in a subsidence rate of about 0.02 cm yr
<sup>1</sup>
; furthermore, Kirk et al. (
<xref ref-type="bibr" rid="CR39">2015</xref>
) used a plant nitrogen uptake efficiency of 50 % derived from the literature for mineral soils and fertilizer and stated that the values could range as high at 70 %. Using the 70 % value in Kirk et al.’s (
<xref ref-type="bibr" rid="CR39">2015</xref>
) calculation and removing the groundwater contribution and the f
<sub>min</sub>
factor denominator resulted in accretion (0.001 cm yr
<sup>1</sup>
) to a small amount of subsidence (0.01 cm yr
<sup>1</sup>
).</p>
<p>Hatala et al. (
<xref ref-type="bibr" rid="CR30">2012</xref>
) estimated the net carbon balance (carbon dioxide sequestered – methane emitted + planted seed – harvested grain) and used data for soil bulk density and soil carbon to estimate subsidence rates in rice during 2009–2011. The range of carbon balance values presented was used here to estimate subsidence values as low as 0.07 cm yr
<sup>1</sup>
. Using data presented in Knox et al. (
<xref ref-type="bibr" rid="CR41">2015</xref>
) for 2013, subsidence rates ranging from 0.02 to 0.13 cm yr
<sup>1</sup>
were estimated. The lower range of estimates of subsidence and accretion based on the data published in Hatala et al. (
<xref ref-type="bibr" rid="CR30">2012</xref>
), Knox et al. (
<xref ref-type="bibr" rid="CR41">2015</xref>
) and Kirk et al. (
<xref ref-type="bibr" rid="CR39">2015</xref>
), are more consistent with measured land-surface elevation changes which indicate a small rate of overall accretion. Qualitatively and in light of uncertainty, the preponderance of evidence summarized here preliminarily indicates that rice cultivation greatly reduces subsidence or may slightly reverse the effects of subsidence. Longer-term data will provide improved quantification of the long-term subsidence mitigation benefit due to rice cultivation. Knox and colleagues (S. Knox, University of California-Berkeley, personal communication, 2016) presented results of 6 years of eddy covariance measurements of CO
<sub>2</sub>
and CH
<sub>4</sub>
fluxes in the Twitchell Island rice field. These data show heretofore unreported substantial annual variability in photosynthesis and methane fluxes driven primarily by variability in soil temperatures and resulted in substantial variability in soil carbon budgets.</p>
</sec>
<sec id="Sec28">
<title>Mitigation</title>
<p>The primary Delta subsidence mitigation tools are rice cultivation and permanently flooded wetlands. As demonstrated here and elsewhere (e.g. Miller et al.
<xref ref-type="bibr" rid="CR45">2008</xref>
; Deverel et al.
<xref ref-type="bibr" rid="CR21">2014</xref>
), both of these land use practices stop, greatly reduce or reverse the effects of subsidence. The work described here and in Deverel et al. (
<xref ref-type="bibr" rid="CR21">2014</xref>
) provide guidance for implementation of these land-use changes. Areas below elevations of –2 m are candidate areas for implementation because there is active subsidence occurring (Fig. 
<xref rid="Fig11" ref-type="fig">11</xref>
). Moreover and consistently, Deverel et al. (
<xref ref-type="bibr" rid="CR22">2015</xref>
) demonstrated that artesian conditions prevail below –2 m and 81 % of wet, non- or marginally farmable areas were at or below –2 m. Implementation of rice and wetlands in these areas will prevent or reduce subsidence and associated consequences.</p>
</sec>
</sec>
<sec id="Sec29" sec-type="conclusion">
<title>Summary and conclusions</title>
<p>Subsidence due primarily to oxidation of soil organic matter in the Sacramento-San Joaquin Delta affects sustainability of California’s water supply system and local agriculture. By reducing the landmass and resistance to hydraulic pressure from adjacent channels, subsidence has contributed to levee failure and island inundation which potentially affects water for use by over 25 million Californians and irrigation of millions of hectares of agricultural land. Since the mid-nineteenth century, oxidation has resulted in up to 8 m of subsidence. Subsidence rates have declined with time due to the disappearance of about 2.5 billion m
<sup>3</sup>
of organic soil and consequent decreases in soil organic matter content and changing management practices. Present-day Delta subsidence rates have not heretofore been extensively recorded or estimated.</p>
<p>Land-surface elevation data were collected to assess present-day subsidence rates and preliminarily evaluate rice as a land use for subsidence mitigation. To depict Delta-wide present-day rates of subsidence, the previously developed and reported SUBCALC model was revisited and calibrated using recent subsidence rates and carbon flux data. The primary inputs to the SUBCALC model include depth to groundwater, soil organic matter content and soil temperatures which were spatially estimated using multiple data sources. These inputs were used to map estimated subsidence rates. Land-surface elevation change data was collected and evaluated relative to indirect estimates of subsidence and accretion using carbon and nitrogen flux data for rice.</p>
<p>Extensometer data in a cornfield on Twitchell Island demonstrate seasonal variations in land-surface elevations associated with groundwater-level fluctuations from 2009 through 2013 and an inelastic subsidence rate of 0.83 cm yr
<sup>–1</sup>
. Leveling data resulted in a similar estimated subsidence rate of 0.7 cm yr
<sup>–1</sup>
from 2000–2012 on Twitchell Island. The Sherman Island extensometer data indicated a rate of 0.52 cm yr
<sup>–1</sup>
where there was lower soil organic matter content. Calibration of the SUBCALC model indicated accuracy of ±0.10 cm yr
<sup>–1</sup>
where depth to groundwater, soil organic matter content and temperature are known, while regional estimates of subsidence based on spatial variations in estimated soil organic matter content, depth to groundwater and soil temperature range from less than 0.3 to over 1.8 cm yr
<sup>–1</sup>
. The primary uncertainty is the distribution of soil organic matter content which results in spatial averaging in the mapping of subsidence rates at the level of tens to hundreds of hectares.</p>
<p>Analysis of leveling and extensometer data in the Twitchell Island rice field resulted in an estimated accretion rate of 0.02 to 0.8 cm yr
<sup>–1</sup>
. Indirect estimates based on measurements of carbon fluxes and nitrogen mineralization resulted in estimates of low subsidence rates to low accretion rates. The preponderance of evidence presented here preliminarily demonstrates that rice will stop or greatly reduce subsidence for most of the Delta. Areas below elevations of –2 m are candidate areas for implementation because there is active subsidence occurring at rates greater than 0.4 cm yr
<sup>–1</sup>
.</p>
</sec>
</body>
<back>
<ack>
<title>Acknowledgements</title>
<p>We gratefully acknowledge funding from the California Department of Water Resources, Department of Land, Air, and Water Resources at the University of California, Davis, the US. Department of Agriculture Grant NIFA No. 2011-67003-30371, Subaward No. 201015552-06 and the Jet Propulsion Laboratory, California Institute of Technology and National Aeronautics and Space Administration. We thank California Department of Water Resources personnel under the direction of Bryan Brock for surveying data on Twitchell Island. We also thank two anonymous reviewers for their helpful comments.</p>
</ack>
<ref-list id="Bib1">
<title>References</title>
<ref id="CR1">
<mixed-citation publication-type="other">Atwater BF (1982) Geologic maps of the Sacramento-San Joaquin Delta. US Geoll Surv Miscell Field Studies Map MF-1401</mixed-citation>
</ref>
<ref id="CR2">
<mixed-citation publication-type="other">Bawden GW, Howle J, Bond S, Shriro M, Buck P (2014) Three-dimensional imaging, change detection, and stability assessment during the centerline trench levee seepage experiment using terrestrial light detection and ranging technology, Twitchell Island, California, 2012. US Geol Surv Open-File Rep 2014-1092</mixed-citation>
</ref>
<ref id="CR3">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Boumans</surname>
<given-names>RMJ</given-names>
</name>
<name>
<surname>Day</surname>
<given-names>JWJ</given-names>
</name>
</person-group>
<article-title>High precision measurements of sediment elevation in shallow coastal areas using a sedimentation-erosion table</article-title>
<source>Estuaries</source>
<year>1993</year>
<volume>16</volume>
<issue>2</issue>
<fpage>375</fpage>
<lpage>380</lpage>
<pub-id pub-id-type="doi">10.2307/1352509</pub-id>
</element-citation>
</ref>
<ref id="CR4">
<mixed-citation publication-type="other">Brooks BA, Bawden G, Manjunath D, Werner C, Knowles N, Foster J (2012) Contemporaneous subsidence and levee overtopping potential, Sacramento-San Joaquin Delta, California, San Franc Estuary Watershed Sci 10(1)</mixed-citation>
</ref>
<ref id="CR5">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Browder</surname>
<given-names>JA</given-names>
</name>
<name>
<surname>Volk</surname>
<given-names>BG</given-names>
</name>
</person-group>
<article-title>Systems model of carbon transformations in soil subsidence</article-title>
<source>Ecol Modeling</source>
<year>1978</year>
<volume>5</volume>
<issue>4</issue>
<fpage>269</fpage>
<lpage>292</lpage>
<pub-id pub-id-type="doi">10.1016/0304-3800(78)90038-8</pub-id>
</element-citation>
</ref>
<ref id="CR6">
<mixed-citation publication-type="other">California Central Valley Flood Control Association (2011) Comments on Flood Risk White Paper, Delta Stewardship Council, January 2011.
<ext-link ext-link-type="uri" xlink:href="http://www.deltacouncil.ca.gov/sites/default/files/documents/files/CVFCA_012011_0.pdf">http://www.deltacouncil.ca.gov/sites/default/files/documents/files/CVFCA_012011_0.pdf</ext-link>
. Accessed 07 July 2015</mixed-citation>
</ref>
<ref id="CR7">
<mixed-citation publication-type="other">California Department of Water Resources (CDWR) (2007) 2007 Sacramento-San Joaquin Delta LiDAR acquisition, California Department of Water Resources. CDWR, Sacramento, CA.
<ext-link ext-link-type="uri" xlink:href="http://www.atlas.ca.gov/downloadhtml%23/casil/imageryBaseMapsLandCover/lidar2009">http://www.atlas.ca.gov/downloadhtml#/casil/imageryBaseMapsLandCover/lidar2009</ext-link>
. Accessed 13 July 2015</mixed-citation>
</ref>
<ref id="CR8">
<mixed-citation publication-type="other">California Department of Water Resources (CDWR) (2009) Delta Risk Management Strategy (DRMS) Phase 1, California Dept. of Water Resources. CDWR, Sacramento, CA.
<ext-link ext-link-type="uri" xlink:href="http://www.water.ca.gov/floodsafe/fessro/levees/drms/phase1_information.cfm">http://www.water.ca.gov/floodsafe/fessro/levees/drms/phase1_information.cfm</ext-link>
. Accessed 13 July 2015</mixed-citation>
</ref>
<ref id="CR9">
<mixed-citation publication-type="other">California Department of Water Resources (CDWR) (2015) California State Water Project Overview, 2015. CDWR, Sacramento, CA.
<ext-link ext-link-type="uri" xlink:href="http://www.water.ca.gov/swp/index.cfm">http://www.water.ca.gov/swp/index.cfm</ext-link>
. Accessed 2 December 2015</mixed-citation>
</ref>
<ref id="CR10">
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Carlton</surname>
<given-names>AB</given-names>
</name>
<name>
<surname>Schultz</surname>
<given-names>HB</given-names>
</name>
</person-group>
<source>Annual statements of progress for project 1686, 1955 to 1966, peat land conservation and peat dust abatement</source>
<year>1966</year>
<publisher-loc>Davis, CA</publisher-loc>
<publisher-name>Agricultural Experiment Station, Dept. of Soils and Plant Nutrition</publisher-name>
</element-citation>
</ref>
<ref id="CR11">
<mixed-citation publication-type="other">Cohen DA, Deverel SJ, Johnson LA (1998) Feasibility study of differential SAR interferometry for subsidence monitoring in the Sacramento-San Joaquin Delta. Proc. IEEE Int. IGARSS 1998, vol 3, Seattle, WA, August 1998, pp 1629–1631. doi:10.1109/IGARSS.1998.691658</mixed-citation>
</ref>
<ref id="CR12">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cook</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Coleman</surname>
<given-names>D</given-names>
</name>
</person-group>
<article-title>Disaster strikes Isleton</article-title>
<source>Reclamation Era</source>
<year>1973</year>
<volume>59</volume>
<fpage>1</fpage>
<lpage>9</lpage>
</element-citation>
</ref>
<ref id="CR13">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Couwenberg</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Hooijer</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Towards robust subsidence-based soil carbon emission factors for peat soils in South-East Asia, with special reference to oil palm plantations</article-title>
<source>Mires Peat</source>
<year>2013</year>
<volume>12</volume>
<issue>1</issue>
<fpage>1</fpage>
<lpage>13</lpage>
</element-citation>
</ref>
<ref id="CR14">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Davidson</surname>
<given-names>EA</given-names>
</name>
<name>
<surname>Sudeep</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Caramori</surname>
<given-names>SS</given-names>
</name>
<name>
<surname>Savage</surname>
<given-names>K</given-names>
</name>
</person-group>
<article-title>The Dual Arrhenius and Michaelis–Menten kinetics model for decomposition of soil organic matter at hourly to seasonal time scales</article-title>
<source>Global Change Biol</source>
<year>2012</year>
<volume>18</volume>
<fpage>371</fpage>
<lpage>384</lpage>
<pub-id pub-id-type="doi">10.1111/j.1365-2486.2011.02546.x</pub-id>
</element-citation>
</ref>
<ref id="CR15">
<mixed-citation publication-type="other">Delta Protection Commission (2012) Economic sustainability plan for the Sacramento-San Joaquin Delta.
<ext-link ext-link-type="uri" xlink:href="http://www.delta.ca.gov/res/docs/ESP/ESP_P2_FINAL.pdf">http://www.delta.ca.gov/res/docs/ESP/ESP_P2_FINAL.pdf</ext-link>
. Accessed 14 April 2015</mixed-citation>
</ref>
<ref id="CR16">
<mixed-citation publication-type="other">Deverel SJ, Leighton DA (2010) Historic, recent, and future subsidence, Sacramento-San Joaquin Delta, California, USA. San Franc Estuary Watershed Sci 8(2). doi:10.15447/sfews.2010v8iss2art1</mixed-citation>
</ref>
<ref id="CR17">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Deverel</surname>
<given-names>SJ</given-names>
</name>
<name>
<surname>Rojstaczer</surname>
<given-names>SA</given-names>
</name>
</person-group>
<article-title>Subsidence of agricultural lands in the Sacramento-San Joaquin Delta, California: role of aqueous and gaseous carbon fluxes</article-title>
<source>Water Resour Res</source>
<year>1996</year>
<volume>32</volume>
<issue>8</issue>
<fpage>2359</fpage>
<lpage>2367</lpage>
<pub-id pub-id-type="doi">10.1029/96WR01338</pub-id>
</element-citation>
</ref>
<ref id="CR18">
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Deverel</surname>
<given-names>SJ</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>B</given-names>
</name>
<name>
<surname>Rojstaczer</surname>
<given-names>SA</given-names>
</name>
</person-group>
<person-group person-group-type="editor">
<name>
<surname>Borchers</surname>
<given-names>JW</given-names>
</name>
</person-group>
<article-title>Subsidence in the Sacramento-San Joaquin Delta</article-title>
<source>Proceedings of the Joseph Poland Subsidence Symposium, Association of Engineering Geologists. Spec. Publ. no. 8</source>
<year>1998</year>
<publisher-loc>Belmont, CA</publisher-loc>
<publisher-name>Star</publisher-name>
<fpage>489</fpage>
<lpage>502</lpage>
</element-citation>
</ref>
<ref id="CR19">
<mixed-citation publication-type="other">Deverel SJ, Leighton DA, Sola–Llonch N (2007a) Appendix C: Evaluation of island drain flow, seepage, and organic carbon loads, Sacramento-San Joaquin Delta—results from the Delta Learning Laboratory Project, objectives 2 and 3. California Department of Water Resources and CALFED Bay, Sacramento, CA</mixed-citation>
</ref>
<ref id="CR20">
<mixed-citation publication-type="other">Deverel SJ, Leighton DA, Finlay MR (2007b) Processes affecting agricultural drainwater quality and organic carbon loads in California’s Sacramento-San Joaquin Delta. San Franc Estuary Watershed Sci 5(2). doi:10.15447/sfews.2007v5iss2art2</mixed-citation>
</ref>
<ref id="CR21">
<mixed-citation publication-type="other">Deverel SJ, Ingrum T, Lucero C, Drexler JZ (2014) Impounded marshes on subsided islands: simulated vertical accretion, processes and effects, Sacramento-San Joaquin Delta, CA USA. San Franc Estuary Watershed Sci 12(2). doi:10.15447/sfews.2014v12iss2art5</mixed-citation>
</ref>
<ref id="CR22">
<mixed-citation publication-type="other">Deverel SJ, Lucero C, Bachand S (2015) Evolution of arability and land use, Sacramento-San Joaquin Delta, California. San Franc Estuary Watershed Sci 13(2). doi:10.15447/sfews.2015v13iss2art4</mixed-citation>
</ref>
<ref id="CR23">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Drexler</surname>
<given-names>JZ</given-names>
</name>
<name>
<surname>de Fontaine</surname>
<given-names>CS</given-names>
</name>
<name>
<surname>Brown</surname>
<given-names>TA</given-names>
</name>
</person-group>
<article-title>Peat accretion histories during the past 6000 years in marshes of the Sacramento-San Joaquin Delta, California, USA</article-title>
<source>Estuar Coasts</source>
<year>2009</year>
<volume>32</volume>
<fpage>871</fpage>
<lpage>892</lpage>
<pub-id pub-id-type="doi">10.1007/s12237-009-9202-8</pub-id>
</element-citation>
</ref>
<ref id="CR24">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Drexler</surname>
<given-names>JZ</given-names>
</name>
<name>
<surname>de Fontaine</surname>
<given-names>CS</given-names>
</name>
<name>
<surname>Deverel</surname>
<given-names>SJ</given-names>
</name>
</person-group>
<article-title>The legacy of wetland drainage on the remaining peat in the Sacramento-San Joaquin Delta, California, USA</article-title>
<source>Wetlands</source>
<year>2009</year>
<volume>29</volume>
<issue>1</issue>
<fpage>372</fpage>
<lpage>386</lpage>
<pub-id pub-id-type="doi">10.1672/08-97.1</pub-id>
</element-citation>
</ref>
<ref id="CR25">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ewing</surname>
<given-names>JM</given-names>
</name>
<name>
<surname>Vepraskas</surname>
<given-names>MJ</given-names>
</name>
</person-group>
<article-title>Estimating primary and secondary subsidence in an organic soil 15, 20, and 30 years after drainage</article-title>
<source>Wetlands</source>
<year>2006</year>
<volume>26</volume>
<issue>1</issue>
<fpage>119</fpage>
<lpage>130</lpage>
<pub-id pub-id-type="doi">10.1672/0277-5212(2006)26[119:EPASSI]2.0.CO;2</pub-id>
</element-citation>
</ref>
<ref id="CR26">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fleck</surname>
<given-names>JA</given-names>
</name>
<name>
<surname>Bossio</surname>
<given-names>DA</given-names>
</name>
<name>
<surname>Fujii</surname>
<given-names>R</given-names>
</name>
</person-group>
<article-title>Dissolved organic carbon and disinfection by-product precursor release from managed peat soils</article-title>
<source>Environ Qual</source>
<year>2004</year>
<volume>33</volume>
<issue>2</issue>
<fpage>465</fpage>
<lpage>475</lpage>
<pub-id pub-id-type="doi">10.2134/jeq2004.4650</pub-id>
</element-citation>
</ref>
<ref id="CR27">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Florsheim</surname>
<given-names>JL</given-names>
</name>
<name>
<surname>Dettinger</surname>
<given-names>MD</given-names>
</name>
</person-group>
<article-title>Climate and floods still govern California levee breaks</article-title>
<source>Geophys Res Lett</source>
<year>2007</year>
<volume>34</volume>
<pub-id pub-id-type="doi">10.1029/2007GL031702</pub-id>
</element-citation>
</ref>
<ref id="CR28">
<mixed-citation publication-type="other">Gaddie V, Mierza M, Marr J (2006) Levee failures in the Sacramento-San Joaquin River Delta. California Department of Water Resources, Sacramento, CA.
<ext-link ext-link-type="uri" xlink:href="http://www.water.ca.gov/floodmgmt/docs/DeltaLeveeFailures_FMA_200709.pdf">http://www.water.ca.gov/floodmgmt/docs/DeltaLeveeFailures_FMA_200709.pdf</ext-link>
. Accessed 13 July 2015</mixed-citation>
</ref>
<ref id="CR29">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gambolati</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Putti</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Teatini</surname>
<given-names>P</given-names>
</name>
<name>
<surname>Gasparetto Stori</surname>
<given-names>G</given-names>
</name>
</person-group>
<article-title>Subsidence due to peat oxidation and impact on drainage infrastructures in a farmland catchment south of the Venice Lagoon</article-title>
<source>Environ Geol</source>
<year>2006</year>
<volume>49</volume>
<issue>6</issue>
<fpage>814</fpage>
<lpage>820</lpage>
<pub-id pub-id-type="doi">10.1007/s00254-006-0176-6</pub-id>
</element-citation>
</ref>
<ref id="CR30">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hatala</surname>
<given-names>JA</given-names>
</name>
<name>
<surname>Detto</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Sonnentag</surname>
<given-names>O</given-names>
</name>
<name>
<surname>Deverel</surname>
<given-names>SJ</given-names>
</name>
<name>
<surname>Verfaillie</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Baldocchi</surname>
<given-names>DD</given-names>
</name>
</person-group>
<article-title>Greenhouse gas (CO2, CH4, H2O) fluxes from drained and flooded agricultural peatlands in the Sacramento-San Joaquin Delta</article-title>
<source>Agric Ecosyst Environ</source>
<year>2012</year>
<volume>150</volume>
<fpage>1</fpage>
<lpage>18</lpage>
<pub-id pub-id-type="doi">10.1016/j.agee.2012.01.009</pub-id>
</element-citation>
</ref>
<ref id="CR31">
<mixed-citation publication-type="other">Heim WA, Deverel S, Ingrum T, Piekarski W, Stephenson M (2009) Assessment of methylmercury contributions from Sacramento-San Joaquin Delta farmed islands. Final report, Central Valley Regional Water Quality Control Board, Fresno, CA</mixed-citation>
</ref>
<ref id="CR32">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hirano</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Segah</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Kusin</surname>
<given-names>K</given-names>
</name>
<name>
<surname>Limin</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Takahashi</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Osaki</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>Effects of disturbances on the carbon balance of tropical peat swamp forests</article-title>
<source>Global Chang Biol</source>
<year>2012</year>
<volume>18</volume>
<issue>11</issue>
<fpage>3410</fpage>
<lpage>3422</lpage>
<pub-id pub-id-type="doi">10.1111/j.1365-2486.2012.02793.x</pub-id>
</element-citation>
</ref>
<ref id="CR33">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hobbs</surname>
<given-names>NB</given-names>
</name>
</person-group>
<article-title>Mire morphology and the properties and behaviour of some British and foreign peats</article-title>
<source>Q J Eng Geol Hydrogeol (London)</source>
<year>1986</year>
<volume>19</volume>
<fpage>7</fpage>
<lpage>80</lpage>
<pub-id pub-id-type="doi">10.1144/GSL.QJEG.1986.019.01.02</pub-id>
</element-citation>
</ref>
<ref id="CR34">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hooijer</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Page</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Jauhiainen</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>WA</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>XX</given-names>
</name>
<name>
<surname>Idris</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Anshari</surname>
<given-names>G</given-names>
</name>
</person-group>
<article-title>Subsidence and carbon loss in drained tropical peatlands</article-title>
<source>Biogeosciences</source>
<year>2012</year>
<volume>9</volume>
<fpage>1053</fpage>
<lpage>1071</lpage>
<pub-id pub-id-type="doi">10.5194/bg-9-1053-2012</pub-id>
</element-citation>
</ref>
<ref id="CR35">
<mixed-citation publication-type="other">Jones C, Bawden G, Deverel S, Dudas J, Hensley S (2011) Characterizing land surface change and levee stability in the Sacramento-San Joaquin Delta using UAVSAR radar imagery. IEEE Int. Symposium, IGARSS 2011, Vancouver, July 2011, pp 1638–1641. doi:10.1109/IGARSS.2011.6049546</mixed-citation>
</ref>
<ref id="CR36">
<mixed-citation publication-type="other">Jones C, Bawden G, Deverel S, Dudas J, Hensley S, Yun SH (2012) Study of movement and seepage along levees using DINSAR and the airborne UAVSAR instrument, Proc. SPIE 8536, SAR Image Analysis, Modeling, and Techniques, vol XII, 85360E, Edinburgh, September 2012. doi:10.1117/12.976885</mixed-citation>
</ref>
<ref id="CR37">
<mixed-citation publication-type="other">Kaat A, Joosten H (2009) Fact book for UNFCCC policies on peat carbon emissions. Wetlands Int., Ede, The Netherlands, 24 pp</mixed-citation>
</ref>
<ref id="CR38">
<mixed-citation publication-type="other">Kerr BD, Wang B, Drexler JZ (2003) Carbon fluxes, water levels, and related environmental data, Twitchell Island, Sacramento-San Joaquin Delta, California, November 1992 through September 1995. US Geol Surv Open-File Rep 03-370</mixed-citation>
</ref>
<ref id="CR39">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kirk</surname>
<given-names>ER</given-names>
</name>
<name>
<surname>van Kessel</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Horwath</surname>
<given-names>WR</given-names>
</name>
<name>
<surname>Linquist</surname>
<given-names>BA</given-names>
</name>
</person-group>
<article-title>Estimating annual soil carbon loss in agricultural peatland soils using a nitrogen budget approach</article-title>
<source>PLoS ONE</source>
<year>2015</year>
<volume>10</volume>
<issue>3</issue>
<fpage>L e0121432</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pone.0121432</pub-id>
</element-citation>
</ref>
<ref id="CR40">
<mixed-citation publication-type="other">Knipling EB, Schroeder VN, Duncan WO (1970) CO
<sub>2</sub>
evolution from Florida organic soils. Soil Crop Sci Soc Fla Proc 30:320–326</mixed-citation>
</ref>
<ref id="CR41">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Knox</surname>
<given-names>SH</given-names>
</name>
<name>
<surname>Sturtevant</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Hatala Matthes</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Koteen</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Verfaillie</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Baldocchi</surname>
<given-names>D</given-names>
</name>
</person-group>
<article-title>Agricultural peatland restoration: effects of land-use change on greenhouse gas (CO
<sub>2</sub>
and CH
<sub>4</sub>
) fluxes in the Sacramento-San Joaquin Delta</article-title>
<source>Glob Chang Biol</source>
<year>2015</year>
<volume>21</volume>
<issue>2</issue>
<fpage>750</fpage>
<lpage>765</lpage>
<pub-id pub-id-type="doi">10.1111/gcb.12745</pub-id>
<pub-id pub-id-type="pmid">25229180</pub-id>
</element-citation>
</ref>
<ref id="CR42">
<mixed-citation publication-type="other">Maurer EP, Wood AW, Adam JC, Lettenmaier DP, Nijssen B (2002) A long-term hydrologically-based data set of land surface fluxes and states for the conterminous United States. J Climate 15:3237–3251. doi:10.1175/1520-0442(2002)015%3C3237:ALTHBD%3E2.0.CO;2. Available at
<ext-link ext-link-type="uri" xlink:href="http://www.engr.scu.edu/%7Eemaurer/gridded_obs/index_gridded_obs.html">http://www.engr.scu.edu/~emaurer/gridded_obs/index_gridded_obs.html</ext-link>
. Accessed 8 July 2015</mixed-citation>
</ref>
<ref id="CR43">
<mixed-citation publication-type="other">McElhiney MA (1992) Soil survey of San Joaquin County, California. USDA, Soil Conservation Service, Agricultural Experiment Station, Davis, CA and California Dept. of Conservation, Sacramento, CA</mixed-citation>
</ref>
<ref id="CR44">
<mixed-citation publication-type="other">Miller RL, Hastings L, Fujii R (2000) Hydrologic treatments affect gaseous carbon loss from organic soils, Twitchell Island, California, October 1995–December 1997. US Geol Surv Water Resour Invest Rep 00–4042</mixed-citation>
</ref>
<ref id="CR45">
<mixed-citation publication-type="other">Miller RL, Fram MS, Wheeler G, Fujii R (2008) Subsidence reversal in a re-established wetland in the Sacramento-San Joaquin Delta, California, USA. San Franc Estuary Watershed Sci 6(3). doi:10.15447/sfews.2008v6iss3art1</mixed-citation>
</ref>
<ref id="CR46">
<mixed-citation publication-type="other">Mount J, Twiss R (2005) Subsidence, sea level rise, seismicity in the Sacramento-San Joaquin Delta. San Franc Estuary Watershed Sci 3(1). doi:10.15447/sfews.2005v3iss1art7</mixed-citation>
</ref>
<ref id="CR47">
<mixed-citation publication-type="other">Pellerin B, Anderson F, Bergamaschi B (2013) Assessing the role of winter flooding on baseline greenhouse gas fluxes from cornfields in the Sacramento-San Joaquin Bay Delta. CEC-500-2014-077, California Energy Commission, Sacramento, CA</mixed-citation>
</ref>
<ref id="CR48">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Prokopovitch</surname>
<given-names>NP</given-names>
</name>
</person-group>
<article-title>Subsidence of peat in California and Florida</article-title>
<source>Bull Assoc Eng Geol</source>
<year>1985</year>
<volume>22</volume>
<issue>4</issue>
<fpage>395</fpage>
<lpage>420</lpage>
</element-citation>
</ref>
<ref id="CR49">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rojstaczer</surname>
<given-names>SA</given-names>
</name>
<name>
<surname>Deverel</surname>
<given-names>SJ</given-names>
</name>
</person-group>
<article-title>Time dependence in atmospheric carbon inputs from drainage of organic soils</article-title>
<source>Geophys Res Lett</source>
<year>1993</year>
<volume>20</volume>
<fpage>1383</fpage>
<lpage>1386</lpage>
<pub-id pub-id-type="doi">10.1029/93GL01339</pub-id>
</element-citation>
</ref>
<ref id="CR50">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rojstaczer</surname>
<given-names>SA</given-names>
</name>
<name>
<surname>Deverel</surname>
<given-names>SJ</given-names>
</name>
</person-group>
<article-title>Land subsidence in drained histosols and highly organic mineral soils of California</article-title>
<source>Soil Sci Soc Am J</source>
<year>1995</year>
<volume>59</volume>
<fpage>1162</fpage>
<lpage>1167</lpage>
<pub-id pub-id-type="doi">10.2136/sssaj1995.03615995005900040031x</pub-id>
</element-citation>
</ref>
<ref id="CR51">
<mixed-citation publication-type="other">Rojstaczer SA, Hamon RE, Deverel SJ, Massey CA (1991) Evaluation of selected data to assess the causes of subsidence in the Sacramento-San Joaquin Delta, California. US Geol Surv Open File Rep 91–193</mixed-citation>
</ref>
<ref id="CR52">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schipper</surname>
<given-names>LA</given-names>
</name>
<name>
<surname>McLeod</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>Subsidence rates and carbon loss in peat soils following conversion to pasture in the Waikato Region New Zealand</article-title>
<source>Soil Use Manag</source>
<year>2002</year>
<volume>18</volume>
<fpage>91</fpage>
<lpage>93</lpage>
<pub-id pub-id-type="doi">10.1111/j.1475-2743.2002.tb00225.x</pub-id>
</element-citation>
</ref>
<ref id="CR53">
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Schorthorst</surname>
<given-names>CJ</given-names>
</name>
</person-group>
<source>Subsidence of low moor peat soil in the western Netherlands</source>
<year>1977</year>
<publisher-loc>Wageningen, The Netherlands</publisher-loc>
<publisher-name>Institute of Land and Water Management Research, Technical Bull 102</publisher-name>
<fpage>265</fpage>
<lpage>291</lpage>
</element-citation>
</ref>
<ref id="CR54">
<mixed-citation publication-type="other">Schultz HB, Carlton AB (1959) Field windbreaks for row crops. California Agriculture: November</mixed-citation>
</ref>
<ref id="CR55">
<mixed-citation publication-type="other">Schultz HB, Carlton AB, Lory F (1963) Interpolating methods for wind erosion protection in San Joaquin asparagus. California Agric September 1963</mixed-citation>
</ref>
<ref id="CR56">
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Shlemon</surname>
<given-names>RJ</given-names>
</name>
<name>
<surname>Begg</surname>
<given-names>EL</given-names>
</name>
</person-group>
<person-group person-group-type="editor">
<name>
<surname>Suggate</surname>
<given-names>RP</given-names>
</name>
<name>
<surname>Cressel</surname>
<given-names>MM</given-names>
</name>
</person-group>
<article-title>Late Quaternary evolution of the Sacramento-San Joaquin Delta, California</article-title>
<source>Quaternary studies. Bull 13</source>
<year>1975</year>
<publisher-loc>Wellington, New Zealand</publisher-loc>
<publisher-name>The Royal Society of New Zealand</publisher-name>
<fpage>259</fpage>
<lpage>266</lpage>
</element-citation>
</ref>
<ref id="CR57">
<mixed-citation publication-type="other">Stephens JC, Stewart EH (1976) Effect of climate on organic soil subsidence. IAHS Publ 121, IAHS, Wallingford, UK, pp 649–655</mixed-citation>
</ref>
<ref id="CR58">
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Stephens</surname>
<given-names>JC</given-names>
</name>
<name>
<surname>Allen</surname>
<given-names>LH</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>E</given-names>
</name>
</person-group>
<person-group person-group-type="editor">
<name>
<surname>Holzer</surname>
<given-names>TL</given-names>
</name>
</person-group>
<article-title>Organic soil subsidence</article-title>
<source>Man-induced land subsidence. Reviews in Engineering Geology, vol VI</source>
<year>1984</year>
<publisher-loc>Boulder, CO</publisher-loc>
<publisher-name>Geological Society of America</publisher-name>
</element-citation>
</ref>
<ref id="CR59">
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Strack</surname>
<given-names>M</given-names>
</name>
</person-group>
<source>Peatlands and climate change</source>
<year>2008</year>
<publisher-loc>Parkano, Finland</publisher-loc>
<publisher-name>International Peat Society</publisher-name>
</element-citation>
</ref>
<ref id="CR60">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Syvitski</surname>
<given-names>JP</given-names>
</name>
<name>
<surname>Kettner</surname>
<given-names>AJ</given-names>
</name>
<name>
<surname>Overeem</surname>
<given-names>I</given-names>
</name>
<name>
<surname>Hutton</surname>
<given-names>EW</given-names>
</name>
<name>
<surname>Hannon</surname>
<given-names>MT</given-names>
</name>
<name>
<surname>Brakenridge</surname>
<given-names>GR</given-names>
</name>
<name>
<surname>Day</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Vörösmarty</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Saito</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Giosan</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Nicholls</surname>
<given-names>RJ</given-names>
</name>
</person-group>
<article-title>Sinking deltas due to human activities</article-title>
<source>Nat Geosci</source>
<year>2009</year>
<volume>2</volume>
<fpage>681</fpage>
<lpage>686</lpage>
<pub-id pub-id-type="doi">10.1038/ngeo629</pub-id>
</element-citation>
</ref>
<ref id="CR61">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tate</surname>
<given-names>RL</given-names>
<suffix>III</suffix>
</name>
</person-group>
<article-title>Effect of flooding on microbial activities in organic soil: carbon metabolism</article-title>
<source>Soil Sci</source>
<year>1979</year>
<volume>128</volume>
<issue>5</issue>
<fpage>267</fpage>
<lpage>273</lpage>
<pub-id pub-id-type="doi">10.1097/00010694-197911000-00002</pub-id>
</element-citation>
</ref>
<ref id="CR62">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tate</surname>
<given-names>RL</given-names>
<suffix>III</suffix>
</name>
</person-group>
<article-title>Microbial oxidation of organic matter of histosols</article-title>
<source>Adv Microb Ecol</source>
<year>1980</year>
<volume>4</volume>
<fpage>169</fpage>
<lpage>201</lpage>
<pub-id pub-id-type="doi">10.1007/978-1-4615-8291-5_5</pub-id>
</element-citation>
</ref>
<ref id="CR63">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tate</surname>
<given-names>RL</given-names>
<suffix>III</suffix>
</name>
</person-group>
<article-title>Effect of several environmental parameters on carbon metabolism in histosols</article-title>
<source>Adv Microb Ecol</source>
<year>1980</year>
<volume>5</volume>
<issue>4</issue>
<fpage>329</fpage>
<lpage>336</lpage>
<pub-id pub-id-type="doi">10.1007/BF02020340</pub-id>
</element-citation>
</ref>
<ref id="CR64">
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
<surname>Terzaghi</surname>
<given-names>K</given-names>
</name>
</person-group>
<source>Principles of soil mechanics. IV: settlement and consolidation of clay. Engineering News Record</source>
<year>1925</year>
<publisher-loc>New York</publisher-loc>
<publisher-name>McGraw–Hill</publisher-name>
</element-citation>
</ref>
<ref id="CR65">
<mixed-citation publication-type="other">Thompson J (1957) The settlement geography of the Sacramento-San Joaquin Delta, California. PhD Thesis, Stanford University, Palo Alto, CA, USA</mixed-citation>
</ref>
<ref id="CR66">
<mixed-citation publication-type="other">Tugel AJ (1993) Soil survey of Sacramento County, California. USDA, Soil Conservation Service, University of California, Agricultural Experiment Station, Davis, CA</mixed-citation>
</ref>
<ref id="CR67">
<mixed-citation publication-type="other">US Department of Agriculture (USDA) (2006) Soil Survey Geographic (SSURGO) database for Solano County, California, USDA, Natural Resources Conservation Service, Davis, CA.
<ext-link ext-link-type="uri" xlink:href="http://websoilsurvey.nrcs.usda.gov/">http://websoilsurvey.nrcs.usda.gov/</ext-link>
. Accessed 21 December 2015</mixed-citation>
</ref>
<ref id="CR68">
<mixed-citation publication-type="other">US Department of Agriculture (USDA) (2007) Soil Survey Geographic (SSURGO) database for Yolo County, California. USDA, Natural Resources Conservation Service, Davis, CA.
<ext-link ext-link-type="uri" xlink:href="http://websoilsurvey.nrcs.usda.gov/">http://websoilsurvey.nrcs.usda.gov/</ext-link>
. Accessed 21 December 2015</mixed-citation>
</ref>
<ref id="CR69">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Volk</surname>
<given-names>BG</given-names>
</name>
</person-group>
<article-title>Everglades histosol subsidence: 1. CO2 evolution as affected by soil type, temperature, and moisture</article-title>
<source>Soil Crop Sci Soc Fla Proc</source>
<year>1973</year>
<volume>32</volume>
<fpage>132</fpage>
<lpage>135</lpage>
</element-citation>
</ref>
<ref id="CR70">
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Weir</surname>
<given-names>WW</given-names>
</name>
</person-group>
<article-title>Subsidence of peat lands of the Sacramento-San Joaquin Delta, California</article-title>
<source>Hilgardia</source>
<year>1950</year>
<volume>20</volume>
<issue>3</issue>
<fpage>37</fpage>
<lpage>56</lpage>
<pub-id pub-id-type="doi">10.3733/hilg.v20n03p037</pub-id>
</element-citation>
</ref>
<ref id="CR71">
<mixed-citation publication-type="other">Welch LE (1977) Soil survey of Contra Costa County, California. USDA, Soil Conservation Service, University of California, Agricultural Experiment Station, Davis, CA</mixed-citation>
</ref>
<ref id="CR72">
<mixed-citation publication-type="other">Zanello F, Teatini P, Putti M, Gambolati G (2011) Long-time peatland subsidence: experimental study and modeling scenarios in the Venice Coastland. J Geophys Res 116(F4). doi:10.1029/2011JF002010</mixed-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>

Pour manipuler ce document sous Unix (Dilib)

EXPLOR_STEP=$WICRI_ROOT/Wicri/Agronomie/explor/SisAgriV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000055 | SxmlIndent | more

Ou

HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000055 | SxmlIndent | more

Pour mettre un lien sur cette page dans le réseau Wicri

{{Explor lien
   |wiki=    Wicri/Agronomie
   |area=    SisAgriV1
   |flux=    Pmc
   |étape=   Corpus
   |type=    RBID
   |clé=     PMC:4944668
   |texte=   Present-day oxidative subsidence of organic soils and mitigation in the Sacramento-San Joaquin Delta, California, USA
}}

Pour générer des pages wiki

HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i   -Sk "pubmed:27471427" \
       | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd   \
       | NlmPubMed2Wicri -a SisAgriV1 

Wicri

This area was generated with Dilib version V0.6.28.
Data generation: Wed Mar 29 00:06:34 2017. Site generation: Tue Mar 12 12:44:16 2024