skip to main content


Title: Wind Stress, Wind Stress Curl, and Upwelling Velocities in the Northwest Atlantic (80-45W, 30-45N) during 1980-2019

This dataset contains three netcdf files that pertain to monthly, seasonal, and annual fields of surface wind stress, wind stress curl, and curl-derived upwelling velocities over the Northwest Atlantic (80-45W, 30-45N) covering a forty year period from 1980 to 2019. Six-hourly surface (10 m) wind speed components from the Japanese 55-year reanalysis (JRA-55; Kobayashi et al., 2015) were processed from 1980 to 2019 over a larger North Atlantic domain of 100W to 10E and 10N to 80N. Wind stress was computed using a modified step-wise formulation, originally based on (Gill, 1982) and a non-linear drag coefficient (Large and Pond, 1981), and later modified for low speeds (Trenberth et al., 1989). See Gifford (2023) for more details.   

After the six-hourly zonal and meridional wind stresses were calculated, the zonal change in meridional stress (curlx) and the negative meridional change in zonal stress (curly) were found using NumPy’s gradient function in Python (Harris et al., 2020) over the larger North Atlantic domain (100W-10E, 10-80N). The curl (curlx + curly) over the study domain (80-45W, 10-80N) is then extracted, which maintain a constant order of computational accuracy in the interior and along the boundaries for the smaller domain in a centered-difference gradient calculation. 

The monthly averages of the 6-hour daily stresses and curls were then computed using the command line suite climate data operators (CDO, Schulzweida, 2022) monmean function. The seasonal (3-month average) and annual averages (12-month average) were calculated in Python using the monthly fields with NumPy (NumPy, Harris et al., 2020). 

Corresponding upwelling velocities at different time-scales were obtained from the respective curl fields and zonal wind stress by using the Ekman pumping equation of the study by Risien and Chelton (2008; page 2393). Please see Gifford (2023) for more details.   

The files each contain nine variables that include longitude, latitude, time, zonal wind stress, meridional wind stress, zonal change in meridional wind stress (curlx), the negative meridional change in zonal wind stress (curly), total curl, and upwelling. Units of time begin in 1980 and are months, seasons (JFM etc.), and years to 2019. The longitude variable extends from 80W to 45W and latitude is 30N to 45N with uniform 1.25 degree resolution.  

Units of stress are in Pascals, units of curl are in Pascals per meter, and upwelling velocity is described by centimeters per day. The spatial grid is a 29 x 13 longitude x latitude array. 

Filenames: 

monthly_windstress_wsc_upwelling.nc: 480 time steps from 80W to 45W and 30N to 45N.

seasonal_windstress_wsc_upwelling.nc: 160 time steps from 80W to 45W and 30N to 45N.

annual_windstress_wsc_upwelling.nc: 40 time steps from 80W to 45W and 30N to 45N.

