Simulation Studies to Quantify the Impact of Demand Side Management on Environmental Footprint

: The increased use of energy leads to increased energy-related emissions. Demand side management (DSM) is a potential means of mitigating these emissions from electric utility generating units. DSM can signiﬁcantly reduce emissions and provide economic and reliability beneﬁts. This work presents some DSM techniques, such as load shifting, energy conservation, and valley ﬁlling. Furthermore, this work explains the most common DSM programs. To quantify the effect of DSM in diminishing carbon footprint, this paper performs power ﬂow analysis on a yearly load proﬁle corresponding to Fort Collins, Colorado, U.S. This work used the IEEE 13-node test system to simulate several scenarios from the multi-criteria decision-making (MCDM) alternatives, both individually and integrated. For the base case, emissions decrease by 16% from the 2005 level. The “energy conservation” option achieved a 20% reduction in emissions, integrating both alternatives increased the emissions mitigation up to 22%. Simulation of the residential sector shows the “communication and intelligence” option reduces emissions about 14% from the 2005 level. A scenario that combines “electric stationary storage” with “communication and intelligence” diminishes the emissions by more than 15%. The last scenario examined all MCDM alternatives combined into one option, resulting in a 20% emissions reduction. We also conducted a cost beneﬁt analysis (CBA) to investigate economic, technical, and environmental costs and beneﬁts associated with each alternative. The economic evaluation shows that “electric stationary storage” is the best option since it charges during lower electricity prices and discharges during peaking demand. The economic analysis presents a trade-off chart, so the decision maker can select the alternative based on their preference.


Introduction
The increased use of energy leads to increased energy-related emissions. It is important to control the greenhouse gas (GHG) portion of these emissions from electric generating units. The previous work presented in [1] discussed several governmental plans to mitigate the impact of GHG. Reference [1] also investigated using MCDM algorithms to pare down the numerous options for demand side management (DSM) by establishing a problem framework corresponding to the Climate Action Plan (CAP) goals of Fort Collins, Colorado, U.S. The findings showed that "communication and intelligence" has the highest impact in meeting the goal followed by "stationary electric storage" and "energy conservation".
This manuscript aims to prioritize the MCDM alternatives based on their environmental impact. This work evaluates the environmental and economic aspects of the MCDM alternatives to find an optimal combination among the DSM options. In that regard, the paper quantifies the impact of DSM in reducing the environmental footprints using the IEEE 13-node test feeder through the Open Distribution System Simulator ® (OpenDSS ® ) This test feeder includes data for lines, transformers, capacitors, spot loads, and distributed loads. The next section explains the findings after simulating the base case on a distribution system emulating the one in Fort Collins. We also present some strategies to test the ability of the MCDM solutions in achieving emissions reductions.

Modeling and Simulation of the Distribution Systems
The need to shape the load profile by controlling demand became apparent in the 1970s [13]. Electric load profiles change daily, seasonally, and annually. Therefore, this simulation considers the yearly load profile from which we can obtain the output of the generators to meet the demand and, in turn, estimate the GHG emissions. To determine the environmental impact, the load profile corresponding to Fort Collins is translated and mapped on the IEEE 13-node system through the OpenDSS ® simulation tool. The IEEE 13node test feeder is designed to evaluate and benchmark algorithms, and this test system provides simple ways to modify the test feeders to include DERs. Thus, it is used to adapt the load of the Fort Collins distribution grid on the small circuit test system. The load curve in 2017 has a peak of 660 MW. This demand curve is scaled down to match the peak of the IEEE 13-node test system, i.e., 3.577 MW. To establish a strategy for GHG reduction, this section simulates the load profile and generation mix for the base case. After that, the analysis considers different scenarios from the MCDM options to reduce the GHG emissions compared to the base case. The analysis also combines several scenarios from among the MCDM alternatives to simulate their emissions reduction and to investigate the effectiveness of each alternative in achieving the goal. In that regard, the work studies the most preferred alternatives from the MCDM ranking list from [1], i.e., communication and intelligence, stationary storage, and energy conservation.

Base Case
To obtain the load profile and generation mix for the year 2017, using OpenDSS ® with the COM interface, the actual load curve is mapped on the test system. Figure 2 illustrates the demand curve and generation mix for the base case on the IEEE 13-node test feeder using Fort Collins load data. This figure indicates that the annual peak demand occurred on July 19th, 2017, at 3 pm. After analyzing the supply and demand curves, environmental

Modeling and Simulation of the Distribution Systems
The need to shape the load profile by controlling demand became apparent in the 1970s [13]. Electric load profiles change daily, seasonally, and annually. Therefore, this simulation considers the yearly load profile from which we can obtain the output of the generators to meet the demand and, in turn, estimate the GHG emissions. To determine the environmental impact, the load profile corresponding to Fort Collins is translated and mapped on the IEEE 13-node system through the OpenDSS ® simulation tool. The IEEE 13-node test feeder is designed to evaluate and benchmark algorithms, and this test system provides simple ways to modify the test feeders to include DERs. Thus, it is used to adapt the load of the Fort Collins distribution grid on the small circuit test system. The load curve in 2017 has a peak of 660 MW. This demand curve is scaled down to match the peak of the IEEE 13-node test system, i.e., 3.577 MW. To establish a strategy for GHG reduction, this section simulates the load profile and generation mix for the base case. After that, the analysis considers different scenarios from the MCDM options to reduce the GHG emissions compared to the base case. The analysis also combines several scenarios from among the MCDM alternatives to simulate their emissions reduction and to investigate the effectiveness of each alternative in achieving the goal. In that regard, the work studies the most preferred alternatives from the MCDM ranking list from [1], i.e., communication and intelligence, stationary storage, and energy conservation.

