Impact on Renewable Design Requirements of Net-Zero Carbon Buildings under Potential Future Climate Scenarios

This paper presents an analysis to foresee renewable design requirement changes of net- zero carbon buildings (NZCBs) under different scenarios of potential future climate scenarios in the U.S. Northeast and Midwest regions. A climate change model is developed in this study using the Gaussian random distribution method with monthly temperature changes over the whole Northeast and Midwest regions, which are predicted based on a high greenhouse gas (GHG) emission scenario (i.e., the representative concentration pathways (RCP) 8.5). To reflect the adoption of NZCBs potential in future, this study also considers two representative future climate scenarios in the 2050s and 2080s of climate change years in the U.S. Northeast and Midwest regions. An office prototype building model integrates with an on-site photovoltaics (PV) power generation system to evaluate NZCB performance under the climate change scenarios with an assumption of a net-metering electricity purchase agreement. Appropriate capacities of the on-site PV system needed to reach NZCB balances are determined based on the building energy consumption impacted by the simulated climate scenarios. Results from this study demonstrated the emission by electricity consumption increases as moving toward the future scenarios of up to about 25 tons of CO2-eq (i.e., about 14% of the total CO2-eq produced by the electricity energy source) and the PV installation capacity to offset the emission account for the electricity consumption increases significantly up to about 40 kWp (i.e., up to more than 10% of total PV installation capacities) as the different climate scenarios are applied. It is concluded that the cooling energy consumption of office building models would significantly impact GHG emission as future climate scenarios are considered. Consequently, designers of NZCBs should consider high performance cooling energy systems in their designs to reduce the renewable energy generation system capacity to achieve net-zero carbon emission goals.


Introduction
The building sector accounts for approximately 40% of the total end-use energy consumption (i.e., primary and secondary energy usages) in the United States [1], including about 21% and 18% for residential and commercial buildings, respectively. In the commercial building sector, about 74% of the energy consumption was captured for the electric power sector to generate and supply electricity retails. The U.S. Energy Information Administration (EIA) also estimated that electricity use in the commercial building sector was expected to increase by about 14% by 2050 compared to levels in 2019 [2]. This growth has led to concerns regarding the future increased energy consumption and human emission of carbon dioxide (CO 2 ) from the building sector. According to the Congressional Research Service reports [3,4], it was pointed out that the electric power sector contributed the second-largest percentage, 35% of U.S. CO 2 emission contributions based on the 2017 U.S. CO 2 emission contribution data by sectors. Within the building sector, the commercial sector presented more than 37% of the electric power sector's CO 2 contributions. Although the amount of electricity power generation had remained flatted after 2010, CO 2 emissions continued a general trend of reduction after 2007 mainly because of increased renewable sources and natural gas use as well as reduced coal use in electricity generation, resulting in a 20% reduction of CO 2 emissions when compared to the CO 2 emissions levels of 2005 [3,4]. However, there are still concerns regarding greenhouse gas (GHG) emission and global warming issues as the global growth rate of GHG emission from human activities has still increased [5][6][7].
As an effort to reduce GHG emission, the concept of net-zero carbon emission buildings (NZCBs) has been developed, and the idea has been gaining significance globally to achieve feasible long-term goals of GHG emission reductions for the coming decades [8]. To appropriately implement the NZCB concept, high energy-efficient appliances and systems must be encouraged [9], and an appropriate selection and design of on-site renewable energy systems must be incorporated to effectively reduce both on-site and central generation of GHG emission [10,11]. Louwen et al. [12] presented a good overview of the impact of GHG emission by PV-based renewable energy generation. They reviewed over 40 years of PV development and analyzed energy demand and GHG emission impacts in terms of PV production. Their results illustrated that when the doubling of installed PV capacity occurred, energy consumption and GHG emission could be decreased by about 13% and 17-24%, respectively. Pinel et al. [13] evaluated the impact of allowing to buy CO 2 compensation to enable the zero emission neighborhoods design. They pointed out that there were large differences in energy system design and GHG emission generation depending on variations in weather conditions and the price of external compensation.
Researchers have conducted numerous studies on climate change impacts on building energy consumption. Wan et al. [14] investigated building energy use variations based on climate change scenarios in subtropical climates. They demonstrated that improving efficiencies in building components, such as lighting and chiller, can provide potential mitigation for increased future energy use due to climate changes in the late 21st century. Nik and Kalagasidis [15] assessed possible climate changes and uncertainties in future energy performance of the residential building stock in Stockholm. Four uncertainty factors (i.e., global climate models, regional climate models, emissions scenarios, and initial conditions) were considered in their study, and they pointed out that all their climate change scenarios up to 2100 showed that the future heating demands would be decreased by about 30% compared to the date before 2011 while the cooling demand would increase mainly because of increased outdoor air temperatures. Kikumoto et al. [16] constructed future standard weather data for use in building design and energy analysis for 2030s (i.e., 2031-2035). Based on their study, it was observed that there would be a substantial impact of climate change on the energy performance of a detached house, including 15% increases in sensible heat loads predicted under their study condition. Shen [17] evaluated the energy consumption variation of residential and commercial buildings by climate change scenarios in four representative cities in the U.S. using EnergyPlus. They showed that climate changes could have potential impacts on residential and office building energy usages during the years 2040-2069, including −1.64% to 14.07% and −3.27 to −0.12% of annual energy changes for residential and office buildings, respectively. Wang et al. [18] presented energy simulation-based analysis to investigate impacts of climate change on building energy consumptions in U.S. climate zones. The EnergyPlus medium office prototype model was used to predict energy consumption with two different climate change models, including Hadley Centre Coupled Model version 3 and NCAR Community Earth System Model version 1. With their generated future weather files, they examined that climate change mitigation measures related to heating, ventilation, and air conditioning (HVAC) operations such as adjustment of thermostat setpoints and reduced HVAC operation hours. In addition, Summa et al. [19] compared energy consumption of a current residential NZEB in Rome, Italy with that expected in 2050 under Representative Concentration Pathway (RCP) 4.5 and 8.5 scenarios. Their results showed that about 50% of increases in cooling energy use were presented, while heating energy use showed significant decreases under the future climate changes.
Although there have been several pieces of research investigating building energy use due to climate changes, there is still a lack of studies understanding the adoption of NZCBs associated with future potential climate changes in a quantitative manner in the literature [20]. To fill this research gap, this paper investigates the impact of climate changes on commercial building energy consumption and carries out an analysis to foresee renewable design requirements of NZCBs under different scenarios of future climate changes in the U.S. Northeast and Midwest regions. The medium-sized office prototype building model, which complies with a widely and globally accepted commercial building energy standard (i.e., ASHRAE Standard 90.1-2019 [21]), is subjected in this study to evaluate NZCB design requirements under different climate change scenarios. The focus of this paper is to provide useful insights into variations in building energy consumption and renewable design requirements when office buildings are subjected to a net zero carbon design under the current typical meteorological year (TMY) data and future climate scenarios. It is important to note that this paper does not investigate different strategies for enabling net zero carbon operations to commercial buildings, e.g., advanced building envelope and HVAC technologies to reduce building energy consumption.

