Cancer mortality and quantitative oil production in the Amazon region of Ecuador, 1990–2010

Purpose Controversy persists over whether cancer risk is increased in communities surrounding oil fields, especially in the Oriente region of Ecuador. This ecologic study uses quantitative exposure data, updated mortality data, and improved statistical methods to study the impact of oil exploration and production activities on cancer mortality rates in the Oriente. Methods Cancer mortality rates in the Oriente in 1990 through 2010 were compared between seven cantons with active oil exploration and production as of 1990 and thirteen cantons with little or no such activities. Poisson regression was used to estimate mortality rate ratios (RRs) adjusted for age and sex. In a two-stage analysis, canton-specific log-RRs were regressed against quantitative estimates of cumulative barrels of oil produced and well-years per canton, adjusting for canton-level demographic and socioeconomic factors. Results Overall and site-specific cancer mortality rates were comparable between oil-producing and non-oil-producing cantons. For overall cancer mortality in males and females combined, the RR comparing oil-producing to non-oil-producing cantons was 0.85 [95 % confidence interval (CI) 0.72–1.00]. For leukemia mortality, the corresponding RR was 0.80 (95 % CI 0.57–1.13). Results also revealed no excess of mortality from acute non-lymphocytic, myeloid, or childhood leukemia. Standardized mortality ratios were consistent with RRs. Canton-specific RRs showed no pattern in relation to oil production volume or well-years. Conclusions Results from this first ecologic study to incorporate quantitative measures of oil exploration and production showed no association between the extent of these activities and cancer mortality, including from cancers associated with benzene exposure. Electronic supplementary material The online version of this article (doi:10.1007/s10552-013-0308-8) contains supplementary material, which is available to authorized users.


Introduction
Little is known about the potential adverse human health impact of oil exploration and production on surrounding communities. In 1989, the International Agency for Research on Cancer (IARC) [1] determined that crude oil is ''not classifiable as to its carcinogenicity in humans,'' based on ''inadequate evidence'' for carcinogenicity in humans and ''limited evidence'' for carcinogenicity in experimental animals. However, questions persist about the health impact of oil exploration and production on surrounding communities. One reason for the paucity of knowledge about the potential environmental health effects of oil production is the difficulty of studying this issue rigorously. Any community health impact of oil production is not readily disentangled from the potential effects of socioeconomic status, sanitation, nutrition, health care access, lifestyle, and other health-related factors that may differ between areas with and without oil fields. Furthermore, many regions with oil fields lack high-quality, population-based data on disease incidence and/or mortality, as well as relevant data on exposure to crude oil or oil-related activities.
To date, the few studies of cancer incidence or mortality in communities with oil exploration and production activities have been ecologic in design and most have been based in the Amazon region of Ecuador, where oil extraction has taken place since 1972. Hurtig and San Sebastián [2] reported excesses in the incidence of overall and several site-specific cancers in four oil-producing cantons, compared with eleven non-oil-producing cantons, in this region in 1985-1998. Incident leukemias, but not other cancers, were also reported to be significantly more common among children in oil-producing cantons [3]. However, in an alternative analysis using cancer mortality data from the same region, Kelsh et al. [4] found no evidence that death from these cancers, or cancer overall, was higher in long-term oil-producing than non-oil-producing cantons. Combined with concerns about data quality and availability, exposure assessment, case ascertainment, population estimation, interpretation of results, and study reproducibility [5,6], the inconsistent cancer incidence and mortality results have failed to resolve the question of whether oil production activities increase the risk of cancer in local populations.
To date, no epidemiologic studies of cancer in communities surrounding oil exploration and production activities have used quantitative information on oil-related activities. Rather, previous studies have broadly classified geographic regions as either active or not active in oil exploration and production, thereby ignoring any variation in the level of activity. To enhance prior findings by capturing the extent of oil-related activities more precisely, we sought to incorporate canton-level data on oil well locations and oil production volumes. In addition, we extended prior studies by using a more flexible and detailed statistical approach, additional years of mortality and population data, and supplemental population data on socioeconomic status, ethnicity, health care access, and residential mobility, to more thoroughly examine cancer mortality in regions with different levels of oil exploration and production activity in the Ecuadorian Amazon region.

