Operational Global Atmospheric Prediction System Research Articles

Abstract Ocean models need over-ocean atmospheric forcing. However, such forcing is not necessarily provided near the land–sea boundary because 1) the atmospheric model grid used for forcing is frequently much coarser than the ocean model grid, and 2) some of the atmospheric model grid over the ocean includes land values near coastal regions. This paper presents a creeping sea-fill methodology to reduce the improper representation of scalar atmospheric forcing variables near coastal regions, a problem that compromises the usefulness of the fields for ocean model simulations and other offshore applications. For demonstration, atmospheric forcing variables from archived coarse-resolution gridded products—the 1.125° × 1.125° 15-yr European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERA-15) and 1.0° × 1.0° Navy Operational Global Atmospheric Prediction System (NOGAPS)—are used here. A fine-resolution [1/25° × 1/25° cos(lat)], (longitude × latitude) (∼3.2 km) eddy-resolving Black Sea Hybrid Coordinate Ocean Model (HYCOM) is then forced with/without sea-filled atmospheric variables from these gridded products to simulate monthly mean climatological sea surface temperature (SST). Using only over-ocean values from atmospheric forcing fields in the ocean model simulations significantly reduces the climatological mean SST bias (by ∼1°–3°C) and rms SST difference over the seasonal cycle (by ∼2°–3°C) in coastal regions. Performance of the creeping sea-fill methodology is also directly evaluated using measurements of wind speed at 10 m above the surface from the SeaWinds scatterometer on the NASA Quick Scatterometer (QuikSCAT) satellite. Comparisons of original monthly mean wind speeds from operational ECMWF and NOGAPS products with those from QuikSCAT give basin-averaged rms differences of 1.6 and 1.4 m s−1, respectively, during 2000–03. Similar comparisons performed with sea-filled monthly mean wind speeds result in a much lower rms difference (0.7 m s−1 for both products) during the same time period, clearly confirming the accuracy of the methodology even on interannual time scales. Most of the unrealistically low wind speeds from ECMWF and NOGAPS near coastal boundaries are appropriately corrected with the use of the creeping sea fill. Wind speed errors for ECWMF and NOGAPS (mean bias of ≥ 2.5 m s−1 with respect to QuikSCAT during 2000–03) are substantially eliminated (e.g., almost no bias) near most of the land–sea boundaries. Finally, ocean, atmosphere, and coupled atmospheric–oceanic modelers need to be aware that the creeping sea fill is a promising methodology in significantly reducing the land contamination resulting from an improper land–sea mask existing in gridded coarse-resolution atmospheric products (e.g., ECMWF).

Abstract The influence of convective momentum transport (CMT) on tropical cyclone (TC) track forecasts is examined in the Navy Operational Global Atmospheric Prediction System (NOGAPS) with the Emanuel cumulus parameterization. Data assimilation and medium-range forecast experiments show that for 35 tropical cyclones during August and September 2004 the inclusion of CMT in the cumulus parameterization significantly improves the TC track forecasts. The tests show that the track forecasts are very sensitive to the magnitude of the Emanuel parameterization’s convective momentum transport parameter, which controls the CMT tendency returned by the parameterization. While the overall effect of this formulation of CMT in NOGAPS data assimilation/medium-range forecasts results in the surface pressure of tropical cyclones being less intense (and more consistent with the analysis), the parameterization is not equivalent to a simple diffusion of winds in the presence of convection. This is demonstrated by two data assimilation/medium-range forecast tests in which a vertical diffusion algorithm replaces the CMT. Two additional data assimilation/medium-range forecast experiments were conducted to test whether the skill increase primarily comes from the CMT in the immediate vicinity of the tropical cyclones. The results show that the inclusion of the CMT calculation in the vicinity of the TC makes the largest contribution to the increase in forecast skill, but the general contribution of CMT away from the TC also plays an important role.

A multistage projection pursuit (PP) approach is applied to the classification of tropical cyclones (TCs) during extratropical transition (ET) using 500-hPa Navy Operational Global Atmospheric Prediction System (NOGAPS) geopotential height analyses. PP algorithms reduce the dimensionality of the high-dimensional data while minimizing the loss of information that discriminates among classes of ET types. In this paper, a prediction system is developed for ET TC classification and is applied to 85 western North Pacific storms during 1997-2004 using NOGAPS geopotential height analyses to group them as either intensifiers or dissipators. A classification is developed based on the 1997-2002 analyses, and the forecasting performance of the technique is tested on the storms from 2003 and 2004 in two different ways. In the first, the technique is applied at individual times utilizing just spatial information. In the second, the change in spatial patterns with time is taken into account with spatiotemporal algorithms. The preliminary classification results are slightly less accurate than those of NOGAPS but are found to be promising, and further possible improvements are discussed

