Global comparisons of barotropic and internal tides generated in an eddy‐resolving ocean circulation model are made with tidal estimates obtained from altimetric sea surface heights and an ...altimetry‐constrained tide model. As far as we know, our Hybrid Coordinate Ocean Model (HYCOM) simulations shown here and in an earlier paper are the only published high‐resolution global simulations to contain barotropic tides, internal tides, the general circulation, and mesoscale eddies concurrently. Comparing the model barotropic tide with a global data‐assimilative shallow water tide model shows that the global tidal elevation differences are approximately evenly split between discrepancies in tidal amplitude and phase. Both the model and observations show strong generation of internal tides at a limited number of “hot spot” regions with propagation of beams of energy for thousands of kilometers away from the sources. The model internal tidal amplitudes compare well with observations near these energetic tidal regions. Averaged over these regions, the model and observation internal tide amplitude estimates agree to approximately 15% for the four largest semidiurnal constituents and 23% for the four largest diurnal constituents. Away from the hot spots, the comparison between the model and altimetric amplitude is not as good due, in part, to two problems, errors in the model barotropic tides and overestimation of the altimetric tides in regions of strong mesoscale eddy activity. Examining the general energy distribution of the simulated internal tide is an important first step in the evaluation of internal tides in HYCOM.
Key Points
Model and observations show generation of internal tides in limited regions
Our global model is able to generate internal waves consistent with observations
Barotropic phase errors are a major source of errors in the model internal tide
This study presents the improvement in ice edge error within the US Navy's operational sea ice forecast systems gained by assimilating high horizontal resolution satellite-derived ice concentration ...products. Since the late 1980's, the ice forecast systems have assimilated near real-time sea ice concentration derived from the Defense Meteorological Satellite Program (DMSP) Special Sensor Microwave/Imager (SSMI and then SSMIS). The resolution of the satellite-derived product was approximately the same as the previous operational ice forecast system (25 km). As the sea ice forecast model resolution increased over time, the need for higher horizontal resolution observational data grew. In 2013, a new Navy sea ice forecast system (Arctic Cap Nowcast/Forecast System – ACNFS) went into operations with a horizontal resolution of ~ 3.5 km at the North Pole. A method of blending ice concentration observations from the Advanced Microwave Scanning Radiometer (AMSR2) along with a sea ice mask produced by the National Ice Center (NIC) has been developed, resulting in an ice concentration product with very high spatial resolution. In this study, ACNFS was initialized with this newly developed high resolution blended ice concentration product. The daily ice edge locations from model hindcast simulations were compared against independent observed ice edge locations. ACNFS initialized using the high resolution blended ice concentration data product decreased predicted ice edge location error compared to the operational system that only assimilated SSMIS data. A second evaluation assimilating the new blended sea ice concentration product into the pre-operational Navy Global Ocean Forecast System 3.1 also showed a substantial improvement in ice edge location over a system using the SSMIS sea ice concentration product alone. This paper describes the technique used to create the blended sea ice concentration product and the significant improvements in ice edge forecasting in both of the Navy's sea ice forecasting systems.
Spatial and temporal variability of the impact of air‐sea stratification on the differences between satellite‐derived 10 m equivalent neutral wind speeds and stability‐dependent (e.g., in situ) 10 m ...wind speeds are quantitatively examined over the global ocean. The influences of stability are compared with three air‐sea flux algorithms, Coupled Ocean‐Atmosphere Response Experiment (version 3.0), Bourassa‐Vincent‐Wood, and Liu‐Katsaros‐Businger. Analyses are first presented at many individual buoy locations and then are extended to the global ocean with the use of rain‐free wind measurements from the SeaWinds scatterometer on the QuikSCAT satellite, gridded at a resolution of 0.25° × 0.25°. Overall, stability‐dependent winds are found to be weaker than equivalent neutral winds by 0.2 m s−1 on the basis of 7619 monthly mean values from 208 buoys during 2000–2005. Differences based on hourly winds can be as large as ±0.5 m s−1. Results remain robust regardless of which air‐sea flux algorithm is used. Monthly rain‐free gridded QuikSCAT measurements, combined with atmospheric stability determined using near‐surface variables from the European Centre for Medium‐Range Weather Forecasts 40‐year reanalysis, demonstrate the effects of stratification on the 10 m winds globally. Differences in stability‐dependent and neutral winds are substantially nonsymmetrical and reveal locations where the former is stronger than the latter. These differences may cause physically significant biases in air‐sea fluxes if they are not properly considered, especially near the Kuroshio and Gulf Stream current systems.
This study introduces exchange coefficients for wind stress (CD ), latent heat flux (CL ), and sensible heat flux (CS ) over the global ocean. They are obtained from the state-of-the-art Coupled ...Ocean-Atmosphere Response Experiment (COARE) bulk algorithm (version 3.0). Using the exchange coefficients from this bulk scheme, CD , CL , and CS are then expressed as simple polynomial functions of air-sea temperature difference (Ta - Ts )-where air temperature (Ta ) is at 10 m, wind speed (Va ) is at 10 m, and relative humidity (RH) is at the air-sea interface-to parameterize stability. The advantage of using polynomial-based exchange coefficients is that they do not require any iterations for stability. In addition, they agree with results from the COARE algorithm but at -5 times lower computation cost, an advantage that is particularly needed for ocean general circulation models (OGCMs) and climate models running at high horizontal resolution and short time steps. The effects of any water vapor flux in calculating the exchange coefficients are taken into account in the polynomial functions, a feature that is especially important at low wind speeds (e.g., Va < 5 m s-1) because air-sea mixing ratio difference can have a major effect on the stability, particularly in tropical regions. Analyses of exchange coefficients demonstrate the fact that water vapor can have substantial impact on air-sea exchange coefficients at low wind speeds. An example application of the exchange coefficients from the polynomial approach is the recalculation of climatological mean wind stress magnitude from 40-yr European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERA-40) data in the North Pacific Ocean over 1979-2002. Using ECMWF 10-m winds and the authors' methodology provides accurate surface stresses while largely eliminating the orographically induced Gibb's waves found in the original ERA-40 surface wind stresses. These can have a large amplitude near mountainous regions and can extend far into the ocean interior. This study introduces exchange coefficients of air-sea fluxes, which are applicable to the wide range of conditions occurring over the global ocean, including the air-sea stability differences across the Gulf Stream and Kuroshio, regions which have been the subject of many climate model studies. This versatility results because CD , CL , and CS are determined for Va values of 1 to 40 m s-1, (Ta - Ts ), intervals of -8' to 7'C, and RH values of 0% to 100%. Exchange coefficients presented here are called the Naval Research Laboratory (NRL) Air-Sea Exchange Coefficients (NASEC) and they are suitable for a wide range of air-sea interaction studies and model applications.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, UILJ, UKNU, UL, UM, UPUK
This article provides an overview of the effort centered on the HYbrid Coordinate Ocean Model (HYCOM) to develop an eddy-resolving, real-time global and basin-scale ocean hindcast, nowcast, and ...prediction system in the context of the Global Ocean Data Assimilation Experiment (GODAE). The main characteristics of HYCOM are first presented, followed by a description and assessment of the present near real-time Atlantic forecasting system. Regional/coastal applications are also discussed since an important attribute of the data assimilative HYCOM simulations is the capability to provide boundary conditions to regional and coastal models. The final section describes the steps taken toward the establishment of the fully global eddy-resolving HYCOM data assimilative system and discusses some of the difficulties associated with advanced data assimilation given the size of the problem.
The accuracy of wind speed at 10 m above the sea surface from two satellite and three numerical weather prediction (NWP) products is investigated over the global ocean. Rain‐free equivalent neutral ...winds from the Quick Scatterometer (QuikSCAT) are converted to stability‐dependent winds to be consistent with those from NWP products and are taken as truth in comparisons to winds from other products. Quantitative statistical analyses presented at each grid point over the global ocean reveal that monthly winds from NWP products have almost perfect skill relative to those from QuikSCAT winds during the 3‐year common period (September 1999 to August 2002). Exceptions occur in tropical regions and high southern latitudes. Wind speeds adjusted to 10 m at many moored buoys located in different regions of the global ocean further confirm the accuracy of monthly NWP winds, giving RMS difference of 1.0 m s−1 based on 1281 monthlong time series. The satellite‐based QuikSCAT winds agree with buoy winds relatively better than NWP products. While there is good agreement among wind products on monthly timescales, large differences (>3 m s−1 and more) in NWP winds are found in comparison to QuikSCAT winds on shorter time intervals at high latitudes. Daily means of sensible and latent heat fluxes based on NWP winds can therefore differ as much as 100 W m−2 in comparison to those based on QuikSCAT winds. In general, NWP wind‐based sensible and latent heat fluxes are more similar to their QuikSCAT wind‐based counterparts in tropical regions and midlatitudes.
Global maps of the mesoscale eddy available potential energy (EAPE) field at a depth of 500 m are created using potential density anomalies in a high‐resolution 1/12.5° global ocean model. Maps made ...from both a free‐running simulation and a data‐assimilative reanalysis of the HYbrid Coordinate Ocean Model (HYCOM) are compared with maps made by other researchers from density anomalies in Argo profiles. The HYCOM and Argo maps display similar features, especially in the dominance of western boundary currents. The reanalysis maps match the Argo maps more closely, demonstrating the added value of data assimilation. Global averages of the simulation, reanalysis, and Argo EAPE all agree to within about 10%. The model and Argo EAPE fields are compared to EAPE computed from temperature anomalies in a data set of “moored historical observations” (MHO) in conjunction with buoyancy frequencies computed from a global climatology. The MHO data set allows for an estimate of the EAPE in high‐frequency motions that is aliased into the Argo EAPE values. At MHO locations, 15–32% of the EAPE in the Argo estimates is due to aliased motions having periods of 10 days or less. Spatial averages of EAPE in HYCOM, Argo, and MHO data agree to within 50% at MHO locations, with both model estimates lying within error bars observations. Analysis of the EAPE field in an idealized model, in conjunction with published theory, suggests that much of the scatter seen in comparisons of different EAPE estimates is to be expected given the chaotic, unpredictable nature of mesoscale eddies.
Key Points
Global maps of the mesoscale eddy available potential energy are made from a HYCOM simulation and reanalysis
Modeled eddy available potential energy compares well to Argo observations globally, and to moored instruments locally
Model‐data comparisons of eddy available potential energy exhibit intrinsic scatter
A set of numerical simulations is used to investigate the Pacific Ocean circulation north of 20°S, with emphasis on the Kuroshio/Oyashio current system. The primitive equation models used for these ...simulations have a free surface and realistic geometry that includes the deep marginal seas, such as the Sea of Japan. Most of the simulations have 1/8° resolution for each variable but range from 1/2°, 1.5‐layer reduced gravity to 1/16°, six layer with realistic bottom topography. These are used to investigate the dynamics of the Kuroshio/Oyashio current system and to identify the processes that contribute most to the realism of the simulations. This is done by model‐data comparisons, by using the modularity of layered ocean models to include/exclude certain dynamical processes, by varying the model geometry and bottom topography, and by varying model parameters, such as horizontal grid resolution, layer structure, and eddy viscosity. In comparison with observational data, the simulations show that the barotropic mode, at least one internal mode, nonlinearity, high “horizontal” resolution (1/8° or finer), the regional bottom topography, and the wind forcing are critical for realistic simulations. The first four are important for baroclinic instability (eddy‐mean energetics actually show mixed barotropic‐baroclinic instability), the wind curl pattern for the formation and basic placement of the current system, and the bottom topography for the distribution of the instability and for influences on the pathways of the mean flow. Both the Hellerman and Rosenstein (1983) (HR) monthly wind stress climatology and 1000‐mbar winds from the European Centre for Medium‐Range Weather Forecasts (ECMWF) have been used to drive the model. East of about 150°E, they give a mean latitude for the Kuroshio Extension that differs by about 3°, approximately 34°N for HR, 37°N for ECMWF, and 35°N observed. The subarctic front is the northern boundary of the subtropical gyre. It is associated with the annual and April–September mean zero wind stress curl lines (which are similar), while the Kuroshio Extension is associated with wintertime zero wind stress curl. This means that part of the flow from the Kuroshio must pass north of the Kuroshio Extension and connect with the Oyashio and subarctic front. Realistic routes for this connection are flow through the Sea of Japan, a nonlinear route separated from the east coast of Japan, and bifurcation of the Kuroshio at the Shatsky Rise. In addition, the six‐layer simulations show a 3‐Sv meridional overturning cell with southward surface flow and northward return flow centered near 400 m depth. Baroclinic instability plays a critical role in coupling the shallow and abyssal layer circulations and in allowing the bottom topography to strongly influence the shallow circulation. By this means, the Izu Ridge and Trench and seamounts upstream and downstream of these have profound influence on (1) the mean path of the Kuroshio and its mean meanders south and east of Japan and (2) on separating the northward flow connecting the Kuroshio and the Oyashio/subarctic front from the east coast of Japan. Without the topographic influence, the models show an unrealistic northward current along the east coast of Japan. In essence, the topography regulates the location and strength of the baroclinic instability. The baroclinic instability gives eddy‐driven deep mean flows that follow the f/h contours (where f is the Coriolis parameter and h is the depth of the water column) of the bottom topography. These abyssal currents then strongly influence the pathway for subtropical gyre flow north of the Kuroshio Extension and steer the mean meanders in the Kuroshio south and east of Japan. This is corroborated by current meter data from the Kuroshio Extension Regional Experiment (World Ocean Circulation Experiment line PCM 7). The meander path south of Japan depends on the occurrence of baroclinic instability west of the Izu Ridge; otherwise, a straight path occurs. The pathway shows little sensitivity to the Tokara Strait transport over the range simulated (36–72 Sv in yearly means). However, interannual increases in wind forcing or Tokara Strait transport give rise to a predominant meander path, while decreases yield a predominant straight path. Resolution of 1/8° in an ocean model is comparable to the 2.5° resolution used in atmospheric forecast models in the early 1980s based on the first internal mode Rossby radius of deformation. Model comparisons at 1/8° and 1/16° resolution and comparisons with current meter data and Geosat altimeter data show that 1/16° resolution is needed for adequate eastward penetration of the high eddy kinetic energy associated with the Kuroshio Extension.
The ocean west of the main Hawaiian Islands is characterized by enhanced eddy kinetic energy arising from the abundance of locally generated mesoscale eddies, most frequently in the area west of the ...island of Hawaii. Two mechanisms of eddy generation in the wake of an island are examined with numerical model experiments. The first, eddy generation and shedding by an oceanic flow around an oceanic barrier, requires the existence of strong westward flows to the north and south of the island of Hawaii. Model solutions show such westward flows and generation of eddies by these flows although the intensity of the eddies and the generation frequency are much lower than that derived from altimetry. As a result, these eddies contribute an insignificant amount of eddy kinetic energy in the region. The second, eddy generation and shedding by an atmospheric flow around an atmospheric barrier, is based on oceanic upwelling and downwelling induced by surface wind shear, effectively introducing sinks and sources to the ocean interior. Previous idealized modeling studies have shown that oceanic eddies can be generated by sufficiently strong forcing (source or sink), providing an explanation why eddy occurrences in the lee of the island of Hawaii coincide with periods of strong trade winds. Eddy generation characteristics in the model experiments are consistent with this mechanism in terms of time of occurrence, strength and the resulting eddy kinetic energy. Major discrepancies are in eddy propagation and therefore eddy distribution in the regional domain due to the complex nature of eddy‐eddy interactions.
Key Points
Eddies in the lee of the Hawaiian Islands are primarily wind‐driven
Eddies can also be induced by oceanic flows, but the resulting EKE is low
Eddy‐eddy interactions greatly influence eddy propagation pathways
Through a comprehensive analysis, reliability of 10 m wind speeds is presented near the land‐sea boundaries over the global ocean. Winds from three numerical weather prediction (NWP) centers and two ...satellite‐based products are analyzed. NWP products are 1.875° × 1.875° National Center Environmental Prediction reanalyses, 1.125° × 1.125° European Centre for Medium‐Range Weather Forecasts 40‐year Reanalysis (ERA‐40), and 1.0° × 1.0° Navy Operational Global Atmospheric Prediction System (NOGAPS) operational product. These are compared to much finer resolution (0.25° × 0.25°) satellite winds, Quick Scatterometer (QSCAT) and Special Sensor Microwave/Imager. Large biases (e.g., >3 m s−1) may exist in NWP products near the land‐sea boundaries, because wind speeds from the uniformly gridded global fields are generally at a spatial scale too coarse to appropriately define the contrast between water and land grid points. This so‐called land contamination of ocean‐only winds varies, and typically depends on the extent of the land‐sea mask. A creeping sea‐fill methodology is introduced to reduce errors in winds. It is based on the elimination of land‐corrupted NWP grid points and replacement by adjacent, purely over‐ocean values. In comparison to winds from many moored buoys, the methodology diminishes RMS errors (from >4 m s−1 to <1 m s−1) for NOGAPS and ERA‐40. The creeping sea‐fill is not advised for NCEP winds which have low contrast between land and sea points, thereby resulting in little impact from the land contamination.