Population data
Most oil exploration and production activity in Ecuador is found in the Oriente (East) region within Napo, Pastaza, Orellana, and Sucumbíos Provinces. The population and mortality data of these four provinces from 1990 through 2010 are analyzed in this study.
Population counts for cantons in the Oriente provinces in 1990, 2001, and 2010 were obtained from the Ecuador National Census (www.inec.gov.ec) [7,8]. We also used the 2001 and 2010 census data on residential locations 5 years previously to estimate population counts in 1996 and 2005. To estimate intercensal population counts, we interpolated between the population counts in 1990,1996,2001,2005, and 2010 by using a Poisson regression model that included 5-year age group, sex, year, and age-sex, age-year, and sex-year interactions to account for age-and sex-specific trends in population growth. The expected population, P ij , in the ith age group and jth sex group in each canton was estimated by the following Poisson model: Because Ecuador's administrative divisions have changed in the past 30 years, statistical adjustments to the census and mortality data were made to conform to the administrative divisions in 2010 (Supplementary Table 1).

Mortality data
Annual mortality data from 1990 through 2010 were obtained from the Ecuador National Census. We examined all cancer-related mortality and 25 site-specific cancer causes of death (Supplementary Table 2), including leukemia, childhood leukemia (ages \ 15 years at diagnosis), acute non-lymphocytic leukemia (ANLL), and acute myeloid leukemia (AML, which comprised 83 % of ANLL). Death rates were analyzed based on the canton of residence at the time of death. Foreign residents who died in Ecuador were excluded. Records with missing age (0.25 %) or without a valid code for province/canton (0.06 %) were also excluded.
Oil well and oil production data To quantify the association between mortality and canton-level oil exploration and production activities, we obtained information on oil wells and oil fields from Empresa Pública Petroecuador (www.eppetroecuador.ec/idc/ groups/public/documents/archivo/001373.pdf and www. eppetroecuador.ec/idc/groups/public/documents/archivo/ 001375.pdf). The locations of these wells and fields were overlaid on the province-canton boundaries to quantify oil exploration and production in each canton (Fig. 1). We calculated ''well-year'' as a measure of the cumulative number of oil wells and total years of existence within each canton (Table 1). For a more direct quantification of oil production activity in a given canton, we compiled the total volume of oil produced from 1972 to 2011 based on Petroecuador's annual reports of oil production [9,10] . Oil production volume was reported at the level of the oil field or ''bloque'' (i.e., geographic area designated for oil exploration and production). When an oil field or bloque crossed over canton boundaries, production data were divided proportionally based on the number of wells within each canton. The cumulative number of well-years and the total amount of oil produced from 1972 to 1990 were used to quantify oil exploration and production in each canton. Alternative cutoff dates were evaluated in sensitivity analyses. Seven cantons in the four provinces had a history of oil exploration and production activities as of 1990, whereas thirteen cantons had little or no oil-related activity (Table 1).

