In this paper, we evaluate several timely, daily air-sea heat flux products (NCEP, NCEP2, ERA-Interim and OAFlux/ISCCP) against observations and present the newly developed TropFlux product. This new product uses bias-corrected ERA-interim and ISCCP data as input parameters to compute air-sea fluxes from the COARE v3.0 algorithm. Wind speed is corrected for mesoscale gustiness. Surface net shortwave radiation is based on corrected ISCCP data. We extend the shortwave radiation time series by using ''near real-time'' SWR estimated from outgoing longwave radiation. All products reproduce consistent intraseasonal surface net heat flux variations associated with the Madden-Julian Oscillation in the Indian Ocean, but display more disparate interannual heat flux variations associated with El Niño in the eastern Pacific. They also exhibit marked differences in mean values and seasonal cycle. Comparison with global tropical moored buoy array data, I-COADS and fully independent mooring data sets shows that the two NCEP products display lowest correlation to mooring turbulent fluxes and significant biases. ERA-interim data captures well temporal variability, but with significant biases. OAFlux and TropFlux perform best. All products have issues in reproducing observed longwave radiation. Shortwave flux is much better captured by ISCCP data than by any of the re-analyses. Our ''near real-time'' shortwave radiation performs better than most re-analyses, but tends to underestimate variability over the cold tongues of the Atlantic and Pacific. Compared to independent mooring data, NCEP and NCEP2 net heat fluxes display *0.78 correlation and [65 W m -2 rms-difference, ERA-I performs better (*0.86 correlation and *48 W m -2 ) while OAFlux and TropFlux perform best (*0.9 correlation and *43 W m -2 ). TropFlux hence provides a useful option for studying flux variability associated with ocean-atmosphere interactions, oceanic heat budgets and climate fluctuations in the tropics.

Turbulent mixing and entrainment at the boundary of a cloud is studied by means of direct numerical simulations that couple the Eulerian description of the turbulent velocity and water vapor fields with a Lagrangian ensemble of cloud water droplets that can grow and shrink by condensation and evaporation, respectively. The focus is on detailed analysis of the relaxation process of the droplet ensemble during the entrainment of subsaturated air, in particular the dependence on turbulence time scales, droplet number density, initial droplet radius and particle inertia. We find that the droplet evolution during the entrainment process is captured best by a phase relaxation time that is based on the droplet number density with respect to the entire simulation domain and the initial droplet radius. Even under conditions favoring homogeneous mixing, the probability density function of supersaturation at droplet locations exhibits initially strong negative skewness, consistent with droplets near the cloud boundary being suddenly mixed into clear air, but rapidly approaches a narrower, symmetric shape. The droplet size distribution, which is initialized as perfectly monodisperse, broadens and also becomes somewhat negatively skewed. Particle inertia and gravitational settling lead to a more rapid initial evaporation, but ultimately only to slight depletion of both tails of the droplet size distribution. The Reynolds number dependence of the mixing process remained weak over the parameter range studied, most probably due to the fact that the inhomogeneous mixing regime could not be fully accessed when phase relaxation times based on global number density are considered.Comment: 17 pages, 10 Postscript figures (figures 3,4,6,7,8 and 10 are in reduced quality), to appear in Theoretical Computational Fluid Dynamic

The entrainment of clear air and its subsequent mixing with a filament of cloudy air, as occurs at the edge of a cloud, is studied in three-dimensional direct numerical simulations that combine the Eulerian description of the turbulent velocity, temperature, and vapor fields with a Lagrangian cloud droplet ensemble. Forced and decaying turbulence is considered, such as when the dynamics around the filament is driven by larger-scale eddies or during the final period of the life cycle of a cloud. The microphysical response depicted in n d 2 hr 3 i space (where n d and r are droplet number density and radius, respectively) shows characteristics of both homogeneous and inhomogeneous mixing, depending on the Damk€ ohler number. The transition from inhomogeneous to homogeneous mixing leads to an offset of the homogeneous mixing curve to larger dilution fractions. The response of the system is governed by the smaller of the single droplet evaporation time scale and the bulk phase relaxation time scale. Variability within the n d 2 hr 3 i space increases with decreasing sample volume, especially during the mixing transients. All of these factors have implications for the interpretation of measurements in clouds. The qualitative mixing behavior changes for forced versus decaying turbulence, with the latter yielding remnant patches of unmixed cloud and stronger fluctuations. Buoyancy due to droplet evaporation is observed to play a minor role in the mixing for the present configuration. Finally, the mixing process leads to the transient formation of a pronounced nearly exponential tail of the probability density function of the Lagrangian supersaturation, and a similar tail emerges in the droplet size distribution under inhomogeneous conditions.

During boreal winter, there is a prominent maximum of intraseasonal sea-surface temperature (SST) variability associated with the Madden-Julian Oscillation (MJO) along a Thermocline Ridge located in the southwestern Indian Ocean (5 o S−10°S, 60°E-90°E; TRIO region). There is an ongoing debate about the relative importance of air-sea heat fluxes and oceanic processes in driving this intraseasonal SST variability. Furthermore, various studies have suggested that interannual variability of the oceanic structure in the TRIO region could modulate the amplitude of the MJOdriven SST response. In this study, we use observations and ocean general circulation model (OGCM) experiments to quantify these two effects over the 1997−2006 period. Observational analysis indicates that Ekman pumping does not contribute significantly (on average) to intraseasonal SST variability. It is, however, difficult to quantify the relative contribution of net heat fluxes and entrainment to SST intraseasonal variability from observations alone. We therefore use a suite of OGCM experiments to isolate the impacts of each process. During 1997−2006, wind stress contributed only 20% of the intraseasonal SST variability (averaged over the TRIO region), while heat fluxes contributed to 70%, with forcing by shortwave radiation (75%) dominating the other flux components (25%). This estimate is consistent with an independent air-sea flux product, which indicates that shortwave radiation contributes 68% of intraseasonal heat flux variability. The time scale of the heat-flux perturbation, in addition to its amplitude, is also important in controlling the intraseasonal SST signature, with longer periods favouring a larger response. There are also strong year-to-year variations in the respective role of heat fluxes and wind stress. Intraseasonal-wind stress dominates the SST signature for one (in 2001) and contributes significantly to another (in 2000) of the 5 strong cooling events identified in both observations and the model (2 in 1999, 1 in 2000, 2001 and 2002). Interannual variations of the subsurface thermal structure associated with the Indian Ocean Dipole or El Niño/La Niña events modulate the MJOdriven SST signature only moderately (by up to 30%), mainly by changing the temperature of water entrained into the mixed layer. Our results therefore suggest that the primary factor that controls year-to-year changes in the amplitude of TRIO, intraseasonal SST anomalies is the amplitude and timescale of the intraseasonal heat-flux perturbations.