Evaluating the Potential for Tidal Phase Diversity to Produce Smoother Power Profiles

Although tidal energy conversion technologies are not yet commercially available or cost-competitive with other renewable energy technologies like wind turbines and solar panels, tides are a highly predictable resource. Tidal energy’s predictability indicates that the resource could introduce less volatility into balancing the electric grid when compared to other renewables, a fundamentally desirable attribute for the electric system. More specifically, tidal energy resources are unique in that they have the potential to produce relatively smoother power profiles over time through aggregation. In order to generate smooth power profiles from tidal resources, sufficient complexity within the timing of tides is necessary within electrical proximity. This study evaluates the concept of aggregating diverse tides for the purpose of reducing periods of no and low energy production and creating smoother power profiles in regions around Alaska and Washington by calculating cross-correlations of tidal current velocity time series. Ultimately, study results show limited potential to exploit the resources for this purpose and describe the institutional mechanisms necessary to realize the benefits in practice.


Introduction
Unlike other renewable energy sources like wind and solar, tides are highly predictable without specialized forecasting. As a result of this predictability, tidal energy conversion technologies may become more competitive as a renewable resource that introduces less uncertainty-and thus, less volatility-into balancing the electric system. Although predictable, tides are still variable due to the constantly changing tidal energy resource. Unlike conventional power plants with controlled fuels, power plants that are fueled by fluctuating renewable resources are termed variable energy resources (VER). The variable nature of wind, Sun, subsurface heat, and river flows are fundamentally different Earth system processes, resulting in volatile power profiles from renewable energy power plants. While VERs represent clean, non-emitting sources of power generation, their natural variation in production due to changing weather conditions and daily cycles requires generating reserves over many time intervals and volumes. In order to generate a smooth power profile from tidal resources, sufficient complexity within the timing of tides is necessary; the resource must be diverse in proximity. If these characteristics are realized, it might be possible to demonstrate that tidal phase diversity can reduce periods of no or low production and generate smoother power profiles; ultimately reducing costs to balance the electric grid compared to other resources. 2 of 17 Several international studies including ones around the United Kingdom, Ireland, and northwestern Europe have evaluated the potential of tidal phase diversity to provide either baseload power (i.e., a flat power profile) or firm power (i.e., continually available power during a designated time of commitment) through aggregation [1][2][3]. Clarke et al. [1] show that a small amount of baseload power can be achieved by aggregating complementary sites around the United Kingdom and that limiting turbine capacity improves capacity factors and project economics. Notably, there is a relationship between the distances required to achieve sufficient tidal phasing for aggregation and the ability to achieve certain grid benefits within reasonable costs. In comparison, Neill, Hashemi, and Lewis [3] and Giorgi and Ringwood [2], respectively, show that there is minimal phase diversity in high tidal stream regions around northwestern Europe indicating the possibility of a limited amount of firm power and that a minimal amount of baseload power can be achieved from diverse tides around Ireland. No such study has been conducted for the United States, but the strength of tidal resources around the country has been assessed.
The Electrical Power Research Institute provided a first evaluation of the most energetic tidal stream sites in the United States. Sites in Alaska [4], California [5], Massachusetts [6], Maine [7], and Washington [5,7] were identified as having potential for energy extraction. Haas et al. [8] provided the first synoptic view of the tidal stream by mapping the resource in all U.S. coastal regions (except Hawaii and territories in the Pacific Ocean). The study is divided into four different regions in the U.S. Pacific Coast: Bristol Bay, Alaska; Cook Inlet and Kodiak, Alaska; Inside Passage, Alaska; and Salish Sea, Washington.
This study surveys areas known to have promise for tidal phase diversity by assessing time series of tidal velocities of representative locations in coastal areas of Alaska and Washington. This paper builds upon the work completed in [9], which analyzed the M 2 tidal constituent for the purpose of providing smoother power output through aggregation around the United States. As the most energetic tidal constituent, the M 2 represented the dominant trend in the tidal velocity patterns over time, offering a meaningful first-order analysis of the phenomenon. The work that follows introduces the full suite of tidal constituents from modeled data and verifies that modeled data to help identify locations where tidal resource aggregation (i.e., the sum of tidal power at a cluster of sites) can benefit the grid by minimizing the duration of time during which there is no power availability (slack tides) and producing smoother power profiles in aggregate.
This study will not analyze the amount of baseload or firm power that can be achieved or optimize deployment or centralized electricity dispatch decisions as that will require assumptions on the efficiency of technologies and intensity of energy extraction. Furthermore, achieving perfect baseload power is not the only way that tidal phasing can benefit the grid. Reducing the variability or fluctuation in power over time is valuable. Thus, through this inquiry, tidal phase diversity as a feasible and meaningful concept is evaluated.
To that end, Section 2 provides background information on the value of tidal energy in the United States' grid and the natural conditions necessary to realize diverse tides within proximity. This is followed by a description of the employed mathematical methods and the data sources used in the analysis in Section 3. Results from the analysis are presented in Section 4, followed by an exploration of the practical implications for achieving smoother power profiles in practice in Section 5 and concluding remarks in Section 6.

