A Mediterranean nocturnal heavy rainfall and tornadic event. Part I: Overview,damage survey and radar analysis

This study presents an analysis of a severe weather case that took place during the early morning of the 2nd of November 2008, when intense convective activity associated with a rapidly evolving low pressure system affected the southern coast of Catalonia (NE Spain). The synoptic framework was dominated by an upper level trough and an associated cold front extending from Gibraltar along the Mediterranean coast of the Iberian Peninsula to SE France, which moved north-eastward. South easterly winds in the north of the Balearic Islands and the coast of Catalonia favoured high values of 0–3 km storm relative helicity which combined with moderate MLCAPE values and high shear favoured the conditions for organized convection. A number of multicell storms and others exhibiting supercell features, as indicated by Doppler radar observations, clustered later in a mesoscale convective system, and moved north-eastwards across Catalonia. They produced ground-level strong damaging wind gusts, an F2 tornado, hail and heavy rainfall. Total lightning activity (intra-cloud and cloud to ground flashes) was also relevant, exhibiting several classical features such as a sudden increased rate before ground level severe damage, as discussed in a companion study. Remarkable surface observations of this event include 24 h precipitation accumulations exceeding 100 mm in four different observatories and 30 minute rainfall amounts up to 40 mm which caused local flash floods. As the convective system evolved northward later that day it also affected SE France causing large hail, ground level damaging wind gusts and heavy rainfall.     

Loop Current Eddy formation and baroclinic instability

The formation of three Loop Current Eddies, Ekman, Franklin, and Hadal, during the period April 2009 through November 2011 was observed by an array of moored current meters and bottom mounted pressure equipped inverted echo sounders. The array design, areal extent nominally 89° W to 85° W, 25° N to 27° N with 30–50 km mesoscale resolution, permits quantitative mapping of the regional circulation at all depths. During Loop Current Eddy detachment and formation events, a marked increase in deep eddy kinetic energy occurs coincident with the growth of a large-scale meander along the northern and eastern parts of the Loop Current. Deep eddies develop in a pattern where the deep fields were offset and leading upper meanders consistent with developing baroclinic instability. The interaction between the upper and deep fields is quantified by evaluating the mean eddy potential energy budget. Largest down-gradient heat fluxes are found along the eastern side of the Loop Current. Where strong, the horizontal down-gradient eddy heat flux (baroclinic conversion rate) nearly balances the vertical down-gradient eddy heat flux indicating that eddies extract available potential energy from the mean field and convert eddy potential energy to eddy kinetic energy.     

Analysis of rainfall characteristics of the Madden–Julian oscillation using TRMM satellite data

Rainfall characteristics of the Madden–Julian oscillation (MJO) are analyzed primarily using tropical rainfall measuring mission (TRMM) precipitation radar (PR), TRMM microwave imager (TMI) and lighting imaging sensor (LIS) data. Latent heating structure is also examined using latent heating data estimated with the spectral latent heating (SLH) algorithm.The zonal structure, time evolution, and characteristic stages of the MJO precipitation system are described. Stratiform rain fraction increases with the cloud activity, and the amplitude of stratiform rain variation associated with the MJO is larger than that of convective rain by a factor of 1.7. Maximum peaks of both convective rain and stratiform rain precede the minimum peak of the outgoing longwave radiation (OLR) anomaly which is often used as a proxy for the MJO convection. Stratiform rain remains longer than convective rain until ∼4000 km behind the peak of the mature phase. The stratiform rain contribution results in the top-heavy heating profile of the MJO.Associated with the MJO, there are tri-pole convective rain top heights (RTH) at 10–11, ∼7 and ∼3 km, corresponding to the dominance of afternoon showers, organized systems, and shallow convections, respectively. The stratiform rain is basically organized with convective rain, having similar but slightly lower RTH and slightly lags the convective rain maximum. It is notable that relatively moderate (∼7 km) RTH is dominant in the mature phase of the MJO, while very tall rainfall with RTH over 10 km and lightning frequency increase in the suppressed phase. The rain-yield-per flash (RPF) varies about 20–100% of the mean value of ∼2–10 × 10⁹ kg fl⁻¹ over the tropical warm ocean and that of ∼2–5 × 10⁹ kg fl⁻¹ over the equatorial Islands, between the convectively suppressed phase and the active phase of MJO, in the manner that RPF is smaller in the suppressed phase and larger in the active phase.     

