Observational Characterization of the Downward Atmospheric Longwave Radiation at the Surface in the City of Sao Paulo

This work describes the seasonal and diurnal variations of downward longwave atmospheric irradiance (LW) at the surface in Sao Paulo, Brazil, using 5-min-averaged values of LW, air temperature, relative humidity, and solar radiation observed continuously and simultaneously from 1997 to 2006 on a micrometeorological platform, located at the top of a 4-story building. An objective procedure, including 2-step filtering and dome emission effect correction, was used to evaluate the quality of the 9-yr-long LW dataset. The comparison between LW values observed and yielded by the Surface Radiation Budget project shows spatial and temporal agreement, indicating that monthly and annual average values of LW observed in one point of Sao Paulo can be used as representative of the entire metropolitan region of Sao Paulo. The maximum monthly averaged value of the LW is observed during summer (389 +/- 14 W m(-2): January), and the minimum is observed during winter (332 +/- 12 W m(-2); July). The effective emissivity follows the LW and shows a maximum in summer (0.907 +/- 0.032; January) and a minimum in winter (0.818 +/- 0.029; June). The mean cloud effect, identified objectively by comparing the monthly averaged values of the LW during clear-sky days and all-sky conditions, intensified the monthly average LW by about 32.0 +/- 3.5 W m(-2) and the atmospheric effective emissivity by about 0.088 +/- 0.024. In August, the driest month of the year in Sao Paulo, the diurnal evolution of the LW shows a minimum (325 +/- 11 W m(-2)) at 0900 LT and a maximum (345 12 W m-2) at 1800 LT, which lags behind (by 4 h) the maximum diurnal variation of the screen temperature. The diurnal evolution of effective emissivity shows a minimum (0.781 +/- 0.027) during daytime and a maximum (0.842 +/- 0.030) during nighttime. The diurnal evolution of all-sky condition and clear-sky day differences in the effective emissivity remain relatively constant (7% +/- 1%), indicating that clouds do not change the emissivity diurnal pattern. The relationship between effective emissivity and screen air temperature and between effective emissivity and water vapor is complex. During the night, when the planetary boundary layer is shallower, the effective emissivity can be estimated by screen parameters. During the day, the relationship between effective emissivity and screen parameters varies from place to place and depends on the planetary boundary layer process. Because the empirical expressions do not contain enough information about the diurnal variation of the vertical stratification of air temperature and moisture in Sao Paulo, they are likely to fail in reproducing the diurnal variation of the surface emissivity. The most accurate way to estimate the LW for clear-sky conditions in Sao Paulo is to use an expression derived from a purely empirical approach.

The Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CATT-BRAMS) - Part 1: Model description and evaluation

We introduce the Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CATT-BRAMS). CATT-BRAMS is an on-line transport model fully consistent with the simulated atmospheric dynamics. Emission sources from biomass burning and urban-industrial-vehicular activities for trace gases and from biomass burning aerosol particles are obtained from several published datasets and remote sensing information. The tracer and aerosol mass concentration prognostics include the effects of sub-grid scale turbulence in the planetary boundary layer, convective transport by shallow and deep moist convection, wet and dry deposition, and plume rise associated with vegetation fires in addition to the grid scale transport. The radiation parameterization takes into account the interaction between the simulated biomass burning aerosol particles and short and long wave radiation. The atmospheric model BRAMS is based on the Regional Atmospheric Modeling System (RAMS), with several improvements associated with cumulus convection representation, soil moisture initialization and surface scheme tuned for the tropics, among others. In this paper the CATT-BRAMS model is used to simulate carbon monoxide and particulate material (PM(2.5)) surface fluxes and atmospheric transport during the 2002 LBA field campaigns, conducted during the transition from the dry to wet season in the southwest Amazon Basin. Model evaluation is addressed with comparisons between model results and near surface, radiosondes and airborne measurements performed during the field campaign, as well as remote sensing derived products. We show the matching of emissions strengths to observed carbon monoxide in the LBA campaign. A relatively good comparison to the MOPITT data, in spite of the fact that MOPITT a priori assumptions imply several difficulties, is also obtained.

Transport of North African dust from the Bod,l, depression to the Amazon Basin: a case study