The Value of Tidal Energy in the United States' Grid
Energy generation sources across the United States are transitioning from conventional energy sources that emit greenhouse gasses to cleaner, renewable energy sources. By the end of 2018, 29 states and Washington D.C. had adopted renewable portfolio standards that dictate the percentage of electricity that must be delivered from specified renewable sources by a given date [10]. These renewable portfolio standards widely vary in their goals, and many have recently updated their policies to increase their renewable supply targets from original percentages [10]. Solar and wind energy, two VERs, are dominating deployment trends, accounting for at least 95% of domestic renewable energy capacity additions [11] and 84% of global renewable energy capacity additions in 2018 [12].
As the amount of VERs increases, balancing the electric grid to ensure supply meets demand becomes increasingly complex and costly. This is particularly true for small or remote grids that are more susceptible to disruptions from variations in production and consumption. Even in a system without a significant VER penetration, operators have some foresight but limited control of electric loads. To compensate for this lack of control, traditional generating resources are designed to be predictable and controllable [13].
Notably, tidal current devices are not presently cost-effective nor have they reached a commercial state; however, given the range of grid conditions across the country and the forward-looking energy policies striving for high penetrations of renewables, it is critical to show that tidal energy resources can provide favorable characteristics from a grid perspective. Enabling decisionmakers to include and embrace these technologies as solutions and contributors in the future is necessary to reach policy goals. Identifying tidal energy's grid-friendly traits is one step towards rethinking our conventional paradigms to deploying and valuing energy sources.
The mechanism and the rate at which renewable energy resources are adopted will differ across the country due to resource availability, favorable policy, but also the characteristics of energy capture technologies and the effects of integrating that renewable resource into the grid. As the electric grid and the generation resources feeding into it are evolving, decisionmakers have begun to modernize the way in which generation sources are incentivized and compensated, and prioritize locations where development is encouraged [14]. Innovative thinking around valuation permits grid planners, developers, and policymakers to better identify the future benefits that emerging technology could provide to a changing grid.

Tidal Conditions for Phase Diversity
Tides are caused by the gravitational interactions between the Earth, Moon, and the Sun. Because of the periodicity of the orbital motions, tides are reliable and predictable with a high level of accuracy. This predictability in the resource translates to predictability in energy production. In many parts of the world, the principal lunar semi-diurnal constituent (M 2 ) is the dominant tidal constituent. This tidal generating force occurs because of the gravitational interaction between the Earth and the Moon with a period of 12 h and 25 min. This period is a result of the 24-h rotation of the Earth and the rotation of the moon around the Earth.
Tides are long waves due to their long periods with respect to the depth of the ocean. Their behavior can be described by the shallow water wave equations. Long waves travel at a speed controlled by the water depth: where c is the phase velocity, L is the wave length, T is the wave period, σ (=2π/T) is the radian frequency, k (=2π/L) is the wave number, g is the gravitational acceleration, and h is the water depth. For instance, for an M 2 tide propagating along a 200 m deep shelf the tidal phase velocity equal to gh = 44.3 m/s, with an associated wave length of T gh = 1980 km. In an inviscid environment and in the absence of interactions (e.g., friction, obstructions), tides propagate as progressive waves such that the water surface elevation (η) and the wave-induced horizontal velocity (u) are in phase with the following relationship: where H is the tidal range. The tidal currents are maximum under the crest and trough of the long waves. Continuing with the shelf example, regions of absolute maximum and zero horizontal velocity at the same time are 495 km apart at L/4. For that idealized tide with only an M 2 component, if two extraction sites are located at 495 km and the device can harvest in both ebb and flood tide, the aggregate theoretical power of the resources at those sites will continually be nonzero. As tides enter coastal environments, they become increasingly affected by bottom friction and by interactions with coastal features such as inlets, fjords, estuaries, etc. These interactions alter the phase relationship between η and u. This might lead to sites in proximity closer than L/4 to experience complementary tidal velocities. To account for this effect, this study considers data from multiple locations along the coasts of interest, which will be further described in Section 3.
In this study sites with high tidal current velocities are identified and the phase shifts of the tides to identify sites where large phase differences can be achieved at spatial scales shorter than L/4 are evaluated. For a topographic example of clustered but variable phasing, consider the extreme case where the ocean is connected to a large inland sea through a narrow inlet. Maximum currents will occur at the mean water level in the coastal sea. This results in a 90 • phase shift between the water surface elevation and the currents and in turn a correlation of 0 between the signals (this is further discussed in Section 3.2). Additionally, tidal wave reflection and advection effects can also be present in a given system. The natural occurrence of these physical processes leads to many solutions to the wave equations where the phase of the tides is non-uniform. Due to the complex tidal circulation in natural estuaries, inlets, bays, etc. this problem cannot be fully analyzed theoretically; instead, numerical models are used in these scenarios.
Tidal currents are the quantity of interest because they are directly proportional to the instantaneous in-stream power per unit area: where ρ is the water density. Realistically, the amount of extracted energy will be less than the theoretical resource and will depend upon energy capture efficiency and intensity of extraction. These can be included in the model by multiplying Equation (3) by the area of extraction (i.e., turbine area and number of turbines) and an efficiency factor.