Below-cloud rain scavenging of atmospheric aerosols for aerosol deposition models

Below-cloud aerosol scavenging is generally estimated from field measurements using advanced instruments that measure changes in aerosol distributions with respect to rainfall. In this study, we discuss various scavenging mechanisms and scavenging coefficients from past laboratory and field measurements. Scavenging coefficients derived from field measurements (representing natural aerosols scavenging) are two orders higher than that of theoretical ones for smaller particles (D_p < 2 μm). Measured size-resolved scavenging coefficients can be served as a better option to the default scavenging coefficient (e.g. a constant of 10^?4 s^?1 for all size of aerosols, as used in the CALPUFF model) for representing below-cloud aerosol scavenging. We propose scavenging correction parameter (C_R) as an exponential function of size-resolved scavenging coefficients, winds and width in the downwind of the source–receptor system. For a wind speed of 3 m s^?1, C_R decrease with the width in the downwind for particles of diameters D_p < 0.1 μm but C_R does not vary much for particles in the accumulation mode (0.1 < D_p < 2 μm). For a typical urban aerosol distribution, assuming 3 m s^?1 air-flow in the source–receptor system, 10 km downwind width, 2.84 mm h^?1 of rainfall and using aerosol size dependent scavenging coefficients in the C_R, scavenging of aerosols is found to be 16% in number and 24% in volume of total aerosols. Using the default scavenging coefficient (10^?4 s^?1) in the CALPUFF model, it is found to be 64% in both number and volume of total aerosols.     

Analysis of quantitative precipitation forecasts using the Dynamic State Index

The German Weather Service (DWD) has two non-hydrostatic operational weather prediction models with different spatial resolution and precipitation parametrisations. The coarser COSMO-EU model has a spatial resolution of 7 km, whereas the higher-resolution COSMO-DE model has a gridspace of 2.8 km and explicitly resolves deep convection. To improve the numerical weather prediction (NWP) models it is necessary to understand precipitation processes. A central goal is the statistical evaluation of precipitation forecasts with dynamic parameters. Here, the Dynamic State Index (DSI) is used as a dynamic threshold parameter. The DSI theoretically describes the change of atmospheric flow fields as deviations from a stationary adiabatic solution of the primitive equations (Névir, 2004). For seasonal area means the DSI shows a remarkably high correlation with the precipitation forecasts provided by the COSMO-DE model. This is especially the case for the summer of 2007. The same analysis has been performed with the COSMO-EU forecast data and the results were compared with those from the COSMO-DE model. Moreover, an independent precipitation analysis, with a resolution corresponding to 7 km and 2.8 km, has been compared with respect to modelled precipitation and the DSI. In addition, correlations between the DSI and modelled as well as observed precipitation as a function of the forecast time for the different grid resolutions are also presented. The results show, that after 12 h, the correlation of the persistence forecast with the DSI reaches two thirds of the initial value. Thus, the DSI offers itself as a new dynamic forecast tool for precipitation events.     

Variability of the Somali Current and eddies during the southwest monsoon regimes

The meso-scale eddies and currents in the Arabian Sea are analyzed using different satellite observations, Simple Oceanic Data Assimilation (SODA) reanalysis, and Ocean Reanalysis System 4 (ORAS4) from 1993 to 2016 to investigate the impacts of Southwest (SW) Monsoon strength on Somali Current (SC) mesoscale circulations such as the Great Whirl (GW), the Socotra Eddy (SE), the Southern Gyre (SG), and smaller eddies. Increased Ekman pumping during stronger SW monsoons strengthens coastal upwelling along the Somali coast. The Arabian Sea basin-wide anticyclonic circulation and presence of the GW form mesoscale circulation patterns favourable to advection of upwelled waters eastward into the central Arabian Sea. In September, after the SW monsoon winds reach peak strength in July and August, a higher number of discrete anticyclonic eddies with higher ( > 20 cm) sea surface height anomalies develop in strong and normal intensity SW monsoon seasons than weaker SW monsoon seasons.     

Assessing 4D-VAR for dynamical mapping of coastal high-frequency radar in San Diego

