Over 20 years of egg sampling data were used to reconstruct the geographical and phenological patterns of walleye pollock (Theragra chalcogramma) spawning aggregations in the Gulf of Alaska (GOA). The analyzed time series (1972, 1978–1979, 1981–2000) included a documented event of climate change (i.e., 1988–1989) and the rise and fall of the GOA pollock population abundance and harvests. We compared results from two generalized additive model (GAM) formulations: one assuming no change of egg distribution and phenology over the examined time series (stationary) and the other admitting such changes (nonstationary) across an epoch determined from the data. Results from both model formulations corroborate the existence of a high egg concentration in Shelikof Strait, historically the primary spawning area of pollock in the GOA. However, model results also highlight the presence of other secondary, and possibly transitory, centers of egg distribution at various locations along the shelf and slope regions of the GOA. In addition, results from the nonstationary (and statistically superior) formulation indicate that the abundance of the non-Shelikof aggregations has increased over time, along with a tendency for earlier occurrence and displacement toward shallower areas of the high egg density regions.