Data Sources
Model-derived tidal currents are the primary source of data because they provide coverage at a high spatial resolution with extensive spatial coverage. Observations are used in the Salish Sea to verify the results obtained using the model-derived database; this is discussed in Appendix A.
In order to evaluate the possible complementary sites that would reduce periods of no or lower production, a consistent dataset that covers the entire area of interest is necessary. The main dataset used in this study comes from Haas, Fritz, French, Smith, and Neary [8] (H11 hereafter), which is a model-based study of tidal streams for the mainland United States, the Pacific coast of Alaska, Puerto Rico, and the U.S. Virgin Islands. Multiple overlapping models based on the Regional Ocean Modeling System (ROMS) were used to hindcast the tidal streams for a 32-day period, of which the first two days were discarded to account for spin-up effects. The models had an average spatial resolution of 350 m, for a total of 3,247,981 grid points. The models were forced using the Q 1 , O 1 , K 1 , S 2 , M 2 , N 2 , and M 4 constituents. The harmonic constituents and mean currents at each output point were computed using harmonic analysis on the 30-day time series. Tidal ellipses at each grid point for the main tidal constituents were available to the authors. The time series were recovered from the ellipses for each tidal constituent using the tidal ellipse toolbox for MATLAB. The final velocity time series were obtained by a adding time series from the different constituents.
Results from two refined models in the Salish Sea [15,16] (YW14 hereafter) and Cook Inlet, AK [17] (WY19 hereafter), are used to complement the results from H11. These are numerical models based on the Finite-Volume, Community Ocean Model (FVCOM) [18,19]. Similar to ROMS, FVCOM solves the Reynolds-averaged Navier-Stokes equations. However, this model implements a finite-volume approach making it suitable for coastal applications with complex shoreline geometries and topography. YW14 has a spatial resolution of~200 m in the Puget Sound and is relaxed to 2000 m at the open boundaries outside of the Puget Sound. The model was forced with observed tides at two tidal gauges located near Neah Bay, Washington and Campbell River, British Columbia, respectively. Similarly, WY19 implements a variable resolution from 200 to 300 m inside the inlet to~1000 m at the open boundary. This model was forced with Oregon State University's TPXO8-atlas tidal database [20], which consists of 13 major harmonic constituents.

