A dynamical analysis of the drag conditions at the sea surface

When applied to a sea surface, shortcomings are noted for the ordinary classification of drag conditions at rigid underlying surfaces according to the Reynolds roughness number Re _s . It is shown that in the case of mobile underlying surfaces, it would be more natural to use the dynamical classification of drag conditions according to the order of magnitude of the ratio ( = /) of the momentum flux toward the waves ( _w) to the viscous momentum flux through the surface ( _w). The relevant estimates of for the main stages of development of the wind waves indicate that the observed values of the drag coefficient of the sea surface correspond to the case of underdeveloped roughness.     

The Askervein Hill Project: Mean wind variations at fixed heights above ground

This is one of a series of papers on the Askervein Hill Project. It presents results on the variations in mean wind speed at fixed heights (z) above the ground from linear arrays of anemometer posts and towers. Most of the data are for z = 10 m but some are for z = 3 m. Selected and directionally grouped data from the 55 Mean Flow runs are presented together with mean flow data from Askervein '83 Turbulence runs. Comparisons are made between the data and guideline estimates of fractional speed-up ratio at hilltop locations and between the data and MS3DJH/3 model predictions along the tower lines. There is good agreement in most cases.     

Examination of ‘global atmospheric temperature monitoring with satellite microwave measurements’: 2. Analysis of satellite data

Using Microwave Sounding Unit (MSU) channel 2 (Ch. 2, 53.74 GHz) data, Spencer and Christy (1992a) determined that the earth exhibits no temperature trend in the period 1979–90, while other authors find a temperature increase of roughly 0.1 K. Based on a theoretical analysis Prabhakara et al. (1995) showed that the information about the global atmospheric temperature deduced from MSU Ch. 2 observations has a small contamination, T ₂, as a result of the attenuation due to hydrometeors in the atmosphere. A method is developed in this study, that utilizes coincident measurements made by MSU in Ch. 1 (50.3 GHz), to estimate this T ₂ over the global oceans. The magnitude of T ₂ is found to be about 1 K over significant parts of the tropical oceanic rain belts and about 0.25 K over minor portions of the mid-latitude oceanic storm tracks. Due to events such as El Niôo, there is variability from year to year in the rain areas and rain intensity leading to significant change in the patterns of T ₂. The patterns of T ₂ derived for March 82 and March 83 reveal such a change. When averaged over the global oceans, from 50° N to 50° S, T ₂ has a value of 0.25 and 0.29 K for March 1982 and 1983, respectively. Due to these reasons the interannual temperature change derived by Spencer and Christy from MSU Ch. 2 will contain a residual hydrometeor effect. Thus in evaluating decadal trend of the global mean temperature of the order of 0.1 K from MSU Ch. 2 data one has to take into account completely the contamination due to hydrometeors.     

Major wet interval in white mountains medieval warm period evidenced inδ 13C of bristlecone pine tree rings

A long ¹³C chronology was developed from bristlecone pine (Pinus longaeva) at the Methuselah Walk site in the White Mountains of California. The chronology represents cellulose from five-year ring groups pooled from multiple radii of multiple trees. The most dramatic isotopic event in the chronology appears from A.D. 1080–1129, when ¹³C values are depressed to levels ~ 2 below the mean for the period A.D. 925–1654. This isotopic excursion appears to represent a real event and is not an artifact of sampling circumstances; in fact, a similar excursion occurs in a previously-reported, independent ¹³C chronology from bristlecone pine. By carbon isotope fractionation models, the shift to low ¹³C values is consistent with abundant soil moisture, permitting leaf stomata to remain open, and allowing ready access of CO₂ from which carbon fixation may discriminate more effectively against¹³C in favor of¹²C. According to this model, the¹³C-depleted 50-yr isotopic excursion represents the wettest period in the White Mountains in the past 1000 yr, during which isotope-reconstructed July Palmer Drought Severity Indices averaged ~ +2.2.     

Tropical Glacier and Ice Core Evidence of Climate Change on Annual to Millennial Time Scales