Description of Climate Change Model
This section describes the methodology to generate climate change data for building energy simulations. Various weather conditions, including dry-bulb temperature, wet-bulb temperature, solar radiation, relative humidity, and wind speed/direction, can impact the building thermal load. Such factors can be considered in comprehensive building energy modeling and prediction to estimate energy consumption in a building. Typical meteorological year 3 (TMY3) weather data [22], which is statistically derived annual climate data, is used in the building energy simulations of this study. Our previous study [23] revealed that the ambient temperature is a dominant contributing factor to the thermal load during both heating and cooling seasons among all weather variables. Therefore, predicted changes in the dry-bulb temperature of outdoor air are assumed to be a single dominant variable for the climate change model in this study. A stochastic method using the normal distribution with a random sampling technique is applied to generate future weather profiles for building energy simulation. The probability density function of the normal distribution is defined as: where µ is the mean value and σ is the standard deviation, and the standard deviation is assumed to be one (1) in this study. As the value of σ increases, the range of potential random errors also increases, therefore, this study considers appropriate distribution patterns by assuming the standard deviation of the predicted dry-bulb temperature to be 1 • C to ensure that the adjacent hourly temperature does not deviate more than 2 • C. The mean values of the predicted dry-bulb temperature are determined based on the study conducted by Byun et al. [24]. The Byun's model [24] was developed based on hydrologic modeling experiments over 20 Midwestern watersheds using the Variable Infiltration Capacity (VIC) model forced by historical observed datasets and future projections from the statistically downscaled Global Climate model (GCM). Their model provides monthly projected changes in hydrometeorolgical variables over the whole Midwest and Northeast regions, including precipitation, outdoor air temperature, evapotranspiration, and soil moisture. The projected monthly temperature changes from the Byun's model are used to generate future Climate 2021, 9, 17 4 of 20 weather files for all selected locations in this study. Figure 1 presents the temperature changes of the Byun's model for two representative scenarios of future study periods (i.e., 2050s and 2080s). The future weather projections in Byun's study are based on two GHG concentration scenarios, including Representative Concentration Pathway (RCP) 4.5 and 8.5 for medium and high cases, respectively. The two scenarios are based on 30-year windows centered on the 2050s (2041-2070) and 2080s (2071-2100).
Climate 2021, 9, x FOR PEER REVIEW 4 of 20 changes in hydrometeorolgical variables over the whole Midwest and Northeast regions, including precipitation, outdoor air temperature, evapotranspiration, and soil moisture. The projected monthly temperature changes from the Byun's model are used to generate future weather files for all selected locations in this study. Figure 1 presents the temperature changes of the Byun's model for two representative scenarios of future study periods (i.e., 2050s and 2080s). The future weather projections in Byun's study are based on two GHG concentration scenarios, including Representative Concentration Pathway (RCP) 4.5 and 8.5 for medium and high cases, respectively. The two scenarios are based on 30-year windows centered on the 2050s (2041-2070) and 2080s (2071-2100).

Figure 1.
Monthly projected changes in the outdoor air temperature [24]. RCP: Representative concentration pathways.
The RCP 8.5 model's projection shown in Figure 1 is selected in this study to generate random distributions of temperature changes with the scenarios with an assumption of a high-emission scenario. Monthly mean values of outdoor air temperature changes are determined from the Byun's RCP 8.5 model shown in Table 1. These mean values are used in the probability density function in Equation (1)  Three locations, including New York, NY, Buffalo, NY, and Rochester, MN, in the U.S. Northeast and Midwest regions are chosen to evaluate the impact of climate change on the building's energy use. Those locations are selected because they can represent three different climate zones in the U.S. Northeast and Midwest regions defined by ASHRAE 90.1-2019, i.e., climate zones 4A, 5A, and 6A, which can show variations in net zero carbon requirement and design within the same region. Hourly temperature changes of the scenario in the 2050s and 2080s are randomly generated separately for each climate location using the mean values in Table 1 and Equation (1). Figures 2 and 3 show monthly distributions of randomly generated hourly temperature changes for the scenario in the 2050s and 2080s, respectively, in climate location 4A New York, NY. Note that the dotted line in each plot indicate the mean value of each month, and each bin size in the histogram is determined based on the number of hourly temperature change data in each month and scenario. Figure 4 shows climate change patterns of hourly temperatures corresponding to each scenario in the selected locations. The RCP 8.5 model's projection shown in Figure 1 is selected in this study to generate random distributions of temperature changes with the scenarios with an assumption of a high-emission scenario. Monthly mean values of outdoor air temperature changes are determined from the Byun's RCP 8.5 model shown in Table 1. These mean values are used in the probability density function in Equation (1) to generate random distributions of hourly temperature variations for each climate change scenario in different climate locations. Table 1 summarizes monthly mean values of outdoor air temperature changes extracted from the Byun's RCP 8.5 model. Three locations, including New York, NY, Buffalo, NY, and Rochester, MN, in the U.S. Northeast and Midwest regions are chosen to evaluate the impact of climate change on the building's energy use. Those locations are selected because they can represent three different climate zones in the U.S. Northeast and Midwest regions defined by ASHRAE 90.1-2019, i.e., climate zones 4A, 5A, and 6A, which can show variations in net zero carbon requirement and design within the same region. Hourly temperature changes of the scenario in the 2050s and 2080s are randomly generated separately for each climate location using the mean values in Table 1 and Equation (1). Figures 2 and 3 show monthly distributions of randomly generated hourly temperature changes for the scenario in the 2050s and 2080s, respectively, in climate location 4A New York, NY. Note that the dotted line in each plot indicate the mean value of each month, and each bin size in the histogram is determined based on the number of hourly temperature change data in each month and scenario. Figure 4 shows climate change patterns of hourly temperatures corresponding to each scenario in the selected locations.       This study uses typical meteorological year version 3 (TMY3) weather data [25] for the hourly building energy simulations. The TMY3 weather dataset was created using data collected over a 1991-2005 period based on the procedures developed by National Renewable National Laboratories, and its dataset covers more than 1000 locations in the U.S. [25]. This study assumes that the TMY3 weather data would represent typical weather conditions in the present time in the selected climate locations, and it is used as a baseline weather scenario for climate change analysis. To generate the weather scenario in the 2050s and 2080s, TMY3 weather datasets are modified to reflect the climate changes shown in Figure 4. More specifically, the hourly outdoor air dry-bulb temperatures in the TMY3 weather data file in each location are modified to determine the weather scenario in the 2050s and 2080s. The modified hourly dry-bulb temperatures (TOA,modified,i) for the weather scenario in the 2050s and 2080s can be determined by: where , , is the hourly outdoor air dry-bulb temperature of current TMY3 weather files and , is the hourly temperature changes, presented in Figure 4.