Abstract Cyclone structure is known to be directly linked to the sensible weather effects produced by the weather system. The extratropical transition (ET) process leads to immense changes in cyclone structure and therefore to changes in the associated weather experienced. Although structure is clearly an important cyclone characteristic, validation of cyclone structure forecasts in operational numerical models has not been previously performed. In this study, short-term (12–36 h) forecasts of cyclone structure from tropical genesis to the completion of ET are validated using fields from the Navy Operational Global Atmospheric Prediction System and the NCEP Aviation model. The cyclone phase space (CPS) is used to quantify differences between forecast and analyzed storm structure, both on a point-by-point basis and through a cyclone-type-based comparison. This cyclone-type comparison exploits a previously defined breakdown of cyclone structure regimes in the CPS. The impacts of synthetic vortex insertion on the ensuing agreement between forecast and analyzed storm structure are explored. While the results show reasonable forecast skill for well-defined (i.e., nonhybrid) systems, cyclones in the process of ET are found to be poorly forecast, emphasizing the need for improved understanding and simulation of the structural changes experienced by ET cyclones.

Abstract The Weber barotropic model (WBAR) was originally developed using predefined 850–200-hPa analyses and forecasts from the NCEP Global Forecasting System. The WBAR tropical cyclone (TC) track forecast performance was found to be competitive with that of more complex numerical weather prediction models in the North Atlantic. As a result, WBAR was revised to incorporate the Navy Operational Global Atmospheric Prediction System (NOGAPS) analyses and forecasts for use at the Joint Typhoon Warning Center (JTWC). The model was also modified to analyze its own storm-dependent deep-layer mean fields from standard NOGAPS pressure levels. Since its operational installation at the JTWC in May 2003, WBAR TC track forecast performance has been competitive with the performance of other more complex NWP models in the western North Pacific. Its TC track forecast performance combined with its high availability rate (93%–95%) has warranted its inclusion in the JTWC operational consensus. The impact of WBAR on consensus TC track forecast performance has been positive and WBAR has added to the consensus forecast availability (i.e., having at least two models to provide a consensus forecast).

Daily inter-annual simulations of SST and MLD using atmospherically forced OGCMs: Model evaluation in comparison to buoy time series

The U. S. Navy Operational Global Atmospheric Prediction System (NOGAPS) daily analysis fields are used to study the African easterly waves to show how they impacted two tropical cyclones. Danielle and Earl (2004) formed approximately at the same time in the Atlantic. A three to eight day band‐pass filter is applied to the daily global analysis fields to reveal the characteristics of the African easterly waves related to the two cyclones. The time‐filtered 850 mb vorticity shows that African easterly waves reaching the formation locations of Danielle and Earl on and about 14 August 2004 bear similar structures, indicating a unique situation in which one train of the easterly wave is responsible for the genesis of two tropical cyclones located one‐wavelength apart. In addition, using piecewise time‐filtered analysis for different latitude bands, it is demonstrated that the movements of both cyclones are significantly modulated by the propagation of the easterly waves, even during Danielle's poleward movement for more than ten‐degrees latitude.

