Published name
Please provide any comments in relation to these provisions.
We applaud the inclusion of information on the electricity grid to which the facility is connected
Please provide any comments in relation to these provisions.
We find that
1) the emissions intensity assigned to grid-based electricity using both the PGO location-based and market-based methods can lead to significant differences in the certified and actual embedded emissions intensity of hydrogen produced by grid-connected electrolysers.
2) when using RECs to certify hydrogen emissions intensity, we find that geographic correlation coupled with less stringent temporal correlation is an effective policy measure to ensure low-emissions hydrogen without driving up cost.
3) under current requirements around the use of renewable nergy certificates, (i.e. without geographic correlation and only yearly temporal correlation), we show that it is possible for hydrogen that has been certified as zero-emissions under the PGO scheme to generate actual grid electricity emissions exceeding 20 kgCO2e/kgH2, similar to that of hydrogen produced via coal gasification.
Please see the submitted documentation for details.
Do you have any additional feedback in relation to PGO rules (Parts 1-3 of the exposure draft) discussed in the consultation paper that you haven’t already provided?
Would you like to upload a document?
Upload a submission
Submission to the Consultation on
the Future Made in Australia
(Guarantee of Origin) Rules 2025
Associate Professor Fiona J Beck, Mr Chengzhe (Kevan) Li, Dr Anna Skobeleva, Dr Michael
Thomas
School of Engineering
The Australian National University
Canberra ACT 2600 Australia www.anu.edu.au
CRICOS Provider No. 00120C
The Australian National University 1
Contents
Contents ........................................................................................................................................................ 2
Executive Summary................................................................................................................................. 3
Details of Analysis ................................................................................................................................... 5
Energy system modelling framework ................................................................................... 6
Results and Discussion.................................................................................................................. 9
References ............................................................................................................................................. 13
The Australian National University 2
Executive Summary
This submission provides an analysis of the accuracy of emissions intensity of hydrogen certified under the Product Guarantee of Origin (PGO) Scheme, and the implications of different policy settings on the cost and emissions intensity of hydrogen produced by electrolysers connected to the grid in different states making up the National Energy Market (NEM). Specifically, we investigate the importance of: geographic correlation - which is the requirement that renewable energy certificates (RECs) used to certify the emissions intensity of hydrogen be generated on the same regional grid as the electricity sourced to power the electrolyser, and temporal correlation - which is the requirement that RECs used to certify the emissions intensity of hydrogen should be generated within the same specified time period as the electricity sourced to power the electrolyser.
Our analysis employs energy system modelling to quantify the difference between the actual emissions and those assigned by the PGO certificate under current rules. We explore how and why proposed emissions accounting frameworks lead to inaccuracies and identify effective policy settings that can keep both emissions and costs low.
Assessing accuracy
We find that the emissions intensity assigned to grid-based electricity using both the PGO location-based and market-based methods can lead to significant differences in the certified and actual embedded emissions intensity of hydrogen produced by grid-connected electrolysers.
We analyse systems consisting of dedicated renewable energy (RE) capacity and electrolysers which can buy and sell electricity from the Australian NEM to minimise operational costs (See
Figure 1 and associated description).
The PGO market-based method tends to overestimate hydrogen emissions intensities when large amounts of electricity are sourced from the grid (by up to 274% in some cases), and overestimates emissions offsets when selling electricity by surplus generation of RECs (up to 245% in some cases). This is because the market-based method uses an emissions factor (RMF) and applicable renewable power percentage (ARPP) that are estimated at the national level[1]. The discrepancies are larger for electrolysers connected to grids with higher RE penetration and lower emissions, disadvantaging producers in these states. We suggest that use of more granular, state-level factors would improve accuracy (see Figure 2 and associated discussion for more information).
The PGO location-based method overestimates hydrogen emissions intensities for grids with relatively high average electricity emissions factors (up to 63% in Victoria) and underestimates them for grids with relatively low hydrogen emissions intensities (up to 42% in South Australia).
This is because the emissions intensity of the grid electricity used by the electrolyser is determined by the ‘marginal’ generator which would respond to additional demand on the local grid, which can differ significantly from the state-level average electricity emissions factor used in the location-based method (see Figure 2 and associated discussion for more information).
While overestimation is consistent with the principle of conservativeness recommended for embedded emissions accounting[2], underestimation of embedded emissions contravenes the conservativeness principle and undermines the scheme’s emissions-reduction goals.
Policy measures to ensure the production of low-emissions hydrogen when relying on RECs
When using RECs to certify hydrogen emissions intensity (using the PGO market-based method), we find that geographic correlation coupled with less stringent temporal correlation is an effective policy measure to ensure low-emissions hydrogen without driving up cost.
However, while yearly temporal correlation is sufficient to ensure that the certified emissions intensities under the PGO market-based method are all below zero, actual emissions intensities
The Australian National University 3
can still be above zero in some cases: as much as 3.6 kgCO2e/kgH2 for the system in NSW, similar to those of natural gas based blue hydrogen systems in NSW [3]. Ensuring that actual emissions from the production of hydrogen are below zero for all states will require temporal correlation intervals somewhere between monthly and yearly durations and will likely change with location and over time, as it is determined by the interplay between regional electricity prices, local renewable energy resources and RE penetration in the local state grid (see Figure 3 and associated discussion for more information).
Discrepancies in certified emissions intensity of hydrogen under the PGO scheme when relying on
RECs
The Future Made in Australia (Guarantee of Origin) Act 2024 imposes a temporal correlation requirement of a minimum of 1 year and 1 hour and a maximum of 2 years under Clauses 50 and
51, but no geographical correlation.
Under these requirements, we show that it is possible for hydrogen that has been certified as zero-emissions under the PGO scheme to generate actual grid electricity emissions exceeding 20 kgCO2e/kgH2, similar to that of hydrogen produced via coal gasification[4]. We also show that while the certified emissions are consistent across states, actual emissions strongly depend on the location of the RE capacity generating the RECs (see Figure 4 and associated discussion for more information).
Since Australia aims to produce relatively large amounts of hydrogen for industry and export, this could lead to large increase in domestic emissions due to increased electricity demand, which would not be accurately attributed.
Further details can be found in the analysis presented below, and in a preprint of this analysis accessible here. We welcome further discussion of these findings and are currently undertaking further work to ascertain optimal temporal correlation settings.
The Australian National University 4
Details of Analysis
Certification and additionality
The certification of hydrogen produced via electrolysis is relatively straightforward when the power is sourced directly from dedicated, off-grid renewable energy (RE), such as wind or solar.
Certification becomes more complex for grid-connected electrolyser systems, as the emissions from grid electricity provided by a dynamic mix of generators vary with time and location[5]. The challenge for certification schemes is to balance the requirement to accurately account for electricity emissions without imposing excessive regulatory burden or driving up production costs[2].
To address this issue, the European Union (EU) has established a methodology in the Renewable
Energy Directive II[6] (RED II) and associated Delegated Regulations[7], [8] to certify electricity emissions associated with grid-connected hydrogen production under its CertifHy scheme[9].
Grid electricity is defined as “fully renewable”, and hence zero emissions, if it meets three criteria:
1. RE is supplied by renewable capacity commissioned no more than 36 months before the
start of hydrogen production (additionality);
2. RE should be generated and consumed by hydrogen production within a certain time
interval (temporal correlation or simultaneity);
3. and RE should be generated within the same interconnected bidding zone as the
electrolyser (geographic correlation).
If these conditions are not met, then the electricity emissions are calculated using the average emissions intensity of the electricity grid[7].
Among the RED II criteria, temporal correlation has proven to be an important and controversial policy lever[10], [11]. While strict time-matching may lower the emissions of hydrogen by ensuring that the system can only use grid electricity when RE is being generated, it also impacts the profitability by limiting the operation of the electrolyser and requiring electricity or hydrogen storage to firm supply[12], [13]. Recognizing the trade-off between economic competitiveness and minimizing emissions, jurisdictions have imposed less restrictive temporal correlation in the short term to enable grid-connected hydrogen production. For example, the EU’s CertifHy has initially adopted temporal correlation over calendar months, and plans to gradually shift to hourly matching by 2030[8]. Likewise, the Inflation Reduction Act (IRA) in the United States starts with the requirement of annual temporal correlation, gradually shifting to an hourly basis to assess qualification for the emissions intensity dependent tiers of the hydrogen production tax credit[14].
Additionality criteria under the Future Made in Australia (Guarantee of Origin) Act 2024
We recognise that the Australian PGO scheme does not seek to define emission thresholds for
“low-emissions” or “renewable” hydrogen, unlike CertifHy. Instead, the emission intensity for each batch of hydrogen is certified using a modular approach designed to meet the requirements of the International Partnership for Hydrogen and Fuel Cells in the Economy’s methodology[15]. We applaud this approach and recognise that accounting for emissions in each part of the supply chain allows interoperability with different international schemes with different definitions of “low-emissions hydrogen” [16].
The PGO scheme proposes two methods to estimate the emissions of electricity purchased from the grid, using methodologies defined under the National Greenhouse and Energy Reporting
(NGER) scheme[17]. The location-based method applies government defined annual average state emissions factors to electricity sourced in different states, while the market-based method allows consumers of grid electricity to offset their emissions by purchasing and surrendering renewable energy certificates (RECs, including Renewable Energy Guarantee of Origin
The Australian National University 5
certificates (REGOs), large-scale generation certificates, and other related certificates) which are generated for each megawatt hour of eligible renewable electricity sold to the grid.
We recognise that information included with each PGO certificate allows flexibility for the holders of registered profiles to choose to meet additionality as well as geographic and temporal correlation requirements in order to be eligible under different certification schemes, or qualify for incentives under green industrial policies[2]. We applaud that the newly developed REGOs will include the commissioning date and location of the power station, and an hourly time stamp recording when the electricity is generated[18]. We also support the fact that the Future Made in
Australia (Guarantee of Origin) Rules 2025 also require PGO certificates to state (Clause
20(1)(f)(g)) (i) whether the [renewable] electricity is from the same grid to which the facility is connected; and (ii) the time period within which the electricity was produced and subsequently used for the production pathway.
However, the Future Made in Australia (Guarantee of Origin) Act 2025 only specifies the following requirements for the use of RECs when calculating the emissions intensity of hydrogen:
Additionality
RECs can be generated by eligible RE systems built before the electrolyser, as long as they began generating renewable electricity after 1997, and do not have to be surrendered as part of
Renewable Energy Target[19], [20], [21].
Geographic Correlation
There is no requirement for geographic correlation. REGOs surrendered as part of a PGO certificate do not need to be sourced from the same grid as the electrolyser.
Temporal correlation
RECs used to claim RE for the production of hydrogen under the PGO must have been generated within one year of the start of production of the batch (Division 4, Subdivision A, Clause 51(8)). In addition, the amount of time taken to produce the batch (Clause 50(1)(h)) must be between 1 hour and 1 year (Clause 50(3)).
Taken together, these rules in effect impose temporal correlation requirement of a minimum 1 year and 1 hour and a maximum of 2 years.
In the following analysis, we assess the implications of these policy settings on the costs and emissions of hydrogen production.
Energy system modelling framework
We establish an energy system model of a grid-connected hydrogen production system, consisting of dedicated wind and solar, an electrolyser, two compressors, hydrogen storage facilities and pipelines, as shown schematically in Figure 1. Emissions primarily stem from upstream scope 2 GHG emissions attributable to electricity imports from the grid. Other emissions such as those from the import of potable water and hydrogen leakage 1 are neglected in this study.
The system is operated with hourly time resolution and delivers a constant supply of hydrogen to industrial end-users throughout the year[23]. 2 Local wind and solar electricity generation are calculated using a standard model developed by NREL[24] based on historical weather data for each local area[25].
To accurately track the grid electricity emissions, we assume that the system can only access
RECs generated by the dedicated RE capacity that makes up the system. This is equivalent to assuming that hydrogen production systems will need to source a dedicated supply of RE and not
1
Noting that hydrogen has a global warming potential of 11[22]
2
Note that the constant hourly supply is to ensure the hydrogen supply reliability required by refineries, ammonia plants, liquefaction facilities, and other industrial operations.
The Australian National University 6
just rely on the RECs market, which we consider reasonable given the relatively tight market for
RECs in Australia[26].
The system can use power directly from the dedicated RE capacity to run the electrolyser and compressors or can choose to buy or sell electricity from the Australian NEM, modelled using historical data on spot prices and emission factors for the five interconnected states[27].
Renewable energy capacity can be directly connected to the electrolyser when co-located, bypassing the grid, or connected through the grid. Purchasing electricity from the grid incurs a
“transmission use of system fee” (TUOS) representing the costs required to maintain the electricity grid set by grid operators[28], [29].
Figure 1. System boundaries for emissions accounting of hybrid grid-connected hydrogen production under the Product Guarantee of Origin scheme. The arrows indicate energy commodity flows and their directions: orange arrows for electricity, and blue arrows for hydrogen. The electrolyser converts electricity into hydrogen with an efficiency of 70%[30].
The linear optimization model is used to determine the capacity of the different components as well as the interaction with the grid to ensure the lowest hydrogen production cost, calculated as the levelized cost (LCOH). The capacity of the system is constrained so as not to perturb the spot price or emission intensity of the grid. All details of the costs can be found here[ref]
Policy settings are incorporated into the model by introducing additional constraints. Geographic correlation can be imposed by requiring that the dedicated RE capacity is connected in the same state as the electrolyser (directly connected to the electrolyser or connected through the grid).
Temporal correlation requires that the total RE sold to the grid by the system must be larger or equal to the total grid-electricity purchased to run the electrolyser within the specified time interval[12].
Calculating emissions from grid connected electricity
The emissions associated with the electricity used for hydrogen production are calculated using three different methods. Two of these methods are defined by the Australian PGO scheme[31], differentiated by whether they include the buying and selling of RECs, known as the market- based method, or use a state based emissions factor defined by the Australian government[17], known as the location-based method.
The Australian National University 7
For the PGO market-based method, electricity emissions ( ) are calculated using the applicable renewable power percentage (), which gives the fraction of the electricity on the grid that is renewable, and the residual mix factor (), which is the estimated average emission factor for the grid that excludes all the electricity generation claimed by RECs[1]. The emissions of electricity used to produce hydrogen is then
= × (1 − ) − × 1000 × , for the general case where is the quantity of grid electricity used and is the quantity of RECs surrendered.
We apply the market-based method in our model from time 1 to 2 as
2
= ( () × (1 − ) − ()) × ,
=1
where () and () are the quantity of electricity bought and sold to the grid at each time step which is then summed over time 1 to 2 . We assume that each MWh of RE sold to the grid
generates one REC such that, = ()/1000, which can be surrendered to cover emissions from using grid electricity. We assume that excess RECs will be traded to other facilities, enabling them to claim reductions in their electricity emissions under the market-based method.
Consequently, one unit of eligible RECs can represent an emissions reduction on the grid of 0.81 tCO2e/MWh.
The location-based method ( ), uses the annual average emissions intensity of the grid in each state ( ), defined by the Australian government as part of the National Greenhouse Gas
Accounting Factors[17]. These emissions factors are applied to the net energy bought or sold to the grid in each time step, summed over time 1 to 2 as
2
= ( () − ()) ×
=1
The third method represents the best estimate for the actual emissions from grid-based electricity as it assigns time dependent marginal emissions factors (MEFs) from historical NEM data for each of the five interconnected states to electricity bought from and sold to the grid in each time step[32], [33], [34]. The marginal emissions factor represents the emissions intensity of the ‘marginal’ generator which would be dispatched to meet the additional electricity demand from the electrolyser. MEFs can also quantify the grid emissions avoided by displacing electricity from the marginal generator with RE sold into the grid. We calculate the actual emissions intensity of electricity used by the electrolyser using the time dependent marginal emission factor in each Australian state, (), applied the net electricity consumption within each time step, before summing over time 1 to 2
2
= (( () − ()) × ()).
=1
The electricity emissions intensity can be negative, representing the case when the systems contribute excess RE to the grid and displace grid emissions.
Comparing the emissions calculated using the PGO methods with the actual emissions calculated using the MEF-based method allows us to quantify how well the proposed certification methodology captures the actual emissions.
The Australian National University 8
Results and Discussion
Assessing the accuracy of the embedded emissions accounting frameworks in the GO scheme
We calculate the emission intensity of hydrogen produced by systems that maintain a constant electrolyser size (equal to that of an optimized off-grid system, around 10-50MW) but have varying levels of dedicated RE capacity. The RE factor is defined as the ratio between RE generators (sum of wind and solar) and electrolyser capacity: an RE factor of 0 is a fully grid- reliant system, while an RE factor of 2 indicates the system has RE capacity two times larger than the electrolyser. The dedicated RE capacity is co-located and directly connected to the electrolyser ensuring geographic correlation and the systems can interact with the grid flexibly without any temporal correlation. Further details are provided in the Figure 2 caption.
(a) Gladstone (QLD) (b) Newcastle (NSW) (c) Geelong (VIC)
1.2 1.2 1.2
Emission intensity (tCO2e/MWhH2)
Emission intensity (tCO2e/MWhH2)
Emission intensity (tCO2e/MWhH2)
1.0 1.0 1.0
0.8 0.8 0.8
0.6 0.6 0.6
0.4 0.4 0.4
0.2 0.2 0.2
0.0 0.0 0.0
−0.2 −0.2 −0.2
−0.4 −0.4 −0.4
−0.6 −0.6 −0.6
−0.8 −0.8 −0.8
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 0.0 0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4 2.7
RE factor RE factor RE factor
(d) Burnie (TAS) (e) Spencer Gulf (SA)
1.0 1.0
Emission intensity (tCO2e/MWhH2)
Emission intensity (tCO2e/MWhH2)
0.8 0.8
0.6 0.6
0.4 0.4
EI_RECs
0.2 0.2 EI_Market
0.0 0.0 EI_MEF
−0.2 −0.2 EI_L
−0.4 −0.4
−0.6 −0.6
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 0.0 0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4 2.7 3.0
RE factor RE factor
Figure 2 (a)-(e) Emission intensities of hydrogen under various RE factors, defined as the ratio between RE generators and electrolyser capacity in five states of the NEM, calculated using three different methods for electricity emissions. Actual emissions are tracked using time and state dependent marginal emission factors to net grid electricity consumption (EI_MEF, orange curve). Certified emissions are calculated using the market-based (EI_Market, blue curve) and location-based (EI_L, green curve) adopted in the PGO scheme.
Negative emissions intensity values indicate that more renewable electricity is sold to the grid than is bought per unit hydrogen produced, resulting in a reduction of grid emissions (for location and MEF methods).
EI_RECs represents the RECs that can be generated for each MWh hydrogen produced under the market- based method, when excess RE is generated and subsequently reduces grid electricity emissions (green bars). The higher heating value of hydrogen is used throughout the manuscript (39.4 kWh/kg H2[35]).
The key result from this analysis is that actual emissions intensities of hydrogen, represented by
EI_MEF (red line), are not always well represent by the PGO emission accounting methods.
The emissions intensity calculated using the PGO market-based method (blue line) is the same for all systems. It overestimates the hydrogen emissions intensity for fully grid reliant systems
(RE Factor = 0) and overestimates the emissions avoided by systems generating excess RE that it sells to the grid (maximum RE Factor) quantified by excess RECs (shown as green bars as negative emission intensities). The largest discrepancies are observed for TAS, with emissions intensities overestimated by up to 274%, and RECs overestimated by up to 245%. Deviations reduce as less grid-electricity is bought and sold.
The Australian National University 9
We attribute these discrepancies to the adoption of nation-wide factors in the PGO market-based method. The proportion of RE on the state grids can be much higher (93% and 71% in 2023 for
TAS and SA respectively) than the national level applicable renewable power percentage (ARPP
= 18.72%). Likewise, the national level residual mix factor (RMF) is calculated by adjusting the average emission factor of all generation on the NEM, South-West Interconnected System, and the Darwin and Katherine Interconnected System[1], by excluding all RE that has generated RECs.
This neglects the energy mix of the state grids on the NEM, resulting in larger discrepancies for states with higher RE penetration.
The PGO location-based method overestimates the emissions intensity in QLD, NSW, and VIC by up to 63% but underestimates emissions for SA by 42% (RE Factor = 0). This is because MEFs for states dominated by coal (QLD, VIC, NSW) tend to be lower than the average state electricity emission factors used in the location-based method, as over 49% of the time the marginal generator will be a lower emissions technology like gas or renewables. Conversely, SA tends to have MEFs higher than average state electricity emission factors as coal is more likely to be the marginal generator when importing electricity.
Impact of temporal correlation
Next, we optimise the configuration (i.e. the capacity of different components) and operation of hybrid systems to achieve lowest cost hydrogen production under different constraints. We compare off-grid systems, fully flexible systems with no constraints, and systems with different temporal correlation requirements, while ensuring geographic correlation as above. Further details are provided in the Figure 3 caption.
Burnie (TAS) Newcastle (NSW) Geelong (VIC)
45 45
7 7 40 7
40 40
6 6 35
6
Emission intensity (kgCO2e/kg H2)
Emission intensity (kgCO2e/kg H2)
Emission intensity (kgCO2e/kg H2)
35 30 35
5
5 30 25 5 30
4
4 25 20 4 25
LCOH ($/kg H2)
LCOH ($/kg H2)
LCOH ($/kg H2)
3
20 15 20
3 2 10 3
15 15
2 1 5 2
10 10
0 0
1 5 1 5
−1 −5
0 0 −2 −10 0 0
−5 −15 −5
−1 −3 −1
−10 −20 −10
−2 −4 −25 −2
−15 −5 −15
−3 −30 −3
−20 −6 −35 −20
id
y
le
id
y
le
id
y
ly
ly
≤0
y
ly
ly
≤0
ly
y
le
y
ly
≤0
rl
rl
rl
nl
nl
nl
ib
ib
ib
ai
ai
ai
th
th
th
gr
gr
gr
EF
EF
EF
ea
ea
ea
o
o
o
D
ex
D
ex
D
ex
on
on
on
ff-
ff-
ff-
ll
ll
ll
Y
Y
Y
_M
_M
_M
Fl
Fl
Fl
Se
Se
Se
M
M
M
O
O
O
EI
EI
EI
Spencer Gulf (SA) Gladstone (QLD)
45 45
7 7
40 40
6 6
Emission intensity (kgCO2e/kg H2)
Emission intensity (kgCO2e/kg H2)
35 35
LCOH
5 30 5 30
25 25
EI_Market
4 4
LCOH ($/kg H2)
LCOH ($/kg H2)
20 20 EI_MEF
3 3
15 15 Grid electricity cost
2 2
10 10 Electrolyser CAPEX+OM
1 5 1 5
Hydrogen Storage CAPEX
0 0 0 0
−5 −5
PV CAPEX+OM
−1 −1
−10 −10 Wind CAPEX+OM
−2 −2
−15 −15
−3 −3
−20 −20
id
ly
y
le
id
le
y
ly
≤0
y
ly
y
ly
≤0
rl
rl
nl
nl
ib
ib
ai
ai
th
th
gr
gr
EF
EF
ea
ea
o
o
D
ex
D
ex
on
on
ff-
ff-
ll
ll
Y
Y
_M
_M
Fl
Fl
Se
Se
M
M
O
O
EI
EI
Figure 3 Total LCOH (white dots), average yearly actual (EI_MEF, yellow triangles) and certified emission intensities of hydrogen (EI_Market, red triangles) calculated using the PGO market-based method, optimized system configurations in five states across different under different constraints including temporal correlation from hourly (sell-only) to yearly. Data for off-grid systems are given for comparison.
EI_MEF≤0 systems are optimised to minimize cost while constraining actual emissions to be equal to or below zero. The emission intensity is given in kgCO2e/kgH2 to allow comparison with existing certification schemes. Bars show LCOH deconvoluted to components: grid electricity cost (sum of the total expenditure and profit from buying and selling electricity to and from the grid, which can be negative), capital expenditure (CAPEX) of hydrogen storage, and CAPEX and operation and maintenance costs (OM) of electrolyser, PV and wind generators. Negative emission intensity values for market-based method
(EI_Market) quantify grid electricity emissions offset by excess RECs generated for each kg hydrogen
The Australian National University 10
produced, while they represent the actual displaced electricity emissions for the MEF-based method
(EI_MEF).
Grid connection reduces costs for all states, with the lowest LCOH occurring for fully flexible scenarios with no additional constraints. There are two distinct trends for optimized system configurations and emission intensity of hydrogen: one that holds for TAS, NSW and VIC, and one that holds for SA and QLD.
Due to the interplay between spot prices and the quality of the RE resources, the cheapest way to produce electrolytic hydrogen in TAS, NSW and VIC is to minimize capital investment and source large amounts of power from the grid. For TAS and VIC, relatively low grid electricity prices make dedicated RE uncompetitive. While NSW has the highest yearly average electricity price in the
NEM, the local RE resources in Newcastle are less abundant (quantified by the RE capacity factor, leading to higher RE costs. When unconstrained by temporal correlation in the flexible scenario, the system builds less dedicated RE capacity in TAS and NSW compared to the off-grid scenario, and none in VIC, and invests in smaller hydrogen storage and electrolyser capacity. Critically, relying on grid power in the lowest cost flexible scenarios can result in very high actual hydrogen emission intensities, with EI_MEF =~25 kgCO2e/kgH2 for VIC and NSW, comparable to hydrogen made from coal gasification[4]. Once again, we see that the PGO market-based method overestimates the certified emissions intensity, particularly for TAS and VIC.
Once yearly temporal correlation is imposed in TAS, NSW and VIC, the system is forced to build RE capacity, and the actual emission intensity falls to below zero in TAS and slightly above zero in
NSW and VIC. However, the LCOH is only very slightly increased.
As the temporal correlation requirements become more restrictive, shifting from yearly to hourly time matching, the system increases the capacity of the electrolyser and storage, and the configuration progressively resembles an off-grid setup. Consequently, LCOH increases by 30%-
38%. Once again, the PGO market-based emission calculation overestimates potential emissions offsets by RECs, becoming more significant when the system sells larger amounts of electricity to the grid.
The trend is different for SA and QLD. All systems have actual emissions intensities close to or below zero, and the flexible and yearly temporal correlation scenarios all result in nearly identical system configurations, with similar and lowest cost LCOH. The cost reduction in these two scenarios compared to the off-grid system is mainly attributable to the systems’ ability to profit from selling surplus electricity to the grid, driven by favourable local RE resources and higher regional electricity prices. The system significantly oversizes the cheapest RE capacity and sells large amounts of electricity to subsidize the cost of hydrogen production. As the temporal correlation requirements become more restrictive, shifting from yearly to hourly time matching, the restrictions on importing electricity from the grid results in the establishment of configurations similar to the off-grid system with PV and larger storage components. The systems are prevented from subsidizing the cost of hydrogen by selling electricity, increasing the LCOH by
78 % in QLD, and 52% in SA.
One way to assess the efficacy of temporal correlation as a policy is to compare with systems optimized to minimize cost while constraining actual emissions to be equal to or below zero
(EI_MEF≤0). Figure 3 shows that such systems are most similar to those optimized under yearly temporal correlation, suggesting it is an efficient policy measure to keep emissions low, at low cost. However, it is important to note that while yearly temporal correlation is sufficient to ensure that the certified emissions intensities under the PGO market-based method are all below zero, actual emissions intensities are still 1.9 kgCO2e/kgH2 for the VIC case and 3.6 kgCO2e/kgH2 for the
NSW case. These emissions intensities are similar to those of natural gas based blue hydrogen systems [3].
The Australian National University 11
Impact of geographic correlation
So far, all hybrid systems have been geographically correlated, with a direct connection between co-located RE and hydrogen electrolysers. Next, we relax that constraint and compare scenarios in which the electrolyser and the dedicated RE generators are located in different states and interact only through grid-connection. In these scenarios, the RE is sold to the local grid where it is located, generating RECs and offsetting actual emissions in that state, while the electrolyser buys electricity from its local grid, generating actual emissions and using the RECs generated by the RE to offset emissions. We focus on an electrolyser located in Newcastle (NSW), which is the site of the announced Hunter Valley Hydrogen Hub, slated to include a 50 MW electrolyser powered by grid-connected RE through purchase and retirement of RECs[36]. We compare systems under the constraint of yearly temporal correlation with dedicated RE located in different states, with a fully grid-powered system and an off-grid system (Figure 4).
The LCOH can be reduced by as much as 21% when the RE is located in QLD due to the superior
RE resources and high electricity costs, as the net cost of grid electricity is determined by buying and selling electricity to and from two different grids and depends on the spot prices in both states.
Hydrogen production system
in Newcastle (NSW)
7.0
(a) LCOH
6.5
6.6
6.5 LCOH with direct line
LCOH with RECs cost
6.0
LCOH ($/kg H2)
5.5
5.1
5.0 5.2
4.8 4.8
4.6
4.5
4.5
4.6
4.0
3.8
3.5
40 (b) 38
35 EI_Market
30 EI_MEF 32
Emission intensity (kgCO2e/kg H2)
25 21
20
15
10
10 7.6
3.7
5
0
−5 −8.7 −8.7 −8.7 −8.7
−3.6
−10
−15
−20 −24
−25
−30
VIC QLD SA TAS NSW Fully-grid Off-grid
(Geelong) (Gladstone) (Spencer Gulf) (Burnie) (Newcastle) (Newcastle) (Newcastle)
RE generators' location
Figure 4 (a) LCOH and (b) emission intensity of hydrogen produced by a grid-connected electrolyser located in Newcastle NSW, with RECs sourced from grid-connected RE capacity located different states, as indicated along the x-axis. Systems are optimised under the constraint of yearly temporal correlation and are compared to off-grid and fully-grid connected systems in NSW. The green line indicates the estimated cost of purchasing RECs from a third party to offset emissions for the fully on-grid case, assuming REC costs of 20-60AUD[26]. Note that the LCOH for a system with RE located in NSW connected to the electrolyser through the grid is slightly higher than that calculated previously (see Fig. 3) for a system with direct connection between co-located RE and electrolyser in NSW (red dot) due to increased transmission use of system fees. The emission intensity of hydrogen is calculated by MEF-based method (yellow triangles) and for the PGO market-based (red triangles).
Results clearly show that yearly temporal correlation is sufficient to ensure that all certified emissions (EI_Market, red triangles) are below zero even without geographic correlation, as each unit of electricity used by the electrolyser is compensated for by a REC generated within the year.
However, the actual emissions (EI_MEF, yellow triangles) vary widely with the location of the RE as the emissions displaced by selling RE will depend on the makeup of its local grid and are not necessarily equal to emissions generated by using electricity in NSW. In the worst case, the emission intensity can reach 21 kgCO2e/kgH2, similar to that of hydrogen produced via coal
The Australian National University 12
gasification[4], as high emissions intensity electricity consumed in NSW is being offset by displacing relatively low emissions intensity electricity in TAS.
In all scenarios so far, we assume a strict additionality requirement: that dedicated RE is deployed concurrently to the electrolyser. The PGO scheme does not require such strict additionality: RECs can be generated by eligible RE systems built before the electrolyser, as long as they began generating renewable electricity after 1997, and do not have to be surrendered as part of
Renewable Energy Target[19], [20], [21]. The fully grid-powered electrolyser scenario in Figure 4 approximates this case, buying electricity from the NSW grid at the spot price and generating the corresponding emissions according to the state-based MEF. Under this scenario, the additional cost for RECs to offset the certified emissions under the market-based method would increase the LCOH by 16-44% assuming RECs costs of AUD 20-60[26], while estimating the actual emissions would require knowing the time and state dependent MEF for each REC.
References
[1] DCCEEW, “Australian National Greenhouse Accounts Factors,” 2024. [Online]. Available:
https://www.dcceew.gov.au/sites/default/files/documents/national-greenhouse-account-
factors-2024.pdf
[2] L. V. White, E. Aisbett, O. Pearce, and W. Cheng, “Principles for embedded emissions
accounting to support trade-related climate policy,” Climate Policy, pp. 1–17, May 2024, doi:
10.1080/14693062.2024.2356803.
[3] T. Longden, F. J. Beck, F. Jotzo, R. Andrews, and M. Prasad, “‘Clean’ hydrogen? – Comparing
the emissions and costs of fossil fuel versus renewable electricity based hydrogen,” Applied
Energy, vol. 306, p. 118145, Jan. 2022, doi: 10.1016/j.apenergy.2021.118145.
[4] L. Eric et al., “Comparison of Commercial, State-of-the-Art, Fossil-Based Hydrogen
Production Technologies,” National Energy Technology Laboratory (NETL), 2022. [Online].
Available: https://doi.org/10.2172/1862910
[5] L. Engstam, L. Janke, C. Sundberg, and Å. Nordberg, “Grid-supported electrolytic hydrogen
production: Cost and climate impact using dynamic emission factors,” Energy Conversion
and Management, vol. 293, p. 117458, Oct. 2023, doi: 10.1016/j.enconman.2023.117458.
[6] European Commission, “Directive (EU) 2018/2001 of the European Parliament and of the
Council of 11 December 2018 on the Promotion of the Use of Energy from Renewable
Sources (Recast) (Text with EEA Relevance),” 2024. Accessed: Nov. 18, 2024. [Online].
Available: https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:02018L2001-
20240716
[7] “Commission Delegated Regulation (EU) 2023/1185 of 10 February 2023 supplementing
Directive (EU) 2018/2001 of the European Parliament and of the Council by establishing a
minimum threshold for greenhouse gas emissions savings of recycled carbon fuels and by
specifying a methodology for assessing greenhouse gas emissions savings from renewable
liquid and gaseous transport fuels of non-biological origin and from recycled carbon fuels,”
Cambridge University Press, 2023. doi: 10.1017/CBO9780511610851.034.
[8] European Commission, “Commission Delegated Regulation (EU) 2023/1184 of 10 February
2023 supplementing Directive (EU) 2018/2001 of the European Parliament and of the
Council by establishing a Union methodology setting out detailed rules for the production of
renewable liquid and gaseous transport fuels of non-biological origin,” 2023. Accessed: Oct.
09, 2024. [Online]. Available: https://eur-lex.europa.eu/legal-
content/EN/TXT/PDF/?uri=CELEX:02023R1184-20240610
[9] CERTIFHY, “CERTIFHY EU RFNBO VOLUNTARY SCHEME.” Accessed: Dec. 14, 2024. [Online].
Available: https://www.certifhy.eu/wp-
content/uploads/2024/09/CertifHy_Scheme_doc_Scheme_document.pdf
[10] M. A. Giovanniello, A. N. Cybulsky, T. Schittekatte, and D. S. Mallapragada, “The influence of
additionality and time-matching requirements on the emissions from grid-connected
The Australian National University 13
hydrogen production,” Nat Energy, vol. 9, no. 2, pp. 197–207, Jan. 2024, doi: 10.1038/s41560-
023-01435-0.
[11] J. Brandt et al., “Cost and competitiveness of green hydrogen and the effects of the
European Union regulatory framework,” Nat Energy, vol. 9, no. 6, pp. 703–713, May 2024, doi:
10.1038/s41560-024-01511-z.
[12] O. Ruhnau and J. Schiele, “Flexible green hydrogen: The effect of relaxing simultaneity
requirements on project design, economics, and power sector emissions,” Energy Policy, vol.
182, p. 113763, Nov. 2023, doi: 10.1016/j.enpol.2023.113763.
[13] D. Schlund and P. Theile, “Simultaneity of green energy and hydrogen production: Analysing
the dispatch of a grid-connected electrolyser,” Energy Policy, vol. 166, p. 113008, Jul. 2022,
doi: 10.1016/j.enpol.2022.113008.
[14] Internal Revenue Service Treasury, “Section 45V Credit for Production of Clean Hydrogen;
Section 48(a)(15) Election To Treat Clean Hydrogen Production Facilities as Energy
Property.” Accessed: Oct. 10, 2024. [Online]. Available:
https://www.govinfo.gov/content/pkg/FR-2023-12-26/pdf/2023-28359.pdf
[15] International Partnership for Hydrogen and Fuel Cells in the Economy, “Methodology for
Determining the Greenhouse Gas Emissions Associated with the Production of Hydrogen,”
2023. Accessed: Oct. 11, 2024. [Online]. Available:
https://www.iphe.net/_files/ugd/45185a_8f9608847cbe46c88c319a75bb85f436.pdf
[16] D. Goodwin, F. Gale, H. Lovell, K. Beasy, H. Murphy, and M. Schoen, “Sustainability
certification for renewable hydrogen: An international survey of energy professionals,”
Energy Policy, vol. 192, p. 114231, Sep. 2024, doi: 10.1016/j.enpol.2024.114231.
[17] National Greenhouse and Energy Reporting Scheme, “National Greenhouse and Energy
Reporting (Measurement) Determination 2008,” 2024. [Online]. Available:
https://www.legislation.gov.au/F2008L02309/latest/text
[18] DCCEEW, “Future Made in Australia (Guarantee of Origin) Act 2024,” 2024. [Online].
Available: https://www.legislation.gov.au/C2024A00121/asmade/2024-12-
10/text/original/pdf
[19] DCCEEW, “Renewable Energy (Electricity) Regulations 2001,” 2024. [Online]. Available:
https://www.legislation.gov.au/F2001B00053/2024-10-14/2024-10-14/text/original/pdf
[20] DCCEEW, “Renewable Energy (Electricity) Act 2000.” [Online]. Available:
https://www.legislation.gov.au/C2004A00767/2025-01-01/2025-01-01/text/original/pdf
[21] Clean Energy Regulator, “Renewable Energy Target.” Accessed: Mar. 24, 2025. [Online].
Available: https://cer.gov.au/schemes/renewable-energy-target
[22] Frazer-Nash Consultancy, “Fugitive hydrogen emissions in a future hydrogen economy,”
2022. [Online]. Available:
https://assets.publishing.service.gov.uk/media/624ec79cd3bf7f600d4055d1/fugitive-
hydrogen-emissions-future-hydrogen-economy.pdf
[23] J. M. Bracci, E. D. Sherwin, N. L. Boness, and A. R. Brandt, “A cost comparison of various
hourly-reliable and net-zero hydrogen production pathways in the United States,” Nat
Commun, vol. 14, no. 1, Art. no. 1, Nov. 2023, doi: 10.1038/s41467-023-43137-x.
[24] PySAM Version 2.2.0. (Sep. 18, 2024). C. National Renewable Energy Laboratory. Accessed:
Oct. 10, 2024. [Online]. Available: https://github.com/NREL/pysam
[25] M. M. Rienecker et al., “MERRA: NASA’s Modern-Era Retrospective Analysis for Research
and Applications,” Journal of Climate, vol. 24, no. 14, pp. 3624–3648, Jul. 2011, doi:
10.1175/JCLI-D-11-00015.1.
[26] Clean Energy Regulator, “Large-scale generation certificate (LGC) reported spot and
forward prices.” Accessed: Mar. 18, 2025. [Online]. Available:
The Australian National University 14
https://cer.gov.au/markets/reports-and-data/quarterly-carbon-market-reports/quarterly-
carbon-market-report-december-quarter-2024/large-scale-generation-certificates-lgcs
[27] AEMO, “NEMweb market data.” Accessed: Oct. 10, 2024. [Online]. Available:
https://aemo.com.au/energy-systems/electricity/national-electricity-market-nem/data-
nem/market-data-nemweb
[28] TasNetworks, “Prescribed Transmission Service Prices for 2024-25.” Accessed: Mar. 07,
2025. [Online]. Available:
https://www.tasnetworks.com.au/config/getattachment/96b170ef-8d44-43b9-97bf-
7ba6b2304bd0/2024-25-prescribed-transmission-service-prices.pdf
[29] N. Salmon and R. Bañares-Alcántara, “Impact of grid connectivity on cost and location of
green ammonia production: Australia as a case study,” Energy Environ. Sci., vol. 14, no. 12, pp.
6655–6671, 2021, doi: 10.1039/D1EE02582A.
[30] AEMO, “2023 costs and technical parameter review AEMO,” 2023. Accessed: Oct. 10, 2024.
[Online]. Available: https://aemo.com.au/-
/media/files/stakeholder_consultation/consultations/nem-consultations/2023/2024-
forecasting-assumptions-update-consultation-page/aurecon---2023-cost-and-technical-
parameters-review.pdf
[31] DCCEEW, “Guarantee of Origin - Emissions Accounting Approach paper,” 2023. [Online].
Available: https://storage.googleapis.com/files-au-climate/climate-
au/p/prj291cc9979281a4ffc59d8/public_assets/Guarantee%20of%20Origin%20-
%20Emissions%20Accounting%20Approach%20paper.pdf
[32] N. A. Ryan, J. X. Johnson, and G. A. Keoleian, “Comparative Assessment of Models and
Methods To Calculate Grid Electricity Emissions,” Environ. Sci. Technol., vol. 50, no. 17, pp.
8937–8953, Sep. 2016, doi: 10.1021/acs.est.5b05216.
[33] A. D. Hawkes, “Estimating marginal CO2 emissions rates for national electricity systems,”
Energy Policy, vol. 38, no. 10, pp. 5977–5987, Oct. 2010, doi: 10.1016/j.enpol.2010.05.053.
[34] A. G. N. Elenes, E. Williams, E. Hittinger, and N. S. Goteti, “How Well Do Emission Factors
Approximate Emission Changes from Electricity System Models?,” Environ. Sci. Technol., vol.
56, no. 20, pp. 14701–14712, Oct. 2022, doi: 10.1021/acs.est.2c02344.
[35] K. W. Harrison, R. Remick, G. D. Martin, and A. Hoskin, “Hydrogen Production: Fundamentals
and Case Study Summaries,” 2010. [Online]. Available:
https://www.nrel.gov/docs/fy10osti/47302.pdf
[36] CSIRO, “Hunter Valley Hydrogen Hub (Updated Project Description),” HyResource.
Accessed: Mar. 04, 2025. [Online]. Available: https://research.csiro.au/hyresource/hunter-
valley-hydrogen-hub/
The Australian National University 15
Upload supporting file
Assessing emission certification schemes for grid-connected
hydrogen in Australia
Chengzhe Li, Lee V. White, Reza Fazeli, Anna Skobeleva, Michael Thomas, Shuang Wang, Fiona J. Beck*
Abstract
Certification of low emissions products is crucial to ensure that they can attract a “green premium” in international markets but need to be designed so that they do not impose undue regulatory burden or drive-up production costs. In this work, we employ energy system modelling to evaluate how different policy choices affect the cost and certified emissions of hydrogen produced by grid-connected electrolysers in Australia. As a major exporter of energy and ores, Australia has ambitions to become a
‘renewable energy superpower’ by providing commodities with low embedded emissions. Our analysis shows that emissions are not always accurately accounted for under the currently proposed Australian methodologies. However, imposing geographic correlation and relatively relaxed temporal correlation requirements on the use of renewable energy certificates could be sufficient to ensure minimal emissions at low cost. Our findings provide evidence to assist in the design of robust and practical certification schemes.
Introduction
Internationally recognized certification of embedded emissions is critical to enable trade in so-called
“green” commodities like hydrogen and its derivatives where production processes and emissions can vary widely1,2, allowing customers to differentiate between otherwise identical products3,4. In the absence of a global price on carbon, such differentiation will be needed for low-emission products to compete in global markets5, and to facilitate the implementation of green industrial policies that incentivize decarbonization6, for example, hydrogen production tax credits in the Inflation Reduction
Act (IRA) of the United States scale inversely with certified emissions intensity7.
In general, the certification of hydrogen produced via electrolysis is relatively straightforward when the power is sourced directly from dedicated, off-grid renewable energy (RE), such as wind or solar.
Certification becomes more complex for grid-connected electrolyser systems, as the emissions from grid electricity provided by a dynamic mix of generators vary with time and location8. The challenge for certification schemes is to balance the requirement to accurately account for electricity emissions without imposing excessive regulatory burden, or driving up production costs6.
1
To address this issue, the European Union (EU) has established a methodology in the Renewable Energy
Directive II9 (RED II) and associated Delegated Regulations10,11 to certify electricity emissions associated with grid-connected hydrogen production under its CertifHy scheme12. Grid electricity is defined as “fully renewable”, and hence zero emissions, if it meets three criteria: RE is supplied by renewable capacity commissioned no more than 36 months before the start of hydrogen production
(additionality); RE should be generated and consumed by hydrogen production within a certain time interval (temporal correlation or simultaneity); and RE should be generated within the same interconnected bidding zone as the electrolyser (geographic correlation). If these conditions are not met, then the electricity emissions are calculated using the average emissions intensity of the electricity grid10.
Among the RED II criteria, temporal correlation has proven to be an important and controversial policy lever7,13. While strict time-matching may lower the emissions of hydrogen by ensuring that the system can only use grid electricity when RE is being generated, it also impacts the profitability by limiting the operation of the electrolyser and requiring electricity or hydrogen storage to firm supply14,15.
Recognizing the trade-off between economic competitiveness and minimizing emissions, jurisdictions have imposed less restrictive temporal correlation in the short term to enable grid-connected hydrogen production. For example, the EU’s CertifHy has initially adopted temporal correlation over calendar months, and plans to gradually shift to hourly matching by 203011. Likewise, the IRA in the US starts with the requirement of annual temporal correlation, gradually shifting to an hourly basis to assess qualification for the emissions intensity dependent tiers of the hydrogen production tax credit16.
The Australian government is developing a Product Guarantee of Origin (PGO) scheme to certify emissions embedded in products such as hydrogen17, to support its ambition to become “renewable energy superpower” by exporting low emission products and commodities18,19. A proposal for the methodology and emissions accounting framework that will underpin the PGO scheme has been released for consultation and the policy settings are under development20–22. Unlike the IRA or CertifHy, the Australian PGO scheme does not define emission thresholds for “low-emissions” or “renewable” hydrogen. Instead, the emission intensity for each batch of hydrogen is certified using a modular approach designed to meet the requirements of the International Partnership for Hydrogen and Fuel
Cells in the Economy’s methodology23. Accounting for emissions in each part of the supply chain allows interoperability with different international schemes with different definitions of “low-emissions hydrogen”4.
The PGO scheme proposes two methods to estimate the emissions of electricity purchased from the grid, using methodologies defined under the National Greenhouse and Energy Reporting (NGER) scheme24. The location-based method applies government defined annual electricity emissions factors
2
for each state in Australia, while the market-based method allows consumers of grid electricity to offset their emissions by purchasing and surrendering renewable energy certificates (RECs) which are generated for each megawatt hour of renewable electricity sold to the grid. Each REC will report the commissioning date and location of the power station, and a hourly time stamp recording when the electricity is generated17. While the information included with each REC grants the PGO scheme flexibility to develop rules to ensure additionality as well as geographic and temporal correlation, optimum settings have not yet been determined.
The Australian National Electricity Market (NEM) is both a wholesale market and a physical power system, providing electricity to 88% of Australia’s population living along the eastern seaboard25. The
NEM is one of the longest transmission networks in the world, connecting five state grids which each have unique mix of electricity generators25. South Australia has rapidly increased the penetration of renewables on its grid over the last 20 years, reaching 71% in 202326, and is on track to be powered by
100% renewable energy by 202727,28. Queensland has exceptional solar resources, but started the transition much later, reaching 28% in 2023 and aiming for 80% by 203529. The different levels of ambition present unique challenges for designing effective national policy settings for the PGO scheme.
While previous analyses have shown the importance of geographic and temporal correlation requirements on costs and emissions in European and US grids7,13–15,30, there has been no analysis of the impact of these policy choices on the competitiveness of hydrogen produced across states in
Australia.
In this context, our study aims to assess the implications of different policy settings in the proposed
PGO certification scheme on the cost and emissions intensity of hydrogen produced by grid-connected hydrogen production systems in Australia. We establish a linear optimization model to quantify the difference between the actual emissions and those assigned by the proposed PGO accounting methods for different cost-optimal system configurations. We explore how and why proposed emissions accounting frameworks under the PGO will lead to inaccuracies and identify effective and efficient policy settings that can keep both emissions and costs low. Lastly, we provide policy suggestions for the future development of the Australian PGO scheme.
Hybrid grid-connected hydrogen production
3
Figure 1. System boundaries for emissions accounting of hybrid grid-connected hydrogen production under the
Product Guarantee of Origin scheme. The arrows indicate energy commodity flows and their directions: orange arrows for electricity, and blue arrows for hydrogen.
The emission accounting system boundary of the Australian PGO scheme split into two main parts20: 1) the production boundary encompassing upstream and direct emissions; and 2) the post-production boundary including emissions from transport and storage of hydrogen. The PGO scheme requires that hydrogen producers define a batch period over which they are certifying their hydrogen production, ranging from hourly to yearly. During this batch period, the producer collects emissions and RECs information and reports it to the PGO scheme to label the certified emission intensity for the batch of hydrogen.
In this study, we evaluate the LCOH and associated actual and certified emissions for an optimized grid- connected hydrogen production system, powered by a mix of renewable and grid electricity as shown in Figure 1. Production emissions primarily stem from upstream scope 2 greenhouse gas (GHG) emissions attributable to electricity imports from the grid. Emissions from the import of potable water are included in the PGO scheme20 but are not included in our model for simplicity, as they would be identical whether hydrogen production is powered by off-grid RE or on-grid electricity. In the post- production boundary, we include scope 2 emissions associated with electricity used for compression, transport and storage. We neglect fugitive emissions of hydrogen which has a GWP of 1131 which could occur at various points during the supply chain.
4
We establish an energy system model of the grid-connected hydrogen production system, based on previously developed models of off-grid systems32 consisting of dedicated wind and solar, an electrolyser, two compressors, hydrogen storage facilities and pipeline (details in Method section). We assume that the system can only access RECs from RE capacity installed at the same time as the electrolyser, which can be directly connected, or accessed through the grid. This assumption allows us to accurately track emissions associated with RECs and reflects the rapid deployment of new capacity to meet the growing demand for RECs in Australia33.
The system is operated with hourly time resolution and delivers a constant supply of hydrogen to industrial end-users throughout the year34. The system can use power directly from the dedicated RE capacity or can choose to buy or sell electricity from the Australian NEM, modelled using historical data on spot prices and emission factors for the five interconnected states35. Purchasing electricity from the grid incurs a “transmission use of system fee” representing the costs required to maintain the electricity grid set by grid operators36,37. A linear optimization model is used to determine the capacity of the different components as well as the interaction with the grid required to ensure the lowest hydrogen production cost, calculated as the levelized cost (LCOH). The capacity of the system is constrained so as not to perturb the spot price or emission intensity of the grid.
We leverage our model to assign time and state dependent emission factors to net grid electricity consumption in each time step. Average emissions factors (AEFs) are computed by averaging emissions generated by all electricity generated across each of the five states, and marginal emissions factors
(MEFs) represent the emissions of the marginal generator, which is the generator that would be dispatched to meet additional electricity demand. MEFs best represent the actual electricity emissions caused by the additional demand on the grid from our system38–40 and quantify the grid emissions avoided by displacing electricity from the marginal generator with RE sold into the grid. Comparing the emissions calculated using the PGO methods described previously with the MEF-based method allows us to quantify how well the proposed certification methodology captures the actual emissions.
Policy settings can be incorporated into the model by introducing additional constraints. Geographic correlation can be imposed by requiring the dedicated RE capacity to be in the same state as the electrolyser. Temporal correlation requires that the total electricity sold to the grid must be larger than or equal to the total purchased electricity within the time interval14.
Evaluating certified emission intensity of hydrogen
34
Note that the constant hourly supply is to ensure the hydrogen supply reliability required by refineries, ammonia plants, liquefaction facilities, and other industrial operations.
5
(a) Gladstone (QLD) (b) Newcastle (NSW) (c) Geelong (VIC)
1.2 1.2 1.2
Emission intensity (tCO2e/MWhH2)
Emission intensity (tCO2e/MWhH2)
Emission intensity (tCO2e/MWhH2)
1.0 1.0 1.0
0.8 0.8 0.8
0.6 0.6 0.6
0.4 0.4 0.4
0.2 0.2 0.2
0.0 0.0 0.0
−0.2 −0.2 −0.2
−0.4 −0.4 −0.4
−0.6 −0.6 −0.6
−0.8 −0.8 −0.8
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 0.0 0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4 2.7
RE factor RE factor RE factor
(d) Burnie (TAS) (e) Spencer Gulf (SA)
1.0 1.0
Emission intensity (tCO2e/MWhH2)
Emission intensity (tCO2e/MWhH2)
0.8 0.8
0.6 0.6
0.4 0.4
EI_RECs
0.2 0.2 EI_Market
0.0 0.0 EI_MEF
−0.2 −0.2 EI_L
−0.4 −0.4
−0.6 −0.6
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 0.0 0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4 2.7 3.0
RE factor RE factor
Figure 2 (a)-(e) Greenhouse gas (GHG) emission intensities of hydrogen under various RE factors, defined as the ratio between RE generators and electrolyser capacity in five states of the NEM, calculated using three different methods for electricity emissions. Actual emissions are tracked using time and state dependent marginal emission factors to net grid electricity consumption (EI_MEF, orange curve). Certified emissions are calculated using the market-based (EI_Market, blue curve) and location-based (EI_L, green curve) adopted in the PGO scheme.
Negative emissions intensity values indicate that more renewable electricity is sold to the grid than is bought per unit hydrogen produced, resulting in a reduction of grid emissions (for location and MEF methods). When the system generates more RECs that results EI_Market = 024 and grid electricity emissions can be reduced by excess
RECs generated for each MWh hydrogen produced (EI_RECs, green bars). The conversion rate between energy content units of hydrogen throughout the manuscript is based on its higher heating value (39.4 kWh/kg H 241).
We calculate the emission intensity of hydrogen produced by hybrid systems in the five states for systems that maintain a constant electrolyser size (equal to that of an optimized off-grid system, around
10-50MW) but have varying levels of dedicated RE capacity (Figure 2). The dedicated RE capacity is co-located ensuring geographic correlation and the systems can interact with the grid flexibly without temporal correlation. The RE factor is defined as the ratio between RE generators (sum of wind and solar) and electrolyser capacity: an RE factor of 0 is a fully grid-reliant system, while an RE factor of 2 indicates the system has RE capacity two times larger than the electrolyser. For each location, the largest
RE factor is set to be 1.5 times greater than the RE factor for the optimised off-grid system. By adjusting
6
the RE factor from 0 to its maximum value, we explore the variation in the emission intensities of hydrogen using the three different accounting methods.
In general, comparison between the emission intensities calculated using the AEF-based and PGO location-based methods, given in Figure 15, show relatively good agreement for all systems investigated, as they both use state-based AEFs, albeit at different temporal resolutions (1 hour for
EI_AEF, 1 year for EI_L)42. However, the location-based method does not reproduce the actual emissions intensity well, calculated with MEF-based method. It overestimates actual emissions in QLD,
NSW, and VIC by up to 63% and underestimates actual emissions for SA by 42% (RE Factor = 0). This is because MEFs for states dominated by coal (QLD, VIC, NSW) tend to be lower than the AEFs used in the location-based method, as over 49% of the time the marginal generator will be a lower emissions technology like gas or renewables (Figure 13). Conversely, SA tends to have MEFs higher than AEFs as coal is more likely to be the marginal generator when importing electricity.
The PGO market-based method overestimates the emissions for fully grid reliant systems (RE Factor =
0) and overestimates the electricity emissions avoided by self-sufficient systems (maximum RE Factor) by generating excess RECs. The largest discrepancies are observed for TAS and SA, with emissions intensities overestimated by up to 274%, and EI_RECs overestimated by up to 245% in TAS. Deviations reduce as less grid-electricity is bought and sold.
We attribute these discrepancies to the adoption of nation-wide factors in the PGO market-based method
(Eq. 18a). The proportion of RE on the state grids can be much higher (93% and 71% in 2023 for TAS and SA respectively, see Table 7) than the national-wide applicable renewable power percentage (ARPP) of 18.72%. Likewise, the emissions intensity used in the calculation (RMF) is calculated by adjusting the average emission factor of the grid at the national level (including the NEM, South-West
Interconnected System, and the Darwin, Katherine Interconnected System43) by excluding all RE that has generated RECs (Eq. (27)). This overlooks the energy mix of regional grids, resulting in larger discrepancies for states with higher RE penetration, overestimating emissions when buying grid electricity.
Impact of temporal correlation
7
Burnie (TAS) Newcastle (NSW) Geelong (VIC)
45 45
7 7 40 7
40 40
6 35
Emission intensity (kgCO2e/kg H2)
Emission intensity (kgCO2e/kg H2)
Emission intensity (kgCO2e/kg H2)
6 35 6 35
5 30
5 30 25 5 30
4
25 20 25
LCOH ($/kg H2)
LCOH ($/kg H2)
LCOH ($/kg H2)
4 3 4
20 15 20
3 2 10 3
15 15
2 1 5 2
10 10
0 0
1 5 1 5
−1 −5
0 0 −2 −10 0 0
−5 −15 −5
−1 −3 −1
−10 −20 −10
−2 −4 −25 −2
−15 −5 −15
−3 −30 −3
−20 −6 −35 −20
id ly
ai
ly hl
y rl
y
ib
le ≤0 id ly
ai
ly ly rl
y
bl
e ≤ 0 id ly
ai
ly hl
y rl
y
bl
e ≤0
gr on nt ea ex EF gr on th ea xi EF gr on nt ea xi EF
ff- ll D o Y Fl ff- ll D on Y le ff- ll D o Y le
O Se M I_M O Se M F _M O Se M F _M
E EI EI
Spencer Gulf (SA) Gladstone (QLD)
45 45
7 7
40 Emission intensity (kgCO2e/kg H2) 40
Emission intensity (kgCO2e/kg H2)
6 35 6 35
LCOH
5 30 5 30
25 25
EI_Market
LCOH ($/kg H2)
LCOH ($/kg H2)
4 4
20 20 EI_MEF
3 3
15 15 Grid electricity cost
2 2
10 10 Electrolyser CAPEX+OM
1 5 1 5
Hydrogen Storage CAPEX
0 0 0 0
−5 −5
PV CAPEX+OM
−1 −1
−10 −10 Wind CAPEX+OM
−2 −2
−15 −15
−3 −3
−20 −20
id nl
y
ai
ly ly rl
y
ib
le ≤0 id nl
y
ai
ly ly rl
y
bl
e ≤0
gr o th ea ex EF gr o th ea ex
i
EF
ff- ll D on Y Fl M ff- ll D on Y Fl
O Se M _ O Se M _M
EI EI
Figure 3 Total LCOH (white dots), average yearly actual (EI_MEF, yellow triangles) and certified emission intensities of hydrogen (EI_Market, red triangles) calculated using the PGO market-based method, optimized system configurations in five states across different under different constraints described in Table 1. The emission intensity is given in kgCO2e/kgH2 to allow comparison with existing certification schemes. Bars show LCOH deconvoluted to components: grid electricity cost (sum of the total expenditure and profit from buying and selling electricity to and from the grid, which can be negative), capital expenditure (CAPEX) of hydrogen storage, and
CAPEX and operation and maintenance costs (OM) of electrolyser, PV and wind generators. Negative emission intensity values for market-based method (EI_Market) quantify grid electricity emissions offset by excess RECs generated for each kg hydrogen produced, while they represent the actual emissions displaced for MEF-based method (EI_MEF).
Table 1 Scenario design
Scenario Grid connection Temporal correlation Emission intensity constraint
Off-grid NO None None
Sell only Yes – but only sell electricity Hourly None
Daily Yes – buy and sell Daily None
Monthly Yes – buy and sell Monthly None
Yearly Yes – buy and sell Yearly None
Flexible Yes – buy and sell None None
8
EI_MEF≤0 Yes – buy and sell None ≤0
In the following, we use our model to optimize the configuration and operation of hybrid grid-connected hydrogen systems to minimize LCOH under constraints, as described in Table 1. We compare the scenarios include the reference off-grid scenario where no grid interaction is permitted; four scenarios requiring temporal correlation over different intervals (hourly - which is equivalent to only allowing the selling of electricity to the grid labelled “Sell only”, daily, monthly, yearly), while ensuring geographic correlation by co-locating the dedicated RE capacity, and a “flexible” scenario with no requirements for temporal correlation. To assess the effectiveness of temporal correlation policies for lowering emissions, we also include a scenario where the system can interact with the grid without temporal correlation, and instead constrain the “actual” emissions intensity of hydrogen to be less than zero, i.e. EI_MEF≤0.
This final scenario is not meant to model a realistic policy setting, but to understand what a system optimised for zero-emissions at lowest cost looks like. We focus on the PGO market-based method, that allows the trading of RECs and explore the impact of different temporal correlation policy settings on reducing actual emissions, and those certified by the PGO scheme.
LCOH is highest for off-grid production systems and reduces with grid connection for all states, with the lowest LCOH occurring for fully flexible scenarios with no additional constraints (Figure 3). There are two distinct trends for optimized system configurations and emission intensity of hydrogen: one that holds for TAS, NSW and VIC, and one that holds for SA and QLD.
Due to the interplay between spot prices and the quality of the RE resources, the cheapest way to produce electrolytic hydrogen in TAS, NSW and VIC is to minimize capital investment and source large amounts of power from the grid. For TAS and VIC, relatively low grid electricity prices make dedicated
RE uncompetitive. While NSW has the highest yearly average electricity price in the NEM, the local
RE resources in Newcastle are less abundant (quantified by the RE capacity factor, Table 6), leading to higher RE costs. When unconstrained by temporal correlation in the flexible scenario, the system builds less dedicated RE capacity in TAS and NSW compared to the off-grid scenario, and none in VIC, and invests in smaller hydrogen storage and electrolyser capacity. Critically, relying on grid power in the lowest cost flexible scenarios can result in very high hydrogen emission intensities, with EI_MEF =~25 kgCO2e/kgH2 for VIC and NSW, comparable to hydrogen made from coal gasification44. Once again, we see that the PGO market-based method overestimates the certified emissions intensity, particularly for TAS and VIC.
Once yearly temporal correlation is imposed in TAS, NSW and VIC, the system is forced to build RE capacity, the LCOH increases, and the actual emission intensity falls to below zero in TAS and slightly
9
above zero in NSW and VIC. As the temporal correlation requirements become more restrictive, shifting from yearly to hourly time matching, the system increases the capacity of the electrolyser and storage, and the configuration progressively resembles an off-grid setup. Consequently, LCOH increases by
30%-38%. Once again, the PGO market-based emission calculation overestimates potential emissions offsets by RECs, becoming more significant when the system sells larger amounts of electricity to the grid.
The trend is different for SA and QLD. All systems have actual emissions intensities close to or below zero, and the flexible and yearly temporal correlation scenarios all result in nearly identical system configurations, with similar and lowest cost LCOH. The cost reduction in these two scenarios compared to the off-grid system is mainly attributable to the systems’ ability to profit from selling surplus electricity to the grid, driven by favourable local RE resources and higher regional electricity prices
(Table 6, 7). The system significantly oversizes the cheapest RE capacity and sells large amounts of electricity to subsidize the cost of hydrogen production. As the temporal correlation requirements become more restrictive, shifting from yearly to hourly time matching, the restrictions on importing electricity from the grid results in the establishment of configurations similar to the off-grid system with
PV and larger storage components. The systems are prevented from subsidizing the cost of hydrogen by selling electricity, increasing the LCOH by 78 % in QLD, and 52% in SA.
One way to assess the efficacy of temporal correlation as a policy is to compare with systems optimized to minimize cost while constraining EI_MEF≤0. Figure 3 shows that such systems are most similar to those optimized under yearly temporal correlation, suggesting it is an efficient policy measure to keep emissions low, at low cost. However, it is important to note that while yearly temporal correlation is sufficient to ensure that the certified emissions intensities under the PGO market-based method are all below zero, actual emissions intensities are still 1.9 kgCO2e/kgH2 for the VIC case and 3.6 kgCO2e/kgH2 for the NSW case. These emissions intensities are similar to those of natural gas based blue hydrogen systems1.
Impact of geographic correlation
10
Hydrogen production system
in Newcastle (NSW)
7.0
(a) LCOH
6.5
6.6
6.5 LCOH with direct line
LCOH with RECs cost
6.0
LCOH ($/kg H2)
5.5
5.1
5.0 5.2
4.8 4.8
4.6
4.5
4.5
4.6
4.0
3.8
3.5
40 (b) 38
35 EI_Market
Emission intensity (kgCO2e/kg H2)
30 EI_MEF 32
25 21
20
15
10
10 7.6
3.7
5
0
−5 −8.7 −8.7 −8.7 −8.7
−3.6
−10
−15
−20 −24
−25
−30
VIC QLD SA TAS NSW Fully-grid Off-grid
(Geelong) (Gladstone) (Spencer Gulf) (Burnie) (Newcastle) (Newcastle) (Newcastle)
RE generators' location
Figure 4 (a) LCOH and (b) emission intensity of hydrogen produced by a grid-connected electrolyser located in
Newcastle NSW, with RECs sourced from grid-connected RE capacity located different states, as indicated along the x-axis. Systems are optimised under the constraint of yearly temporal correlation and are compared to off-grid and fully-grid connected systems in NSW. The green line indicates the estimated cost of purchasing RECs from a third party to offset emissions for the fully on-grid case, assuming REC costs of 20-60AUD33. Note that the LCOH for a system with RE located in NSW connected to the electrolyser through the grid is slightly higher than that calculated previously (see Figure 3) for a system with direct connection between co-located RE and electrolyser in NSW (red dot) due to increased transmission use of system fees. The emission intensity of hydrogen is calculated by MEF-based method (yellow triangles) and for the PGO market-based (red triangles).
So far, all hybrid systems have been constrained by geographic correlation, with a direct connection between co-located RE and hydrogen electrolysers. Next, we relax that constraint and compare scenarios in which the electrolyser and the dedicated RE generators are located in different states and
11
interact only through grid-connection, as would occur under certain pre-purchase agreement39. In these scenarios, the RE is sold to its local grid, generating RECs and offsetting actual emissions in that state, while the electrolyser buys electricity from its local grid, generating actual emissions and using the
RECs generated by the RE to offset emissions (see Modelling geographic correlation in the Method section). We focus on an electrolyser located in Newcastle (NSW), which is the site of the announced
Hunter Valley Hydrogen Hub, slated to include a 50 MW electrolyser powered by grid-connected RE through purchase and retirement of RECs45. We compare systems under the constraint of yearly temporal correlation with dedicated RE located in different states, with a fully grid-powered system and an off-grid system (Figure 4). The LCOH can be reduced by as much as 21% when the RE is located in QLD due to the superior RE resources and high electricity costs, as the net cost of grid electricity is determined by buying and selling electricity to and from two different grids and depends on the spot prices in both states.
Results show that yearly temporal correlation is sufficient to ensure that all certified emissions
(EI_Market, red triangles) are below zero even without geographic correlation, as each unit of electricity used by the electrolyser is compensated for by a REC generated within the year. However, the actual emissions (EI_MEF, yellow triangles) vary widely with the location of the RE as the emissions displaced by selling RE will depend on the makeup of its local grid and are not necessarily equal to emissions generated using electricity in NSW. In the worst case, the emission intensity can reach 21 kgCO2e/kgH2, similar to that of hydrogen produced via coal gasification44, as high emissions intensity electricity consumed in NSW is being offset by displacing relatively low emissions intensity electricity in TAS.
In all scenarios so far, we assume a strict additionality requirement: that dedicated RE is deployed concurrently to the electrolyser. The PGO scheme does not require such strict additionality as RECs can be generated by eligible RE systems built before the electrolyser, as long as they began generating renewable electricity after 1997, and do not have to be surrendered as part of Renewable Energy
Target46–48 or any other scheme. The fully grid-powered electrolyser scenario in Figure 4 approximates this case, buying electricity from the NSW grid at the spot price and generating the corresponding emissions according to the state-based MEF. Under this scenario, the additional cost for RECs to offset the certified emissions under the market-based method would increase the LCOH by 16-44% assuming
RECs costs of AUD 20-6033, while estimating the actual emissions would require the time and state dependent MEF for each REC.
45
The LCOH of the NSW system connected through the grid is slightly higher than when it is directly connected with RE, as plotted in Figure 3, since electricity purchases from the grid and incurs a “transmission use of system fee” (see Table 3)
12
Conclusion
Australia has the potential and the ambition to be a producer and exporter of large volumes of low- emissions hydrogen. This paper evaluates the efficacy of proposed emissions accounting methodologies and certification scheme rules for hydrogen production in Australia’s PGO certification scheme, which is currently under development.
We employ an energy system model to quantify emissions and costs of hybrid grid-connected hydrogen production systems located across Australian NEM. By imposing different constraints, we assess the accuracy of the PGO certification scheme methodologies and the efficacy of different policy settings in minimizing emissions while maintaining low costs.
We find that PGO emissions accounting methods can lead to significant discrepancies compared to the actual emissions. The differences vary across states within the NEM due to differences in state grid generator mix and renewable resource availability. The PGO location-based method overestimates emissions when AEFs are larger than MEFs (e.g. 63% for fully grid reliant in VIC), and underestimates when they are smaller (e.g. 42% for fully grid reliant in SA). This occurs because actual emissions are determined by the marginal generator, not the average grid emissions intensity. The high frequency of hydro as a marginal generator leads to overestimation in heavily coal-reliant states with high AEFs
(QLD, NSW, VIC), while coal results in underestimation in high renewable penetration states with low
AEFs (SA).
The PGO market-based method tends to overestimate emissions in all states when large amounts of electricity is bought (e.g. 274% for fully grid-reliant systems in TAS). The discrepancies are larger for states with more RE penetration as the market-based method uses factors relating to RE penetration and residual emissions (i.e. ARPP and RMF) that are currently estimated using the RECs information aggregated at the national level43. It is likely that the use of more granular, state-level factors would improve accuracy.
We note that the overestimation is consistent with the principle of conservativeness recommended for embedded emissions accounting6, and hence is not necessarily a design flaw, underestimation of embedded emissions contravenes the conservativeness principle and undermines the scheme’s emissions-reduction goals.
We further investigate the effect of PGO policy settings on the cost as well as the real and certified emissions when trading renewable energy certificates (i.e. PGO market-based method). By ensuring that the systems only use RE from dedicated, co-located wind and solar capacity, we are implicitly applying additionality and geographic correlation criteria. We then examine different temporal correlation settings. Despite differences in cost optimal system configurations across the five states, our
13
main finding is robust and in line with previous work7,13: allowing less onerous temporal correlation requirements can reduce LCOH, while still limiting emission intensities. Specifically, optimizing systems operating under yearly temporal correlation results in similar system configurations as those optimized to keep certified emissions below zero at lowest cost (i.e. constraint of EI_MEF≤0). This is an important outcome that suggests that less stringent temporal correlation is an effective and efficient policy measure to ensure near-zero emissions systems at low cost and low regulatory overhead.
However, in certain circumstances, discrepancies between the actual and certified emissions calculated under the PGO market-based method could lead to systems that still generate hydrogen with actual emissions under yearly correlation, even though they are certified as zero emissions, particularly when the geographic correlation is not met (e.g. 21 kgCO2e/kgH2 in Figure 4). Ensuring actual emissions are below zero for all states will require geographic correlation coupled with temporal correlation intervals somewhere between monthly and yearly and will change with location and over time, as it is determined by the interplay between regional electricity prices, local renewable resources and RE penetration in the local state grid.
Our findings provide timely evidence to aid the design of the Australian PGO scheme as well as evidence to support the development of robust and practical certification schemes internationally.
Method
Model description
To quantify the cost and emission intensity of hydrogen production by the hybrid grid-connected hydrogen production system, we establish a linear optimization model to determine the optimal capacity of the system components consisting of dedicated wind and solar, an electrolyser, two compressors, hydrogen storage facilities and pipeline, and optimal interaction with the grid to ensure the lowest hydrogen production cost. The model is operated with hourly time resolution over a one-year period and is constrained to deliver a constant supply of hydrogen to industrial end-users throughout the year.
The modelled system can use the dedicated RE capacity to supply power for electrolysis and subsequent compression processes, (directly or through the grid) or can choose to buy or sell electricity from the
Australian NEM. The system can exchange unlimited power with the grid however RE generation and electrolyser capacities are kept within the range of 10MW to 50MW so that the effect on the price and emission intensity of grid electricity can be neglected. These values were chosen to reflect the capacity of proposed hydrogen production projects in Australia49, and the average instantaneous electricity demand of the smallest regional grid in the NEM, Tasmania, which is approximately 1200MW50. Local wind and solar electricity generation are calculated using NREL’s PySAM engine51 based on historical local weather data from NASA’s MERRA-2 dataset52. The electrolyser converts electricity into
14
hydrogen with an efficiency of 70%53. We assume that the hydrogen produced by the electrolyser is initially at a pressure of 30 bar and compressor 1 pressurizes hydrogen to 100 bar for direct pipeline transport to the load, while compressor 2 increases the pressure to 150 bar for storage in underground lined rock caverns (LRC) if needed. The H2 storage acts as a buffer to the load, allowing it to be stockpiled when local RE generation is abundant and grid electricity is cheap, and to be supplied when local RE generation is scarce, and electricity is expensive. The system could choose to store electricity in batteries; however, due to the high capital expenditure of batteries compared to hydrogen storage, the model optimization process favours investing in hydrogen storage. The system can also choose to import and export electricity from the state grid (bidding zone) in which it is located, defined by time and location dependent historical electricity prices and emissions factors, described in detail in section
Modelling the Australian National Electricity Market below.
Hybrid grid-connected hydrogen production optimization model
Figure 5 Overview of the hybrid grid-connected hydrogen production optimization model. Ovals represent nodes, which are system components responsible for transformation of energy commodities. Lines represent edges which indicate the flow of an energy commodity and arrows represent the flow directions. Orange lines indicate electricity flow (kW), while green lines represent hydrogen flow (kg/hour). The electricity connection point serves as the node ensuring electricity balance, where electricity consumption aligns with electricity generation. The grid node represents the Australian NEM, and the hydrogen production system can participate in the Australian NEM via the transmission line, either as a generator selling electricity or as a consumer purchasing electricity.
The linear optimization model is shown schematically in Figure 5, where ovals correspond to nodes representing system components for converting between energy commodities, lines representing edges which indicate the energy commodity flows, with arrows indicating the direction of these flows. The variables and indicate flows of electricity in kilowatts and hydrogen in kilograms per hour between different component, respectively, where the superscript specifies the component, and
15
indicates the flow direction. The production system simulates operations by ensuring the equilibrium of electricity and hydrogen flows, as described in the following.
Electricity flow balance
The production system is powered by electricity which can be sourced from local RE generation or
()) imported from the grid. At each time , local wind generation ( ()) and solar generation ( are produced based on the local weather profile. Renewable energy (RE) can either be used to (i) power
the electrolyser for hydrogen production ( ()) and the compressor for hydrogen transmission into
1 2 either the pipeline ( ()) or storage facilities ( ()), (ii) exported to the grid ( ()), or
(iii) curtailed ( ()), where the variables refer to power flows shown in Figure 5. If the local RE is insufficient to support the system’s operation, the system can choose to import electricity from the grid
( ()). The decision of electricity dispatching is made by solving the flow balance equation at the
Electricity Connection Point given by
1 2 () ()
() + () + () + () + () = + + () ∀ (1)
Hydrogen flow balance
The system converts electricity into hydrogen through electrolysis and then either pumps the hydrogen into the pipeline for direct transportation to the hydrogen load or stores it in the hydrogen storage equipment, which serves as a reserve for the load, storing hydrogen to supply when production by the electrolyser is insufficient to meet the demand. At each time , hydrogen is produced by the electrolyser with efficiency, = 70%53, as following:
()×
()
=
∀ (2)
where =39.4 kWh/kgH241 is the higher heating value of hydrogen and () represents the hydrogen produced by electrolyser. This hydrogen is subsequently allocated by the system to either
1 compressor 1 ( ()), which pumps it into the pipeline for direct transportation to end-users at a
2 pressure of 100 bar, or to compressor 2 ( ()), where it is pressurized to 150 bar and then stored in the hydrogen storage:
() 1 2
= () + () ∀ (3)
The hydrogen load in each time ( () = 180) is met either by hydrogen flow directly from
the electrolyser via pipeline and/or hydrogen from storage ( ()) according to:
1
() = () + () ∀ (4)
16
Wind and solar generation
Local wind and PV generation at each time depends on the optimized capacity of the generators. To
enable the optimization, we define reference capacities for the wind farm ( = 320 MW)54 and PV
field ( = 1 MW), and calculate the hourly reference generation (
and
()) using
NREL’s PySAM engine51 and historical local weather data from MERRA-2 dataset52. The actual generation is then calculated through linearly scaling the reference output based on the ratio of actual capacity to reference capacity:
() = ×
() ∀ (5)
() = × () ∀ (6)
Where () and () represent the actual hourly generation of the wind farm and the PV field, respectively. and indicate the capacity for the wind farm and PV field, respectively.
Grid node
We introduce the grid node to simulate the interaction between production system and the grid. The grid is defined by historical electricity prices (spot prices) and emissions factors (average and marginal) data sourced from the NEM. The modelling approach is introduced in the following section Modelling the
Australian National Electricity Market. Detailed information on the NEM data is provided in following section NEM Data.
Compressor 1 and compressor 2
The electricity consumption of compressor 1, which pressurizes hydrogen produced by the electrolyser for pipeline injection, and compressor 2, which pressurizes hydrogen for storage is given by
1
1
() = () × 1
∀, (7)
2
2
() = () × 2
∀, (8)
1 2 where and are the electricity consumption for each unit of hydrogen compressed by compressor 1 and compressor 2, respectively.
Hydrogen storage
The hydrogen storage level at each time ( _ ()) is updated according to:
54
Note that a wind farm requires a layout of multiple turbines to account for the wake effect between them. Here, we assume that 320 MW capacity represents a sufficiently large plant to consider a turbine layout.
17
2
() = 0 + ∑0 (−
()
+ ()) ∀, (9)
where () represents the hydrogen sent out by the hydrogen storage to meet the load and __0 is the initial level of hydrogen storage as we assume the system has been in trial operation for a period.
The hydrogen storage level in the last time point should be equal to the initial level to ensure that all the hydrogen sent out is produced within the modelled year.
Optimization
The objective of the optimization model is to minimize the LCOH, defined as55:
×+&+
= 2
, (10) where 2 is the total hydrogen produced annually, is the capital expenditure including installation and equipment cost53, O&M is the operation and maintenance costs, is the capital recovery factor, and is the electricity cost over the year. The CAPEX and O&M are given by the sum of the costs associated with individual system components:
= ∑∈ , (11)
& = ∑∈ ( + ), (12) where represents the installed capacity of a component, is investment needed per unit installed capacity, and are the fixed and variable operation and maintenance cost of component , where is the component set which includes wind, PV, electrolyser and 2 storage facilities. The can be calculated according to the plant lifetime and the interest rate as:
(1+)
= (1+) −1, (13)
The electricity cost over one year includes costs associated with buying electricity from the grid and negative “costs” from selling electricity to the grid in each time and is given by,
= ∑8759
=0 ( () × ( () + ) − () × ()), (14) where () is the electricity spot price in state at each time . In addition, the system will need to pay an additional transmission use of system fee when importing electricity37. The electricity into the electrolyser and the hydrogen storage level are both limited by the capacity of electrolyser , and hydrogen storage , respectively:
() ≤ ∀, (15)
() ≤ ∀, (16)
18
In addition to the constraints for each component and the flow balances established above, the system can incorporate additional constraints such as temporal correlation and limits on emissions intensity.
The temporal correlation constraint requires electricity consumed by the electrolyser to equal dedicated
RE generation over certain time intervals, following previous work14. In the model, this constraint mandates that the total electricity sold to the grid must be larger than or equal to the total purchased electricity within each time interval:
∑ +(
() − ()) ≥ 0 ∀ ∈ , (17) where indicates the length of the time interval required, and marks the beginning of the temporal correlation interval, which falls within the set of time steps tailored to meet the specified time interval.
In this work, ranges from one hour to one year. For example, for a monthly correlation requirement,
will include the initial time points of each financial month.
Other policy settings can be incorporated into the model by introducing additional constraints.
Additionality is ensured as dedicated RE capacity is installed at the same time as the electrolyser, and geographic correlation is applied by assuming that electrolysers are in the same state as the dedicated
RE capacity (i.e. weather data and grid data are state specific). The system can also constrain the emission intensity for hydrogen to below zero to assess the impact on LCOH.
Constraints are also placed on the size of the RE capacity to avoid unphysical results. For example, in states with high electricity prices, the system may grossly oversize the local RE generation capacity to profit from selling electricity to the grid, in essence creating a second business to subsidize hydrogen production. This phenomenon does not align with the practical constraints as the RE generation capacity is usually limited by the available land, project budget and the risk of investment 37. Thus, we also constrain the of the on-grid system such that it does not exceed the of the optimized off-grid system without grid-connection.
Methods for calculating the electricity emissions
The emissions associated with the electricity used for hydrogen production are calculated using four different methods. Two of these methods are defined by the Australian PGO scheme20, differentiated by whether they include the buying and selling of RECs, known as the market-based method, or use a location dependent emissions factor defined by the Australian government24, known as the location- based method. The other two methods take advantage of the ability of our model to assign time and location dependent MEFs or AEFs to the net electricity consumption at each point in time.
For the PGO market-based method, electricity emissions ( ) from time 1 to 2 are calculated using the applicable renewable power percentage (), which gives the fraction of the electricity taken
19
from the grid that comes from RE source, and the residual mix factor (), which is the estimated average emission factor for the grid that excludes all the electricity generation claimed by RECs43. The emissions of electricity used to produce hydrogen is then
= ( × (1 − ) − × 1000) × . (18a)
2
= (∑= 1
( () × (1 − ) − ()) × . (18b)
Eq. 18a shows the general case in which the residual electricity, calculated as the total purchased or acquired electricity that is not classified as renewable minus the eligible RECs, is multiplied by the
RMF43, where is the quantity of grid electricity used and is the quantity of RECs surrendered, and Eq. 18b shows how we apply the market-based method in our model from time 1 to
2 . We assume that each MWh of RE sold to the grid generates one RECs such that, =
()/1000, which can be surrendered to cover emissions from using grid electricity. If more eligible RECs are surrendered to estimate electricity emissions than the total required to reach zero emissions, is equal to zero24. Here, we assume that excess RECs will be traded to other facilities, enabling them to claim reductions in their electricity emissions under the market-based method.
Consequently, one unit of eligible RECs can represent an emissions reduction on the grid of 0.81 tCO2e in the market-base method.
To calculate the electricity emissions using the location-based method ( ), the location dependent emissions factor for each state ( ) is defined by the Australian government as part of the National
Greenhouse Gas Accounting Factors24. The location-based emissions factor is based on the AEFs in each state averaged over the year43 to calculate the emissions intensity of the electricity used to produce hydrogen:
2
= ∑= 1
( () − ()) × (19)
If the amount of () < () , then the system contributes to the reduction of grid GHG emissions.
We define two other emissions calculation methods, using the time and location dependent MEFs or
AEFs and the net electricity consumption at each point in time. The MEF-based and AEF-based method keep track of the electricity emissions ( and ) from buying and selling electricity to the grid from time 1 to 2 and are defined as
2
= ∑= 1
(( () − ()) × ()) (20)
2
= ∑= 1
(( () − ()) × ()) (21)
20
Where the MEFs and AEFs of the grid in Australian states connected to the NEM at each time are represented as () and () , respectively. In both cases, the values can be negative, representing the case when more RE is sold to the grid than is bought in each time, and interaction of the system grid with the system is reducing the grid emissions.
Modelling the Australian National Electricity Market
The NEM is made up of networks in five regions (roughly corresponding to the states) interconnected by a transmission network, including Queensland (QLD), New South Wales (NSW), South Australia
(SA), Victoria (VIC) and Tasmania (TAS). Electricity is traded within and between states by the central dispatch engine run by the Australian Energy Market Operator (AEMO)56. The spot prices of electricity in each state are determined by the highest bid accepted by the AEMO from a generator to fulfill demand for each 5-minute interval. The historical data on spot prices and emission factors used in our study is from 2023 sourced from AEMO’s dataset NEMweb at five-minute intervals35, and subsequently averaged to one-hour timesteps. Here, a representative historical time series is shown in Figure 6a, depicting hourly AEFs and MEFs for Queensland grid profile on 18 Feb 2023. The associated hourly spot prices for the same time and location are shown in Figure 6b, along with a bar chart indicating how frequently each type of generator (coal, gas, renewables) becomes the marginal generator over the five-minute intervals in each hourly dataset.
21
1×105 1.0
1×10 5 (a) Solar Wind Hydro Fossil MEF AEF
0.9
5
1×10
0.8
9×104
Emission factor (kgCO2e/kWh)
Electricity generation (MW)
0.7
8×104
0.6
7×104
0.5
6×104
5×104 0.4
4×104 0.3
4
3×10 0.2
4
2×10
0.1
1×104
0.0
Frequency of marginal generator occurrence
(b) Black coal Brown coal Coal seam methane Hydro Natural Gas Solar Wind Battery Storage 0.28
18
Electricity Price 0.24
Electricity Price (AUD/kWh)
16
14 0.20
12
0.16
10
0.12
8
6 0.08
4 0.04
2 0.00
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23
Hour
Figure 6. The example of the grid profile in QLD on 18 Feb 2023. (a) shows marginal emissions factors (MEFs)
(purple curve) and average emissions factors (AEFs) (green curve) in relation to the electricity generation of the grid (bars) in each hour. (b) shows the frequency of each type of fuel generator becoming marginal (bars), along with electricity prices (blue curve) for each hour.
The graphs show that AEFs reflect the generation mix of the grid, with high AEFs during the early morning and evening peaks due to the reliance on fossil fuel-generated electricity, and lower AEFs during midday when the proportion of RE generation on the grid increases. In contrast, the MEF values are determined by the emissions intensity of the marginal generators and so can change rapidly across the day. The two emission factors can vary widely over some time periods. For example, in hour 18 the
AEF is 0.82 kgCO2e/kWh due to the generation mix including RE and fossil sources. At the same time, the MEF plummets to zero, as hydro is the marginal generator and will be used to meet any additional demand. In this work, we consider MEFs to best represent the actual electricity emissions caused by the additional demand on the grid from our system. For example, if our system purchases electricity at hour 18, the generators will be hydro as shown in Figure 6b, resulting in no emissions from our purchase. In addition, it tells us which generators would be displaced by RE sold from our system into the grid (and thus quantifies avoided emissions). For instance, if the system sells 1 kWh of renewable
22
electricity at hour 7, the displaced generators will be black coal, reducing grid emissions by 0.88 kgCO2e. Further discussion about the comparison between MEFs and AEFs is given in the NEM data section below.
It is important to understand how the generator mix and emission profiles of the state grids differ, and how this is reflected in the AEF and MEF values. Yearly average AEFs and MEFs for each state are given in Table 7, and the energy generation mix as quantified by the percentage contribution of different generators to annual generation is shown in Figure 13. Local grids in QLD, VIC and NSW are dominated by coal power (~60%)26, with RE meeting a significant portion of the remaining demand,
(also given in Table 7 as 28%, 42%, 31% respectively) leading to average AEFs of 0.69, 0.73, and 0.63 kgCO2e/kWh in 2023. In contrast, 71% of SA electricity was provided by renewables in 2023, with the rest coming from gas power and imports from other states26, resulting in a much lower average AEF of
0.23 kgCO2e/kWh. The average AEF in TAS was lower still, at 0.12 kgCO2e/kWh, as over 90% of electricity generation was renewable, mostly from hydro (~70%), followed by wind (~20%).
The average MEF values vary much less between the states, from 0.4-0.52 kgCO2e/kWh, with the exception of TAS, which has a much lower average MEF of 0.19 kgCO2e/kWh. This is because the marginal generators that respond most frequently to increased demand in each state are similar: coal
(33-50% of the time), followed by hydro (25-32%), as shown in Figure 13. Since TAS has abundant hydro, it is the marginal generator over 75% of the time. The high frequency of hydro as the marginal generator leads to MEFs that are lower than AEF in heavily coal-reliant states (QLD, NSW, VIC), while the fossil fuel technologies are marginal generators 45% of the time in SA, resulting in larger MEF than
AEF.
Modelling geographic correlation
Figure 7 Overview of hybrid grid-connected hydrogen production optimization model without the geographic correlation constraints. Ovals represent nodes, which are system components responsible for transformation of
23
energy commodities. Lines represent edges which indicate the flow of an energy commodity and arrows represent the flow directions. Orange lines indicate electricity flow (kW), while green lines represent hydrogen flow
(kg/hour), and pink lines show renewable energy certificates (RECs) associated with renewable electricity integrated into the grid. The dashed lines indicate the boundary between two states. The Sell Grid node locates in
State 1, where the renewable generators (wind farm and PV field) integrate electricity to the grid, while the hydrogen production system and the Buy Grid node are situated in State 2. State 1 and State 2 are interconnected two bidding zones via a transmission line.
Apart from the scenario in which renewable generators (the wind farm and solar field) are co-located with the electrolyser, meeting the geographic correlation requirement as shown in Figure 5, our study uses the same model but with modified constraints to study a scenario where the electrolyser is not co- located with renewable generators (without the direct line connection) to investigate the policy impact of geographic correlation on the emissions and cost of hydrogen production as shown in Figure 7.
Under this scenario, we assume that the hydrogen producer invests in self-owned renewable generators located in states different from where the hydrogen production system is situated 57. All renewable electricity generated will be sold to the regional grid (the Sell Grid node), generating profit and producing RECs in the process. The hydrogen production system is grid-connected with the regional grid (the Buy Grid node), leveraging the NEM’s transmission network to match its electricity consumption with the self-invested renewable generators using RECs information.
Distinct from Eq. (14) above, the electricity cost over one year includes costs associated with buying electricity from the Buy Grid and negative “costs” from selling electricity to the Sell Grid in each time
and is given by,
= ∑8759
=0 ( () × ( () + ) − () × ()), (22) where () and () are the electricity spot price in the Buy Grid and Sell Grid at each time , respectively. For example, the hydrogen producer establishes a hydrogen production system in
Newcastle (NSW), including a grid-connected electrolyser, compressors, hydrogen storage, and transportation pipeline facilities, while developing a wind farm in Burnie (TAS). The electricity purchase cost will be determined by the Buy Grid Electricity price in NSW ( ()), while the revenue from electricity sales will be based on the Sold Grid Electricity price in TAS ( ()). The RECs generated by the Burnie wind farm will be used to certify the electricity consumed by the hydrogen production system in Newcastle.
Data input
Selected locations
24
Figure 8 The five locations selected as the case study.
Table 2 The latitude and longitude of the five locations
State Location Latitude Longitude
Queensland Gladstone -24.08 151.28
South Australia Spencer Gulf -32.94 137.45
Tasmania Burnie -40.89 145.25
Victoria Geelong -38.07 144.39
New South Wales Newcastle -32.92 151.75
In this study, we select one location for a hydrogen production plant in each of the five states participating in the Australian National Electricity Market (NEM) as the example projects using
Australian hydrogen projects dataset as the guidance49. The database provides a comprehensive overview of all hydrogen projects across Australia, including detailed information such as the energy sources used for hydrogen production, the production methods employed, and the annual hydrogen production amounts for each project. We focus on projects that source their energy from the grid and select those that are also participating in the current trial projects under the PGO scheme58.
The geographical locations of the five cases we study, all of which have available grid connection conditions, are shown in Figure 8, where the purple dots represent the locations of the hydrogen
25
production systems, and the red line indicates the grid electricity transmission line sourced from
Geoscience Australia dataset59. The specific latitude and longitude are listed in Table 2.
Model input technical and economic parameters
The input technical and economic parameters for the optimization model are chosen from best available data sources in the literature and listed in Table 3. In particular, we use data from the latest Gencost report as extensively as possible. This economic report estimates the costs of building new electricity generation and hydrogen production facilities specific to Australia60. Costs are given for 2023 and converted to US dollars with an exchange rate of 0.7 US dollars per Australian dollar which is in line with the exchange rate adopted in latest Gencost report60.
Table 3 Model input parameters
Parameters Expression Units Current Value
Capex of electrolyser USD/kW 1343.360
Capax of wind USD/kW 2126.660
Capex of PV USD/kW 1068.260
FOM of electrolyser USD/kW 37.461
FOM of wind USD/kW 17.560
FOM of PV USD/kW 11.960
VOM of electrolyser USD/kgH2 0.0262
Electricity consumption of pumping unit 1 kWh/kgH2 0.83
hydrogen into pipeline
Efficiency of electrolyser 70%53
Higher heating value of hydrogen kWh/kg H2 39.441
Discount rate 6%
Project lifetime 25
Transmission use of system fee USD/kWh 0.00763
In this study, we explore two options of hydrogen storage equipment: pipeline storage and the Lined
Rock Cavern (LRC) underground storage. The relationship between the total capex for both the pipeline
26
and LRC underground storage systems and their storage capacity are calculated using the following Eqs.
(23)-(24)64:
10 ( ) = −0.028510 (
1000
) + 2.7853 (23)
2
10 ( ) = 0.217956 (10 (
1000
)) − 1.57520910 (
1000
) + 4.463930 (24)
Where unit cost of pipeline storage and LRC underground storage are in USD per kg of H2, and their capacity is in kg of H2.
Figure 9 The unit capex of hydrogen storage as a function of the capacity.
Here, we utilize the piecewise function to linearize these two cost functions to reduce the solving complexity, and the curve of unit capex of those two types of hydrogen storage based on their capacity is shown in Figure 9. Noticeably, LRC underground storage capex decreases significantly and become more economical than pipeline storage for larger storage sizes. Hence, we assume the system will opt for pipeline storage if the hydrogen storage size is under 21742kg and the electricity consumption of
2 pumping unit hydrogen into the storage is 0.83 kWh/kgH2. Otherwise, the system will choose
2
LRC underground storage and is 1.24 kWh/kgH2, due to the higher-pressure requirement of underground storage (150 bar) compared to the pipeline (100 bar).
27
Renewable generation
Table 4 PV field model parameters.
Parameter Value
Rated Capacity 1MW
Array type 2-axis tracking
Tilt 0°
Azimuth 180°
System losses 15%
Table 5 Wind farm model parameters.
Category Parameter Value
Turbine Model 2020 ATB NREL Ref. 4 MW
Hub height 150m
Number of turbines 80
Rated Capacity 320MW
Availability losses Turbine 3.58%
Balance of plant 0.5%
Grid 0.5%
Electric losses Efficiency 3%
Turbine performance losses Sub-optimal performance 1.1%
Generic power curve adj. 1.7%
Site-specific power curve adj. 0.81%
High wind hysteresis 0.4%
Environmental losses Degradation 0.1%
Environmental 0.4%
For PV generation, we import the PVWatt module from the NREL’s PySAM package51 to calculate the reference electricity generation of a 1MW PV field for the selected locations. The system’s parameters are provided in Table 4. For wind power generation, we select the NREL’s 2020 ATB Reference 4MW
28
wind turbine with the system parameters provided in Table 5. The reference electricity generation of a
320MW wind farm is calculated using NREL PySAM’s windpower module51.
Table 6 RE capacity factor of five locations in 2021 and 2023.
Location Year RE capacity factor
Gladstone (QLD) 2021 0.28
2023 0.29
Spencer Gulf (SA) 2021 0.35
2023 0.34
Burnie (TAS) 2021 0.55
2023 0.54
Geelong (VIC) 2021 0.42
2023 0.39
Newcastle (NSW) 2021 0.24
2023 0.23
Table 6 provides RE capacity factor of five locations for the years 2021 and 2023. The RE capacity factor ( ) of each location is calculated using the optimized off-grid system’s configuration to evaluate the local RE resources as shown in the following equation:
∑8759
=0 ( ()+ ())
= (25)
( + )×8760
() ()
Where and are the wind and PV generation at each time , and and represent the capacity of wind and PV, respectively.
29
Burnie (TAS) Newcastle (NSW) Geelong (VIC)
0.75 0.35 0.60
Monthly CF - 2021 Monthly CF - 2021 Monthly CF - 2021
0.70 Monthly CF - 2023 Monthly CF - 2023 0.55 Monthly CF - 2023
0.30
0.65 0.50
Capacity factor (wind+PV)
Capacity factor (wind+PV)
Capacity factor (wind+PV)
0.60 0.45
0.25
0.55 0.40
0.20
0.50 0.35
0.45 0.30
0.15
0.40 0.25
0.35 0.10 0.20
Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
Month Month Month
Spencer Gulf (SA) Gladstone (QLD)
0.50 0.45
Monthly CF - 2021 Monthly CF - 2021
Monthly CF - 2023 Monthly CF - 2023
0.45 0.40
Capacity factor (wind+PV)
Capacity factor (wind+PV)
0.40 0.35
0.35 0.30
0.30 0.25
0.25 0.20
0.20 0.15
Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
Month Month
Figure 10 Monthly average capacity factor in five locations in 2021 and 2023.
Figure 10 visualizes the monthly average RE capacity factors for the five cases studied across five states in the NEM, based on the optimized off-grid system configuration. The different coloured curves represent the monthly average RE capacity factors for the year 2021 and 2023.
NEM data
Table 7 Regional grid information of five locations in 2021 and 2023.
Location Year Regional yearly Regional grid RE Yearly average MEF Yearly average AEF of
average electricity penetration of regional grid regional grid
prices (AUD/MWh) (kgCO2e/kWh) (kgCO2e/kWh)
Gladstone (QLD) 2021 89 19% 0.63 0.76
2023 95 28% 0.47 0.69
Spencer Gulf (SA) 2021 59 63% 0.51 0.3
2023 92 71% 0.4 0.23
Burnie (TAS) 2021 35 102% 0.26 0.09
2023 57 93% 0.19 0.12
Geelong (VIC) 2021 48 34% 0.54 0.8
2023 64 42% 0.43 0.73
Newcastle (NSW) 2021 73 23% 0.61 0.74
2023 98 31% 0.52 0.63
30
Table 7 provides regional grid profiles of five locations for the years 2021 and 2023. The data of RE penetration in the regional grid are sourced from Open Electricity’s Tracker26 and yearly electricity prices, MEFs, and AEFs for the regional grid are sourced from the NEMweb dataset, which provides data at 5-minute intervals, and are averaged across the year35.
It is worth noting that the grid profile during 2021 and 2023 for each of the five states reveals significant electricity prices variations across different years. For instance, in 2023, SA yearly average electricity price surged to 92 AUD/MWh, much higher compared to 2021 (59 AUD/MWh). The yearly average
MEFs in five states are all decreasing, and similar decreasing trend can be seen in the yearly average
AEFs due to with higher grid RE penetration in NSW, VIC, SA and QLD except for TAS, where RE penetration is reducing from 102% in 2021 to 93% in 2023 and the yearly average AEFs increase from
0.09 kgCO2e/kWh in 2021 to 0.12 kgCO2e/kWh in 2023.
QLD NSW SA
0.95 0.45 0.75 0.30
0.90
MEF AEF Spot Price MEF AEF Spot Price 0.30 0.70 MEF AEF Spot Price
Emission Factor of Electricity (kgCO2e/kWh)
Emission Factor of Electricity (kgCO2e/kWh)
Emission Factor of Electricity (kgCO2e/kWh)
0.90
0.40 0.85
0.85 0.65 0.25
0.80
0.80 0.35 0.25 0.60
0.75
0.55
Spot Price (AUD/kWh)
Spot Price (AUD/kWh)
Spot Price (AUD/kWh)
0.75
0.30 0.70 0.20
0.70 0.50
0.65 0.20
0.65 0.25 0.45
0.60 0.40 0.15
0.60
0.20 0.55 0.15 0.35
0.55
0.50 0.30
0.50 0.15 0.10
0.45 0.25
0.45 0.10
0.10 0.40 0.20
0.40 0.05
0.35 0.15
0.35 0.05 0.05
0.30 0.10
0.30
0.00 0.25 0.05 0.00
0 2 4 6 8 10 12 14 16 18 20 22 0 2 4 6 8 10 12 14 16 18 20 22 0 2 4 6 8 10 12 14 16 18 20 22
Hour Hour Hour
TAS VIC
0.08 0.20
0.40 0.95 MEF AEF Spot Price
MEF AEF Spot Price
Emission Factor of Electricity (kgCO2e/kWh)
Emission Factor of Electricity (kgCO2e/kWh)
0.90 0.18
0.35 0.07 0.85
0.16
0.80
0.30
Spot Price (AUD/kWh)
Spot Price (AUD/kWh)
0.75 0.14
0.06
0.70
0.25 0.12
0.65
0.20 0.05 0.60 0.10
0.55
0.08
0.15 0.50
0.04
0.45 0.06
0.10 0.40
0.04
0.03 0.35
0.05 0.30 0.02
0.25
0.00 0.02 0.00
0 2 4 6 8 10 12 14 16 18 20 22 0 2 4 6 8 10 12 14 16 18 20 22
Hour Hour
Figure 11 Hourly average of MEFs vs AEFs and spot prices of the regional grid in five states in 2021.
31
QLD NSW SA
0.30
0.90 0.35 0.90 0.60
MEF AEF Spot Price MEF AEF Spot Price 0.30
Emission Factor of Electricity (kgCO2e/kWh)
Emission Factor of Electricity (kgCO2e/kWh)
Emission Factor of Electricity (kgCO2e/kWh)
MEF AEF Spot Price
0.85 0.85
0.55
0.80 0.30 0.80 0.25
0.25 0.50
0.75 0.75
Spot Price (AUD/kWh)
Spot Price (AUD/kWh)
Spot Price (AUD/kWh)
0.25 0.70 0.45
0.70 0.20
0.65 0.20 0.40
0.65 0.20
0.60 0.60 0.35 0.15
0.55 0.15 0.55 0.15 0.30
0.50 0.50
0.25 0.10
0.10 0.45
0.45 0.10 0.20
0.40
0.40 0.05 0.15
0.35 0.05
0.35 0.05
0.00 0.30 0.10
0.30
0.25 0.05 0.00
0 2 4 6 8 10 12 14 16 18 20 22 0 2 4 6 8 10 12 14 16 18 20 22 0 2 4 6 8 10 12 14 16 18 20 22
Hour Hour Hour
TAS VIC
0.35 0.10 0.90 0.20
MEF AEF Spot Price MEF AEF Spot Price
Emission Factor of Electricity (kgCO2e/kWh)
Emission Factor of Electricity (kgCO2e/kWh)
0.85
0.09 0.18
0.30 0.80
0.75 0.16
0.08
Spot Price (AUD/kWh)
Spot Price (AUD/kWh)
0.25 0.70 0.14
0.07 0.65
0.12
0.20 0.60
0.06 0.55 0.10
0.15 0.50 0.08
0.05
0.45
0.06
0.04 0.40
0.10
0.35 0.04
0.03 0.30 0.02
0.05
0.25
0.02 0.00
0 2 4 6 8 10 12 14 16 18 20 22 0 2 4 6 8 10 12 14 16 18 20 22
Hour Hour
Figure 12 Hourly average of MEFs vs AEFs and spot prices of the regional grid in five states in 2023.
We simulate the interaction between the hydrogen production system and the grid via modelling the grid using the historical electricity price, MEF and AEF data published by Australian Energy Market
Operator (AEMO). With the help of the python package NEMOSIS65, we get access to hourly electricity prices of five regions in the NEM. We obtain MEF utilising python package NEMED66 by identifying the marginally dispatched generators from AEMO’s price setter files and associating each marginal generator with its respective carbon dioxide equivalent intensity index published in AEMO’s NEMweb dataset for each time interval35, AEF is computed by the energy sent out consumed from the grid
(including the inter-region electricity trade) and its associated emissions. Emissions data is obtained using NEMED, and data of energy sent out and the inter-region electricity trade is acquired via
NEMOSIS.
Through comparing the hourly average of MEFs vs AEFs in 2021 and 2023 as shown in Figure 11 and
Figure 12, it is noted that the hourly average MEFs are quite different with AEFs. In this work, we adopt the two methods: AEF-based method and MEF-based method, to compare the electricity emissions calculated by two PGO methodologies. The two emissions factors, AEFs and MEFs, answer distinct questions. If we assume that the electricity demand in question represents a portion of the existing demand, then AEFs are more appropriate for answering the question of what share of grid emissions this electricity demand is responsible for. Conversely, if the electricity demand under analysis represents an incremental increase or decrease, MEFs are preferred to determine the impact of the newly changed demand on grid emissions38. This is because changes in demand are not distributed proportionally across all generators within the electricity system. In real operation, specific generators respond to these electricity demand change, and it is the emission intensity of these generators (MEFs)
32
that drives the resulting emission variations39. Thus, we consider the MEF-based method to best represent actual emissions in this work because we use historical grid data to assess emission variations arising from the additional electricity demand of the hybrid grid-connected hydrogen production systems and study how policy interventions like temporal and geographic correlations requirements may impact emissions.
Previous research has emphasised that MEFs are more accurate than AEFs for calculating emissions resulting from shifts in electricity demand and for studying the effects of policy interventions. For example, Sengupta et al. (2022) pointed out MEFs provide accurate assessments of emissions, particularly when analysts or policymakers evaluate interventions that result in changes to electricity demand67. Holland et al. (2022) discovered that using average rather than marginal emissions to predict the impacts significantly overestimates the emission benefits of the policy intervention of Biden administration’s 2030 EV target. They also highlighted it is estimates of marginal emissions that are needed to accurately evaluate the impacts of policies or behaviours that cause changes in the demand or supply of electricity68. Elenes et al. (2022) found that average emission factors have lower accuracy when estimating emissions from demand shifts. In contrast, incremental marginal emission factors
(emissions in response to a small shift in electricity demand) can reproduce the emission changes of a power grid model under many testing conditions and scenarios69. Li et al. (2017) indicated that for evaluating the emission implications of policy and technology interventions, such as electric vehicle tax credits and energy storage solutions, MEFs are a more suitable metric than AEFs70. Siler-Evans et al.
(2012) revealed that AEFs can substantially misestimate the avoided emissions resulting from an intervention, in contrast to MEFs71.
Table 8 Scope 2 emission factors and residual mix factors for consumption of electricity and applicable renewable power percentage.
State Scope 2 emission factor (kgCO2e/kWh) Residual mix factor (kgCO2e/kWh) ARPP
Queensland 0.71 0.81 18.72%
South Australia 0.23 0.81 18.72%
Tasmania 0.15 0.81 18.72%
Victoria 0.77 0.81 18.72%
New South Wales 0.66 0.81 18.72%
Table 8 provides the scope 2 emission factors of five states used in the PGO location-based method, and the residual mix factor and applicable renewable power percentage applied in the PGO market- based method, sourced from the latest National Greenhouse and Energy Reporting (Measurement)
33
Determination 200824. The scope 2 emission factors are state-based, calculated as the average emission factor for all electricity consumed from the grid within a specific state, territory, or electricity grid. The applicable renewable power percentage () is calculated by72,
= ( + )/( − ) (26)
Where is renewable electricity required for the year, is the cumulative adjustment, is the estimated relevant acquisitions of electricity for the year and represents estimated exemptions.
Following the example provided in Australian National Greenhouse Accounts Factors 202443, we adopt the average value for the calendar years 2023 and 2024 in this work72, which is 18.72%. Another key parameter used in the PGO market-based method is the residual mix factor ( ), which is calculated nationally by adjusting the national location-based scope 2 emission factor to exclude the emissions benefits of all claimable renewable generation through LGCs43, calculated by,
= / ( – ( × 1000)) (27)
Where is the residual mix factor for the financial year. represents the emissions from on-grid generation based on NGER reports for the preceding reporting year. is the electricity generation sent out to the grid. represent renewable energy certificates created for renewable energy generation which occurred in the same financial year. The spot price of RECs is estimated to range from 20 AUD to 60 AUD per unit, based on spot price variations observed during 2023-202433.
Generation mix 65.7% 7.3% 20.6%
QLD
Marginal generator 46.7% 8.6% 24.7% 14.2%
Generation mix 23.7% 28.1% 47.4%
SA
Marginal generator 33.0% 12.0% 27.5% 5.7% 15.8%
Generation mix 78.1% 18.1%
TAS
Marginal generator 17.8% 75.6%
Generation mix 60.7% 11.8% 21.2%
VIC
Marginal generator 40.3% 7.4% 31.8% 5.7% 11.2% Coal
Gas
Hydro
Generation mix 63.8% 20.6% 8.4%
Solar
NSW
Marginal generator 50.3% 7.9% 30.2% 5.8%
Wind
Other
0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100%
Figure 13 Electricity generation mix and frequency of marginal generator fuel types across 5-min intervals of five states in 2023
34
Time series illustrating the modelled system operation
Hour
0 2 4 6 8 10 12 14 16 18 20 22
20 Operating power of electrolyser
(I) Electricity import
18
Electricity export
16 Wind generation
Electricity flow (MW)
14 PV generation
12
10
8
6
4
2
0
1.4
1.3 (II) Hydrogen storage level
1.2 Hydrogen storage input
Direct delivery to load
1.1
Hydrogen flow (tonnes)
Hydrogen storage output
1.0
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0.0
0 2 4 6 8 10 12 14 16 18 20 22
Hour
Figure 14 Example of time resolved output of the model over 24 hours in the case of Gladstone (QLD) on 9 Sep
2023. (I) Electricity flow [MW] from different components and (II) hydrogen flow [in tonnes].
Figure 14 gives a representative time series illustrating the modelled system operation on 9 Sep 2023 in Gladstone (QLD). By tracing the electrical flow from various components in Figure 14I and the hydrogen flow in of Figure 14II, we can monitor the operation of the electrolyser, whether power is provided by local RE or the grid, and the hydrogen storage level.
At time 0 (12 am), local RE generation is minimal and the hydrogen storage does not have sufficient inventory to meet the demand. Thus, the system imports electricity from the grid to power the electrolyser, producing hydrogen to meet demand. Once the local RE generation is increasing from 3 am to 5 am, the system gradually reduces its import of grid electricity, instead using local RE to operate the electrolyser. From 7 am to 10 am, as wind and PV generation rises, the system not only utilises local
RE to produce hydrogen for meeting demand and replenishing hydrogen storage but also exports excess electricity to the grid for profit. The system resumes importing the electricity at 11 am to maintain the electrolyser at nominal power as wind generation decreases and continues importing until 16 pm due to low electricity prices during this period. Then, during the high electricity prices of the day in the evening from 6 pm to 11 pm, the system chooses to cease hydrogen production from the electrolyser and export
35
all locally generated electricity to the grid. At this period, the hydrogen stored in hydrogen storage is transported to meet demand.
Limitations
Firstly, our model optimizes the system configuration over a single year with perfect foresight, and we do not investigate how natural variation in weather patterns and the dynamic changes in grid spot price and emissions over the years would affect optimal sizing or operation. Secondly, we use fixed estimates for component costs and do not include location dependent cost factors. Thirdly, we do not consider additional revenue stream for the producer from selling the RECs, due to the high uncertainty of the future fluctuations of the RECs prices73. Instead, we focus on the comparative emissions and costs of systems in different states, under different constraints to assess the accuracy of PGO certification methodologies and the efficacy of different policy settings in minimizing emissions while maintaining low costs.
Additional Results
AEF-based method vs Location-based method
EI_AEF-EI_L (2021) EI_AEF-EI_L (2023)
0.6 0.6
0.5 (a) 0.5
(b)
△Emission intensity (tCO2e/MWhH2)
0.4
△Emission intensity (tCO2e/MWhH2)
0.4
0.3 0.3
0.2 0.2
0.1 0.1
0.0 0.0
−0.1 −0.1
−0.2 −0.2
QLD QLD
−0.3 TAS −0.3
TAS
−0.4 SA −0.4 SA
NSW NSW
−0.5 −0.5
VIC VIC
−0.6 −0.6
0.0 0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4 2.7 3.0 3.3 3.6 0.0 0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4 2.7 3.0 3.3 3.6
RE factor RE factor
Figure 15 (a)-(b) The difference in emission intensities calculated by AEF-based method and PGO location-based method for different RE factors across the five states in 2021-2023.
Figure 15 (a)-(b) shows the differences in emission intensities during 2021 and 2023, calculated by subtracting the PGO location-based method from the AEF-based method. The 2021 results are included here to validate the findings. A good alignment between the PGO location-based method and the AEF- based method can be observed across all systems investigated during 2021 and 2023, with no distinct
36
trends. The differences in emission intensities (EI_AEF - EI_L) remain within a narrow range of -0.15 to 0.1 tCO2e/MWhH2. This demonstrates that PGO location-based method accurately reflects the average GHG emissions associated with electricity consumption from a specific regional grid, and the temporal resolution of the average emission factors does not significantly affect the emissions estimates calculated by this approach.
Table 9 Selected numerical results of Figure 2
Location RE EI_MEF EI_Market EI_RECs EI_L DL DMarket
factor (kgCO2e/kgH2) (kgCO2e/kgH2) (kgCO2e/kgH2) (kgCO2e/kgH2)
Gladstone 0 27.61 37.6 0 40.55 -47% -36%
(QLD) 3.33 -13.62 0 -26.02 -21.41 -57% 91%
Spencer 0 22.55 37.6 0 13.14 42% -67%
Gulf (SA) 2.91 -9.15 0 -20.28 -5.36 41% 122%
Burnie 0 10.06 37.6 0 8.57 15% -274%
(TAS) 1.73 -6.54 0 -22.57 -3.82 42% 245%
Geelong 0 26.9 37.6 0 43.98 -63% -40%
(VIC) 2.52 -9.65 0 -20 -17.54 -82% 107%
Newcastle 0 32.42 37.6 0 37.7 -16% -16%
(NSW) 3.43 -14.07 0 -27.03 -20.59 -46% 92%
The conversion between (kgCO2e/kgH2 and tCO2e/MWhH2 is based on the higher heating value of hydrogen (39.4 kWh/kgH241). The difference percentage between EI_MEF and EI_Market ( ), between EI_MEF and EI_L ( ) is calculated by:
(_−_)
= |_|
(28a)
(_−_)
= |_|
(28b)
(_−_)
= |_|
(29)
When the system generates more RECs that result in the emission intensity calculated by the market- based method falling below zero, EI_Market is assumed to be zero24. Then we use the emission offset by RECs (EI_RECs) as the emission intensity of hydrogen calculated by the market-based method to compare the normalised difference scale with the EI_MEF ( ).
MEF-based method vs PGO methods
37
EI_MEF-EI_Market EI_MEF-EI_L
0.5
0.3
0.4
0.3 0.2
△Emission intensity (tCO2e/MWhH2)
△Emission intensity (tCO2e/MWhH2)
0.2
0.1
0.1
PGO scheme underestimate PGO scheme underestimate
0.0 0.0
overestimate overestimate
−0.1
−0.1
−0.2
−0.3 −0.2
−0.4 QLD QLD
−0.3
TAS TAS
−0.5
SA SA
−0.6 −0.4
NSW NSW
−0.7
VIC VIC
−0.5
0.0 0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4 2.7 3.0 3.3 3.6 0.0 0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4 2.7 3.0 3.3 3.6
RE factor RE factor
Figure 16.The difference in emission intensities calculated by subtracting the PGO market-based method (left) and PGO location-based method (right) from the MEF-based method for different RE factors across the five states in 2023.
B.2 Verification
38
(a) Gladstone (QLD) (b) Newcastle (NSW) (c) Geelong (VIC)
1.2 1.2 1.2
1.0 1.0
Emission intensity (tCO2e/MWhH2)
Emission intensity (tCO2e/MWhH2)
1.0
Emission intensity (tCO2e/MWhH2)
0.8 0.8 0.8
0.6 0.6 0.6
0.4 0.4 0.4
0.2 0.2 0.2
0.0 0.0 0.0
−0.2 −0.2 −0.2
−0.4 −0.4 −0.4
−0.6 −0.6 −0.6
−0.8 −0.8 −0.8
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 0.0 0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4 2.7 3.0 3.3 3.6 3.9 0.0 0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4
RE factor RE factor RE factor
(d) Burnie (TAS) (e) Spencer Gulf (SA)
1.0 1.0
Emission intensity (tCO2e/MWhH2)
Emission intensity (tCO2e/MWhH2)
0.8 0.8
0.6 0.6 EI_RECs
0.4 0.4 EI_Market
0.2 0.2
EI_MEF
EI_L
0.0 0.0
−0.2 −0.2
−0.4 −0.4
−0.6 −0.6
−0.8 −0.8
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 0.0 0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4 2.7 3.0
RE factor RE factor
Figure 17 (a)-(e) GHG emission intensities of hydrogen under various RE factors in 2021, calculated using three different methods for electricity emissions.
Figure 17 (a)-(e) replicate Figure 2, utilizing grid data inputs from 2021. These figures illustrate the
GHG emission intensities of hydrogen under various RE factors, calculated using three different methods for electricity emissions, along with the differences in emission intensities between the market- based method and the MEF-based method.
39
Burnie (TAS) Newcastle (NSW) Geelong (VIC)
7 45 7 45 7 45
6 40 6 40 6 40
Emission intensity (kgCO2e/kg H2)
Emission intensity (kgCO2e/kg H2)
Emission intensity (kgCO2e/kg H2)
35 35 35
5 5 5
30 30 30
4 25 4 25 4 25
LCOH ($/kg H2)
LCOH ($/kg H2)
LCOH ($/kg H2)
3 20 3 20 3 20
15 15 15
2 2 2
10 10 10
1 5 1 5 1 5
0 0 0 0 0 0
−1 −5 −1 −5 −1 −5
−10 −10 −10
−2 −2 −2
−15 −15 −15
−3 −20 −3 −20 −3 −20
ly
ly
ly
y
y
y
e
≤0
le
≤0
le
≤0
id
y
ly
id
y
ly
id
y
ly
bl
rl
rl
rl
nl
nl
nl
ib
ib
ai
ai
ai
th
th
th
gr
gr
gr
EF
EF
EF
ea
ea
ea
i
o
o
o
D
D
D
ex
ex
ex
on
on
on
ff-
ff-
ff-
ll
ll
ll
Y
Y
Y
_M
_M
_M
Fl
Fl
Fl
Se
Se
Se
M
M
M
O
O
O
EI
EI
EI
Spencer Gulf (SA) Gladstone (QLD)
7 45 7 45
6 40 6 40
Emission intensity (kgCO2e/kg H2)
Emission intensity (kgCO2e/kg H2)
35 35
5 5 LCOH
30 30
4 25 4 25
EI_Market
LCOH ($/kg H2)
LCOH ($/kg H2)
3 20 3 20 EI_MEF
2
15
2
15 Grid electricity cost
10 10
1 1
Electrolyser CAPEX+OM
5 5
0 0 0 0
Hydrogen Storage CAPEX
−1 −5 −1 −5 PV CAPEX+OM
−2
−10
−2
−10 Wind CAPEX+OM
−15 −15
−3 −20 −3 −20
ly
ly
y
y
le
≤0
le
≤0
id
y
ly
id
y
ly
rl
rl
nl
nl
ib
ib
ai
ai
th
th
gr
gr
EF
EF
ea
ea
o
o
D
D
ex
ex
on
on
ff-
ff-
ll
ll
Y
Y
_M
_M
Fl
Fl
Se
Se
M
M
O
O
EI
EI
Figure 18 The LCOH (white dot points) and yearly average emission intensities (red and yellow triangles) of hydrogen calculated for five states in 2021 under different configurations which have been optimized under different scenarios to minimize LCOH.
Figure 18 depicts the LCOH (white dot points) and yearly average emission intensities (red and yellow triangles) of hydrogen for five states calculated under different configurations optimized to minimize
LCOH. Figure 18 replicates Figure 3 using the grid data in 2021 for optimization. It is worth noting that the case of Spencer Gulf (SA) in 2021 exhibits different optimization patterns compared to the system optimized using 2023 data, as shown in Figure 3. Although the RE capacity factor in Spencer
Gulf (SA) is higher in 2021 (0.35) than in 2023 (0.34), indicating more favourable renewable resources, the significantly lower electricity prices in 2021 (59 AUD/MWh) compared to 2023 (92 AUD/MWh) result in the optimized systems under different scenarios in Spencer Gulf (SA) following similar patterns to those observed in TAS, VIC, and NSW. Although the yearly average RE capacity factors of these five cases do not vary significantly (as shown in Table 6), the differences in the optimized system results between 2023 (as presented in the Figure 3) and 2021 are highly related to variations in electricity prices.
Conflicts of interest
There are no conflicts to declare.
Data availability
The data supporting this article have been included in the manuscript.
40
Author CRediT statement
Chengzhe Li: Conceptualisation, Investigation, Methodology, Visualization, Writing-Original Draft,
Writing-Review and Editing. Lee V. White: Conceptualisation, Investigation, Writing-Review and
Editing. Reza Fazeli: Conceptualisation, Investigation, Writing-Review and Editing. Anna Skobeleva:
Methodology, Writing-Review and Editing, Michael Thomas: Methodology, Writing-Review and
Editing, Shuang Wang: Methodology, Writing-Review and Editing. Fiona J Beck: Conceptualisation,
Methodology, Writing-Original Draft, Writing-Review and Editing, Supervision.
Acknowledgements
This work has been supported by the Heavy Industry Low-carbon Transition Cooperative Research
Centre (HILT CRC) whose activities are funded by industry, research and government partners along with the Australian Government’s Cooperative Research Centre Program. C.L. acknowledges support from the ANU Centre for Energy Systems.
References
1. Longden, T., Beck, F. J., Jotzo, F., Andrews, R. & Prasad, M. ‘Clean’ hydrogen? – Comparing the
emissions and costs of fossil fuel versus renewable electricity based hydrogen. Applied Energy 306,
118145 (2022).
2. Vallejo, V., Nguyen, Q. & Ravikumar, A. P. Geospatial variation in carbon accounting of hydrogen
production and implications for the US Inflation Reduction Act. Nat Energy (2024)
doi:10.1038/s41560-024-01653-0.
3. White, L. V. et al. Towards emissions certification systems for international trade in hydrogen: The
policy challenge of defining boundaries for emissions accounting. Energy 215, 119139 (2021).
4. Goodwin, D. et al. Sustainability certification for renewable hydrogen: An international survey of
energy professionals. Energy Policy 192, 114231 (2024).
5. Köveker, T. et al. Green premiums are a challenge and an opportunity for climate policy design.
Nat. Clim. Chang. 13, 592–595 (2023).
41
6. White, L. V., Aisbett, E., Pearce, O. & Cheng, W. Principles for embedded emissions accounting to
support trade-related climate policy. Climate Policy 1–17 (2024)
doi:10.1080/14693062.2024.2356803.
7. Giovanniello, M. A., Cybulsky, A. N., Schittekatte, T. & Mallapragada, D. S. The influence of
additionality and time-matching requirements on the emissions from grid-connected hydrogen
production. Nat Energy 9, 197–207 (2024).
8. Engstam, L., Janke, L., Sundberg, C. & Nordberg, Å. Grid-supported electrolytic hydrogen
production: Cost and climate impact using dynamic emission factors. Energy Conversion and
Management 293, 117458 (2023).
9. European Commission. Directive (EU) 2018/2001 of the European Parliament and of the Council
of 11 December 2018 on the Promotion of the Use of Energy from Renewable Sources (Recast)
(Text with EEA Relevance). https://eur-lex.europa.eu/legal-
content/EN/TXT/PDF/?uri=CELEX:02018L2001-20240716 (2024).
10. Commission Delegated Regulation (EU) 2023/1185 of 10 February 2023 Supplementing Directive
(EU) 2018/2001 of the European Parliament and of the Council by Establishing a Minimum
Threshold for Greenhouse Gas Emissions Savings of Recycled Carbon Fuels and by Specifying a
Methodology for Assessing Greenhouse Gas Emissions Savings from Renewable Liquid and
Gaseous Transport Fuels of Non-Biological Origin and from Recycled Carbon Fuels. 507–531
https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:32023R1185 (2023)
doi:10.1017/CBO9780511610851.034.
11. European Commission. Commission Delegated Regulation (EU) 2023/1184 of 10 February 2023
Supplementing Directive (EU) 2018/2001 of the European Parliament and of the Council by
Establishing a Union Methodology Setting out Detailed Rules for the Production of Renewable
Liquid and Gaseous Transport Fuels of Non-Biological Origin. https://eur-lex.europa.eu/legal-
content/EN/TXT/PDF/?uri=CELEX:02023R1184-20240610 (2023).
12. CERTIFHY. CERTIFHY EU RFNBO VOLUNTARY SCHEME. https://www.certifhy.eu/wp-
content/uploads/2024/09/CertifHy_Scheme_doc_Scheme_document.pdf (2024).
42
13. Brandt, J. et al. Cost and competitiveness of green hydrogen and the effects of the European Union
regulatory framework. Nat Energy 9, 703–713 (2024).
14. Ruhnau, O. & Schiele, J. Flexible green hydrogen: The effect of relaxing simultaneity requirements
on project design, economics, and power sector emissions. Energy Policy 182, 113763 (2023).
15. Schlund, D. & Theile, P. Simultaneity of green energy and hydrogen production: Analysing the
dispatch of a grid-connected electrolyser. Energy Policy 166, 113008 (2022).
16. Internal Revenue Service Treasury. Section 45V Credit for Production of Clean Hydrogen; Section
48(a)(15) Election To Treat Clean Hydrogen Production Facilities as Energy Property.
https://www.govinfo.gov/content/pkg/FR-2023-12-26/pdf/2023-28359.pdf (2023).
17. DCCEEW. Future Made in Australia (Guarantee of Origin) Act 2024.
https://www.legislation.gov.au/C2024A00121/asmade/2024-12-10/text/original/pdf (2024).
18. DCCEEW. National Hydrogen Strategy 2024.
https://www.dcceew.gov.au/sites/default/files/documents/national-hydrogen-strategy-2024.pdf
(2024).
19. The Treasury. Establishing a ‘Front Door’ for Major, Transformational Projects - Consultation
Paper. https://treasury.gov.au/sites/default/files/2024-09/c2024-571335-consultation-paper.pdf
(2024).
20. DCCEEW. Guarantee of Origin - Emissions Accounting Approach Paper.
https://storage.googleapis.com/files-au-climate/climate-
au/p/prj291cc9979281a4ffc59d8/public_assets/Guarantee%20of%20Origin%20-
%20Emissions%20Accounting%20Approach%20paper.pdf (2023).
21. DCCEEW. Guarantee of Origin Scheme Design Paper. https://storage.googleapis.com/files-au-
climate/climate-
au/p/prj291cc9979281a4ffc59d8/public_assets/Guarantee%20of%20Origin%20Scheme%20desig
n%20paper.pdf (2023).
22. DCCEEW. Guarantee of Origin scheme. https://www.dcceew.gov.au/energy/renewable/guarantee-
of-origin-scheme (2025).
43
23. International Partnership for Hydrogen and Fuel Cells in the Economy. Methodology for
Determining the Greenhouse Gas Emissions Associated with the Production of Hydrogen.
https://www.iphe.net/_files/ugd/45185a_8f9608847cbe46c88c319a75bb85f436.pdf (2023).
24. National Greenhouse and Energy Reporting Scheme. National Greenhouse and Energy Reporting
(Measurement) Determination 2008. https://www.legislation.gov.au/F2008L02309/latest/text
(2024).
25. AEMO. national-electricity-market-fact-sheet. https://aemo.com.au/-
/media/files/electricity/nem/national-electricity-market-fact-sheet.pdf (2024).
26. The Superpower Institute. Open Electricity.
https://explore.openelectricity.org.au/energy/nem/?range=7d&interval=30m&view=discrete-
time&group=Detailed (2025).
27. Giles Parkinson. From zero to 100 pct renewables in just 20 years: South Australia’s remarkable
energy transition. RenewEconomy https://reneweconomy.com.au/from-zero-to-100-pct-
renewables-in-just-20-years-south-australias-remarkable-energy-transition/ (2024).
28. Government of South Australia. New target for renewables. Government of South Australia
https://www.premier.sa.gov.au/media-releases/news-archive/new-target-for-renewables (2024).
29. Queensland Government. Queensland’s renewable energy targets.
https://www.energyandclimate.qld.gov.au/energy/reasons-for-renewable/renewable-energy-
targets.
30. Langenmayr, U. & Ruppert, M. Renewable origin, additionality, temporal and geographical
correlation – eFuels production in Germany under the RED II regime. Energy Policy 183, 113830
(2023).
31. Frazer-Nash Consultancy. Fugitive Hydrogen Emissions in a Future Hydrogen Economy.
https://assets.publishing.service.gov.uk/media/624ec79cd3bf7f600d4055d1/fugitive-hydrogen-
emissions-future-hydrogen-economy.pdf (2022).
32. Mojiri, A., Wang, Y., Rahbari, A., Pye, J. & Coventry, J. Current and future cost of large-scale green
hydrogen generation. SSRN Scholarly Paper at https://doi.org/10.2139/ssrn.4477925 (2023).
44
33. Clean Energy Regulator. Large-scale generation certificate (LGC) reported spot and forward prices.
https://cer.gov.au/markets/reports-and-data/quarterly-carbon-market-reports/quarterly-carbon-
market-report-december-quarter-2024/large-scale-generation-certificates-lgcs (2024).
34. Bracci, J. M., Sherwin, E. D., Boness, N. L. & Brandt, A. R. A cost comparison of various hourly-
reliable and net-zero hydrogen production pathways in the United States. Nat Commun 14, 7391
(2023).
35. AEMO. NEMweb market data. https://aemo.com.au/energy-systems/electricity/national-
electricity-market-nem/data-nem/market-data-nemweb (2024).
36. TasNetworks. Prescribed Transmission Service Prices for 2024-25.
https://www.tasnetworks.com.au/config/getattachment/96b170ef-8d44-43b9-97bf-
7ba6b2304bd0/2024-25-prescribed-transmission-service-prices.pdf.
37. Salmon, N. & Bañares-Alcántara, R. Impact of grid connectivity on cost and location of green
ammonia production: Australia as a case study. Energy Environ. Sci. 14, 6655–6671 (2021).
38. Ryan, N. A., Johnson, J. X. & Keoleian, G. A. Comparative Assessment of Models and Methods
To Calculate Grid Electricity Emissions. Environ. Sci. Technol. 50, 8937–8953 (2016).
39. Hawkes, A. D. Estimating marginal CO2 emissions rates for national electricity systems. Energy
Policy 38, 5977–5987 (2010).
40. Elenes, A. G. N., Williams, E., Hittinger, E. & Goteti, N. S. How Well Do Emission Factors
Approximate Emission Changes from Electricity System Models? Environ. Sci. Technol. 56,
14701–14712 (2022).
41. Harrison, K. W., Remick, R., Martin, G. D. & Hoskin, A. Hydrogen Production: Fundamentals and
Case Study Summaries. https://www.nrel.gov/docs/fy10osti/47302.pdf (2010).
42. Schäfer, M., Cerdas, F. & Herrmann, C. Towards standardized grid emission factors:
methodological insights and best practices. Energy Environ. Sci. 17, 2776–2786 (2024).
43. DCCEEW. Australian National Greenhouse Accounts Factors.
https://www.dcceew.gov.au/sites/default/files/documents/national-greenhouse-account-factors-
2024.pdf (2024).
45
44. Eric, L. et al. Comparison of Commercial, State-of-the-Art, Fossil-Based Hydrogen Production
Technologies. https://doi.org/10.2172/1862910 (2022).
45. CSIRO. Hunter Valley Hydrogen Hub (Updated Project Description). HyResource
https://research.csiro.au/hyresource/hunter-valley-hydrogen-hub/ (2024).
46. DCCEEW. Renewable Energy (Electricity) Regulations 2001.
https://www.legislation.gov.au/F2001B00053/2024-10-14/2024-10-14/text/original/pdf (2024).
47. DCCEEW. Renewable Energy (Electricity) Act 2000.
https://www.legislation.gov.au/C2004A00767/2025-01-01/2025-01-01/text/original/pdf.
48. Clean Energy Regulator. Renewable Energy Target. https://cer.gov.au/schemes/renewable-energy-
target (2025).
49. Rees S., Grubnik P., Easton L., & Feitz A. Australian Hydrogen Projects Dataset (June 2024).
http://pid.geoscience.gov.au/dataset/ga/146425 (2024).
50. Australian Energy Regulator. Annual electricity consumption - NEM.
https://www.aer.gov.au/industry/registers/charts/annual-electricity-consumption-nem (2024).
51. PySAM Version 2.2.0. National Renewable Energy Laboratory (2024).
52. Rienecker, M. M. et al. MERRA: NASA’s Modern-Era Retrospective Analysis for Research and
Applications. Journal of Climate 24, 3624–3648 (2011).
53. AEMO. 2023 Costs and Technical Parameter Review AEMO. https://aemo.com.au/-
/media/files/stakeholder_consultation/consultations/nem-consultations/2023/2024-forecasting-
assumptions-update-consultation-page/aurecon---2023-cost-and-technical-parameters-review.pdf
(2023).
54. Freeman, J., Jorgenson, J., Gilman, P. & Ferguson, T. Reference Manual for the System Advisor
Model’s Wind Power Performance Model. NREL/TP-6A20-60570, 1150800
http://www.osti.gov/servlets/purl/1150800/ (2014) doi:10.2172/1150800.
55. Umlaut & Agora Industry. Levelised Cost of Hydrogen: Making the Application of the LCOH
Concept More Consistent and More Useful. https://www.agora-
46
industry.org/fileadmin/Projekte/2022/2022-12-10_Trans4Real/A-EW_301_LCOH_WEB.pdf
(2023).
56. Australian Energy Regulator. State of the Energy Market 2021.
https://www.aer.gov.au/system/files/State%20of%20the%20energy%20market%202021%20-
%20Chapter%202%20-%20National%20Electricity%20Market_0.pdf (2021).
57. ARENA. Corporate Renewable Power Purchase Agreements in Australia – State of the Market
2019. https://arena.gov.au/assets/2018/10/corporate-renewable-power-purchase-agreements-in-
australia-state-of-the-market-2019.pdf (2019).
58. Clean Energy Regulator. Designing the Guarantee of Origin. https://cer.gov.au/schemes/guarantee-
origin-scheme/designing-guarantee-origin (2024).
59. Geoscience Australia. Electricity Transmission Lines. https://digital.atlas.gov.au/maps/electricity-
transmission-lines/explore (2024).
60. Graham, P., Hayward, J, & Foster J. GenCost 2023‐24: Final Report. https://www.csiro.au/-
/media/Energy/GenCost/GenCost2023-24Final_20240522.pdf (2024).
61. Rosenstiel, A. et al. Electrochemical Hydrogen Production Powered by PV/CSP Hybrid Power
Plants: A Modelling Approach for Cost Optimal System Design. Energies 14, 3437 (2021).
62. Terlouw, T., Bauer, C., McKenna, R. & Mazzotti, M. Large-scale hydrogen production via water
electrolysis: a techno-economic and environmental assessment. Energy Environ. Sci. 15, 3583–
3602 (2022).
63. Salmon, N. & Bañares-Alcántara, R. Impact of grid connectivity on cost and location of green
ammonia production: Australia as a case study. Energy Environ. Sci. 14, 6655–6671 (2021).
64. Papadias, D. D. & Ahluwalia, R. K. Bulk storage of hydrogen. International Journal of Hydrogen
Energy 46, 34527–34541 (2021).
65. Gorman, N. NEMOSIS, a program for accessing and compiling data on the Australian National
Electricity Market. https://github.com/UNSW-CEEM/NEMOSIS (2018).
66. Heim D & Naderi S. NEMED. https://github.com/UNSW-CEEM/NEMED (2022).
47
67. Sengupta, S. et al. Current and Future Estimates of Marginal Emission Factors for Indian Power
Generation. Environ. Sci. Technol. 56, 9237–9250 (2022).
68. Holland, S. P., Kotchen, M. J., Mansur, E. T. & Yates, A. J. Why marginal CO 2 emissions are not
decreasing for US electricity: Estimates and implications for climate policy. Proc. Natl. Acad. Sci.
U.S.A. 119, e2116632119 (2022).
69. Elenes, A. G. N., Williams, E., Hittinger, E. & Goteti, N. S. How Well Do Emission Factors
Approximate Emission Changes from Electricity System Models? Environ. Sci. Technol. 56,
14701–14712 (2022).
70. Li, M., Smith, T. M., Yang, Y. & Wilson, E. J. Marginal Emission Factors Considering Renewables:
A Case Study of the U.S. Midcontinent Independent System Operator (MISO) System. Environ.
Sci. Technol. 51, 11215–11223 (2017).
71. Siler-Evans, K., Azevedo, I. L. & Morgan, M. G. Marginal Emissions Factors for the U.S.
Electricity System. Environ. Sci. Technol. 46, 4742–4748 (2012).
72. Clean Energy Regulator. Renewable power percentage. https://cer.gov.au/schemes/renewable-
energy-target/renewable-energy-target-liability-and-exemptions/renewable-power-percentage
(2024).
73. Energy Action. Understanding LGC Prices for Energy Investments.
https://energyaction.com.au/understanding-lgc-prices-energy-investments/ (2024).
48