Turbulent heat fluxes by profile and inertial dissipation methods: analysis of the atmospheric surface layer from shipboard measurements during the SOFIA/ASTEX and SEMAPHORE experiments

Heat flux estimates obtained using the inertial dissipation method, and the profile method applied to radiosonde soundings, are assessed with emphasis on the parameterization of the roughness lengths for temperature and specific humidity. Results from the inertial dissipation method show a decrease of the temperature and humidity roughness lengths for increasing neutral wind speed, in agreement with previous studies. The sensible heat flux estimates were obtained using the temperature estimated from the speed of sound determined by a sonic anemometer. This method seems very attractive for estimating heat fluxes over the ocean. However allowance must be made in the inertial dissipation method for non-neutral stratification. The SOFIA/ASTEX and SEMAPHORE results show that, in unstable stratification, a term due to the transport terms in the turbulent kinetic energy budget, has to be included in order to determine the friction velocity with better accuracy. Using the profile method with radiosonde data, the roughness length values showed large scatter. A reliable estimate of the temperature roughness length could not be obtained. The humidity roughness length values were compatible with those found using the inertial dissipation method.     

Vertical velocity field in the convective boundary layer as observed with an acoustic Doppler sodar

The morning development of the daytime convective boundary layer (CBL) during fine weather has been observed with an acoustic Doppler sodar of the C.R.P.E. In particular, the vertical profile of the vertical velocity third-order statistic W* ³ has been obtained. This quantity is a maximum near 0.3z _I where z _I, is the height of the CBL. The histogram of vertical velocity in the CBL shows a relationship between W ³ and the convective velocity W _* and is useful for convective plume determination.     

Density calculations for silicate liquids. I. Revised method for aluminosilicate compositions

Equations are developed for calculating the density of aluminosilicate liquids as a function of composition and temperature. The mean molar volume at reference temperature T_r, is given by $\text{V}{}_{\mathrm{r}}\text{= ∑X}{}_{\mathrm{i}}\text{V}\text{?}{}^{\mathrm{o}}{}_{\mathrm{i}}\text{+ X}{}_{\mathrm{A}}\text{V}\text{?}{}^{\mathrm{o}}{}_{\mathrm{A}}$, where the summation is taken over all oxide components except A1₂O₃, X stands for mole fraction, terms are constants derived independently from an analysis of volume-composition relations in alumina-free silicate liquids, and $\text{V}\text{?}{}^{\mathrm{o}}{}_{\mathrm{A}}$ is the composition-dependent apparent partial molar volume of Al₂O₃. The thermal expansion coefficient of aluminosilicate liquids is given by $\text{α = ∑X}{}_{\mathrm{i}}\text{}\text{?}\text{ga}{}_{\mathrm{i}}{}^{\mathrm{o}}\text{+ X}{}_{\mathrm{A}}\text{}\text{?}\text{ga}{}_{\mathrm{A}}{}^{\mathrm{o}}$, where $\text{}\text{?}\text{ga}{}_{\mathrm{i}}{}^{\mathrm{o}}$ terms are constants independent of temperature and composition, and $\text{}\text{?}\text{ga}{}^{\mathrm{o}}{}_{\mathrm{A}}$ is a composition-dependent term representing the effect of Al₂O₃ on the thermal expansion. Parameters necessary to calculate the volume of silicate liquids at any temperature T according to V(T) = V_rexp[α(T-T_r)], where T_r = 1400°C have been evaluated by least-square analysis of selected density measurements in aluminosilicate melts. Mean molar volumes of aluminosilicate liquids calculated according to the model equation conform to experimentally measured volumes with a root mean square difference of 0.28 $\text{cc}\text{mole}$ and an average absolute difference of 0.90% for 248 experimental observations. The compositional dependence of $\text{V}\text{?}{}^{\mathrm{o}}{}_{\mathrm{A}}$ is discussed in terms of several possible interpretations of the structural role of Al³⁺ in aluminosilicate melts.     

Convective plumes in the atmospheric boundary layer as observed with an acoustic Doppler sodar

