Savage, A. C., Arbic, B. K., Alford, M. H., Ansong, J. K., Farrar, J. T., Menemenlis, D., et al. (2017). Spectral decomposition of internal gravity wave sea surface height in global models: INTERNAL GRAVITY WAVE SEA SURFACE HEIGHT. J. Geophys. Res. Oceans , 122 (10), 7803–7821.
Abstract: Two global ocean models ranging in horizontal resolution from 1/128 to 1/488 are used to study the space and time scales of sea surface height (SSH) signals associated with internal gravity waves (IGWs). Frequency-horizontal wavenumber SSH spectral densities are computed over seven regions of the world ocean from two simulations of the HYbrid Coordinate Ocean Model (HYCOM) and three simulations of the Massachusetts Institute of Technology general circulation model (MITgcm). High wavenumber, high-frequency SSH variance follows the predicted IGW linear dispersion curves. The realism of high-frequency motions (>0:87 cpd) in the models is tested through comparison of the frequency spectral density of dynamic height variance computed from the highest-resolution runs of each model (1/258 HYCOM and 1/488 MITgcm) with dynamic height variance frequency spectral density computed from nine in situ profiling instruments. These high-frequency motions are of particular interest because of their contributions to the small-scale SSH variability that will be observed on a global scale in the upcoming Surface Water and Ocean Topography (SWOT) satellite altimetry mission. The variance at supertidal frequencies can be comparable to the tidal and low-frequency variance for high wavenumbers (length scales smaller than 50 km), especially in the higher-resolution simulations. In the highest-resolution simulations, the high-frequency variance can be greater than the low-frequency variance at these scales.
Luecke, C. A., Arbic, B. K., Bassette, S. L., Richman, J. G., Shriver, J. F., Alford, M. H., et al. (2017). The Global Mesoscale Eddy Available Potential Energy Field in Models and Observations: GLOBAL LOW-FREQUENCY EDDY APE. J. Geophys. Res. Oceans , 122 (11), 9126–9143.
Abstract: 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.
Yu, B., Seed, A., Pu, L., & Malone, T. (2019). Integration of weather radar data into a raster GIS framework for improved flood estimation. Atmos. Sci. Lett. , 6 (1).
Abstract: We present in this paper the interannual variability of seasonal temperature and rainfall in the Indian meteorological subdivisions (IMS) for boreal winter and summer seasons that take in to account the varying length of the seasons.Our study reveals that accounting for the variations in the length of the sea-sons produces stronger teleconnections between the seasonal anomalies of surface temperature and rainfall over India with corresponding sea surface temperature anomalies of the tropical Oceans (especially over the northern Indian and the equatorial Pacific Oceans) compared to the same teleconnections from fixed length seasons over the IMS. It should be noted that the IMS show significant spatial heterogeneity in these teleconnections
Smith, S. R., Briggs, K., Bourassa, M. A., Elya, J., & Paver, C. R. (2018). Shipboard automated meteorological and oceanographic system data archive: 2005-2017. Geosci Data J , 5 (2), 73–86.
Abstract: Since 2005, the Shipboard Automated Meteorological and Oceanographic System (SAMOS) initiative has been collecting, quality-evaluating, distributing, and archiving underway navigational, meteorological, and oceanographic observations from research vessels. Herein we describe the procedures for acquiring ship and instrumental metadata and the one-minute interval observations from 44 research vessels that have contributed to the SAMOS initiative from 2005 to 2017. The overall data processing workflow and quality control procedures are documented along with data file formats and version control procedures. The SAMOS data are disseminated to the user community via web, FTP, and Thematic Real-time Environmental Distributed Data Services from both the Marine Data Center at the Florida State University and the National Centers for Environmental Information, which serves as the long-term archive for the SAMOS initiative. They have been used to address topics ranging from air-sea interaction studies, the calibration, evaluation, and development of satellite observational products, the evaluation of numerical atmospheric and ocean models, and the development of new tools and techniques for geospatial data analysis in the informatics community. Maps provide users the geospatial coverage within the SAMOS dataset, with a focus on the Essential Climate/Ocean Variables, and recommendations are made regarding which versions of the dataset should be accessed by different user communities.
Hoffman, R. N., Privé, N., & Bourassa, M. (2017). Comments on “Reanalyses and Observations: What's the Difference?”. Bull. Amer. Meteor. Soc. , 98 (11), 2455–2459.
Abstract: Are there important differences between reanalysis data and familiar observations and measurements? If so, what are they? This essay evaluates four possible answers that relate to: the role of inference, reliance on forecasts, the need to solve an ill-posed inverse problem, and understanding of errors and uncertainties. The last of these is argued to be most significant. The importance of characterizing uncertainties associated with results—whether those results are observations or measurements, analyses or reanalyses, or forecasts—is emphasized.
Maloney, E. D., Gettelman, A., Ming, Y., Neelin, J. D., Barrie, D., Mariotti, A., et al. (2019). Process-Oriented Evaluation of Climate and Weather Forecasting Models. Bull. Amer. Meteor. Soc. , 100 (9), 1665–1686.
Abstract: Realistic climate and weather prediction models are necessary to produce confidence in projections of future climate over many decades and predictions for days to seasons. These models must be physically justified and validated for multiple weather and climate processes. A key opportunity to accelerate model improvement is greater incorporation of process-oriented diagnostics (PODs) into standard packages that can be applied during the model development process, allowing the application of diagnostics to be repeatable across multiple model versions and used as a benchmark for model improvement. A POD characterizes a specific physical process or emergent behavior that is related to the ability to simulate an observed phenomenon. This paper describes the outcomes of activities by the Model Diagnostics Task Force (MDTF) under the NOAA Climate Program Office (CPO) Modeling, Analysis, Predictions and Projections (MAPP) program to promote development of PODs and their application to climate and weather prediction models. MDTF and modeling center perspectives on the need for expanded process-oriented diagnosis of models are presented. Multiple PODs developed by the MDTF are summarized, and an open-source software framework developed by the MDTF to aid application of PODs to centers' model development is presented in the context of other relevant community activities. The paper closes by discussing paths forward for the MDTF effort and for community process-oriented diagnosis.
Liu, Y., Tan, Z. - M., & Wu, Z. (2019). Noninstantaneous Wave-CISK for the Interaction between Convective Heating and Low-Level Moisture Convergence in the Tropics. J. Atmos. Sci. , 76 (7), 2083–2101.
Abstract: The interaction between tropical convective heating and thermally forced circulation is investigated using a global dry primitive-equation model with the parameterization of wave-conditional instability of the second kind (CISK). It is demonstrated that deep convective heating can hardly sustain itself through the moisture convergence at low levels regardless of the fraction of immediate consumption of converged moisture. In contrast, when the fraction is large, shallow convective heating and its forced circulation exhibit preferred growth of small scales. As the “CISK catastrophe” mainly comes from the instantaneous characters of moisture-convection feedback in the conventional wave-CISK, a noninstantaneous wave-CISK is proposed, which highlights the accumulation-consumption (AC) time scale for the convective heating accumulation and/or the converged moisture consumption. In the new wave-CISK, once moisture is converged, the release of latent heat takes place gradually within an AC time scale. In this sense, convective heating is not only related to the instantaneous moisture convergence at the current time, but also to that which occurred in the past period of the AC time scale. The noninstantaneous wave-CISK could guarantee the occurrence of convective heating and/or moisture convergence at larger scales, and then favor the growth of long waves, and thus solve the problem of CISK catastrophe. With the new wave-CISK and AC time scale of 2 days, the simulated convective heating-driven system bears a large similarity to that of the observed convectively coupled Kelvin wave.
Holbach, H. M., Uhlhorn, E. W., & Bourassa, M. A. (2018). Off-Nadir SFMR Brightness Temperature Measurements in High-Wind Conditions. J. Atmos. Oceanic Technol. , 35 (9), 1865–1879.
Abstract: Wind and wave-breaking directions are investigated as potential sources of an asymmetry identified in off-nadir remotely sensed measurements of ocean surface brightness temperatures obtained by the Stepped Frequency Microwave Radiometer (SFMR) in high-wind conditions, including in tropical cyclones. Surface wind speed, which dynamically couples the atmosphere and ocean, can be inferred from SFMR ocean surface brightness temperature measurements using a radiative transfer model and an inversion algorithm. The accuracy of the ocean surface brightness temperature to wind speed calibration relies on accurate knowledge of the surface variables that are influencing the ocean surface brightness temperature. Previous studies have identified wind direction signals in horizontally polarized radiometer measurements in low to moderate (0�20 m s−1) wind conditions over a wide range of incidence angles. This study finds that the azimuthal asymmetry in the off-nadir SFMR brightness temperature measurements is also likely a function of wind direction and extends the results of these previous studies to high-wind conditions. The off-nadir measurements from the SFMR provide critical data for improving the understanding of the relationships between brightness temperature, surface wave�breaking direction, and surface wind vectors at various incidence angles, which is extremely useful for the development of geophysical model functions for instruments like the Hurricane Imaging Radiometer (HIRAD).
Misra, V., & Bhardwaj, A. (2019). Defining the Northeast Monsoon of India. Mon. Wea. Rev. , 147 (3), 791–807.
Abstract: This study introduces an objective definition for onset and demise of the Northeast Indian Monsoon (NEM). The definition is based on the land surface temperature analysis over the Indian subcontinent. It is diagnosed from the inflection points in the daily anomaly cumulative curve of the area-averaged surface temperature over the provinces of Andhra Pradesh, Rayalseema, and Tamil Nadu located in the southeastern part of India. Per this definition, the climatological onset and demise dates of the NEM season are 6 November and 13 March, respectively. The composite evolution of the seasonal cycle of 850hPa winds, surface wind stress, surface ocean currents, and upper ocean heat content suggest a seasonal shift around the time of the diagnosed onset and demise dates of the NEM season. The interannual variations indicate onset date variations have a larger impact than demise date variations on the seasonal length, seasonal anomalies of rainfall, and surface temperature of the NEM. Furthermore, it is shown that warm El Niño�Southern Oscillation (ENSO) episodes are associated with excess seasonal rainfall, warm seasonal land surface temperature anomalies, and reduced lengths of the NEM season. Likewise, cold ENSO episodes are likely to be related to seasonal deficit rainfall anomalies, cold land surface temperature anomalies, and increased lengths of the NEM season.
Ahern, K., Bourassa, M. A., Hart, R. E., Zhang, J. A., & Rogers, R. F. (2019). Observed Kinematic and Thermodynamic Structure in the Hurricane Boundary Layer During Intensity Change. Mon. Wea. Rev. , .
Abstract: The axisymmetric structure of the inner-core hurricane boundary layer (BL) during intensification [IN; intensity tendency ≥ 20 kt (24 h)−1], weakening [WE; intensity tendency < −10 kt (24 h)−1], and steady-state [SS; the remainder] periods are analyzed using composites of GPS dropwindsondes from reconnaissance missions between 1998 and 2015. A total of 3,091 dropsondes were composited for analysis below 2.5 km elevation—1,086 during IN, 1,042 during WE, and 963 during SS. In non-intensifying hurricanes, the lowlevel tangential wind is greater outside the radius of maximum wind (RMW) than for intensifying hurricanes, implying higher inertial stability (I) at those radii for non-intensifying hurricanes. Differences in tangential wind structure (and I) between the groups also imply differences in secondary circulation. The IN radial inflow layer is of nearly equal or greater thickness than nonintensifying groups, and all groups show an inflow maximum just outside the RMW. Non-intensifying hurricanes have stronger inflow outside the eyewall region, likely associated with frictionally forced ascent out of the BL and enhanced subsidence into the BL at radii outside the RMW. Equivalent potential temperatures (θe) and conditional stability are highest inside the RMW of non-intensifying storms, which is potentially related to TC intensity. At greater radii, inflow layer θe is lowest in WE hurricanes, suggesting greater subsidence or more convective downdrafts at those radii compared to IN and SS hurricanes. Comparisons of prior observational and theoretical studies are highlighted, especially those relating BL structure to large-scale vortex structure, convection, and intensity.