Introduction

Greenhouse gases and global warming

Solar energy heats the Earth’s surface where some of it is absorbed and the remaining is radiated back toward space as infrared radiation (IR). For 2000 years, the atmospheric concentrations of GHGs have remained relatively constant, but since the 1700s, the concentration of atmospheric carbon dioxide (CO2) has risen about 30%, and the Earth’s average temperature has risen 1.6 ± 0.16 °C (IPCC 2014; Houghton 2005; Keller 2009). Both increases are concomitant and greater than would be expected from natural phenomenon such as volcanism. High atmospheric concentrations of greenhouse gases (GHGs) such as CO2 and methane (CH4), due to the large-scale burning of fossil fuels such as coal and oil for energy, are known to be responsible for global warming. Some other gases such as nitrous oxide (N2O) and refrigerants and dielectric gases are emitted in much smaller quantities than CO2, but are important because they absorb IR radiation 25–22,800 times more efficiently and have atmospheric lifetimes of 125–22,880 years (EPA 2017; IPCC 2014). In the USA, CO2 accounts for 81% of GHG emissions, while CH4, N2O, and refrigerants and dielectric gases account for 10.0%, 5.1%, and 2.8% respectively (EPA 2017).

In 2015, the total US GHG emissions and sinks were calculated to be 5827.7 MMTCO2e, 17.3% of the global emissions of 33,733 MMTCO2e (EPA 2017). Electricity production emitted 1900.7 MMTCO2e, and landfills released 115.7 MMTCO2e, mainly from CH4, produced by anaerobic digestion of organic matter by microbes (EPA 2017). Global warming results from increasing atmospheric GHG concentrations. The ramifications to society are climate change, which is a reorganization of Earth’s weather patterns in response to excess atmospheric energy (Houghton 2005; Held and Soden 2006; Keller 2009; Hansen et al. 2016). Most projections show that climate change will be costly to society, and therefore it is advisable that GHG emissions must be reduced. In the 2015 United Nations global conference on climate change, held in Paris, the USA is committed to reducing its GHG emissions in the range of 17% by 2020, and 26–28% by 2025, relative to 2005 emission levels of 7429 MMTCO2e (IPCC 2014; UN 2016). This commitment was later withdrawn in 2017 (Mohan and MacDonald 2016a b).

Impact of materials

This research has set up mathematical models for computing lifetime GHG emissions associated with materials for service lives of 50 and 100 years, and for three end-of-life treatments of demolition debris: incineration, landfill, and deconstruction. These models are generic and can be applied to any material or assembly of materials; their application in this paper has been demonstrated on lumber from an example residential building.

This paper introduces a new method for measuring lifetime GHG emissions associated with materials termed “Global Warming Impact of Materials (GWIM).” The GWIM is subdivided into two portions: (i) productive portion (GWIMp) that includes the materials’ emissions until the service life of the facility and (ii) non-productive portion (GWIMnp) which includes the materials’ GHG emissions beyond the service life until they are eliminated from the atmosphere. In place of the current, static, EC measurements (kgCO2e or MTCO2e), this model reports the GWIMs in units of kgCO2e-years or MTCO2e-years, which includes the effects of “time of use” of a facility.

The research presented in this paper focuses on the global warming impact of materials, specifically from the GHG emissions occurring during their manufacture, and from GHG emissions arising from disposal of construction and demolition (C&D) debris at the end of their service life. Service life, defined as the time of a structure’s productive use, has been included in the models to calculate the productive and unproductive GHG emissions.

Commercial and residential buildings in the USA use 39.4% of all the energy and 67.9% of all the electricity produced, and are responsible for 38.1% of all the CO2 produced (AIA 2010). The largest emissions occur during the operation phase of a building’s lifetime where more than 50% of the total CO2 emissions occur as a result of heating, cooling, lighting, and powering appliances (DOE 2012). Many agencies are working to improve energy efficiencies for this phase. Generally, the global warming potential (GWP) of building materials is measured by the embodied carbon (EC) of the material. This is the amount of GHG emitted during the manufacture of the material. The EC is a static measure of GWP that does not factor in the length of use of the material or the GHGs emitted after demolition when the C&D debris are disposed of. This work extends EC emissions by including the end-of-life treatments, and models GHG reductions via the use of deconstruction instead of demolition when the building becomes obsolete. This paper has introduced a new concept for measuring the lifetime GHG emissions associated with materials, and termed it Global Warming Impact of Materials (GWIM). For demonstrating the concept, lifetime GHG emissions from building materials have been mathematically modeled for two service lives of 50 years and 100 years, and for three end-of-life treatments at demolition: incineration, landfill disposal, and deconstruction.

Currently, the GHG emissions are measured in kgCO2e or MMTCO2e. This paper has suggested that GHG emissions be measured in kgCO2e-years to include the effect of “time of use,” or service life of a facility. GHG emissions have been computed for two service lives of 50 and 100 years using mathematical models for an example building to demonstrate the effect of service life.

Materials and methods

Modeling lifetime GHG emissions or global warming impact of materials