Abstract A high-altitude version of the Navy Operational Global Atmospheric Prediction System (NOGAPS) spectral forecast model is used to simulate the unusual September 2002 Southern Hemisphere stratospheric major warming. Designated as NOGAPS-Advanced Level Physics and High Altitude (NOGAPS-ALPHA), this model extends from the surface to 0.005 hPa (∼85 km altitude) and includes modifications to multiple components of the operational NOGAPS system, including a new radiative heating scheme, middle-atmosphere gravity wave drag parameterizations, hybrid vertical coordinate, upper-level meteorological initialization, and radiatively active prognostic ozone with parameterized photochemistry. NOGAPS-ALPHA forecasts (hindcasts) out to 6 days capture the main features of the major warming, such as the zonal mean wind reversal, planetary-scale wave amplification, large upward Eliassen–Palm (EP) fluxes, and splitting of the polar vortex in the middle stratosphere. Forecasts beyond 6 days have reduced upward EP flux in the lower stratosphere, reduced amplitude of zonal wavenumbers 2 and 3, and a middle stratospheric vortex that does not split. Three-dimensional EP-flux diagnostics in the troposphere reveal that the longer forecasts underestimate upward-propagating planetary wave energy emanating from a significant blocking pattern over the South Atlantic that played a large role in forcing the major warming. Forecasts of less than 6 days are initialized with the blocking in place, and therefore are not required to predict the blocking onset. For a more thorough skill assessment, NOGAPS-ALPHA forecasts over 3 weeks during September–October 2002 are compared with operational NOGAPS 5-day forecasts made at the time. NOGAPS-ALPHA forecasts initialized with 2002 operational NOGAPS analyses show a modest improvement in skill over the NOGAPS operational forecasts. An additional, larger improvement is obtained when NOGAPS-ALPHA is initialized with reanalyzed 2002 fields produced with the currently operational (as of October 2003) Naval Research Laboratory (NRL) Atmospheric Variational Data Assimilation System (NAVDAS). Thus the combination of higher model top, better physical parameterizations, and better initial conditions all yield improved forecasting skill over the NOGAPS forecasts issued operationally at the time.

Abstract A 34-member ensemble-mean trajectory through the cyclone phase space (CPS) is calculated using Navy Operational Global Atmospheric Prediction System (NOGAPS) analyses for North Atlantic tropical cyclones (TCs) undergoing extratropical transition (ET). Synoptic composites at four ET milestones are examined: 24 h prior to the beginning of ET (TB − 24), the beginning of ET (TB), the end of ET (TE), and 24 h after the end of ET (TE + 24). While the extratropically transitioning TC structure is tightly constrained in its tropical phase, it has a variety of evolutions after TE. Partitioning the ensemble based upon post-ET intensity change or structure discriminates among statistically significant ET precursor conditions. Compositing the various post-ET intensity regimes provides insight into the important environmental factors governing post-ET development. A TC that intensifies (weakens) after TE begins transition (t = TB) with a negatively (positively) tilted trough 1000 km (1500 km) upstream. The negative tilt permits a contraction and intensification of the eddy potential vorticity (PV) flux, while the positive trough tilt prevents contraction and intensification of the forcing. In 6 of the 34 cases, the posttropical cold-core cyclone develops a warm-seclusion structure, rather than remaining cold core. Anticipation of this warm-seclusion evolution is critical since it represents a dramatically increased risk of middle- to high-latitude wind and wave damage. The warm-seclusion evolution is most favored when the scale of the interacting trough closely matches the scale of the transitioned TC, focusing the eddy PV flux in the outflow layer of the transitioning TC. The sensitivity of structural evolution prior to and after TE illustrated here gives insight into the degradation of global model midlatitude forecast accuracy during a pending ET event. Eliassen–Palm flux cross sections suggest that ET is primarily driven by the eddy angular momentum flux of the trough, rather than the eddy heat flux associated with the trough. The response of the transitioning TC to the eddy angular momentum forcing is to produce adiabatic ascent and cooling radially inward and beneath the region of the forcing to restore thermal wind balance. In the case of ET, the forcing is maximized lower in the atmosphere, and spread over much greater depth, than in the case of trough-induced TC intensification. Only after TE is the eddy heat flux forcing as significant as the eddy angular momentum forcing, further supporting a physical foundation for the CPS description of cyclone evolution.