Base Case
To obtain the load profile and generation mix for the year 2017, using OpenDSS ® with the COM interface, the actual load curve is mapped on the test system. Figure 2 illustrates the demand curve and generation mix for the base case on the IEEE 13-node test feeder using Fort Collins load data. This figure indicates that the annual peak demand occurred on 19 July 2017, at 3 pm. After analyzing the supply and demand curves, environmental assessment calculates the emissions generated from the conventional generating units for meeting the demand. Four generators use fossil fuels (and consequently emit CO 2 ) in the system: (a) coal 1, (b) coal 2, (c) coal 3, and (d) natural gas [14]. 13,692 tons of CO2 per year. A GHG equivalence calculator equates this to the emissions produced by burning more than 7484.29 tons of coal per year and the captured emissions from about 16,115 acres of the U.S. forests in one year [15]. As the city's energy environmental indicator shows a 16% emissions reduction in 2017 from 2005 level due to the utilization of more renewable generation, the simulation follows with the same reduction in emissions from 2005 level of 16.26% for the same year [16].

ESS
Storage is one of the top priorities of MCDM technologies [1]. ESS can be used during peak hours to shave the load using the energy stored in ESS during off-peak periods [17]. The base case results in Figure 3 show that the emissions from the electricity sector is 13,692 tons of CO 2 per year. A GHG equivalence calculator equates this to the emissions produced by burning more than 7484.29 tons of coal per year and the captured emissions from about 16,115 acres of the U.S. forests in one year [15]. As the city's energy environmental indicator shows a 16% emissions reduction in 2017 from 2005 level due to the utilization of more renewable generation, the simulation follows with the same reduction in emissions from 2005 level of 16.26% for the same year [16].  Figure 3 show that the emissions from the electricity sector is 13,692 tons of CO2 per year. A GHG equivalence calculator equates this to the emissions produced by burning more than 7484.29 tons of coal per year and the captured emissions from about 16,115 acres of the U.S. forests in one year [15]. As the city's energy environmental indicator shows a 16% emissions reduction in 2017 from 2005 level due to the utilization of more renewable generation, the simulation follows with the same reduction in emissions from 2005 level of 16.26% for the same year [16].

ESS
Storage is one of the top priorities of MCDM technologies [1]. ESS can be used during peak hours to shave the load using the energy stored in ESS during off-peak periods [17].

ESS
Storage is one of the top priorities of MCDM technologies [1]. ESS can be used during peak hours to shave the load using the energy stored in ESS during off-peak periods [17]. This, in turn, minimizes the need for powering on generators that emit high volumes of CO 2 during times of high demand. Platte River Power Authority (PRPA), an electric utility that generates and delivers electricity in Fort Collins, proposes using a Li-ion battery since it is the second most used technology (after Pb-acid) in the stationary storage market. Furthermore, the Li-ion technology can operate over more and deeper cycles than a lead acid battery, resulting in a lower cost per cycle. PRPA proposes batteries that can supply a peak power or 50 MW for four hours, totaling 200 MWh of energy storage.
This simulation scales down the storage parameters to fit the test system. Therefore, the ESS provides about 271 kW for the duration of four hours (i.e., 1084 kWh) to shave the peak load. After looking at the 2017 load curve, it is noteworthy that the peak hourly load usually exceeds the average hourly load by 20%. Thus, the simulation is designed to enable the ESS to shave 10% of the peak load when the load at the specified hour exceeds the average load by 20%. Coincident peak refers to the highest user demand that occurs one hour a month when the system demand is at its highest [18]. During the coincident peak, a storage system can provide aggressive reductions in peaking load, say, up to 15%. The storage system remains inactive when the load is lower than 120% of the average load. The system charges the energy storage during off-peak hours when the load is at its lowest. Figure 4 shows the load and generation mix from the test system for 2017, while Figure 5 represents a simulation of the load profile on 5 September 2019. The results show that ESS, operated according to the charge-discharge technique described above, mitigates the total emissions from the 2005 level by 18.13%, equivalent to 13,385 tons of CO 2 per year. Figure 6 explains the emissions reduction per source over a year after integrating ESS. This mitigation comes from diminishing the emissions from the coal 1 power plant. This, in turn, minimizes the need for powering on generators that emit high volumes of CO2 during times of high demand. Platte River Power Authority (PRPA), an electric utility that generates and delivers electricity in Fort Collins, proposes using a Li-ion battery since it is the second most used technology (after Pb-acid) in the stationary storage market. Furthermore, the Li-ion technology can operate over more and deeper cycles than a lead acid battery, resulting in a lower cost per cycle. PRPA proposes batteries that can supply a peak power or 50 MW for four hours, totaling 200 MWh of energy storage. This simulation scales down the storage parameters to fit the test system. Therefore, the ESS provides about 271 kW for the duration of four hours (i.e., 1084 kWh) to shave the peak load. After looking at the 2017 load curve, it is noteworthy that the peak hourly load usually exceeds the average hourly load by 20%. Thus, the simulation is designed to enable the ESS to shave 10% of the peak load when the load at the specified hour exceeds the average load by 20%. Coincident peak refers to the highest user demand that occurs one hour a month when the system demand is at its highest [18]. During the coincident peak, a storage system can provide aggressive reductions in peaking load, say, up to 15%. The storage system remains inactive when the load is lower than 120% of the average load. The system charges the energy storage during off-peak hours when the load is at its lowest.