The vertical velocity field and the convective plumes in the atmospheric boundary layer have been observed during morning hours with the acoustic Doppler sounder of the C.R.P.E. A method for plume determination using acoustic soundings in the well-mixed layer is presented. Using Telford's 1970 and Manton's 1975 models, a comparison is made between the predictions of the models and the plume properties as observed by the Doppler sodar. The mean plume velocity is found to be parabolic. It is shown, restricting Monin and Obukhov similarity to conditions inside plumes and using only vertical velocity within plumes, that the observed convective plumes carry nearly sixty percent of the sensible heat flux at the top of the surface layer.     

Study of the air-sea interactions at the mesoscale: the SEMAPHORE experiment

The SEMAPHORE (Structure des Echanges Mer-Atmosphère, Propriétés des Hétérogénéités Océaniques: Recherche Expérimentale) experiment has been conducted from June to November 1993 in the Northeast Atlantic between the Azores and Madeira. It was centered on the study of the mesoscale ocean circulation and air-sea interactions. The experimental investigation was achieved at the mesoscale using moorings, floats, and ship hydrological survey, and at a smaller scale by one dedicated ship, two instrumented aircraft, and surface drifting buoys, for one and a half month in October-November (IOP: intense observing period). Observations from meteorological operational satellites as well as spaceborne microwave sensors were used in complement. The main studies undertaken concern the mesoscale ocean, the upper ocean, the atmospheric boundary layer, and the sea surface, and first results are presented for the various topics. From data analysis and model simulations, the main characteristics of the ocean circulation were deduced, showing the close relationship between the Azores front meander and the occurrence of Mediterranean water lenses (meddies), and the shift between the Azores current frontal signature at the surface and within the thermocline. Using drifting buoys and ship data in the upper ocean, the gap between the scales of the atmospheric forcing and the oceanic variability was made evident. A 2 °C decrease and a 40-m deepening of the mixed layer were measured within the IOP, associated with a heating loss of about 100 W m⁻². This evolution was shown to be strongly connected to the occurrence of storms at the beginning and the end of October. Above the surface, turbulent measurements from ship and aircraft were analyzed across the surface thermal front, showing a 30% difference in heat fluxes between both sides during a 4-day period, and the respective contributions of the wind and the surface temperature were evaluated. The classical momentum flux bulk parameterization was found to fail in low wind and unstable conditions. Finally, the sea surface was investigated using airborne and satellite radars and wave buoys. A wave model, operationally used, was found to get better results compared with radar and wave-buoy measurements, when initialized using an improved wind field, obtained by assimilating satellite and buoy wind data in a meteorological model. A detailed analysis of a 2-day period showed that the swell component, propagating from a far source area, is underestimated in the wave model. A data base has been created, containing all experimental measurements. It will allow us to pursue the interpretation of observations and to test model simulations in the ocean, at the surface and in the atmospheric boundary layer, and to investigate the ocean-atmosphere coupling at the local and mesoscales.     

Precipitation of auroral particles in the central polar cap during a sudden commencement magnetic storm

All-sky camera observations from two stations in the inner (northern) polar cap and an auroral zone station are combined with photometer records from the polar cap station Nord in a study of the brilliant auroral display following the ssc of the storm of 7 November 1970. This display is the large, poleward expanding bulge of a substorm triggered by the ssc. It is composed of brilliant discrete forms embedded in low-intensity diffuse electron and proton aurora. The poleward edge of the diffuse electron aurora is 5° north of the discrete auroras and 3° north of the proton aurora. The intensity of the discrete aurora varies as the strength of the auroral electrojet as shown by magnetograms from auroral zone stations. Succeeding the retreating display a subvisible low-energy electron precipitation, which may be identified as the polar squall (Winningham and Heikkila, 1974) is observed over the polar cap during the main phase of the storm.In the early morning sector already existing diffuse auroras broaden towards the equator from the time of the ssc and at least during the following half hour.Ssc-triggered displays have been found (Feldstein, 1959) to withdraw from the inner polar cap as the initial (positive H) phase of the storm ends. A comparison of the records from seven low-latitude stations shows that during this particular storm the positive phase appears to be composed by two overlapping disturbances, i.e. the proper initial phase, which is generally thought to be due to compression of the inner magnetosphere and a series of positive bays accompanying the negative bays in auroral latitudes. These positive bays are observable over a great range of longitudes with a maximum of amplitude near midnight. As judged from the dayside magnetograms the initial (compression) phase ends at an early stage of the substorm. The observed coincidence between the withdrawal of the display and the cessation of the positive H phase of the storm is a consequence of the fact that the second component—the positive bays—and the auroral display over the polar cap are both signatures of the substorm activity.     