Description of Building Model
The Department of Energy (DOE)'s flagship building energy modeling software, En-ergyPlus version 9.0, is used to model office buildings in the selected climate locations [22]. In this study, the DOE's medium-sized office prototype building model is used to represent all office building types. The office building type is chosen in this study because the building model has the highest construction weights among all other building types (i.e., construction weight of 14.98% based on 2003-2007 surveyed data in the U.S.) and complies with the minimum energy code requirements prescribed in ASHRAE Standard 90.1-2019 [25,26]. Figure 5 shows a three-dimensional (3D) perspective view and a floor plan of the DOE's medium-sized office prototype model. This study uses typical meteorological year version 3 (TMY3) weather data [25] for the hourly building energy simulations. The TMY3 weather dataset was created using data collected over a 1991-2005 period based on the procedures developed by National Renewable National Laboratories, and its dataset covers more than 1000 locations in the U.S. [25]. This study assumes that the TMY3 weather data would represent typical weather conditions in the present time in the selected climate locations, and it is used as a baseline weather scenario for climate change analysis. To generate the weather scenario in the 2050s and 2080s, TMY3 weather datasets are modified to reflect the climate changes shown in Figure 4. More specifically, the hourly outdoor air dry-bulb temperatures in the TMY3 weather data file in each location are modified to determine the weather scenario in the 2050s and 2080s. The modified hourly dry-bulb temperatures (T OA,modified,i ) for the weather scenario in the 2050s and 2080s can be determined by: where T OA,original,i is the hourly outdoor air dry-bulb temperature of current TMY3 weather files and T change,i is the hourly temperature changes, presented in Figure 4.

Description of Building Model
The Department of Energy (DOE)'s flagship building energy modeling software, Ener-gyPlus version 9.0, is used to model office buildings in the selected climate locations [22]. In this study, the DOE's medium-sized office prototype building model is used to represent all office building types. The office building type is chosen in this study because the building model has the highest construction weights among all other building types (i.e., construction weight of 14.98% based on 2003-2007 surveyed data in the U.S.) and complies with the minimum energy code requirements prescribed in ASHRAE Standard 90.1-2019 [25,26]. Figure 5 shows a three-dimensional (3D) perspective view and a floor plan of the DOE's medium-sized office prototype model. The prototype office building is configured as a three-story building with a rectangular floor plan, as shown in Figure 5. Each floor has five (5) thermal zones, including four (4) perimeter zones and one (1) core zone, with a total floor area of 4982 m². Its construction consists of steel-framed exterior walls, an insulation layer entirely above deck for a flat roof, and slab-on-grade floors. All the construction details and internal heat gains comply with the minimum energy code requirements, i.e., ASHRAE 90.1-2019 [21]. Table  2 summaries representative input parameter and its characteristics of the office model. In terms of a HVAC system, the original prototype office model employs a rooftop unit (RTU) system with variable air volume (VAV) and electric reheat box to meet required heating and cooling loads in each zone of a building, as shown in Figure 5. With a RTU-VAV system, natural gas and electricity are used for a central heating coil and a reheat coil of a VAV reheat box, respectively. For cooling, a direct expansion (DX) cooling coil with two-speed fan operation is used, and the cooling coil requires two sets of performance data to calculate electricity consumption based on the fan speed. All the same performance curves that the original prototype model provides are used in this study. Table  3 summarizes detailed information of the HVAC systems used in the prototype building model.  The prototype office building is configured as a three-story building with a rectangular floor plan, as shown in Figure 5. Each floor has five (5) thermal zones, including four (4) perimeter zones and one (1) core zone, with a total floor area of 4982 m 2 . Its construction consists of steel-framed exterior walls, an insulation layer entirely above deck for a flat roof, and slab-on-grade floors. All the construction details and internal heat gains comply with the minimum energy code requirements, i.e., ASHRAE 90.1-2019 [21]. Table 2 summaries representative input parameter and its characteristics of the office model. In terms of a HVAC system, the original prototype office model employs a rooftop unit (RTU) system with variable air volume (VAV) and electric reheat box to meet required heating and cooling loads in each zone of a building, as shown in Figure 5. With a RTU-VAV system, natural gas and electricity are used for a central heating coil and a reheat coil of a VAV reheat box, respectively. For cooling, a direct expansion (DX) cooling coil with two-speed fan operation is used, and the cooling coil requires two sets of performance data to calculate electricity consumption based on the fan speed. All the same performance curves that the original prototype model provides are used in this study. Table 3 summarizes detailed information of the HVAC systems used in the prototype building model.  Modeling the energy impacts of the building's internal heat gains in EnergyPlus requires assumptions about the internal heat gain intensity and operation schedules. All input values and schedules regarding internal heat gains are directly obtained from the original version of a prototype office model without modifications. Figure 6 shows fraction schedules of representative internal heat gains [26,27] (i.e., lighting, plug, and occupancy), which are used to determine total internal heat gain values inside by simply multiplying the internal heat gain density times its faction value in the time horizon. In addition to the internal heat gains, ASHRAE Standard 90.1-2019 [21] requires additional plug loads consumed by the office building, such as exterior lighting and electricity used by an elevator. All the parts of building energy consumption are considered for NZCB design in this study.  Modeling the energy impacts of the building's internal heat gains in EnergyPlus requires assumptions about the internal heat gain intensity and operation schedules. All input values and schedules regarding internal heat gains are directly obtained from the original version of a prototype office model without modifications. Figure 6 shows fraction schedules of representative internal heat gains [26,27] (i.e., lighting, plug, and occupancy), which are used to determine total internal heat gain values inside by simply multiplying the internal heat gain density times its faction value in the time horizon. In addition to the internal heat gains, ASHRAE Standard 90.1-2019 [21] requires additional plug loads consumed by the office building, such as exterior lighting and electricity used by an elevator. All the parts of building energy consumption are considered for NZCB design in this study.