The problem of dynamically mapping high-frequency (HF) radar radial velocity observations is investigated using a three-dimensional hydrodynamic model of the San Diego coastal region and an adjoint-based assimilation method. The HF radar provides near-real-time radial velocities from three sites covering the region offshore of San Diego Bay. The hydrodynamical model is the Massachusetts Institute of Technology general circulation model (MITgcm) with 1 km horizontal resolution and 40 vertical layers. The domain is centered on Point Loma, extending 117 km offshore and 120 km alongshore. The reference run (before adjustment) is initialized from a single profile of T and S and is forced with wind data from a single shore station and with zero heat and fresh water fluxes. The adjoint of the model is used to adjust initial temperature, salinity, and velocity, hourly temperature, salinity and horizontal velocities at the open boundaries, and hourly surface fluxes of momentum, heat and freshwater so that the model reproduces hourly HF radar radial velocity observations. Results from a small number of experiments suggest that the adjoint method can be successfully used over 10-day windows at coastal model resolution. It produces a dynamically consistent model run that fits HF radar data with errors near the specified uncertainties. In a test of the forecasting capability of the San Diego model after adjustment, the forecast skill was shown to exceed persistence for up to 20 h.     

Water on an urban planet: Urbanization and the reach of urban water infrastructure

Urban growth is increasing the demand for freshwater resources, yet surprisingly the water sources of the world's large cities have never been globally assessed, hampering efforts to assess the distribution and causes of urban water stress. We conducted the first global survey of the large cities’ water sources, and show that previous global hydrologic models that ignored urban water infrastructure significantly overestimated urban water stress. Large cities obtain 78 ± 3% of their water from surface sources, some of which are far away: cumulatively, large cities moved 504 billion liters a day (184 km³ yr⁻¹) a distance of 27,000 ± 3800 km, and the upstream contributing area of urban water sources is 41% of the global land surface. Despite this infrastructure, one in four cities, containing $4.8 ± 0.7 trillion in economic activity, remain water stressed due to geographical and financial limitations. The strategic management of these cities’ water sources is therefore important for the future of the global economy.     

A method for determining thermophysical properties of organic material in aqueous solutions: Succinic acid

