Publications

Export 19 results:
Sort by: Author [ Title  (Asc)] Type Year
A B C D E F G H I J K L M N O P Q R [S] T U V W X Y Z   [Show ALL]
A
Valero, FPJ, Bucholtz A, Bush BC, Pope SK, Collins WD, Flatau P, Strawa A, Gore WJY.  1997.  Atmospheric Radiation Measurements Enhanced Shortwave Experiment (ARESE): Experimental and data details. Journal of Geophysical Research-Atmospheres. 102:29929-29937.   10.1029/97jd02434   AbstractWebsite

Atmospheric Radiation Measurements Enhanced Shortwave Experiment (ARESE) was conducted to study the magnitude and spectral characteristics of the absorption Of solar radiation by the clear and cloudy atmosphere. Three aircraft platforms, a Grob Egrett, a NASA ER-2, and a Twin Otter, were used during ARESE in conjunction with the Atmospheric Radiation Measurements (ARM) central and extended facilities in north central Oklahoma. The aircraft were coordinated to simultaneously measure solar irradiances in the total spectral broadband (0.224-3.91 mu m), near infrared broadband (0.678-3.3 mu m), and in seven narrow band-pass (similar to 10 nm width) channels centered at 0.500, 0.862, 1.064, 1.249, 1.501, 1.651, and 1.750 mu m. Instrumental calibration issues are discussed in some detail, in particular radiometric power, angular, and spectral responses. The data discussed in this paper are available at the ARM ARESE data archive via anonymous FTP to ftp.arm.gov.

C
Markowicz, KM, Flatau PJ, Kardas AE, Remiszewska J, Stelmaszczyk K, Woeste L.  2008.  Ceilometer retrieval of the boundary layer vertical aerosol extinction structure. Journal of Atmospheric and Oceanic Technology. 25:928-944.   10.1175/2007jtecha1016.1   AbstractWebsite

The CT25K ceilometer is a general-purpose cloud height sensor employing lidar technology for detection of clouds. In this paper it is shown that it can also be used to retrieve aerosol optical properties in the boundary layer. The authors present a comparison of the CT25K instrument with the aerosol lidar system and discuss its good overall agreement for both the range-corrected signals and the retrieved extinction coefficient profiles. The CT25K aerosol profiling is mostly limited to the boundary layer, but it is capable of detecting events in the lower atmosphere such as mineral dust events between 1 and 3 km. Assumptions needed for the estimation of the aerosol extinction profiles are discussed. It is shown that, when a significant part of the aerosol layer is in the boundary layer, knowledge of the aerosol optical depth from a sun photometer allows inversion of the lidar signal. In other cases, surface observations of the aerosol optical properties are used. It is demonstrated that additional information from a nephelometer and aethalometer allows definition of the lidar ratio. Extinction retrievals based on spherical and randomly oriented spheroid assumptions are performed. It is shown, by comparison with the field measurements during the United Arab Emirates Unified Aerosol Experiment, that an assumption about specific particle shape is important for the extinction profile inversions. The authors indicate that this limitation of detection is a result of the relatively small sensitivity of this instrument in comparison to more sophisticated aerosol lidars. However, in many cases this does not play a significant role because globally only about 20% of the aerosol optical depth is above the boundary layer.

Schmidt, JM, Flatau PJ, Yates RD.  2014.  Convective cells in altocumulus observed with a high-resolution radar. Journal of the Atmospheric Sciences. 71:2130-2154.   10.1175/jas-d-13-0172.1   AbstractWebsite