Please contact igifford@earth.miami.edu for any queries. {"references": ["Gifford, I.H., 2023. The Synchronicity of the Gulf Stream Free Jet and the Wind Induced Cyclonic Vorticity Pool. MS Thesis, University of Massachusetts Dartmouth. 75pp.", "Gill, A. E. (1982). Atmosphere-ocean dynamics (Vol. 30). Academic Press.", "Harris, C.R., Millman, K.J., van der Walt, S.J. et al. Array programming with NumPy. Nature 585, 357\u2013362 (2020). DOI: 10.1038/s41586-020-2649-2.", "Japan Meteorological Agency/Japan (2013), JRA-55: Japanese 55-year Reanalysis, Daily 3-Hourly and 6-Hourly Data, https://doi.org/10.5065/D6HH6H41, Research Data Archive at the National Center for Atmospheric Research, Computational and Information Systems Laboratory, Boulder, Colo. (Updated monthly.)", "Kobayashi, S., Ota, Y., Harada, Y., Ebita, A., Moriya, M., Onoda, H., Onogi, K., Kamahori, H., Kobayashi, C., Endo, H. and Miyaoka, K., 2015. The JRA-55 reanalysis: General specifications and basic characteristics.\u202fJournal of the Meteorological Society of Japan. Ser. II,\u202f93(1), pp.5-48.", "Large, W.G. and Pond, S., 1981. Open ocean momentum flux measurements in moderate to strong winds.\u202fJournal of physical oceanography,\u202f11(3), pp.324-336.", "Risien, C.M. and Chelton, D.B., 2008. A global climatology of surface wind and wind stress fields from eight years of QuikSCAT scatterometer data.\u202fJournal of Physical Oceanography,\u202f38(11), pp.2379-2413.", "Schulzweida, Uwe. (2022). CDO User Guide (2.1.0). Zenodo. https://doi.org/10.5281/zenodo.7112925.", "Trenberth, K.E., Large, W.G. and Olson, J.G., 1989. The effective drag coefficient for evaluating wind stress over the oceans.\u202fJournal of Climate,\u202f2(12), pp.1507-1516."]} 
more » « less
Award ID(s):
2123283 1851242 2122726
NSF-PAR ID:
10448886
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Zenodo
Date Published:
Edition / Version:
1.0.0
Subject(s) / Keyword(s):
wind stress wind stress curl upwelling JRA-55 Northwest Atlantic Gulf Stream North Atlantic Python Computation CDO Biology
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. This dataset includes multiple fields: (i) files for monthly and annual fields for the max curl line and the zero curl line at 0.1 degree longitudinal resolutions; (ii) files for monthly and annual GS path obtained from Altimetry and originally processed by Andres (2016) at 0.1 degree longitudinal resolution. The maximum curl line (MCL) and the zero curl line (ZCL) calculations are briefly described here and are based on the original wind data (at 1.25 x 1.25 degree) provided by the Japanese reanalysis (JRA-55; Kobayashi et al., 2015) and available at https://zenodo.org/record/8200832 (Gifford et al. 2023). For details see Gifford, 2023. 

    The wind stress curl (WSC) fields used for the MCL and ZCL calculations extend from 80W to 45W and 30N to 45N at the 1.25 by 1.25-degree resolution.  The MCL is defined as the maximum WSC values greater than zero within the domain per 1.25 degree longitude. As such, it is a function of longitude and is not a constant WSC value unlike the zero contour. High wind stress curl values that occurred near the coast were not included within this calculation. After MCL at the 1.25 resolution was obtained the line was smoothed with a gaussian smoothing and interpolated on to a 0.1 longitudinal resolution. The smoothed MCL lines at 0.1 degree resolution are provided in separate files for monthly and annual averages (2 files). Similarly, 2 other files (monthly and annual) are provided for the ZCL.    

    Like the MCL, the ZCL is a line derived from 1.25 degree longitude throughout the domain under the condition that it's the line of zero WSC. The ZCL is constant at 0 and does not vary spatially like the MCL. If there are more than one location of zero curl for a given longitude the first location south of the MCL is selected. Similar to the MCL, the ZCL was smoothed with a gaussian smoothing and interpolated on to a 0.1 longitudinal resolution.   

    The above files span the years from 1980 through 2019. So, the monthly files have 480 months starting January 1980, and the annual files have 40 years of data. The files are organized with each row being a new time step and each column being a different longitude. Therefore, the monthly MCL and ZCL files are each 480 x 351 for the 0.1 resolution data. Similarly, the annual files are 40 x 351 for the 0.1 degree resolution data.  

    Note that the monthly MCLs and ZCLs are obtained from the monthly wind-stress curl fields. The annual MCLs and ZCLs are obtained from the annual wind-stress curl fields.

    Since the monthly curl fields preserves more atmospheric mesoscales than the annual curl fields, the 12-month average of the monthly MCLs and ZCLs will not match with the annual MCLs and ZCLs derived from the annual curl field.  The annual MCLs and ZCLs provided here are obtained from the annual curl fields and representative metrics of the wind forcing on an annual time-scale. 

    Furthermore, the monthly Gulf Stream axis path (25 cm isoheight from Altimeter, reprocessed by Andres (2016) technique) from 1993 through 2019 have been made available here. A total of 324 monthly paths of the Gulf Stream are tabulated. In addition, the annual GS paths for these 27 years (1993-2019) of altimetry era have been put together for ease of use. The monthly Gulf Stream paths have been resampled and reprocessed for uniqueness at every 0.1 degree longitude from 75W to 50W and smoothed with a 100 km (10 point) running average via matlab. The uniqueness has been achieved by using Consolidator algorithm (D’Errico, 2023). 

    Each monthly or annual GS path has 251 points between 75W to 50W at 0.1 degree resolution.  

    Please contact igifford@earth.miami.edu for any queries. {"references": ["Andres, M., 2016. On the recent destabilization of the Gulf Stream path downstream of Cape Hatteras. Geophysical Research Letters, 43(18), 9836-9842.", "D'Errico, J., 2023. Consolidator (https://www.mathworks.com/matlabcentral/fileexchange/ 8354-consolidator), MATLAB Central File Exchange. Retrieved June 17, 2023.", "Gifford, Ian. H., 2023. The Synchronicity of the Gulf Stream Free Jet and the Wind Induced Cyclonic Vorticity Pool. MS Thesis, University of Massachusetts Dartmouth. 75pp.", "Gifford, Ian, H., Avijit Gangopadhyay, Magdalena Andres, Glen Gawarkiewicz, Hilde Oliver, Adrienne Silver, 2023. Wind Stress, Wind Stress Curl, and Upwelling Velocities in the Northwest Atlantic (80-45W, 30-45N) during 1980-2019, https://zenodo.org/record/8200832.", "Kobayashi, S., Ota, Y., Harada, Y., Ebita, A., Moriya, M., Onoda, H., Onogi, K., Kamahori, H., Kobayashi, C., Endo, H. and Miyaoka, K., 2015. The JRA-55 reanalysis: General specifications and basic characteristics.\u202fJournal of the Meteorological Society of Japan. Ser. II,\u202f93(1), pp.5-48. Kobayashi, S., Ota, Y., Harada, Y., Ebita, A., Moriya, M., Onoda, H., Onogi, K., Kamahori, H., Kobayashi, C., Endo, H. and Miyaoka, K., 2015. The JRA-55 reanalysis: General specifications and basic characteristics.\u202fJournal of the Meteorological Society of Japan. Ser. II,\u202f93(1), pp.5-48."]} 
    more » « less
  2. null (Ed.)
    Abstract. This study quantifies differences among four widely usedatmospheric reanalysis datasets (ERA5, JRA-55, MERRA-2, and CFSR) in theirrepresentation of the dynamical changes induced by springtime polarstratospheric ozone depletion in the Southern Hemisphere from 1980 to 2001.The intercomparison is undertaken as part of the SPARC(Stratosphere–troposphere Processes and their Role in Climate) ReanalysisIntercomparison Project (S-RIP). The reanalyses are generally in goodagreement in their representation of the strengthening of the lowerstratospheric polar vortex during the austral spring–summer season,associated with reduced radiative heating due to ozone loss, as well as thedescent of anomalously strong westerly winds into the troposphere duringsummer and the subsequent poleward displacement and intensification of thepolar front jet. Differences in the trends in zonal wind between thereanalyses are generally small compared to the mean trends. The exception isCFSR, which exhibits greater disagreement compared to the other threereanalysis datasets, with stronger westerly winds in the lower stratospherein spring and a larger poleward displacement of the tropospheric westerlyjet in summer. The dynamical changes associated with the ozone hole are examined byinvestigating the momentum budget and then the eddy heat and momentumfluxes in terms of planetary- and synoptic-scale Rossby wave contributions.The dynamical changes are consistently represented across the reanalysesand support our dynamical understanding of the response of the coupledstratosphere–troposphere system to the ozone hole. Although our resultssuggest a high degree of consistency across the four reanalysis datasets inthe representation of these dynamical changes, there are larger differencesin the wave forcing, residual circulation, and eddy propagation changes compared to the zonal wind trends. In particular, there is a noticeabledisparity in these trends in CFSR compared to the other three reanalyses,while the best agreement is found between ERA5 and JRA-55. Greateruncertainty in the components of the momentum budget, as opposed to meancirculation, suggests that the zonal wind is better constrained by theassimilation of observations compared to the wave forcing, residualcirculation, and eddy momentum and heat fluxes, which are more dependent onthe model-based forecasts that can differ between reanalyses. Lookingforward, however, these findings give us confidence that reanalysis datasetscan be used to assess changes associated with the ongoing recovery ofstratospheric ozone. 
    more » « less
  3. Abstract

    Upwelling along ocean eastern boundaries is expected to intensify due to coastal wind strengthening driven by increasing land-sea contrast according to the Bakun hypothesis. Here, the latest high-resolution climate simulations that exhibit drastic improvements of upwelling processes reveal far more complex future upwelling changes. The Southern Hemisphere upwelling systems show a future strengthening in coastal winds with a rapid coastal warming, whereas the Northern Hemisphere coastal winds show a decrease with a comparable warming trend. The Bakun mechanism cannot explain these changes. Heat budget analysis indicates that temperature change in the upwelling region is not simply controlled by vertical Ekman upwelling, but also influenced by horizontal heat advection driven by strong near-coast wind stress curl that is neglected in the Bakun hypothesis and poorly represented by the low-resolution models in the Coupled Model Intercomparison Project. The high-resolution climate simulations also reveal a strong spatial variation in future upwelling changes, which is missing in the low-resolution simulations.

     
    more » « less
  4. Changes in background zonal wind in the tropical Pacific are often invoked to explain changes in ENSO properties. However, the sensitivity of ENSO to mean zonal winds has been thoroughly explored only in intermediate coupled models (following Zebiak and Cane), not in coupled GCMs. The role of mean meridional winds has received even less attention. Accordingly, the goal of this study is to examine systematically the effects of both zonal (equatorial) and meridional (cross-equatorial) background winds on ENSO using targeted experiments with a comprehensive climate model (CESM). Changes in the mean winds are generated by imposing heat flux forcing in two confined regions at a sufficient distance north and south of the equator. We find that the strengthening of either wind component reduces ENSO amplitude, especially eastern Pacific SST variability, and inhibits meridional swings of the intertropical convergence zone (ITCZ). The effect of zonal winds is generally stronger than that of meridional winds. A stability analysis reveals that the strengthening of zonal and meridional winds weakens the ENSO key positive feedbacks, specifically the zonal advection and thermocline feedbacks, which explains these changes. Zonal wind enhancement also intensifies mean upwelling and hence dynamical damping, leading to a further weakening of El Niño events. Ultimately, this study argues that the zonal and, to a lesser extent, meridional wind strengthening of the past decades may have contributed to the observed shift of El Niño characteristics after the year 2000.

     
    more » « less
  5. Cool-season occurrences of blocks, extratropical cyclones that undergo explosive cyclogenesis, and tropical cyclones (TCs) that undergo extratropical transition (ET) from 1980 to 2015 are analyzed using the National Aeronautics and Space Administration’s Modern-Era Retrospective Analysis for Research and Applications, version 2, dataset. These synoptic events are first examined in a climatological analysis that includes identifying consecutive synoptic events, namely, blocks that follow bombs or ET events as well as extratropical cyclones that follow ET events. These synoptic events are then analyzed with respect to three tropical modes of variability: the Madden–Julian oscillation (MJO), El Niño–Southern Oscillation, and the stratospheric quasi-biennial oscillation (QBO). The QBO was considered from both a momentum and thermal point of view, using the equatorial 30-hPa zonal-mean wind and the equatorial zonal wind shear between 30 and 50 hPa, respectively. The results show that in the seven days prior to cool-season blocks and ET events, there is a statistically significant frequency minimum in MJO phases 7 and 3, respectively. With respect to the QBO, there is a statistically significant frequency maximum in neutral QBO conditions during bomb onset and a frequency minimum during ET onset. When stratifying bombs by latitude, there is a significant reduction in Arctic (i.e., poleward of 55°N) bomb onset during easterly QBO conditions. The results show that both tropospheric and stratospheric tropical modes of variability can modulate the frequency of extratropical synoptic events to a similar degree.

     
    more » « less