Abstract Magnetotelluric (MT) impedances from 62 sites in southern South Island of Aotearoa New Zealand have been used to model geomagnetically induced currents (GIC) in four transformers during two solar storms. Induced electric fields during the storms are calculated from the MT impedances using the magnetic fields measured at the Eyrewell (EYR) geomagnetic observatory, approximately 200 km north of the study area. Calculated GIC during the sudden storm commencements (SSC) give a generally good match to GIC measured by the network operator, Transpower New Zealand. Long period GIC (periods longer than about 10,000 s) are less well modeled. Calculations based on thin‐sheet modeling, which has restrictions on the shortest period of variation which can be modeled, perform less well for the GIC associated with SSC, but are equally good, if not better, at modeling longer period GIC. Consistent underestimation of large GIC at one transformer (HWBT4) near Dunedin are likely to be the result of uncertainty in the assumed values of line, transformer, and earthing resistances. The assumption of a spatially uniform magnetic field across the study area, which is implied by use of the magnetic field measured at EYR as a basis for calculation, may also lead to incorrect calculation of GIC. For one storm use of magnetic field data from a magnetometer within the study area leads to much improved modeling of the observed GIC. This study compares modeled and measured GIC using specifically measured MT impedance data.

Abstract Reducing the impact of Geomagnetically induced currents (GICs) on electrical power networks is an essential step to protect network assets and maintain reliable power transmission during and after storm events. In this study, multiple mitigation strategies are tested during worst‐case extreme storm scenarios in order to investigate their effectiveness for the New Zealand transmission network. By working directly with our industry partners, Transpower New Zealand Ltd, a mitigation strategy in the form of targeted line disconnections has been developed. This mitigation strategy proved more effective than previous strategies at reducing GIC magnitudes and durations at transformers at most risk to GIC while still maintaining the continuous supply of power throughout New Zealand. Under this mitigation plan, the average 60‐min mean GIC decreased for 27 of the top 30 at‐risk transformers, and the total network GIC was reduced by 16%. This updated mitigation has been adopted as an operational procedure in the New Zealand national control room to manage GIC. In addition, simulations show that the installation of 14 capacitor blocking devices at specific transformers reduces the total GIC sum in the network by an additional 16%. As a result of this study Transpower is considering further mitigation in the form of capacitor blockers. We strongly recommend collaborating with the relevant power network providers to develop effective mitigation strategies that reduce GIC and have a minimal impact on power distribution.

AbstractGas pipelines can experience elevated pipe to soil potentials (PSPs) during geomagnetic disturbances due to the induced geoelectric field. Gas pipeline operators use cathodic protection to keep PSPs between −0.85 and −1.2 V to prevent corrosion of the steel pipes and disbondment of the protective coating from the pipes. We have developed a model of the gas pipelines in the North Island of New Zealand to identify whether a hazard exists to these pipelines and how big this hazard is. We used a transmission line representation to model the pipelines and a nodal admittance matrix method to calculate the PSPs at nodes up to 5 km apart along the pipelines. We used this model to calculate PSPs resulting from an idealized 100 mVkm−1 electric field, initially to the north and east. The calculated PSPs are highest are at the ends of the pipelines in the direction of the applied electric field vector. The calculated PSP follows a characteristic curve along the length of the pipelines that matches theory, with deviations due to branchlines and changes in pipeline direction. The modeling shows that the PSP magnitudes are sensitive to the branchline coating conductance with higher coating conductances decreasing the PSPs at most locations. Enhanced PSPs produce the highest risk of disbondment and corrosion occurring, and hence this modeling provides insights into the network locations most at risk.

Geomagnetically induced currents (GICs) produced during geomagnetic disturbances pose a risk to the safe operation of electrical power networks. One route to determine the hazard of large and extreme geomagnetic disturbances to national electrical networks is a validated model to predict GIC across the entire network. In this study, we improve upon an earlier model for New Zealand, expanding the approach to cover transformers nationwide by making use of multiple storms to develop national scaling factors. We exploit GIC observations which have been made and archived by Transpower New Zealand Ltd, the national grid operator. For some transformers the GIC observations span nearly 2 decades, while for others only a few years. GICs can vary wildly between transformers, particularly due to differences in the electrical network characteristics, transformer properties, and ground conductivity. Modeling these individual transformers is required if an accurate representation of the GIC distribution throughout the network is to be produced. Here we model the GIC during 25 disturbed periods, ranging from large geomagnetic storms to weakly active periods. We adopt the approach of scaling model output using observed GIC power spectra, finding that it improves the correlations between the maximum model and observed GIC by between 10% and 40% depending on the transformer. The modeled GIC at the 73 transformers which have measured GIC are analyzed to create local and national scaling curves. These are used to allow modeling for transformers without in‐situ GIC. We present approaches to utilize this technique for future storms, including non‐monitored transformers.