Methodology
With the goal of identifying sites with diverse tides that result in complementary resources, the tidal velocity time series were initially filtered for feasibility. Following similar studies estimating first-generation tidal technologies like [3], locations with modeled data with mean tidal velocities under 50 cm/s or a minimum water depth under 10 m were discarded.
Once the modeled data points below the threshold were discarded, representative locations were selected within each of the regions, one per feature (i.e., inlet, strait, headland). A cross-correlation at lag zero was calculated for the time series between each pair of locations within a region, such that where r xy is the correlation coefficient, x(t) and y(t) are the tidal velocity time series, and n x and n y are the mean values of their respective time series. The cross-correlation indicates a measure of the time series' movement with respect to one another. Positive correlation coefficients indicate temporal similarity between the time series patterns while negative correlation coefficients indicate dissimilarity. For the purpose of identifying complementary sites, negative correlations are ideal. The correlation coefficient does not reflect the magnitude of the sources creating a perfect offset of one another but represents the pattern of the fluctuations between the resources offsetting each other over time. Thus, a strong anticorrelation indicates decreased times of no or lower power availability when the resources are aggregated and the potential to achieve a smoother power profile. A matrix of correlation coefficients is presented for the locations within each region.
One alternate way to analyze sites for complementary phases is to compare the phases of the tidal components, e.g., [3]. This method can be implemented if looking at one tidal constituent only. Since the present study considers all the tidal constituents, this method cannot be implemented and thus the correlation method is chosen.

