Comparison of ADCP and ECCOv4r4 Currents in the Pacific Equatorial UndercurrentAbstract
The Pacific Equatorial Undercurrent (EUC) flows eastward across the Pacific at the equator in the thermocline. Its variability is related to El Niño. Moored acoustic Doppler current profiler (ADCP) measurements recorded at four widely separated sites along the equator in the EUC were compared to currents generated by version 4 release 4 of the Estimating the Circulation and Climate of the Ocean (ECCOv4r4) global model–data synthesis product. We are interested to learn how well ECCOv4r4 currents could complement sparse in situ current measurements. ADCP measurements were not assimilated in ECCOv4r4. Comparisons occurred at 5-m depth intervals at 165°E, 170°W, 140°W, and 110°W over time intervals of 10–14 years from 1995 to 2010. Hourly values of ECCOv4r4 and ADCP EUC core speeds were strongly correlated, similar for the EUC transport per unit width (TPUW). Correlations were substantially weaker at 110°W. Although we expected means and standard deviations of ECCOv4r4 currents to be smaller than ADCP values because of ECCOv4r4’s grid representation error, the large differences were unforeseen. The appearance of ECCOv4r4 diurnal-period current oscillations was surprising. As the EUC moved eastward from 170° to 140°W, the ECCOv4r4 TPUW exhibited a much smaller increase compared to the ADCP TPUW. A consequence of smaller ECCOv4r4 EUC core speeds was significantly fewer instances of gradient Richardson number (Ri) less than 1/4 above and below the depth of the core speed compared to Ri computed with ADCP observations. We present linear regression analyses to use monthly-mean ECCOv4r4 EUC core speeds and TPUWs as proxies for ADCP measurements.
Significance Statement
Hundreds of scientific papers have used ECCO data products generated with a continually evolving state-of-the-art ocean-model–data synthesis system. We ask, How representative is the latest version of ECCO equatorial ocean currents? We use independent in situ current measurements as the reference dataset to establish the accuracy of ECCO currents in the tropical Pacific. Attention is focused on the Pacific Equatorial Undercurrent (EUC) because it contributes to the formation of El Niño and La Niña events. ECCO EUC core speeds were smaller in magnitude and less variable in time compared to observations. As a consequence, ECCO currents generated smaller vertical mixing in the EUC compared to that inferred from current measurements. We developed a linear regression model to improve representation of monthly-mean ECCO currents.
Abstract
The Pacific Equatorial Undercurrent (EUC) flows eastward across the Pacific at the equator in the thermocline. Its variability is related to El Niño. Moored acoustic Doppler current profiler (ADCP) measurements recorded at four widely separated sites along the equator in the EUC were compared to currents generated by version 4 release 4 of the Estimating the Circulation and Climate of the Ocean (ECCOv4r4) global model–data synthesis product. We are interested to learn how well ECCOv4r4 currents could complement sparse in situ current measurements. ADCP measurements were not assimilated in ECCOv4r4. Comparisons occurred at 5-m depth intervals at 165°E, 170°W, 140°W, and 110°W over time intervals of 10–14 years from 1995 to 2010. Hourly values of ECCOv4r4 and ADCP EUC core speeds were strongly correlated, similar for the EUC transport per unit width (TPUW). Correlations were substantially weaker at 110°W. Although we expected means and standard deviations of ECCOv4r4 currents to be smaller than ADCP values because of ECCOv4r4’s grid representation error, the large differences were unforeseen. The appearance of ECCOv4r4 diurnal-period current oscillations was surprising. As the EUC moved eastward from 170° to 140°W, the ECCOv4r4 TPUW exhibited a much smaller increase compared to the ADCP TPUW. A consequence of smaller ECCOv4r4 EUC core speeds was significantly fewer instances of gradient Richardson number (Ri) less than 1/4 above and below the depth of the core speed compared to Ri computed with ADCP observations. We present linear regression analyses to use monthly-mean ECCOv4r4 EUC core speeds and TPUWs as proxies for ADCP measurements.
Significance Statement
Hundreds of scientific papers have used ECCO data products generated with a continually evolving state-of-the-art ocean-model–data synthesis system. We ask, How representative is the latest version of ECCO equatorial ocean currents? We use independent in situ current measurements as the reference dataset to establish the accuracy of ECCO currents in the tropical Pacific. Attention is focused on the Pacific Equatorial Undercurrent (EUC) because it contributes to the formation of El Niño and La Niña events. ECCO EUC core speeds were smaller in magnitude and less variable in time compared to observations. As a consequence, ECCO currents generated smaller vertical mixing in the EUC compared to that inferred from current measurements. We developed a linear regression model to improve representation of monthly-mean ECCO currents.
Abstract
The influence of meridional shift of the oceanic subtropical front (STF) on the Agulhas Current (AC) regime shifts is studied using satellite altimeter data and a 1.5-layer ocean model. The satellite observations suggest the northward shift of the STF leads to the AC leaping across the gap with little Agulhas leakage, and the southward shift of the STF mainly results in the AC intruding into the Atlantic Ocean in the forms of a loop current and an eddy-shedding path, while there are three flow patterns of AC for moderate latitude of the STF. The ocean model results suggest no hysteresis (associated with multiple equilibrium states) exists in the AC system. The model reproduces similar AC regimes depending on different gap widths as in the observations, and model results can be used to explain the observed Agulhas leakage well. We also present the parameter space of the critical AC strength that results in different AC flow patterns as a function of the gap width. The vorticity dynamics of the AC regime shift suggests that the β term is mainly balanced by the viscosity term for the AC in the leaping and loop current paths, while the β and instantaneous vorticity terms are mainly balanced by the advection and viscosity terms for the AC in the eddy-shedding path. These findings help explain the dynamics of the AC flowing across the gateway beyond the tip of Africa affected by the north–south shift of the STF in the leaping regime or penetrating regime.
Abstract
The influence of meridional shift of the oceanic subtropical front (STF) on the Agulhas Current (AC) regime shifts is studied using satellite altimeter data and a 1.5-layer ocean model. The satellite observations suggest the northward shift of the STF leads to the AC leaping across the gap with little Agulhas leakage, and the southward shift of the STF mainly results in the AC intruding into the Atlantic Ocean in the forms of a loop current and an eddy-shedding path, while there are three flow patterns of AC for moderate latitude of the STF. The ocean model results suggest no hysteresis (associated with multiple equilibrium states) exists in the AC system. The model reproduces similar AC regimes depending on different gap widths as in the observations, and model results can be used to explain the observed Agulhas leakage well. We also present the parameter space of the critical AC strength that results in different AC flow patterns as a function of the gap width. The vorticity dynamics of the AC regime shift suggests that the β term is mainly balanced by the viscosity term for the AC in the leaping and loop current paths, while the β and instantaneous vorticity terms are mainly balanced by the advection and viscosity terms for the AC in the eddy-shedding path. These findings help explain the dynamics of the AC flowing across the gateway beyond the tip of Africa affected by the north–south shift of the STF in the leaping regime or penetrating regime.
Abstract
Streamwise vorticity currents (SVCs) have been hypothesized to enhance low-level mesocyclones within supercell thunderstorms and perhaps increase the likelihood of tornadogenesis. Recent observational studies have confirmed the existence of SVCs in supercells and numerical simulations have allowed for further investigation of SVCs. A suite of 19 idealized supercell simulations with varying midlevel shear orientations is analyzed to determine how SVC formation and characteristics may differ between storms. In our simulations, SVCs develop on the cold side of left-flank convergence boundaries and their updraft-relative positions are partially dependent on downdraft location. The magnitude, duration, and mean depth of SVCs do not differ significantly between simulations or between SVCs that precede tornado-like vortices (TLVs) and those that do not. Trajectories initialized within SVCs reveal two primary airstreams, one that flows through an SVC for the majority of its length, and another that originates in the modified inflow in the forward flank and then merges with the SVC. Vorticity budgets calculated along trajectories reveal that the first airstream exhibits significantly greater maximum streamwise vorticity magnitudes than the second airstream. The vorticity budgets also indicate that stretching of horizontal streamwise vorticity is the dominant contributor to the large values of streamwise vorticity within the SVCs. TLV formation does not require the development of an SVC beforehand; 44% of TLVs in the simulations are preceded by SVCs. When an SVC occurs, it is followed by a TLV 53% of the time, indicating not all SVCs lead to TLV formation.
Significance Statement
Streamwise vorticity currents (SVCs) are features within thunderstorms hypothesized to strengthen updraft rotation and increase the likelihood of tornado formation. SVCs in a suite of 19 thunderstorm simulations are analyzed to investigate how they develop, if their characteristics differ between storms, and how often they precede tornado production. The rotation in an SVC is amplified as air accelerates toward the updraft, which is the main process contributing to SVC formation. The likelihood of SVCs may vary with differences in the winds 3–6 km above the ground. These findings may aid in developing strategies for better observing SVCs.
Abstract
Streamwise vorticity currents (SVCs) have been hypothesized to enhance low-level mesocyclones within supercell thunderstorms and perhaps increase the likelihood of tornadogenesis. Recent observational studies have confirmed the existence of SVCs in supercells and numerical simulations have allowed for further investigation of SVCs. A suite of 19 idealized supercell simulations with varying midlevel shear orientations is analyzed to determine how SVC formation and characteristics may differ between storms. In our simulations, SVCs develop on the cold side of left-flank convergence boundaries and their updraft-relative positions are partially dependent on downdraft location. The magnitude, duration, and mean depth of SVCs do not differ significantly between simulations or between SVCs that precede tornado-like vortices (TLVs) and those that do not. Trajectories initialized within SVCs reveal two primary airstreams, one that flows through an SVC for the majority of its length, and another that originates in the modified inflow in the forward flank and then merges with the SVC. Vorticity budgets calculated along trajectories reveal that the first airstream exhibits significantly greater maximum streamwise vorticity magnitudes than the second airstream. The vorticity budgets also indicate that stretching of horizontal streamwise vorticity is the dominant contributor to the large values of streamwise vorticity within the SVCs. TLV formation does not require the development of an SVC beforehand; 44% of TLVs in the simulations are preceded by SVCs. When an SVC occurs, it is followed by a TLV 53% of the time, indicating not all SVCs lead to TLV formation.
Significance Statement
Streamwise vorticity currents (SVCs) are features within thunderstorms hypothesized to strengthen updraft rotation and increase the likelihood of tornado formation. SVCs in a suite of 19 thunderstorm simulations are analyzed to investigate how they develop, if their characteristics differ between storms, and how often they precede tornado production. The rotation in an SVC is amplified as air accelerates toward the updraft, which is the main process contributing to SVC formation. The likelihood of SVCs may vary with differences in the winds 3–6 km above the ground. These findings may aid in developing strategies for better observing SVCs.
Abstract
Doppler-lidar wind-profile measurements at three sites were used to evaluate NWP model errors from two versions of NOAA’s 3-km-grid HRRR model, to see whether updates in the latest version 4 reduced errors when compared against the original version 1. Nested (750-m grid) versions of each were also tested to see how grid spacing affected forecast skill. The measurements were part of the field phase of the Second Wind Forecasting Improvement Project (WFIP2), an 18-month deployment into central Oregon–Washington, a major wind-energy-producing region. This study focuses on errors in simulating marine intrusions, a summertime, 600–800-m-deep, regional sea-breeze flow found to generate large errors. HRRR errors proved to be complex and site dependent. The most prominent error resulted from a premature drop in modeled marine-intrusion wind speeds after local midnight, when lidar-measured winds of greater than 8 m s−1 persisted through the next morning. These large negative errors were offset at low levels by positive errors due to excessive mixing, complicating the interpretation of model “improvement,” such that the updates to the full-scale versions produced mixed results, sometimes enhancing but sometimes degrading model skill. Nesting consistently improved model performance, with version 1’s nest producing the smallest errors overall. HRRR’s ability to represent the stages of sea-breeze forcing was evaluated using radiation budget, surface-energy balance, and near-surface temperature measurements available during WFIP2. The significant site-to-site differences in model error and the complex nature of these errors mean that field-measurement campaigns having dense arrays of profiling sensors are necessary to properly diagnose and characterize model errors, as part of a systematic approach to NWP model improvement.
Significance Statement
Dramatic increases in NWP model skill will be required over the coming decades. This paper describes the role of major deployments of accurate profiling sensors in achieving that goal and presents an example from the Second Wind Forecast Improvement Program (WFIP2). Wind-profile data from scanning Doppler lidars were used to evaluate two versions of HRRR, the original and an updated version, and nested versions of each. This study focuses on the ability of updated HRRR versions to improve upon predicting a regional sea-breeze flow, which was found to generate large errors by the original HRRR. Updates to the full-scale HRRR versions produced mixed results, but the finer-mesh versions consistently reduced model errors.
Abstract
Doppler-lidar wind-profile measurements at three sites were used to evaluate NWP model errors from two versions of NOAA’s 3-km-grid HRRR model, to see whether updates in the latest version 4 reduced errors when compared against the original version 1. Nested (750-m grid) versions of each were also tested to see how grid spacing affected forecast skill. The measurements were part of the field phase of the Second Wind Forecasting Improvement Project (WFIP2), an 18-month deployment into central Oregon–Washington, a major wind-energy-producing region. This study focuses on errors in simulating marine intrusions, a summertime, 600–800-m-deep, regional sea-breeze flow found to generate large errors. HRRR errors proved to be complex and site dependent. The most prominent error resulted from a premature drop in modeled marine-intrusion wind speeds after local midnight, when lidar-measured winds of greater than 8 m s−1 persisted through the next morning. These large negative errors were offset at low levels by positive errors due to excessive mixing, complicating the interpretation of model “improvement,” such that the updates to the full-scale versions produced mixed results, sometimes enhancing but sometimes degrading model skill. Nesting consistently improved model performance, with version 1’s nest producing the smallest errors overall. HRRR’s ability to represent the stages of sea-breeze forcing was evaluated using radiation budget, surface-energy balance, and near-surface temperature measurements available during WFIP2. The significant site-to-site differences in model error and the complex nature of these errors mean that field-measurement campaigns having dense arrays of profiling sensors are necessary to properly diagnose and characterize model errors, as part of a systematic approach to NWP model improvement.
Significance Statement
Dramatic increases in NWP model skill will be required over the coming decades. This paper describes the role of major deployments of accurate profiling sensors in achieving that goal and presents an example from the Second Wind Forecast Improvement Program (WFIP2). Wind-profile data from scanning Doppler lidars were used to evaluate two versions of HRRR, the original and an updated version, and nested versions of each. This study focuses on the ability of updated HRRR versions to improve upon predicting a regional sea-breeze flow, which was found to generate large errors by the original HRRR. Updates to the full-scale HRRR versions produced mixed results, but the finer-mesh versions consistently reduced model errors.
Abstract
Air–sea coupling system in the southwestern Indian Ocean (SWIO; 35°–55°S, 40°–75°E) exhibits predominant multidecadal variability that is the strongest during austral summer. It is characterized by an equivalent barotropic atmospheric high (low) pressure over warm (cold) sea surface temperature (SST) anomalies and a poleward (equatorward) shift of the westerlies during the positive (negative) phase. In this study, physical processes of this multidecadal variability are investigated by using observations/reanalysis and CMIP6 model simulations. Results suggest that the multidecadal fluctuation can be explained by the modulation of the Atlantic meridional overturning circulation (AMOC) and the local air–sea positive feedback in the SWIO. In both observations/reanalysis and CMIP6 model simulations, the AMOC fluctuation presents a significantly negative correlation with the multidecadal SST variation in the SWIO when the AMOC is leading by about a decade. The mechanisms are that the preceding AMOC variation can cause an interhemispheric dipolar pattern of SST anomalies in the Atlantic Ocean. Subsequently, the SST anomalies in the midlatitudes of the South Atlantic can propagate to the SWIO by the oceanic Rossby wave under the influence of the Antarctic Circumpolar Current (ACC). Once the SST anomalies reach the SWIO, these SST anomalies in the oceanic front can affect the baroclinicity in the lower troposphere to influence the synoptic transient eddy and then cause the atmospheric circulation anomaly via the eddy–mean flow interaction. Subsequently, the anomalous atmospheric circulation over the SWIO can significantly strengthen the SST anomalies through modifying the oceanic meridional temperature advection and latent and sensible heat flux.
Abstract
Air–sea coupling system in the southwestern Indian Ocean (SWIO; 35°–55°S, 40°–75°E) exhibits predominant multidecadal variability that is the strongest during austral summer. It is characterized by an equivalent barotropic atmospheric high (low) pressure over warm (cold) sea surface temperature (SST) anomalies and a poleward (equatorward) shift of the westerlies during the positive (negative) phase. In this study, physical processes of this multidecadal variability are investigated by using observations/reanalysis and CMIP6 model simulations. Results suggest that the multidecadal fluctuation can be explained by the modulation of the Atlantic meridional overturning circulation (AMOC) and the local air–sea positive feedback in the SWIO. In both observations/reanalysis and CMIP6 model simulations, the AMOC fluctuation presents a significantly negative correlation with the multidecadal SST variation in the SWIO when the AMOC is leading by about a decade. The mechanisms are that the preceding AMOC variation can cause an interhemispheric dipolar pattern of SST anomalies in the Atlantic Ocean. Subsequently, the SST anomalies in the midlatitudes of the South Atlantic can propagate to the SWIO by the oceanic Rossby wave under the influence of the Antarctic Circumpolar Current (ACC). Once the SST anomalies reach the SWIO, these SST anomalies in the oceanic front can affect the baroclinicity in the lower troposphere to influence the synoptic transient eddy and then cause the atmospheric circulation anomaly via the eddy–mean flow interaction. Subsequently, the anomalous atmospheric circulation over the SWIO can significantly strengthen the SST anomalies through modifying the oceanic meridional temperature advection and latent and sensible heat flux.
Abstract
Recurving tropical cyclones (TCs) in the western North Pacific often cause heavy rainfall events (HREs) in East Asia. However, how their interactions with midlatitude flows alter the characteristics of HREs remains unclear. The present study examines the synoptic–dynamic characteristics of HREs directly resulting from TCs in South Korea with a focus on the role of midlatitude baroclinic condition. The HREs are categorized into two clusters based on midlatitude tropopause patterns: strongly (C1) and weakly (C2) baroclinic conditions. C1, which is common in late summer, is characterized by a well-defined trough–ridge couplet and jet streak at the tropopause. As TCs approach, the trough–ridge couplet amplifies, but is anchored by divergent TC outflow. This leads to phase locking of the upstream trough with TCs and thereby prompts substantial structural changes of TCs reminiscent of extratropical transition. The synergistic TC–midlatitude flow interactions allow for widely enhanced quasigeostrophic forcing for ascent to the north of the TC center. This allows HREs to occur even before TC landfall with more inland rainfall than C2 HREs. In contrast, C2, which is mainly observed in midsummer, does not accompany the undulating tropopause. In the absence of strong interactions with midlatitude flows, TCs rapidly dissipate after HREs while maintaining their tropical features. The upward motion is confined to the inherent TC convection, and thus HREs occur only when TCs are located in the vicinity of the country. These findings suggest that midlatitude baroclinic condition determines the spatial extent of TC rainfall and the timing of TC-induced HREs in South Korea.
Significance Statement
This study suggests that the midlatitude flows can substantially modulate heavy rainfall events directly caused by tropical cyclones. By analyzing the 42-yr tropical cyclone–induced heavy rainfall events in South Korea, it is found that tropical cyclones and midlatitude flows strongly interact with each other, especially when the midlatitude flows meander in conjunction with a strong jet stream. Their synergistic interactions result in a poleward expansion of the tropical cyclones’ precipitation shields, leading to heavy rainfall events even before they make landfall in the country. Consequently, it is advisable to carefully monitor the midlatitude conditions as well as tropical cyclones themselves as earlier heavy rainfall warnings may be necessary depending on the former.
Abstract
Recurving tropical cyclones (TCs) in the western North Pacific often cause heavy rainfall events (HREs) in East Asia. However, how their interactions with midlatitude flows alter the characteristics of HREs remains unclear. The present study examines the synoptic–dynamic characteristics of HREs directly resulting from TCs in South Korea with a focus on the role of midlatitude baroclinic condition. The HREs are categorized into two clusters based on midlatitude tropopause patterns: strongly (C1) and weakly (C2) baroclinic conditions. C1, which is common in late summer, is characterized by a well-defined trough–ridge couplet and jet streak at the tropopause. As TCs approach, the trough–ridge couplet amplifies, but is anchored by divergent TC outflow. This leads to phase locking of the upstream trough with TCs and thereby prompts substantial structural changes of TCs reminiscent of extratropical transition. The synergistic TC–midlatitude flow interactions allow for widely enhanced quasigeostrophic forcing for ascent to the north of the TC center. This allows HREs to occur even before TC landfall with more inland rainfall than C2 HREs. In contrast, C2, which is mainly observed in midsummer, does not accompany the undulating tropopause. In the absence of strong interactions with midlatitude flows, TCs rapidly dissipate after HREs while maintaining their tropical features. The upward motion is confined to the inherent TC convection, and thus HREs occur only when TCs are located in the vicinity of the country. These findings suggest that midlatitude baroclinic condition determines the spatial extent of TC rainfall and the timing of TC-induced HREs in South Korea.
Significance Statement
This study suggests that the midlatitude flows can substantially modulate heavy rainfall events directly caused by tropical cyclones. By analyzing the 42-yr tropical cyclone–induced heavy rainfall events in South Korea, it is found that tropical cyclones and midlatitude flows strongly interact with each other, especially when the midlatitude flows meander in conjunction with a strong jet stream. Their synergistic interactions result in a poleward expansion of the tropical cyclones’ precipitation shields, leading to heavy rainfall events even before they make landfall in the country. Consequently, it is advisable to carefully monitor the midlatitude conditions as well as tropical cyclones themselves as earlier heavy rainfall warnings may be necessary depending on the former.
Abstract
Radiative feedbacks over interannual time scales can be potentially useful for global warming estimation. However, the diversity of the lead–lag relationships in global mean surface temperature (GMST) and net radiation flux at the top of the atmosphere (GMTOA) create uncertainty during the estimation of radiative feedbacks. In this study, key physical processes controlling lead–lag relationships were elucidated by categorizing preindustrial control simulations of CMIP6 into three groups based on cross correlation values of GMTOA against GMST at lag 0 and lag +1 year. The diversity in the lead–lag relationships was primarily caused by the climatological state difference of the atmosphere over the equatorial Pacific, which modulated the strength of convective activity and sensitivity of low-level clouds. Diminished atmospheric stability caused enhanced convective activity, more efficient energy release, and smaller lags. In addition, enhanced stability in the lower atmosphere rendered the low-level clouds more sensitive to sea surface temperature changes and considerably delayed the radiative response. The climatological state difference of the atmosphere resulted from model-inherent atmospheric conditions. These findings suggest that the diversity of lead–lag relationships of GMST and GMTOA over interannual time scales could represent the characteristics of general atmospheric circulation models and possible solutions of the actual atmosphere, which could also affect long-term feedback features.
Abstract
Radiative feedbacks over interannual time scales can be potentially useful for global warming estimation. However, the diversity of the lead–lag relationships in global mean surface temperature (GMST) and net radiation flux at the top of the atmosphere (GMTOA) create uncertainty during the estimation of radiative feedbacks. In this study, key physical processes controlling lead–lag relationships were elucidated by categorizing preindustrial control simulations of CMIP6 into three groups based on cross correlation values of GMTOA against GMST at lag 0 and lag +1 year. The diversity in the lead–lag relationships was primarily caused by the climatological state difference of the atmosphere over the equatorial Pacific, which modulated the strength of convective activity and sensitivity of low-level clouds. Diminished atmospheric stability caused enhanced convective activity, more efficient energy release, and smaller lags. In addition, enhanced stability in the lower atmosphere rendered the low-level clouds more sensitive to sea surface temperature changes and considerably delayed the radiative response. The climatological state difference of the atmosphere resulted from model-inherent atmospheric conditions. These findings suggest that the diversity of lead–lag relationships of GMST and GMTOA over interannual time scales could represent the characteristics of general atmospheric circulation models and possible solutions of the actual atmosphere, which could also affect long-term feedback features.
Abstract
Two global atmospheric circulation datasets (ERA5 and NCEP-FNL) with horizontal resolutions of 0.25°×0.25° are investigated in terms of kinetic energy (KE) spectra at 200 hPa (roughly between 11 and 12 km). The horizontal KE (HKE) in NCEP-FNL is larger and flatter than that in ERA5 at subsynoptic scales and mesoscales. Restoring the energy of this wavenumber range to the physical space shows that the HKE in NCEP-FNL is larger than that in ERA5 over most areas but smaller mainly in the Indo-Pacific warm pool. The spectral budgets show that at these scales, the positive contribution from net vertical flux in ERA5 is stronger than that in NCEP-FNL, while the negative contribution from available potential energy (APE) conversion is smaller; assuming that the atmosphere is in a quasi-stationary state, more dissipation is found in ERA5 than in NCEP-FNL, which should be responsible for the HKE spectrum in ERA5 to be steeper and weaker than that in NCEP-FNL. Our formulation shows that the APE conversion and net vertical flux are related to the pressure vertical velocity (PVV). The APE conversion and net vertical flux differences between the two datasets, like the PVV difference, are mainly from the tropical region. At large scales, the vertical motion in ERA5 is larger than that in NCEP-FNL. The amplitude differences of the PVV spectra between two datasets are consistent with those of the large-scale precipitation spectra associated with microphysics parameterizations. These results support that vertical motion is a key dynamical factor explaining energy discrepancies at mesoscales.
Abstract
Two global atmospheric circulation datasets (ERA5 and NCEP-FNL) with horizontal resolutions of 0.25°×0.25° are investigated in terms of kinetic energy (KE) spectra at 200 hPa (roughly between 11 and 12 km). The horizontal KE (HKE) in NCEP-FNL is larger and flatter than that in ERA5 at subsynoptic scales and mesoscales. Restoring the energy of this wavenumber range to the physical space shows that the HKE in NCEP-FNL is larger than that in ERA5 over most areas but smaller mainly in the Indo-Pacific warm pool. The spectral budgets show that at these scales, the positive contribution from net vertical flux in ERA5 is stronger than that in NCEP-FNL, while the negative contribution from available potential energy (APE) conversion is smaller; assuming that the atmosphere is in a quasi-stationary state, more dissipation is found in ERA5 than in NCEP-FNL, which should be responsible for the HKE spectrum in ERA5 to be steeper and weaker than that in NCEP-FNL. Our formulation shows that the APE conversion and net vertical flux are related to the pressure vertical velocity (PVV). The APE conversion and net vertical flux differences between the two datasets, like the PVV difference, are mainly from the tropical region. At large scales, the vertical motion in ERA5 is larger than that in NCEP-FNL. The amplitude differences of the PVV spectra between two datasets are consistent with those of the large-scale precipitation spectra associated with microphysics parameterizations. These results support that vertical motion is a key dynamical factor explaining energy discrepancies at mesoscales.
Abstract
Large-amplitude internal solitary waves were recently observed in a coastal plain estuary and were hypothesized to evolve from an internal lee wave generated at the channel-shoal interface. To test this mechanism, a 3D nonhydrostatic model with nested domains and adaptive grids was used to investigate the generation of the internal solitary waves and their subsequent nonlinear evolution. A complex sequence of wave propagation and transformation was documented and interpreted using the nonlinear wave theory based on the Korteweg-de Vries equation. During the ebb tide a mode-2 internal lee wave is generated by the interaction between lateral flows and channel-shoal topography. This mode-2 lee wave subsequently propagates onto the shallow shoal and transforms into a mode-1 wave of elevation as strong mixing on the flood tide erases stratification in the bottom boundary layer and the lower branch of the mode-2 wave. The mode-1 wave of elevation evolves into an internal solitary wave due to nonlinear steepening and spatial changes in the wave phase speed. As the solitary wave of elevation continues to propagate over the shoaling bottom, the leading edge moves ahead as a rarefaction wave while the trailing edge steepens and disintegrates into a train of rank-ordered internal solitary waves, due to the combined effects of shoaling and dispersion. Strong turbulence in the bottom boundary layer dissipates wave energy and causes the eventual destruction of the solitary waves. In the meantime, the internal solitary waves can generate elevated shear and dissipation rate in local regions.
Abstract
Large-amplitude internal solitary waves were recently observed in a coastal plain estuary and were hypothesized to evolve from an internal lee wave generated at the channel-shoal interface. To test this mechanism, a 3D nonhydrostatic model with nested domains and adaptive grids was used to investigate the generation of the internal solitary waves and their subsequent nonlinear evolution. A complex sequence of wave propagation and transformation was documented and interpreted using the nonlinear wave theory based on the Korteweg-de Vries equation. During the ebb tide a mode-2 internal lee wave is generated by the interaction between lateral flows and channel-shoal topography. This mode-2 lee wave subsequently propagates onto the shallow shoal and transforms into a mode-1 wave of elevation as strong mixing on the flood tide erases stratification in the bottom boundary layer and the lower branch of the mode-2 wave. The mode-1 wave of elevation evolves into an internal solitary wave due to nonlinear steepening and spatial changes in the wave phase speed. As the solitary wave of elevation continues to propagate over the shoaling bottom, the leading edge moves ahead as a rarefaction wave while the trailing edge steepens and disintegrates into a train of rank-ordered internal solitary waves, due to the combined effects of shoaling and dispersion. Strong turbulence in the bottom boundary layer dissipates wave energy and causes the eventual destruction of the solitary waves. In the meantime, the internal solitary waves can generate elevated shear and dissipation rate in local regions.