Statistical analysis
Comparisons of overall and site-specific cancer mortality between the seven oil-producing cantons and the thirteen non-oil-producing cantons were conducted using both Poisson regression and indirect standardization. Following the age-cohort method proposed by Breslow and Day [11,12], we used Poisson regression to model overall and sitespecific cancer mortality rates as a function of age, sex, and canton-level oil exploration and production activities. The expected number of deaths, D ijk , was calculated from the multiplicative contributions of the ith age group (with ages \ 35 years combined for some analyses), the jth sex group, the activity level of the kth canton, and the age-, sex-, and canton-specific person-years, PY ijk , and was estimated by the following Poisson model: The Oil Activity k factor was equal to 1 if the kth canton was active in oil exploration and production, and 0 otherwise. The parameter associated with this factor provided an estimate of the mortality rate ratio (RR) comparing oilproducing with non-oil-producing cantons. For comparability to prior publications that reported standardized incidence and mortality ratios [2][3][4], we used the indirect standardization method to estimate standardized mortality ratios (SMRs) comparing the observed with the expected number of deaths in the seven oil-producing cantons. The expected number of deaths was calculated using age-and sex-specific mortality rates from the thirteen non-oil-producing cantons and applying those rates to the person-years from the seven oil-producing cantons. For SMR analyses including males and females, the expected number of deaths was calculated as follows: where R ij was the mortality rate for the ith five-year age group and jth sex group in non-oil-producing cantons, and PY ij was the corresponding age-and sex-specific personyears in oil-producing cantons. In SMR analyses of males and females considered separately, the expected number of deaths was summed over age-specific mortality rates and the corresponding age-specific person-years for each sex. We used the method suggested by Rothman and Boice [13] to estimate confidence intervals (CIs) and associated p values for the SMRs.
To further understand the variation in cancer mortality rates among the study cantons, the Poisson model was used to estimate cancer-specific mortality RRs for each of the 20 cantons, without designating particular cantons as active or inactive in oil exploration and production. We used Lago Agrio Canton in Sucumbíos Province as the reference because it had the largest population in the study area; use of a different reference group would not affect the overall results. Scatterplots were created to examine the patterns of association between the RR estimates and oil production metrics, with a nonparametric Loess regression line added to facilitate detection of any trends. To estimate the strength of association more quantitatively, we treated the Poisson regression as the first stage in the regression analysis and, as a second stage, regressed the canton-specific Poisson log-RRs as the dependent variable against cantonlevel oil production volume, well-years, and censusderived data on the proportion of adults who had completed high school, indigenous fraction in the population, availability of health care facilities per capita, and residential mobility in the previous 5 years. Although oil exploration and production began in the 1970s in many areas, we performed sensitivity analyses allowing for an additional 10-year induction period by relating oil production volume or well-years as of 1990 to cancer mortality in 2000-2010.
Statistical analyses were performed with SAS v9.3.