In this work, the atmospheric concentrations of GHGs associated with building materials, during manufacturing, and end-of-life treatments were simulated with a continuously stirred tank reactor (CSTR), as is done in chemical engineering. Figure 1 shows the schematic of the CSTR model, where the Earth’s atmosphere is represented by the volume of the CSTR. In this model, G is used as a generic variable representing the atmospheric concentration of any GHG. The inlet on the left side of the CSTR represents lifetime GHG emissions into the Earth’s atmosphere. These may come from the manufacture of building materials and from their disposal after the structure is demolished. The GHG emissions associated with manufacturing and from incineration of construction and demolition (C&D) debris are considered instantaneous pulse inputs into the CSTR, as their production times are extremely short. Biogenic carbon is the measure of the amount of GHG emitted when a material is incinerated or used for fuel. In this work, it represents the stored carbon in biological materials such as wood and in organic materials such as plastic. The emission of landfill gas (LFG) from decomposing building materials is a time-dependent source function that gradually releases GHGs into the atmosphere, which depends on the decay rate (k) and methane production potential (Lo) of landfilled materials.

Fig. 1
figure 1

CSTR model of greenhouse gas emissions associated with building materials

The outlet represents the Earth’s physical, chemical, and biological mechanisms that remove CO2 and other GHGs from the atmosphere (Trussell and Hand 2005). The removal mechanisms include partitioning of CO2 into the oceans, CO2 uptake by trees and plants, and oxidation of CH4 by OH-radicals in the atmosphere. The removal mechanisms are a function of the atmospheric residence lifetime of the individual GHGs. Residence lifetime is a concept used in chemical engineering that describes the average amount of time a molecule spends in a reactor before exiting (Trussell and Hand 2005). The atmospheric residence lifetimes of CO2 and CH4 are estimated to be 50–200 years and 12 years respectively (EPA 2017).

The CSTR model is also adaptable to determine the GWP of other GHGs. Refrigerants (e.g., HFCs, CFCs) and dielectric gases (SF6) used in electrical distribution equipment also migrate to the stratosphere and are broken down photochemically by sunlight to produce halogen radicals that catalyze the reduction of O3 to O2.

GHG concentrations in the atmosphere

The derivation of the equation for the atmospheric concentration of a GHG due to an emission input (manufacturing, landfill disposal, incinerations, fugitive emissions, etc.) is shown below. Symbols for the variables and constants used in the GHG emissions model are defined in Table 1. Equation 1 represents the mass flow of a GHG (G) into the atmosphere and a first-order GHG removal mechanism (−QG) that is dependent on the atmospheric concentration (G) and the rate of removal (Q) due to ocean partitioning, plant uptake, etc.

$$ {V}_{\mathrm{atm}}\frac{dG}{dt}=- QG $$
(1)
Table 1 Variables and constants used in modeling GHG emissions

Dividing Eq. 1 by the volume of Earth’s atmosphere (Vatm) gives a differential equation for the concentration of G in terms of the GHG residence time (θ) in the atmosphere:

$$ \frac{dG}{dt}=-\frac{1}{\theta }G $$
(2)

The solution for G(t) from Eq. 2 can be given by:

$$ G(t)={G}_0{e}^{-t/\theta } $$
(3)

with the initial condition G(0) = G0 at t = 0,

where, G0 is the initial concentration of GHG pulse emitted.

Modeling pulse GHG emissions

Equation 3 represents the concentration of the GHG in the atmosphere, emitted as pulse during manufacturing or incineration or as a fugitive emission, and subsequently removed from the atmosphere over time. If the source of the GHG pulse is combustion of the material during incineration, then G0 is given by the mass of the material combusted (m) multiplied by the carbon intensity coefficient, which is termed ECbio in the Inventory of Carbon and Energy (ICE) tables compiled by Hammond et al. (Hammond 2011). Hammond et al. use the subscript “bio” in the ECbio coefficient to refer to the biogenic carbon stored in wood that is released as CO2 upon combustion. In this work, the same coefficient is used to represent the CO2 released from any incinerated, flammable organic material, (e.g., plastic and rubber) and not just wood. Equation 4 gives the concentration of GHG emissions associated with the incineration of a mass (m) of flammable organic material.

$$ {G}_{\mathrm{bio}}(t)=\mathrm{m}\bullet {\mathrm{EC}}_{\mathrm{bio}}{e}^{-t/\theta } $$
(4)

Landfill gas production function

Landfill gas (LFG) production has been previously described by Mohan et al. and EPA (Mohan and MacDonald 2016a b; De La Cruz and Barlaz 2010). LFG production is a complex biological phenomenon carried out by microorganisms. The rate of production and amount of LFG produced by microbes depends on the methane generating potential (Lo) of the material and the microbial decay rate (k) of the material. Methane and CO2 production rates from decomposing building materials were modeled by a modified solution to a first-order landfill decay equation and are given by Eq. 5 (De La Cruz and Barlaz 2010). In Eq. 5, the physical constants are combined together and collectively termed B.

$$ \left[{\mathrm{CH}}_4\right](t)=\left[\frac{M_{\mathrm{o}}{L}_{\mathrm{o}}k\bullet p\bullet \mathrm{mw}}{V_{\mathrm{atm}} RT}\right]{e}^{- kt}=\frac{B}{V_{\mathrm{atm}}}{e}^{- kt} $$
(5)
$$ \mathrm{where}\ B=\left[\frac{M_{\mathrm{o}}{L}_{\mathrm{o}}k\bullet p\bullet \mathrm{mw}}{RT}\right] $$

The factors contained in B are the physical constants that affect the rate of CH4 production in a landfill. The constant M0 is the mass of biodegradable material disposed of in a landfill. The factors mw and R represent the molecular weight of CH4 (16 g/mol) and the ideal gas law constant (0.0821 L atm/mol K) respectively. The factors Lo and k represent the methane production potential of a particular material and the decay rate of biodegradable material. These constants can only be determined empirically, from both laboratory methods and field surveys, and vary widely (Bogner and Spokas 1993; Hansen et al. 2004; Schirmer et al. 2014). This work utilizes the constants Lo = 170 m3/MT and k = 0.05 1/year which are used as parameters for a conventional landfill, in the EPA’s Landfill Gas Emissions Model (LandGEM). Factors p and T represent the barometric pressure of CH4 within the interstices and the temperature within the landfill. For this work, conservative values of p = 1 atm and T = 20 °C (293 K) are used. Pressure and temperature change with landfill depth and season, but range around 1 atm and 20 °C (Yesiller and Hanson 2003).