Abstract Geomagnetically induced currents produced during geomagnetic storms present a potential risk to gas pipeline systems by disturbing the cathodic protection (CP) systems used to protect the pipeline from corrosion. In this study we look at CP monitoring data from a number of sites on the gas pipeline network in the North Island of New Zealand. We focus on variations during geomagnetic storms in three aspects of the CP system: (a) the output voltage of the constant current rectifiers, (b) the potential between the pipe and a Cu/CuSO4 reference cell, and (c) the potential between an installed metal coupon and the reference cell. The industry standard for suitable CP is that the latter potential should lie between −0.85 and −1.2 V when the rectifier is momentarily turned off. Three monitoring sites in particular are identified as showing very large variations in rectifier output voltage during geomagnetic storms, suggesting a possible risk to the system at these sites. Additionally, one of the three sites has coupon potentials which appear to rise above the −0.85 V level. Magnetotelluric impedance data are used to assess how effects at the monitoring sites relate to the magnitude and direction of induced electric fields. It clear that the response at any monitoring site is related to effects on the pipeline network as a whole. We also note the risks of disbonding and hydrogen induced cracking, issues we do not believe have been widely recognized in the research community.

Geomagnetically induced currents (GICs) in the North Island New Zealand power transmission network during two large magnetic storms are calculated from both magnetotelluric (MT) data and a thin‐sheet conductance model of New Zealand previously used to study GIC in the South Island. We focus on the 2015 St. Patrick's Day magnetic storm and the storm of 20 November 2003. Lack of MT data in the northwestern part of the Island means that the transmission network in this region is represented by an equivalent circuit. Lack of GIC observations in the North Island means that results cannot be directly compared with measured GIC. However, our calculation of GIC shows that substations and individual transformers in the lower part of the Island with significant currents are generally the same as those where total harmonic distortion has been observed during periods of enhanced geomagnetic activity. MT data in the period range 2–30 min are used to predict GIC associated with the sudden storm commencement and rapid variations in the magnetic field. In contrast, the thin‐sheet modeling approach shows that GIC may be expected to occur in conjunction with longer‐period variations. Calculations for the 2003 storm suggest that at some locations GIC in excess of 10 A may persist for long periods of time and may produce significant harmonic distortion which could lead to localized transformer heating. It is concluded that despite its relatively low latitude the North Island power network is potentially at risk from significant GIC during extreme storms.

Abstract Geomagnetically induced currents (GICs) during a space weather event have previously caused transformer damage in New Zealand. During the 2015 St. Patrick's Day Storm, Transpower NZ Ltd has reliable GIC measurements at 23 different transformers across New Zealand's South Island. These observed GICs show large variability, spatially and within a substation. We compare these GICs with those calculated from a modeled geolectric field using a network model of the transmission network with industry‐provided line, earthing, and transformer resistances. We calculate the modeled geoelectric field from the spectra of magnetic field variations interpolated from measurements during this storm and ground conductance using a thin‐sheet model. Modeled and observed GIC spectra are similar, and coherence exceeds the 95% confidence threshold, for most valid frequencies at 18 of the 23 transformers. Sensitivity analysis shows that modeled GICs are most sensitive to variation in magnetic field input, followed by the variation in land conductivity. The assumption that transmission lines follow straight lines or getting the network resistances exactly right is less significant. Comparing modeled and measured GIC time series highlights that this modeling approach is useful for reconstructing the timing, duration, and relative magnitude of GIC peaks during sudden commencement and substorms. However, the model significantly underestimates the magnitude of these peaks, even for a transformer with good spectral match. This is because of the limited range of frequencies for which the thin‐sheet model is valid and severely limits the usefulness of this modeling approach for accurate prediction of peak GICs.

AbstractLarge geomagnetic storms are a known space weather hazard to power transmission networks due to the effects of geomagnetically induced currents (GICs). However, research in this area has been hampered by a lack of GIC observations. Previous studies have noted that New Zealand is unusually fortunate in having a comparatively dense, high quality, set of GIC measurements, spanning >60 transformers in >20 substations. However, due to operational reasons these observations are clustered in the mid and lower South Island. In this paper we analyze space weather‐induced GIC impact patterns over the entire country by using a different set of sensors that monitor levels of harmonic distortion, with even and odd harmonics measured separately. GICs lead to half cycle transformer saturation and is one of the few ways in which even harmonics are produced in a well‐run power transmission network. We make use of harmonic distortion measurements at 377 circuit breakers made at 126 separate locations. Focusing on the intense geomagnetic storm activity during 6–9 September 2017, we show how the even harmonic distortion observations provide a useful new picture of GIC‐stressed transformers. These observations demonstrate how GIC effects can be monitored by using even harmonic distortion in locations where no GIC measurements are present (e.g., the most of the North Island). We understand that harmonic distortion measurements are fairly common in electrical networks and could provide a new tool for space weather researchers.