Through long-range transport of dust, the North-African desert supplies essential minerals to the Amazon rain forest. Since North African dust reaches South America mostly during the Northern Hemisphere winter, the dust sources active during winter are the main contributors to the forest. Given that the Bod,l, depression area in southwestern Chad is the main winter dust source, a close link is expected between the Bod,l, emission patterns and volumes and the mineral supply flux to the Amazon. Until now, the particular link between the Bod,l, and the Amazon forest was based on sparse satellite measurements and modeling studies. In this study, we combine a detailed analysis of space-borne and ground data with reanalysis model data and surface measurements taken in the central Amazon during the Amazonian Aerosol Characterization Experiment (AMAZE-08) in order to explore the validity and the nature of the proposed link between the Bod,l, depression and the Amazon forest. This case study follows the dust events of 11-16 and 18-27 February 2008, from the emission in the Bod,l, over West Africa (most likely with contribution from other dust sources in the region) the crossing of the Atlantic Ocean, to the observed effects above the Amazon canopy about 10 days after the emission. The dust was lifted by surface winds stronger than 14 m s(-1), usually starting early in the morning. The lofted dust, mixed with biomass burning aerosols over Nigeria, was transported over the Atlantic Ocean, and arrived over the South American continent. The top of the aerosol layer reached above 3 km, and the bottom merged with the boundary layer. The arrival of the dusty air parcel over the Amazon forest increased the average concentration of aerosol crustal elements by an order of magnitude.

Evaluation of the carbon content of aerosols from the burning of biomass in the Brazilian Amazon using thermal, optical and thermal-optical analysis methods

Aerosol samples were collected at a pasture site in the Amazon Basin as part of the project LBA-SMOCC-2002 (Large-Scale Biosphere-Atmosphere Experiment in Amazonia - Smoke Aerosols, Clouds, Rainfall and Climate: Aerosols from Biomass Burning Perturb Global and Regional Climate). Sampling was conducted during the late dry season, when the aerosol composition was dominated by biomass burning emissions, especially in the submicron fraction. A 13-stage Dekati low-pressure impactor (DLPI) was used to collect particles with nominal aerodynamic diameters (D(p)) ranging from 0.03 to 0.10 mu m. Gravimetric analyses of the DLPI substrates and filters were performed to obtain aerosol mass concentrations. The concentrations of total, apparent elemental, and organic carbon (TC, EC(a), and OC) were determined using thermal and thermal-optical analysis (TOA) methods. A light transmission method (LTM) was used to determine the concentration of equivalent black carbon (BC(e)) or the absorbing fraction at 880 nm for the size-resolved samples. During the dry period, due to the pervasive presence of fires in the region upwind of the sampling site, concentrations of fine aerosols (D(p) < 2.5 mu m: average 59.8 mu g m(-3)) were higher than coarse aerosols (D(p) > 2.5 mu m: 4.1 mu g m(-3)). Carbonaceous matter, estimated as the sum of the particulate organic matter (i.e., OC x 1.8) plus BC(e), comprised more than 90% to the total aerosol mass. Concentrations of EC(a) (estimated by thermal analysis with a correction for charring) and BC(e) (estimated by LTM) averaged 5.2 +/- 1.3 and 3.1 +/- 0.8 mu g m(-3), respectively. The determination of EC was improved by extracting water-soluble organic material from the samples, which reduced the average light absorption Angstrom exponent of particles in the size range of 0.1 to 1.0 mu m from >2.0 to approximately 1.2. The size-resolved BC(e) measured by the LTM showed a clear maximum between 0.4 and 0.6 mu m in diameter. The concentrations of OC and BC(e) varied diurnally during the dry period, and this variation is related to diurnal changes in boundary layer thickness and in fire frequency.

Hybrid modeling of convective laminar flow in a permeable tube associated with the cross-flow process

The confined flows in tubes with permeable surfaces arc associated to tangential filtration processes (microfiltration or ultrafiltration). The complexity of the phenomena do not allow for the development of exact analytical solutions, however, approximate solutions are of great interest for the calculation of the transmembrane outflow and estimate of the concentration, polarization phenomenon. In the present work, the generalized integral transform technique (GITT) was employed in solving the laminar and permanent flow in permeable tubes of Newtonian and incompressible fluid. The mathematical formulation employed the parabolic differential equation of chemical species conservation (convective-diffusive equation). The velocity profiles for the entrance region flow, which are found in the connective terms of the equation, were assessed by solutions obtained from literature. The velocity at the permeable wall was considered uniform, with the concentration at the tube wall regarded as variable with an axial position. A computational methodology using global error control was applied to determine the concentration in the wall and concentration boundary layer thickness. The results obtained for the local transmembrane flux and the concentration boundary layer thickness were compared against others in literature. (C) 2007 Elsevier B.V. All rights reserved.