A method for determining evaporation rates and thermodynamic properties of aqueous solution droplets is introduced. The method combines evaporation rate measurements using modified TDMA technique with data evaluation using an accurate evaporation model. The first set of data has been collected and evaluated for succinic acid aqueous solution droplets.Evaporation rates of succinic acid solution droplets have been measured using a TDMA system at controlled relative humidity (65%) and temperature (298 K). A temperature-dependent expression for the saturation vapour pressure of pure liquid phase succinic acid at atmospheric temperatures has been derived by analysing the evaporation rate data with a numerical model. The obtained saturation vapour pressure of liquid phase succinic acid is ln(p) = 118.41 − 16204.8/T − 12.452ln(T). The vapour pressure is in unit of Pascal and the temperature in Kelvin. A linear expression for the enthalpy of vaporization for liquid state succinic acid is also presented.According to the results presented in the following, a literature expression for the vapour pressure of liquid phase succinic acid defined for temperatures higher than 461 K [Yaws, C.L., 2003. Yaws' Handbook of Thermodynamic and Physical Properties of Chemical Compounds, Knovel] can be extrapolated to atmospheric temperatures with very good accuracy. The results also suggest that at 298 K the mass accommodation coefficient of succinic acid is unity or very close to unity.     

A review of ground-based,mobile, W-band Doppler-radar observations of tornadoes and dust devils

A ground-based, mobile, W-band Doppler-radar has been used in the U.S. during the last decade to obtain high-spatial resolution maps of the radar reflectivity and wind fields in tornadoes and dust devils. This radar is one of the best tools available for studying the substructure of intense, small-scale vortices in the boundary layer. The most significant findings to date are summarized.In one case, it was found that just prior to tornadogenesis in a supercell, a 100–200 m scale cyclonic vortex formed at the leading edge of a bulge in the rear–flank gust front. This vortex appeared to interact with a larger-scale (500 m to 1 km wide) cyclonic vortex, just as the tornado formed. Other small-scale cyclonic vortices were present along the rear–flank gust, but they did not develop into tornadoes. The mature tornado-vortex was dominated by quasi-stationary wavenumber-two disturbances, while the mean vortex resembled a two-celled, Rankine combined vortex. The diameter of the mean vortex narrowed as it intensified and widened as it weakened, even though the tornado condensation funnel narrowed as the tornado was dissipating. Evidence was also found of short-term, inertial-like oscillations in vortex diameter and intensity. Spiral bands and eyes were ubiquitous. The eye in one well-documented case was broader in the lowest few hundred meters than it was aloft. Multiple vortices and “umbilical” cords of very narrow bands of reflectivity have also been found.Both cyclonic and anticyclonic dust devils have been documented. Some dust devils resemble a relatively narrow, Rankine combined vortex, while others are wider and have a broad, calm eye and a narrow annulus of intense vorticity just within the radius of maximum wind (RMW), and rising motion just inside the RMW and sinking motion well inside the RMW. Multiple-vortex structure, Rossby-like wave motion, and the Fujiwhara effect have also been documented.     

A Loop Current experiment: Field and remote measurements

An overview of a new comprehensive observational study of the Loop Current (LC) in the eastern Gulf of Mexico that encompassed full-depth and near-bottom moorings, pressure-equipped inverted echo sounders (PIES) and remote sensing is presented. The study array was designed to encompass the LC from the Campeche Bank to the west Florida escarpment. This overview centers about principal findings as they pertain to mesoscale dynamics. Two companion papers provide in-depth analyses. Three LC anticyclonic eddy separation events were observed with good 3D spatial coverage over the 2½ year extent of the field study; the three separations exhibited similar processes after the LC had extended into the eastern Gulf. Large scale (∼300 km wavelength, 40–60 day periods) southward propagating meanders developed on the eastern side of the LC over deep (∼3000 m) water that were the result of baroclinic instability between the upper layer meandering jet and lower layer cyclones and anticyclones. The lower layer was only highly energetic during relatively short (∼2–3 months) intervals just prior to or during eddy detachments because of baroclinic instability. The steepening of the meanders lead to a pinch-off of LC eddies. The deep lower-layer eddies, constrained by the closed topography of the southeastern Gulf, propagated westward across the detachment zone and appear to assist in achieving separation. Small scale (∼50–100 km, periods ∼10 days) frontal eddies, observed on the western side of the LC along the Campeche Bank slope, decay over the deep water of the northern part of an extended LC, and have little influence on lower layer eddies, the east side meanders and the eddy detachment processes.     

The 183-WSL fast rain rate retrieval algorithm: Part I: Retrieval design

The Water vapour Strong Lines at 183 GHz (183-WSL) fast retrieval method retrieves rain rates and classifies precipitation types for applications in nowcasting and weather monitoring. The retrieval scheme consists of two fast algorithms, over land and over ocean, that use the water vapour absorption lines at 183.31 GHz corresponding to the channels 3 (183.31 ± 1 GHz), 4 (183.31 ± 3 GHz) and 5 (183.31 ± 7 GHz) of the Advanced Microwave Sounding Unit module B (AMSU-B) and of the Microwave Humidity Sounder (MHS) flying on NOAA-15-18 and Metop-A satellite series, respectively.The method retrieves rain rates by exploiting the extinction of radiation due to rain drops following four subsequent steps. After ingesting the satellite data stream, the window channels at 89 and 150 GHz are used to compute scattering-based thresholds and the 183-WSLW module for rainfall area discrimination and precipitation type classification as stratiform or convective on the basis of the thresholds calculated for land/mixed and sea surfaces. The thresholds are based on the brightness temperature difference Δ_win = T_B89 ? T_B150 and are different over land (L) and over sea (S): cloud droplets and water vapour (Δ_win < 3 K L; Δ_win < 0 K S), stratiform rain (3 K < Δ_win < 10 K L; 0 K < Δ_win < 10 K S), and convective rain (Δ_win > 10 K L and S). The thresholds, initially empirically derived from observations, are corroborated by the simulations of the RTTOV radiative transfer model applied to 20000 ECMWF atmospheric profiles at midlatitudes and the use of data from the Nimrod radar network. A snow cover mask and a digital elevation model are used to eliminate false rain area attribution, especially over elevated terrain. A probability of detection logistic function is also applied in the transition region from no-rain to rain adjacent to the clouds to ensure continuity of the rainfall field. Finally, the last step is dedicated to the rain rate retrieval with the modules 183-WSLS (stratiform) and 183WSLC (convective), and the module 183-WSL for total rainfall intensity derivation.A comparison with rainfall retrievals from the Goddard Profiling (GPROF) TRMM 2A12 algorithm is done with good results on a stratiform and hurricane case studies. A comparison is also conducted with the MSG-based Precipitation Index (PI) and the Scattering Index (SI) for a convective-stratiform event showing good agreement with the 183-WSLC retrieval. A complete validation of the product is the subject of Part II of the paper.     

Climate change impacts on ecosystems and the terrestrial carbon sink: a new assessment

Climate output from the UK Hadley Centre's HadCM2 and HadCM3 experiments for the period 1860 to 2100, with IS92a greenhouse gas forcing, together with predicted patterns of N deposition and increasing CO₂, were input (offline) to the dynamic vegetation model, Hybrid v4.1 (Friend et al., 1997; Friend and White, 1999). This model represents biogeochemical, biophysical and biogeographical processes, coupling the carbon, nitrogen and water cycles on a sub-daily timestep, simulating potential vegetation and transient changes in annual growth and competition between eight generalized plant types in response to climate.Global vegetation carbon was predicted to rise from about 600 to 800 PgC (or to 650 PgC for HadCM3) while the soil carbon pool of about 1100 PgC decreased by about 8%. By the 2080s, climate change caused a partial loss of Amazonian rainforest, C₄ grasslands and temperate forest in areas of southern Europe and eastern USA, but an expansion in the boreal forest area. These changes were accompanied by a decrease in net primary productivity (NPP) of vegetation in many tropical areas, southern Europe and eastern USA (in response to warming and a decrease in rainfall), but an increase in NPP of boreal forests. Global NPP increased from 45 to 50 PgC y⁻¹ in the 1990s to about 65 PgC y⁻¹ in the 2080s (about 58 PgC y⁻¹ for HadCM3). Global net ecosystem productivity (NEP) increased from about 1.3 PgC y⁻¹ in the 1990s to about 3.6 PgC y⁻¹ in the 2030s and then declined to zero by 2100 owing to a loss of carbon from declining forests in the tropics and at warm temperate latitudes — despite strengthening of the carbon sink at northern high latitudes. HadCM3 gave a more erratic temporal evolution of NEP than HadCM2, with a dramatic collapse in NEP in the 2050s.     

Pic 2005, a field campaign to investigate low-tropospheric ozone variability in the Pyrenees

The Pic 2005 field campaign took place from 13 June to 7 July 2005 close to the high-altitude permanent atmospheric observatory Pic-du-Midi (PDM), situated at 2875 m asl in the French Pyrenees. The experimental set-up combined in situ ground-based observations at PDM with ozone lidar measurements at two lower sites in close vicinity (600 m asl/28 km away, and 2380 m asl/500 m away). Such an experimental configuration is appropriate to address the question of the vertical layering of the chemical atmosphere in a mountain area and above the plain nearby, and how this influences measurements conducted on a mountain summit under the influence of horizontal transport at regional scale, and vertical transport at local scale. Forecast tools made it possible to plan and carry out 6 one-day Intensive Observation Periods (IOPs), mostly in anticyclonic conditions favoring local thermally induced circulations, with and without local pollution in the lower troposphere.It was thus possible to document i) ozone diurnal variations at PDM; ii) correlation between ozone measurements at PDM and their counterparts at the same altitude in the free troposphere; iii) ozone variability in the vicinity of PDM.The field campaign provided direct experimental evidence that at daytime in the encountered conditions (mostly anticyclonic), PDM failed in a large extent to be representative of the troposphere above the surrounding flat areas at similar altitude. First, ozone daily averages at PDM were found lower than their free-tropospheric counterpart. Thermally induced circulations and convection pumping clean air from the rural boundary layer can account qualitatively for ozone depletion observed at PDM during daytime. However the surface measurements do not support the hypothesis of direct lifting of near-surface air masses up to PDM. Thus, mixing with free-tropospheric air, photochemistry and surface deposition in the valleys appear to be needed ingredients to account quantitatively for the observed variations (in proportions that further studies should determine). Second, ozone variability was found to be much lower at PDM than in the free troposphere—again an indication of atmospheric mixing. In particular at daytime, the PDM observatory did not allow for detection of ozone-rich layers simultaneously visible above the plain. Beyond these first results, the data set presented here paves way to detailed studies of the IOPs.     

