skip to main content


Title: Monthly and Annual contour lines of the zero and the positive maximum of the Wind Stress Curl over Western North Atlantic during 1980-2019 and the Gulf Stream path during 1993-2019.

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
Award ID(s):
2123283 1851242
NSF-PAR ID:
10448887
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Zenodo
Date Published:
Edition / Version:
1.0.0
Subject(s) / Keyword(s):
["Ocean-Air Interaction","Gulf Stream","Force and Response","Ocean Circulation","Wind Stress Curl","Altimetry","Zero Wind Stress Curl","Maximum Wind Stress Curl","Gulf Stream Path"]
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. This dataset contains monthly average output files from the iCAM6 simulations used in the manuscript "Enhancing understanding of the hydrological cycle via pairing of process-oriented and isotope ratio tracers," in review at the Journal of Advances in Modeling Earth Systems. A file corresponding to each of the tagged and isotopic variables used in this manuscript is included. Files are at 0.9° latitude x 1.25° longitude, and are in NetCDF format. Data from two simulations are included: 1) a simulation where the atmospheric model was "nudged" to ERA5 wind and surface pressure fields, by adding an additional tendency (see section 3.1 of associated manuscript), and 2) a simulation where the atmospheric state was allowed to freely evolve, using only boundary conditions imposed at the surface and top of atmosphere. Specific information about each of the variables provided is located in the "usage notes" section below. Associated article abstract: The hydrologic cycle couples the Earth's energy and carbon budgets through evaporation, moisture transport, and precipitation. Despite a wealth of observations and models, fundamental limitations remain in our capacity to deduce even the most basic properties of the hydrological cycle, including the spatial pattern of the residence time (RT) of water in the atmosphere and the mean distance traveled from evaporation sources to precipitation sinks. Meanwhile, geochemical tracers such as stable water isotope ratios provide a tool to probe hydrological processes, yet their interpretation remains equivocal despite several decades of use. As a result, there is a need for new mechanistic tools that link variations in water isotope ratios to underlying hydrological processes. Here we present a new suite of “process-oriented tags,” which we use to explicitly trace hydrological processes within the isotopically enabled Community Atmosphere Model, version 6 (iCAM6). Using these tags, we test the hypotheses that precipitation isotope ratios respond to parcel rainout, variations in atmospheric RT, and preserve information regarding meteorological conditions during evaporation. We present results for a historical simulation from 1980 to 2004, forced with winds from the ERA5 reanalysis. We find strong evidence that precipitation isotope ratios record information about atmospheric rainout and meteorological conditions during evaporation, but little evidence that precipitation isotope ratios vary with water vapor RT. These new tracer methods will enable more robust linkages between observations of isotope ratios in the modern hydrologic cycle or proxies of past terrestrial environments and the environmental processes underlying these observations.   Details about the simulation setup can be found in section 3 of the associated open-source manuscript, "Enhancing understanding of the hydrological cycle via pairing of process‐oriented and isotope ratio tracers." In brief, we conducted two simulations of the atmosphere from 1980-2004 using the isotope-enabled version of the Community Atmosphere Model 6 (iCAM6) at 0.9x1.25° horizontal resolution, and with 30 vertical hybrid layers spanning from the surface to ~3 hPa. In the first simulation, wind and surface pressure fields were "nudged" toward the ERA5 reanalysis dataset by adding a nudging tendency, preventing the model from diverging from observed/reanalysis wind fields. In the second simulation, no additional nudging tendency was included, and the model was allowed to evolve 'freely' with only boundary conditions provided at the top (e.g., incoming solar radiation) and bottom (e.g., observed sea surface temperatures) of the model. In addition to the isotopic variables, our simulation included a suite of 'process-oriented tracers,' which we describe in section 2 of the manuscript. These variables are meant to track a property of water associated with evaporation, condensation, or atmospheric transport. Metadata are provided about each of the files below; moreover, since the attached files are NetCDF data - this information is also provided with the data files. NetCDF metadata can be accessed using standard tools (e.g., ncdump). Each file has 4 variables: the tagged quantity, and the associated coordinate variables (time, latitude, longitude). The latter three are identical across all files, only the tagged quantity changes. Twelve files are provided for the nudged simulation, and an additional three are provided for the free simulations: Nudged simulation files iCAM6_nudged_1980-2004_mon_RHevap: Mass-weighted mean evaporation source property: RH (%) with respect to surface temperature. iCAM6_nudged_1980-2004_mon_Tevap: Mass-weighted mean evaporation source property: surface temperature in Kelvin iCAM6_nudged_1980-2004_mon_Tcond: Mass-weighted mean condensation property: temperature (K) iCAM6_nudged_1980-2004_mon_columnQ: Total (vertically integrated) precipitable water (kg/m2).  Not a tagged quantity, but necessary to calculate depletion times in section 4.3 (e.g., Fig. 11 and 12). iCAM6_nudged_1980-2004_mon_d18O: Precipitation d18O (‰ VSMOW) iCAM6_nudged_1980-2004_mon_d18Oevap_0: Mass-weighted mean evaporation source property - d18O of the evaporative flux (e.g., the 'initial' isotope ratio prior to condensation), (‰ VSMOW) iCAM6_nudged_1980-2004_mon_dxs: Precipitation deuterium excess (‰ VSMOW) - note that precipitation d2H can be calculated from this file and the precipitation d18O as d2H = d-excess - 8*d18O. iCAM6_nudged_1980-2004_mon_dexevap_0: Mass-weighted mean evaporation source property - deuterium excess of the evaporative flux iCAM6_nudged_1980-2004_mon_lnf: Integrated property - ln(f) calculated from the constant-fractionation d18O tracer (see section 3.2). iCAM6_nudged_1980-2004_mon_precip: Total precipitation rate in m/s. Note there is an error in the metadata in this file - it is total precipitation, not just convective precipitation. iCAM6_nudged_1980-2004_mon_residencetime: Mean atmospheric water residence time (in days). iCAM6_nudged_1980-2004_mon_transportdistance: Mean atmospheric water transport distance (in km). Free simulation files iCAM6_free_1980-2004_mon_d18O: Precipitation d18O (‰ VSMOW) iCAM6_free_1980-2004_mon_dxs: Precipitation deuterium excess (‰ VSMOW) - note that precipitation d2H can be calculated from this file and the precipitation d18O as d2H = d-excess - 8*d18O. iCAM6_free_1980-2004_mon_precip: Total precipitation rate in m/s. Note there is an error in the metadata in this file - it is total precipitation, not just convective precipitation. 
    more » « less
  3. This dataset consists of weekly trajectory information of Gulf Stream Warm Core Rings from 2000-2010. This work builds upon Silver et al. (2022a) ( https://doi.org/10.5281/zenodo.6436380) which contained Warm Core Ring trajectory information from 2011 to 2020. Combining the two datasets a total of 21 years of weekly Warm Core Ring trajectories can be obtained. An example of how to use such a dataset can be found in Silver et al. (2022b).

    The format of the dataset is similar to that of  Silver et al. (2022a), and the following description is adapted from their dataset. This dataset is comprised of individual files containing each ring’s weekly center location and its area for 374 WCRs present between January 1, 2000 and December 31, 2010. Each Warm Core Ring is identified by a unique alphanumeric code 'WEyyyymmddA', where 'WE' represents a Warm Eddy (as identified in the analysis charts); 'yyyymmdd' is the year, month and day of formation; and the last character 'A' represents the sequential sighting of the eddies in a particular year. Continuity of a ring which passes from one year to the next is maintained by the same character in the first sighting.  For example, the first ring in 2002 having a trailing alphabet of 'F' indicates that five rings were carried over from 2001 which were still observed on January 1, 2002. Each ring has its own netCDF (.nc) filename following its alphanumeric code. Each file contains 4 variables, “Lon”- the ring center’s weekly longitude, “Lat”- the ring center’s weekly latitude, “Area” - the rings weekly size in km2, and “Date” in days - representing the days since Jan 01, 0000. 

    The process of creating the WCR tracking dataset follows the same methodology of the previously generated WCR census (Gangopadhyay et al., 2019, 2020). The Jenifer Clark’s Gulf Stream Charts used to create this dataset are 2-3 times a week from 2000-2010. Thus, we used approximately 1560 Charts for the 10 years of analysis. All of these charts were reanalyzed between 75° and 55°W using QGIS 2.18.16 (2016) and geo-referenced on a WGS84 coordinate system (Decker, 1986). 

     

    Silver, A., Gangopadhyay, A, & Gawarkiewicz, G. (2022a). Warm Core Ring Trajectories in the Northwest Atlantic Slope Sea (2011-2020) (1.0.0) [Data set]. Zenodo. https://doi.org/10.5281/zenodo.6436380

    Silver, A., Gangopadhyay, A., Gawarkiewicz, G., Andres, M., Flierl, G., & Clark, J. (2022b). Spatial Variability of Movement, Structure, and Formation of Warm Core Rings in the Northwest Atlantic Slope Sea. Journal of Geophysical Research: Oceans127(8), e2022JC018737. https://doi.org/10.1029/2022JC018737 

    Gangopadhyay, A., G. Gawarkiewicz, N. Etige, M. Monim and J. Clark, 2019. An Observed Regime Shift in the Formation of Warm Core Rings from the Gulf Stream, Nature - Scientific Reports, https://doi.org/10.1038/s41598-019-48661-9. www.nature.com/articles/s41598-019-48661-9.

    Gangopadhyay, A., N. Etige, G. Gawarkiewicz, A. M. Silver, M. Monim and J. Clark, 2020.  A Census of the Warm Core Rings of the Gulf Stream (1980-2017). Journal of Geophysical Research, Oceans, 125, e2019JC016033. https://doi.org/10.1029/2019JC016033.

    QGIS Development Team. QGIS Geographic Information System (2016).

    Decker, B. L. World Geodetic System 1984. World geodetic system 1984 (1986).

     

    Funded by two NSF US grants OCE-1851242, OCE-212328 {"references": ["Silver, A., Gangopadhyay, A, & Gawarkiewicz, G. (2022). Warm Core Ring Trajectories in the Northwest Atlantic Slope Sea (2011-2020) (1.0.0) [Data set]. Zenodo. https://doi.org/10.5281/zenodo.6436380", "Silver, A., Gangopadhyay, A., Gawarkiewicz, G., Andres, M., Flierl, G., & Clark, J. (2022b). Spatial Variability of Movement, Structure, and Formation of Warm Core Rings in the Northwest Atlantic Slope Sea.\u00a0Journal of Geophysical Research: Oceans,\u00a0127(8), e2022JC018737.\u00a0https://doi.org/10.1029/2022JC018737", "Gangopadhyay, A., G. Gawarkiewicz, N. Etige, M. Monim and J. Clark, 2019. An Observed Regime Shift in the Formation of Warm Core Rings from the Gulf Stream, Nature - Scientific Reports, https://doi.org/10.1038/s41598-019-48661-9. www.nature.com/articles/s41598-019-48661-9.", "Gangopadhyay, A., N. Etige, G. Gawarkiewicz, A. M. Silver, M. Monim and J. Clark, 2020. A Census of the Warm Core Rings of the Gulf Stream (1980-2017). Journal of Geophysical Research, Oceans, 125, e2019JC016033. https://doi.org/10.1029/2019JC016033.", "QGIS Development Team. QGIS Geographic Information System (2016).", "Decker, B. L. World Geodetic System 1984. World geodetic system 1984 (1986)."]} 
    more » « less
  4. Abstract

    Global warming is causing an unprecedented loss of species and habitats worldwide. This is particularly apparent for tropical coral reefs, with an increasing number of reefs experiencing mass bleaching and mortality on an annual basis. As such, there is a growing need for a standardized experimental approach to rapidly assess the thermal limits of corals and predict the survival of coral species across reefs and regions. Using a portable experimental system, the Coral Bleaching Automated Stress System (CBASS), we conducted standardized 18 h acute thermal stress assays to quantitively determine the upper thermal limits of four coral species across the length of the Red Sea coastline, from the Gulf of Aqaba (GoA) to Djibouti (~ 2100 km). We measured dark-acclimated photosynthetic efficiency (Fv/Fm), algal symbiont density, chlorophyll a, and visual bleaching intensity following heat stress.Fv/Fmwas the most precise response variable assessed, advancing theFv/Fmeffective dose 50 (ED50, i.e., the temperature at which 50% of the initialFv/Fmis measured) as an empirically derived proxy for thermal tolerance. ED50 thermal thresholds from the central/southern Red Sea and Djibouti populations were consistently higher forAcropora hemprichii, Pocillopora verrucosa,andStylophora pistillata(0.1–1.8 °C above GoA corals, respectively), in line with prevailing warmer maximum monthly means (MMMs), though were lower than GoA corals relative to site MMMs (1.5–3.0 °C).P. verrucosahad the lowest thresholds overall. Despite coming from the hottest site, thresholds were lowest forPorites lobatain the southern Red Sea, suggesting long-term physiological damage or ongoing recovery from a severe, prior bleaching event. Altogether, the CBASS resolved historical, taxonomic, and possibly recent environmental drivers of variation in coral thermal thresholds, highlighting the potential for a standardized, short-term thermal assay as a universal approach for assessing ecological and evolutionary variation in the upper thermal limits of corals.

     
    more » « less
  5. This dataset contains the atmospheric river catalogues and the associated precipitation and temperature data for the Preindustrial and Last Glacial Maximum CESM2 simulations presented in the GRL manuscript:  Atmospheric river contributions to ice sheet hydro climate at the Last Glacial Maximum. The atmospheric river catalogue files (zipped) are in netcdf format and organized by year. There are 100 years of data for both simulations.  The Preindustrial simulation catalogue begins in model year 41 and ends in model year 140.  The LGM simulation catalogue begins in model year 1 and ends in year 100. Each yearly file has a temporal resolution of 6 hours (1460 time steps each file) and a spatial resolution of 0.9° x 1.25° (the native resolution of the CESM simulation). A variable in the file called "ar_binary_tag" indicates whether an atmospheric river is present at each grid cell and each tilmestep: 1 indicates an atmospheric river is present; 0 indicates an atmospheric river is not present.  The precipitation and temperature files are 100-year annual or 100-year seasonal averages of atmospheric river precipitation/temperature. See the Methods section of the article for more details on the atmospheric river detection algorithm and precipitation/temperature calculations.

    Associated article abstract:

    Atmospheric rivers (ARs) are an important driver of surface mass balance over today’s Greenland and Antarctic ice sheets. Using paleoclimate simulations with the Community Earth System Model, we find ARs also had a key influence on the extensive ice sheets of the Last Glacial Maximum (LGM). ARs provide up to 53% of total precipitation along the margins of the eastern Laurentide ice sheet and up to 22-27% of precipitation along the margins of the Patagonian, western Cordilleran, and western Fennoscandian ice sheets. Despite overall cold conditions at the LGM, surface temperatures during AR events are often above freezing, resulting in more rain than snow along ice sheet margins and conditions that promote surface melt. The results suggest  ARs may have had an important role in ice sheet growth and melt during previous glacial periods and may have accelerated ice sheet retreat following the LGM.

     
    more » « less