Very-high-resolution Doppler radar observations are used together with aircraft measurements to document the dynamic and thermodynamic structure of a dissipating altocumulus cloud system associated with a deep virga layer. The cloud layer circulation is shown to consist of shallow vertical velocity couplets near cloud top and a series of subkilometer-scale Rayleigh-Benard-like cells that extend vertically through the depth of the cloud layer. The subcloud layer was observed to contain a number of narrow virga fall streaks that developed below the more dominant Rayleigh-Benard updraft circulations in the cloud layer. These features were discovered to be associated with kilometer-scale horizontally orientated rotor circulations that formed along the lateral flanks of the streaks collocated downdraft circulation. The Doppler analysis further reveals that a layer mean descent was present throughout both the cloud and subcloud layers. This characteristic of the circulation is analyzed with regard to the diabatic and radiative forcing on horizontal length scales ranging from the Rayleigh-Benard circulations to the overall cloud layer width. In particular, linear analytical results indicate that a deep and broad mesoscale region of subsidence is quickly established in middle-level cloud layers of finite width when a layer-wide horizontal gradient in the cloud-top radiative cooling rate is present. A conceptual model summarizing the primary observed and inferred circulation features of the altocumulus layer is presented.

D
Flatau, MK, Flatau PJ, Schmidt J, Kiladis GN.  2003.  Delayed onset of the 2002 Indian monsoon. Geophysical Research Letters. 30   10.1029/2003gl017434   AbstractWebsite

[1] We show that there is a set of dynamical predictors, which facilitate forecasting of a delayed monsoon onset. The main dynamical contributor is the early May propagation of the "bogus onset Intraseasonal Oscillation'' which triggers a set of events precluding the climatological monsoon onset. We analyze in detail the 2002 monsoon onset and show that it followed a pattern described in our previous study. We notice that the 2003 monsoon onset followed very similar pattern and was delayed.

E
Schmidt, JM, Flatau PJ, Harasti PR.  2017.  Evidence for a nimbostratus uncinus in a convectively generated mixed-phase stratiform cloud shield. Journal of the Atmospheric Sciences. 74:4093-4116.   10.1175/jas-d-17-0074.1   AbstractWebsite

The structure of a melting layer associated with a mesoconvective system is examined using a combination of in situ aircraft measurements and a unique Doppler radar operated by the U.S. Navy that has a range resolution as fine as 0.5 m. Interest in this case was motivated by ground-based all-sky camera images that captured the transient development of midlevel billow cloud structures within a precipitating trailing stratiform cloud shield associated with a passing deep convective system. A sequence of high-fidelity time-height radar measurements taken of this storm system reveal that the movement of the billow cloud structure over the radar site corresponded with abrupt transitions in the observed low-level precipitation structure. Of particular note is an observed transition from stratiform to more periodic and vertically slanted rain shaft structures that both radar and aircraft measurements indicate have the same temporal periodicity determined to arise visually between successive billow cloud bands. Doppler, balloon, and aircraft measurements reveal these transient bands are associated with a shallow circulation field that resides just above the melting level in a layer of moist neutral stability and strong negative vertical wind shear. The nature of these circulations and their impact on the evolving precipitation field are described in the context of known nimbostratus cloud types.

G
Lelieveld, J, Berresheim H, Borrmann S, Crutzen PJ, Dentener FJ, Fischer H, Feichter J, Flatau PJ, Heland J, Holzinger R, Korrmann R, Lawrence MG, Levin Z, Markowicz KM, Mihalopoulos N, Minikin A, Ramanathan V, de Reus M, Roelofs GJ, Scheeren HA, Sciare J, Schlager H, Schultz M, Siegmund P, Steil B, Stephanou EG, Stier P, Traub M, Warneke C, Williams J, Ziereis H.  2002.  Global air pollution crossroads over the Mediterranean. Science. 298:794-799.   10.1126/science.1075457   AbstractWebsite

The Mediterranean Intensive Oxidant Study, performed in the summer of 2001, uncovered air pollution layers from the surface to an altitude of 15 kilometers. In the boundary layer, air pollution standards are exceeded throughout the region, caused by West and East European pollution from the north. Aerosol particles also reduce solar radiation penetration to the surface, which can suppress precipitation. In the middle troposphere, Asian and to a lesser extent North American pollution is transported from the west. Additional Asian pollution from the east, transported from the monsoon in the upper troposphere, crosses the Mediterranean tropopause, which pollutes the lower stratosphere at middle latitudes.