This paper examines the potential of the stable isotopic ratios, ¹⁸O/¹⁶O ( ¹⁸O_ice)and ²H/¹H ( D_ice), preserved in mid to low latitude glaciers as a toolfor paleoclimate reconstruction. Ice cores are particularly valuable as they contain additional data, such as dust concentrations, aerosol chemistry, and accumulation rates, that can be combined with the isotopic information to assist with inferences about the regional climate conditions prevailing at the time of deposition. We use a collection of multi-proxy ice core histories to explore the ¹⁸O-climate relationship over the last 25,000 years that includes both Late Glacial Stage (LGS) and Holocene climate conditions. These results suggest that on centennial to millennial time scales atmospheric temperature is the principal control on the ¹⁸O_ice of the snowfall that sustains these high mountainice fields.Decadally averaged ¹⁸O_ice records from threeAndean and three Tibetan ice cores are composited to produce a low latitude ¹⁸O_ice history for the last millennium. Comparison ofthis ice core composite with the Northern Hemisphere proxy record (1000–2000A.D.) reconstructed by Mann et al. (1999) and measured temperatures(1856–2000) reported by Jones et al. (1999) suggests the ice cores have captured the decadal scale variability in the global temperature trends. These ice cores show a 20th century isotopic enrichment that suggests a large scale warming is underway at low latitudes. The rate of this isotopically inferred warming is amplified at higher elevations over the Tibetan Plateau while amplification in the Andes is latitude dependent with enrichment (warming) increasing equatorward. In concert with this apparent warming, in situobservations reveal that tropical glaciers are currently disappearing. A brief overview of the loss of these tropical data archives over the last 30 years is presented along with evaluation of recent changes in mean ¹⁸O_ice composition. The isotopic composition of precipitation should be viewed not only as a powerful proxy indicator of climate change, but also as an additional parameter to aid our understanding of the linkages between changes in the hydrologic cycle and global climate.     

The wind-induced shaping and migration of an isolated dune: A numerical experiment

Numerical experiments are carried out to simulate the development and migration of a barchan dune starting with a conical pile of sand. Such an experiment is done in three steps: (1) computation of the steady-state wind field over and around a barchan using the numerical meso-scale simulation model FITNAH, whereby the horizontal variation of the friction velocity is also calculated; (2) computation of the sand transport using the friction velocity in the transport formula by Lettau and Lettau (1978); (3) computation of the erosion and deposition rates as the divergence of the sand transport, where a special treatment is used for the slip-face of the barchan dune. Adding these rates to the height field leads to a different shape of the dune after a time step t _h . Then this procedure has to be repeated for the next time step t _h .The results are in good agreement with observations: the initial pile of sand develops wings (horns) and a slip-face between them. In addition, flow separation over the lee-side can be simulated. Finally, the tendency to form a barchan in equilibrium is considered.     

MCC and gull flight behavior

Wintertime observations of mesoscale cellular convection (MCC) over the East China Sea have resulted in criteria that have a remarkable similarity to those reported by Woodcock (1975) in the study of thermals and gull flight behavior. It has been determined that the surface wind speed (V) and the air-sea temperature difference (T) prescribe unique and compatible conditons for both the occurrence of MCC and soaring by sea gulls. Specifically, the onset of MCC when V is between 5 and 9 m s^–1 is inversely proportional to T in the range 5 to 7 °C. Elsewhere, the onset of MCC occurs under conditions of direct proportionality between V and T. Necessary conditions for the occurrence of MCC due to heating from below are T 5 °C and V 5 m s^–1. The boundaries of the convective regime for MCC are discussed and interpreted in accordance with the regime for sea-gull soaring and similarity concepts.     

A reversing temperature-difference measurement system for Bowen ratio determination

The design and performance of a reversing temperature difference measurement system are reported. This system employs five-junction copper-constantan thermopiles for the measurement of T and T _w, while a linearized thermistor is used to measure T _w. Field performance has been checked against a precision lysimeter as well as against a second temperature difference measurement system in which diodes are used for temperature measurement. In both cases, the agreement between the systems is satisfactory for the measurement of hourly values of the Bowen ratio.     

Numerical integration errors in calculated tropospheric photodissociation rate coefficients