Energy Conservation
This alternative means reducing the demand and, thus, the energy needed to supply it. Energy conservation can be implemented either by deploying more energy-efficient loads or educating users to change their usage behavior to reduce the work requested of the system. In that regard, this scenario uses the ideal case for the simulation based on a 5% reduction of the entire load. Therefore, a one-year simulation is implemented after reducing the demand by 5% at all time steps. The red curve in Figure 7 illustrates such a reduction in the total demand curve for 2017. This energy saving successfully minimized the total emissions by more than 20% compared to the 2005 level. The amount of emissions generated from the test system after reducing the energy demand is equal to 13,016 tons of CO 2 per year, compared to 16,350 tons of CO 2 in 2005. Figure

Energy Conservation
This alternative means reducing the demand and, thus, the energy needed to supply it. Energy conservation can be implemented either by deploying more energy-efficient loads or educating users to change their usage behavior to reduce the work requested of the system. In that regard, this scenario uses the ideal case for the simulation based on a 5% reduction of the entire load. Therefore, a one-year simulation is implemented after reducing the demand by 5% at all time steps. The red curve in Figure 7 illustrates such a

ESS with Energy Conservation
This scenario studies the load behavior and emissions amount after combining the two options: ESS and energy conservation. The aim of considering such a scenario is to investigate the potential of achieving greater reductions in emissions than deploying just

ESS with Energy Conservation
This scenario studies the load behavior and emissions amount after combining the two options: ESS and energy conservation. The aim of considering such a scenario is to investigate the potential of achieving greater reductions in emissions than deploying just

ESS with Energy Conservation
This scenario studies the load behavior and emissions amount after combining the two options: ESS and energy conservation. The aim of considering such a scenario is to investigate the potential of achieving greater reductions in emissions than deploying just one option. This scenario follows the same simulation strategies for ESS and conservation programs that are explained in Sections 3.2 and 3.3. Figure 9 illustrates that such a combination affects the total energy demand. This scenario can achieve a 22% reduction in emissions, which is more than the city's goal of reaching a 20% emissions reduction by Sustainability 2021, 13, 9504 8 of 24 2020. As stated, the base case generates 16,350 CO 2 tons of emissions while this option can minimize it to 12,725 CO 2 tons. Figure 10 presents the generated emissions from each emitting source. Unsurprisingly, the ESS displaced the energy needed from Coal 1. one option. This scenario follows the same simulation strategies for ESS and conservation programs that are explained in Sections 3.2 and 3.3. Figure 9 illustrates that such a combination affects the total energy demand. This scenario can achieve a 22% reduction in emissions, which is more than the city's goal of reaching a 20% emissions reduction by 2020. As stated, the base case generates 16,350 CO2 tons of emissions while this option can minimize it to 12,725 CO2 tons. Figure 10 presents the generated emissions from each emitting source. Unsurprisingly, the ESS displaced the energy needed from Coal 1.   one option. This scenario follows the same simulation strategies for ESS and conservation programs that are explained in Sections 3.2 and 3.3. Figure 9 illustrates that such a combination affects the total energy demand. This scenario can achieve a 22% reduction in emissions, which is more than the city's goal of reaching a 20% emissions reduction by 2020. As stated, the base case generates 16,350 CO2 tons of emissions while this option can minimize it to 12,725 CO2 tons. Figure 10 presents the generated emissions from each emitting source. Unsurprisingly, the ESS displaced the energy needed from Coal 1.

Communication and Intelligence
The aim of using communication and intelligence is to improve the matching of load to the availability of variable clean generating sources and to reduce peak demand. This option is a means of implementing residential DR encompassing the Internet of Things, smart meters and AMI, real-time pricing, and smart appliances with controls Sustainability 2021, 13, 9504 9 of 24 for automating energy conservation. Therefore, the simulation only applies DR to the residential load curve.
Fort Collins Utility's electric rates in the summer season (May-September) are USD 0.07 during off-peak and USD 0.24 from 2 p.m. to 7 p.m. The rates in non-summer seasons (October-April) remain the same, USD 0.07, during off-peak hours, but the prices during peak demand hours (5 p.m. to 9 p.m.) are USD 0.22 [19,20].
In DR, as with most other commodities, the demand changes as the prices change according to the price elasticity of demand, ε.
where ∆Q D is the change in demand and ∆P is the change in prices [21]. In [22], there are two values for price elasticity of demand: long-term and short-term. Since this work conducts a one-year simulation, we consider the price elasticity for the short-term demand, −0.02. After performing the analysis, we notice that the Time-of-Day rates help in reducing the load during peak hours. Figure 11 illustrates the residential load profile after applying DR in the residential sector. Figures 12 and 13 demonstrate the change in the demand with DR during the summer season and the non-summer season, respectively.
While the goal of deploying DR is shaving the peak load, the total emissions from the residential sector is reduced by about 14% from the residential emissions in 2005 level. The pollution from powering the residential sector was 5068 CO 2 tons in 2005 levels and then DR dropped the emissions to 4363 CO 2 tons. Figure 14 explains the emissions curve for each source.