L
Flatau, PJ, Stephens GL, Draine BT.  1990.  Light-Scattering by Rectangular Solids in the Discrete-Dipole Approximation - a New Algorithm Exploiting the Block-Toeplitz Structure. Journal of the Optical Society of America a-Optics Image Science and Vision. 7:593-600.   10.1364/josaa.7.000593   AbstractWebsite

The discrete-dipole approximation is used to study the problem of light scattering by homogeneous rectangular particles. The structure of the discrete-dipole approximation is investigated, and the matrix formed by this approximation is identified to be a symmetric, block-Toeplitz matrix. Special properties of block-Toeplitz arrays are explored, and an efficient algorithm to solve the dipole scattering problem is provided. Timings for conjugate gradient, Linpack, and block-Toeplitz solvers are given; the results indicate the advantages of the block-Toeplitz algorithm. A practical test of the algorithm was performed on a system of 1400 dipoles, which corresponds to direct inversion of an 8400 × 8400 real matrix. A short discussion of the limitations of the discrete-dipole approximation is provided, and some results for cubes and parallelepipeds are given. We briefly consider how the algorithm may be improved further.

M
Welton, EJ, Voss KJ, Quinn PK, Flatau PJ, Markowicz K, Campbell JR, Spinhirne JD, Gordon HR, Johnson JE.  2002.  Measurements of aerosol vertical profiles and optical properties during INDOEX 1999 using micropulse lidars. Journal of Geophysical Research-Atmospheres. 107   10.1029/2000jd000038   AbstractWebsite

[1] Micropulse lidar (MPL) systems were used to measure aerosol properties during the Indian Ocean Experiment (INDOEX) 1999 field phase. Measurements were made from two platforms: the NOAA ship R/V Ronald H. Brown, and the Kaashidhoo Climate Observatory (KCO) in the Maldives. Sun photometers were used to provide aerosol optical depths (AOD) needed to calibrate the MPL. This study focuses on the height distribution and optical properties (at 523 nm) of aerosols observed during the campaign. The height of the highest aerosols (top height) was calculated and found to be below 4 km for most of the cruise. The marine boundary layer (MBL) top was calculated and found to be less than 1 km. MPL results were combined with air mass trajectories, radiosonde profiles of temperature and humidity, and aerosol concentration and optical measurements. Humidity varied from approximately 80% near the surface to 50% near the top height during the entire cruise. The average value and standard deviation of aerosol optical parameters were determined for characteristic air mass regimes. Marine aerosols in the absence of any continental influence were found to have an AOD of 0.05+/-0.03, an extinction-to-backscatter ratio (S ratio) of 33+/-6 sr, and peak extinction values around 0.05 km(-1) (near the MBL top). The marine results are shown to be in agreement with previously measured and expected values. Polluted marine areas over the Indian Ocean, influenced by continental aerosols, had AOD values in excess of 0.2, S ratios well above 40 sr, and peak extinction values approximately 0.20 km(-1) (near the MBL top). The polluted marine results are shown to be similar to previously published values for continental aerosols. Comparisons between MPL derived extinction near the ship (75 m) and extinction calculated at ship level using scattering measured by a nephelometer and absorption using a particle soot absorption photometer were conducted. The comparisons indicated that the MPL algorithm (using a constant S ratio throughout the lower troposphere) calculates extinction near the surface in agreement with the ship-level measurements only when the MBL aerosols are well mixed with aerosols above. Finally, a review of the MPL extinction profiles showed that the model of aerosol vertical extinction developed during an earlier INDOEX field campaign (at the Maldives) did not correctly describe the true vertical distribution over the greater Indian Ocean region. Using the average extinction profile and AOD obtained during marine conditions, a new model of aerosol vertical extinction was determined for marine atmospheres over the Indian Ocean. A new model of aerosol vertical extinction for polluted marine atmospheres was also developed using the average extinction profile and AOD obtained during marine conditions influenced by continental aerosols.