CONCENTRATION FIELDS NEAR AIR-WATER INTERFACES DURING INTERFACIAL MASS TRANSPORT: OXYGEN TRANSPORT AND RANDOM SQUARE WAVE ANALYSIS

Mass transfer across a gas-liquid interface was studied theoretically and experimentally, using transfer of oxygen into water as the gas-liquid system. The experimental results support the conclusions of a theoretical description of the concentration field that uses random square waves approximations. The effect of diffusion over the concentration records was quantified. It is shown that the peak of the normalized rills concentration fluctuation profiles must be lower than 0.5, and that the position of the peak of the rms value is an adequate measure of the thickness of the diffusive layer. The position of the peak is the boundary between the regions more subject to molecular diffusion or to turbulent transport of dissolved mass.

Enhanced flow visualisation using near-resonant holographic interferometry

Holographic interferometry measurements have been performed on high-speed, high-temperature gas flows with a laser output tuned near a resonant sodium transition. The technique allows the detection and quantification of the sodium concentration in the flow. By controlling the laser detuning and seeded sodium concentration, we performed flow visualization in low-density flows that are not normally detectable with standard interferometry. The technique was also successfully used to estimate the temperature in the boundary layer of the flow over a flat plate.

Turbulence in Small Sub-Tropical Estuary with Semi-Diurnal Tides

In natural estuaries, contaminant transport is driven by the turbulent momentum mixing. The predictions of scalar dispersion can rarely be predicted accurately because of a lack of fundamental understanding of the turbulence structure in estuaries. Herein detailed turbulence field measurements were conducted at high frequency and continuously for up to 50 hours per investigation in a small subtropical estuary with semi-diurnal tides. Acoustic Doppler velocimetry was deemed the most appropriate measurement technique for such small estuarine systems with shallow water depths (less than 0.5 m at low tides), and a thorough post-processing technique was applied. The estuarine flow is always a fluctuating process. The bulk flow parameters fluctuated with periods comparable to tidal cycles and other large-scale processes. But turbulence properties depended upon the instantaneous local flow properties. They were little affected by the flow history, but their structure and temporal variability were influenced by a variety of mechanisms. This resulted in behaviour which deviated from that for equilibrium turbulent boundary layer induced by velocity shear only. A striking feature of the data sets is the large fluctuations in all turbulence characteristics during the tidal cycle. This feature was rarely documented, but an important difference between the data sets used in this study from earlier reported measurements is that the present data were collected continuously at high frequency during relatively long periods. The findings bring new lights in the fluctuating nature of momentum exchange coefficients and integral time and length scales. These turbulent properties should not be assumed constant.

Characteristics of undular hydraulic jumps: Experiments and analysis

The writers measured velocity, pressure and energy distributions, wavelengths, and wave amplitudes along undular jumps in a smooth rectangular channel 0.25 m wide. In each case the upstream flow was a fully developed shear flow. Analysis of the data shows that the jump has strong three-dimensional features and that the aspect ratio of the channel is an important parameter. Energy dissipation on the centerline is far from negligible and is largely constrained to the reach between the start of the lateral shock waves and the first wave crest of the jump, in which the boundary layer develops under a strong adverse pressure gradient. A Boussinesq-type solution of the free-surface profile, velocity, and energy and pressure distributions is developed and compared with the data. Limitations of the two-dimensional analysis are discussed.

Longevity and stability of cratonic lithosphere: Insights from numerical simulations of coupled mantle convection and continental tectonics