Net-Zero Carbon Building (NZCB) Design
Net-zero GHG emission building criteria can be achieved when the annual site GHG emission is equal or less than zero. The net GHG emission (EMnet) can be calculated by: where , and , indicate CO emission values of the whole building electricity and national gas end-use, respectively, and i varies from 1 to 8760 [hr] for an annual simulation.
, is the saved CO emission value due to the on-site electricity generation. Each emission value can be calculated by multiplying the building electricity and natural gas end-uses times each conversion factor as follows:

Net-Zero Carbon Building (NZCB) Design
Net-zero GHG emission building criteria can be achieved when the annual site GHG emission is equal or less than zero. The net GHG emission (EM net ) can be calculated by: where EM elec−used,i and EM gas−used,i indicate CO 2 emission values of the whole building electricity and national gas end-use, respectively, and i varies from 1 to 8760 [hr] for an annual simulation. EM elec−gen,i is the saved CO 2 emission value due to the on-site Climate 2021, 9, 17 9 of 20 electricity generation. Each emission value can be calculated by multiplying the building electricity and natural gas end-uses times each conversion factor as follows: EM elec−gen = E PV−gen × CF elec (6) where E elec−used and E gas−used are total building electricity and natural gas end-uses, respectively. CF elec and CF gas are emission conversion factors for electricity and natural gas, respectively, in each state. Table 4 shows the U.S. CO 2 emission conversion factors by fuel used for this study. A solar photovoltaic (PV) system is one of the most promising distributed power technologies and can be used on building footprints and building site areas for on-site electricity generation. Although there are many ways to integrate PV systems with a building, this study only considers open area and roof mount types for PV installation with a fixed tilt angle because this comparative study focuses only on climate changes and GHG emission reductions. To enable a NZCB balance, both roof and parking areas are used for solar PV installation locations. This study uses the PVWatts module in EnergyPlus to calculate on-site electricity generation. The PVWatts model is originally developed by the National Renewable Energy Laboratory (NREL), and is a web application and module in NREL's system advisor model software that estimates the electricity generation of a grid-connected PV system for buildings based on a few simple inputs. EnergyPlus uses the same calculation algorithm without any additional needs to call the web application during the calculation process. The hourly available PV production is defined [29] as: where P dc,i and P dc0 indicate hourly generated DC power from the plane-of-array and specified nameplate DC rate, I tr,i is the transmitted irradiance on a cell surface, γ is the temperature coefficient, and T cell,i and T re f are the computed cell temperature at each time step using a first-principles heat transfer energy balance model developed by Fuentes [30] and the reference cell temperature, respectively. The reference cell temperature is 25 • C, and reference irradiance is 1000 W/m 2 [29,31]. The input values of the PVWatts module used in this study are listed in Table 5.  This study chooses the "standard" option because this option represents typical ployor mono-crystalline silicon modules, with efficiencies in the range of 14-17%, which is commonly used for roof mounted PV installation. Losses in the PV system include the impacts of soiling, shading, snow cover, mismatch, wiring, connections, light-induced degradation, nameplate rating, system age, and operational availability [31,32]. This study uses 14% of the default total system loss, which is provided by the EnergyPlus PVWatts model. The inverter model is based on an analysis of California Energy Commission (CEC) inverter performance data [33]. Hourly AC power outputs can be determined from the DC outputs [29] as: where η is the inverter efficiency, P dc,i is the hourly generated DC power outputs from the plane-of-array, P dc0 is the specified nameplate DC rate, and P ac0 is the AC nameplate rating determined from the DC rating of the system and the DC-to-AC ratio. P ac0 can be calculated by simply multiplying P dc0 times η 96% of the default nominal efficiency used for the calculation in this study. The electric load center transformer in EnergyPlus allows for feeding the surplus power from on-site power generators to the electrical grid. The load center object matches the voltage between on-site power generation and total electricity used by a building facility. All surplus electricity produced by the PV power generators can be fed back into the electrical grid through the main panel operation that is connected to both the utility grid and on-site power generators of the NZCB. For the transformer performance, the nominal efficiency is set to 97.7%, which is directly obtained from EnergyPlus prototype building's default value [32]. Figure 7 depicts the on-site PV power generators with net-metering operation for the NZCB model. This study chooses the "standard" option because this option represents typical ployor mono-crystalline silicon modules, with efficiencies in the range of 14-17%, which is commonly used for roof mounted PV installation. Losses in the PV system include the impacts of soiling, shading, snow cover, mismatch, wiring, connections, light-induced degradation, nameplate rating, system age, and operational availability [31,32]. This study uses 14% of the default total system loss, which is provided by the EnergyPlus PVWatts model. The inverter model is based on an analysis of California Energy Commission (CEC) inverter performance data [33]. Hourly AC power outputs can be determined from the DC outputs [29] as: = 0 where is the inverter efficiency, , is the hourly generated DC power outputs from the plane-of-array, is the specified nameplate DC rate, and is the AC nameplate rating determined from the DC rating of the system and the DC-to-AC ratio.
can be calculated by simply multiplying times 96% of the default nominal efficiency used for the calculation in this study.
The electric load center transformer in EnergyPlus allows for feeding the surplus power from on-site power generators to the electrical grid. The load center object matches the voltage between on-site power generation and total electricity used by a building facility. All surplus electricity produced by the PV power generators can be fed back into the electrical grid through the main panel operation that is connected to both the utility grid and on-site power generators of the NZCB. For the transformer performance, the nominal efficiency is set to 97.7%, which is directly obtained from EnergyPlus prototype building's default value [32]. Figure 7 depicts the on-site PV power generators with netmetering operation for the NZCB model.