Conant, WC, Seinfeld JH, Wang J, Carmichael GR, Tang YH, Uno I, Flatau PJ, Markowicz KM, Quinn PK.  2003.  A model for the radiative forcing during ACE-Asia derived from CIRPAS Twin Otter and R/V Ronald H. Brown data and comparison with observations. Journal of Geophysical Research-Atmospheres. 108   10.1029/2002jd003260   AbstractWebsite

Vertical profiles of aerosol size, composition, and hygroscopic behavior from Center for Interdisciplinary Remotely Piloted Aircraft Studies (CIRPAS) Twin Otter and National Oceanic and Atmospheric Administration R/V Ronald H. Brown observations are used to construct a generic optical model of the Asian Pacific Regional Aerosol Characterization Experiment (ACE-Asia) aerosol. The model accounts for sulfate, black carbon, organic carbon, sea salt, and mineral dust. The effects of relative humidity and mixing assumptions (internal versus external, coating of dust by pollutants) are explicitly accounted for. The aerosol model is integrated with a Monte Carlo radiative transfer model to compute direct radiative forcing in the solar spectrum. The predicted regional average surface aerosol forcing efficiency (change in clear-sky radiative flux per unit aerosol optical depth at 500 nm) during the ACE-Asia intensive period is -65 Wm(-2) for pure dust and -60 Wm(-2) for pure pollution aerosol (clear skies). A three-dimensional atmospheric chemical transport model (Chemical Weather Forecast System (CFORS)) is used with the radiative transfer model to derive regional radiative forcing during ACE-Asia in clear and cloudy skies. Net regional solar direct radiative forcing during the 5-15 April 2001 dust storm period is -3 Wm(-2) at the top of the atmosphere and -17 Wm(-2) at the surface for the region from 20degreesN to 50degreesN and 100degreesE to 150degreesE when the effects of clouds on the direct forcing are included. The model fluxes and forcing efficiencies are found to be in good agreement with surface radiometric observations made aboard the R. H. Brown. Mean cloud conditions are found to moderate the top of atmosphere (TOA) radiative forcing by a factor of similar to3 compared to clear-sky calculations, but atmospheric absorption by aerosol is not strongly affected by clouds in this study. The regional aerosol effect at the TOA ("climate forcing") of -3 Wm(-2) is comparable in magnitude, but of opposite sign, to present-day anthropogenic greenhouse gas forcing. The forcing observed during ACE-Asia is similar in character to that seen during other major field experiments downwind of industrial and biomass black carbon sources (e.g., the Indian Ocean Experiment (INDOEX)), insofar as the primary effect of aerosol is to redistribute solar heating from the surface to the atmosphere.

Baranowski, DB, Flatau MK, Flatau PJ, Schmidt JM.  2017.  Multiple and spin off initiation of atmospheric convectively coupled Kelvin waves. Climate Dynamics. 49:2991-3009.   10.1007/s00382-016-3487-7   AbstractWebsite

A novel atmospheric convectively coupled Kelvin wave trajectories database, derived from Tropical Rainfall Measuring Mission precipitation data, is used to investigate initiation of sequential Kelvin wave events. Based on the analysis of beginnings of trajectories from years 1998-2012 it is shown that sequential event initiations can be divided into two distinct categories: multiple initiations and spin off initiations, both of which involve interactions with ocean surface and upper ocean temperature variability. The results of composite analysis of the 83 multiple Kelvin wave initiations show that the local thermodynamic forcing related to the diurnal sea surface temperature variability is responsible for sequential Kelvin wave development. The composite analysis of 91 spin off Kelvin wave initiations shows that the dynamic forcing is a dominant effect and the local thermodynamic forcing is secondary. Detail case studies of both multiple and spin off initiations confirm statistical analysis. A multiple initiation occurs in the presence of the high upper ocean diurnal cycle and a spin off initiation results from both dynamic and local thermodynamic processes. The dynamic forcing is related to increased wind speed and latent heat flux likely associated with an off equatorial circulation. In addition a theoretical study of the sequential Kelvin waves is performed using a shallow water model. Finally, conceptual models of these two types of initiations are proposed.