Multiplying both sides of Eq. 5 by Vatm provides the mass of CH4 (instead of concentration [CH4]) at time t, emitted from landfilled material:

$$ {\mathrm{CH}}_4(t)=\left[\frac{M_{\mathrm{o}}{L}_{\mathrm{o}}k\bullet p\bullet \mathrm{mw}}{RT}\right]{e}^{- kt}=B{e}^{- kt} $$
(6)

Note that the symbols with square brackets represent concentration in mass/volume, while symbols with no brackets represent mass of CH4.

Landfill gas is generally considered to be a mixture of 50% CH4 and 50% CO2. The mass of CO2 produced by landfills is derived as a proportion of the mass of CH4 produced (EPA 2005; De La Cruz and Barlaz 2010). The amount of CO2 in LFG is given as a percentage of CH4 (PCH4) as shown below:

$$ \left(1-{P}_{\mathrm{CH}4}\right)\bullet \mathrm{LFG}(t)=\left(1-{P}_{\mathrm{CH}4}\right)\bullet \left[{\mathrm{CH}}_4(t)+{\mathrm{CO}}_2(t)\right]={\mathrm{CO}}_2(t) $$
(7)

Rearranging Eq. 7 for CO2(t) gives an equation for the mass of CO2 produced as a proportion of the mass of CH4 generated:

$$ {\mathrm{CO}}_2(t)=\left[\frac{1}{P_{\mathrm{CH}4}}-1\right]{\mathrm{CH}}_4(t) $$
(8)

The LFG(t) produced by the landfill is the sum of the masses of CH4(t) and CO2(t) as below. To report the GWP (in CO2 equivalents) of LFG(t) from landfilled materials, CH4 is multiplied by 25 to account for its 25-fold increased heat trapping efficiency of methane.

$$ \mathrm{LFG}(t)={\mathrm{CO}}_2(t)+25{\mathrm{CH}}_4(t) $$
(9)

Substituting Eqs. 7 and 8 into Eq. 9 gives Eq. 10, which provides an expression for GWP (in CO2 equivalents) from LFG production in terms of CH4 and time.

$$ \mathrm{LFG}(t)=\left[\frac{1}{P_{\mathrm{CH}4}}-1\right]{\mathrm{CH}}_4(t)+25{\mathrm{CH}}_4(t)=\left(\frac{1}{P_{\mathrm{CH}4}}+24\right){\mathrm{CH}}_4(t) $$
(10)

Substituting Eq. 6, for the mass of CH4 produced by a mass of landfilled material (CH4(t)), in Eq. 10, provides the expression for the time-dependent production of LFG (Eq. 11), based on the mass of the landfilled material and physical constants.

$$ \mathrm{LFG}(t)=\left(\frac{1}{P_{\mathrm{CH}4}}+24\right)B{e}^{- kt} $$
(11)

Removal of GHG emissions associated with landfilled building materials over time

When considering the time-dependent removal of GHGs from the material disposed of in a landfill, Eq. 1 (for a pulse input of GHG into the atmosphere) must be modified to include the time-dependent input of LFG into the CSTR as well as the time-dependent removal of the GHGs by Earth’s GHG removal mechanisms (sinks). Methane and CO2 emissions are modeled separately and added together at the end, since each gas has a different atmospheric residence time. Equation 12 represents the mass flow of CH4 produced by a mass of material in a landfill (Bekt) flowing into the atmosphere and removal of CH4 by natural mechanisms (−QCH4[CH4]).

$$ {V}_{\mathrm{atm}}\frac{d\left[{\mathrm{CH}}_4\right]}{dt}=B{e}^{- kt}-{Q}_{\mathrm{CH}4}\left[{\mathrm{CH}}_4\right] $$
(12)

Equation 12 gives the production rate of CH4 from anaerobically digested biodegradable organic material. Dividing through by Vatm and rearranging Eq. 12 gives the non-homogeneous differential equation:

$$ \frac{d\left[{\mathrm{CH}}_4\right]}{dt}+\frac{1}{\theta_{\mathrm{CH}4}}\left[{\mathrm{CH}}_4\right]=\frac{B}{V_{\mathrm{atm}}}{e}^{- kt} $$
(13)

where θCH4 is the atmospheric residence time of CH4.

With the initial condition [CH4] = 0, at t = 0, the solution to Eq. 13 is an explicit, time-dependent expression for the atmospheric concentration of CH4.

$$ \left[{\mathrm{CH}}_4\right](t)=\frac{B}{V_{\mathrm{atm}}}\left[\frac{1}{\frac{1}{\theta_{\mathrm{CH}4}}-k}\right]\left[{e}^{- kt}-{e}^{-\frac{t}{\theta_{\mathrm{CH}4}}}\right] $$
(14)

Multiplying through by Vatm gives Eq. 15, which is the mass of atmospheric CH4 at each moment in time.

$$ {\mathrm{CH}}_4(t)=B\left[\frac{1}{\frac{1}{\theta_{\mathrm{CH}4}}-k}\right]\left[{e}^{- kt}-{e}^{-\frac{t}{\theta_{\mathrm{CH}4}}}\right] $$
(15)