17 During space weather events, geomagnetically induced currents (GICs) can be induced 18 in high voltage transmission networks, damaging individual transformers within substations. 19 A common approach to modeling a transmission network has been to assume that every sub20 station can be represented by a single resistance to Earth. We have extended that model by 21 building a transformer-level network representation of New Zealand’s South Island trans22 mission network. We represent every transformer winding at each earthed substation in the 23 network by its known DC resistance. Using this network representation significantly changes 24 the GIC hazard assessment, compared to assessments based on the earlier assumption. Fur25 ther, we have calculated the GIC flowing through a single phase of every individual trans26 former winding in the network. These transformer-level GIC calculations show variation in 27 GICs between transformers within a substation due to transformer characteristics and con28 nections. The transformer-level GIC calculations alter the hazard assessment by up to an 29 order of magnitude in some places. In most cases the calculated GIC variations match mea30 sured variations in GIC flowing through the same transformers. This comparison with an 31 extensive set of observations demonstrates the importance of transformer-level GIC calcula32 tions in models used for hazard assessment . 33

During space weather events, geomagnetically induced currents (GICs) can be induced in high‐voltage transmission networks, damaging individual transformers within substations. A common approach to modeling a transmission network has been to assume that every substation can be represented by a single resistance to Earth. We have extended that model by building a transformer‐level network representation of New Zealand's South Island transmission network. We represent every transformer winding at each earthed substation in the network by its known direct current resistance. Using this network representation significantly changes the GIC hazard assessment, compared to assessments based on the earlier assumption. Further, we have calculated the GIC flowing through a single phase of every individual transformer winding in the network. These transformer‐level GIC calculations show variation in GICs between transformers within a substation due to transformer characteristics and connections. The transformer‐level GIC calculations alter the hazard assessment by up to an order of magnitude in some places. In most cases the calculated GIC variations match measured variations in GIC flowing through the same transformers. This comparison with an extensive set of observations demonstrates the importance of transformer‐level GIC calculations in models used for hazard assessment.

Transfer functions are calculated for periods between 2 and 1,000 min between geomagnetically induced currents (GIC) measured at three transformers in the South Island of New Zealand and variations in the horizontal components of the geomagnetic field measured at the Eyrewell Observatory near Christchurch. Using an inverse Fourier transform, the transfer functions allow the GIC expected in these transformers to be estimated for any variation of the inducing magnetic field. Comparison of the predicted GIC with measured GIC for individual geomagnetic storms shows remarkable agreement, although the lack of high‐frequency measurements of GIC and the need for interpolation of the measurements lead to a degree of underestimation of the peak GIC magnitude. An approximate correction for this is suggested. Calculation of the GIC for a magnetic storm in November 2001 that led to the failure of a transformer in Dunedin suggests that peak GIC were as large as about 80 A. Use of spectral scaling to estimate the likely GIC associated with a geomagnetic storm of the magnitude of the 1859 Carrington Event indicates that GIC of at least 10 times this magnitude may occur at some locations. Although the impact of changes to the transmission network on calculated transfer functions remains to be explored, it is suggested that the use of this technique may provide a useful check on estimates of GIC produced by other methods such as thin sheet modeling.

Geomagnetically induced current (GIC) observations made in New Zealand over 14 years show induction effects associated with a rapidly varying horizontal magnetic field (dBH/dt) during geomagnetic storms. This study analyzes the GIC observations in order to estimate the impact of extreme storms as a hazard to the power system in New Zealand. Analysis is undertaken of GIC in transformer number six in Islington, Christchurch (ISL M6), which had the highest observed currents during the 6 November 2001 storm. Using previously published values of 3,000 nT/min as a representation of an extreme storm with 100 year return period, induced currents of ~455 A were estimated for Islington (with the 95% confidence interval range being ~155–605 A). For 200 year return periods using 5,000 nT/min, current estimates reach ~755 A (confidence interval range 155–910 A). GIC measurements from the much shorter data set collected at transformer number 4 in Halfway Bush, Dunedin, (HWB T4), found induced currents to be consistently a factor of 3 higher than at Islington, suggesting equivalent extreme storm effects of ~460–1,815 A (100 year return) and ~460–2,720 A (200 year return). An estimate was undertaken of likely failure levels for single‐phase transformers, such as HWB T4 when it failed during the 6 November 2001 geomagnetic storm, identifying that induced currents of ~100 A can put such transformer types at risk of damage. Detailed modeling of the New Zealand power system is therefore required to put this regional analysis into a global context.

