- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Volume 2013 (2013), Article ID 359584, 8 pages
Using the Contingent Grouping Method to Value Forest Attributes
1Department of Applied Economics, Autonomous University of Barcelona, 08193 Bellaterra, Spain
2Department of Economics, Rovira i Virgili University, Avinguda Universitat 1, 43204 Reus, Spain
3Department of Economics, Pablo de Olavide University, Ctrretera de Utrera km 1, 41013 Sevilla, Spain
Received 29 November 2012; Accepted 4 February 2013
Academic Editors: P. Newton and M. Pensa
Copyright © 2013 Pere Riera et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This paper presents the first application of a recently proposed stated preference valuation method called contingent grouping. The method is an alternative to other choice modeling methods such as contingent choice or contingent ranking. It was applied to an afforestation program in the northeast of Spain. The attributes included (and the marginal values estimated per individual) were allowing picnicking in the new forests (€2.47), sequestering 1000 tons of CO2 (€0.04), delaying the loss of land productivity by 100 years, due to erosion in the new forests area (€0.783), and allowing four-wheel driving (€6.5), which is perceived as a welfare loss.
Forest ecosystems generate a wide variety of goods and services not only for the forest owners but also for society at large. They provide a number of public goods, like enjoyment from recreational opportunities, nontimber products (e.g., mushrooms, berries, or aromatic herbs), carbon sequestration, erosion prevention and biodiversity preservation, among others. In order to make sound decisions for the whole society, forest planning and management ought to take into account the value of forests for both the landowner and the other affected persons. The field of economics helps in this process by being able to estimate the value, in monetary units, of the forest at stake. Their estimation could constitute a significant source of information for further forest policy design and the development of financial instruments.
Forest valuation is often undertaken from choice modeling techniques. They involve surveying people and asking them to state their preferences among a set of alternatives characterized by attributes fixed at different levels . These preferences may be stated, for example, selecting the most preferred alternative from a choice set (named choice experiment; see, e.g., Louviere et al. ) or ranking the alternatives included in the choice set (named contingent ranking; see, e.g., Chapman and Staelin ) according to their preferences. The different choice modeling variants, like the aforementioned contingent choice and contingent ranking, are able to obtain separate social values for different forest goods and services.
Recently, Brey et al.  proposed a variant named contingent grouping (CG). It requests individuals to classify alternatives included in a choice set as “better than” or “worse than” a status quo or reference situation. The purpose of this paper is to illustrate an application of CG in order to determine how Catalan people choose among potential afforestation programs described by six attributes. The marginal economic values of several forest functions are estimated using a survey with a CG elicitation format undertaken in Catalonia, a region in the northeast of Spain. Following Hensher and Greene  and Hensher et al. , the WTP values are presented in three different ways: a point estimates, the confidence interval, and the whole distribution. The latter can help to better understand the results obtained.
The structure of this paper is as follows. Section 2 introduces the CG method. Section 3 describes the application. Section 4 reports and discusses the results. Finally, Section 5 summarizes the main conclusions.
Within the different valuation methods at hand, stated preference methods are most often used when valuing changes in forests or forest-related goods and services. The contingent valuation method (CVM) tends to be applied when a holistic approach is required, whereas attribute-based valuation methods (ABVM) can isolate the value of different forest attributes.
Typically, CVM analyzes the individual tradeoff between the provision of a good and a payment. The elicitation question can take different formats . For instance, in the closed-ended dichotomous choice variant [8, 9], people may be asked whether they would accept the provision of a nonmarket good at a given price (usually called bid) that varies among subsamples. The allowable answer is therefore closed and dichotomous: yes or no, although a fraction of the population may not know or may not the answer. From the yes or no responses to the different bid amounts, the researcher estimates a response probability function, out of which the mean or median of the individual maximum willingness to pay (WTP) is computed as economic value. Some CVM applications to forests are, for example, [10–16].
An ABVM contains a set, or several sets, of alternatives (choice set) defined by options with different attribute levels, varying across the sample. Usually, the choice sets include a status quo or “business-as-usual” (BAU) option. Individuals express their preferences for the alternatives making some kind of choices, like picking the most preferred-called choice experiment or contingent choice (CC), or ranking the alternatives of the choice set from best to worse, or vice versa, labeled contingent ranking (CR). If one of the attributes is the money that a person would have to pay to secure the change, it is possible to generate estimates of the marginal value of changes in each attribute. In a sense, these methods are able to provide more information than CVM, but at the expense of a more demanding statistical treatment, and maybe a higher burden on the individuals.
CC seems to be relatively easier for respondents to answer than CR, since it resembles the kind of choices that individuals face in actual markets. From a statistical point of view, though, CC provides less information per choice set and individual than CR . The increased complexity of CR may lead some individuals to rank arbitrarily or to engage in strategies of ranking different from utility maximization, which could lead to unreliable welfare measures [17–19]. CC exercises have been applied to forest ecosystems [11, 20–23]. Examples of CR applications in the field of forest economics can be found [24, 25].
A new variant of the ABVM family for modeling preferences for goods, where goods are described in terms of their attributes and levels, is the contingent grouping method (CG) [4, 26]. This variant requests individuals to classify alternatives included in the choice set as “better than” or “worse than” a status quo or reference alternative. This elicitation method provides welfare measures conforming to standard consumer theory. Contingent grouping improves upon CC by generally collecting more information without the cognitive demands of CR . In a way it can be regarded as an appealing alternative to CC and CR, when balancing the amount of statistical information and burden on respondents.
The CG method can be introduced as follows. Consider an individual , , facing a set of five alternatives, as will be used in the empirical application. Denote these alternatives as , , , , and , where represents the status quo or BAU alternative. The utility provided by an alternative can be expressed with a deterministic and a stochastic part, as in where is a vector of observable variables describing the alternative for individual is a vector of unobserved coefficients representing the tastes of individual , and is an unobserved stochastic component independent of and .
Next, consider that the individual is requested to classify the alternatives as “better than” or “worse than” the status quo alternative. Under the assumption of rational behavior, the individual will group the different alternatives according to whether they provide more or less utility than . Therefore, there are 16 possible complete grouping results (see Table 1). In other words, assuming that the probability for an individual to select an alternative from the choice set is given by the probability that its utility is greater than the sum of the utilities of all the other alternatives included in the choice set and that the utilities of these alternatives are independent (i.e., we have independent RUM models ), information about the individual preferences underlying the sixteen different groupings can be obtained as shown in Table 1.
Depending on the BAU position in the choice set in terms of the utility provided to the individual, the information obtained differs. Thus, a BAU dummy variable is to be included in the utility function model to obtain consistent estimates. This variable informs of where option stands in relation to the other alternatives in the preferences of the individual. Using conditional logit models to estimate these probabilities, Brey et al.  showed by numerical simulation that CG generally provides better welfare measure estimates than CC but worse than CR.
The case study application took place in Catalonia, a region in the northeast of Spain. About 40% is occupied by forests. Having a Mediterranean climate in most of the region, Aleppo Pine (Pinus halepensis) is the dominant species; coniferous trees account for 73% of the total forest area. Deciduous trees occupy 14% of the forest land, being Holm oaks (Quercus ilex) the most abundant species. Three quarters of the forest are in private hands, whereas the rest is publicly own, mostly municipal . However, commercial forests are only marginal, with 2% of the agrarian production in Catalonia . The main reason is the low profitability of its timber. Nevertheless, forests provide the Catalan society with many goods and services that seem to increase due to the economic development in the last decades .
Fire is probably the principal disturbance of the Catalan forest. It plays a relevant role in determining the landscape structure and plant community composition, but also in the amount of carbon sequestration and soil erosion.
An afforestation program was proposed in the valuation exercise, increasing the forest surface from 40% to 50% of Catalonia. The expansion would come at the expense of marginal agricultural land. The proposal was in line with the government’s policy. Social and economic changes that occur within developed societies lead up to a situation in which large pieces of agricultural land are being left abandoned as many rural areas become depopulated. Afforestation may be an attractive way of managing fallow lands. Numerous European countries implement subsidized afforestation programs, which provide land owners with financial support for afforestation and management of the planted forest [31, 32].
3.1. CG Design
The program contemplated five different nonmonetary attributes, allowing for picnics in the new forests, four-wheel drive access, picking mushrooms, different amounts of C sequestration, and deferring erosion over time and a required payment. The attributes and levels are summarized in Table 2. The selection of these attributes and levels come from a combination of expert opinions, focus groups, and pilot tests.
A total of 512 combinations were obtained out of the levels and attributes (). A fractional factorial design was then applied to obtain 16 alternatives that were included in four choice sets of four alternatives each.
The first part of the questionnaire described the most relevant aspects of the current forest situation in Catalonia and the afforestation program, with its likely effects. The central part contained the description of the market institution and the value elicitation questions. Each individual was faced with one of the predetermined choice sets and was asked to classify the alternatives as better or worse than the BAU option. The final part of the questionnaire collected several debriefing questions and the request of personal and economic data from the individuals.
The questionnaire went through different focus group sessions and two pilot tests. Focus group participants found the program to be very credible and its description in the questionnaire is understandable for the general public. More than 90% of surveyed people perceived the afforestation program as an overall positive or very positive initiative. Figure 1 reproduces the map included in the questionnaire showing the forested areas before and after the proposed afforestation.
3.3. Data Collection
The questionnaire was administered face-to-face in personal home interviews. A stratified random sample of 800 individuals was selected in Catalonia, from which 732 valid responses were obtained. The strata used were age, gender, and size of the town of residence. The interviews involved 25 randomly selected locations and administered to the population of 18 years of age or older proportionally to the population of each location. In each location, the questionnaires were distributed using random survey routes, with the sample being stratified to include ten respondents selected in terms of gender and age. No significant problems were appreciated in applying the survey. The response rate (91.5%) has been considerably high, suggesting that the cognitive burden of the grouping task may be not too strong.
4. Results and Discussion
The parameters were estimated by maximizing the log-likelihood function , where denotes the –xv expression for the corresponding event associated to individual are the weights associated to the individual, and is the number of sampled individuals. Weights were included in the log-likelihood function to give the same importance to each individual. Expressions –xv were estimated using random parameter logit models [5, 33]. Thus, CG is formulated here as a random parameter logit model with one or repeated choices with stable tastes from the sampled individuals.
The random variables were determined combining two approaches. One was the classical procedure based on estimating different models and using the likelihood ratio tests to select from them. The other procedure was based on the inclusion of artificial variables, as suggested by McFadden and Train . The application of both approaches determined variables and Payment as random. The estimated spreads for these variables were different from zero at 5% significance level. Out of a number of different distributions tested, the triangular showed the best fit. A triangular density function has the shape of a triangle . One of the main advantages of this distribution is that it has relatively short tails, avoiding the strong influence of outliers on the mean estimation . On the other hand, it may provide difficult-to-explain sign changes for some variables. This could be the case, for example, with variable Payment, where the coefficient is always expected to be negative. To avoid this problem, Hensher and Greene  propose making the standard deviation or spread of this distribution a function of its mean, thus ensuring that the distribution provides only positive or negative values. However, it was not necessary to introduce this constraint since the estimated values of the spread for the random parameters were significantly lower than their mean.
All the estimated coefficients, except for Mushrooms, are significant at 10% significance level, and all but Mushrooms and Picnic are significant at 5%. The model contains multiplicative interactions of the status quo dummy variable with some dummies representing different age intervals, since respondents seemed to be more or less prone to select the “business as usual” alternative depending on their age. These dummy variables , , and are set to 1 for individuals between 30 and 44 years of age, 45 and 64, or more than 60-year-old, respectively. The interval age between 18 and 29-year-old is then considered as the reference interval. The results indicate that the probability of choosing the status quo, everything else held constant, increases with the age of the respondents.
4.1. Marginal WTP Estimates
Following Hensher and Greene  and Hensher et al. , the mean WTP, the confidence interval, and its entire distribution were estimated. The simplest approach estimates coefficients as fixed points. Assuming a linear utility model, the WTP is estimated as the negative of the ratio between the mean coefficients of the corresponding attribute and the coefficient of the Payment attribute. Figure 2 shows the WTP estimates.
To estimate the sampling variance, the Krinsky-Robb procedure was adopted . This procedure uses random draws from the estimated asymptotic normal distribution of parameter estimates to calculate different WTP estimates. Figure 2 shows WTP point estimate and the 95% confidence intervals calculated with 1000 replications.
The estimation of the whole marginal WTP distribution is also based on the Krinsky-Robb procedure. For each replication, marginal WTP values were obtained considering the triangular distribution of the random variables according to the following expressions: for attributes Picnic, Drive, Mushrooms, and Erosion , and for attribute , where , and denote random variables following a triangular distribution with centre at zero and spread of one.
The values obtained from the previous expressions can be used to apply some nonparametric density estimators, allowing the description of the distribution of the marginal WTP. This is done without assuming any underlying analytical distribution. Two different estimators of the density function will be reported here: a histogram, which constitutes the simplest representation, and the kernel density estimator, which can be seen as a generalization and improvement over histograms . Figure 2 shows the histograms and the Kernel density estimates constructed from 5000 simulated observations of the marginal WTP, for the five attributes.
The values of 2.47, −1.978, and −6.505 of Picnic, Mushrooms, and Drive, respectively, correspond to the maximum amount of euros (at 2005 values) that on average an individual would pay annually for a discrete change, from not being allowed to picnic, pick mushrooms, or drive cars in the new forests, to being allowed to do those recreational activities. The marginal WTP of 2.735 euros reflects the value that a new forest provides to society by sequestering a quantity equivalent to the emissions of that, on average, a Catalan city of 100.000 inhabitants generates annually in production and consumption activities, that is, approximately equivalent to 68000 tone of per year . Moreover, in the case of the sequestered annually by the new forests, the distribution of the marginal WTP has presented a long tail indicating that a considerable part of the population is willing to pay amounts significantly larger than the estimated mean. Finally, 0.78 euros is the marginal WTP for delaying the loss of land productivity for a hundred years, or in other words, not avoiding erosion for that 100 year period would imply a land productivity cost of 0.783 euros per person and year.
Respondents valued positively the availability of new forests where picnicking is allowed, whereas they seem to be more indifferent toward the possibility of picking mushrooms. The nonsignificance of the Mushrooms coefficient might be due to the fact that picking mushrooms is a leisure activity requiring some specialized knowledge, and thus the number of people involved in this activity being more reduced than in other leisure activities in the forests.
The marginal WTP of Drive indicates that people tend to assess negatively this attribute. Furthermore, the estimated density function never reaches positive values, with a relatively long tail to the left. Thus, the four-wheel driving would not be part of the benefits that individuals perceive from a forest area. Allowing driving may be interpreted by individuals as a source of pollution, confronting with some of the environmental values associated to forests. Four-wheel driving may bring to forests one thing that many respondents might be trying to escape from.
Results also suggest that people tend to assess positively the environmental attributes considered, reinforcing the importance of forest environmental services to society.
This paper has illustrated the use of a recently proposed stated preference variant, the contingent grouping (CG) method, to analyze the way in which individuals regard some attributes representing both recreational and environmental functions of forests in Catalonia, Spain. CG requests individuals to classify alternatives as “better than” or “worse than” a status quo or reference alternative. This elicitation format has only recently been proposed and its feasibility in the empirical ground has to be tested. The results obtained seem to be in line with expectations regarding the sign for the attributes and no particular problems were detected in the CG application, when exposed to respondents. The response rate (91.5%) has been considerably high, suggesting that the cognitive burden of the grouping task may be not too strong.
In this paper, the CG estimations have been formulated as a random parameter logit model with one or repeated choices with stable tastes from the sampled individuals. The WTP estimates have been presented in three different ways: point estimates, confidence intervals, and their entire estimated distribution. On average, in 2005 values, a Catalan over 18-year-old would pay 2.47 euros per year for being allowed to picnic in the new forests and 2.73 euros per year for the sequestered annually by the new forests, which is equivalent to the emissions from a Catalan city of 100.000 inhabitants. The relatively large tail of the estimated distribution seems to indicate that a portion of the population has large WTP for this attribute. Respondents value negatively the possibility that driving through the new forests is allowed; this implying an individual cost of 6.5 euros per year. Likewise, they would be willing to pay 0.783 euros per year to delay one 100-year period of the loss of land productivity caused by erosion.
The results may be relevant for policy design in the afforestation practices in Catalonia. For example, the WTP estimates that could be used by forest planners in a social cost-benefit analysis to assess different variants of an afforestation program are the best to implement. Likewise, policy makes facing mushroom picking legislation ought to be aware of the social effects of such regulations. Management plans for public forests can benefit from the estimates if social welfare is taken into consideration. The results could also be used in forest damage assessments and the estimation of compensations.
CG does not seem to have notably different problems from other choice modeling variants. However, more applications would be needed to further test the applicability of CG.
The authors greatly acknowledge the partial financial support from the Ministry of the Environment of the Spanish and Catalan governments, and from the Forest Technological Center of Catalonia through its MEDFOREX program and the institutional support of the Ministerio de Educación y Cultura of the Spanish government (grant ECO2010-17728).
- N. Hanley, S. Mourato, and R. E. Wright, “Choice modelling approaches: a superior alternative for environmental valuation?” Journal of Economic Surveys, vol. 15, no. 3, pp. 435–462, 2001.
- J. J. Louviere, D. A. Hensher, and J. D. Swait, Stated Choice Method. Analysis and Application, Cambridge University Press, Cambridge, UK, 2000.
- R. G. Chapman and R. Staelin, “Exploiting rank ordered choice set data within the stochastic utility model,” Journal of Marketing Research, vol. 19, no. 3, pp. 288–301, 1982.
- R. Brey, O. Bergland, and P. Riera, “A contingent grouping approach for stated preferences,” Resource and Energy Economics, vol. 33, no. 3, pp. 745–755, 2011.
- D. A. Hensher and W. H. Greene, “The mixed logit model: the state of practice,” Transportation, vol. 30, no. 2, pp. 133–176, 2003.
- D. A. Hensher, J. M. Rose, and W. H. Greene, Applied Choice Analysis. A Primer, Cambridge University Press, New York, NY, USA, 2005.
- R. C. Mitchell and R. T. Carson, Using Surveys to Value Public Goods: The Contingent Valuation Method, Resource for the Future, Washington, DC, USA, 1989.
- R. C. Bishop and T. A. Heberlein, “Measuring values of extramarket goods: are indirect measures biased?” American Journal of Agricultural Economics, vol. 61, no. 5, pp. 926–930, 1979.
- W. M. Hanemann, “Welfare evaluations in contingent valuation experiments with discrete responses,” American Journal of Agricultural Economics, vol. 66, no. 3, pp. 332–341, 1984.
- R. A. Kramer and D. E. Mercer, “Valuing a global environmental good: U.S. residents willingness to pay to protect tropical rain forests,” Land Economics, vol. 73, no. 2, pp. 196–210, 1997.
- E. Lehtonen, J. Kuuluvainen, E. Pouta, M. Rekola, and C. Z. Li, “Non-market benefits of forest conservation in southern Finland,” Environmental Science and Policy, vol. 6, no. 3, pp. 195–204, 2003.
- M. Lockwood, J. Loomis, and T. DeLacy, “A contingent valuation survey and benefit-cost analysis of forest preservation in East Gippsland, Australia,” Journal of Environmental Management, vol. 38, no. 3, pp. 233–243, 1993.
- M. Lockwood, J. Loomis, and T. DeLacy, “The relative unimportance of a nonmarket willingness to pay for timber harvesting,” Ecological Economics, vol. 9, no. 2, pp. 145–152, 1994.
- J. B. Loomis and A. Gonzalez-Caban, “A willingness-to-pay function for protecting acres of spotted owl habitat from fire,” Ecological Economics, vol. 25, no. 3, pp. 315–322, 1998.
- M. Rekola and E. Pouta, “Public preferences for uncertain regeneration cuttings: a contingent valuation experiment involving Finnish private forests,” Forest Policy and Economics, vol. 7, no. 4, pp. 635–649, 2005.
- R. Scarpa, S. M. Chilton, W. G. Hutchinson, and J. Buongiorno, “Valuing the recreational benefits from the creation of nature reserves in Irish forests,” Ecological Economics, vol. 33, no. 2, pp. 237–250, 2000.
- M. Ben-Akiva, T. Morikawa, and F. Shiroishi, “Analysis of the reliability of preference ranking data,” Journal of Business Research, vol. 24, no. 2, pp. 149–164, 1991.
- V. Foster and S. Mourato, “Testing for consistency in contingent ranking experiments,” Journal of Environmental Economics and Management, vol. 44, no. 2, pp. 309–328, 2002.
- J. A. Hausman and P. A. Ruud, “Specifying and testing econometric models for rank-ordered data,” Journal of Econometrics, vol. 34, no. 1-2, pp. 83–104, 1987.
- W. Adamowicz, J. Swait, P. Boxall, J. Louviere, and M. Williams, “Perceptions versus objective measures of environmental quality in combined revealed and stated preference models of environmental valuation,” Journal of Environmental Economics and Management, vol. 32, no. 1, pp. 65–84, 1997.
- P. C. Boxall and B. MacNab, “Exploring the preferences of wildlife recreationists for features of boreal forest management: a choice experiment approach,” Canadian Journal of Forest Research, vol. 30, no. 12, pp. 1931–1941, 2000.
- P. Horne, P. C. Boxall, and W. L. Adamowicz, “Multiple-use management of forest recreation sites: a spatially explicit choice experiment,” Forest Ecology and Management, vol. 207, no. 1-2, pp. 189–199, 2005.
- J. Rolfe, J. Bennett, and J. Louviere, “Choice modelling and its potential application to tropical rainforest preservation,” Ecological Economics, vol. 35, no. 2, pp. 289–302, 2000.
- G. D. Garrod and K. G. Willis, “The non-use benefits of enhancing forest biodiversity: a contingent ranking study,” Ecological Economics, vol. 21, no. 1, pp. 45–61, 1997.
- J. Siikamaki and D. Layton, “Pooled models for contingent valuation and contingent raking data: valuing benefits from biodiversity conservation,” Working Paper, Department of Agricultural and Resource Economics, University of California, Davis, Calif, USA, 2001.
- P. Riera, M. Giergiczny, J. Peñuelas, and P. A. Mahieu, “A choice modelling case study on climate change involving two-way interactions,” Journal of Forest Economics, vol. 18, no. 4, pp. 345–354, 2012.
- R. D. Luce and P. Suppes, “Preference, utility and subjective probabiblity,” in Handbook of Mathematical Psychology, R. D. Luce, R. R. Bush, and E. Galanter, Eds., J. Wiley and Sons, New York, NY, USA, 1965.
- Departament d’Agricultura, RamaderIa I Pesca, DARP, Pla General de Política Forestal. Generalitat de Catalunya, Barcelona, Spain, 1994.
- C. García, Estimació de les macromaginituds agràries de les comarques de Catalunya, 1993, Serveis de Publicacions, Universitat de Lleida, Lleida, Spain, 1997.
- M. Merlo and E. Rojas, “Policy instruments for promoting positive externalities of Mediterranean forests,” in Proceedings of the European Forest Institute, Annual Conference, Chartreuse Ittingen, Switzerland, 1999.
- Council Regulation (EC) No 2080/92 Community aid scheme for forestry measures in agriculture, 30 June 1992.
- Council Regulation (EC) No 1257/1999 Support for rural development from the European Agricultural Guidance and Guarantee Fund (EAGGF) and amending and repealing certain Regulations, 17 May 1999.
- K. Train, Discrete Choice Methods with Simulation, Cambridge University Press, New York, NY, USA, 2003.
- W. H. Greene, NLogit Version 3. 0 Reference Guide, Econometric Software, Plainview, NY, USA, 2002.
- D. McFadden and K. Train, “Mixed MNL models for discrete response,” Journal of Applied Econometrics, vol. 15, no. 5, pp. 447–470, 2000.
- M. Evans, N. Hastings, and B. Peacock, Statistical Distributions, John Wiley and Sons, New York, NY, USA, 3rd edition, 2000.
- I. Krinsky and A. L. Robb, “On approximating the statistical properties of elasticities,” The Review of Economics and Statistics, vol. 68, no. 4, pp. 715–719, 1986.
- B. W. Silverman, Density Estimation for Statistics and Data Analysis, Chapman and Hall, New York, NY, USA, 1986.
- Departament de Medi Ambient, Les emissions a l’atmosfera a Catalunya. Una aproximació quantitativa. Quaderns de Medi Ambient, 5. Generalitat de Catalunya, Barcelona, Spain, 1996.