Tropospheric photodissociation rate coefficients (J values) were calculated for NO₂, O₃, HNO₂, CH₂O, and CH₃CHO using high spectral resolution (0.1 mm wavelength increments), and compared to the J values obtained with numerically degraded resolution (=1, 2, 4, 6, 8, and 10 nm, and several commonly used nonuniform grids). Depending on the molecule, substantial errors can be introduced by the larger increments. Thus for =10 nm, errors are less than 1% for NO₂, less than 2% for HNO₂, +6.5% to -16% for CH₂O, -6.9% to +24% for CH₃CHO, and -24% to +110% for O₃. The errors for CH₂O arise from the fine structure of its absorption spectrum, and are prevalently negative (underestimate of J). The errors for O₃, and to a lesser extent for CH₃CHO, arise mainly from under-resolving the overlap of the molecular action spectrum and the tropospheric actinic flux in the wavelength region of stratospheric ozone attenuation. The sign of those errors depends on whether the actinic flux is averaged onto the grid before or after the radiative transfer calculation. In all cases studied, grids with 2 nm produced errors no larger than 5%.     

The Influence of O3, Relative Humidity, NO and NO2 on the Oxidation of α-Pinene and Δ3-Carene

Upto 13% of -pinene and ³-carene had reacted after 213 s in this dark experimental set-up, where O₃, NO and NO₂ were mixed with terpenes at different relative humidities (RHs). The different experiments were planned according to an experimental design, where O₃, NO₂, NO, RH and reaction time were varied between high and low settings (25 and 75 ppb, 15 and 42%, 44 and 213 s). An increased amount of -pinene and ³-carene reacted in the chamber was observed, when the level of O₃, NO and reaction time was increased and RH was decreased. In the study, it was found that different interactions affected the amount of terpene reacted as well. These interactions were between O₃ and NO, O₃ and reaction time, NO and RH, and between NO and reaction time.     

The carbon isotopic ratio of atmospheric carbon dioxide at Tsukuba,Japan

To find out the secular and seasonal trends of the ¹³C value and CO₂ concentration in the surface air and the determination of the ¹³C in the atmospheric CO₂ collected at Tsukuba Science City was carried out during the period from July 1981 to October 1983. The monthly average of the ¹³C value of CO₂ in the surface air collected at 1400 LMT ranged from -7.52 to \s-8.45 with an average of -7.96±0.25 and the CO₂ concentration in the air varied from 334.5 l 1^-1 to 359 l 1^-1 with an average of 347.2±6.3 l 1^-1. The ¹³C value is high in summer and low in winter and is negatively correlated with the CO₂ concentration. In general, the relationship between the ¹³C and the CO₂ concentration is explainable by a simple mixing model of two different constant carbon isotopic species but the relationship does not always follow the model. The correlation between the ¹³C value and the CO₂ concentration is low during the plant growth season and high at other times. The observed negative deviation of the ¹³C value from the simple mixing model in the plant growth season is partly due to the isotopic fractionation process which takes place in the land biota.     

Statistics of Surface-Layer Turbulence Over Terrain with Metre-Scale Heterogeneity

Refuge has patchy vegetation in sandy soil. During midday and at night, the surface sources and sinks for heat and moisture may thus be different. Although the Sevilleta is broad and level, its metre-scale heterogeneity could therefore violate an assumption on which Monin-Obukhov similarity theory (MOST) relies. To test the applicability of MOST in such a setting, we measured the standard deviations of vertical (_w) and longitudinal velocity (_u), temperature (_t), and humidity (_q), the temperature-humidity covariance (¯tq), and the temperature skewness (S_t). Dividing the former five quantities by the appropriate flux scales (u_*, _*, and q_*) yielded the nondimensional statistics _w/u_*, _u/u_*, _t/|t_*|, _q/|q_*|, and ¯tq/t_*q_*. _w/u_*, _t/|t_*|, and S_t have magnitudes and variations with stability similar to those reported in the literature and, thus, seem to obey MOST. Though _u/u_* is often presumed not to obey MOST, our _u/u_* data also agree with MOST scaling arguments. While _q/|q_*| has the same dependence on stability as _t/|t_*|, its magnitude is 28% larger. When we ignore ¯tq/t_*q_* values measured during sunrise and sunset transitions – when MOST is not expected to apply – this statistic has essentially the same magnitude and stability dependence as (_t/t_*)². In a flow that truly obeys MOST, (_t/t_*)², (_q/q_*)², and ¯tq/t_*q_* should all have the same functional form. That (_q/q_*)² differs from the other two suggests that the Sevilleta has an interesting surface not compatible with MOST. The sources of humidity reflect the patchiness while, despite the patchiness, the sources of heat seem uniformly distributed.     