O
Vogelmann, AM, Flatau PJ, Szczodrak M, Markowicz KM, Minnett PJ.  2003.  Observations of large aerosol infrared forcing at the surface. Geophysical Research Letters. 30   10.1029/2002gl016829   AbstractWebsite

Studies of aerosol effects on the Earth's energy budget usually consider only the cooling effects at short (solar) wavelengths, but we demonstrate that they also have important warming effects at thermal infrared (IR) wavelengths that have rarely been observed and are commonly ignored in climate models. We use high-resolution spectra to obtain the IR radiative forcing at the surface for aerosols encountered in the outflow from northeastern Asia. The spectra were measured by the Marine-Atmospheric Emitted Radiance Interferometer (M-AERI) from the NOAA Ship Ronald H. Brown during the Aerosol Characterization Experiment-Asia (ACE-Asia). We show that the daytime surface IR forcing are often a few Wm(-2) and can reach almost 10 Wm(-2) for large aerosol loadings. Thus, even the smaller aerosol IR forcing observed here are comparable to or greater than the 1 to 2 Wm(-2) IR surface enhancement from increases in greenhouse gases. These results highlight the importance of aerosol IR forcing which should be included in climate model simulations.

Flatau, PJ, Stephens GL.  1988.  On the Fundamental Solution of the Radiative-Transfer Equation. Journal of Geophysical Research-Atmospheres. 93:11037-11050.   10.1029/JD093iD09p11037   AbstractWebsite

This paper outlines the general solution of the one-dimensional, azimuthally averaged radiative transfer equation in terms of a matrix exponential. The link between this fundamental solution and those more commonly used in radiative transfer is established. The formulation is developed for a general vertically inhomogeneous atmosphere with sources. Several new concepts, based on properties of the matrix exponentials, are described in the context of radiative transfer, including the use of the commutator and product integrals. It is also demonstrated how the matrix exponential formulation provides for new insights, not only into improvements of the numerical efficiency and stability of the solution, but also into the understanding of radiative transfer through a layered atmosphere. The various concepts introduced in this paper are illustrated throughout by the two-stream simplification of the general radiative transfer equation.

Maslowska, A, Flatau PJ, Stephens GL.  1994.  On the Validity of the Anomalous Diffraction Theory to Light-Scattering by Cubes. Optics Communications. 107:35-40.   10.1016/0030-4018(94)90099-x   AbstractWebsite

The extinction and absorption efficiencies of a cube at light incidence normal to its four-fold symmetry axis are calculated using the anomalous diffraction theory (ADT). The results are compared with those based on the discrete dipole approximation (DDA). It is shown that for certain cases of the orientation of a cube relative to the direction of the incident light the extinction efficiency calculated using DDA and ADT do not agree. However, the ADT-based absorption efficiencies for the cases studied are dependent on a particle volume and exhibit smaller errors. Hence the validity of the ADT for cubes is not as good as for spheres.

Stramski, D, Wozniak SB, Flatau PJ.  2004.  Optical properties of Asian mineral dust suspended in seawater. Limnology and Oceanography. 49:749-755. AbstractWebsite