Results
Demographic characteristics and cancer mortality rates of populations residing in the four northern Amazon provinces are summarized in Supplementary Tables 3 and 4. Results from the Poisson regression analysis of cancer mortality in oil-producing versus non-oil-producing cantons among males and females analyzed together and separately are shown in Fig. 2. The corresponding numerical results from both the Poisson regression and SMR analyses are shown in Table 2. For males and females combined, the RR for all cancer-related deaths was 0.85 (95% CI 0.72-1.00) comparing the seven oil-producing cantons with the thirteen non-oil-producing cantons. When males and females were analyzed separately, the RRs showed a similar deficit. We found few consistent elevations in the mortality rate of any site-specific cancer in oil-producing versus non-oil-producing cantons based on either RRs or SMRs in males and females together or separately. Ten or fewer deaths were identified in the oilproducing cantons for each of the following cancers, resulting in imprecise RR estimates: lip/mouth/pharynx (the only cancer for which RR estimates were [1.0 in males, females, and both sexes combined), testis, skin, thyroid, kidney, bladder, and multiple myeloma. Mortality from leukemia was not elevated in oil-producing compared with non-oil-producing cantons (Fig. 2).
Likewise, mortality from ANLL or AML was not higher in cantons that were active in oil exploration and production, although results were based on small numbers of deaths. Leukemia-related mortality among children up to age 14 years also was not associated with the presence of oilrelated activities. Of all the specific cancer sites examined, only mortality from cancer of the lip, mouth, and pharynx was elevated among both males and females in the oilproducing cantons, but estimates were statistically unstable. Classification of canton-level oil production status according to the system used by [2] yielded no substantial difference (data not shown).
The 20 canton-specific RRs from the Poisson regression analysis of overall and site-specific cancer mortality (with Lago Agrio as the reference, RR = 1), adjusted for age and sex, are shown in Fig. 3. The figure reveals no apparent association between oil-related activity in each canton and the RRs for mortality from overall cancer, overall leukemia, childhood leukemia, ANLL, AML, or lymphoma. Rank-ordering of the RRs showed no apparent patterns to suggest increased RRs in the oil-producing cantons. Among the oil-producing cantons, the magnitude of the RR bore no relationship with the amount of oil produced, as represented by the size of the markers in Fig. 3. Due to sparse numbers for AML and ANLL, RRs could not be estimated in some cantons with insufficient data. Scatterplots of the RRs for mortality from overall cancer and other major cancer sites against total volume of oil produced or total well-years also showed no consistent differences in cancer mortality according to level of oil production activity (Fig. 4).
To further examine the association between cancer mortality and the extent of oil exploration and production activities, the age-and sex-adjusted canton-specific log-RRs were regressed against canton-level barrels of oil produced (per 100 million) and well-years (per 1,000) as of 1990, 2000, and 2010, with or without adjustment for educational attainment, indigenous fraction, and health care facilities per capita. Residential mobility was not associated with overall or site-specific cancer mortality and therefore was not included. After multivariate adjustment, no consistent, stable positive associations were observed between the two metrics of oil production and exploration activities and site-specific cancer mortality rates (Table 3). When the analysis was not adjusted for canton-level educational attainment, indigenous fraction, or health care facilities (data not shown), most of the results were not meaningfully different. Well-years but not oil production was positively associated with mortality from cancer of the lip, mouth, and pharynx, whereas inverse associations were detected with mortality from multiple myeloma and cancers of the pancreas, testis, thyroid, and bladder and other urinary organs. Sensitivity analyses allowing for a 10-year induction period revealed no consistent positive associations with overall or site-specific cancer mortality (data not shown).

Discussion
In this ecologic study, we found no evidence of increased overall or site-specific cancer mortality in association with increased level of oil exploration and production activities in the Oriente region of Ecuador. Whether oil-related activity was classified broadly or more finely based on well-years or volume of oil produced, and whether using the traditional SMR approach or the more flexible and detailed Poisson regression approach, we observed no apparent excess of cancer mortality in cantons with more oil exploration and production.
On the contrary, for several cancer sites, mortality was markedly lower in oil-producing than non-oil-producing cantons. If oil-producing cantons have more complete and accurate reporting of cause of death than non-oil-producing cantons due to greater access to mainstream health services as a result of oil-related economic activity, then differential outcome classification would be likely to result in overestimated-not underestimated-RRs. Given that the proportion of death certificates signed by a physician was similar between oil-producing (65 %) and non-oil-producing (58 %) cantons, and census measures of access to health care were also comparable between regions, it is improbable that information bias due to poorer vital statistics reporting in oil-producing than non-oil-producing cantons accounts for the absence of an observed association with cancer mortality. Instead, another plausible explanation for the observed deficits of cause-specific mortality in oil-producing cantons may be unmeasured differences in behavioral, social, cultural, and/or structural factors, rather than a direct beneficial effect of oil-related activities. This explanation is consistent with the fact that most associations were weaker in magnitude after adjustment for canton-level education attainment, indigenous fraction, and density of health care facilities.
Information potentially relevant to the evaluation of the health effects of crude oil exposure can be derived from occupational health studies of oil exploration and production workers. To our knowledge, five cohort studies (each with multiple publications) [14][15][16][17][18][19] and five case-control studies [6,[20][21][22][23] have evaluated cause-specific mortality and/or cancer incidence among oil exploration and  production (i.e., ''upstream'') workers. Overall, no clear picture of excess risk of cancer incidence or mortality has emerged from these studies, and no cancers occurred in significant excess in the majority of studies. Kelsh et al. [4] found that liver cancer mortality in the Ecuadorian Amazon region in 1990-2005 was elevated in cantons with oil production activities. None of the occupational studies described above detected an excess of liver cancer incidence or mortality among upstream petroleum industry workers, nor did we detect such an excess in our updated analysis. Hepatitis B virus (HBV), the cause of the majority of liver cancer worldwide [24], is endemic in the Amazon basin, where 2-14 % of the population is chronically infected, with differences in the prevalence of chronic infection by ethnicity and geographic region [25,26]. Thus, chronic HBV infection may be a major determinant of regional differences in liver cancer incidence and mortality in the Amazon region of Ecuador.  In the study by Hurtig and San Sebastián [2], the two malignancies that accounted for the greatest proportion of the excess overall cancer risk were stomach cancer in men and cervical cancer in women. Neither of these malignancies was consistently positively associated with oilrelated activities in our study, the study by Kelsh et al. [4], or the occupational studies of upstream petroleum industry workers described above. The primary causes of these cancers are also infectious agents, namely Helicobacter pylori as the leading cause of stomach cancer (and a minority of lymphomas), and oncogenic human papillomaviruses (HPV) as the leading cause of cervical cancer (and a substantial proportion of anogenital and oropharyngeal cancers) [24]. Both of these infections are common worldwide, including in Ecuador [27,28]. The prevalence of H. pylori infection varies by individual-and area-level Fig. 4 Scatterplots of relative risk for overall and site-specific cancer mortality versus cumulative well-years (per 1,000) and oil produced (per 100 million barrels) as of 1990, 2000, and 2010, with nonparametric Loess regression lines. AML = acute myeloid leukemia; ANLL = acute nonlymphocytic leukemia; leuk, 0-14 = leukemia in children aged 0-14 years socioeconomic status, urbanization, sanitation, water quality, health care access, ethnicity, and birthplace [29,30], while the prevalence of HPV infection varies by sexual behavior, which in turn depends on population migration and social, cultural, religious, and economic factors [31]. Any of these determinants could explain the observed geographic differences in stomach and cervical cancer incidence in the Ecuadorian Amazon region. Furthermore, disparities in cervical cancer incidence and mortality in Latin American countries have been attributed to differential access to cervical cancer screening and treatment [32].
The other findings of Hurtig and San Sebastián [2,3], including excesses of incident cancers of the rectum, connective/soft tissue, kidney, uterine cervix, and lymph nodes and childhood leukemia, were not confirmed by our study or most studies of oil exploration and production workers. However, Yang and Zhang [33] observed an excess of leukemia around oil fields in China, and Gazdek et al. [34] reported a significant excess of certain hematopoietic malignancies, albeit not lymphomas or all leukemias combined, in Croatian populations living near oil and natural gas fields.
A key methodological difference that may explain part of the inconsistency in results across studies is our use of cancer mortality rather than incidence data. Hurtig and San Sebastián [2,3], Yang and Zhang [33], and Gazdek et al. [34] used cancer incidence data, which more accurately reflect the risk of developing disease than cancer mortality data, especially for cancer types with relatively high survival. However, in regions that lack mandatory populationbased cancer surveillance, incident cases may be missed and those that are reported may represent a biased sample of total incident cases. For example, Hurtig and San Table 3 Associations of age-and sex-adjusted canton-specific log-relative risks of cancer mortality with barrels of oil produced or well-years in 1990, adjusted for canton-level educational attainment, percent indigenous, and health care facilities per capita, northern Amazon provinces of Ecuador, 1990Ecuador, -2010  Sebastián were able to include only incident cancer cases diagnosed in Quito and reported to the National Tumor Registry with a permanent residence in the Amazon region [2,3,35]. Suspected cancer cases in the Amazon region are referred to the capitol city of Quito for diagnosis and treatment, but the long distance-requiring as much as a 12-hour bus ride-and cultural differences between the Amazon and Quito most likely pose a substantial barrier to many residents of the study area. Therefore, cancer incidence among residents of the Amazon region may be grossly underreported, and cancer cases identified in the National Tumor Registry may differ considerably from unreported cases in terms of disease characteristics, patient attributes, and exposures. For example, it is conceivable that cancer cases in oil-producing areas, compared with those in non-oil-producing areas, have better access to navigable roads and/or transportation, enabling them to travel to Quito for diagnosis and treatment and resulting in overestimated relative risks. By contrast, our study used mortality data abstracted from death certificates. Mortality rate ratios are unbiased estimates of incidence rate ratios if the exposure of interest does not affect disease survival or reporting. Currently, no evidence shows that proximity to oil exploration and production activities influences cancer survival. Compared with cancer incidence data in the Ecuadorian Amazon region, mortality data are likely to be more complete and less systematically biased. However, the use of death certificate data, especially in developing regions, entails important limitations in data quality and population coverage. First, the accuracy of the recorded cause of death depends on the diagnostic abilities of the responsible medical facility and/or physician. Second, mortality data may be deficient due to incomplete coverage of the civil registration system, leading to under-registration of deaths by an estimated 13.5 % [36] to 30 % [37], and such underregistration may be unequal between oil-producing and non-oil-producing regions. Third, all death records had a cause of death listed, but 25 % of deaths in oil-producing cantons and 28 % in non-oil-producing cantons had ''symptoms and ill-defined conditions'' identified as the cause of death. Thus, misclassification of causes of death was undoubtedly present, with potential variation across cantons in accordance with degree of development, access to medical care, and other demographic and socioeconomic factors, leading to an unknown impact on the results.
Other differences between our analysis and those of Hurtig and San Sebastián include the definitions of areas with or without oil exploration and production, and methods for estimating the annual population at risk in the study area. We used information on well and oil field locations, drilling dates, and oil field production volumes to characterize the extent of oil exploration and production activities. By contrast, the sources and methods used by Hurtig and San Sebastián were not clearly specified and resulted in different classifications than ours [2,3], although our results were similar when using their classification. We used data from the 1990, 2001, and 2010 national censuses and imputed canton-, age-, and sex-specific population denominators for each intercensal year. Hurtig and San Sebastián used population projections for 1992 or 1993, based on the 1990 national census, as denominators for cancer incidence rates in 1985-1998 or 1985-2000. The latter approach almost certainly underestimated post-1990 populations. Based on national census data showing that oil-producing populations grew more quickly than non-oil-producing populations between 1990 and 2001, this approach would have resulted in overestimated relative risks.
Our study and previous investigations of cancer in communities surrounding oil fields [2-4, 33-35, 38] are all limited by their ecologic design, in which exposure status is assigned at the community level. By assuming that all individuals within a given community have the same exposure status, such studies introduce an unknown degree and direction of bias, as associations observed at the community level may not apply at the individual level. Furthermore, ecologic studies such as these have limited information on potential confounders that may explain observed differences in disease rates between populations. Most studies, including ours, have modest numbers of most site-specific cancers, with correspondingly limited statistical power and ability to control for confounding. An additional concern is that none of these studies can account fully for residential migration and therefore could not classify individuals according to the canton in which they resided for the longest duration or at a biologically relevant latency period prior to death. In fact, biologically plausible latency periods between exposure to oil-related contaminants and cancer diagnosis or death are not established. Finally, we were unable to assess the validity and completeness of oil well and oil production data or cancer mortality data. Even if exposure and outcome misclassification were random, however, the resulting bias would not necessarily attenuate the estimated RRs [39][40][41]. Given these substantial limitations, the reported associations cannot be interpreted as definitively establishing or refuting a causal effect of crude oil on cancer incidence or mortality.
Ideally, studies of the association between residence near oil fields and risk of cancer should use individual-level data on exposure to crude oil and its waste products, as well as abundant data on potential confounders. However, no studies of environmental exposure to oil exploration and production activities have collected such information, and unbiased prospective collection of such data is now virtually impossible, in the aftermath of immense public scrutiny and controversy concerning oil exploration and production in the Ecuadorian Amazon region, along with close involvement of local organizations in setting the agenda of research in the region [42].
Despite these caveats, our study offers several advantages over previously published studies of health outcomes in the Ecuadorian Amazon region. These strengths include more years of follow-up, allowing for longer latency periods and larger sample sizes; more detailed, quantitative information on oil exploration and production; a more refined approach to data analysis; and adjustment for potential confounding by demographic and socioeconomic factors. In particular, a key advantage of the Poisson regression method over the more conventional SMR method is the ability to accommodate a multi-category or continuous rather than binary exposure variable.
In conclusion, in this extended and enhanced analysis of cancer mortality in the Oriente region of Ecuador from 1990 through 2010, we observed no apparent excess of death from any or all cancers in areas with oil exploration and production activities, compared with areas that had little or no oil-related activity. Given the methodological limitations of this study, our findings do not necessarily indicate that exposure to crude oil and oil-related activities is causally unrelated to any form of cancer. However, our findings provide no evidence to support such a relationship and further demonstrate that in the Ecuadorian Amazon region, residing near oil fields appears not to adversely affect cancer mortality.