<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article
  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article" dtd-version="3.0" xml:lang="en">
<front>
<journal-meta>
<journal-id journal-id-type="nlm-ta">PLoS ONE</journal-id>
<journal-id journal-id-type="publisher-id">plos</journal-id>
<journal-id journal-id-type="pmc">plosone</journal-id><journal-title-group>
<journal-title>PLoS ONE</journal-title></journal-title-group>
<issn pub-type="epub">1932-6203</issn>
<publisher>
<publisher-name>Public Library of Science</publisher-name>
<publisher-loc>San Francisco, USA</publisher-loc></publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">PONE-D-13-32183</article-id>
<article-id pub-id-type="doi">10.1371/journal.pone.0088463</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Agriculture</subject><subj-group><subject>Agricultural production</subject><subj-group><subject>Environmental impacts</subject></subj-group></subj-group><subj-group><subject>Crops</subject><subj-group><subject>Crop management</subject></subj-group></subj-group><subj-group><subject>Sustainable agriculture</subject></subj-group></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Biology</subject><subj-group><subject>Ecology</subject><subj-group><subject>Agroecology</subject><subject>Ecological economics</subject></subj-group></subj-group><subj-group><subject>Population biology</subject><subj-group><subject>Epidemiology</subject><subj-group><subject>Economic epidemiology</subject></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Earth sciences</subject><subj-group><subject>Geography</subject><subj-group><subject>Economic geography</subject></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Social and behavioral sciences</subject><subj-group><subject>Economics</subject><subj-group><subject>Human capital</subject><subj-group><subject>Economics of training and education</subject></subj-group></subj-group><subj-group><subject>Agricultural economics</subject><subject>Economic models</subject><subject>Economic geography</subject><subject>Environmental economics</subject></subj-group></subj-group><subj-group><subject>Sociology</subject><subj-group><subject>Social research</subject></subj-group></subj-group></subj-group></article-categories>
<title-group>
<article-title>An Integrated Framework for Assessing Vulnerability to Climate Change and Developing Adaptation Strategies for Coffee Growing Families in Mesoamerica</article-title>
<alt-title alt-title-type="running-head">Vulnerability to Climate Change in Coffee Growing</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Baca</surname><given-names>María</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib>
<contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Läderach</surname><given-names>Peter</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib>
<contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Haggar</surname><given-names>Jeremy</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib>
<contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Schroth</surname><given-names>Götz</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib>
<contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Ovalle</surname><given-names>Oriana</given-names></name><xref ref-type="aff" rid="aff4"><sup>4</sup></xref></contrib>
</contrib-group>
<aff id="aff1"><label>1</label><addr-line>International Center for Tropical Agriculture (CIAT), Managua, Nicaragua</addr-line></aff>
<aff id="aff2"><label>2</label><addr-line>Natural Resources Institute (NRI), University of Greenwich, Chatham Maritime, Kent, United Kingdom</addr-line></aff>
<aff id="aff3"><label>3</label><addr-line>Rainforest Alliance, Wageningen, The Netherlands</addr-line></aff>
<aff id="aff4"><label>4</label><addr-line>International Center for Tropical Agriculture (CIAT), Cali, Colombia</addr-line></aff>
<contrib-group>
<contrib contrib-type="editor" xlink:type="simple"><name name-style="western"><surname>Bond-Lamberty</surname><given-names>Ben</given-names></name>
<role>Editor</role>
<xref ref-type="aff" rid="edit1"/></contrib>
</contrib-group>
<aff id="edit1"><addr-line>DOE Pacific Northwest National Laboratory, United States of America</addr-line></aff>
<author-notes>
<corresp id="cor1">* E-mail: <email xlink:type="simple">magubago@gmail.com</email></corresp>
<fn fn-type="conflict"><p>Funding from Green Mountain Coffee Roasters (GMCR). Funds were used for employment and development of research. There are no further patents, products in development or marketed products to declare. This does not alter the authors′ adherence to all the PLOS ONE policies on sharing data and materials, as detailed online in the guide for authors.</p></fn>
<fn fn-type="con"><p>Conceived and designed the experiments: MB PL JH. Performed the experiments: MB PL JH OO. Analyzed the data: MB PL JH GS. Contributed reagents/materials/analysis tools: MB PL JH OO. Wrote the paper: MB GS PL JH.</p></fn>
</author-notes>
<pub-date pub-type="collection"><year>2014</year></pub-date>
<pub-date pub-type="epub"><day>26</day><month>2</month><year>2014</year></pub-date>
<volume>9</volume>
<issue>2</issue>
<elocation-id>e88463</elocation-id>
<history>
<date date-type="received"><day>11</day><month>7</month><year>2013</year></date>
<date date-type="accepted"><day>7</day><month>1</month><year>2014</year></date>
</history>
<permissions>
<copyright-year>2014</copyright-year>
<copyright-holder>Baca et al</copyright-holder><license xlink:type="simple"><license-p>This is an open-access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/" xlink:type="simple">Creative Commons Attribution License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p></license></permissions>
<abstract>
<p>The Mesoamerican region is considered to be one of the areas in the world most vulnerable to climate change. We developed a framework for quantifying the vulnerability of the livelihoods of coffee growers in Mesoamerica at regional and local levels and identify adaptation strategies. Following the Intergovernmental Panel on Climate Change (IPCC) concepts, vulnerability was defined as the combination of exposure, sensitivity and adaptive capacity. To quantify exposure, changes in the climatic suitability for coffee and other crops were predicted through niche modelling based on historical climate data and locations of coffee growing areas from Mexico, Guatemala, El Salvador and Nicaragua. Future climate projections were generated from 19 Global Circulation Models. Focus groups were used to identify nine indicators of sensitivity and eleven indicators of adaptive capacity, which were evaluated through semi-structured interviews with 558 coffee producers. Exposure, sensitivity and adaptive capacity were then condensed into an index of vulnerability, and adaptation strategies were identified in participatory workshops. Models predict that all target countries will experience a decrease in climatic suitability for growing Arabica coffee, with highest suitability loss for El Salvador and lowest loss for Mexico. High vulnerability resulted from loss in climatic suitability for coffee production and high sensitivity through variability of yields and out-migration of the work force. This was combined with low adaptation capacity as evidenced by poor post harvest infrastructure and in some cases poor access to credit and low levels of social organization. Nevertheless, the specific contributors to vulnerability varied strongly among countries, municipalities and families making general trends difficult to identify. Flexible strategies for adaption are therefore needed. Families need the support of government and institutions specialized in impacts of climate change and strengthening of farmer organizations to enable the adjustment of adaptation strategies to local needs and conditions.</p>
</abstract>
<funding-group><funding-statement>This research was conducted under the CGIAR Research Program on Climate Change, Agriculture and Food Security (CCAFS) in collaboration with Catholic Relief Services (CRS) and additional funding from Green Mountain Coffee Roasters (GMCR). Funds were used for employment and development of research. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. There are no further patents, products in development or marketed products to declare.</funding-statement></funding-group><counts><page-count count="11"/></counts></article-meta>
</front>
<body><sec id="s1">
<title>Introduction</title>
<p>Climate change represents a serious threat for Mesoamerican countries due to the multiple impacts predicted to directly affect the population as well as various sectors of the economy <xref ref-type="bibr" rid="pone.0088463-DFID1">[1]</xref>, <xref ref-type="bibr" rid="pone.0088463-IPCC1">[2]</xref>. The Vulnerability-Resilience Indicators developed by Yohe suggested high exposure to climate change for the Mesoamerican and Caribbean Region <xref ref-type="bibr" rid="pone.0088463-Yohe1">[14]</xref>. Climate projections indicate that increases in temperature will reduce crop yields in general <xref ref-type="bibr" rid="pone.0088463-IPCC2">[3]</xref> and particularly those of Arabica coffee, one of the region's major exports. Arabica coffee responds strongly to seasonal temperature patterns and coffee of high quality requires relatively stable temperatures within a fairly narrow range <xref ref-type="bibr" rid="pone.0088463-Gay1">[4]</xref>, <xref ref-type="bibr" rid="pone.0088463-Baker1">[5]</xref>. Läderach and others <xref ref-type="bibr" rid="pone.0088463-Lderach1">[6]</xref> predicted that optimal conditions for growing Arabica coffee in Mesoamerica will move from currently 800 to 1,400 m.a.s.l upwards to 1,200 to 1,600 m.a.s.l by 2050. Studies in Ethiopia and Kenya have similarly foreseen significant impacts of climate change on the distribution of wild coffee and coffee pests <xref ref-type="bibr" rid="pone.0088463-Davis1">[7]</xref>, <xref ref-type="bibr" rid="pone.0088463-Jaramillo1">[8]</xref>, in the later case extending the areas affected by coffee berry-borer.</p>
<p>During the last 40 years, agriculture has contributed 10% of the GDP in Latin American countries and is a major export earner. It is an important sector in the regional economy, since it employs 30% to 40% of the economically active population and is essential for the food security of the poorest segment of society <xref ref-type="bibr" rid="pone.0088463-IPCC2">[3]</xref>, <xref ref-type="bibr" rid="pone.0088463-Tucker1">[9]</xref>. Across Mexico and Central America, over 4 million people depend directly on coffee production for their livelihoods <xref ref-type="bibr" rid="pone.0088463-CEPAL1">[10]</xref>. According to CEPAL <xref ref-type="bibr" rid="pone.0088463-Tucker1">[9]</xref>, coffee production, purchasing, and processing employ an estimated 8.5 million people in the region. Employment and income generation from coffee are particularly significant for many indigenous peoples in Mexico and Guatemala. The environmental services generated by shade coffee farms, including carbon sequestration, watershed services and the conservation of biodiversity, have also been highlighted by many authors <xref ref-type="bibr" rid="pone.0088463-Moguel1">[11]</xref>, <xref ref-type="bibr" rid="pone.0088463-Escamilla1">[12]</xref>, while Haggar and others <xref ref-type="bibr" rid="pone.0088463-Haggar1">[13]</xref> have shown how land-use change can affect the provision of environmental services when shade coffee is replaced by other land-uses <xref ref-type="bibr" rid="pone.0088463-Haggar1">[13]</xref>.</p>
<p>The Intergovernmental Panel on Climate Change (IPCC) presents an integrated concept where “vulnerability to climate change is the degree by which a system is susceptible or unable to face the adverse effects of climate change, including climate extremes and variability. Moreover, vulnerability depends on the nature, magnitude, and rate of climate change, as well as the variation to which a system is exposed, its sensitivity and its capacity for adaptation”. Exposure is the nature and extent of changes that a place’s climate is subjected to with regard to variables such as temperature, precipitation, and extreme weather events. Sensitivity is a measure of how systems could be affected by the change in climate (e.g. how much crop yields change or how much human health might be affected). In contrast, adaptive capacity is defined as a system’s ability to adjust to climate change in order to reduce or mitigate possible damage <xref ref-type="bibr" rid="pone.0088463-IPCC2">[3]</xref>. Adaptive capacity is dynamic, and depends partly on the societý productive base, such as: natural and artificial assets, social benefits and networks, human capital and institutions, governance, national income, health and technology <xref ref-type="bibr" rid="pone.0088463-IPCC1">[2]</xref>, and how much capability a society has to adapt to the changes so as to maintain, minimize loss of, or maximize gain in welfare.</p>
<p>The current study was conducted to evaluate the vulnerability of coffee farming communities in El Salvador, Guatemala, Mexico and Nicaragua and to identify adaptation strategies to climate change <xref ref-type="bibr" rid="pone.0088463-DFID1">[1]</xref>, <xref ref-type="bibr" rid="pone.0088463-IPCC2">[3]</xref>.</p>
</sec><sec id="s2" sec-type="materials|methods">
<title>Materials and Methods</title>
<p>In order to assess the vulnerability to climate change and define appropriate adaptation strategies, we adapted the IPCC's definition of vulnerability and applied it to small coffee producers <xref ref-type="bibr" rid="pone.0088463-IPCC2">[3]</xref>. For our methodology, vulnerability is defined as changes in climate variables that affect agricultural and natural systems over a timeframe. The vulnerability in the livelihoods of small coffee farmers is a function of three factors: exposure, sensitivity and adaptive capacity.</p>
<p>These factors are related to the interaction between climate change and access and availability of resources to farming families. Exposure is quantified by modelled coffee crop suitability change comparing current and future climates, representing how familieś livelihoods will be impacted by climate change. Sensitivity and adaptive capacity are measured by indicators based on family resources-such as-natural, human, social, physical and financial capital <xref ref-type="bibr" rid="pone.0088463-DFID1">[1]</xref>. We quantified vulnerability levels by combining exposure, sensitivity and adaptive capacity. Then, we identified adaptation strategies based on vulnerability levels applying participatory methods with coffee producing communities and organizations (<xref ref-type="fig" rid="pone-0088463-g001">Figure 1</xref>). The communities and organizations included were from four countries. The study was part of a development project seeking to facilitate adaptation to climate change among coffee producers that was implemented in Mexico, Nicaragua, Guatemala and El Salvador. These countries represent a range of economic and social development in the region from Nicaragua the poorest to Mexico the richest country.</p>
<fig id="pone-0088463-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0088463.g001</object-id><label>Figure 1</label><caption>
<title>Framework to assess the vulnerability of coffee communities and to identify strategies for adaptation to climate change.</title>
</caption><graphic mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0088463.g001" position="float" xlink:type="simple"/></fig><sec id="s2a">
<title>Exposure</title>
<p>To quantify exposure to climate change, crop suitability models predicting future changes of climatic suitability of coffee were used for the four countries. The methodology combined current climate data with future climate change predictions. To map current climatic suitability, the historical climate database WorldClim (<ext-link ext-link-type="uri" xlink:href="http://www.worldclim.org" xlink:type="simple">www.worldclim.org</ext-link>) was used. The variables included a total of 19 bioclimatic variables derived from monthly precipitation, monthly median temperature, minimum and maximum temperature <xref ref-type="bibr" rid="pone.0088463-Hijmans1">[15]</xref>. Bioclimatic variables represent annual trends, seasonality, and extreme conditions.</p>
<p>To predict future climate, the SRES-A2a scenario 19 IPCC Global Circulation Models were used. The Delta method was used to down-scale the climate change data, based on the sum of the anomalies interpolated with the WorldClim monthly high-resolution surfaces <xref ref-type="bibr" rid="pone.0088463-Hijmans1">[15]</xref>. The method produces a softened surface (interpolation) of climate change (deltas or anomalies). It implies that changes in climate are only relevant at coarse-scale and that the interactions between variables are maintained in the future <xref ref-type="bibr" rid="pone.0088463-Ramrez1">[16]</xref>.</p>
<p>The Maximum entropy (MAXENT) method, a general-purpose method for making predictions or inferences based on incomplete information <xref ref-type="bibr" rid="pone.0088463-Phillips1">[17]</xref>, was used to predict the future climatic suitability for coffee. The model requires calibration with climate data for current coffee production areas, which is provided by GPS coordinates. The model assumes that a certain future climate at a given site is as suitable or unsuitable for the crop as is the same climate at another site in the present. This assumption is reasonable as long as crop genetics and cropping systems do not significantly change. It thus predicts what will happen in terms of relative climatic suitability for a crop if these factors do not change and helps identify those sites where adaptations in crops and cropping systems are necessary in order to avoid the consequences of a predicted decline in climatic suitability. This approach has previously been used for coffee <xref ref-type="bibr" rid="pone.0088463-Lderach1">[6]</xref>, <xref ref-type="bibr" rid="pone.0088463-Schroth1">[18]</xref>.</p>
<p>Two measures of uncertainty were calculated: (1) the agreement of calculated models as a percentage of models that predict changes in the same direction and (2) the coefficient of variation (CV) among models.</p>
</sec><sec id="s2b">
<title>Sensitivity and adaptive capacity</title>
<p>Indicators of the sensitivity to climate change and adaptive capacity were devised in collaboration with organizations and experts from the region using an expert panel, focus groups, and semi-structured interviews. For the expert panel, semi-structured individual interviews were conducted with 17 key informants of the coffee sector in Nicaragua, including technicians, farmers and researchers. It included questions about the most important factors affecting coffee production. Four focus groups were carried out in Nicaragua and three groups in each of the remaining countries (El Salvador, Guatemala and Mexico). Participants discussed and assessed the significance of climate change over time and identified key indicators for coffee livelihoods. The list of key indicators was structured according to the five community capitals (natural, human, social, physical and financial) of the Livelihoods Approach <xref ref-type="bibr" rid="pone.0088463-DFID1">[1]</xref>.</p>
<p>Parameters were then constructed to evaluate each indicator as shown in <xref ref-type="supplementary-material" rid="pone.0088463.s001">Table S1 in File S1</xref>. To quantify the parameters scales from 1 to 5 were applied or a binary scale of 0 and 1, depending on the nature of the parameter. The final values for each indicator were calculated by averaging all the parameters and then transformed to a 0-1 continuous variable scale, with 0 being low and 1 being high sensitivity and adaptive capacity. For example, access to and availability of water is an important natural resource for familieś livelihoods and coffee production. To measure the water access and availability indicator we considered the parameters source, distance, quality and quantity of water. Water availability, for example, is measured on a scale of 1 to 5, 1 being least sensitive and 5 being most sensitive. A value of 1 means there is never sufficient water and a value of 5 means there is an abundance of water all year. Then we developed a semi-structured interview by adapting qualitative tools <xref ref-type="bibr" rid="pone.0088463-Geilfus1">[19]</xref> and validated the tool with six coffee families.</p>
<p>The indicators were used to assess the vulnerability of coffee farms in each country. From a population of 7,000 farmer members from 15 organizations across the four countries, 558 farmers were interviewed. The farmers may be considered representative of small-scale organized farmers, but should not be considered representative of the coffee farmers as a whole in each country. The sample size was defined using the formula for finite populations <xref ref-type="bibr" rid="pone.0088463-Surez1">[20]</xref> and then individual farmers were selected randomly, stratified according to exposure level and country by 2050 (<xref ref-type="table" rid="pone-0088463-t001">Table 1</xref>).</p>
<table-wrap id="pone-0088463-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0088463.t001</object-id><label>Table 1</label><caption>
<title>Number of interviewed families by country and exposure level<xref ref-type="table-fn" rid="nt101">a</xref>.</title>
</caption><alternatives><graphic id="pone-0088463-t001-1" position="float" mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0088463.t001" xlink:type="simple"/>
<table><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup>
<thead>
<tr>
<td align="left" rowspan="1" colspan="1">Country</td>
<td align="left" rowspan="1" colspan="1">Department or State</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Exposure level</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Total</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">High</td>
<td align="left" rowspan="1" colspan="1">Medium</td>
<td align="left" rowspan="1" colspan="1">Low</td>
<td align="left" rowspan="1" colspan="1"/>
</tr>
</thead>
<tbody>
<tr>
<td align="left" rowspan="1" colspan="1">Nicaragua</td>
<td align="left" rowspan="1" colspan="1">Jinotega</td>
<td align="left" rowspan="1" colspan="1">12</td>
<td align="left" rowspan="1" colspan="1">14</td>
<td align="left" rowspan="1" colspan="1">15</td>
<td align="left" rowspan="1" colspan="1">41</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Matagalpa</td>
<td align="left" rowspan="1" colspan="1">36</td>
<td align="left" rowspan="1" colspan="1">20</td>
<td align="left" rowspan="1" colspan="1">5</td>
<td align="left" rowspan="1" colspan="1">61</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Madriz</td>
<td align="left" rowspan="1" colspan="1">4</td>
<td align="left" rowspan="1" colspan="1">14</td>
<td align="left" rowspan="1" colspan="1">14</td>
<td align="left" rowspan="1" colspan="1">32</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Nueva Segovia</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">16</td>
<td align="left" rowspan="1" colspan="1">16</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Total</td>
<td align="left" rowspan="1" colspan="1">50</td>
<td align="left" rowspan="1" colspan="1">50</td>
<td align="left" rowspan="1" colspan="1">50</td>
<td align="left" rowspan="1" colspan="1">150</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">El Salvador</td>
<td align="left" rowspan="1" colspan="1">Usulután</td>
<td align="left" rowspan="1" colspan="1">11</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">9</td>
<td align="left" rowspan="1" colspan="1">20</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Santa Ana</td>
<td align="left" rowspan="1" colspan="1">8</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">8</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">La Libertad</td>
<td align="left" rowspan="1" colspan="1">4</td>
<td align="left" rowspan="1" colspan="1">11</td>
<td align="left" rowspan="1" colspan="1">14</td>
<td align="left" rowspan="1" colspan="1">29</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Ahuachapán</td>
<td align="left" rowspan="1" colspan="1">20</td>
<td align="left" rowspan="1" colspan="1">32</td>
<td align="left" rowspan="1" colspan="1">20</td>
<td align="left" rowspan="1" colspan="1">72</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Total</td>
<td align="left" rowspan="1" colspan="1">43</td>
<td align="left" rowspan="1" colspan="1">43</td>
<td align="left" rowspan="1" colspan="1">43</td>
<td align="left" rowspan="1" colspan="1">129</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Guatemala</td>
<td align="left" rowspan="1" colspan="1">Chiquimula</td>
<td align="left" rowspan="1" colspan="1">14</td>
<td align="left" rowspan="1" colspan="1">23</td>
<td align="left" rowspan="1" colspan="1">10</td>
<td align="left" rowspan="1" colspan="1">47</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Sololá</td>
<td align="left" rowspan="1" colspan="1">17</td>
<td align="left" rowspan="1" colspan="1">2</td>
<td align="left" rowspan="1" colspan="1">4</td>
<td align="left" rowspan="1" colspan="1">23</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Chimaltenango</td>
<td align="left" rowspan="1" colspan="1">4</td>
<td align="left" rowspan="1" colspan="1">5</td>
<td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">10</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">San Marcos</td>
<td align="left" rowspan="1" colspan="1">8</td>
<td align="left" rowspan="1" colspan="1">13</td>
<td align="left" rowspan="1" colspan="1">28</td>
<td align="left" rowspan="1" colspan="1">49</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Total</td>
<td align="left" rowspan="1" colspan="1">43</td>
<td align="left" rowspan="1" colspan="1">43</td>
<td align="left" rowspan="1" colspan="1">43</td>
<td align="left" rowspan="1" colspan="1">129</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">México</td>
<td align="left" rowspan="1" colspan="1">Chiapas</td>
<td align="left" rowspan="1" colspan="1">4</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">32</td>
<td align="left" rowspan="1" colspan="1">36</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Oaxaca</td>
<td align="left" rowspan="1" colspan="1">45</td>
<td align="left" rowspan="1" colspan="1">30</td>
<td align="left" rowspan="1" colspan="1">11</td>
<td align="left" rowspan="1" colspan="1">86</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">San Luis Potosí</td>
<td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">20</td>
<td align="left" rowspan="1" colspan="1">7</td>
<td align="left" rowspan="1" colspan="1">28</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Total</td>
<td align="left" rowspan="1" colspan="1">50</td>
<td align="left" rowspan="1" colspan="1">50</td>
<td align="left" rowspan="1" colspan="1">50</td>
<td align="left" rowspan="1" colspan="1">150</td>
</tr>
</tbody>
</table>
</alternatives><table-wrap-foot><fn id="nt101"><label>a</label><p>For the definition of exposure levels see section 2.3.</p></fn></table-wrap-foot></table-wrap></sec><sec id="s2c">
<title>Vulnerability</title>
<p>For exposure, the relative decreases in climatic suitability according to the MAXENT model were divided into three classes of suitability loss (low, medium, high). For sensitivity and adaptive capacity, indicators were identified and quantified through interviews with the farming families.</p>
<p>A cluster analysis was carried out for each indicator of sensitivity and adaptive capacity based on the score of each family using the Ward method with Euclidean distance. Then an Analysis of Variance (ANOVA) was applied using the LSD-Fisher test to compare the averages for each indicator by cluster. The indicators in each cluster that obtained significantly different sample averages were classified in three levels on a scale of 0 to 1 (0–0.33 = low, 0.34–0.66 = medium, 0.67–1 = high). Clusters with the greatest number of indicators with high, medium or low averages were classified as having high, medium or low sensitivity and adaptive capacity <xref ref-type="bibr" rid="pone.0088463-DiRienzo1">[21]</xref>.</p>
<p>Each factor (exposure, sensitivity and adaptive capacity), as previously explained, and was classified into three levels (high, medium, low). To calculate the vulnerability equation we assigned each level a quantitative value: low = 1, medium = 2, high = 3. With three factors and three levels per factor, we obtained 27 possible combinations. After applying the equation we obtained 7 values (–1,0,1,2,3,4,5), which we used to define low (–1,0), medium (1,2,3,) and high (4,5) levels of vulnerability (<xref ref-type="fig" rid="pone-0088463-g001">Figure 1</xref>). A Principal Components Analysis (PCA) was carried out to identify the indicators that most contribute to the sensitivity or adaptive capacity of families in different municipalities.</p>
</sec><sec id="s2d">
<title>Identifying adaptation strategies</title>
<p>Workshops were carried out to identify possible adaptation strategies. Participants were families of farmers, technicians, and presidents of different cooperatives. Workshops began by presenting the results of the vulnerability analysis by state, department, or municipality. After a general discussion, the following question was asked to groups: Given that the conditions for coffee production will change by 2050, what can be done to maintain the level of production? This was first discussed in the plenary and then in subgroups of 3 to 5 people, followed by the presentation of results. Each participant then noted the three most important ideas from the discussion on separate cards and assigned them to one of the five types of resources (natural, human, social, physical and financial) where he felt they most closely linked. Then subgroups were formed for each resource type to organize the ideas and build an outline of a climate change adaptation strategy. In this strategy, key actors, roles, resource availability, and time needed to implement the strategy were identified. The subgroups presented their results to the entire group and received feedback, until a general consensus was reached for each adaptation strategy.</p>
</sec></sec><sec id="s3">
<title>Results</title>
<sec id="s3a">
<title>Exposure</title>
<p>According to the climate change models, total annual precipitation in all countries is predicted to decrease by 2050. Nicaragua, already the driest of the four countries, is predicted to have a 5% precipitation decrease by 2050. The mean annual temperature will increase by 2.3°C, while the maximum temperature will increase in all countries by between 2.2°C and 2.4°C (<xref ref-type="supplementary-material" rid="pone.0088463.s001">Table S2 in File S1</xref>). The most decisive climatic variables for the predicted decrease in climatic suitability for coffee were the increase in maximum mean temperature of the hottest month. Mexico is predicted to reach a mean maximum temperature of 36°C by 2050.</p>
<p>In interviews the producers confirmed the trend of the climate models (<xref ref-type="supplementary-material" rid="pone.0088463.s001">Table S2 in File S1</xref>). They mentioned changes in climate seasonality and predictability, including hotter and longer dry seasons (from three-four months to four-six months) and shorter and drier wet seasons. Also, they perceived an increase in extreme temperatures, drought and wind as well as changes in the intensity and distribution of precipitation. They also highlighted an increase in extreme weather events, such as cold periods, hail, drought, and hurricanes.</p>
<p>The climatic suitability for Arabica coffee has been predicted to decrease in the lowest altitude areas as a result of the increase in temperature to which coffee quality is sensitive <xref ref-type="bibr" rid="pone.0088463-Castro1">[22]</xref>. Changes in temperature and rainfall will decrease the area suitable for coffee and effectively displace coffee up the altitudinal gradient to cooler climates. On a national level, El Salvador and Nicaragua have the highest percentage of land affected by decreases in suitability of 40% or greater (<xref ref-type="table" rid="pone-0088463-t002">Table 2</xref>). Models predicted that, in Central America as a whole, the optimal coffee-growing elevation will shift from 1,200 m.a.s.l currently to 1,600 m.a.s.l by 2050 <xref ref-type="bibr" rid="pone.0088463-Lderach1">[6]</xref>, <xref ref-type="bibr" rid="pone.0088463-Lderach2">[23]</xref>.</p>
<table-wrap id="pone-0088463-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0088463.t002</object-id><label>Table 2</label><caption>
<title>Projected changes in overall suitability for coffee production and altitudinal range suitable for production in Mesoamerica by 2050.</title>
</caption><alternatives><graphic id="pone-0088463-t002-2" position="float" mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0088463.t002" xlink:type="simple"/>
<table><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup>
<thead>
<tr>
<td align="left" rowspan="1" colspan="1">Country</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Changes in overall suitability for coffee production</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Altitude suitable for production in meters above sea level</td>
<td align="left" rowspan="1" colspan="1"/>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">–40% or more</td>
<td align="left" rowspan="1" colspan="1">–40% to –20%</td>
<td align="left" rowspan="1" colspan="1">–20% to 0%</td>
<td align="left" rowspan="1" colspan="1">&lt;0%</td>
<td align="left" rowspan="1" colspan="1">Current model</td>
<td align="left" rowspan="1" colspan="1">Future model</td>
</tr>
</thead>
<tbody>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">(1950–2005)</td>
<td align="left" rowspan="1" colspan="1">(2050)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">El Salvador</td>
<td align="left" rowspan="1" colspan="1">45.5</td>
<td align="left" rowspan="1" colspan="1">43.7</td>
<td align="left" rowspan="1" colspan="1">10.9</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">700 to 1700</td>
<td align="left" rowspan="1" colspan="1">1000 to 1700</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Guatemala</td>
<td align="left" rowspan="1" colspan="1">12.9</td>
<td align="left" rowspan="1" colspan="1">25.5</td>
<td align="left" rowspan="1" colspan="1">54.2</td>
<td align="left" rowspan="1" colspan="1">7.4</td>
<td align="left" rowspan="1" colspan="1">600 to 1800</td>
<td align="left" rowspan="1" colspan="1">1200 to 2200</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Mexico</td>
<td align="left" rowspan="1" colspan="1">18.2</td>
<td align="left" rowspan="1" colspan="1">34.6</td>
<td align="left" rowspan="1" colspan="1">46.9</td>
<td align="left" rowspan="1" colspan="1">0.3</td>
<td align="left" rowspan="1" colspan="1">500 to 2000</td>
<td align="left" rowspan="1" colspan="1">1200 to 2300</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Nicaragua</td>
<td align="left" rowspan="1" colspan="1">35.3</td>
<td align="left" rowspan="1" colspan="1">32.1</td>
<td align="left" rowspan="1" colspan="1">32.5</td>
<td align="left" rowspan="1" colspan="1">0.1</td>
<td align="left" rowspan="1" colspan="1">700 to 1500</td>
<td align="left" rowspan="1" colspan="1">1000 to 1600</td>
</tr>
</tbody>
</table>
</alternatives><table-wrap-foot><fn id="nt102"><label/><p>Adapted from Läderach et al. (2010b).</p></fn></table-wrap-foot></table-wrap>
<p>In general, the climatic suitability for Arabica coffee will be maintained at the highest altitudes. In Nicaragua, the area with the greatest loss in suitability is located in the Pacific zone in the departments of Carazo and Managua, while the highlands around Apanás Lake in Jinotega will maintain a high suitability for coffee quality. In El Salvador, the area with loss of suitability will be the departments of San Miguel and Usulután, while the department of Ahuachapán will maintain high suitability. In Guatemala, greatest loss of suitability will be the southern slope of the Pacific volcanic chain, northern Zacapa and eastern Chiquimula, while Chimaltenango will maintain high suitability. In Mexico, areas with greatest loss of suitability will be northern Chiapas, Veracruz and Tabasco, whereas high suitability will be maintained in the Chiapas highlands (<xref ref-type="fig" rid="pone-0088463-g002">Figure 2</xref>).</p>
<fig id="pone-0088463-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0088463.g002</object-id><label>Figure 2</label><caption>
<title>Prediction of the relative climatic suitability for Arabica coffee production in Mexico, Guatemala, El Salvador and Nicaragua in 2010 and 2050 (large maps), coefficient of variation (CV; small map to the left), and consistency between models (small map to the mid-right).</title>
</caption><graphic mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0088463.g002" position="float" xlink:type="simple"/></fig></sec><sec id="s3b">
<title>Sensitivity and adaptive capacity</title>
<sec id="s3b1">
<title>Sensitivity</title>
<p>The sampled families from each country were divided into three sensitivity levels using cluster analysis; on average 33% of the families fell into the highly sensitive level. The families of El Salvador and Guatemala were more likely to cluster in the high sensitivity level, while those of Nicaragua and Mexico were more likely to cluster in the medium sensitivity class. Guatemala had the highest percentage of families with high sensitivity (49%) and El Salvador the highest percentage with low sensitivity (34%) (<xref ref-type="table" rid="pone-0088463-t003">Table 3</xref>).</p>
<table-wrap id="pone-0088463-t003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0088463.t003</object-id><label>Table 3</label><caption>
<title>Distribution of families per level sensitivity between countries.</title>
</caption><alternatives><graphic id="pone-0088463-t003-3" position="float" mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0088463.t003" xlink:type="simple"/>
<table><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup>
<thead>
<tr>
<td align="left" rowspan="1" colspan="1">Country</td>
<td align="left" rowspan="1" colspan="1">High sensitivity (%)</td>
<td align="left" rowspan="1" colspan="1">Medium sensitivity (%)</td>
<td align="left" rowspan="1" colspan="1">Low sensitivity (%)</td>
<td align="left" rowspan="1" colspan="1">Total (%)</td>
</tr>
</thead>
<tbody>
<tr>
<td align="left" rowspan="1" colspan="1">Nicaragua</td>
<td align="left" rowspan="1" colspan="1">22</td>
<td align="left" rowspan="1" colspan="1">61</td>
<td align="left" rowspan="1" colspan="1">17</td>
<td align="left" rowspan="1" colspan="1">100</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">El Salvador</td>
<td align="left" rowspan="1" colspan="1">40</td>
<td align="left" rowspan="1" colspan="1">26</td>
<td align="left" rowspan="1" colspan="1">34</td>
<td align="left" rowspan="1" colspan="1">100</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Guatemala</td>
<td align="left" rowspan="1" colspan="1">49</td>
<td align="left" rowspan="1" colspan="1">37</td>
<td align="left" rowspan="1" colspan="1">14</td>
<td align="left" rowspan="1" colspan="1">100</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Mexico</td>
<td align="left" rowspan="1" colspan="1">23</td>
<td align="left" rowspan="1" colspan="1">46</td>
<td align="left" rowspan="1" colspan="1">31</td>
<td align="left" rowspan="1" colspan="1">100</td>
</tr>
</tbody>
</table>
</alternatives></table-wrap>
<p>An ANOVA test indicated that there are significant differences (p&lt;0.001) between high, medium and low clusters for the sensitivity indicators for each country. <xref ref-type="fig" rid="pone-0088463-g003">Figure 3</xref> shows some indicators that were significantly different in the high sensitivity level compared by cluster for each country. The families in the high sensitivity level in each country were characterized by high yield variability. In El Salvador, Guatemala and Nicaragua this indicator was significantly higher for high sensitivity than for medium or low sensitivity clusters (p&lt;0.001). In contrast, in Mexico the indicator was not significantly different (p = 0.090) between sensitivity clusters (<xref ref-type="fig" rid="pone-0088463-g003">Figure 3</xref>). Yield variability leads to frequent reductions in income and the ability of the families to respond to external stresses such as climate change. Furthermore, in Nicaragua, Mexico and Guatemala the migration indicator was also significantly greater for the high sensitivity cluster in each country (p&lt;0.001). Migration, being the temporary or permanent work-related move of one or more family members to a foreign country, reduces the availability of family labour and, thus, the resilience to respond to climate change impacts. The resulting lack of labour for the coffee harvest was frequently mentioned during the interviews.</p>
<fig id="pone-0088463-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0088463.g003</object-id><label>Figure 3</label><caption>
<title>Sensitivity indicators in the livelihoods of small coffee producers to climate change in four countries of Mesoamerica (a high value equals high sensitivity).</title>
</caption><graphic mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0088463.g003" position="float" xlink:type="simple"/></fig>
<p>In El Salvador, a notable contributor to high sensitivity was lack of conservation practices (p&lt;0.001), including lack of maintenance of vegetation protecting water sources and communally or individually managed forest areas. In parts of the country, farm sizes are small (between 1–2 hectares) and therefore there is little space for leaving areas of natural vegetation. On the other hand, producers in some communities in the departments of Ahuachapán and La Libertad in El Salvador are organized in groups of 40–50 families that collectively own 600–1100 hectares of land on which they maintain water sources and forest areas. In Mexico, high sensitivity was particularly related to difficulties with transportation of products (p&lt;0.001) (<xref ref-type="fig" rid="pone-0088463-g003">Figure 3</xref>).</p>
<p>The Principal Components Analysis showed that key indicators of vulnerability differed among municipalities even within the same country. The value of the first axis and second axis per country is shown in the <xref ref-type="fig" rid="pone-0088463-g004">Figure 4</xref>. The variables that were highly correlated with axis 1 included yield variability (r = 0.48) in Nicaragua, road type (r = 0.48) in El Salvador, transport of products (r = 0.46) in Guatemala and, health and nutrition (r = 0.46) in Mexico. For example, poor transport and bus connections restricting market access was a key factor of vulnerability in San Lucas Tolimán in Sololá, Guatemala. Farmers grow coffee among the volcanoes surrounding Lake Atitlan and had to walk between 3–4 hours from their farms to the nearest market town, transporting their produce either on mules or on their backs. Similar conditions prevail in many coffee communities in remote mountain areas across the region.</p>
<fig id="pone-0088463-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0088463.g004</object-id><label>Figure 4</label><caption>
<title>Principal components analysis of association of sensitivity indicators with different municipalities.</title>
</caption><graphic mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0088463.g004" position="float" xlink:type="simple"/></fig>
<p>In El Tuma-La Dalia in Matagalpa the coffee yield of producer farms had high variability between years. Ranges in yield were between +/– 25% on average for four years, with a variability of +/– 50% reported in some cases; the average farmer yield was 331 Kg/ha, in contrast with the Matagalpan department average of 538 Kg/ha. High variability of yields and thus income between years limited the access and availability of food, nutrition, education and health of families. Sometimes when the yield was very low, in Nicaragua the families migrate to other places in search of casual labour. Depending on the municipality or community, migration was an important factor in each country and included both seasonal and permanent migrants (e.g. in Matagalpa, Nicaragua, ten out of eleven interviewed families had some of its members living abroad as migrants) (<xref ref-type="fig" rid="pone-0088463-g004">Figure 4</xref>).</p>
</sec><sec id="s3b2">
<title>Adaptive capacity</title>
<p>The sampled families from each of the four countries were divided into three adaptive capacity levels using the cluster analysis; on average, 32% of families fell into the low adaptive capacity level. The families of Nicaragua, El Salvador and Guatemala were more likely to cluster in the low adaptive capacity level (<xref ref-type="table" rid="pone-0088463-t004">Table 4</xref>).</p>
<table-wrap id="pone-0088463-t004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0088463.t004</object-id><label>Table 4</label><caption>
<title>Distribution of families per level adaptive capacity between countries.</title>
</caption><alternatives><graphic id="pone-0088463-t004-4" position="float" mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0088463.t004" xlink:type="simple"/>
<table><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup>
<thead>
<tr>
<td align="left" rowspan="1" colspan="1">Country</td>
<td align="left" rowspan="1" colspan="1">High adaptive capacity (%)</td>
<td align="left" rowspan="1" colspan="1">Medium adaptive capacity (%)</td>
<td align="left" rowspan="1" colspan="1">Low adaptive capacity (%)</td>
<td align="left" rowspan="1" colspan="1">Total (%)</td>
</tr>
</thead>
<tbody>
<tr>
<td align="left" rowspan="1" colspan="1">Nicaragua</td>
<td align="left" rowspan="1" colspan="1">41</td>
<td align="left" rowspan="1" colspan="1">22</td>
<td align="left" rowspan="1" colspan="1">37</td>
<td align="left" rowspan="1" colspan="1">100</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">El Salvador</td>
<td align="left" rowspan="1" colspan="1">50</td>
<td align="left" rowspan="1" colspan="1">15</td>
<td align="left" rowspan="1" colspan="1">35</td>
<td align="left" rowspan="1" colspan="1">100</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Guatemala</td>
<td align="left" rowspan="1" colspan="1">53</td>
<td align="left" rowspan="1" colspan="1">13</td>
<td align="left" rowspan="1" colspan="1">34</td>
<td align="left" rowspan="1" colspan="1">100</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Mexico</td>
<td align="left" rowspan="1" colspan="1">38</td>
<td align="left" rowspan="1" colspan="1">38</td>
<td align="left" rowspan="1" colspan="1">24</td>
<td align="left" rowspan="1" colspan="1">100</td>
</tr>
</tbody>
</table>
</alternatives></table-wrap>
<p>The ANOVA test indicated significant differences between the high, medium and low adaptive capacity clusters in each country, but the indicators varied across countries. <xref ref-type="fig" rid="pone-0088463-g005">Figure 5</xref> shows indicators that were significantly different in low adaptive capacity level compared to clusters for each country. The families in the low adaptive capacity cluster, in Mexico, Guatemala, El Salvador and Nicaragua were characterized by low viability of post harvest infrastructure for drying coffee. Low adaptive capacity families in Mexico and El Salvador had poorer post harvest infrastructure for drying coffee (p&lt;0.001), than higher adaptive capacity families. Between 50 and 80% of families interviewed had only a single technique for drying coffee irrespective of weather conditions. While in Mexico, El Salvador and Guatemala producers had access to drying patios, in Nicaragua producers dry their coffee using drying tables, drying patios or plastic tarps and, in Mexico producer dry coffee using sacks, or on the floors or house roofs. In very humid regions, these drying techniques were not suitable due to lack of adequate sun; as a result, the drying quality is poor. Some cooperatives in Nicaragua and Guatemala with economic resources transport the coffee to drier areas, while other cooperatives use a mechanical drying process. In both cases this increases the cost of processing.</p>
<fig id="pone-0088463-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0088463.g005</object-id><label>Figure 5</label><caption>
<title>Adaptive capacity indicators in the livelihoods of small coffee producers to climate change in four countries of Mesoamerica (a low value equals low adaptive capacity).</title>
</caption><graphic mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0088463.g005" position="float" xlink:type="simple"/></fig>
<p>Furthermore, low adaptive capacity farmers in Nicaragua presented lower income diversification (p&lt;0.005).Families depended on coffee sales for between 50 to 60% of their yearly income. Some families diversify their income with basic grains, bananas, oranges and avocados, but the intermediaries pay very little and the markets are distant. In all countries farmers received credit from cooperatives. However, in Guatemala low adaptive capacity farmers had more limited access to credit (p&lt;0.001) than higher adaptive capacity farmers, with 27% of 129 farmers interviewed having had no access to credits. In Mexico low adaptive capacity farmers had lower knowledge levels regarding coffee sector policies and environmental and land use laws (p&lt;0.001), than higher adaptive capacity famers. In general, interviewed families had knowledge of only one to three coffee sector or environmental policies, and they did not have active participation in the application of these laws. Additionally, in Mexico low adaptive capacity families had low access to alternative technologies, but this result was not significantly different from the other adaptive capacity levels (p = 0.079). For example, they did not collect water for their own consumption or for crops, and they did not use drip irrigation (<xref ref-type="fig" rid="pone-0088463-g005">Figure 5</xref>).</p>
<p>There was a close association of some municipalities with certain indicators in the Principal Components Analysis. The value of the first axis and second axis per country is shown in <xref ref-type="fig" rid="pone-0088463-g006">Figure 6</xref>. The variables that were least correlated with axis 1 included post harvest infrastructure viability (r = –0.09) in Nicaragua and (r = –0.159) Mexico, organization (r = –0.43) in El Salvador and, education (r = –0.31) in Guatemala. For example, in Xilitla in San Luis de Potosí, Mexico interviewed families had little access to post-harvest infrastructure; farmers dry their coffee on the floor, in sacks, on plastic or other forms, and some producers sell their coffee as unprocessed cherries with corresponding price reductions (<xref ref-type="fig" rid="pone-0088463-g006">Figure 6</xref>).</p>
<fig id="pone-0088463-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0088463.g006</object-id><label>Figure 6</label><caption>
<title>Principal components analysis of the association of adaptive capacity indicators to different municipalities.</title>
</caption><graphic mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0088463.g006" position="float" xlink:type="simple"/></fig>
<p>In San Lucas Tolimán in Sololá, Guatemala, the families had little access to formal and informal education. This is because the heads of household are of Mayan origin and speak their local language, rather than Spanish. They had little access to education, and there is no technical assistance in their own language. Currently children have access to primary school and Spanish language instruction.</p>
</sec></sec><sec id="s3c">
<title>Vulnerability and adaptation strategies</title>
<p>The results of the vulnerability equation indicate that in Nicaragua, of the 143 families, 18% had a high level of vulnerability; in El Salvador, 14% out of 129 interviewed families showed high vulnerability; in Guatemala, 22% out of 129 families showed high vulnerability; and in Mexico 9% out of 150 families showed high vulnerability (<xref ref-type="table" rid="pone-0088463-t005">Table 5</xref>).</p>
<table-wrap id="pone-0088463-t005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0088463.t005</object-id><label>Table 5</label><caption>
<title>Percentage of families by vulnerability level in each country.</title>
</caption><alternatives><graphic id="pone-0088463-t005-5" position="float" mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0088463.t005" xlink:type="simple"/>
<table><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup>
<thead>
<tr>
<td align="left" rowspan="1" colspan="1">Country</td>
<td align="left" rowspan="1" colspan="1">Department or state</td>
<td align="left" rowspan="1" colspan="1">Vulnerability level</td>
<td align="left" rowspan="1" colspan="1">Total (%)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">High</td>
<td align="left" rowspan="1" colspan="1">Medium</td>
<td align="left" rowspan="1" colspan="1">Low</td>
<td align="left" rowspan="1" colspan="1"/>
</tr>
</thead>
<tbody>
<tr>
<td align="left" rowspan="1" colspan="1">Nicaragua</td>
<td align="left" rowspan="1" colspan="1">Jinotega</td>
<td align="left" rowspan="1" colspan="1">6</td>
<td align="left" rowspan="1" colspan="1">13</td>
<td align="left" rowspan="1" colspan="1">10</td>
<td align="left" rowspan="1" colspan="1">29</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Matagalpa</td>
<td align="left" rowspan="1" colspan="1">9</td>
<td align="left" rowspan="1" colspan="1">22</td>
<td align="left" rowspan="1" colspan="1">6</td>
<td align="left" rowspan="1" colspan="1">38</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Madriz</td>
<td align="left" rowspan="1" colspan="1">3</td>
<td align="left" rowspan="1" colspan="1">10</td>
<td align="left" rowspan="1" colspan="1">9</td>
<td align="left" rowspan="1" colspan="1">22</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Nueva Segovia</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">6</td>
<td align="left" rowspan="1" colspan="1">5</td>
<td align="left" rowspan="1" colspan="1">11</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Total (%)</td>
<td align="left" rowspan="1" colspan="1">18</td>
<td align="left" rowspan="1" colspan="1">51</td>
<td align="left" rowspan="1" colspan="1">31</td>
<td align="left" rowspan="1" colspan="1">100</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">El Salvador</td>
<td align="left" rowspan="1" colspan="1">Usulután</td>
<td align="left" rowspan="1" colspan="1">5</td>
<td align="left" rowspan="1" colspan="1">10</td>
<td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">16</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Santa Ana</td>
<td align="left" rowspan="1" colspan="1">2</td>
<td align="left" rowspan="1" colspan="1">5</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">6</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">La Libertad</td>
<td align="left" rowspan="1" colspan="1">5</td>
<td align="left" rowspan="1" colspan="1">10</td>
<td align="left" rowspan="1" colspan="1">8</td>
<td align="left" rowspan="1" colspan="1">22</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Ahuachapán</td>
<td align="left" rowspan="1" colspan="1">3</td>
<td align="left" rowspan="1" colspan="1">43</td>
<td align="left" rowspan="1" colspan="1">9</td>
<td align="left" rowspan="1" colspan="1">56</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Total (%)</td>
<td align="left" rowspan="1" colspan="1">14</td>
<td align="left" rowspan="1" colspan="1">68</td>
<td align="left" rowspan="1" colspan="1">18</td>
<td align="left" rowspan="1" colspan="1">100</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Guatemala</td>
<td align="left" rowspan="1" colspan="1">Chiquimula</td>
<td align="left" rowspan="1" colspan="1">13</td>
<td align="left" rowspan="1" colspan="1">22</td>
<td align="left" rowspan="1" colspan="1">2</td>
<td align="left" rowspan="1" colspan="1">36</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Sololá</td>
<td align="left" rowspan="1" colspan="1">9</td>
<td align="left" rowspan="1" colspan="1">9</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">18</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Chimaltenango</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">7</td>
<td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">8</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">San Marcos</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">26</td>
<td align="left" rowspan="1" colspan="1">12</td>
<td align="left" rowspan="1" colspan="1">38</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Total (%)</td>
<td align="left" rowspan="1" colspan="1">22</td>
<td align="left" rowspan="1" colspan="1">64</td>
<td align="left" rowspan="1" colspan="1">15</td>
<td align="left" rowspan="1" colspan="1">100</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">México</td>
<td align="left" rowspan="1" colspan="1">Chiapas</td>
<td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">15</td>
<td align="left" rowspan="1" colspan="1">9</td>
<td align="left" rowspan="1" colspan="1">24</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Oaxaca</td>
<td align="left" rowspan="1" colspan="1">9</td>
<td align="left" rowspan="1" colspan="1">40</td>
<td align="left" rowspan="1" colspan="1">9</td>
<td align="left" rowspan="1" colspan="1">57</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">San Luis Potosí</td>
<td align="left" rowspan="1" colspan="1">0</td>
<td align="left" rowspan="1" colspan="1">18</td>
<td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">19</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Total (%)</td>
<td align="left" rowspan="1" colspan="1">9</td>
<td align="left" rowspan="1" colspan="1">73</td>
<td align="left" rowspan="1" colspan="1">18</td>
<td align="left" rowspan="1" colspan="1">100</td>
</tr>
</tbody>
</table>
</alternatives></table-wrap>
<p>High vulnerability was related to high and medium exposure, which was represented by a loss of climatic suitability for coffee production but in Mexico this tended to be less severe than in the three countries further south. The families’ farms located in the high exposure level will not have optimal conditions for production quality coffee by 2050 (<xref ref-type="fig" rid="pone-0088463-g007">Figure 7</xref>).</p>
<fig id="pone-0088463-g007" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0088463.g007</object-id><label>Figure 7</label><caption>
<title>Small-scale variability of vulnerability to climate change among coffee producing communities in four countries of Mesoamerica.</title>
</caption><graphic mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0088463.g007" position="float" xlink:type="simple"/></fig>
<p>In addition high vulnerability was related to high and medium sensitivity, which was due to the high variability in productivity levels. The variability of production is very important for the families because it represents their principal income (on average 50 to 65%) and they depend on this income for food security, health and education. Furthermore, Nicaragua, Guatemala and Mexico had high sensitivity caused by migration. Additionally, high vulnerability related to low adaptive capacity resulted in all four countries from poor post-harvest infrastructure and limited access to credit and alternative technologies. Often the individual producers did not have access to machines for pulping coffee, drying infrastructure, solar dryers, drip irrigation or water harvesting. The low level of organization was due to the lack of participation of many families in joint activities, projects, training and exchanges. The low level of income diversification was because many families depend mainly on coffee for cash to buy food, healthcare, education and transportation or to invest in the farms (<xref ref-type="table" rid="pone-0088463-t006">Table 6</xref>).</p>
<table-wrap id="pone-0088463-t006" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0088463.t006</object-id><label>Table 6</label><caption>
<title>Vulnerability indicators in relation to adaptation strategies and their specific adaption options.</title>
</caption><alternatives><graphic id="pone-0088463-t006-6" position="float" mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0088463.t006" xlink:type="simple"/>
<table><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup>
<thead>
<tr>
<td align="left" rowspan="1" colspan="1">Vulnerability indicators</td>
<td align="left" rowspan="1" colspan="1">Adaptation strategies</td>
<td align="left" rowspan="1" colspan="1">Specific adaptations options</td>
</tr>
</thead>
<tbody>
<tr>
<td align="left" rowspan="1" colspan="1">Decrease of suitability for coffee production</td>
<td align="left" rowspan="1" colspan="1">Programs of research, validation, transfer and adoption of agricultural technologies that adapt coffee to changing climate</td>
<td align="left" rowspan="1" colspan="1">Drip irrigation, water harvesting and management of available water</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">High variability of annual productivity</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Management of shade, fertility, crop residues, pest and diseases</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Low soil fertility and forest conservation</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Conservation of soil, water and natural forest</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Low income diversification</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Improved varieties and hybrids</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Diversification with other crops where loss of suitability for coffee production</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Poor health and nutrition</td>
<td align="left" rowspan="1" colspan="1">Integral programs with Institutional support improving human and social resources</td>
<td align="left" rowspan="1" colspan="1">Improved environmental education (schools, organizations, committees)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Low level of organizational capacity</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Implementing food and health security programs</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Low level of knowledge of polices of coffee sector and local laws</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Provide cooperatives with social experts to improve the level of participation of producers</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">High migration rate</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Empowering families in policies and laws of their environment sector and to improve implementation</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Low access to credit</td>
<td align="left" rowspan="1" colspan="1">Implementation of long term financial rural programs</td>
<td align="left" rowspan="1" colspan="1">Financial education</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Low viability of post harvest infrastructure</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Planning for the investment of resources</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Low access to technologies</td>
<td align="left" rowspan="1" colspan="1"/>
<td align="left" rowspan="1" colspan="1">Planning of long-term credits (in cash, tools, supplies and others) with technical assistance</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Low access to transport and types of roads</td>
<td align="left" rowspan="1" colspan="1">Implementation of investment programs to improve road infrastructure, quality housing and basic services</td>
<td align="left" rowspan="1" colspan="1">Planning with municipalities, private sector, international cooperation</td>
</tr>
</tbody>
</table>
</alternatives></table-wrap>
<p>Families identified as general strategies for adaptation to climate change the need to develop or improve technologies such as drip irrigation in areas with high risk of drought, shade management, soil fertility management, pest and disease control, conservation of soil and ground water, and adoption of new crops to adapt to future conditions.</p>
</sec></sec><sec id="s4">
<title>Discussion</title>
<p>The strategy used to determine and integrate estimates of exposure, sensitivity and capacity to adapt to climate change proved effective to differentiate between high, medium and low vulnerability families and to identify the livelihood characteristics that contribute to those states. Nevertheless, factors that contribute to vulnerability are distinct between departments, municipalities and families. This is consistent with Cutteŕs analysis that “it is place that forms the fundamental unit of analysis” for vulnerability <xref ref-type="bibr" rid="pone.0088463-Cutter1">[24]</xref>. Thus it is not unexpected that the nature of vulnerability is very site or even family specific <xref ref-type="bibr" rid="pone.0088463-Adger1">[25]</xref>.</p>
<p>The farms located in areas with high vulnerability level will not have suitable conditions for quality coffee production by 2050. These conditions include changing climate factors (e.g. temperature, precipitation) by 2050, high variability of coffee production and high levels of migration in some communities, low adaptive capacity in post harvest infrastructure and in Guatemala and Mexico low access to credit. Some of these are similar to those identified by Eakin <xref ref-type="bibr" rid="pone.0088463-Eakin1">[26]</xref> as the important drivers of change, which include: torrential rainfall, credit, declining soil fertility, new market opportunities, and declining international coffee price. Tucker found that coffee farmers in Central America primarily adapted to global change through changes in crops and crop practices, supporting the importance of diversification options. Nevertheless, these changes were not associated with perceptions of climate risk but rather with market demands <xref ref-type="bibr" rid="pone.0088463-CEPAL1">[10]</xref>. We found that some communities in Mexico, Nicaragua and Guatemala had high migration rates as an additional characteristic of high vulnerability. While migration is definitely on the rise in many parts of Mesoamerica, it appeared that migration in coffee households was difficult to link specifically to climate change drivers <xref ref-type="bibr" rid="pone.0088463-Eakin1">[26]</xref>.</p>
<p>In areas that will remain suitable for coffee growing, but with some reductions in suitability, better agronomic management could lessen the impacts of climate change <xref ref-type="bibr" rid="pone.0088463-Lderach1">[6]</xref>, while in those where a low suitability for coffee growing has been predicted, farmers will have to identify crop alternatives. Solving the problem of variable yield is crucial to the survival of farmers who live in marginalized environments where agro-climatic conditions have always been a challenge. Diversification is, therefore, an important strategy for production risk management in small farming systems. In general, traditional agro-ecosystems are less vulnerable to catastrophic loss because the wide variety of crops and various spatial and temporal arrangements show compensation in case of loss <xref ref-type="bibr" rid="pone.0088463-Altieri1">[27]</xref>.</p>
<p>The adaptation of smallholders also relates to the global supply chain, which is governed by traders and industries that determine the market price and requirements, considering that an upgrading of the coffee chain could help reduce its economic vulnerability <xref ref-type="bibr" rid="pone.0088463-Daz1">[28]</xref>. But ultimately it is the farmers who decide what to farm and how. Strategic decisions must be made with the uncertainty of climate conditions, pricing, costs, government programs, and others factors <xref ref-type="bibr" rid="pone.0088463-Smith1">[29]</xref>.</p>
<p>Adaptation solutions should focus on the development of new infrastructure, policies, and support institutions that facilitate, coordinate, and maximize the benefits of the new systems of land management and use. This can be accomplished by improving governance, ensuring that development programs take climate change into account; increasing investments in irrigation infrastructure and technologies that would increase water use efficiency; creating appropriate transportation and storage infrastructure; reviewing the agrarian property regime and establishing accessible, efficient markets for products, assets and financial services including insurance <xref ref-type="bibr" rid="pone.0088463-IPCC3">[30]</xref>.Finance is critical: access to microloans and formal credit for farm-level investments will help households strategically invest in coffee varieties, complementary crops and livelihood enhancements that effectively reduce risk and improve social welfare <xref ref-type="bibr" rid="pone.0088463-Eakin1">[26]</xref>.</p>
<p>Nevertheless, the solutions require the support of social organizations (civil society groups, cooperatives and small-business organizations) to enable rural households to access the resources and knowledge necessary for adaptation, while empowering communities to shape the direction of the coffee sector to meet their diverse development needs <xref ref-type="bibr" rid="pone.0088463-Eakin1">[26]</xref>. In the case of coping with weather-related hazards, social networks play a primary role in adaptation and recovery. Social and institutional diversity itself promotes resilience <xref ref-type="bibr" rid="pone.0088463-Adger1">[25]</xref>. The decision of a farmer to participate in such an organization is thought to provide improved access to resources and knowledge, as well as to provide a social network that could facilitate recovery and re-establishment of activities following severe shocks (e.g. climatic and price changes). Facilitating adaptation will involve renewed attention to helping households acquire information about markets and new technology. In the face of declining public investment in agriculture, it is clear that public support is needed for research on low-cost, low-input strategies to manage climatic extremes <xref ref-type="bibr" rid="pone.0088463-Eakin1">[26]</xref>. Nevertheless, investment is also needed in the social organization to enable the adaptation of these strategies to local community and individual family vulnerability characteristics.</p>
<p>The role of the state remains important for planned adaptation and sustainable development. Governance is vital in managing global environmental risks and in promoting sustainable technologies. Also, it is necessary for the implementation of human and social programs to have government support and participation of organizations, social networks, education centres and the international cooperation <xref ref-type="bibr" rid="pone.0088463-Adger2">[31]</xref>.</p>
</sec><sec id="s5">
<title>Supporting Information</title>
<supplementary-material id="pone.0088463.s001" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" xlink:href="info:doi/10.1371/journal.pone.0088463.s001" position="float" xlink:type="simple"><label>File S1</label><caption>
<p><bold>Table S1, Key sustainable livelihood indicators of sensitivity and adaptive capacity of coffee producing families in Mesoamerica. Table S2. Current and future values of bioclimatic variables with (averages of coffee growing area in each country) their coefficients of variation (CV) of 19 GCM models for Nicaragua, El Salvador, Guatemala and Mexico.</bold> The CV indicates the variation between models (i.e. the greater the CV the more variability between GCMs).</p>
<p>(DOCX)</p>
</caption></supplementary-material></sec></body>
<back>
<ack>
<p>The authors give special thanks to the coffee growing families and the organizations of each country. Advice and guidance by Dr. Tamara Benjamin, Carolina Backer and Wilfredo Chávez are gratefully acknowledged. We also acknowledge the support of colleagues at CIAT such as Stephania Carmona, Lorena Gómez, Samuel Ocón, Samuel García and Dr. Carlos Zealaya.</p>
</ack>
<ref-list>
<title>References</title>
<ref id="pone.0088463-DFID1"><label>1</label>
<mixed-citation publication-type="other" xlink:type="simple">DFID (Departament for Internacional Development UK) (1999) Hojas Orientativas sobre los Medios de Vida Sostenibles. Available: <ext-link ext-link-type="uri" xlink:href="http://community.eldis.org/.59c21877/SP-GS1.pdf" xlink:type="simple">http://community.eldis.org/.59c21877/SP-GS1.pdf</ext-link> and <ext-link ext-link-type="uri" xlink:href="http://community.eldis.org/.59c21877/SP-GS2.pdf" xlink:type="simple">http://community.eldis.org/.59c21877/SP-GS2.pdf</ext-link>. Accessed 15 October 2009.</mixed-citation>
</ref>
<ref id="pone.0088463-IPCC1"><label>2</label>
<mixed-citation publication-type="other" xlink:type="simple">IPCC (Intergovermental Panel on Climate Change) (2007) Climate Change 2007: Synthesis report. Available: <ext-link ext-link-type="uri" xlink:href="http://www.ipcc.ch/pdf/assessment-report/ar4/syr/ar4_syr.pdf" xlink:type="simple">http://www.ipcc.ch/pdf/assessment-report/ar4/syr/ar4_syr.pdf</ext-link>. Accessed 2009 Oct 15.</mixed-citation>
</ref>
<ref id="pone.0088463-IPCC2"><label>3</label>
<mixed-citation publication-type="other" xlink:type="simple">IPCC (Intergovermental Panel on Climate Change) (2001) Climate Change 2001: Working Group II: Impacts, Adaptation and Vulnerability. Summary for Policy makers.OMM-PNUMA. Ginebra, Suiza. 95p.</mixed-citation>
</ref>
<ref id="pone.0088463-Gay1"><label>4</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Gay</surname><given-names>C</given-names></name>, <name name-style="western"><surname>Estrada</surname><given-names>F</given-names></name>, <name name-style="western"><surname>Conde</surname><given-names>C</given-names></name>, <name name-style="western"><surname>Eakin</surname><given-names>H</given-names></name>, <name name-style="western"><surname>Villers</surname><given-names>L</given-names></name> (<year>2006</year>) <article-title>Potential impacts of climate change on agriculture: A case of study of coffee production in Veracruz, Mexico</article-title>. <source>Climatic Change</source> <volume>79</volume>: <fpage>259</fpage>–<lpage>288</lpage> <comment>DOI:<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1007/s10584-006-9066-x" xlink:type="simple">10.1007/s10584-006-9066-x</ext-link></comment></mixed-citation>
</ref>
<ref id="pone.0088463-Baker1"><label>5</label>
<mixed-citation publication-type="other" xlink:type="simple">Baker P, Haggar J (2007) Global Warming: the impact on global coffee; presentation handout. SCAA.</mixed-citation>
</ref>
<ref id="pone.0088463-Lderach1"><label>6</label>
<mixed-citation publication-type="book" xlink:type="simple">Läderach P, Lundy M, Jarvis A, Ramírez J, Pérez PE, et al. (2010a) Predicted impact of climate change on coffee-supply chains. <italic>In</italic> Leal Filho, W. (ed) The Economic, social and Political Elements of Climate Change, Springer Verlag, Berlin, DE. 19 p.</mixed-citation>
</ref>
<ref id="pone.0088463-Davis1"><label>7</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Davis</surname><given-names>AP</given-names></name>, <name name-style="western"><surname>Gole</surname><given-names>TW</given-names></name>, <name name-style="western"><surname>Baena</surname><given-names>S</given-names></name>, <name name-style="western"><surname>Moat</surname><given-names>J</given-names></name> (<year>2012</year>) <article-title>The Impact of Climate Change on Indigenous Arabica Coffee (<italic>Coffea arabica</italic>): Predicting Future Trends and Identifying Priorities</article-title>. <source>PLoS ONE</source> <volume>7(11)</volume>: <fpage>e47981</fpage> <comment>doi:<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1371/journal.pone.0047981" xlink:type="simple">10.1371/journal.pone.0047981</ext-link></comment></mixed-citation>
</ref>
<ref id="pone.0088463-Jaramillo1"><label>8</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Jaramillo</surname><given-names>J</given-names></name>, <name name-style="western"><surname>Muchugu</surname><given-names>E</given-names></name>, <name name-style="western"><surname>Vega</surname><given-names>FE</given-names></name>, <name name-style="western"><surname>Davis</surname><given-names>A</given-names></name>, <name name-style="western"><surname>Borgemeister</surname><given-names>C</given-names></name>, <etal>et al</etal>. (<year>2011</year>) <article-title>Some Like It Hot: The Influence and Implications of Climate Change on Coffee Berry Borer (<italic>Hypothenemus hampei)</italic> and Coffee Production in East Africa</article-title>. <source>PLoS ONE</source> <volume>6(9)</volume>: <fpage>e24528</fpage>.</mixed-citation>
</ref>
<ref id="pone.0088463-Tucker1"><label>9</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Tucker</surname><given-names>C</given-names></name>, <name name-style="western"><surname>Eakin</surname><given-names>H</given-names></name>, <name name-style="western"><surname>Castellanos</surname><given-names>E</given-names></name> (<year>2010</year>) <article-title>Perceptions of risk and adaptation: Coffee producers, market shocks, and extreme weather in Central America and Mexico</article-title>. <source>Glob Environ Change</source> <volume>20</volume>: <fpage>23</fpage>–<lpage>32</lpage>.</mixed-citation>
</ref>
<ref id="pone.0088463-CEPAL1"><label>10</label>
<mixed-citation publication-type="other" xlink:type="simple">CEPAL (Comisión Económica para América Latina y el Caribe) (2002) Centroamérica: El Impacto de la Caída de los Precios del Café en el 2001. Available: http:// <ext-link ext-link-type="uri" xlink:href="http://www.fondominkachorlavi.org/cafe/docs/cepal2002.pdf" xlink:type="simple">www.fondominkachorlavi.org/cafe/docs/cepal2002.pdf</ext-link>. Accessed 2009 Oct 15.</mixed-citation>
</ref>
<ref id="pone.0088463-Moguel1"><label>11</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Moguel</surname><given-names>P</given-names></name>, <name name-style="western"><surname>Toledo</surname><given-names>V</given-names></name> (<year>1999</year>) <article-title>Biodiversity conservation in traditional coffee systems of Mexico</article-title>. <source>Conservation Biology 13</source> <volume>(1)</volume>: <fpage>11</fpage>–<lpage>21</lpage>.</mixed-citation>
</ref>
<ref id="pone.0088463-Escamilla1"><label>12</label>
<mixed-citation publication-type="other" xlink:type="simple">Escamilla E, Ruiz O, Díaz G, Landeros C, Platas D, et al. (2005) El agroecosistema café orgánico en México. Manejo Integrado de Plagas y Agroecología (Costa Rica) Nro.76p.5–16, 2005.</mixed-citation>
</ref>
<ref id="pone.0088463-Haggar1"><label>13</label>
<mixed-citation publication-type="other" xlink:type="simple">Haggar J, Medina B, Aguilar R, Munoz C (2013) Land use change on coffee farms in southern Guatemala and its Environmental Consequences. Environmental Management. DOI<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1007/s00267-013-0019-7" xlink:type="simple">10.1007/s00267-013-0019-7</ext-link>.</mixed-citation>
</ref>
<ref id="pone.0088463-Yohe1"><label>14</label>
<mixed-citation publication-type="other" xlink:type="simple">Yohe G, Malone E, Brenkert A, Schlesinger M, Meij H, et al. (2006) “A Synthetic Assessment of the Global Distribution of Vulnerability to Climate Change from the IPCC Perspective that Reflects Exposure and Adaptive Capacity.” Palisades, New York: CIESIN (Center for International Earth Science Information Network), Columbia University. Available: <ext-link ext-link-type="uri" xlink:href="http://sedac.ciesin.columbia.edu/mva/ccv/" xlink:type="simple">http://sedac.ciesin.columbia.edu/mva/ccv/</ext-link>.</mixed-citation>
</ref>
<ref id="pone.0088463-Hijmans1"><label>15</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Hijmans</surname><given-names>RJ</given-names></name>, <name name-style="western"><surname>Cameron</surname><given-names>SE</given-names></name>, <name name-style="western"><surname>Parra</surname><given-names>JL</given-names></name>, <name name-style="western"><surname>Jones</surname><given-names>PG</given-names></name>, <name name-style="western"><surname>Jarvis</surname><given-names>A</given-names></name> (<year>2005</year>) <article-title>Very high resolution interpolated climate surfaces for global land areas</article-title>. <source>Int JClimatol</source> <volume>25</volume>: <fpage>19651978</fpage> <comment>DOI:<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1002/joc.1276" xlink:type="simple">10.1002/joc.1276</ext-link></comment></mixed-citation>
</ref>
<ref id="pone.0088463-Ramrez1"><label>16</label>
<mixed-citation publication-type="other" xlink:type="simple">Ramírez J, Jarvis A (2010) Disaggregation of Global Circulation Model Outputs. Available: <ext-link ext-link-type="uri" xlink:href="http:www.ccafs-climate.org/dowlands/docs/Disaggregation-WP-02.pdf" xlink:type="simple">http:www.ccafs-climate.org/dowlands/docs/Disaggregation-WP-02.pdf</ext-link></mixed-citation>
</ref>
<ref id="pone.0088463-Phillips1"><label>17</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Phillips</surname><given-names>SJ</given-names></name>, <name name-style="western"><surname>Anderson</surname><given-names>RP</given-names></name>, <name name-style="western"><surname>Schapire</surname><given-names>RE</given-names></name> (<year>2006</year>) <article-title>Maximum entropy modeling of species geographic distributions</article-title>. <source>Ecol Model</source> <volume>190</volume>: <fpage>231</fpage>–<lpage>259</lpage>.</mixed-citation>
</ref>
<ref id="pone.0088463-Schroth1"><label>18</label>
<mixed-citation publication-type="other" xlink:type="simple">Schroth G, Ruf F (2013) Farmer strategies for tree crop diversification in the humid tropics. A review. Agronomy for Sustainable Development, in press.</mixed-citation>
</ref>
<ref id="pone.0088463-Geilfus1"><label>19</label>
<mixed-citation publication-type="book" xlink:type="simple">Geilfus F (1997) 80 Herramientas para el Desarrollo Participativo: Diagnóstico, planificación, monitoreo y evaluación. IICA, SAGAR Press. México. 210p.</mixed-citation>
</ref>
<ref id="pone.0088463-Surez1"><label>20</label>
<mixed-citation publication-type="book" xlink:type="simple">Suárez Fabio (1999) Fundamentos de Estadística aplicado al sector agropecuario. En: Rojas Eberhard, editores. Santa Fe de Bogotá. 426p.</mixed-citation>
</ref>
<ref id="pone.0088463-DiRienzo1"><label>21</label>
<mixed-citation publication-type="book" xlink:type="simple">Di Rienzo JA, Casanoves F, Balzarini MG, Gonzalez L, Tablada M, et al. (2008) InfoStat version 2008. FCA, Universidad Nacional de Córdoba Press. Argentina.</mixed-citation>
</ref>
<ref id="pone.0088463-Castro1"><label>22</label>
<mixed-citation publication-type="other" xlink:type="simple">Castro F, Montes E, Raine M (2004) Centroamerica la crisis cafetalera: Efectos y estrategias para hacerle frente. The World Bank, Latin America and Caribbean Region, Sustainable Development Working Paper 23.</mixed-citation>
</ref>
<ref id="pone.0088463-Lderach2"><label>23</label>
<mixed-citation publication-type="other" xlink:type="simple">Läderach P, Ovalle O, Lau C, Haggar J, Eitzinger A, et al. (2012b) Climate Change at Mesoamerican Origins. In: Oberthür T, Läderach P, Pohlan J. and Cock J, editors Specialty Coffee: Managing Quality. International Plant Nutrition Institute (IPNI). pp. 39–49.</mixed-citation>
</ref>
<ref id="pone.0088463-Cutter1"><label>24</label>
<mixed-citation publication-type="other" xlink:type="simple">Cutter S (1996) Vulnerability to environmental hazards. Prog Hum Geogr. 20,4 (1996) pp. 529–539.</mixed-citation>
</ref>
<ref id="pone.0088463-Adger1"><label>25</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Adger</surname><given-names>WN</given-names></name>, <name name-style="western"><surname>Kelly</surname><given-names>PM</given-names></name> (<year>1999</year>) <article-title>Social Vulnerability to Climate Change and the Architecture to Entitlements</article-title>. <source> Mitigation and Adaptation Strategies for Global Change</source>. Vol. 47.Nro. <volume>4</volume>: <fpage>253</fpage>–<lpage>266</lpage>.</mixed-citation>
</ref>
<ref id="pone.0088463-Eakin1"><label>26</label>
<mixed-citation publication-type="other" xlink:type="simple">Eakin H, Bojórquez-Tapia LA, Monterde DR, Castellanos E, Haggar J (2011) Adaptive Capacity and Social-Environmental Change:Theoretical and Operational Modeling of Smallholder Coffee Systems Response in Mesoamerican Pacific Rim. Environmental Management (2011) 47:352–367. DOI<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1007/s00267-010-9603-2" xlink:type="simple">10.1007/s00267-010-9603-2</ext-link>.</mixed-citation>
</ref>
<ref id="pone.0088463-Altieri1"><label>27</label>
<mixed-citation publication-type="book" xlink:type="simple">Altieri M, Nicholls C (2009) Cambio Climático y Agricultura Campesina: Impactos y respuestas adaptativas. LEISA. Vol. 24. Universidad de California, Berkeley USA. 4p.</mixed-citation>
</ref>
<ref id="pone.0088463-Daz1"><label>28</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Díaz</surname><given-names>R</given-names></name>, <name name-style="western"><surname>Heakin</surname><given-names>H</given-names></name>, <name name-style="western"><surname>Castellanos</surname><given-names>E</given-names></name>, <name name-style="western"><surname>Jiménez</surname><given-names>G</given-names></name> (<year>2009</year>) <article-title>Condiciones para la adaptación de pequeños productores de café ante presiones económicas mediante procesos de “upgrading” en la cadena productiva. </article-title>. <source>Revista Iberoamericana de la Red de Economía Ecológica. Vol</source>. <volume>10</volume>: <fpage>61</fpage>–<lpage>72</lpage>.</mixed-citation>
</ref>
<ref id="pone.0088463-Smith1"><label>29</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Smith</surname><given-names>B</given-names></name>, <name name-style="western"><surname>Mcnabb</surname><given-names>D</given-names></name>, <name name-style="western"><surname>Smithers</surname><given-names>J</given-names></name> (<year>1996</year>) <article-title>Agricultural adaptation to climatic variation</article-title>. <source>Clim Change </source> 33: 7–29 <volume>1996</volume>: <fpage>23p</fpage>.</mixed-citation>
</ref>
<ref id="pone.0088463-IPCC3"><label>30</label>
<mixed-citation publication-type="book" xlink:type="simple">IPCC (Grupo Intergubernamental de Expertos sobre el Cambio Climático) (2008) El Cambio Climático y el Agua. Documento Técnico VI.OMM-PNUMA. Ginebra, Suiza. 224 p.</mixed-citation>
</ref>
<ref id="pone.0088463-Adger2"><label>31</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Adger</surname><given-names>WN</given-names></name> (<year>2003</year>) <article-title>Social Capital, Collective Action, and Adaptation to Climate Change</article-title>. <source>Econ Geogr</source> <volume>79</volume><issue>(4)</issue>: <fpage>387</fpage>–<lpage>404, 2003</lpage>.</mixed-citation>
</ref>
</ref-list></back>
</article>