Analytical solutions for the Ekman layer

The PBL equation that governs the transition from the constant-stress surface layer to the geostrophic wind in a neutrally stratified atmosphere for which the eddy viscosityK(z) is assumed to vary smoothly from the surface-layer value U _*z (0.4,U _*=friction velocity,z=elevation) to the geostrophic asymptoteK _GU _*d forzd is solved through an expansion in fd/U _*1 (f=Coriolis parameter). The resulting solution is separated into Ekman's constant-K solution an inner component that reduces to the classical logarithmic form forzd and isO() relative to the Ekman component forzd. The approximationKU _*d is supported by the solution of Nee and Kovasznay's phenomenological transport equation forK(z), which yieldsK–U _*d exp(–z/d), where is an empirical constant for which observation implies, 1. The parametersA andB in Kazanskii and Monin's similarity relation forG/U _* (G=geostrophic velocity) are determined as functions of . The predicted values ofG/U _* and the turning angle are in agreement with the observed values for the Leipzig wind profile. The predicted value ofB based on the assumption of asymptotically constantK is 4.5, while that based on the Nee-Kovasznay model is 5.1; these compare with the observed value of 4.7 for the Leipzig profile. A thermal wind correction, an asymptotic solution for arbitraryK(z) and 1, and an exact (unrestricted ) solution forK(z)=U _*d[1–exp(–z/d)] are developed in appendices.     

Nearly steady convection and the boundary-layer budgets of water vapor and sensible heat

The budgets of water vapor and sensible heat in the convective atmospheric boundary (mixed) layer are analyzed by means of a simple slab approach adapted to steady large-scale advective conditions with radiation and cloud activity. The entrainment flux for sensible heat is assumed to be a linear function of the surface flux. The flux of water vapor at the top of the mixed layer is parameterized by extending the first-order Betts-Deardorff approach, i.e., by adopting linear changes for both the specific humidity and the flux across the mixed layer and across the inversion layer of finite thickness. In this way the dissimilarity of sensible heat and water vapor transport in the mixed layer can be taken into account. The experimental data were obtained from the Air Mass Transformation Experiment (AMTEX). The entrainment constant for sensible heat at the top of the mixed layer was found to have values similar to those observed in other weakly convective situations, i.e., around 0.4 to 0.6. This appears to indicate that the effect of mechanical turbulence was not negligible; however, the inclusion of this effect in the formulation did not improve the correlation. In contrast to the first-order approach, the zero-order approach, i. e., the jump equation commonly used for the flux of a scalar at the inversion, (ovwc ) _h = w_e c (where w _e is the entrainment velocity and c the concentration jump across the inversion), was found to be invalid and incapable of describing the data.     

Markov-chain simulation of particle dispersion in inhomogeneous flows: The mean drift velocity induced by a gradient in Eulerian velocity variance

The Langevin equation is used to derive the Markov equation for the vertical velocity of a fluid particle moving in turbulent flow. It is shown that if the Eulerian velocity variance _wE is not constant with height, there is an associated vertical pressure gradient which appears as a force-like term in the Markov equation. The correct form of the Markov equation is: w(t + t) = aw(t) + b _wE + (1 – a)T _L ( _wE ²)/z, where w(t) is the vertical velocity at time t, a random number from a Gaussian distribution with zero mean and unit variance, T _L the Lagrangian integral time scale for vertical velocity, a = exp(–t/T _L), and b = (1 – a ²)^1/2. This equation can be used for inhomogeneous turbulence in which the mean wind speed, _wE and T _L vary with height. A two-dimensional numerical simulation shows that when this equation is used, an initially uniform distribution of tracer remains uniform.     

Observed variations of σθ and σΦ in the nocturnal drainage flow in a deep valley