[1] The physical conditions required to provide for the tectonic stability of cratonic crust and for the relative longevity of deep cratonic lithosphere within a dynamic, convecting mantle are explored through a suite of numerical simulations. The simulations allow chemically distinct continents to reside within the upper thermal boundary layer of a thermally convecting mantle layer. A rheologic formulation, which models both brittle and ductile behavior, is incorporated to allow for plate-like behavior and the associated subduction of oceanic lithosphere. Several mechanisms that may stabilize cratons are considered. The two most often invoked mechanisms, chemical buoyancy and/or high viscosity of cratonic root material, are found to be relatively ineffective if cratons come into contact with subduction zones. High root viscosity can provide for stability and longevity but only within a thick root limit in which the thickness of chemically distinct, high-viscosity cratonic lithosphere exceeds the thickness of old oceanic lithosphere by at least a factor of 2. This end-member implies a very thick mechanical lithosphere for cratons. A high brittle yield stress for cratonic lithosphere as a whole, relative to oceanic lithosphere, is found to be an effective and robust means for providing stability and lithospheric longevity. This mode does not require exceedingly deep strength within cratons. A high yield stress for only the crustal or mantle component of the cratonic lithosphere is found to be less effective as detachment zones can then form at the crust-mantle interface which decreases the longevity potential of cratonic roots. The degree of yield stress variations between cratonic and oceanic lithosphere required for stability and longevity can be decreased if cratons are bordered by continental lithosphere that has a relatively low yield stress, i.e., mobile belts. Simulations that combine all the mechanisms can lead to crustal stability and deep root longevity for model cratons over several mantle overturn times, but the dominant stabilizing factor remains a relatively high brittle yield stress for cratonic lithosphere.

Satellite observation and climate system model simulation of the St. Lawrence Island polynya

The St. Lawrence Island polynya (SLIP) is a commonly occurring winter phenomenon in the Bering Sea, in which dense saline water produced during new ice formation is thought to flow northward through the Bering Strait to help maintain the Arctic Ocean halocline. Winter darkness and inclement weather conditions have made continuous in situ and remote observation of this polynya difficult. However, imagery acquired from the European Space Agency ERS-1 Synthetic Aperture Radar (SAR) has allowed observation of the St. Lawrence Island polynya using both the imagery and derived ice displacement products. With the development of ARCSyM, a high resolution regional model of the Arctic atmosphere/sea ice system, simulation of the SLIP in a climate model is now possible. Intercomparisons between remotely sensed products and simulations can lead to additional insight into the SLIP formation process. Low resolution SAR, SSM/I and AVHRR infrared imagery for the St. Lawrence Island region are compared with the results of a model simulation for the period of 24-27 February 1992. The imagery illustrates a polynya event (polynya opening). With the northerly winds strong and consistent over several days, the coupled model captures the SLIP event with moderate accuracy. However, the introduction of a stability dependent atmosphere-ice drag coefficient, which allows feedbacks between atmospheric stability, open water, and air-ice drag, produces a more accurate simulation of the SLIP in comparison to satellite imagery. Model experiments show that the polynya event is forced primarily by changes in atmospheric circulation followed by persistent favorable conditions: ocean surface currents are found to have a small but positive impact on the simulation which is enhanced when wind forcing is weak or variable.

Infiltration effects on sediment mobility under waves