The possible role of Brazilian promontory in Little Ice Age

The Gulf Stream, one of the strongest currents in the world, transports approximately 31 Sv of water (Kelly and Gille, 1990, Baringer and Larsen, 2001, Leaman et al., 1995) and 1.3 × 10¹⁵ W (Larsen, 1992) of heat into the Atlantic Ocean, and warms the vast European continent. Thus any change of the Gulf Stream could lead to the climate change in the European continent, and even worldwide (Bryden et al., 2005). Past studies have revealed a diminished Gulf Stream and oceanic heat transport that was possibly associated with a southward migration of intertropical convergence zone (ITCZ) and may have contributed to Little Ice Age (AD ∼1200 to 1850) in the North Atlantic (Lund et al., 2006). However, the causations of the Gulf Stream weakening due to the southward migration of the ITCZ remain uncertain. Here we use satellite observation data and employ a model (oceanic general circulation model – OGCM) to demonstrate that the Brazilian promontory in the east coast of South America may have played a crucial role in allocating the equatorial currents, while the mean position of the equatorial currents migrates between northern and southern hemisphere in the Atlantic Ocean. Northward migrations of the equatorial currents in the Atlantic Ocean have little influence on the Gulf Stream. Nevertheless, southward migrations, especially abrupt large southward migrations of the equatorial currents, can lead to the increase of the Brazil Current and the significant decrease of the North Brazil Current, in turn the weakening of the Gulf Stream. The results from the model simulations suggest the mean position of the equatorial currents in the Atlantic Ocean shifted at least 180–260 km southwards of its present-day position during the Little Ice Age based on the calculations of simple linear equations and the OGCM simulations.     

