Data for "Increasing Frequency and Persistence of the Summertime Greenland High Regime Not Captured by a Seasonal Prediction Model Very Large Ensemble"

These data support the work presented in Lee & Polvani (2025) "Increasing Frequency and Persistence of the Summertime Greenland High Regime Not Captured by a Seasonal Prediction Model Very Large Ensemble", published in Geophysical Research Letters: https://doi.org/10.1029/2025GL119421This work is based on the year-round weather regime classification introduced by Lee et al. (2023) and Lee & Messori (2024), with slight modifications as outlined in Lee & Polvani (2025).era5_1981_2024_year_round_north_america_weather_regimes_running_mean.ncThe ERA5 daily weather regime classification for 1981–2024 and associated parameters. The weather regime index (WRI) is also supplied here following Lee & Messori (2024). Dimensions:latitude = 61longitude = 151time = 16060eof_number = 12regime = 5trend_params = 2doy = 365Variables:latitude(latitude)longitude(longitude)time(time) (units = "days since 1970-01-01")z_anoms_detrend(time, latitude, longitude); units = "gpm" ; long_name = "Detrended Z500 anomalies 1981-2024 ERA5 climate"z_norm(time, latitude, longitude) ; units = "std" ; long_name = "Z500 anomalies normalized by seasonal cycle of area-averaged standard deviation"z_climo(doy, latitude, longitude) ; units = "gpm" ; long_name = "Z500 climatology 1981-2024 ERA5"trend(doy, trend_params) ; units = "m per day since 1981-01-01" ; long_name = "60d smoothed linear trend slope & intercept,1981-2024"norm_factor(doy) ; long_name = "Variance normalization factor"pc(time, eof_number) ; units = "raw pc" ; long_name = "pc time series"eof(eof_number, latitude, longitude) ; long_name = "raw EOFs"eigs(eof_number) ; long_name = "eigenvalues"var_frac(eof_number) ; long_name = "variance fraction"regime(time) ; long_name = "Daily regime attribution"WRI(regime, time) ; units = "std" ; long_name = "Weather regime index"WRI_mean(regime) ; long_name = "WRI mean"WRI_std(regime) ; long_name = "WRI standard deviation"cluster_mean_norm(regime, latitude, longitude) ; units = "std" ; long_name = "Cluster-mean normalized Z500 anomalies"double cluster_mean(regime, latitude, longitude) ; units = "gpm" ; long_name = "Cluster-mean Z500 anomalies"double cluster_centroids(regime, eof_number) ; long_name = "Cluster centroids in PC space"double dist_to_regime_centroid(time) ; long_name = "Distance to centroid of assigned regime"seas5_may_1981_2024_jja_year_round_north_america_weather_regimes_monte_carlo_10000_member.ncThe weather regime classification for 1 June–31 August using a Monte Carlo/random sampling of ECMWF's SEAS5 hindcasts & forecasts from 1981–2024. Variables are defined similarly to the ERA5 dataset above.Dimensions:lead_time = 92year = 44ens_member = 10000cluster = 5pc = 12trend_parameter = 2Variables:year(year)lead_time(lead_time) (units = "days since 1 May")ens_member(ens_member) regime(ens_member, year, lead_time) # Daily regime attributionwri(ens_member, year, lead_time, cluster) # Weather regime indexpcs(ens_member, year, lead_time, pc) # PCs obtained from projecting onto ERA5 EOFsnorm_factor(ens_member, lead_time)  # variance normalisation factortrend(ens_member, lead_time, trend_parameter) # Z500 trend in each ensemble member (slope + intercept)member_id(ens_member, year) # this corresponds to the true ensemble member selected in each random sampling