Sediment mobility measurements with a horizontal sand bed under non-breaking waves are reported. Conditions include no seepage and steady downward seepage corresponding to head gradients up to 2.5. The results indicate that infiltration tends to inhibit sediment mobility for a horizontal bcd of 0.2 mm quartz sand exposed to moderated wave induced bed shear stresses. The effect is weak for the parameter range of the present study. The two opposing effects of shear stress increase due to boundary layer thinning and the stabilizing downward drag are successfully accounted for through the modified Shields parameter of Nielsen [Nielsen, P., 1997. Coastal groundwater dynamics. Proc. Coastal Dynamics '97, Plymouth, ASCE, Dp, 546-555] using coefficients derived from independent studies. That is, from the shear stress experiments of Conley [Conley, D.C., 1993. Ventilated oscillatory boundary layers. PhD Thesis, University of California, San Diego, 74 pp.] and the slope stability experiments of Martin and Aral [Martin, C.S. and M.M. Aral, 1971. Seepage force on interfacial bed particles. J. Hydraulics Div., proc. ASCE, Vol. 97, No. Hy7, pp. 1081-1100]. (C) 2001 Elsevier Science B.V. All rights reserved.

Generalized quasi mode theory of macroscopic canonical quantization in cavity quantum electrodynamics and quantum optics II. Application to reflection and refraction at a dielectric interface

The generalization of the quasi mode theory of macroscopic quantization in quantum optics and cavity QED presented in the previous paper, is applied to provide a fully quantum theoretic derivation of the laws of reflection and refraction at a boundary. The quasi mode picture of this process involves the annihilation of a photon travelling in the incident region quasi mode, and the subsequent creation of a photon in either the incident region or transmitted region quasi modes. The derivation of the laws of reflection and refraction is achieved through the dual application of the quasi mode theory and a quantum scattering theory based on the Heisenberg picture. Formal expressions from scattering theory are given for the reflection and transmission coefficients. The behaviour of the intensity for a localized one photon wave packet coming in at time minus infinity from the incident direction is examined and it is shown that at time plus infinity, the light intensity is only significant where the classical laws of reflection and refraction predict. The occurrence of both refraction and reflection is dependent upon the quasi mode theory coupling constants between incident and transmitted region quasi modes being nonzero, and it is seen that the contributions to such coupling constants come from the overlap of the mode functions in the boundary layer region, as might be expected from a microscopic theory.

Shear stress and sediment transport calculations for swash zone modelling

It is shown that the observed difference in sediment transporting efficiency by the swash uprush, compared with the downrush, could be mainly due to greater bed shear stress for a given velocity in the more abruptly accelerated uprush. The bed shear stress generated by an arbitrary free stream velocity time series is modelled in terms of usual wave boundary layer models plus a phase lead (phi(tau) of the bed shear stress compared with the free stream velocity at the peak frequency. With this approach, the total transport amounts in uprush and downrush can be modelled satisfactorily with the same sediment transport formula, without the need for different uprush and downrush coefficients. While the adaptation of sediment transport formulae from steady flow can thus lead to the right total amounts of sediment moved by this method, the timing of the instantaneous sediment transport rates are probably not accurately modelled due to the highly unsteady nature of the swash and the presence of pre-suspended sediment in the uprush. Nevertheless, the proposed method is a useful intermediate step before we have a complete understanding of sediment transport under very rapid accelerations and of the relative contribution of pre-suspended sediment to the onshore sediment transport in swash zones. (C) 2002 Published by Elsevier Science B.V.

Shear stress and sediment transport calculations for sheet flow under waves

A simple method is provided for calculating transport rates of not too fine (d(50) greater than or equal to 0.20 mm) sand under sheet flow conditions. The method consists of a Meyer-Peter-type transport formula operating on a time-varying Shields parameter, which accounts for both acceleration-asymmetry and boundary layer streaming. While velocity moment formulae, e.g.., = Constant x calibrated against U-tube measurements, fail spectacularly under some real waves (Ribberink, J.S., Dohmen-Janssen, C.M., Hanes, D.M., McLean, S.R., Vincent, C., 2000. Near-bed sand transport mechanisms under waves. Proc. 27th Int. Conf. Coastal Engineering, Sydney, ASCE, New York, pp. 3263-3276, Fig. 12), the new method predicts the real wave observations equally well. The reason that the velocity moment formulae fail under these waves is partly the presence of boundary layer streaming and partly the saw-tooth asymmetry, i.e., the front of the waves being steeper than the back. Waves with saw-tooth asymmetry may generate a net landward sediment transport even if = 0, because of the more abrupt acceleration under the steep front. More abrupt accelerations are associated with thinner boundary layers and greater pressure gradients for a given velocity magnitude. The two real wave effects are incorporated in a model of the form Q(s)(t) = Q(s)[theta(t)] rather than Q(S)(t) = Q(S)[u(infinity)(t)], i.e., by expressing the transport rate in terms of an instantaneous Shields parameter rather than in terms of the free stream velocity, and accounting for both streaming and accelerations in the 0(t) calculations. The instantaneous friction velocities u(*)(t) and subsequently theta(t) are calculated as follows. Firstly, a linear filter incorporating the grain roughness friction factor f(2.5) and a phase angle phi(tau) is applied to u(infinity)(t). This delivers u(*)(t) which is used to calculate an instantaneous grain roughness Shields parameter theta(2.5)(t). Secondly, a constant bed shear stress is added which corresponds to the streaming related bed shear stress -rho ($) over bar((u) over tilde(w) over tilde)(infinity) . The method can be applied to any u(infinity)(t) time series, but further experimental validation is recommended before application to conditions that differ strongly from the ones considered below. The method is not recommended for rippled beds or for sheet flow with typical prototype wave periods and d(50) < 0.20 turn. In such scenarios, time lags related to vertical sediment movement become important, and these are not considered by the present model. (C) 2002 Elsevier Science B.V. All rights reserved.