Lunar polymict breccia 14321: a petrographic study

Seven petrographic thin sections of lunar rock sample 14321, ‘Big Bertha’, have been examined. It is a complex rock incorporating diverse lithic and single crystal fragments and represents a sampling of the heterogeneous Fra Mauro formation, considered by the writers to be lithified debris from the Imbrium impact event. Electron probe microanalysis and microscopic study of textures reveal the assembly history of this breccia which in turn allows some interpretation of the nature of the pre-Imbrium crust and the effect of the Imbrium impact and the subsequent transportation to the Apollo 14 site. The present-day polymict breccia 14321 is composed of basaltic clasts originating from the fragmentation of a single or closely related set of lava cooling units, a set of fragmental clasts designated as microbreccia 3 (themselves polymict microbreccias), and a light colored matrix which formed rock 14321 by cementing the two major groups of clasts. The light colored matrix material is derived from the fragmentation and mutual abrasion of the basalt and microbreccia 3. On the basis of consistent textural relations two older sets of microbreccias have been identified within microbreccia 3. Microbreccia 1 clasts are well-rounded, relatively light colored, and noritic. They are always completely enclosed within microbreccia 3, most often forming the central cores of rounded accretionary lapilli structures which we have designated as microbreccia 2. Microbreccias 1, 2, 3, and macrobreccia 14321 represent a chronological series of fragmentation and lithification events. Each of these events involved some thermal and/or shock metamorphism as evidenced by mineralogical and textural criteria, and the chronological order of formation of the breccias also corresponds to a decreasing intensity of associated thermal effects. The petrology and mineralogy of 14321 are described in detail in this paper. A more general interpretation of the combined petrographic and chemical data is given in Duncanet al. (1975a).     

Accuracy and efficiency of time integration methods for 1D diffusive wave equation

The diffusive wave approximation of the Saint-Venant equations is commonly used in hydrological models to describe surface flow processes. Numerous numerical approaches can be used to solve this highly nonlinear equation. Nonlinear time integration schemes—also called methods of lines (MOL)—were proven very efficient to solve other nonlinear problems in geosciences but were never considered to deal with surface flow modeling with the diffusive wave equation. In this paper, we study the relative performance of different time and space integration schemes by comparing the results obtained with classical approaches and with nonlinear time integration approaches. The results show that (i) the integration method with a higher order in space shows high accuracy regarding an integrated indicator such as the global mass balance error but is less accurate regarding local indicators, and (ii) nonlinear time integration techniques perform better than classical ones. Overall, it seems that integration techniques combining nonlinear time integration and a low spatial order need to be considered when developing hydrological modeling tools owing to their simplicity of implementation and very good performance.     

A low-dimensional subsurface model for saturated and unsaturated flow processes: ability to address heterogeneity

A low-dimensional model that describes both saturated and unsaturated flow processes in a single equation is presented. Subsurface flow processes in the groundwater, the vadose zone, and the capillary fringe are accounted for through the computation of aggregated hydrodynamic parameters that result from the integration of the governing flow equations from the bedrock to the land surface. The three-dimensional subsurface flow dynamics are thus described by a two-dimensional equation, allowing for a drastic reduction of model unknowns and simplification of the model parameterizations. This approach is compared with a full resolution of the Richards equation in different synthetic test cases. Because the model reduction stems from the vertical integration of the flow equations, the test cases all use different configurations of heterogeneity for vertical cross-sections of a soil-aquifer system. The low-dimensional flow model shows strong consistency with results from a complete resolution of the Richards equation for both the water table and fluxes. The proposed approach is therefore well suited to the accurate reproduction of complex subsurface flow processes.