Communication and Intelligence
The aim of using communication and intelligence is to improve the matching of load to the availability of variable clean generating sources and to reduce peak demand. This option is a means of implementing residential DR encompassing the Internet of Things, smart meters and AMI, real-time pricing, and smart appliances with controls for automating energy conservation. Therefore, the simulation only applies DR to the residential load curve.
Fort Collins Utility's electric rates in the summer season (May-September) are USD 0.07 during off-peak and USD 0.24 from 2 p.m. to 7 p.m. The rates in non-summer seasons (October-April) remain the same, USD 0.07, during off-peak hours, but the prices during peak demand hours (5 p.m. to 9 p.m.) are USD 0.22 [19,20].
In DR, as with most other commodities, the demand changes as the prices change according to the price elasticity of demand, .
where ∆ is the change in demand and ∆ is the change in prices [21]. In [22], there are two values for price elasticity of demand: long-term and short-term. Since this work conducts a one-year simulation, we consider the price elasticity for the short-term demand, −0.02.
After performing the analysis, we notice that the Time-of-Day rates help in reducing the load during peak hours. Figure 11 illustrates the residential load profile after applying DR in the residential sector. Figures 12 and 13 demonstrate the change in the demand with DR during the summer season and the non-summer season, respectively.
While the goal of deploying DR is shaving the peak load, the total emissions from the residential sector is reduced by about 14% from the residential emissions in 2005 level. The pollution from powering the residential sector was 5068 CO2 tons in 2005 levels and then DR dropped the emissions to 4363 CO2 tons. Figure 14 explains the emissions curve for each source.      Emissions from the electricity residential sector from the IEEE test system for 2017, after using residential DR.

Communication and Intelligence with ESS
This subsection investigates the combination of residential DR and ESS in diminishing the emissions from powering the residential electric sector. In such a simulation, we apply the same above-described approaches for residential DR and ESS. Since both these options aim to shave the peaking load, Figure 15 demonstrates the reduction in the peak load. Integrating ESS helped in saving emissions from the residential sector by about 2.48%. As a result, this scenario helped in minimizing the total emissions to almost 15.68% compared to the 2005 level. The total emissions obtained from the IEEE 13-node system dropped to 4273 CO2 tons for this scenario. Figure 16 illustrates the generated emissions after applying this scenario. Figure 14. Emissions from the electricity residential sector from the IEEE test system for 2017, after using residential DR.

Communication and Intelligence with ESS
This subsection investigates the combination of residential DR and ESS in diminishing the emissions from powering the residential electric sector. In such a simulation, we apply the same above-described approaches for residential DR and ESS. Since both these options aim to shave the peaking load, Figure 15 demonstrates the reduction in the peak load. Integrating ESS helped in saving emissions from the residential sector by about 2.48%. As a result, this scenario helped in minimizing the total emissions to almost 15.68% compared to the 2005 level. The total emissions obtained from the IEEE 13-node system dropped to 4273 CO 2 tons for this scenario. Figure 16 illustrates the generated emissions after applying this scenario. Emissions from the electricity residential sector from the IEEE test system for 2017, after using residential DR.

Communication and Intelligence with ESS
This subsection investigates the combination of residential DR and ESS in diminishing the emissions from powering the residential electric sector. In such a simulation, we apply the same above-described approaches for residential DR and ESS. Since both these options aim to shave the peaking load, Figure 15 demonstrates the reduction in the peak load. Integrating ESS helped in saving emissions from the residential sector by about 2.48%. As a result, this scenario helped in minimizing the total emissions to almost 15.68% compared to the 2005 level. The total emissions obtained from the IEEE 13-node system dropped to 4273 CO2 tons for this scenario. Figure 16 illustrates the generated emissions after applying this scenario.   Figure 15. Residential load profile and generation mix from the IEEE test system for 2017, after using residential DR with ESS. Figure 16. Emissions from the electricity residential sector from the IEEE test system for 2017, after using residential DR with ESS.

Communication and Intelligence with Energy Conservation
This scenario considers merging residential DR with energy conservation programs to examine all the potential options that could reduce the carbon footprints from the electric sector. In fact, energy conservation can play a major role in attaining the goal when it is integrated with residential DR. In Figure 17, we realize a reduction in the demand curve due to the aggressive curtailment in energy use. Figure 18 depicts the emissions from the electric residential sector for 2017. It shows that this combination resulted in more than 18% emission reduction, dropping from 5068.67 CO2 tons in 2005 to 4151.43 CO2 tons. Figure 16. Emissions from the electricity residential sector from the IEEE test system for 2017, after using residential DR with ESS.

Communication and Intelligence with Energy Conservation
This scenario considers merging residential DR with energy conservation programs to examine all the potential options that could reduce the carbon footprints from the electric sector. In fact, energy conservation can play a major role in attaining the goal when it is integrated with residential DR. In Figure 17, we realize a reduction in the demand curve due to the aggressive curtailment in energy use. Figure 18 depicts the emissions from the electric residential sector for 2017. It shows that this combination resulted in more than 18% emission reduction, dropping from 5068.67 CO 2 tons in 2005 to 4151.43 CO 2 tons. Sustainability 2021, 13, x FOR PEER REVIEW 12 of 25 Figure 15. Residential load profile and generation mix from the IEEE test system for 2017, after using residential DR with ESS. Figure 16. Emissions from the electricity residential sector from the IEEE test system for 2017, after using residential DR with ESS.

Communication and Intelligence with Energy Conservation
This scenario considers merging residential DR with energy conservation programs to examine all the potential options that could reduce the carbon footprints from the electric sector. In fact, energy conservation can play a major role in attaining the goal when it is integrated with residential DR. In Figure 17, we realize a reduction in the demand curve due to the aggressive curtailment in energy use. Figure 18 depicts the emissions from the electric residential sector for 2017. It shows that this combination resulted in more than 18% emission reduction, dropping from 5068.67 CO2 tons in 2005 to 4151.43 CO2 tons. Figure 17. Residential load profile and generation mix from the IEEE test system for 2017, after using residential DR with energy conservation.  Figure 17. Residential load profile and generation mix from the IEEE test system for 2017, after using residential DR with energy conservation. Figure 18. Emissions from the electricity residential sector from the IEEE test system for 2017, after using residential DR with energy conservation.