To model CO2 production from landfilled materials, Eq. 8 gives the yield of CO2 as a percentage of CH4 produced in a landfill, and Eq. 16 gives the CO2 mass input from a landfill (\( \left[\frac{1}{P_{\mathrm{CH}4}}-1\right]B{e}^{- kt} \)), into the atmosphere with an exit (−QCO2[CO2]) representing Earth’s CO2 removal mechanisms.

$$ {V}_{\mathrm{atm}}\frac{d\left[{\mathrm{CO}}_2\right]}{dt}=\left[\frac{1}{P_{\mathrm{CH}4}}-1\right]B{e}^{- kt}-{Q}_{\mathrm{CO}2}\left[{\mathrm{CO}}_2\right] $$
(16)

Dividing by Vatm and rearranging Equation 16 gives a non-homogeneous differential, where (\( {\theta}_{\mathrm{CO}2}=\frac{V_{\mathrm{atm}}}{Q} \)) is the atmospheric residence time of CO2.

$$ \frac{d\left[{\mathrm{CO}}_2\right]}{dt}+\frac{1}{\theta_{\mathrm{CO}2}}\left[{\mathrm{CO}}_2\right]=\left[\frac{1}{P_{\mathrm{CH}4}}-1\right]\frac{B}{V_{\mathrm{atm}}}{e}^{-\mathrm{k}t} $$
(17)

The solution to Eq. 17, with the initial condition, [CO2] = 0, at t = 0, is an expression for the atmospheric concentration of CO2 from a mass of landfilled material, as a function of time.

$$ \left[{\mathrm{CO}}_2\right](t)=\frac{B}{V_{\mathrm{atm}}}\left[\frac{1}{\frac{1}{\theta_{\mathrm{CO}2}}-k}\right]\left[{e}^{- kt}-{e}^{-\frac{t}{\theta_{\mathrm{CO}2}}}\right] $$
(18)

Multiplying through by Vatm provides the mass of CO2 produced from the landfilled material.

$$ {\mathrm{CO}}_2(t)=B\left[\frac{1}{\frac{1}{\theta_{\mathrm{CO}2}}-k}\right]\left[{e}^{- kt}-{e}^{-\frac{t}{\theta_{\mathrm{CO}2}}}\right] $$
(19)

The total time-dependent GWP of the landfill gas emitted by material (in kgCO2e) is given by Eq. 9, where the mass of CH4 produced is multiplied by 25 to account for methane’s 25-fold greater heat trapping capacity. Substituting Eqs. 15 and 19 into the CH4 and CO2 masses in Eq. 9 gives Eq. 20, which is the GWP due to CH4 and CO2 (LFG) produced in time by a mass of landfilled material.

$$ \mathrm{LFG}(t)=B\left\{\left[\frac{25}{\frac{1}{\theta_{\mathrm{CH}4}}-k}\right]\left[{e}^{- kt}-{e}^{-\frac{t}{\theta_{\mathrm{CH}4}}}\right]+\left[\frac{1}{P_{\mathrm{CH}4}}-1\right]\left[\frac{1}{\frac{1}{\theta_{\mathrm{CO}2}}-k}\right]\left[{e}^{- kt}-{e}^{-\frac{t}{\theta_{\mathrm{CO}2}}}\right]\right\} $$
(20)

Lifetime embodied carbon GHG emissions associated with building materials

The GWP for the lifetime of a building material is the sum of the GHG emitted during manufacture (Eq. 3) and the end-of-life treatment: (i) incineration (Eq. 4) or (ii) landfill disposal (Eq. 20). The time of demolition (TDemo) of a building is the time span from its construction to its demolition. The amount of atmospheric GHG, due to the manufacture, demolition at TDemo, with incineration disposal of the material is given by Eq. 21.