Measurements and simulation of the M-component current and simultaneous electromagnetic fields at 60 m and 550 m

Simultaneous measurements of the M-component current (surges superimposed on lightning continuing currents) and the corresponding electromagnetic fields at 60 m and 550 m from the lightning channel are analyzed and simulated with a two-wave model. The measured results reveal that the M-component current at the bottom of the channel exhibits a V-shape character with a leading edge of 78 μs and a trailing edge of 194 μs, while the electric field pulses at 60 m and 550 m have trailing edges faster than leading edges. The peak of the M-component current lags behind the electric field peak by tens of microseconds, when the distance increases to 550 m, the disparity of the time shift increases as well. However, the waveshape of the M-component current is similar to that of the magnetic field pulse. The M-component electric fields at 60 m and 550 m are 1.16 kV/m and 0.17 kV/m, respectively, and exhibit a logarithmic distance dependence which implies that the M-component charge density increases with height. Additionally, a two-wave model is used to examine the sensitivity of the predicted electric and magnetic fields to the speed and current reflection coefficient variations of the M-component. The simulated results show that the effects are different for the electric and magnetic fields. The M-component speed essentially controls the electric field, but has little effect on the magnetic field. Larger reflection coefficient results in a larger magnetic field, but a smaller electric field.     

Kinematic metrics of the Loop Current in the Gulf of Mexico from satellite altimetry

We analyzed a 20-year time series (January 1st, 1993 through December 31st, 2012) of Loop Current (LC) surface area derived from satellite altimetry in the eastern Gulf of Mexico to estimate kinematical metrics of this potent flow. On average the LC intrudes to its maximum northward position about 216 ± 126 days after the previous eddy separation; and ∼30 ± 31 days later sheds a large anticyclonic eddy. When the northern extent of the LC intrusion following the previous eddy separation is greater than 27°N, the current retreats very quickly until it sheds another eddy with the entire separation process occurring on the order of 30 days. To first order the change in areal extent of the LC during intrusion into the Gulf occurs at an average rate of 225 km² day⁻¹, which corresponds to an intrusion velocity of 1.7 cm s⁻¹ of the LC front, and adds Caribbean water to the Gulf at a rate of 2.6 ± 0.7 Sv.     

Concentrations and sources of PAHs at three stations in Istanbul,Turkey

The chemical mass balance model was applied to atmospheric Polycyclic Aromatic Hydrocarbons (PAHs) in Istanbul, Turkey. A total of 326 airborne samples were collected and analyzed for 16 PAHs and Total Suspended Particles (TSP) in the September 2006–December 2007 period at three monitoring stations: Yildiz, DMO (urban sites) and Kilyos (rural site). The total average PAH concentrations were 100.66 ± 61.26, 84.63 ± 46.66 and 25.12 ± 13.34 ng m^?3 and the TSP concentrations were 101.16 ± 53.22, 152.31 ± 99.12, 49.84 ± 18.58 μg m^?3 for Yildiz, DMO and Kilyos stations respectively. At all the sites, the lighter compounds were the most abundant, notably Nap, AcPy and PA. The average correlation values between TSP and total heavier PAH were greater than 0.5 for Yildiz and DMO stations. The patterns of PAH and TSP concentrations showed spatial and temporal variations. PAH concentrations were evaluated for the PAH contribution from four sources (diesel engines, gasoline engines, natural gas combustion, and coal + wood burning). Vehicle emissions appear to be the major source with contributions of 61.2%, 63.3% and 54.1% for Yildiz, DMO and Kilyos stations respectively. Seasonal and yearly variations had different trends for all sites.     

