Analytical approach
This study aims to explore the impacts and implications of system deployment pathways on deploying cost-effective green hydrogen energy systems for different urban communities. Therefore, energy system design should cover all potential community energy input scenarios subject to multiple uncertainties to ensure robust and optimal outcomes. However, climate change exhibits a distinct probabilistic nature, i.e., extreme climate conditions are significantly less likely to occur than typical climate conditions. The energy input scenarios should reflect this variability in weather probabilities to ensure reasonable realism. Therefore, we use stochastic optimization to formulate and integrate the objective function into the energy model31. In this process, each energy scenario is treated independently, and then their results are combined based on the probability of occurrence to obtain the mathematical expectation of the outcome.
The overview of modelling multi-source uncertainties
Energy models are subject to inherent input parameter uncertainties derived from complex real-world communities. It is crucial for energy models to transparently disclose the impacts of community uncertainties on energy system inputs and provide reliable energy input scenarios for subsequent evaluation. Parameterizing the compound impacts of these uncertainties on energy input scenarios for the energy model is challenging due to the complexity of sources and the differences in their influence pathways. Thus, we distinguish the community uncertainties into two levels: climate-human uncertainties at the household level and community design layer uncertainties at the community level (Supplementary Notes 1). As a bottom-up energy model, we quantify uncertainties layer by layer from the household level to the community level and superimpose them to quantify compound impacts.
Modelling climate-human uncertainties
To parameterize the climate-human uncertainties at the household level, we distinguish climate and human systems into two treatments based on their different sources. We then associate them into compound scenarios to account for their superimposed impacts on energy input. Climate uncertainty stems from the unpredictable nature of climate period paths and their corresponding extreme weather events, leading to the loss of potential risk information on weather-dependent energy production and demand in defined climate scenarios32. Current climate change further increases uncertainty, exacerbating climate risk to energy systems due to the increasing frequency of extreme climate events and ongoing shifts in climate patterns33. Thus, we developed a method to synthesize future representative weather datasets based on regional climate models (RCMs) data to capture potential climate uncertainty information related to weather patterns and extreme weather events (Supplementary Note 2). This method can derive an arbitrary number of stochastic hourly-resolution weather data to decrease the simulation time without destroying month-scale weather features. Thus, this method maximizes our ability to generate enough weather data to measure climate uncertainty information and then trims these data to representative weather scenario sets with occurrence probabilities, which facilitates compatibility with stochastic energy system optimization.
Human uncertainty comes from the complex stochastic nature of human behavior and their interaction with energy systems, in which the greater control freedom for the residential occupants makes the uncertainty of energy demand profiles even more significant. Deterministic occupant behavior variables not only lead to overestimating or underestimating energy demand profiles at the household level but also decrease the diversity of energy demand profiles at the community level. Thus, in human systems, we sample human-related operational information from assumed probability distributions and parameter ranges based on urban statistical data, ignoring the behavioral logic of the occupants17. We explore seven aspects of human behavior uncertainties that affect building energy demand profiles. These features encompass the full range of human-involved operational parameters in building energy simulations. Due to the aleatory nature and weak relationship of uncertainty from the approach, we specified independent probability functions and then drew random samples from these features using Latin Hypercube Sampling (Supplementary Note 3)34,35. We specify uniform distribution as probability functions to equal probabilities of each interval in uncertainty variables because there is no explicit reason to value one probability distribution over another. Meanwhile, we extracted the boundary of these variables from the urban residential energy surveys36. Every random human behavior scenario was then associated with each climate scenario to create a compound scenario pool to treat the uncertainties in both climate and humans properly.
In this study, we created eight sets of composite climate-human uncertainty scenario pools in North America (climate zones 1 ~ 7), where each scenario pool yielded a total of 50,000 household-level production and demand profiles, i.e., 50 climate scenarios and 1,000 human scenarios.
Modelling regional household benchmark buildings
Besides climate-human scenarios, an urban residential archetype is crucial in generating household energy scenarios using dynamic energy simulations. We extracted physics-based benchmark building models based on local urban information: urban building survey data and building standards of climate zones (Supplementary Note 4). In this study, we created a set of benchmark bungalows and two-story houses in each city for aggregating artificial urban communities. We combined building information modeling (BIM) and building energy modeling (BEM) to generate household energy scenarios systematically. Specifically, BIM crafts benchmark archetypes with Level of Detail-2 (LoD2), incorporating fundamental construction information such as building typologies, floors, roofs, external walls, and windows, all sourced from urban building surveys. Concurrently, BEM analyzes building energy output by associating benchmark archetypes and climate-human scenarios. This process incorporates energy characteristic factors like envelope construction, internal loads, HVAC systems, and renewable energy systems, with specifications adhering to building standards within the local climate zone. Once completed, the EnergyPlus engine executes energy simulations to analyze hourly-resolution dynamic energy performance.
Modelling community design layer uncertainties
Parameterizing uncertainties within the community design layer in bottom-up energy models entails creating detailed community models that encompass the full array of urban community features. Extracting community samples from cities to parameterize uncertainties falls into a dilemma between scenario scopes and feature completeness due to the intricate nature of urban community features. Due to the lack of boundary samples, the representative samples derived from regional-average community features may not provide complete uncertainty information for extrapolation to all urban communities. Thus, we create a set of artificial communities using community characteristic matrix: household scale and prosumer scale. The household scale quantifies uncertainties in community composition using boundary and mean urban community features, while the prosumer scale accounts for uncertainties in renewable energy using linear scale expansion. These matrix parameters define the boundary and node samples to mimic community feature information, which provides fully quantitative information on community design layer uncertainties. In this study, the community scale was assumed to be 50 households. Each community characteristic matrix for the local climate zone considered three types of household scales and five sets of prosumer scales.
The community design layer uncertainties were parameterized in two sequential processes. Initially, we clustered aggregate household energy scenarios into designed community energy scenarios. Subsequently, we assigned household roles as either prosumers or consumers within these community energy scenarios. In the first process, we translated each type of household energy scenario into corresponding shares of community energy scenarios based on household scale using k-means clustering. Subsequently, in the second step, to mitigate the impact of the direction of renewable energy penetration, we arranged households within the community in descending order of energy sharing levels. This assumed that households with higher levels of sharing were more inclined to adopt renewable energy systems. Next, we categorized households into prosumer and consumer roles based on prosumer scale. Finally, this pipeline of uncertainty propagation returned a total of 750 community energy scenarios, i.e., 50 community scenarios under per-community features. The community energy profiles of all communities were aggregated to return 15 profiles for the entire city.
Modelling system deployment pathways
From a government planning perspective, we model three system deployment pathways that coordinate households into energy community, including household distributed programs, household centralized programs, and community centralized programs37. These potential community collaboration programs govern energy allocation and trades to ensure effective and equitable energy flows among households in their local communities. This requires that energy models include energy dispatch strategies and pricing mechanisms under different system deployment pathways. To distinguish the roles of households clearly in energy models, we defined intra-household layers and intra-community layers based on control logic relationships rather than physical relationships. The intra-household layers refer to energy interactions between modules within the system, while the intra-community layers refer to energy interactions between different systems. Consequently, energy dispatch strategies are present in both layers, whereas pricing mechanisms are exclusive to the intra-community layers.
We detail energy dispatch strategies using a two-stage logical framework which defines energy flow structure in intra-household layers and energy sharing structures in intra-community layers from the rational household perspective (Supplementary Note 5). In this study, our analysis of multiple deployment pathways focuses on intra-community layers that reflect collaborative solutions in local community energy markets. Thus, each deployment pathway features the same energy flow structure but differs in energy sharing structure. The energy flow structure in three deployment pathways identifies twelve operating states of the intra-household layer to manage energy flow hierarchy in system modules under varying module capacity and power states, with the goal of eliminating module interference. To exploit the storage advantages of the green hydrogen systems in large storage capacity and long timescale, we managed the battery module as a short-term energy storage option and the hydrogen module as a long-term energy storage option. Energy sharing structures are used to govern the energy trade methods between prosumers and consumers in the intra-community layer. These structures distinguish trade entities based on ownership of energy systems and trade paths based on internal energy sharing pricing. Household distributed programs describe prosumer-owned deployment pathways which specify each prosumer as an energy trade entity with a three-layer trade path. Household centralized programs describe prosumer-grouped deployment pathways which combine all prosumers as energy trade entities with a two-layer trade path. Community centralized programs describe community-owned deployment pathways that coordinate all households as energy trade entities with a one-layer trade path.
The pricing model is essential for determining energy sharing pricing between buyers and sellers in intra-community layers to facilitate energy sharing. We utilize the terms "buyers" and "sellers" to delineate household roles in the pricing model, as prosumers may adapt their behaviors to function as either sellers or buyers based on their net power profiles. We assume that all sellers have equal privilege and are equally influential participants within the energy communities, implying that all sellers should collectively determine energy sharing prices. In accordance with the basic principles of economics, disparities in demand response can result in price fluctuations. This signifies that the market price for buyers favors the grid output price, whereas the market price for sellers favors the grid input price. Thus, we use a dynamic internal pricing model that uses local feed-in tariffs to define trade prices using the production and demand ratio (SDR) of shared energy and price boundaries. Further details can be found in ref.38. To maintain consistency with the time resolution of the energy model, the internal prices of the pricing model are adjusted in tandem with the hourly SDR within the boundary limitations. Moreover, the predicted time horizon of the pricing model is set to one hour ahead, with local electricity prices acting as the benchmark prices39,40. In this study, electricity sales prices were assumed to be 70% of electricity purchase prices. Notably, we disregarded the energy price fluctuation and local government price incentives during the operating cycle of the energy model.
Outline of green hydrogen systems in urban communities
As in refs.41,42, energy systems, long-term climate resilience systems designed for urban communities, must be considered cross-seasonal and large-capacity energy storage due to the high homogeneity of households in both energy production and demand. Compared to conventional battery storage systems43, green hydrogen systems provide higher storage and operation flexibility, which can schedule multi-storage resources in the most efficient way to address uncertainty, especially during extreme energy events. Specifically, green hydrogen systems can achieve synergistic benefits from both the superior efficiency of battery storage and the high energy density and low leakage rate of hydrogen storage. From a technical perspective, the design of green hydrogen systems for urban communities must prioritize safety and portability. Safety is paramount, given the risk of hydrogen leakage and spontaneous combustion44. Portability calls for compact equipment and installation methods that are both straightforward and minimally invasive. To meet these criteria, we have chosen photovoltaic panels as the energy production modules and paired them with lithium-ion battery packages and Liquid Organic Hydrogen Carrier (LOHC) equipment for energy storage modules (Supplementary Note 6).
Modelling design optimization for community green hydrogen systems
Design optimization in energy models involves sizing optimal system components of green hydrogen systems to inform evaluation results for metrics of interest. As discussed earlier, multi-source uncertainties should be incorporated to ensure the robustness and stability of evaluation outcomes. Scenario pools are utilized as part of stochastic optimization in design optimization, capturing compound impacts of uncertainty factors. These scenario pools comprise a set of community energy scenarios with occurrence probabilities. Once the energy input scenarios are established, design optimization maps decision space variables onto the objective space through cyclic simulation under inherent system constraints. Specifically, decision variables encompass the rated capacity of the battery package, along with the rated power and capacity of both the fuel cell and electrolyzer. Meanwhile, the objective variables consist of system affordability and independence, quantified as annual life cycle costs (LCC) and grid interaction level (GIL), respectively. In this study, the multi-objective particle swarm algorithm was used to calculate the design optimization part, as it has demonstrated effectiveness in handling multi-scenarios and multi-objectives simultaneously. To enhance the efficiency of the search process and prevent convergence to local minima, we incorporated mechanisms like self-adaptive adjustment of inertia weights45 and wavelet mutation46.
Formulating objective functions
The annual LCC represents the equivalent annual cost of green hydrogen systems over their lifespan in urban communities, including equivalent investment costs, equivalent operation and maintenance costs, carbon tax, and community trade costs, as formulated in Eq. (1). Notably, it is imperative to note that all the techno-economic data deployed in this study are derivatives of prevailing current social and technical conditions. The price uncertainties linked with technological upgrades and capital inflation are not considered in this study. The details of objective functions are shown in Supplementary Note 7, including each term of LCC and techno-economic data of green hydrogen systems.
$$LCC=\sum _{\forall s\in {N}_{s}}{\rho }_{s}N\left({C}_{s,inv}+{C}_{s,om}\right)CRF+\sum _{\forall n\in N}\sum _{\forall s\in {N}_{s}}{\rho }_{s}\sum _{\forall t\in T}\left({C}_{n,s,t,tax}+{C}_{n,s,t,community}\right),\forall n\in N,\forall s\in {N}_{s},\forall t\in T$$
1
where, n, s, and t denote the system number, expected community energy scenarios and time series; N, Ns, and T denote the total system number (prosumer number in HDP, 1 at HCP and CCP), total scenario number, and total simulation time; ρs is the probability of scenario s; CRF denotes the capital recovery factor; Cinv and Com denote the initial investment cost and operation and maintenance cost per system; Cn,s,t,tax, and Cn,s,t,community denote the carbon tax, and community interaction cost at time t in scenario s for system n.
GIL quantifies the level of independence of distributed energy systems by evaluating interactions with the power grid. The objective for urban residential communities is to minimize their interactions with the grid. This approach not only maximizes the role and benefits of communities but also helps to improve the stability of the urban power grid. This study evaluates GIL based on the amount of electricity imported from and exported to the grid, as shown in Eq. (2).
\(GIL=\sum _{\forall n\in N}\sum _{\forall s\in {N}_{s}}{\rho }_{s}\sum _{\forall t\in T}\left({P}_{n,s,t,gridsell}+{P}_{n,s,t,gridbuy}\right),\) \(\forall n\in N,\forall s\in {N}_{s},\forall t\in T\) (2)
where, Pn,s,t,gridsell t and Pn,s,t,gridbuy denote the amount of electricity purchased from the grid and the amount of electricity sold to the grid at time t in scenario s.