The spectral optical properties of Asian mineral dust suspended in seawater exhibit significant variability associated with the origin (and, hence, the chemistry and mineralogy) and particle size distribution of the samples. The measurements of dust samples from different locations show that the mass-specific absorption coefficient of particles, a(p)*, at a wavelength of light lambda = 440 nm, varies from about 0.028 m(2) g(-1) for the soil dust from Chinese desert Pnear Dunhuang to 0.15 m(2) g(-1) for the soil dust of volcanic origin in Cheju Island (South Korea). At lambda = 400 nm, this range is 0.05-0.23 m(2) g(-1). The aerosol sample collected in the Sea of Japan during a massive dust storm in East Asia shows a(p)*(lambda) > 0.1 m(2) g(-1) for lambda < 425 nm. The mass-specific scattering coefficient, b(p)*(lambda), ranges from about 0.8 to 1.5 m(2) g(-1) at blue and green wavelengths for the samples examined. The single scattering albedo, omega(0) increases with wavelength. For lambda > 400 nm, omega(0) was > 0.78 for the sample from Cheju Island and > 0.9 for other samples. In the near-infrared region (750-850 nm), where absorption by dust particles is small or undetectable, omega(0), was close to 1.

R
Schmidt, JM, Flatau PJ, Harasti PR, Yates RD, Littleton R, Pritchard MS, Fischer JM, Fischer EJ, Kohri WJ, Vetter JR, Richman S, Baranowski DB, Anderson MJ, Fletcher E, Lando DW.  2012.  Radar observations of individual rain drops in the free atmosphere. Proceedings of the National Academy of Sciences of the United States of America. 109:9293-9298.   10.1073/pnas.1117776109   AbstractWebsite

Atmospheric remote sensing has played a pivotal role in the increasingly sophisticated representation of clouds in the numerical models used to assess global and regional climate change. This has been accomplished because the underlying bulk cloud properties can be derived from a statistical analysis of the returned microwave signals scattered by a diverse ensemble comprised of numerous cloud hydrometeors. A new Doppler radar, previously used to track small debris particles shed from the NASA space shuttle during launch, is shown to also have the capacity to detect individual cloud hydrometeors in the free atmosphere. Similar to the traces left behind on film by subatomic particles, larger cloud particles were observed to leave a well-defined radar signature (or streak), which could be analyzed to infer the underlying particle properties. We examine the unique radar and environmental conditions leading to the formation of the radar streaks and develop a theoretical framework which reveals the regulating role of the background radar reflectivity on their observed characteristics. This main expectation from theory is examined through an analysis of the drop properties inferred from radar and in situ aircraft measurements obtained in two contrasting regions of an observed multicellular storm system. The observations are placed in context of the parent storm circulation through the use of the radar's unique high-resolution waveforms, which allow the bulk and individual hydrometeor properties to be inferred at the same time.

Stephens, GL, Tsay SC, Stackhouse PW, Flatau PJ.  1990.  The Relevance of the Microphysical and Radiative Properties of Cirrus Clouds to Climate and Climatic Feedback. Journal of the Atmospheric Sciences. 47:1742-1753.   10.1175/1520-0469(1990)047<1742:trotma>2.0.co;2   AbstractWebsite

This paper examines the effects of the relationship between cirrus cloud ice water content and cloud temperature on climate change. A simple mechanistic climate model is used to study the feedback between ice water content and temperature. The central question studied in this paper concerns the extent to which both the radiative and microphysical properties of cirrus cloud influence such a feedback. To address this question, a parameterization of the albedo and emissivity of clouds is introduced. Observations that relate the ice water content to cloud temperature are incorporated in the parameterization to introduce a temperature dependence to both albedo and emittance. The cloud properties relevant to the cloud feedback are expressed as functions of particles size re, asymmetry parameter g and cloud temperature and analyses of aircraft measurements, lidar and ground based radiometer data are used to select re and g. It was shown that scattering calculations assuming spherical particles with a distribution described by re = 16 μm reasonably matched the lidar and radiometer data. However, comparison of cloud radiation properties measured from aircraft to those parameterized in this study required values of g significantly smaller than those derived for spheres but consistent with our understanding of nonspherical particle scattering.The climate simulations revealed that the influence of cirrus cloud on climate was strongly affected by the choice of re and g: parameters that are both poorly known for cirrus. It was further shown that the effect of ice water feedback on a CO2 warming simulation could be either positive or negative depending on the value of re assumed. Based on these results, it was concluded that prediction of cirrus cloud feedback on climate is both premature and limited by our lack of understanding of the relationship between size and shape of ice crystals and the gross radiative properties of cirrus.