Abstract This study describes atmospheric forcing parameters constructed from different global climatologies, applied to the Black Sea, and investigates the sensitivity of Hybrid Coordinate Ocean Model (HYCOM) simulations to these products. Significant discussion is devoted to construction of these parameters before using them in the eddy-resolving (≈3.2-km resolution) HYCOM simulations. The main goal is to answer how the model dynamics can be substantially affected by different atmospheric forcing products in the Black Sea. Eight wind forcing products are used: four obtained from observation-based climatologies, including one based on measurements from the SeaWinds scatterometer on the Quick Scatterometer (QuikSCAT) satellite, and the rest formed from operational model products. Thermal forcing parameters, including solar radiation, are formed from two operational models: the European Centre for Medium-Range Weather Forecasts (ECMWF) and the Fleet Numerical Meteorology and Oceanography Center (FNMOC) Navy Operational Global Atmospheric Prediction System (NOGAPS). Climatologically forced Black Sea HYCOM simulations (without ocean data assimilation) are then performed to assess the accuracy and sensitivity of the model sea surface temperature (SST) and sea surface circulation to these wind and thermal forcing products. Results demonstrate that the model-simulated SST structure is quite sensitive to the wind and thermal forcing products, especially near coastal regions. Despite this sensitivity, several robust features are found in the model SST in comparison to a monthly 9.3-km-resolution satellite-based Pathfinder SST climatology. Annual mean HYCOM SST usually agreed to within ≈±0.2° of the climatology in the interior of the Black Sea for any of the wind and thermal forcing products used. The fine-resolution (0.25° × 0.25°) wind forcing from the scatterometer data along with thermal forcing from NOGAPS gave the best SST simulation with a basin-averaged rms difference value of 1.21°C, especially improving model results near coastal regions. Specifically, atmospherically forced model simulations with no assimilation of any ocean data suggest that the basin-averaged rms SST differences with respect to the Pathfinder SST climatology can vary from 1.21° to 2.15°C depending on the wind and thermal forcing product. The latter rms SST difference value is obtained when using wind forcing from the National Centers for Environmental Prediction (NCEP), a product that has a too-coarse grid resolution of 1.875° × 1.875° for a small ocean basin such as the Black Sea. This paper also highlights the importance of using high-frequency (hybrid) wind forcing as opposed to monthly mean wind forcing in the model simulations. Finally, there are large variations in the annual mean surface circulation simulated using the different wind sets, with general agreement between those forced by the model-based products (vector correlation is usually >0.7). Three of the observation-based climatologies generally yield unrealistic circulation features and currents that are too weak.

Abstract An adjoint-based method is used to calculate the impact of observation data on a measure of short-range forecast error in the Navy Operational Global Atmospheric Prediction System (NOGAPS) during November and December 2003. The evaluated observations include all regular satellite and in situ data assimilated in the Naval Research Laboratory (NRL) Atmospheric Variational Data Assimilation System (NAVDAS) at 1800 UTC, and also targeted dropsonde profiles provided by the North Atlantic Observing-System Research and Predictability Experiment (THORPEX) Regional Campaign (NA-TReC) field program. Commerical aircraft observations account for 46% of the total forecast error reduction by observations in the NA-TReC domain, which includes the North Atlantic and adjacent regions of North America and Europe. Targeted dropsonde data have high impact per observation, but the impact of all dropsonde data is less than 2% of the total during the 2-month study period. Eight of 12 targeted dropsonde cases reduce forecast error. The percent of total impact for other observations assimilated at 1800 UTC in the NA-TReC domain is as follows: Advanced Microwave Sounding Unit-A (AMSU-A) radiances (16%), satellite winds (14%), land surface data (9%), radiosondes (8%), and ship-surface data (5%). If observations over the entire global domain are evaluated, the largest impact of data provided at 1800 UTC during November and December 2003 is provided by AMSU-A radiance data (48% of total).

Abstract The Kim–Arakawa orographic gravity‐wave drag parametrization scheme, which is a component of the US Navy's NOGAPS ALPHA (Navy Operational Global Atmospheric Prediction System, Advanced‐Level Physics and High Altitude), is extended to include the effects of orographic anisotropy and low‐level flow blocking. The algorithms to calculate the orographic statistics needed for the parametrization are also revised. The extended scheme is evaluated against mountain waves explicitly simulated with COAMPS®† (Coupled Ocean/Atmosphere Mesoscale Prediction System) of NRL (Naval Research Laboratory).Mountain‐wave simulations over Boulder, Colorado, USA, are used for representing realistic situations of different wave activity including severe downslope windstorms. The simulations are area‐averaged and interpolated to the vertical grid of NOGAPS, and are used as the input to the extended Kim–Arakawa scheme. The scheme is calibrated by comparing the parametrized vertical distribution of the momentum fluxes with the counterpart obtained from the explicit mesoscale simulations. Overall, the calibrated scheme successfully represents the simulated magnitudes and vertical divergences of the momentum fluxes. A flow regime diagram is constructed utilizing a time series of the simulations to further evaluate the parametrization. The robustness of the orographic statistics, together with an approximate method to improve it, are also addressed. Copyright © 2005 Royal Meteorological Society.