Communication and Intelligence with ESS and Energy Conservation
The last scenario examines all three options (i.e., residential DR, ESS, and energy conservation) in one scenario and then analyzes their emissions reduction. Each scenario is applied with its designed methodology that is explained in the above subsections. This scenario provides more flexible options during the year and, in turn, achieves more reduction in emissions. The results obtained show significant energy saving after integrating these options. It is notable from Figure 19 that this scenario reduced the total energy used in 2017 from the residential sector from 5.36 GWh for the base case to 4.97 GWh. Therefore, this scenario successfully reduces emissions from the 2005 level by about 19.72%, as detailed in Figure 20. Figure 18. Emissions from the electricity residential sector from the IEEE test system for 2017, after using residential DR with energy conservation.

Communication and Intelligence with ESS and Energy Conservation
The last scenario examines all three options (i.e., residential DR, ESS, and energy conservation) in one scenario and then analyzes their emissions reduction. Each scenario is applied with its designed methodology that is explained in the above subsections. This scenario provides more flexible options during the year and, in turn, achieves more reduction in emissions. The results obtained show significant energy saving after integrating these options. It is notable from Figure 19 that this scenario reduced the total energy used in 2017 from the residential sector from 5.36 GWh for the base case to 4.97 GWh. Therefore, this scenario successfully reduces emissions from the 2005 level by about 19.72%, as detailed in Figure 20.  Figure 17. Residential load profile and generation mix from the IEEE test system for 2017, after using residential DR with energy conservation. Figure 18. Emissions from the electricity residential sector from the IEEE test system for 2017, after using residential DR with energy conservation.

Communication and Intelligence with ESS and Energy Conservation
The last scenario examines all three options (i.e., residential DR, ESS, and energy conservation) in one scenario and then analyzes their emissions reduction. Each scenario is applied with its designed methodology that is explained in the above subsections. This scenario provides more flexible options during the year and, in turn, achieves more reduction in emissions. The results obtained show significant energy saving after integrating these options. It is notable from Figure 19 that this scenario reduced the total energy used in 2017 from the residential sector from 5.36 GWh for the base case to 4.97 GWh. Therefore, this scenario successfully reduces emissions from the 2005 level by about 19.72%, as detailed in Figure 20.

Results and Discussion
The figures here compared all the scenarios and quantified the progress of each scenario in achieving the goal. In that regard, Figure 21 shows the demand curve from the IEEE test system for 2017 after examining the first three scenarios: ESS, energy conservation, and the scenario that combines ESS with energy conservation. Figure 22 illustrates the residential demand curve from the test system after considering communication and intelligence alternative and merging it with other scenarios.

Results and Discussion
The figures here compared all the scenarios and quantified the progress of each scenario in achieving the goal. In that regard, Figure 21 shows the demand curve from the IEEE test system for 2017 after examining the first three scenarios: ESS, energy conservation, and the scenario that combines ESS with energy conservation. Figure 22 illustrates the residential demand curve from the test system after considering communication and intelligence alternative and merging it with other scenarios.

Results and Discussion
The figures here compared all the scenarios and quantified the progress of each scenario in achieving the goal. In that regard, Figure 21 shows the demand curve from the IEEE test system for 2017 after examining the first three scenarios: ESS, energy conservation, and the scenario that combines ESS with energy conservation. Figure 22 illustrates the residential demand curve from the test system after considering communication and intelligence alternative and merging it with other scenarios.   The results also showed the environmental footprint was different for each scenario. Figure 23 demonstrates the emissions from the emitted sources after integrating ESS and energy conservation. Table 1 explained the percentage of emissions reduction from each scenario. ESS successfully reduced more than 18% of emissions. However, another study found that combining ESS with solar PV could only provide an emissions reduction of 16% [23]. These findings pointed out the fact that several factors could affect the penetration of ESS in reducing emissions. The values indicated that the MCDM options can help in a meaningful reduction in total emissions. Merging ESS with energy conservation provided a potential improvement in environmental footprint.  The results also showed the environmental footprint was different for each scenario. Figure 23 demonstrates the emissions from the emitted sources after integrating ESS and energy conservation. Table 1 explained the percentage of emissions reduction from each scenario. ESS successfully reduced more than 18% of emissions. However, another study found that combining ESS with solar PV could only provide an emissions reduction of 16% [23]. These findings pointed out the fact that several factors could affect the penetration of ESS in reducing emissions. The values indicated that the MCDM options can help in a meaningful reduction in total emissions. Merging ESS with energy conservation provided a potential improvement in environmental footprint. The results also showed the environmental footprint was different for each scenario. Figure 23 demonstrates the emissions from the emitted sources after integrating ESS and energy conservation. Table 1 explained the percentage of emissions reduction from each scenario. ESS successfully reduced more than 18% of emissions. However, another study found that combining ESS with solar PV could only provide an emissions reduction of 16% [23]. These findings pointed out the fact that several factors could affect the penetration of ESS in reducing emissions. The values indicated that the MCDM options can help in a meaningful reduction in total emissions. Merging ESS with energy conservation provided a potential improvement in environmental footprint.   By including residential DR as a representative substitute for communication and intelligence, the residential sector saved up to 20% of the 2005 electric residential sector emissions. Specifically, the residential DR scenario achieved a reduction of 13.91% in carbon emissions, while the work in [24] found that DR diminished emissions by up to 14%. A study has been conducted in Europe to investigate the impact of intelligent DSM on controlling and operating electrical devices. The study used three scenarios: low, medium, and high levels. The findings show that GHG emissions can be minimized by up to 10% [25]. Thus, energy savings and emissions reduction might vary depending on the mechanism and the DR techniques. Figure 24 explains the emissions per source from the electric residential sector.  By including residential DR as a representative substitute for communication and intelligence, the residential sector saved up to 20% of the 2005 electric residential sector emissions. Specifically, the residential DR scenario achieved a reduction of 13.91% in carbon emissions, while the work in [24] found that DR diminished emissions by up to 14%. A study has been conducted in Europe to investigate the impact of intelligent DSM on controlling and operating electrical devices. The study used three scenarios: low, medium, and high levels. The findings show that GHG emissions can be minimized by up to 10% [25]. Thus, energy savings and emissions reduction might vary depending on the mechanism and the DR techniques. Figure 24 explains the emissions per source from the electric residential sector. The analysis on the residential sector concluded that merging all the MCDM options can help Fort Collins in meeting its 2020 goal and providing the highest reduction in emissions. Table 2 presents carbon emissions reduction from the DR scenarios. The analysis on the residential sector concluded that merging all the MCDM options can help Fort Collins in meeting its 2020 goal and providing the highest reduction in emissions. Table 2 presents carbon emissions reduction from the DR scenarios. However, this analysis was implemented based on an environmental evaluation while there are other factors that can affect this ranking. In that regard, the next section will perform an economic evaluation for all alternatives using a CBA technique.