Transformers in New Zealand's South Island electrical transmission network have been impacted by geomagnetically induced currents (GIC) during geomagnetic storms. We explore the impact of GIC on this network by developing a thin‐sheet conductance (TSC) model for the region, a geoelectric field model, and a GIC network model. (The TSC is composed of a thin‐sheet conductance map with underlying layered resistivity structure.) Using modeling approaches that have been successfully used in the United Kingdom and Ireland, we applied a thin‐sheet model to calculate the electric field as a function of magnetic field and ground conductance. We developed a TSC model based on magnetotelluric surveys, geology, and bathymetry, modified to account for offshore sediments. Using this representation, the thin sheet model gave good agreement with measured impedance vectors. Driven by a spatially uniform magnetic field variation, the thin‐sheet model results in electric fields dominated by the ocean‐land boundary with effects due to the deep ocean and steep terrain. There is a strong tendency for the electric field to align northwest‐southeast, irrespective of the direction of the magnetic field. Applying this electric field to a GIC network model, we show that modeled GIC are dominated by northwest‐southeast transmission lines rather than east‐west lines usually assumed to dominate.

Transpower New Zealand Limited has measured DC currents in transformer neutrals in the New Zealand electrical network at multiple South Island locations. Near‐continuous archived DC current data exist since 2001, starting with 12 different substations and expanding from 2009 to include 17 substations. From 2001 to 2015 up to 58 individual transformers were simultaneously monitored. Primarily, the measurements were intended to monitor the impact of the high‐voltage DC system linking the North and South Islands when it is operating in “Earth return” mode. However, after correcting for Earth return operation, as described here, the New Zealand measurements provide an unusually long and spatially detailed set of geomagnetically induced current (GIC) measurements. We examine the peak GIC magnitudes observed from these observations during two large geomagnetic storms on 6 November 2001 and 2 October 2013. Currents of ~30–50 A are observed, depending on the measurement location. There are large spatial variations in the GIC observations over comparatively small distances, which likely depend upon network layout and ground conductivity. We then go on to examine the GIC in transformers throughout the South Island during more than 151 h of geomagnetic storm conditions. We compare the GIC to the various magnitude and rate of change components of the magnetic field. Our results show that there is a strong correlation between the magnitude of the GIC and the rate of change of the horizontal magnetic field ( H′ ). This correlation is particularly clear for transformers that show large GIC current during magnetic storms.

Large arrays of tidal turbines are critical to realise the potential of tidal current power. This study is a systematic exploration of large tidal array optimisation in channels with numerically modelled array layouts in 2-D. Crucially, flow along channels is driven by head loss leading to significantly more realistic results than previous models which assume constant velocity. The 2-D adaptive mesh approach bridges the gap between large- and small-scale array models. Hundreds of layouts and turbine tunings have been simulated using LES of turbulent flow in tidally reversing currents to explore channel-scale optimisation and tuning of large arrays. Simulations show that total power capture increases as rows are added to the array although there are diminishing returns on additional turbines. Each turbine in 1 (7), optimally blocked row in a small channel captures 2.5× (0.5×) the power of an isolated turbine. There is an optimum blockage for maximum power per turbine which decreases linearly from 1.0 as the number of rows increases. As array size increases individual turbine wakes become less important than stepped head loss across each row. Free-stream velocity reduces linearly with total power capture, with the gradient increasing with channel size.

Much of the global tidal current energy resource lies in the accelerated flows along narrow channels. These channels have the potential to produce 10–1000s of MW of electricity. However, realizing 100MW of a channel׳s potential is much more complex than installing 100 1-MW turbines because large scale power extraction reduces tidal currents throughout the channel, changing the resource. This synthesis and review gives an overview of the issues and compromises in designing the layout of the large tidal turbine arrays required to realize this potential. The paper focuses on macro- and micro-design of arrays. Macro-design relates to the total number of turbines and their gross arrangement into rows, while micro-design adjusts the relative positions of the turbines within a grid and the spacing between rows. Interdependent macro-design compromises balance the total number of turbines, array power output, the power output of each turbine, the loads the turbines experience, turbine construction costs, maintaining navigability along the channel and any environmental impacts due to flow reduction. A strong emphasis is placed on providing physical insights about how “channel-scale dynamics” and the “duct-effect” impact on the compromises in array design. This work is relevant to the design of any “large” array which blocks more than 2–5% of a channel׳s cross-section, be it 2 turbines in a small channel or 100 turbines in a large channel.