Flatau, M, Schubert WH, Stevens DE.  1994.  The Role of Baroclinic Processes in Tropical Cyclone Motion - the Influence of Vertical Tilt. Journal of the Atmospheric Sciences. 51:2589-2601.   10.1175/1520-0469(1994)051<2589:trobpi>2.0.co;2   AbstractWebsite

The numerical study presented here focuses on baroclinic processes that contribute to tropical cyclone (TC) propagation. A three-dimensional, semispectral, primitive equation model of baroclinic vortex was developed to study TC motion. In a tiled vortex, interaction between upper- and lower-level vorticity anomalies leads to vortex propagation relative to the steering flow. On a beta plane, with no environmental flow, the vortex is tilted toward the south and the interaction between the layers reduces the westward movement of the vortex. The vortex tilting can also occur due to the vertical shear in the environmental wind. On an f plane, the interaction between the layers causes the northward movement of the vortex in westerly linear shear, and southward movement in easterly linear shear, with a meridional velocity of about 1 ms-1. This velocity increases with increasing vortex intensity and vertical motion.

Flatau, M, Stevens DE.  1993.  The Role of Outflow-Layer Instabilities in Tropical Cyclone Motion. Journal of the Atmospheric Sciences. 50:1721-1733.   10.1175/1520-0469(1993)050<1721:trooli>2.0.co;2   AbstractWebsite

The paper examines the role of the development of outflow-layer instabilities on the motion of tropical cyclones. The influence of barotropic instability is examined by comparing the time changes in the storm tracks with the frequencies of free, unstable barotropic modes. For intense vortices barotropic instability is shown to contribute to the slow (periods of a few days) trochoidal motion of a cyclone. The development of instability depends on the horizontal distribution and frequency of environmental forcing. The strongest response occurs when the frequency of the forcing matches the frequency of an unstable mode.

S
Matthews, AJ, Baranowski DB, Heywood KJ, Flatau PJ, Schmidtko S.  2014.  The surface diurnal warm layer in the Indian Ocean during CINDY/DYNAMO. Journal of Climate. 27:9101-9122.   10.1175/jcli-d-14-00222.1   AbstractWebsite

A surface diurnal warm layer is diagnosed from Seaglider observations and develops on half of the days in the Cooperative Indian Ocean Experiment on Intraseasonal Variability/Dynamics of the Madden-Julian Oscillation (CINDY/DYNAMO) Indian Ocean experiment. The diurnal warm layer occurs on days of high solar radiation flux (>80 W m(-2)) and low wind speed (<6 ms(-1)) and preferentially in the inactive stage of the Madden-Julian oscillation. Its diurnal harmonic has an exponential vertical structure with a depth scale of 4-5m (dependent on chlorophyll concentration), consistent with forcing by absorption of solar radiation. The effective sea surface temperature (SST) anomaly due to the diurnal warm layer often reaches 0.8 degrees C in the afternoon, with a daily mean of 0.2 degrees C, rectifying the diurnal cycle onto longer time scales. This SST anomaly drives an anomalous flux of 4Wm(-2) that cools the ocean. Alternatively, in a climate model where this process is unresolved, this represents an erroneous flux that warms the ocean. A simple model predicts a diurnal warm layer to occur on 30%-50% of days across the tropical warm pool. On the remaining days, with low solar radiation and high wind speeds, a residual diurnal cycle is observed by the Seaglider, with a diurnal harmonic of temperature that decreases linearly with depth. As wind speed increases, this already weak temperature gradient decreases further, tending toward isothermal conditions.