CBA
CBA quantifies the net benefit to the system of each scenario. This analysis evaluates the costs and benefits associated with each solution proposed to reduce emissions. The analysis quantifies environmental, technical, and economic costs and benefits from customer and utility perspectives. This evaluation should explain all costs associated with each alternative, such as fixed costs, operating and maintenance costs, and customer dropout and removal costs. Benefits to the utility can be reductions in the costs of generation, transmission, and distribution. Additionally, there can be customer benefits such as reducing the cost of electricity bills or utility benefits, lowering the cost of services or improving operation and efficiency. Social benefits such as reducing environmental degradation, conserving resources, or protecting the global environment are also considered. A key issue is the monetization of environmental and social impacts. A social cost of carbon emissions and a carbon price resulting from climate policies are two proposed approaches to establish monetary values for emissions [26,27].
CBA follows several steps to select the optimal alternative. The first step is to specify the set of alternative options. The next step is to define the boundaries of the analysis to explain which benefits and costs are included.
Step three is to select indicators before measuring all costs and benefits of the selected measurements. Then, we monetize all costs and benefits before applying a discount rate to calculate the NPV, the IRR, the payback period, and BCR to select the most economically beneficial alternative for reducing emissions in this particular distribution system. NPV computes all expected cash flows associated with the project and subtracts them from the capital cost [28][29][30].
where T is the period of the project, r is the discount rate, CF 0 is investment cost, CF t is the net cash flow at time t. The project is accepted if NPV > 0 while the decision maker rejects the project when NPV < 0. IRR is another metric that measures the profitability of the project. So, the IRR is the discount rate that makes the NPV equal to zero. The payback period is the time required to recover the capital investment in a project. The cash flow is summed until it equals the initial investment in the project. BCR is another useful measure that summarizes the relationship between the related costs and benefits of the project. The benefits exceed the costs when BCR > 1 and the project should be accepted. When BCR = 1, the project indicates that the costs and benefits are equal and the project can be accepted with little viability. If BCR < 1, it means the costs are higher than the benefits and the project should be rejected [28][29][30].

BCR = Present value o f bene f its Present value o f costs
The next section investigates the benefits and costs of integrating DSM alternatives programs into the distribution system. Such analysis quantifies all the associated costs and measures the estimated benefits for each scenario. Since each scenario has its cost parameters, the costs vary based on the specified alternative. Further, the benefits obtained from every scenario consider the revenue from the avoided operating and maintenance costs or bill reduction. It also monetizes the impact of every scenario on power quality, reliability, environmental collateral, and socioeconomic equity. Results will rank each scenario, based on a weighted value, thus allowing a decision maker to consider tradeoffs.

ESS
As explained in Section 3.2, ESS is dispatched during peak hours to shave the demand and avoid running resources such as coal. While ESS has a capacity of 50 MW, 200 MWh, and serving for a 4 h peak load, this capacity is adjusted to match the IEEE 13-node system capacity. The new size of ESS is scaled down to 270.98 kW for 4 h of peak shaving (i.e., 1083.93 kWh). In that regard, PRPA proposes the estimated cost range of ESS for peak shaving. The battery cost ranges between USD 340 and USD 450 per kWh. Power conversion system costs USD 150-350 per kW, while power control system starts from USD 80-120 per kW. Balance of plant could cost USD 90-120 per kW and procurement construction costs USD 150-180 per kWh [17]. Further, this study adjusts the cost of recycling to the size of the ESS used in the test system [31]. Since this is an economic evaluation for new technology in the system, we consider the worst case scenario. Hence, the analysis considers the highest price for each parameter. A research study in [32] estimates the discount rate for a storage system project is 5.09%, while PRPA expects a 10 year lifetime for the ESS project. Further, the EIA estimates the average load growth from 2018 to 2050 is 0.2%, which is the load growth assumed in this calculation [33].
The benefits of installing ESS are savings in fuel costs, network support, and environmental benefits. Reference [32] concludes that integrating a storage system in the grid would improve DG integration by increasing its utilization. EIA determines a levelized avoided cost for resources such as coal and natural gas. This calculation uses USD 0.082 and USD 0.080 per kWh for coal and natural gas, respectively, to estimate the fuel cost savings [34]. There is a social cost of using fossil-fueled generating units in which the environmental cost is USD 11 per kg of CO 2 and the cost will increase to reach USD 26 per kg of CO 2 by 2050. Thus, as a result of enabling ESS, the amount of avoided emissions is multiplied by USD 11 to quantify the environmental benefits [35]. The economic benefit of utilizing DG is USD 52.28 per MWh, which is avoided from conventional resources [32]. Maintaining acceptable limits of power quality and reliability could save up to USD 62.71 per kW. Table 3 summarizes the findings from the economic evaluation. The project is beneficial since BCR is greater than 1 and the project will pay its capital cost after about 4 years.