Iron mines in Minnesota are ideally located to assess the accuracy of available atmospheric profiles used in infrasound modeling. These mines are located approximately 400 km away to the southeast (142) of the Lac-Du-Bonnet infrasound station, IS-10. Infrasound data from June 1999 to March 2004 was analyzed to assess the effects of explosion size and atmospheric conditions on observations. IS-10 recorded a suite of events from this time period resulting in well constrained ground truth. This ground truth allows for the comparison of ray trace and PE (Parabolic Equation) modeling to the observed arrivals. The tele-infrasonic distance (greater than 250 km) produces ray paths that turn in the upper atmosphere, the thermosphere, at approximately 120 km to 140 km. Modeling based upon MSIS/HWM (Mass Spectrometer Incoherent Scatter/Horizontal Wind Model) and the NOGAPS (Navy Operational Global Atmospheric Prediction System) and NRL-GS2 (Naval Research Laboratory Ground to Space) augmented profiles are used to interpret the observed arrivals.

Abstract This paper examines the sensitivity of sea surface temperature (SST) to water turbidity in the Black Sea using the eddy-resolving (∼3.2-km resolution) Hybrid Coordinate Ocean Model (HYCOM), which includes a nonslab K-profile parameterization (KPP) mixed layer model. The KPP model uses a diffusive attenuation coefficient of photosynthetically active radiation (kPAR) processed from a remotely sensed dataset to take water turbidity into account. Six model experiments (expt) are performed with no assimilation of any ocean data and wind/thermal forcing from two sources: 1) European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis (ERA) and 2) Fleet Numerical Meteorology and Oceanography Center (FNMOC) Navy Operational Global Atmospheric Prediction System (NOGAPS). Forced with ECMWF, experiment 1 uses spatially and monthly varying kPAR values over the Black Sea, experiment 2 assumes all of the solar radiation is absorbed at the sea surface, and experiment 3 uses a constant kPAR value of 0.06 m−1, representing clear-water constant solar attenuation depth of 16.7 m. Experiments 4, 5, and 6 are twins of 1, 2, and 3 but forced with NOGAPS. The monthly averaged model SSTs resulting from all experiments are then compared with a fine-resolution (∼9 km) satellite-based monthly SST climatology (the Pathfinder climatology). Because of the high turbidity in the Black Sea, it is found that a clear-water constant attenuation depth (i.e., expts 3 and 6) results in SST bias as large as 3°C in comparison with standard simulations (expts 1 and 4) over most of the Black Sea in summer. In particular, when using the clear-water constant attenuation depth as opposed to using spatial and temporal kPAR, basin-averaged rms SST difference with respect to the Pathfinder SST climatology increases ∼46% (from 1.41°C in expt 1 to 2.06°C in expt 3) in the ECMWF forcing case. Similarly, basin-averaged rms SST difference increases ∼36% (from 1.39°C in expt 4 to 1.89°C in expt 6) in the NOGAPS forcing case. The standard HYCOM simulations (expts 1 and 4) have a very high basin-averaged skill score of 0.95, showing overall model success in predicting climatological SST, even with no assimilation of any SST data. In general, the use of spatially and temporally varying turbidity fields is necessary for the Black Sea OGCM studies because there is strong seasonal cycle and large spatial variation in the solar attenuation coefficient, and an additional simulation using a constant kPAR value of 0.19 m−1, the Sea-Viewing Wide Field-of-View Sensor (SeaWiFS) space–time mean for the Black Sea, did not yield as accurate SST results as experiments 1 and 4. Model–data comparisons also revealed that relatively large HYCOM SST errors close to the coastal boundaries can be attributed to the misrepresentation of land– sea mask in the ECMWF and NOGAPS products. With the relatively accurate mask used in NOGAPS, HYCOM demonstrated the ability to simulate accurate SSTs in shallow water over the broad northwest shelf in the Black Sea, a region of large errors using the inaccurate mask in ECMWF. A linear relationship is found between changes in SST and changes in heat flux below the mixed layer. Specifically, a change of ∼50 W m−2 in sub-mixed-layer heat flux results in a SST change of ∼3.0°C, a value that occurs when using clear-water constant attenuation depth rather than monthly varying kPAR in the model simulations, clearly demonstrating potential impact of penetrating solar radiation on SST simulations.