Buoyancy in tropical cyclones and other rapidly rotating atmospheric vortices

Motivated primarily by its application to understanding tropical-cyclone intensification and maintenance, we re-examine the concept of buoyancy in rapidly rotating vortices, distinguishing between the buoyancy of the symmetric balanced vortex or system buoyancy, and the local buoyancy associated with cloud dynamics. The conventional definition of buoyancy is contrasted with a generalized form applicable to a vortex, which has a radial as well as a vertical component. If, for the special case of axisymmetric motions, the balanced density and pressure distribution of a rapidly rotating vortex are used as the reference state, the buoyancy field then characterizes the unbalanced density perturbations, i.e. the local buoyancy. We show how to determine such a reference state without approximation.The generation of the toroidal circulation of a vortex, which is necessary for vortex amplification, is characterized in the vorticity equation by the baroclinicity vector. This vector depends, inter-alia, on the horizontal (or radial) gradient of buoyancy evaluated along isobaric surfaces. We show that for a tropical-cyclone-scale vortex, the buoyancy so calculated is significantly different from that calculated at constant height or on surfaces of constant σ (σ = (p − p^*)/(p_s − p^*), where p is the actual pressure, p^* some reference pressure and p_s is the surface pressure). Since many tropical-cyclone models are formulated using σ-coordinates, we examine the calculation of buoyancy on σ-surfaces and derive an expression for the baroclinicity vector in σ-coordinates. The baroclinic forcing term in the azimuthal vorticity equation for an axisymmetric vortex is shown to be approximately equal to the azimuthal component of the curl of the generalized buoyancy. A scale analysis indicates that the vertical gradient of the radial component of generalized buoyancy makes a comparatively small contribution to the generation of toroidal vorticity in a tropical cyclone, but may be important in tornadoes and possibly also in dust devils.We derive also a form of the Sawyer–Eliassen equation from which the toroidal (or secondary) circulation of a balanced vortex may be determined. The equation is shown to be the time derivative of the toroidal vorticity equation in which the time rate-of-change of the material derivative of potential toroidal vorticity is set to zero. In analogy with the general case, the diabatic forcing term in the Sawyer–Eliassen equation is shown to be approximately equal to the time rate-of-change of the azimuthal component of the curl of generalized buoyancy.Finally, we discuss the generation of buoyancy in tropical cyclones and contrast the definitions of buoyancy that have been used in recent studies of tropical cyclones. We emphasize the non-uniqueness of the buoyancy force, which depends on the choice of a reference density and pressure, and note that different, but equivalent interpretations of the flow dynamics may be expected to arise if different reference quantities are chosen.     

Distribution of inorganic nitrogen-containing species in atmospheric particles from an island in the Yellow Sea

In this study, total suspended particles (TSP) and size-segregated atmospheric aerosol samples were measured on Qianliyan Island in the Yellow Sea in spring (April–May), summer (July–August) and fall (October–November) of 2006 and in water (January–February) of 2007. The mass concentration of the TSP varied from 75.6 to 132.0 μg/m³. The average concentration were 9.37 ± 7.56 μg/m³ and 5.32 ± 4.25 μg/m³ for nitrate and ammonium in the TSP, respectively. TSP concentration showed a significant correlation with those of nitrate (n = 27, r = 0.73) and ammonium (n = 27, r = 0.60). The mass-size distribution of atmospheric particles exhibited two modes with an accumulation mode at 0.43–1.1 μm and a coarse mode at 3.3–4.7 μm throughout the sampling months. A bi-modal size distribution of nitrate in concentration occurred in the April–May, October–November and January–February, but a uni-modal size distribution occurred in the August. The uni-modal size distribution of ammonium at 0.43–0.65 μm was observed throughout the sampling months. The average of inorganic nitrogen in mass concentration accounted for 4.0% of the total mass of aerosol particles while ammonium-N was the dominant fraction of TIN (Total Inorganic Nitrogen), contributing to 62–71% of the TIN.