Energy Conservation
According to the DOE, the cost of implementing a conservation program in Fort Collins includes installing about 85,328 smart meters, 2347 programmable communicating thermostats, and 1710 direct load control devices [36]. The investment cost of this energy conservation project is adjusted to the test system. Further, there is an operating cost of USD 0.035 per kWh saved to deploy such a project [37]. Implementing energy conservation programs can cause inconvenience for the participants. Therefore, this analysis defines socioeconomic cost as a societal cost. EPA estimates socioeconomic cost as USD 0.214 per kWh [38]. Moreover, Fort Collins has a rebate program to replace low efficiency equipment with higher efficiency appliances [39]. The benefits obtained from reducing the energy demand and mitigating environmental footprints are calculated as explained in Section 4.1. Reference [36] shows the avoided and deferred costs after implementing the conservation project while the socioeconomic benefit is obtained from [38], after excluding avoided carbon cost. Table 4 demonstrates the final calculation for the energy conservation project. Although the project is acceptable, the project takes a longer time than ESS to recover its expenses.

ESS with Energy Conservation
The economic analysis of this scenario includes all the expected costs from ESS and energy conservation. We notice that the final cost after 10 years is higher than expected. This is reasonable since this scenario has a capital cost and higher operating costs. After calculating the outcomes of this scenario, it saves energy more than expected since ESS and efficient appliances are displacing the energy needed from other resources. The results are not surprising since this scenario incorporates more cost. The project needs more than 6.5 years to pay its investment cost. However, it is still economically acceptable. Table 5 gives an economic summary about merging ESS and conservation programs.

Communication and Intelligence
As previously explained, the residential DR program represents communication and intelligence since this alternative leads to a change in the end-use demand curve. The economic evaluation of this option is implemented according to the DR model in [29]. Since this alternative is deployed in the residential sector, the capital cost is excluded and scaled down from [36]. The energy cost, energy sales, and peak demand cost for the base case are calculated in (4) [29].
Energy cost = Q * π r where Q is the energy consumption in kWh and π r is the retail price in USD per kWh. PRPA determines the energy charges as USD 0.04282 per kWh for the summer season and USD 0.04109 per kWh for the winter season [40]. The energy sales before DR can be calculated in (5). Energy sales = D * π wh (5) where D is the demand in kW and π wh is the wholesale price in USD per kWh. The utility charges the customers for energy usage based on the old electric rates for residential energy use, before DR. Table 6 shows the usage charge per kWh, regenerated from [41]. Monthly peak demand is the user's demand during the hour that coincides with the system's monthly peak. A charge of 11.56 USD per kW is applied as the demand charge for the summer months. The tariff decreases to USD 8.81 per kW for non-summer months [40].
The utility applies a 60 min charge on coincident demand as shown in (6). where PCD is the peak demand charge in USD/h, P m is the coincident peak demand in kWh, and π DC is the demand charge is USD. The energy cost, energy sales, and peak demand cost for the DR are calculated for the modified energy consumption in kWh, Q .
Energy cost = Q * π r The electricity sold after DR, given in (8), depends on the new demand curve in kWh, D , and the pricing mechanism, explained in Section 3.1.
Energy sales = D * π wh (8) While the peak demand charge, PCD , changes according to the new coincident peak demand, Pm , the monthly peak demand rates remain the same during coincident demand.
The cost of incorporating such a technique requires evaluating the capital cost and the variable costs. The operating and maintenance costs of applying DR are adjusted as 28 USD per kW, to match the size of the test system [42]. Equation (10) shows all the associated costs of DR where DRinv is the investment cost of DR and CO&M is the annual operating and maintenance cost.
Total cost = DRinv + CO&M (10) This alternative incurs a yearly financial benefit by obtaining the difference between the energy sales and the peak demand charges.
The analysis shows several costs and benefits from the program. As smart meters and thermostat devices are the enablers of residential DR, the analysis adjusted the expected benefits, such as reduced cost and investment deferral, from [36]. The economic evaluation shows the project is economically beneficial to reduce the environmental impact. Table 7 explains the economic outcomes of this alternative. Although the project requires annual expenses to implement the program, it takes less time than energy efficiency to recover the capital cost with less benefit.

Communication and Intelligence with ESS
This subsection investigates the economic approach of combining residential DR and ESS. All the costs and benefits follow the methodologies explained in Sections 3.2 and 3.5. Since ESS is integrated into the residential electricity sector, the investment cost is extracted for the residential sector only, and the DR framework model is the same. The results explain how DR and ESS work together to achieve such benefits. The results are not surprising since the main goal of these alternatives is shaving the peak during high electricity demand. Even though this scenario takes a longer time than the residential DR project to pay the capital cost, Table 8 shows this option is more beneficial than the previous alternative.

Communication and Intelligence with Energy Conservation
This scenario applies conservation programs along with residential DR. The conservation program uses the model in Section 3.3. The cost of implementing this might be higher since there is an effect on the convenience of the participants. After specifying all the costs and benefits obtained from this partnership between conservation programs and DR, the findings demonstrate this project has the longest time to recover the cost and to become profitable. Table 9 indicates that this project is accepted based on its BCR value but with less benefit than the other alternatives.