Simulation Setting Description Based on the Scenarios
The detailed steps of the simulation setting based on the two scenarios are depicted in Figure 8. The original EnergyPlus prototype office models of each climate location are first simulated with modified weather files of the different scenarios to evaluate its electricity and national gas consumption. CO2 emission values are then calculated based on each building's electricity-and natural gas-based emissions.

Simulation Setting Description Based on the Scenarios
The detailed steps of the simulation setting based on the two scenarios are depicted in Figure 8. The original EnergyPlus prototype office models of each climate location are first simulated with modified weather files of the different scenarios to evaluate its electricity The total maximum PV power outputs installed on both the roof and parking areas of the simulated office model are presented in Table 6 to enable NZCB performance. About 83% of the entire roof area (1661 m²) is used for the roof-mounted PV installation. Additional PV panels, which are about 3986.4 m², are installed on parking areas until each simulated model shows satisfactory annual NZCB performance corresponding to each scenario.
It should also be noted that only on-site PV systems are applied to the existing prototype office building model with no technological advance changes to enable NZCBs because the objective of this study is only to investigate the future climate change's impacts on electricity and natural gas energy consumption and thus the potential need of PV capacity for NZCB performance. This simplified way can give precise insight into how the current office building model itself reacts to the future climate changes for the potential enabling of a zero carbon emission building with the scenarios over the coming decades.

