rheilmayr / chile_subsidies

Replication code for: Forest subsidies reduced native forest extent, carbon sequestration and biodiversity in Chile
HomePage
1 stars 0 forks source link

Request access to carbonModel.py: Summarizes carbon data by region and land use. #1

Open midraed opened 4 years ago

midraed commented 4 years ago

Hello @rheilmayr, is it possible to have access to carbonModel.py? I would like to understand how the data on co2_metrics.csv was generated.

Thank you!

rheilmayr commented 4 years ago

Estimado Guillermo,

Thanks so much for your interest in our research. I'm facing several important deadlines in the coming ~2 weeks so don't think I'll be able to prepare the source files and clean up and properly annotate the carbon calculation script immediately. However, I'm happy to answer questions you might have.

Cristian told me that you may be concerned that we are using the average biomass values as reported here https://www.odepa.gob.cl/publicaciones/articulos/estimacion-del-carbono-capturado-en-las-plantaciones-de-pino-radiata-y-eucaliptos-relacionadas-con-el-dl-701-de-1974-3. That's not quite the approach we took. Instead, we are using the tables reported at the end of the full report (attached here). We start with the distribution of plantation areas under the PFR scenario (e.g. Cuadro A.1.6

Saludos cordiales Robert

Plantations Plantation timber volumes were calculated using data from ODEPA (2010). The report tabulates Chile’s 2005 distribution of plantations by species (i.e., Pinus radiata, Eucalyptus globulus and E. nitens), age, growing region (i.e., macroregions 1, 2 and 3), site quality (i.e., high, medium and low quality classes), management scheme (i.e., pulp production, multipurpose production, intensive management) and ownership type (i.e., large companies or other). In addition, the report presents the expected timber volume for each combination of these categories as generated by the EUCASIM and RADIATA forest growth simulators. Since our observations of land use change span multiple decades, and since the age distribution of plantations is constantly changing, we chose to use the “regulated forest estate” (Patrimonio Forestal Regulado) scenario. This scenario imposes even-aged distributions and constant 9 rotation ages within each species - growing region - site class - management scheme - ownership type combination. In effect, this assumption yields the average carbon density of a plantation over the full span of its rotation. However, this assumption temporarily over-estimates timber volume during periods of plantation expansion. As with native forests, timber volumes were converted to carbon density using Equation 7. Local studies determined the expansion factors for pine (1.567) and eucalyptus (1.77), biomass density for pine (0.3846) and eucalyptus (0.5432) and the carbon share of pine (0.453) and eucalyptus (0.4286) (Gayoso, Guerra, and Alarcón 2002; Gayoso 2002). Aggregate regional carbon densities for each time period were calculated by weighting the carbon densities for pine and eucalyptus by the share of pine as a percent of all plantations(INFOR 1986; 2007; 2008).

On Tue, Jun 23, 2020 at 6:52 PM Guillermo Federico Olmedo < notifications@github.com> wrote:

Hello @rheilmayr https://github.com/rheilmayr, is it possible to have access to carbonModel.py? I would like to understand how the data on co2_metrics.csv was generated.

Thank you!

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub https://github.com/rheilmayr/chile_subsidies/issues/1, or unsubscribe https://github.com/notifications/unsubscribe-auth/ACASICLVBHKON55DETXKZPLRYFL5BANCNFSM4OGGP65Q .

midraed commented 4 years ago

Dear Robert,

Thank you for the quick reply. No worries, I can wait. I just saw that you were willing to share the code, so, I wanted to take a look. I already know the document and tables and the notes of the article supplementary information, but thank you for attaching them here. And please don't invest too much time annotating the code, I use python and have some experience with this kind of data. And for sure I'll rise any question I might have. Best regards,

Guillermo