Communication and Intelligence with ESS and Energy Conservation
The last scenario investigates the economics of the approach of combining all the proposed alternatives on the electricity residential sector. Environmentally, this option successfully achieved about a 20% reduction in total emissions from the electricity residential sector. However, this subsection studies the economic viability of this approach. In that regard, the analysis takes into consideration the methodologies in Sections 3.2, 3.3 and 3.5. The results demonstrate that this scenario performs better than the base case. Integrating all the available options has increased the potential benefits as well as the convenience level for the participants since ESS discharges the stored energy during peaking demand. The economic results show the improvement after integrating ESS with the two remaining alternatives. The NPV increased as well as BCR while the payback period decreased, as stated in Table 10.

Results and Discussion
The pertinent results to the above-studied scenarios indicated that all the projects were economically accepted. The available options provided positive NPVs, paid their costs during the life of the project, and the BCRs were greater than one. The analysis created a combination between the environmental impact and the economic outputs for the alternatives. The findings showed that ESS has the highest rank among other alternatives in the case of economic preference. This is reasonable since ESS charges during lower electricity prices and discharges during peak demand hours. Thus, the customers can avoid the high electricity charges and the utility does not need to run more generating units. The second-best option is residential DR combined with ESS. This alternative boasted flexibility in shifting the loads to off-peak periods. Thus, the results showed more than 15% emissions reduction and it pays back after 4.09 years. Another study in [43] matched our environmental and economic analysis when aggregating DR with a battery system provided up to 14% GHG emission reduction while the payback period was 4.84 years. Combining all the MCDM alternatives in one option was the third-ranked scenario. This alternative provided more emissions reduction than the previous ones. Although the energy conservation option and ESS with energy conservation project provided less economic impact than residential DR, those two options ranked in fourth and fifth place, respectively, due to their environmental impact. The penultimate alternative was residential DR because it was designed to shift the peak load and it had a socioeconomic cost. The last alternative was combining residential DR with energy conservation. Although this option performed environmentally better than residential DR alone, its socioeconomic cost played a major role in selecting this alternative.
However, priorities might change according to the participants' choices. One may prefer environmental impact over economic output or vice versa. Therefore, a trade-off between the options might lead to a strategic decision that takes into consideration the advantages and disadvantages of each project. In this analysis, two factors can affect the prioritization: environmental impact and economic output. Thus, the result is presented in Figure 25 in a trade-off setup so decision makers can compare the alternatives to select the most suitable option.
alternatives. The findings showed that ESS has the highest rank among other alternatives in the case of economic preference. This is reasonable since ESS charges during lower electricity prices and discharges during peak demand hours. Thus, the customers can avoid the high electricity charges and the utility does not need to run more generating units. The second-best option is residential DR combined with ESS. This alternative boasted flexibility in shifting the loads to off-peak periods. Thus, the results showed more than 15% emissions reduction and it pays back after 4.09 years. Another study in [43] matched our environmental and economic analysis when aggregating DR with a battery system provided up to 14% GHG emission reduction while the payback period was 4.84 years. Combining all the MCDM alternatives in one option was the third-ranked scenario. This alternative provided more emissions reduction than the previous ones. Although the energy conservation option and ESS with energy conservation project provided less economic impact than residential DR, those two options ranked in fourth and fifth place, respectively, due to their environmental impact. The penultimate alternative was residential DR because it was designed to shift the peak load and it had a socioeconomic cost. The last alternative was combining residential DR with energy conservation. Although this option performed environmentally better than residential DR alone, its socioeconomic cost played a major role in selecting this alternative.
However, priorities might change according to the participants' choices. One may prefer environmental impact over economic output or vice versa. Therefore, a trade-off between the options might lead to a strategic decision that takes into consideration the advantages and disadvantages of each project. In this analysis, two factors can affect the prioritization: environmental impact and economic output. Thus, the result is presented in Figure 25 in a trade-off setup so decision makers can compare the alternatives to select the most suitable option.

Conclusions and Future Work
This paper presented an environmental analysis by creating several scenarios. The simulation studied each option individually and then merged an option with other available options. The IEEE 13-node test system was used through the OpenDSS simulation tool to analyze 2017 hourly yearly load data, corresponding to the Fort Collins, Colorado, area. Each scenario shows different CO 2 curves based on the change in supply and demand. The base case simulation for the 2017 data pointed out that the emissions were reduced by 16.26% from the 2005 level. Thus, this complied with the results released by the city. However, studying this case from a test system would not be enough to rank the potential solutions for minimizing their environmental footprints through DSM.
This work investigated the economic consequences of the MCDM options. It studied each alternative based on its associated costs and the expected benefits. Therefore, the economic evaluation obtains the NPV, the IRR, the payback period, and the BCR for every scenario. Since this work examines the potential emissions reduction from the proposed alternatives, it ensured those solutions are economically accepted. The preliminary result shows that investing in ESS is the best option, followed by combining residential DR with ESS.
The results presented in this dissertation rank the alternatives based on the simulation from the IEEE 13-node distribution system. However, studying this case from a test system would not be sufficient in ranking the potential solutions for minimizing environmental footprints through DSM. The future path should consider the real electric distribution system of Fort Collins. As the city has a plan to reduce its dependence on conventional generation and to increase the penetration of renewable energy, analyzing the real system can obtain results that are more accurate. The future use of the test system does not constrain the results of the current studies as the future work will investigate reliability, power quality, and sizing and siting of ESS in the electric distribution system. This can increase the benefits of ESS by decreasing the system loss and increasing reliability and power quality. Our future path is also considering coordinating energy storage charging with available renewable energy generation. This dispatch will increase the utility utilization by smoothing energy output from intermittent resources and will increase the reliability and resilience of the system. It can also reduce electricity bills, defer investment, and offset the emissions from dirty generating units.