The variations of and in the drainage flow in the Brush Creek valley of western Colorado are investigated using data from Doppler acoustic sodars and instrumented towers. The data were obtained on two experimental nights during the 1984 ASCOT field study. There is good agreement between the variations derived from low-level observations of the sodars and those derived from the towers located throughout the valley. The observed hourly average and in the nocturnal drainage flow are about 20 ° to 25 ° and 5 °, respectively; these values are much larger than those generally observed over flat terrain during nighttime stable conditions. After sunrise (about 0600 MST), as the valley warms and the flow direction changes to up-valley, these parameters increase sharply to their peak values at about 0800 MST and then decrease to their normal daytime values after about two hours.In the drainage flow, the hourly average varies inversely with wind speed according to the relation u 0.7ms^-1. The vertical standard deviation is much less enhanced by complex terrain than the horizontal standard deviation. The observed values are predicted fairly well by the local similarity theory.Oak Ridge Associated Universities (ORAU) Summer Research Participant at ATDD in 1987 andOak Ridge Associated Universities (ORAU) Summer Research Participant at ATDD in 1987 and     

The dependence of the measured cool skin of the ocean on wind stress and total heat flux

The temperature drop T between the ocean surface and the 5-cm depth was recorded during GATE, Phase III. With measured values of the total heat flux Q and an assumption about the thickness of the viscous boundary layer of the ocean, the wind-speed dependence of the factor of proportionality between T and Q is determined. This factor depends on the deviations of the thickness of the conductive layer from the thickness of the viscous layer and possibly partially on the wind stress. A further assumption about the thickness of the conductive layer leads to a wind-speed dependence of the ratio between total wind stress and its wave supporting part of it. This ratio increases from a value 1.5 at 1 m s^–1 to 9 at 10 m s^–1, which is in agreement with existing estimates.     

Observation of lake breeze penetration and subsequent development of the thermal internal boundary layer for the Nanticoke II Shoreline Diffusion Experiment

The turbulent structure of the lake breeze penetration and subsequent development of the thermal internal boundary layer (TIBL) was observed using a kytoon-mounted ultrasonic anemometer-thermometer. The lake breeze penetrated with an upward rolling motion associated with the upward flow near the lake breeze front. After the lake breeze front passed, the behaviors of the velocity and temperature at the top of the lake breeze layer were similar to those found in convective boundary layers (CBL). Comparing _gq/_*, _u /w _* and _w /w _* between the present observation of TIBL development after the passage of the lake breeze front and CBL data from the literature, the /_* values showed reasonable agreement; however, _u /w _* and _w /W_* had smaller values in the TIBL than in the CBL at higher altitudes. This is due to the differences in the mean velocity profiles. While the CBL has a uniform velocity profile, the TIBL has a peak at lower elevation due to the lake breeze penetration; the velocity then decreases with height.Present address: The Institute of Behavioral Science, 1-35-7 Yoyogi, Tokyo 151, Japan.     

Equilibrium evaporation beneath a growing convective boundary layer

Expressions for the equilibrium surface Bowen ratio ( _s ) and equilibrium evaporation are derived for a growing convective boundary layer (CBL) in terms of the Bowen ratio at the top of the mixed layer _i and the entrainment parameter A _R . If A_R is put equal to zero, the solution for _s becomes-that previously obtained for the zero entrainment or closed box model. The Priestley-Taylor parameter is also calculated and plotted in terms ofA _R and _i . Realistic combinations of the atmospheric parameters give values of in the range 1.1 to 1.4.     

One-dimensional checkerboards and blending heights

We analyze the checkerboard problem of many alternating surfaces with different properties, on scales up to (say) 3,000 m. Power-law representations of the vertical profiles of mean wind speed and eddy diffusivity lead to solutions in terms of Kelvin and trigonometric functions.These solutions are used to determine blending heights (_*), where deviations from the mean of concentration, or of vertical flux density, fall to some small fraction, , of their value at the surface. Values of _*are important for regional and larger-scale meteorological models. In smaller scale micrometeorological studies, they may serve also as the top levels of surface boundary layers.An important result for both theoretical and experimental contexts is that deviations of flux persist with elevation much more strongly than those of concentration, so that, in general, _* should be based on flux rather than concentration. Representative values of _*, for = 0.05, are of order 5 and 30 m for surface pattern wavelengths of 10² and 10³ m, respectively. Values of _* are robust to changes in adopted power-law indices, and are independent of wind speed. Surface roughness has a mild but calculable effect.