Abstract The extratropical transition (ET) of tropical cyclones in the North Atlantic basin is explored using the framework of the Cyclone Phase Space (CPS). Analyses from the U.S. Navy Operational Global Atmospheric Prediction System (NOGAPS) and National Centers for Environmental Prediction Aviation Model (AVN) are used to diagnose CPS locations of 19 Atlantic tropical cyclones that underwent ET in the years 1998–2002. A nonhierarchical cluster analysis technique, the k-means method, is employed to determine physically meaningful clusters of cyclone structure in this CPS dataset. The cluster analysis successfully partitions observations into subsets representing distinct cyclone structures. It is also found that, while current operational cyclone classification schemes are successful at defining some storm structures, cluster analysis assists in defining a “gray area” of transitioning tropical cyclones, an area where operational schemes overlap. A mean path through the CPS of tropical cyclones undergoin...

Abstract A ⅛° global version of the Navy Coastal Ocean Model (NCOM), operational at the Naval Oceanographic Office (NAVOCEANO), is used for prediction of sea surface height (SSH) on daily and monthly time scales during 1998–2001. Model simulations that use 3-hourly wind and thermal forcing obtained from the Navy Operational Global Atmospheric Prediction System (NOGAPS) are performed with/without data assimilation to examine indirect/direct effects of atmospheric forcing in predicting SSH. Model–data evaluations are performed using the extensive database of daily averaged SSH values from tide gauges in the Atlantic, Pacific, and Indian Oceans obtained from the Joint Archive for Sea Level (JASL) center during 1998–2001. Model–data comparisons are based on observations from 282 tide gauge locations. An inverse barometer correction was applied to SSH time series from tide gauges for model–data comparisons, and a sensitivity study is undertaken to assess the impact of the inverse barometer correction on the SSH validation. A set of statistical metrics that includes conditional bias (Bcond), root-mean-square (rms) difference, correlation coefficient (R), and nondimensional skill score (SS) is used to evaluate the model performance. It is shown that global NCOM has skill in representing SSH even in a free-running simulation, with general improvement when SSH from satellite altimetry and sea surface temperature (SST) from satellite IR are assimilated via synthetic temperature and salinity profiles derived from climatological correlations. When the model was run from 1998 to 2001 with NOGAPS forcing, daily model SSH comparisons from 612 yearlong daily tide gauge time series gave a median rms difference of 5.98 cm (5.77 cm), an R value of 0.72 (0.76), and an SS value of 0.45 (0.51) for the ⅛° free-running (assimilative) NCOM. Similarly, error statistics based on the 30-day running averages of SSH time series for 591 yearlong daily tide gauge time series over the time frame 1998–2001 give a median rms difference of 3.63 cm (3.36 cm), an R value of 0.83 (0.85), and an SS value of 0.60 (0.64) for the ⅛° free-running (assimilated) NCOM. Model– data comparisons show that skill in 30-day running average SSH time series is as much as 30% higher than skill for daily SSH. Finally, SSH predictions from the free-running and assimilative ⅛° NCOM simulations are validated against sea level data from the tide gauges in two different ways: 1) using original detided sea level time series from tide gauges and 2) using the detided data with an inverse barometer correction derived using daily mean sea level pressure extracted from NOGAPS at each location. Based on comparisons with 612 yearlong daily tide gauge time series during 1998–2001, the inverse barometer correction lowered the median rms difference by about 1 cm (15%–20%). Results presented in this paper reveal that NCOM is able to predict SSH with reasonable accuracies, as evidenced by model simulations performed during 1998–2001. In an extension of the validation over broader ocean regions, the authors find good agreement in amplitude and distribution of SSH variability between NCOM and other operational model products.

Various drag mechanisms are currently parametrized in numerical models of the atmosphere. For global models that include the middle atmosphere in particular, these mechanisms profoundly affect weather forecast as well as climate simulation.We have developed an extended-top version of the Navy Operational Global Atmospheric Prediction System (NOGAPS) to include the middle atmosphere by modifying some physical parametrizations and the vertical coordinate. We performed a series of ensemble simulations corresponding to January 2000 for investigating the response of the model to various drag mechanisms, such as mountain drag, orographic gravity wave drag, surface friction drag, and artificial model top drag. Based on the monthly mean fields obtained from the simulations, we first investigate the effect of gravity wave drag due to its direct impact through planetary wave activity as well as indirect impact through induced meridional circulation.We discuss the difficulties in partitioning between the mountain drag due to resolved orography and the gravity wave drag due to unresolved orography, first using conventional diagnostic measures. From analyses of the atmospheric angular momentum budget, we show that various model drag mechanisms when modified interact with one another by redistributing their drag while conserving the total amount. In particular, an overestimation of mountain drag is accompanied by an underestimation of gravity wave drag in the Northern Hemisphere mid-latitudes to conserve the total amount of drag in the model while likely breaking an optimal balance among the mechanisms. Under such a condition, the inclusion of a gravity wave drag parametrization — even if the drag amount itself is reasonable – does not necessarily improve the performance of the model. Diagnosis of this type of imbalance is not clear by conventional monthly mean fields of variables. In this paper, we argue that the budget of atmospheric angular momentum is a useful measure to diagnose impact of such changes in model physics with regard to the partition and balance among drag mechanisms. We also discuss the experimental results that led to the replacement of silhouette orography by mean orography in our model.

