Center for Ocean-Atmospheric Prediction Studies

Abstract: Often the parameterization of the moisture roughness length is not seen as being important, as long as the parameterization seems reasonable; that is, it is within the rather considerable bounds of error for the data sets used to determine the parameterization. However, the choice of parameterization does influence height adjustments of humidity and calculations of turbulent heat fluxes. This paper focuses on the calculation of the turbulent heat fluxes using different parameterizations of roughness length. Five roughness length parameterizations are examined herein. These parameterizations include wall theory; the Clayson, Fairall, Curry parameterization; the Liu, Katsaros, Businger parameterization; Zilitinkevich et al. parameterization; and the COARE3.0 parameterization. Turbulent heat fluxes are calculated from each parameterization of the roughness length and are compared to observed turbulent heat flux values. The bulk latent heat flux estimates have a much better signal to noise ratio than the sensible heat fluxes, and are therefore the focus of the comparison to observations. This comparison indicates how to improve the proportionality in the above roughness length parameterizations, which are causing modeled turbulent heat flux magnitudes to be too large in four of the five parameterizations. The modeled turbulent heat fluxes are evaluated again after the modification of the parameterizations. Significant improvements in both the bias and the root mean square error (RMSE) are seen. Three parameterizations see roughly the same improvements of around 17Wm^-2 in the bias and roughly 10Wm^-2 in the RMSE. The largest improvements are in the Liu, Katsaros, Businger parameterization with bias improvements of over 45Wm^-2 and a RMSE reduction of nearly 32Wm^-2.

Abstract: Mexico and Central America (Mesoamerica) are situated in a complex and unique geographical position with the Caribbean Sea to the East and the tropical Eastern Pacific Ocean to the West. The weather patterns of this region are driven by winds, temperatures, moisture, and orography of several mountain ranges. This study finds the dates of the onset and demise of rainfall regimes on a specific day using NASA’s Tropical Rainfall Measuring Mission (TRMM) rainfall for years 1998–2012, area-averaged over land. Using NASA’s MERRA-2 Reanalysis data, we also look at the phenomenology of the triggers of the rainy season onset and demise on the daily time-scale instead of the monthly scales used by previous studies.
We find that the Mesoamerican Rainy Season can be distinguished into two parts: the Early Spring Rainfall (ESR) associated with light rains and the Late Spring Rainfall (LSR) associated with heavy rains. Two algorithms are used to obtain these rainy season distinctions. A new algorithm was developed during this study, called the SLOPE algorithm, to calculate when the rain rates first start to increase. In the second method, the daily cumulative anomalies of rainfall are compared to the climatological rainfall to find the time of onset of the heavy rains, called the MINCA algorithm. To better understand the phenomenology associated with the timing of the rainfall, we look at the monsoon trough, moisture flux convergence, moist static energy anomalies, and the weakening/strengthening of the winds associated with the Caribbean Low-Level Jet and Panama Jet.
The light rain rates begin, on average, in mid-March, approximately one month after the peak of the winter Caribbean Low-Level Jet and the Panama Jet. The ramp-up between the light rains and heavy rains is associated with a significant weakening of both jets and the northward progression of a monsoon trough off the western coast of Central America. The heavy rain rates begin, on average, in mid-May, and are associated with the timing when the Panama Jet goes to near zero magnitude and a strong monsoon trough in the eastern Pacific. At the demise of the rainfall, approximately in mid-November, the Panama Jet strengthens again, the total moisture flux convergence decreases significantly, and the monsoon trough retreats southward and eastward. The results of this study have positive implications in agriculture and water resources for Mesoamerica, as this information may help resource managers better plan and adapt to climate variability.

Abstract: Despite the fact that latent heating in cloud systems drives many atmospheric circulations, including tropical cyclones, little is known of its magnitude and structure due in large part to inadequate observations. In this work, a reasonably high-resolution (2 km), four-dimensional airborne Doppler radar retrieval of the latent heat of condensation is presented for rapidly intensifying Hurricane Guillermo (1997). Several advancements in the retrieval algorithm are shown including: (1) analyzing the scheme within the dynamically consistent framework of a numerical model, (2) identifying algorithm sensitivities through the use of ancillary data sources and (3) developing a precipitation budget storage term parameterization. The determination of the saturation state is shown to be an important part of the algorithm for updrafts of ~ 5 m s-1 or less. The uncertainties in the magnitude of the retrieved heating are dominated by errors in the vertical velocity. Using a combination of error propagation and Monte Carlo uncertainty techniques, biases were found to be small, and randomly distributed errors in the heating magnitude were ~16 % for updrafts greater than 5 m s-1 and ~156 % for updrafts of 1 m s- 1. The impact of the retrievals is assessed by inserting the heating into realistic numerical simulations at 2 km resolution and comparing the generated wind structure to the Doppler radar observations of Guillermo. Results show that using the latent heat retrievals outperforms a simulation that relies on a state-of-the-art microphysics scheme (Reisner and Jeffery 2009), in terms of wind speed root-mean-square errors, explained variance and eye/eyewall structure. The incorrect transport of water vapor (a function of the sub-grid model and the numerical approximations to advection) and the restrictions on the magnitude of heat release that ensure the present model's stability are suggested as sources of error in the simulation without the retrievals. Motivated by the latent heat retrievals, the dynamics of vortex axisymmetrization from the perspective of thermal anomalies is investigated using an idealized, non-linear atmospheric model (HIGRAD). Attempts at reproducing the results of previous work (Nolan and Grasso 2003; NG03) revealed a discrepancy with the impacts of purely asymmetric forcing. While NG03 found that purely asymmetric heating led to a negligible, largely negative impact on the vortex intensification, in the present study the impacts of asymmetries are found to have an important, largely positive role. Absolute angular momentum budgets revealed that the essential difference between the present work and that of NG03 was the existence of a significant, axisymmetric secondary circulation in the basic-state vortex used in the HIGRAD simulations. This secondary circulation was larger than that present in NG03's simulations. The spin-up of the vortex caused by the asymmetric thermal anomalies was dominated by the axisymmetric fluxes of angular momentum at all times, indicating fundamentally different evolution of asymmetries in the presence of radial flow. Radial momentum budgets were performed to elucidate the mechanisms responsible for the formation of the physically significant secondary circulation. Results show that explicit (sub-grid) diffusion in the model was producing a gradient wind imbalance, which drives a radial inflow and associated secondary circulation in an attempt to re-gain balance. In addition, the production of vorticity anomalies from the asymmetric heating was found to be sensitive to the eddy diffusivity, with large differences between HIGRAD and the widely used WRF model for the exact same value of this uncertain parameter.