Results and Discussion
This section presents and discusses results obtained from building energy simulations of NZCBs under the different climate change scenarios in the U.S. Northeast and Midwest regions. There are two climate change scenarios updated from TMY3 weather data based on future changes obtained from the proposed climate change model. Those scenarios are applied in a conventional office building model (i.e., DOE' medium office The total maximum PV power outputs installed on both the roof and parking areas of the simulated office model are presented in Table 6 to enable NZCB performance. About 83% of the entire roof area (1661 m 2 ) is used for the roof-mounted PV installation. Additional PV panels, which are about 3986.4 m 2 , are installed on parking areas until each simulated model shows satisfactory annual NZCB performance corresponding to each scenario.
It should also be noted that only on-site PV systems are applied to the existing prototype office building model with no technological advance changes to enable NZCBs because the objective of this study is only to investigate the future climate change's impacts on electricity and natural gas energy consumption and thus the potential need of PV capacity for NZCB performance. This simplified way can give precise insight into how the current office building model itself reacts to the future climate changes for the potential enabling of a zero carbon emission building with the scenarios over the coming decades.

Results and Discussion
This section presents and discusses results obtained from building energy simulations of NZCBs under the different climate change scenarios in the U.S. Northeast and Midwest regions. There are two climate change scenarios updated from TMY3 weather data based on future changes obtained from the proposed climate change model. Those scenarios are applied in a conventional office building model (i.e., DOE' medium office prototype building model) to estimate energy consumption and GHG emission as a baseline. Then the baseline office building model is designed to be the NZCB with the PV-based renewable energy generation system. A comparative analysis is conducted to determine the capacity of the on-site PV system to enable net zero carbon balance. The following sections present the climate data obtained from the proposed climate change model, energy performance, GHG emission results for climate change scenarios determined from the prototype office building models, and the potential need for on-site PV system capacities to enable the NZCBs with the climate change scenarios over a yearly cycle.

Climate Change Data
In this section, the profiles of generated hourly dry-bulb temperatures are presented. Each climate location has three different temperature profiles, including the original TMY3. The hourly temperature profiles of the two scenarios are developed based on Equation (2), Figure 9 depicts the hourly profiles of new generated outdoor air dry-bulb temperatures for two scenarios and three climate locations. The temperature rises in hourly temperature profiles throughout all scenarios and climate locations can be recognized because of the rises in monthly projected outdoor air temperatures, as shown in Figure 1. Those temperature profiles, as shown in Figure 9, generated by the random distribution are used to create new TMY3 weather files and the weather flies are then applied to building energy simulation for the comparative analysis of NZCBs. Note that the impact of other parameters, such as solar radiation and wind speed, in the TMY weather file is not reflected in this study when creating new weather files because the try bulb temperature is the most influential parameter factor to evaluate building energy consumption, especially when HVAC energy usages are focused.
Climate 2021, 9, x FOR PEER REVIEW 12 of 20 prototype building model) to estimate energy consumption and GHG emission as a baseline. Then the baseline office building model is designed to be the NZCB with the PVbased renewable energy generation system. A comparative analysis is conducted to determine the capacity of the on-site PV system to enable net zero carbon balance. The following sections present the climate data obtained from the proposed climate change model, energy performance, GHG emission results for climate change scenarios determined from the prototype office building models, and the potential need for on-site PV system capacities to enable the NZCBs with the climate change scenarios over a yearly cycle.

Climate Change Data
In this section, the profiles of generated hourly dry-bulb temperatures are presented. Each climate location has three different temperature profiles, including the original TMY3. The hourly temperature profiles of the two scenarios are developed based on Equation (2), Figure 9 depicts the hourly profiles of new generated outdoor air dry-bulb temperatures for two scenarios and three climate locations. The temperature rises in hourly temperature profiles throughout all scenarios and climate locations can be recognized because of the rises in monthly projected outdoor air temperatures, as shown in Figure 1. Those temperature profiles, as shown in Figure 9, generated by the random distribution are used to create new TMY3 weather files and the weather flies are then applied to building energy simulation for the comparative analysis of NZCBs. Note that the impact of other parameters, such as solar radiation and wind speed, in the TMY weather file is not reflected in this study when creating new weather files because the try bulb temperature is the most influential parameter factor to evaluate building energy consumption, especially when HVAC energy usages are focused.  Based on the dry-bulb temperatures shown in Figure 9, the design day weather data were derived by determining the maximum and minimum temperatures of the scenario in the 2050s and 2080s in design-day months. The design day information is used for sizing the HVAC systems in Section 2.2. The design day information for the baseline and scenario in the 2050s and 2080s for the selected locations are provided in Table 7.

Energy Use Intensity (EUI) of Baseline Model
This section compares the annual whole building EUIs of conventional baseline building models for the original weather file versus the two scenarios in the three U.S. climate locations. There are two energy source types in the prototype office building model, including electricity and natural gas. Space heating and domestic hot water (DHW) systems use natural gas as an energy source, while all other parts use electricity in the model. Figure 10 shows the detailed breakdown of annual building EUIs with the unit of kWh/m 2 -year. As expected, EUIs from HVAC systems (i.e., heating and cooling) vary significantly depending on the climate changes while EUIs for the other part remains almost the same. For example, in the 5A climate condition, increases of the cooling EUIs are 64% and 122% for the scenario in the 2050s and 2080s, respectively, whereas the heating EUIs show 31% and 47% of reductions corresponding to each scenario when compared to the original weather condition. The results from the other two climate locations (i.e., 4A and 6A) also indicate that similar trends are followed. The results indicate that there are significant impacts of climate changes on the HVAC energy consumption, thus it could be observed that GHG emission would vary primarily depending on HVAC energy usages. However, it should be noted that the total annual building EUIs do not show large differences when compared to the original case as shown in Figure 10 because the cooling EUI increases as the heating EUI decreases. Also note that the heating EUI consists of two fuel types, including electricity and natural gas, from VAV reheat coils and central heating coils, respectively. Since it is expected that changes in GHG emission generated by buildings primarily depend on heating and cooling energy usages, GHG emission variation by energy source types is also analyzed in this study. The annual EUIs of baseline models by the energy source type are presented in Figure 11 with the unit of kWh/m 2 -year. The results of all the climate locations indicate that the annual electricity usages tend to increase, as expected, while the annual gas usages show a decreased trend corresponding to climate change scenarios. Climate 2021, 9, x FOR PEER REVIEW 14 of 20 Figure 10. Annual whole building Energy Use Intensity (EUI) of baseline models. Figure 11. Annual electricity and natural gas energy use intensity of baseline models.

GHG Emission of Baseline Model by Energy Source Type
Annual GHG emission values of each energy source type for the scenarios and climate locations in the building are determined by Equation (4) and Equation (5), respectively. Figures 12 and 13 present annual GHG emission values of baseline building models by energy source types and HVAC parts. In Figure 12, the GHG emission values from the electricity consumption show increased patterns and that the natural gas consumption shows reduced values as the scenario of climate changes moves further towards the 2050s and 2080s, whereas the total building EUIs do not vary significantly. For example, in climate 4A, increases in the GHG emission values from the electricity are 6.7 and 13.0 tons of CO2-eq for the 2050s and 2080s scenario, respectively, when compared to the original TMY3 weather condition. Reductions in those from the natural gas are 6.8 and 9.92 tons of CO2-eq for the 2050s and 2080s scenarios, respectively, when compared to the original weather condition. The results from other two climate locations (i.e., 5A and 6A) also show similar patterns of increases and decreases in GHG emission values, but the amount of GHG emission values of the 6A location is much higher than the others due to the emission conversion factor of electricity. Figure 13 depicts annual GHG emission values of baseline models by HVAC components for each scenario and location with the unit of ton CO2-eq. From this figure, changes in the GHG emission by HVAC fan energy usages are shown to be relatively low when compared to the original TMY3. However, the changes of GHG emission in HVAC heating and cooling energy usages are significant. In climate 5A, increases in the GHG emission values from cooling energy use are 5.3 and 10.0 tons of CO2-eq for scenarios in the

GHG Emission of Baseline Model by Energy Source Type
Annual GHG emission values of each energy source type for the scenarios and climate locations in the building are determined by Equation (4) and Equation (5), respectively. Figures 12 and 13 present annual GHG emission values of baseline building models by energy source types and HVAC parts. In Figure 12, the GHG emission values from the electricity consumption show increased patterns and that the natural gas consumption shows reduced values as the scenario of climate changes moves further towards the 2050s and 2080s, whereas the total building EUIs do not vary significantly. For example, in climate 4A, increases in the GHG emission values from the electricity are 6.7 and 13.0 tons of CO2-eq for the 2050s and 2080s scenario, respectively, when compared to the original TMY3 weather condition. Reductions in those from the natural gas are 6.8 and 9.92 tons of CO2-eq for the 2050s and 2080s scenarios, respectively, when compared to the original weather condition. The results from other two climate locations (i.e., 5A and 6A) also show similar patterns of increases and decreases in GHG emission values, but the amount of GHG emission values of the 6A location is much higher than the others due to the emission conversion factor of electricity. Figure 13 depicts annual GHG emission values of baseline models by HVAC components for each scenario and location with the unit of ton CO2-eq. From this figure, changes in the GHG emission by HVAC fan energy usages are shown to be relatively low when compared to the original TMY3. However, the changes of GHG emission in HVAC heating and cooling energy usages are significant. In climate 5A, increases in the GHG emission values from cooling energy use are 5.3 and 10.0 tons of CO2-eq for scenarios in the

GHG Emission of Baseline Model by Energy Source Type
Annual GHG emission values of each energy source type for the scenarios and climate locations in the building are determined by Equation (4) and Equation (5), respectively. Figures 12 and 13 present annual GHG emission values of baseline building models by energy source types and HVAC parts. In Figure 12, the GHG emission values from the electricity consumption show increased patterns and that the natural gas consumption shows reduced values as the scenario of climate changes moves further towards the 2050s and 2080s, whereas the total building EUIs do not vary significantly. For example, in climate 4A, increases in the GHG emission values from the electricity are 6.7 and 13.0 tons of CO 2 -eq for the 2050s and 2080s scenario, respectively, when compared to the original TMY3 weather condition. Reductions in those from the natural gas are 6.8 and 9.92 tons of CO 2 -eq for the 2050s and 2080s scenarios, respectively, when compared to the original weather condition. The results from other two climate locations (i.e., 5A and 6A) also show similar patterns of increases and decreases in GHG emission values, but the amount of GHG emission values of the 6A location is much higher than the others due to the emission conversion factor of electricity.
2050s and 2080s, respectively, when compared to the original case. In terms of GHG emission from HVAC fan energy use in 5A, it indicates that there are 0.3 and 0.5 tons of CO2eq for scenarios in the 2050s and 2080s, respectively, when compared to the original case. Note that the GHG emission values of electricity and national gas energy usages presented in this section are used as the baseline to further develop to NZCBs.

GHG Emission and PV Installation for NZCB
This section discusses variations of GHG emission values when the prototype office building models achieve the performance of NZCBs corresponding to each scenario and climate location. The PV installation capacities needed to enable NZCB goals is also presented in this section. Figure 14 depicts the scatter plots of each NZCB performance corresponding to each scenario and location with the unit of ton CO2-eq. The annual GHG emission from the baseline building energy end-use is compared against that from on-site power generation to investigate carbon emission balances for each scenario and location. The ratio of y-axis to x-axis should be one or greater to reach zero carbon emission performance, as shown in Figure 14. It is observed that all the cases of the building model present acceptable net zero carbon balances. The GHG emissions of NZCBs are in the ranges of 80-120 tons for climate zones 4A and 5A and 220-260 tons for climate zone 6A, respectively. Since the GHG emission of the whole building's energy end-use is determined 2050s and 2080s, respectively, when compared to the original case. In terms of GHG emission from HVAC fan energy use in 5A, it indicates that there are 0.3 and 0.5 tons of CO2eq for scenarios in the 2050s and 2080s, respectively, when compared to the original case. Note that the GHG emission values of electricity and national gas energy usages presented in this section are used as the baseline to further develop to NZCBs.

GHG Emission and PV Installation for NZCB
This section discusses variations of GHG emission values when the prototype office building models achieve the performance of NZCBs corresponding to each scenario and climate location. The PV installation capacities needed to enable NZCB goals is also presented in this section. Figure 14 depicts the scatter plots of each NZCB performance corresponding to each scenario and location with the unit of ton CO2-eq. The annual GHG emission from the baseline building energy end-use is compared against that from on-site power generation to investigate carbon emission balances for each scenario and location. The ratio of y-axis to x-axis should be one or greater to reach zero carbon emission performance, as shown in Figure 14. It is observed that all the cases of the building model present acceptable net zero carbon balances. The GHG emissions of NZCBs are in the ranges of 80-120 tons for climate zones 4A and 5A and 220-260 tons for climate zone 6A, respectively. Since the GHG emission of the whole building's energy end-use is determined Figure 13. Annual GHG emission of baseline models by HVAC components. Figure 13 depicts annual GHG emission values of baseline models by HVAC components for each scenario and location with the unit of ton CO 2 -eq. From this figure, changes in the GHG emission by HVAC fan energy usages are shown to be relatively low when compared to the original TMY3. However, the changes of GHG emission in HVAC heating and cooling energy usages are significant. In climate 5A, increases in the GHG emission values from cooling energy use are 5.3 and 10.0 tons of CO 2 -eq for scenarios in the 2050s and 2080s, respectively, when compared to the original case. In terms of GHG emission from HVAC fan energy use in 5A, it indicates that there are 0.3 and 0.5 tons of CO 2 -eq for scenarios in the 2050s and 2080s, respectively, when compared to the original case. Note that the GHG emission values of electricity and national gas energy usages presented in this section are used as the baseline to further develop to NZCBs.

GHG Emission and PV Installation for NZCB
This section discusses variations of GHG emission values when the prototype office building models achieve the performance of NZCBs corresponding to each scenario and climate location. The PV installation capacities needed to enable NZCB goals is also presented in this section. Figure 14 depicts the scatter plots of each NZCB performance corresponding to each scenario and location with the unit of ton CO 2 -eq. The annual GHG emission from the baseline building energy end-use is compared against that from on-site power generation to investigate carbon emission balances for each scenario and location. The ratio of y-axis to x-axis should be one or greater to reach zero carbon emission performance, as shown in Figure 14. It is observed that all the cases of the building model present acceptable net zero carbon balances. The GHG emissions of NZCBs are in the ranges of 80-120 tons for climate zones 4A and 5A and 220-260 tons for climate zone 6A, respectively. Since the GHG emission of the whole building's energy end-use is determined based on electricity and natural gas usages, both electricity and gas consumed in a building need to be considered to calculate acceptable NZCB balances. Figure 14 also presents the near NZCB balances that considers only electricity consumption. based on electricity and natural gas usages, both electricity and gas consumed in a building need to be considered to calculate acceptable NZCB balances. Figure 14 also presents the near NZCB balances that considers only electricity consumption. Equation (3) is used to determine the total PV installation capacities needed for the NZCB goals by offsetting all GHG emissions from the electricity and natural gas consumption of the building model. Figure 15 depicts the capacities of on-site PV installation and its GHG emission intensities to enable the zero carbon emission balances. This figure shows that the capacity of PV installation to offset the GHG emission from the electricity consumption is significantly higher than that from the natural gas consumption in all climate locations. Figure 15. PV installation capacity and its GHG emission intensity to enable NZCBs.
As expected, the PV capacities of covering the electricity's GHG emission tend to increase gradually according to the climate change scenarios mainly due to rises in outdoor air temperatures. For example, in Buffalo (climate zone 4A), 314, 327, and 341 kWp of PV Equation (3) is used to determine the total PV installation capacities needed for the NZCB goals by offsetting all GHG emissions from the electricity and natural gas consumption of the building model. Figure 15 depicts the capacities of on-site PV installation and its GHG emission intensities to enable the zero carbon emission balances. This figure shows that the capacity of PV installation to offset the GHG emission from the electricity consumption is significantly higher than that from the natural gas consumption in all climate locations. based on electricity and natural gas usages, both electricity and gas consumed in a building need to be considered to calculate acceptable NZCB balances. Figure 14 also presents the near NZCB balances that considers only electricity consumption. Equation (3) is used to determine the total PV installation capacities needed for the NZCB goals by offsetting all GHG emissions from the electricity and natural gas consumption of the building model. Figure 15 depicts the capacities of on-site PV installation and its GHG emission intensities to enable the zero carbon emission balances. This figure shows that the capacity of PV installation to offset the GHG emission from the electricity consumption is significantly higher than that from the natural gas consumption in all climate locations. Figure 15. PV installation capacity and its GHG emission intensity to enable NZCBs.
As expected, the PV capacities of covering the electricity's GHG emission tend to increase gradually according to the climate change scenarios mainly due to rises in outdoor air temperatures. For example, in Buffalo (climate zone 4A), 314, 327, and 341 kWp of PV Figure 15. PV installation capacity and its GHG emission intensity to enable NZCBs.
As expected, the PV capacities of covering the electricity's GHG emission tend to increase gradually according to the climate change scenarios mainly due to rises in outdoor air temperatures. For example, in Buffalo (climate zone 4A), 314, 327, and 341 kWp of PV installation capacities are shown for the climate scenarios of original TMY3, 2050s, and 2080s, respectively. Conversely, for the PV capacities to cover GHG emission from the national gas usage, it is shown that the 2080s scenarios present the lowest capacities, which are relatively higher reductions when compared to other cases. For example, in the climate 4A case, 46 kWp of PV installation capacities is needed to offset the GHG emission from natural gas, whereas about 60 kWp of PV installation capacity is needed for the other scenarios (i.e., TMY3 and 2050s) in the location. The results from the other two climate locations (i.e., 5A and 6A) also show similar patterns of increases and decreases in the PV installation capacities to offset GHG emission values. In addition, Figure 15 presents GHG emission generated from PV installations. A total of 50 g CO 2 -eq/kWh of GHG intensity obtained from [34,35] is used to calculate the GHG emission of PV installation in this study. These results indicate that the emission of PV installation is about one order of magnitude lower than the emission produced by the building energy consumption.

Conclusions
This study investigated the impacts of potential future climate changes on annual building energy use intensities and GHG emission in the U.S. Northeast and Midwest climate locations. It also analyzed renewable design requirements to enable NZCBs under future climate scenarios using an on-site PV system. The Gaussian random distribution method was used to develop climate change scenarios with monthly temperature change values over the whole Northeast and Midwest regions, which were developed based on a high GHG emission scenario (i.e., the representative concentration pathways (RCP) 8.5). A simulation-based study of a prototype office building was conducted to evaluate the annual electricity and natural gas usages of the building models in different climate locations in the U.S. Northeast and Midwest regions. The PV-based on-site power generation system was considered and applied to the simulated building model to enable NZCB performance by assuming net-metering operation under future climate scenarios. Appropriate capacities of the on-site PV power system for reaching NZCB balances were required to be determined according to the energy source types of building energy consumption. The amount of GHG emission from the electricity or natural gas usages of the office building model could vary significantly depending on future climate scenarios. The capacities of the on-site PV systems for NZCBs could also vary based on a combined effect of annual building energy usages in each location and climate scenario, the emission conversion factors, and solar irradiation availability in each climate location. Key findings from this study are as follows:

•
The simulated office building models that comply with the minimum energy code requirements of ASHRAE Standard 90.1-2019 presented around 79-84 kWh/m 2 -year and 13-26 kWh/m 2 -year for the annual electricity and natural gas energy usages, respectively, in U.S. Northeast and Midwest climate locations. As the future climate scenarios of 2050s and 2080s were applied, the amount of annual electricity consumption tended to increase gradually from the baseline year climate conditions, whereas that of annual natural gas use decreased gradually toward the future years. This is primarily because of increases in outdoor air temperatures of climate change data used in the building simulations, which led to increases in the cooling energy consumption in the building, decreasing in heating energy usages; • The GHG emission trends in each climate followed the same pattern as the trends of the energy source type energy consumption (i.e., electricity and natural gas) changes toward the future climate scenarios. The emission by electricity consumption increased by, moving toward the future scenarios, up to about 25 tons of CO 2 -eq (about 14% of the total CO 2 produced by the electricity energy source) while the emission produced by the natural gas consumption in the building decreased by up to about 20 tons of CO 2 -eq (about 30% of the total CO 2 -eq produced by the natural gas energy source).
The total GHG emission of the building in each climate location remained relatively even in all climate scenarios, although each climate location showed a different level of the emission rate. Depending on the climate conditions and electricity emission conversion factors in different climate locations, the same building can generate more emissions as the climate scenarios toward 2050s and 2080s as shown in the results of the New York case (climate zones 4A and 5A) in Figures 13 and 14; • The capacities of on-site PV systems required to enable NZCBs can differ depending on the local climate conditions (i.e., temperature and irradiation) and climate change scenarios as shown in shown in Figure 15 although the total PV capacity responsible by all energy source types did not change significantly (i.e., less than about a 15 kWp increase of the total PV capacity to offset the GHG emissions in the different climate scenarios). The building models of climate zone 4A used around 345-360 kWp for the zero-carbon emission target, while the other two climates (i.e., 5A and 6A) showed similar trends of PV capacities, including 376-387 kWp and 380-396 kWp for climate 5A and 6A, respectively. The results in Figure 15 also showed that the PV installation capacity to offset the emission account for the electricity consumption increases significantly up to about 40 kWp (i.e., up to more than 10% of total PV installation capacities) as the different climate change scenarios were applied.
Such results from this study could provide useful insights into not only GHG emission generated from the electricity and national gas consumption of office buildings under different future climate scenarios in the selected climate locations, but also renewable design requirements (i.e., the capacity of on-site PV systems) when the office buildings were considered for the NZCB design. From the results, it can be concluded that the cooling energy consumption would significantly impact GHG emissions as future potential climate scenarios are considered. Consequently, the designers of NZCBs should consider highly efficient cooling energy systems in their designs to reduce the renewable energy generation system capacity to achieve net-zero carbon emission goals.
It should be noted that there are limitations in this study, including that (a) only outdoor air dry-bulb temperature was used as a weather data parameter for developing the climate change model, and (b) the studied climate regions were limited within the U.S. Northeastern and Midwestern areas mainly due to the limited information of future climate change data in other regions of the U.S. Based on this study, future work can be extended to include various climate parameters in the whole building simulation analysis to reflect climate change scenarios in higher resolutions and improve the estimation of building energy consumptions and renewable energy requirements for NZCBs. Furthermore, all U.S. climate locations could be considered for a national scale analysis of NZCBs.

Data Availability Statement:
The data presented in this study are available within this article.