A hybrid Lagrangian–Eulerian model for calculating the trajectories of near-surface drifters in the ocean is developed in this study. The model employs climatological, near-surface currents computed from a spline fit of all available drifter velocities observed in the Pacific Ocean between 1988 and 1996. It also incorporates contemporaneous wind fields calculated by either the U.S. Navy [the Navy Operational Global Atmospheric Prediction System (NOGAPS)] or the European Centre for Medium-Range Weather Forecasts (ECMWF). The model was applied to 30 drifters launched in the tropical Pacific Ocean in three clusters during 1990, 1993, and 1994. For 10-day-long trajectories the forecasts computed by the hybrid model are up to 164% closer to the observed trajectories compared to the trajectories obtained by advecting the drifters with the climatological currents only. The best-fitting trajectories are computed with ECMWF fields that have a temporal resolution of 6 h. The average improvement over all 30 drifters of the hybrid model trajectories relative to advection by the climatological currents is 21%, but in the open-ocean clusters (1990 and 1993) the improvement is 42% with ECMWF winds (34% with NOGAPS winds). This difference between the open-ocean and coastal clusters is due to the fact that the model does not presently include the effect of horizontal boundaries (coastlines). For zero initial velocities the trajectories generated by the hybrid model are significantly more accurate than advection by the mean currents on time scales of 5–15 days. For 3-day-long trajectories significant improvement is achieved if the drifter's initial velocity is known, in which case the model-generated trajectories are about 2 times closer to observations than persistence. The model's success in providing more accurate trajectories indicates that drifters' motion can deviate significantly from the climatological current and that the instantaneous winds are more relevant to their trajectories than the mean surface currents. It also demonstrates the importance of an accurate initial velocity, especially for short trajectories on the order of 1–3 days. A possible interpretation of these results is that winds affect drifter motion more than the water velocity since drifters do not obey continuity.

Abstract The convective parameterization of Emanuel has been employed in the forecast model of the Navy Operational Global Atmospheric Prediction System (NOGAPS) since 2000, when it replaced a version of the relaxed Arakawa–Schubert scheme. Although in long-period data assimilation forecast tests the Emanuel scheme has been found to perform quite well in NOGAPS, particularly for tropical cyclones, some weaknesses have also become apparent. These weaknesses include underprediction of heavy-precipitation events, too much light precipitation, and unrealistic heating at upper levels. Recent research efforts have resulted in modifications of the scheme that are designed to reduce such problems. One change described here involves the partitioning of the cloud-base mass flux into mixing cloud mass flux at individual levels. The new treatment significantly reduces a heating anomaly near the tropopause that is associated with a large amount of mixing cloud mass flux ascribed to that region in the original Emanuel ...

We performed a systematic comparison of wind speed measurements from the SeaWinds QuikSCAT scatterometer and wind speeds computed from RADARSAT-1 synthetic aperture radar (SAR) normalized radar cross section measurements. These comparisons were made over in the Gulf of Alaska and extended over a two-year period, 2000 and 2001. The SAR wind speed estimates require a wind direction to initialize the retrieval. Here, we first used wind directions from the Navy Operational Global Atmospheric Prediction System (NOGAPS) model. For these retrievals, the standard deviation between the resulting SAR and QuikSCAT wind speed measurements was 1.78 m/s. When we used the QuikSCAT-measured wind directions to initialize the inversion, comparisons improve to a standard deviation of 1.36 m/s. We used these SAR-scatterometer comparisons to generate a new C-band horizontal polarization model function. With this new model function, the wind speed inversion improves to a standard deviation of 1.24 m/s with no mean bias. These results strongly suggest that SAR and QuikSCAT measurements can be combined to make better high-resolution wind measurements than either instrument could alone in coastal areas.