Salish Sea
The Salish Sea is a large inland sea located in the Pacific Northwest Coast, bounded by the coasts of the State of Washington, USA and British Columbia, Canada. The Salish Sea is identified as a region with high potential for in-stream tidal energy extraction [21]. Previous studies have identified regions for potential tidal energy harvesting: Admiralty Inlet [16,21], Tacoma Narrows [16], Agate Pass, and Rich Passage [15]. In addition, sites in the Puget Sound with mean currents above 50 cm/s include Deception Pass, Dana Passage, Pickering Passage, and Port Washington Narrows. The flow around the San Juan Islands is also known to be energetic. To sample all channels with mean flows exceeding the 50 cm/s threshold, 20 representative points are evaluated as shown in Figure 1. In Figure 1b, cross-correlations are calculated and displayed with respect to the anchor point (0) to broadly depict diversity in the region. The cross-correlations between all sites are shown in Table 1.  Figure A1.
Most of the points show correlations above 0, and only Deception Pass (Points 10 and 11) shows negative correlations with multiple points around the region (Table 1). As indicated by the positive correlations, there is no ideal offset generated by the tides in this region apart from combining some points with Deception Pass.   Figure A1.
Most of the points show correlations above 0, and only Deception Pass (Points 10 and 11) shows negative correlations with multiple points around the region (Table 1). As indicated by the positive correlations, there is no ideal offset generated by the tides in this region apart from combining some points with Deception Pass.
Two distinct regions are identified: San Juan Islands to the north and the Puget Sound to the south of Admiralty Inlet. In the San Juan Islands, multiple channels allow tidal flows to propagate between the Strait of Juan de Fuca and the Strait of Georgia. In the Puget Sound, most of the flow enters through Admiralty Inlet and the rest through Deception Pass [21].
In the San Juan Islands (Points 1-9), high correlations in the tides suggest minimal obstruction for propagation between the Strait of Juan de Fuca and the Strait of Georgia. Correlations are generally higher with one other and with Admiralty Inlet (Point 0), showing unfavorable conditions for aggregation to benefit the grid.
In Puget Sound (Points 0 and [12][13][14][15][16][17][18][19], the correlations decrease with distance from Admiralty Inlet (Point 0). This is a result of the tide progressing through the sound. Correlations are above 0.60 for the main channels indicating the tides flow through the main channels quickly. Maximum correlation between Admiralty Inlet and the South Sound is achieved by shifting one of the time series by one hour, indicating the time of travel of the tides between these points. This is consistent with Yang, Wang, Copping, and Geerlofs [16], who showed that the phase difference between Pickering Passage (and the southern end of the Puget Sound) is about 40 min behind Tacoma Narrows for the M 2 tide. The lowest 7 of 17 correlation occurs with Point 14 at Port Washington Narrows. Tides entering the Puget Sound cross Point 0, then 12 and 13, before reaching Dyes Inlet through the Port Washington Narrows (Point 14). The high correlation between 0 and 12 and 13 suggests the tidal waves are slowed down near Dyes Inlet, thus increasing the lag between these points. Deception Pass (Points 10 and 11) shows negative correlations with many of the points in the Salish Sea, including Points 3, 6-7, 12-17, and 19. These results indicate Deception Pass plays a complementary role to many other channels in the region but to a limited degree. It can become an important location for energy harvesting due to this characteristic. Although to varying degrees, having the potential to complement several locations from a timing perspective offers flexibility in planning when trying to achieve smoother power profiles through aggregation.

Alaska
Alaska contains many sites with energetic tidal flows. Due to its size, three subregions are independently analyzed: Bristol Bay, Kodiak and Cook Inlet, and Inside Passage.

Bristol Bay
Bristol Bay is the easternmost reach of the Bering Sea, it is known for having large tides in the upper reaches. The M 2 constituent is the most energetic one in this region and contains an amphidromic point near the north of the Bay [22]. Thus, tides follow a rotational motion around the bay, enhancing the phase diversity relative to a bay where tides propagate perpendicular to the uppermost reaches. In addition, there are several sites that exceed the velocity and depth thresholds set in the methodology. Ten representative sites were selected in this area as shown in Figure 2. uppermost reaches. In addition, there are several sites that exceed the velocity and depth thresholds set in the methodology. Ten representative sites were selected in this area as shown in Figure 2. Starting from the anchor point the correlations decrease clockwise up to Point 9 at the western shore of Unimak Island, off the southwestern end of the Alaska Peninsula (Table 2). This trend indicates that the clockwise rotation identified for the M2 constituent for water surface elevation is dominant when currents from the full array of constituents are considered. Many points show anticorrelations that indicate the possibility of providing smoother power profiles and reducing times of no and low production. Points 7, 8, and 9 are anticorrelated with the highest number of other sites, and those three points largely contain the strongest anticorrelations as well. Thus, not only is there flexibility in the resources that they might aggregate with for energy production, they show the highest probability of achieving smoother power profiles in the region.

Kodiak and Cook Inlet
Tides in the Gulf of Alaska travel from east to west [23,24] around the amphidromic point that lies between Hawaii and California. The open ocean phase difference of the M2 tide is in the order of 20° (or 40 min) between the edge of the Kenai peninsula to the southwestern end of Kodiak Island [25]. This small phase difference is enhanced when the tide flows into Cook Inlet and the channels Starting from the anchor point the correlations decrease clockwise up to Point 9 at the western shore of Unimak Island, off the southwestern end of the Alaska Peninsula (Table 2). This trend indicates that the clockwise rotation identified for the M 2 constituent for water surface elevation is dominant when currents from the full array of constituents are considered. Many points show anticorrelations that indicate the possibility of providing smoother power profiles and reducing times of no and low production. Points 7, 8, and 9 are anticorrelated with the highest number of other sites, and those three points largely contain the strongest anticorrelations as well. Thus, not only is there flexibility in the resources that they might aggregate with for energy production, they show the highest probability of achieving smoother power profiles in the region.

Kodiak and Cook Inlet
Tides in the Gulf of Alaska travel from east to west [23,24] around the amphidromic point that lies between Hawaii and California. The open ocean phase difference of the M 2 tide is in the order of 20 • (or 40 min) between the edge of the Kenai peninsula to the southwestern end of Kodiak Island [25]. This small phase difference is enhanced when the tide flows into Cook Inlet and the channels and banks around the islands. Twelve representative points were selected in this area as shown in Figure 3. and banks around the islands. Twelve representative points were selected in this area as shown in Figure 3. In the channels southwest of Kodiak Islands, Points 0, 2, and 3 are mostly in phase (Table 3), indicating minor tidal rotation around Sitkinak Island. Point 1 located at the Olga Narrows, a narrow strait connecting Olga Bay with the open ocean, shows anticorrelation with the neighboring points. The cross-correlations between the velocity and water surface elevation time series for Points 0 through 3 are 0.17, 0.94, −0.29, and −0.19. This shows that the flow and water surface elevation along the Olga Narrows are in phase; that is not the case for other locations around the islands. Tidal flows around Afognak Island (Points 4, 5, and 6) are also in close phase to each other and to Points 0, 2, and 3. They are in a very close phase with the entrance of Cook Inlet at Point 7. Achieving smoother power profiles within the Kodiak and Afognak system does not appear feasible based on this analysis.
Cook Inlet is known for its large tidal range and fast currents. The M2 tide is again the most energetic component. The high tide associated with this constituent takes more than 4 h to propagate from Point 7 to Point 11 [26]. The points selected to investigate the tidal currents sample locations along the inlet. Tidal currents show a similar picture as shown by the gradients in correlations along the inlet (Figure 3b). The correlation drops monotonically from Point 8 to 11 (Table 3). In the channels southwest of Kodiak Islands, Points 0, 2, and 3 are mostly in phase (Table 3), indicating minor tidal rotation around Sitkinak Island. Point 1 located at the Olga Narrows, a narrow strait connecting Olga Bay with the open ocean, shows anticorrelation with the neighboring points. The cross-correlations between the velocity and water surface elevation time series for Points 0 through 3 are 0.17, 0.94, −0.29, and −0.19. This shows that the flow and water surface elevation along the Olga Narrows are in phase; that is not the case for other locations around the islands. Tidal flows around Afognak Island (Points 4, 5, and 6) are also in close phase to each other and to Points 0, 2, and 3. They are in a very close phase with the entrance of Cook Inlet at Point 7. Achieving smoother power profiles within the Kodiak and Afognak system does not appear feasible based on this analysis. Cook Inlet is known for its large tidal range and fast currents. The M 2 tide is again the most energetic component. The high tide associated with this constituent takes more than 4 h to propagate from Point 7 to Point 11 [26]. The points selected to investigate the tidal currents sample locations along the inlet. Tidal currents show a similar picture as shown by the gradients in correlations along the inlet (Figure 3b). The correlation drops monotonically from Point 8 to 11 (Table 3).
At the upstream end of Cook Inlet, Knik and Turnagain arms might extend the phase diversity of the system but are not included in this analysis because this area is not permanently submerged (Figure 4b)

Inside Passage
There are several locations in the Inside Passage of Southeast Alaska where the mean tidal velocity exceeds 50 cm/s. Some of these sites are in very shallow waters or narrow channels that, due to the model resolution, cannot be accurately modeled with a 350 m resolution. Thus, only seven points are selected for analyzing this region as shown in Figure 5.  As confirmation of the correlations, data from Wang and Yang [17] (WY19) are used. The model extent and correlations are shown in Figure 4. Table 4 shows the correlations of the points inside the model domain. Overall the results are consistent between both models providing confidence in the results. Contrary to the analysis with the H11 dataset, these correlations were derived directly from the hourly model output. This indicates that non-tidal residuals are present, but results show the same behavior (Table 4). The large time lags in M 2 water level propagation [17,26] can be seen from the cross-correlation analysis of the complete time series as well (Figure 4a). Correlations are taken from a point in the deepest part of the offshore boundary condition. The cross-correlation decreases monotonically upstream of the inlet as the tide propagates inward. This, however, is not reflected in the tidal currents as they fall out of phase from the water surface elevation (Figure 4b). The currents decrease monotonically upstream of Point 9. There is a small level of diversity in the upper reaches of Cook Inlet, with differences in correlations of 0.5. To enhance this diversity, points outside of the inlet must be considered.

Inside Passage
There are several locations in the Inside Passage of Southeast Alaska where the mean tidal velocity exceeds 50 cm/s. Some of these sites are in very shallow waters or narrow channels that, due to the model resolution, cannot be accurately modeled with a 350 m resolution. Thus, only seven points are selected for analyzing this region as shown in Figure 5.

Inside Passage
There are several locations in the Inside Passage of Southeast Alaska where the mean tidal velocity exceeds 50 cm/s. Some of these sites are in very shallow waters or narrow channels that, due to the model resolution, cannot be accurately modeled with a 350 m resolution. Thus, only seven points are selected for analyzing this region as shown in Figure 5. In this region, the tide propagates northward and the phase lag for the M2 tidal constituent, the most energetic, is around 25° [22]. However, there is a significant phase difference when the tidal currents are evaluated. Multiple sites show negative correlations as shown in Table 5. Tidal currents north and south of Chichagof Island are out of phase where Point 3 shows correlations between −0.61 and −0.4 with points in the north. Overall, the potential to aggregate resources to reduce periods of In this region, the tide propagates northward and the phase lag for the M 2 tidal constituent, the most energetic, is around 25 • [22]. However, there is a significant phase difference when the tidal currents are evaluated. Multiple sites show negative correlations as shown in Table 5. Tidal currents north and south of Chichagof Island are out of phase where Point 3 shows correlations between −0.61 and −0.4 with points in the north. Overall, the potential to aggregate resources to reduce periods of no or low power production and achieve smoother power profiles in the region exists at Points 3, 5, and 7. Point 3 is also anticorrelated with Points 0 and 1, which are all around Chichagof Island. Local effects play a dominant role since the water surface elevation between Points 1 and 3 are highly correlated at 0.95 (not shown). Further, cross-correlation analysis of water surface elevation is at least 0.89 for all points in the region. Similarly, Point 7 is not correlated with the northern tidal flows. This case is opposite to Cook Inlet where there is more phase diversity in the water surface elevation than in the tidal currents.

Discussion
Of the evaluated regions, pairs of points in Cook Inlet and Bristol Bay indicate the highest probability of reducing times of no power availability and smoother power output through aggregation. The extent to which the resource can be exploited to realize these benefits depends on development decisions, device controls, project economics, infrastructure, and institutional limitations. It is important to note that the analysis considers the correlation of the velocity time series and indicates patterns between the resources, not actual power produced from the specific extraction of the resources. Therefore, even if a strong anticorrelation exists between the resources at two locations, the magnitude of those resources is not explicitly considered. That is to say that the general pattern and fluctuation of the resources complement one another at the evaluated time scale, but one resource could have a much higher average velocity at one location compared to the other and would necessitate development decisions that capitalize on the behaviors that enable smoother power profiles in aggregate.
Results with strong anticorrelations indicate a decrease in times of no or low power availability and thus, the possibility of producing smoother power profiles in aggregate. effect of this asymmetry can be large and must be analyzed in detail at the time of selecting the technology and intensity of extraction. Siting tidal energy converters in ebb and flood dominated parts of a strait has been shown to reduce the variability in peak power production [27,28]. Moreover, the rectilinear misalignments between the dominant flood and ebb directions will have additional impacts on actual power production. Thus, all these factors should be carefully considered in future, finer-scale case studies. There are a number of ways in which smoother power output can be achieved and optimized in practice. They are dependent upon site-specific needs, including those dictated by the grid, and those needs vary by location. Aggregating diverse tides generates an optimization question of how and the extent to which the resources should be harnessed. Part of this optimization is not simply evaluating the hour by hour comparison between the resources but also analyzing the aggregate power over the spring-neap cycle. Weaker resources at one location might need to be extracted at a higher intensity than complementary resources at a location with stronger current velocities to reduce periods of no and low production. To produce smoother power from complementary resources, limiting turbine capacity is an option to cap the maximum production from machines and decrease the range in power output, as also noted in [1]. There are economic implications of limiting power for the purpose of The intensity at which the resources are extracted is also important, as that contributes to the potential for smoother power output. By scaling the theoretical power time series to have the same maximum power value (i.e., mimicking disproportionate development to generate similar amounts of power by increasing the capture area), Figure 6b is generated. In comparison to Figure 6a, the aggregate power does not drop below 200, the maximum aggregate power only increases trivially, and the continual steepness of the ramping has decreased. With fewer dips in available power, there is the opportunity to generate smoother power by limiting the upper capacity of production. It is important to note that these estimations do not consider the efficiency of the devices cut-in or cut-out speeds associated with turbines, but rather focus on the theoretical power. Further, the effects of tidal asymmetry due to the superposition of various tidal constituents can be seen in the example aggregated power (Figure 6) where consecutive cycles can have variations of more than 100% in peak power. When multiple sites are considered in order to achieve continuous power, the aggregated effect of this asymmetry can be large and must be analyzed in detail at the time of selecting the technology and intensity of extraction. Siting tidal energy converters in ebb and flood dominated parts of a strait has been shown to reduce the variability in peak power production [27,28]. Moreover, the rectilinear misalignments between the dominant flood and ebb directions will have additional impacts on actual power production. Thus, all these factors should be carefully considered in future, finer-scale case studies.
There are a number of ways in which smoother power output can be achieved and optimized in practice. They are dependent upon site-specific needs, including those dictated by the grid, and those needs vary by location. Aggregating diverse tides generates an optimization question of how and the extent to which the resources should be harnessed. Part of this optimization is not simply evaluating the hour by hour comparison between the resources but also analyzing the aggregate power over the spring-neap cycle. Weaker resources at one location might need to be extracted at a higher intensity than complementary resources at a location with stronger current velocities to reduce periods of no and low production. To produce smoother power from complementary resources, limiting turbine capacity is an option to cap the maximum production from machines and decrease the range in power output, as also noted in [1]. There are economic implications of limiting power for the purpose of providing smoother power such that decreasing turbine capacity can improve a project's capacity factor and increase project economics [1], but the amount of energy produced over time will be reduced. However, even if the tidal resources are present, institutional mechanisms must be in place for the grid to realize the complementarity benefits and to produce favorable economic results.
In the analysis, reasonable spatial boundaries to "cluster" tidal profiles were applied. Conceptually, each individual tidal project or turbine would not be directly or operationally connected to one another in the water; rather, the smoother net profile or reduction in periods of no or low production is achieved by grid-side aggregation. Aggregator is a term commonly applied to a third party or activity that compiles distributed energy resources into one virtual power plant in order to have a single point of economic transaction. While aggregation is sometimes useful for small generators that are not large enough to enter into complex markets, for tidal resources, a greater value may be achieved through aggregating profiles into a predictable block of energy.
There are limits to aggregation, however. This is principally linked to whether the facilities are located in the same electricity management footprint, rather than whether they are physically nearby. In this electricity footprint, one operator balances electric load and generation in real-time. Electricity boundaries are not necessarily related to political or geographical boundaries, and in some cases, the footprint can be very large. Almost all of the state of California is managed by one operator, the California Independent System Operator (CAISO), which offers options for distributed energy resource aggregators [29].
Therefore, in order to capture the grid benefits of tidal phase diversity for a reduction in periods of no or low production or a smoother power profile, facilities do not need to be close physically, but they need to be close electrically. Once tidal phase diversity is identified, the next steps are to establish points of interconnection and evaluate the ability to aggregate and reach a single off-taker, whether a load-serving utility or a market. In environments with many small isolated electrical systems-e.g., as is common in Alaska-aggregation may not be feasible. Similarly, even perfectly phased tidal regimes, if spread over a very large geographic range, would likely be economically impractical if not under a single system operator. Understanding these requirements in locations where aggregating tidal resources may provide grid benefits is key to moving from theory to reality.

Conclusions
This paper shows there is limited potential in areas of Alaska and Washington for tidal phase diversity to reduce periods of no or low power availability and generate smoother power profiles through aggregation. Locations around Kodiak Island and Cook Inlet show the greatest potential out of the evaluated regions with each of the sampled locations there anticorrelated with at least one other location and several sites showed strong anticorrelations. Locations at Bristol Bay and Inside Passage also present the opportunity to exploit resources for this purpose, but to a more limited degree and over a larger geographic footprint. Finally, locations around Washington are significantly more restricted in their feasibility based on the available resources.
Data from [16] were verified with measurements and other modeled data in the Salish Sea [15,16] and also compared to modeled data from [17] in Alaska. Practical considerations for the phenomenon to be realized in practice were outlined, and it was established that while the concept of tidal phase diversity has limited potential to feasibly benefit the grid around Washington and Alaska, site-specific analyses would be required to determine the practicality of developing the resources in any identified locations for this purpose.
Current Profilers (ADCP). For consistency with the modeled data, the measured currents are depthaveraged before analyzing. Simultaneous measurements in the different study regions with at least 15 days of overlapping data are identified and downloaded. The stations used in the analysis are identified in Figure A1 and Table A1. The Salish Sea has many channels and islands giving it a complex flow pattern. Since time series with multiple constituents are used, it is desirable to cross-compare the results obtained with the H11 model to another model and measurements. The National Oceanic and Atmospheric Administration conducted measurement campaigns in the Puget Sound between 2015 and 2018. The Salish Sea was sampled in sub-regions, thus an analysis of the 20 stations cannot be replicated, but the correlations between groups can be verified. The model-derived data does not overlap temporally with the measurements; thus, direct validation of the models cannot be provided. Cross-correlations are computed between stations that were measured simultaneously as shown in Table A1. Locations for the verification stations are shown in Figure 1a. Results from the YW14 model are sampled at the same stations as the H11 model.  The Salish Sea has many channels and islands giving it a complex flow pattern. Since time series with multiple constituents are used, it is desirable to cross-compare the results obtained with the H11 model to another model and measurements. The National Oceanic and Atmospheric Administration conducted measurement campaigns in the Puget Sound between 2015 and 2018. The Salish Sea was sampled in sub-regions, thus an analysis of the 20 stations cannot be replicated, but the correlations between groups can be verified. The model-derived data does not overlap temporally with the measurements; thus, direct validation of the models cannot be provided. Cross-correlations are computed between stations that were measured simultaneously as shown in Table A1. Locations for the verification stations are shown in Figure 1a. Results from the YW14 model are sampled at the same stations as the H11 model. Table A1. Measurements for model verification. Location of these stations are shown in Figure 1.

ID
Model ID Name Deployment the complexities of the study site. It is worth noting that non-tidal currents are not removed from the measurements and thus the information considered is different. The anticorrelation of Deception Pass with the main channels is confirmed by the measurements (Table A2).