$$ {G}_{\mathrm{I}}(t)=G(t)+{G}_{\mathrm{bio}}(t)={G}_0{e}^{-\frac{t}{\theta_{\mathrm{CO}2}}}+{G}_{\mathrm{bio}}(t)\left\{\begin{array}{c}0\kern9.25em \mathrm{for}\ t<{T}_{\mathrm{Demo}}\\ {}m\bullet {\mathrm{EC}}_{\mathrm{bio}}{e}^{-t/{\theta}_{\mathrm{CO}2}}\kern0.5em \mathrm{for}\ t\ge {T}_{\mathrm{Demo}}\end{array}\right. $$
(21)

Equation 22 represents the amount of atmospheric GHG remaining from the manufacture and landfill disposal of C&D debris at TDemo (GL(t) ).

$$ {G}_L(t)=G(t)+\mathrm{LFG}(t)={G}_0(t){e}^{-\frac{t}{\theta_{\mathrm{CO}2}}}+\mathrm{LFG}(t)\left\{\begin{array}{c}0\kern3.25em \mathrm{for}\ t<{T}_{\mathrm{Demo}}\\ {}\mathrm{LFG}(t)\kern0.5em \mathrm{for}\ t\ge {T}_{\mathrm{Demo}}\end{array}\right. $$
(22)

Embodied carbon GHG emissions associated with entire building, over time

Equation 23 provides the time-dependent (G(t)build), the sum of the G(t) values for manufacture of each building component plus the end-of-life treatments (incineration or landfill:

$$ {G}_{\mathrm{build}}(t)=\sum G(t)+\left\{\begin{array}{c}0\kern3.75em \mathrm{for}\ t<{T}_{\mathrm{Demo}}\\ {}\sum {G}_I(t)\kern1em \mathrm{for}\ t\ge {T}_{\mathrm{Demo}}\\ {}\sum {G}_L(t)\kern1em \mathrm{for}\ t\ge {T}_{\mathrm{Demo}}\end{array}\right. $$
(23)

where, G(t) = GHG emissions from manufacture of all materials in the building.

GI(t) = GHG emissions from combustion of incinerable materials.

GL(t) = GHG emissions from landfill disposal.

Global warming impact due to materials

In this work, a new term, Global Warming Impact of Materials (GWIM), is introduced and defined as below:

“Global Warming Impact of Materials (GWIM) includes all GHG emissions associated with a material including its manufacture (EC), and end-of-life treatment of demolition debris, over a time until they are reduced to zero.”

The GWIM of a material or an assembly of materials is determined from the area under the GHG emissions vs time profiles for materials or assemblies of materials (Eqs. 21, 22, and 23) and is given as Eq. 24.

$$ \mathrm{GWIM}={\int}_0^{T_{\infty }}G(t) dt $$
(24)

In the CSTR model, the atmospheric residence time of molecules is a statistical distribution of the lifetime of molecules in the reactor. The mean of this distribution is the residence time for each GHG (Trussell and Hand 2005). Because this is a statistical distribution, some of the GHG molecules will remain forever in the reactor. The mathematical exponential lifetime GHG emission models show that at 700 years, there are less than 1% of the original GHG molecules left in the atmosphere. Therefore, for purposes of calculating GWIM values, it is assumed that T = 700 years.

Productive and non-productive impact of GHG emissions associated with materials

Any material that is manufactured has an impact on the Earth’s environment. Even the most environmentally friendly materials and manufacturing processes consume raw materials and require energy to produce. In this research, “productive GWIM” (GWIMp) is defined as the GWIM of materials during the time that the material or building is being utilized. Conversely, a non-productive GWIM (GWIMnp) is defined as the GWIM after a building is demolished at TDemo, at the end of its service life. After a building is demolished, some of the embodied carbon from manufacturing of the materials remains in the atmosphere. The remaining emissions can be considered a squandering of the original GHG emissions associated with the materials. These emissions provide no positive value and only negatively affect society by causing global warming. Greenhouse gas emissions from incineration or landfilling of debris are all considered non-productive as they only enhance global warming, but provide no use of the facility. It should be clarified that productive and non-productive do not imply any goodness or badness to the building, instead productive simply implies the material or building is being used and non-productive means it is not. Equations 25 and 26 give the productive and non-productive GWIMs:

$$ {\mathrm{GWIM}}_{\mathrm{p}}={\int}_0^{{\mathrm{T}}_{\mathrm{Demo}}}G(t) dt $$
(25)
$$ {\mathrm{GWIM}}_{\mathrm{np}}={\int}_{T_{\mathrm{Demo}}}^{T_{\infty }}G(t) dt $$
(26)

Added together, GWIMp and GWIMnp give the total GWIM for a material or assembly of materials (Eq. 24). The GWIMp and GWIMnp are of particular importance when describing materials such as cement or glass that have outstanding emissions from their manufacture, but do not emit GHGs after they are disposed of in landfills, against wood or plastics that emit GHGs after they are demolished and either incinerated or landfilled.

Simulations of lifetime GHG emissions for various service lives and various end-of-life treatments for a residential building

The simulations were performed using Microsoft® EXCEL spreadsheet. Integration of the GWIM curves was performed using the trapezoidal rule with the time step set to 0.25 years. In this paper, the lumber portion of the building materials was modeled for demonstrating the simulations. Modeling and computing lifetime GHG emissions, two service lives were selected: 50 years and 100 years, and three most used end-of-life treatments were selected for computations and for making objective comparisons between the treatments: incineration, landfill disposal, and deconstruction.

Simulation results

Lifetime greenhouse gas emissions resulting from manufacturing (EC), service lives, and end-of-life treatments

The structure of the example building used

The example structure, in this paper, is modeled on a residential house in Buffalo, NY. The house is a balloon-framed, wooden structure built about 1900 AD. The plan of the building is a 22 ft × 50 ft (6.7 × 15.2 m) rectangle, 2 stories tall with a 12:12 pitch roof. The house retains most of its original wood siding and moldings. The structure is composed of old-growth Eastern Hemlock (Tsuga canadensis) that was cut to true dimensions, e.g., a 2 × 4 stud (5.08 cm × 10.2 cm) is 2 in. (5.08 cm) thick by 4 in. (10.2 cm) wide unlike modern lumber where a 2 × 4 stud is 1½ (3.8 cm) in. thick by 3½ (8.9 cm) wide. The house has 14 interior doors and 3 exterior doors made of Southern yellow or longleaf pine (Pinus palustris). It has 14, wooden-sash framed exterior windows, made from Eastern white pine (Pinus strobus), jack pine (Pinus banksiana), loblolly pine (Pinus taeda), or similar wood. The foundation is 8 ft (2.4 m) high and composed of large, rough cut blocks of Onondaga limestone. Their lengths and widths vary but their thickness averages 1 ft (0.304 m). The basement floor is 5 in. (12.7 cm) poured cement concrete, thought to have been added in the 1920s. The flooring is ¾ in. (1.09 cm) thick by 2 in. (5.08 cm) wide strips, red oak over a subfloor of 1 in. (2.54 cm) thick by 6 in. (15.2 cm) wide pine (Eastern, jack, loblolly, or other) or hemlock planking. The baseboard, window, and door trim are made from Southern yellow pine. The interior walls of the house are composed of the original sand-cement plaster which averages ½ in. (1.27 cm) thick and is attached to ½ in. (1.27 cm) thick by 2 in. (5.08 cm) wide hemlock lath strips. The exterior siding is original, 6 in. (15.2 cm) wide pine or cedar clapboard. Finally, there are two chimneys made of common, cut, red clay brick, 2 in. (5.08 cm) thick, 4 in. (10.2 cm) wide by 8 in. (20.4 cm) long. The front chimney, built into a fireplace, is irregularly shaped and is more massive than the rear chimney which is 1½ ft (0.457 m) by 1½ ft (0.457 m) square. Both chimneys start in the basement and rise 3 ft (0.91 m) above the ridge line of the roof. In this work, the GWIM values for the lumber component of the house were determined for two service lives of 50 and 100 years, i.e., TDemo denotes the time of demolition. Additionally, the effect of three end-of-life treatments was included: incineration, landfill disposal, and deconstruction. Deconstruction was considered as an alternative to conventional demolition and disposal. The lumber left after deconstruction and recovery was either incinerated or disposed of in a landfill.

Example material—lumber

The example building contains approximately 45,620 boardfeet (BF) (107.7 m3) of softwood lumber and 3300 BF (7.79 m3) of hardwood lumber (total 48,920 BF [115.5 m3]). Pine and hemlock softwoods have a density range of 25–35 lb/ft3 (400 to 550 kg/m3) and an average density of 30 lbs/ft3 (475 kg/m3). Red oak has an average density of about 45 lbs/ft3 (720 kg/m3). Using these densities, the masses of the softwood and hardwood lumber, in the example building are 81,723 kg and 5613 kg respectively. The embodied carbon (EC) and ECbio values used in this research were taken from the Inventory of Carbon and Energy (ICE) v 2.0 tables published by Hammond et al. (Hammond 2011). Using these tables, sawn softwoods have an EC value of 0.2 kgCO2e/kg Hardwood has an EC value of 0.24 kgCO2e/kg and ECbio values of 0.63. The atmospheric lifetime of CO2 is set to 125 years, which is the mean of the range 50–200 years reported by the EPA (EPA 2017). The lifetime GHG emissions have been computed for the 48,920 BF of lumber, in this paper, using mathematical models presented in the earlier sections.

Plotting lifetime GHG emissions associated with lumber

Lifetime GHG emissions from embodied carbon for service lives of 50 and 100 years

The GHG emissions, over time, associated with the manufacture of lumber from the example building are shown in Fig. 2a for TDemo = 50 years and in Fig. 2b for TDemo = 100 years. If demolition was not performed (the “do nothing” scenario) and the material was utilized indefinitely, the GHG emissions remaining in the atmosphere will decrease over time, as shown by the exponential lines in Fig. 2a, b.

Fig. 2
figure 2

a Lifetime GHG emissions associated with building materials for embodied carbon (EC) (TDemo = 50 years). b Lifetime GHG emissions associated with building materials for embodied carbon (EC) (TDemo = 100 years)

The exponential GHG decay curves, line ABE in Fig. 2a and line FGJ in Fig. 2b, are based on decay rates calculated by Eq. 3. For the lumber in the example structure, the total GWIM is calculated to be 21.9 × 105 kgCO2e-years \( \left[\mathrm{GWIM}={\int}_0^{\infty \mathrm{years}}{G}_0{e}^{-t/\left(125\ \mathrm{years}\right)} dt\right] \) where Go is EC coefficients for softwood (0.2 kgCO2e/kg) and for hardwood (0.24 kgCO2e/kg) multiplied by the masses of the softwood (81,723 kg) and hardwood (5613 kg), respectively, in the example structure.

In Fig. 2a, the unshaded region (ABCD) represents the productive portion of GWIM or GWIMp, and the shaded region (area BEC) represents the non-productive portion of GWIM or GWIMnp of 48,920 BF of lumber from the example house, if TDemo = 50 years. The GWIM in this figure includes only GHG emissions from the manufacture (EC) and GHG reduction during the service life of the material. Emissions from end-of-life treatments are not included in Fig. 2a, b. The GWIM values are computed as below.

For TDemo = 50 years:

GWIM = 21.9 × 105 kgCO2e-years

GWIMp = 7.2 × 105 kgCO2e-years, and

GWIMnp = 14.6 × 105 kgCO2e-years.

In Fig. 2b, the unshaded region (area FGHI) represents the GWIMp and the shaded region (area GJH) represents the GWIMnp of 48,920 BF of lumber for TDemo = 100 years. The GWIM values are computed as below:

For TDemo = 100 years:

GWIM = 21.9 × 105 kgCO2e-years

GWIMp = 12.1 × 105 kgCO2e-years, and

GWIMnp = 9.8 × 105 kgCO2e-years.

Extending the service life of the structure from 50 to 100 years, increased the productive GWIM (GWIMP) by 68%, from 7.2 × 105 to 12.1 × 105 kgCO2e-years, and GWIMnp reduced 33% from 14.6 × 105 to 9.8 × 105 kgCO2e-years.

Lifetime GHG emissions from embodied carbon and incineration as end-of life treatment for service lives of 50 and 100 years

Figures 3, 4, and 5 show the effect of the EC, the service lives of the material, and the end-of-life treatment of the demolition debris after the building becomes obsolete and is demolished. The three possible end-of-life treatments, (i) incineration, (ii) landfill disposal, or (iii) deconstruction, have been considered.

Fig. 3
figure 3

a Lifetime GHG emissions associated with building materials for EC and incineration as end-of-life treatment of demolition debris (TDemo = 50) years). b Lifetime GHG emissions associated with building materials for EC and incineration as end-of-life treatment of demolition debris (TDemo = 100 years)

Fig. 4
figure 4

a Lifetime GHG emissions associated with building materials for EC and landfill disposal as end-of-life treatment of demolition debris (TDemo = 50 years). b Lifetime GHG emissions associated with building materials for EC and landfill disposal as end-of-life treatment of demolition debris (TDemo = 100 years) years

Fig. 5
figure 5

a The Global Warming Impact of Materials (GWIM) for 48,920 BF of lumber from the example house. Deconstruction with 30%, 50%, and 70% recovery (TDemo = 50 years). b The Global Warming Impact of Materials (GWIM) for 48,920 BF of lumber from the example house. Deconstruction with 30%, 50%, and 70% recovery (TDemo = 100 years)

In Fig. 3a, the building is demolished at 50 years (TDemo = 50) and all of the construction and demolition (C&D) debris is incinerated. The steep, vertical steps in Fig. 3a, b represent moment of incineration, which is an infinitely short process compared with the atmospheric lifetime of the GHGs. The GWIM values, in Fig. 3a, are calculated from Eq. 21. The total GWIM for this structure, at TDemo = 50 years and incineration as the end-of-life treatment, is 66.2 × 105 kgCO2e-years. The unshaded region (area ABCD), 7.2 × 105 kgCO2e-years represents the productive portion, GWIMp for the material, and the shaded region (area EFC), 58.9 × 105 kgCO2e-years, represents the GWIMnp of the lumber from the example house. If the houses were aged to 100 years at the time of demolition and the lumber were incinerated as the end-of-life treatment, the total GWIM is 66.0 × 105 kgCO2e-years. The unshaded region in Fig. 3b (area GHIJ), 12.1 × 105 kgCO2e-years, represents the GWIMp, and the shaded region (area KLI), 54.0 × 105 kgCO2e-years, represents the GWIMnp. The GWIM values are summarized below:

For TDemo = 50 years:

GWIM = 66.2 × 105 kgCO2e-years

GWIMp = 7.2 × 105 kgCO2e-years, and

GWIMnp = 58.9 × 105 kgCO2e-years.

For TDemo = 100 years:

GWIM = 66.0 × 105 kgCO2e-years

GWIMp = 12.1 × 105 kgCO2e-years, and

GWIMnp = 54.0 × 105 kgCO2e-years.

Lifetime GHG emissions from embodied carbon and landfill as end-of life treatment for service lives of 50 and 100 years

Figure 4a, b represents the lifetime GHG emissions of the lumber from the example house, if the demolition debris is disposed of in landfills, at the end of the building’s service life. Two service lives, TDemo = 50 years and TDemo = 100 years, are used in this demonstration. In Fig. 4a, the GHG emissions associated with manufacture of the lumber and landfill disposal at TDemo = 50 years are shown. The exponential decay curve and GWIM values were calculated using Eq. 22. The total GWIM for this structure with TDemo = 50 years and landfill disposal as the end-of-life treatment is 84.3 × 105 kgCO2e-years. The unshaded region (area ABCD) represents the GWIMp (7.2 × 105 kgCO2e-years) whereas the shaded region (area BEFC) represents the GWIMnp (77.1 × 105 kgCO2e-years). If the house is aged out to 100 years, the GWIM, GWIMp, and GWIMnp values are 84.2. 12.1, and 72.1 × 105 kgCO2e-years respectively. The results are summarizes below.

For TDemo = 50 years:

GWIM = 84.3 × 105 kgCO2e-years

GWIMp = 7.2 × 105 kgCO2e-years, and

GWIMnp = 77.1 × 105 kgCO2e-years.

For TDemo = 100 years:

GWIM = 84.3 × 105 kgCO2e-years

GWIMp = 12.1 × 105 kgCO2e-years, and

GWIMnp = 72.1 × 105 kgCO2e-years.

Deconstruction at the end-of-life demolition

Deconstruction is a process that recovers usable materials at the time of demolition for reuse, recycle, or repurpose. It has been demonstrated that 30% to 70% of the material can be recovered by deconstruction. If the recovered material is reused in new constructions or renovations, the EC of the reused material is used in place of EC for new material, thus lowering the EC of the new structure. The remaining demolition debris, left after deconstruction, is incinerated or landfilled.

Figure 5a presents the GWIM, GWIMp, and GWIMnp for the service life of 50 years, and Fig. 5b for 100 years. Figure 5a gives the various GWIM values for three recovery rates via deconstruction: 30%, 50%, and 70%. The numbers given above the bars in Fig. 5a, b are tabulated in Table 2 for comparing them to the recovery rate of 0% recovery when no deconstruction is done. Table 2 concludes the following:

  1. (i)

    For TDemo = 50 years, GWIM increases from 21.9 × 105 to 66.2 × 105 kgCO2e-years (+ 202%) if the end-of-life debris is incinerated, and the GWIM increases from 21.9 × 105 to 84.3 × 105 kgCO2e-years (+ 285%) if landfill disposal is used.

  2. (ii)

    A 303% increase occurs in the non-productive (GWIMnp) portion of GWIM when incineration is used for demolition debris, and 428% if the demolition debris is disposed of in landfills.

Table 2 Global warming impact of materials in kgCO2e-years (×105) for (i) three end-of-life treatments, (ii) service lives (TDemo) of 50 and 100 years, and (iii) deconstruction with material recovery at 30, 50, and 70%

Deconstruction saves on GWIM values as follows:

  1. (i)

    For 30% recovery, the GWIM is reduced by 20% for incineration and 22.2% if the debris is disposed of in landfills.

  2. (ii)

    Fifty percent material recovery reduces the GWIM by 33.5% for incineration, and 37% for landfill disposal of the demolition debris. These numbers increase to 46.8% for incineration and 51.8% for landfill disposal, for 70% recovery of materials using deconstruction.

  3. (iii)

    For the productive portion of GWIM (GWIMp), 30% material recovery adds 61.1% GWIMp for incineration and landfill disposal. For 50% recovery, the GWIMp for both treatments increases by 103%, and for 70% recovery, both treatments increase the GWIMp by 143%. In the case of non-productive portion of GWIM (GWIMnp), 30% material recovery reduces the GWIMnp by 30% for both incineration and landfill disposal of demolition debris. Fifty percent recovery reduces the GWIMnp by 50% for both end-of-life treatments, and 70% material recovery reduces the GWIMnp by 70%. Similar patterns in GWIM, GWIMp, and GWIMnp follow for TDemo = 100 years, for 30%, 50%, and 70% recovery in dconstruction. These numbers are described in Table 2.

  4. (iv)

    Figure 6a, b displays the reductions in GWIM, GWIMp, and GWIMnp values for 30%, 50%, and 70% material recovery rates after deconstruction. It can be seen that the GWIM values reduce with increasing material recovery rates, in the cases of incineration or landfill disposal of the remaining demolition debris. Also, the non-productive portion (GWIMnp) reduced with increasing recovery percentage. However, the productive portion (GWIMp) increased with increasing recovery rate. All of the above trends follow the same patterns for both service lives of 50 years and 100 years.

Fig. 6
figure 6

a GWIM, GWIMp, and GWIMnp values as functions of recovery percentages after deconstruction. Lumber was incinerated, disposed of in a landfill, or deconstructed with 30, 50, or 70% recovery of materials (TDemo = 50 years). b GWIM, GWIMp, and GWIMnp values as functions of recovery percentages after deconstruction. Lumber was incinerated, disposed of in a landfill or deconstructed with 30, 50, or 70% recovery of materials (TDemo = 100 years)

Summary and conclusions

This research has developed mathematical models for computing lifetime greenhouse gas (GHG) emissions associated with materials. The models include embodied carbon (EC) emissions from the manufacture of materials, and GHG emissions from incineration, or landfill gas (LFG) production from disposal of the material after their service lives. These models are applicable for estimating the global warming impact of any material; however, in this work, they have been used to estimate the global warming impact of lumber from an example residential house. The example structure was modeled with 50- and 100-year service lives to determine the global warming impact of using buildings longer. The global warming impact of the materials, due to their end-of-life treatments (incineration, landfill disposal, or reclamation after deconstruction), was also measured.

This paper introduces a new metric for lifetime GHG emissions associated with materials termed “Global Warming Impact of Materials (GWIM).” The GWIM includes the global warming impact of the manufacturing emissions before the service life and the end-of-life emissions after the service life of the building. The GWIM is subdivided into two portions which represent (i) productive portion (GWIMp) that includes the materials’ emissions before the service life of the facility and (ii) non-productive portion (GWIMnp) which includes the materials’ GHG emissions after the service life until they are reduced to zero. Currently, the global warming impact of materials is given as a static number (kgCO2e or MTCO2e) based on the EC of a material. The models developed in this report give the GWIM in units of kgCO2e-years or MTCO2e-years, which includes the time of use of the facility.

Using the models to calculate the GWIMs for an example residential house, this research has determined that increasing the service life of a building from 50 to 100 years (assuming 100% incineration or landfill disposal of the debris) increases the GWIMp 68.1% and decreases the GWIMnp 8.32% and 6.5% for incineration and landfill respectively.

The models have also been used to compute the effects of deconstructing the example residential house and recovering 30%, 50%, and 70% of the material during demolition. After a 50-year service life, a 30% recovery (assuming the remaining 70% was incinerated or landfilled) affected a GWIM savings of 20% and 22.2% for diversion from incineration and landfill respectively. If at a service life of 50 years deconstruction was performed with a 50% recovery, these values decrease to 33.5% and 37% compared with no deconstruction (0% recovery). If deconstruction is performed with 70% recovery, the GWIM is reduced by 46.8% and 51.8% for incineration and landfill disposal respectively.

In conlcusion, mathematical models can be developed to accurately measure GHG emissions associated with materials at any point in the life of the material, and for any end-of-life treatment: incineration, landfill, or deconstruction. Applied to several alternative materials, these models can help in the decisions to select carbon-efficient materials at the planning stage of the project. These models can also help in determining decisions to rebuild or renovate a project.

The GHG savings due to deconstruction have not been compared and reported in the literature. Models developed in this research have been computed and plotted GHG savings by deconstruction for material recoveries of 30%, 50%, and 70%. The results show the following:

  1. (i)

    Thirty percent material recovery saves approximately 20% GHG savings for both service lives of 50- and 100 years for incineration end-of-life treatment and 22% for landfill disposal.

  2. (ii)

    Fifty percent material recovery saves 34% for both service lives of 50 and 100 years, for incineration end-of-life treatment of the remaining material, and 37% for landfill disposal.

  3. (iii)

    For 70% material recovery, the GHG savings are 47% for both service lives of 50 and 100 years, for incineration and 52% for landfill disposal of the remaining material after material recovery.

The benefits of reclaiming materials, or deconstruction, are not yet well established. The models can compute GHG savings for different material recovery rates. Also, the cost of deconstruction can be estimated and optimal decisions for deconstruction can be made to save on costs and/or on carbon savings.