Turbulent structure and scaling of the inertial subrange in a stratocumulus-topped boundary layer observed by a Doppler lidar

The turbulent structure of a stratocumulus-topped marine boundary layer over a 2-day period is observed with a Doppler lidar at Mace Head in Ireland. Using profiles of vertical velocity statistics, the bulk of the mixing is identified as cloud driven. This is supported by the pertinent feature of negative vertical velocity skewness in the sub-cloud layer which extends, on occasion, almost to the surface. Both coupled and decoupled turbulence characteristics are observed. The length and timescales related to the cloud-driven mixing are investigated and shown to provide additional information about the structure and the source of the mixing inside the boundary layer. They are also shown to place constraints on the length of the sampling periods used to derive products, such as the turbulent dissipation rate, from lidar measurements. For this, the maximum wavelengths that belong to the inertial subrange are studied through spectral analysis of the vertical velocity. The maximum wavelength of the inertial subrange in the cloud-driven layer scales relatively well with the corresponding layer depth during pronounced decoupled structure identified from the vertical velocity skewness. However, on many occasions, combining the analysis of the inertial subrange and vertical velocity statistics suggests higher decoupling height than expected from the skewness profiles. Our results show that investigation of the length scales related to the inertial subrange significantly complements the analysis of the vertical velocity statistics and enables a more confident interpretation of complex boundary layer structures using measurements from a Doppler lidar.

Clouds, aerosols, and precipitation in the marine boundary layer: an ARM Mobile Facility Deployment

The Clouds, Aerosol, and Precipitation in the Marine Boundary Layer (CAP-MBL) deployment at Graciosa Island in the Azores generated a 21-month (April 2009–December 2010) comprehensive dataset documenting clouds, aerosols, and precipitation using the Atmospheric Radiation Measurement Program (ARM) Mobile Facility (AMF). The scientific aim of the deployment is to gain improved understanding of the interactions of clouds, aerosols, and precipitation in the marine boundary layer. Graciosa Island straddles the boundary between the subtropics and midlatitudes in the northeast Atlantic Ocean and consequently experiences a great diversity of meteorological and cloudiness conditions. Low clouds are the dominant cloud type, with stratocumulus and cumulus occurring regularly. Approximately half of all clouds contained precipitation detectable as radar echoes below the cloud base. Radar and satellite observations show that clouds with tops from 1 to 11 km contribute more or less equally to surface-measured precipitation at Graciosa. A wide range of aerosol conditions was sampled during the deployment consistent with the diversity of sources as indicated by back-trajectory analysis. Preliminary findings suggest important two-way interactions between aerosols and clouds at Graciosa, with aerosols affecting light precipitation and cloud radiative properties while being controlled in part by precipitation scavenging. The data from Graciosa are being compared with short-range forecasts made with a variety of models. A pilot analysis with two climate and two weather forecast models shows that they reproduce the observed time-varying vertical structure of lower-tropospheric cloud fairly well but the cloud-nucleating aerosol concentrations less well. The Graciosa site has been chosen to be a permanent fixed ARM site that became operational in October 2013.

Advances in understanding mineral dust and boundary layer processes over the Sahara from Fennec aircraft observations

The Fennec climate program aims to improve understanding of the Saharan climate system through a synergy of observations and modelling. We present a description of the Fennec airborne observations during 2011 and 2012 over the remote Sahara (Mauritania and Mali) and the advances in the understanding of mineral dust and boundary layer processes they have provided. Aircraft instrumentation aboard the UK FAAM BAe146 and French SAFIRE Falcon 20 is described, with specific focus on instrumentation specially developed and relevant to Saharan meteorology and dust. Flight locations, aims and associated meteorology are described. Examples and applications of aircraft measurements from the Fennec flights are presented, highlighting new scientific results delivered using a synergy of different instruments and aircraft. These include: (1) the first airborne measurement of dust particles sized up to 300 microns and associated dust fluxes in the Saharan atmospheric boundary layer (SABL), (2) dust uplift from the breakdown of the nocturnal low-level jet before becoming visible in SEVIRI satellite imagery, (3) vertical profiles of the unique vertical structure of turbulent fluxes in the SABL, (4) in-situ observations of processes in SABL clouds showing dust acting as CCN and IN at −15 °C, (5) dual-aircraft observations of the SABL dynamics, thermodynamics and composition in the Saharan heat low region (SHL), (6) airborne observations of a dust storm associated with a cold-pool (haboob) issued from deep convection over the Atlas, (7) the first airborne chemical composition measurements of dust in the SHL region with differing composition, sources (determined using Lagrangian backward trajectory calculations) and absorption properties between 2011 and 2012, (8) coincident ozone and dust surface area measurements suggest coarser particles provide a route for ozone depletion, (9) discrepancies between airborne coarse mode size distributions and AERONET sunphotometer retrievals under light dust loadings. These results provide insights into boundary layer and dust processes in the SHL region – a region of substantial global climatic importance.

Joint retrievals of cloud and drizzle in marine boundary layer clouds using ground-based radar, lidar and zenith radiances

Active remote sensing of marine boundary-layer clouds is challenging as drizzle drops often dominate the observed radar reflectivity. We present a new method to simultaneously retrieve cloud and drizzle vertical profiles in drizzling boundary-layer clouds using surface-based observations of radar reflectivity, lidar attenuated backscatter, and zenith radiances under conditions when precipitation does not reach the surface. Specifically, the vertical structure of droplet size and water content of both cloud and drizzle is characterised throughout the cloud. An ensemble optimal estimation approach provides full error statistics given the uncertainty in the observations. To evaluate the new method, we first perform retrievals using synthetic measurements from large-eddy simulation snapshots of cumulus under stratocumulus, where cloud water path is retrieved with an error of 31 g m−2 . The method also performs well in non-drizzling clouds where no assumption of the cloud profile is required. We then apply the method to observations of marine stratocumulus obtained during the Atmospheric Radiation Measurement MAGIC deployment in the Northeast Pacific. Here, retrieved cloud water path agrees well with independent three-channel microwave radiometer retrievals, with a root mean square difference of 10–20 g m−2.

Improving a global model from the boundary layer: total turbulent energy and the neutral limit Prandtl number

Model intercomparisons have identified important deficits in the representation of the stable boundary layer by turbulence parametrizations used in current weather and climate models. However, detrimental impacts of more realistic schemes on the large-scale flow have hindered progress in this area. Here we implement a total turbulent energy scheme into the climate model ECHAM6. The total turbulent energy scheme considers the effects of Earth’s rotation and static stability on the turbulence length scale. In contrast to the previously used turbulence scheme, the TTE scheme also implicitly represents entrainment flux in a dry convective boundary layer. Reducing the previously exaggerated surface drag in stable boundary layers indeed causes an increase in southern hemispheric zonal winds and large-scale pressure gradients beyond observed values. These biases can be largely removed by increasing the parametrized orographic drag. Reducing the neutral limit turbulent Prandtl number warms and moistens low-latitude boundary layers and acts to reduce longstanding radiation biases in the stratocumulus regions, the Southern Ocean and the equatorial cold tongue that are common to many climate models.

Impact of Laptev Sea flaw polynyas on the atmospheric boundary layer and ice production using idealized mesoscale simulations

The interaction between polynyas and the atmospheric boundary layer is examined in the Laptev Sea using the regional, non-hydrostatic Consortium for Small-scale Modelling (COSMO) atmosphere model. A thermodynamic sea-ice model is used to consider the response of sea-ice surface temperature to idealized atmospheric forcing. The idealized regimes represent atmospheric conditions that are typical for the Laptev Sea region. Cold wintertime conditions are investigated with sea-ice–ocean temperature differences of up to 40 K. The Laptev Sea flaw polynyas strongly modify the atmospheric boundary layer. Convectively mixed layers reach heights of up to 1200 m above the polynyas with temperature anomalies of more than 5 K. Horizontal transport of heat expands to areas more than 500 km downstream of the polynyas. Strong wind regimes lead to a more shallow mixed layer with strong near-surface modifications, while weaker wind regimes show a deeper, well-mixed convective boundary layer. Shallow mesoscale circulations occur in the vicinity of ice-free and thin-ice covered polynyas. They are forced by large turbulent and radiative heat fluxes from the surface of up to 789 W m−2, strong low-level thermally induced convergence and cold air flow from the orographic structure of the Taimyr Peninsula in the western Laptev Sea region. Based on the surface energy balance we derive potential sea-ice production rates between 8 and 25 cm d−1. These production rates are mainly determined by whether the polynyas are ice-free or covered by thin ice and by the wind strength.

The BALTIMOS (BALTEX Integrated Model System) fi eld experiments: A comprehensive atmospheric boundary layer data set for model validation over the open and ice- covered Baltic Sea

A comprehensive atmospheric boundary layer (ABL) data set was collected in eight fi eld experiments (two during each season) over open water and sea ice in the Baltic Sea during 1998–2001 with the primary objective to validate the coupled atmospheric- ice-ocean-land surface model BALTIMOS (BALTEX Integrated Model System). Measurements were taken by aircraft, ships and surface stations and cover the mean and turbulent structure of the ABL including turbulent fl uxes, radiation fl uxes, and cloud conditions. Measurement examples of the spatial variability of the ABL over the ice edge zone and of the stable ABL over open water demonstrate the wide range of ABL conditions collected and the strength of the data set which can also be used to validate other regional models.

Seasonality of submesoscale flows in the ocean surface boundary layer

A signature of submesoscale flows in the upper ocean is skewness in the distribution of relative vorticity. Expected to result for high Rossby-number flows, such skewness has implications for mixing, dissipation and stratification within the upper ocean. An array of moorings deployed in the Northeast Atlantic for one year as part of the OSMOSIS experiment reveals that relative vorticity is positively skewed during winter even though the scale of the Rossby number is less than 0.5. Furthermore, this skewness is reduced to zero during spring and autumn. There is also evidence of modest seasonal variations in the gradient Rossby number. The proposed mechanism by which relative vorticity is skewed is that the ratio of lateral to vertical buoyancy gradients, as summarized by the inverse gradient Richardson number, restricts its range during winter but less so at other times of the year. These results support recent observations and model simulations suggesting the upper ocean is host to a seasonal cycle in submesoscale turbulence.

Numerical simulation of viscous three-dimensional flow over aerospace vehicles using Euler/boundary-layer zonal approach

This paper presents a viscous three-dimensional simulations coupling Euler and boundary layer codes for calculating flows over arbitrary surfaces. The governing equations are written in a general non orthogonal coordinate system. The Levy-Lees transformation generalized to three-dimensional flows is utilized. The inviscid properties are obtained from the Euler equations using the Beam and Warming implicit approximate factorization scheme. The resulting equations are discretized and approximated by a two-point fmitedifference numerical scheme. The code developed is validated and applied to the simulation of the flowfield over aerospace vehicle configurations. The results present good correlation with the available data.

Suppression of vortex-induced vibration using moving surface boundary-layer control

Experimental results of flow around a circular cylinder with moving surface boundary-layer control (MSBC) are presented. Two small rotating cylinders strategically located inject momentum in the boundary layer of the cylinder, which delays the separation of the boundary layer. As a consequence, the wake becomes narrower and the fluctuating transverse velocity is reduced, resulting in a recirculation free region that prevents the vortex formation. The control parameter is the ratio between the tangential velocity of the moving surface and the flow velocity (U-c/U). The main advantage of the MSBC is the possibility of combining the suppression of vortex-induced vibration (VIV) and drag reduction. The experimental tests are preformed at a circulating water channel facility and the circular cylinders are mounted on a low-damping air bearing base with one degree-of-freedom in the transverse direction of the channel flow. The mass ratio is 1.8. The Reynolds number ranges from 1600 to 7500, the reduced velocity varies up to 17, and the control parameter interval is U-c/U = 5-10. A significant decreasing in the maximum amplitude of oscillation for the cylinder with MSBC is observed. Drag measurements are obtained for statically mounted cylinders with and without MSBC. The use of the flow control results in a mean drag reduction at U-c/U = 5 of almost 60% compared to the plain cylinder. PIV velocity fields of the wake of static cylinders are measured at Re = 3000. The results show that the wake is highly organized and narrower compared to the one observed in cylinders without control. The calculation of the total variance of the fluctuating transverse velocity in the wake region allows the introduction of an active closed-loop control. The experimental results are in good agreement with the numerical simulation studies conducted by other researchers for cylinders with MSBC. (C) 2012 Elsevier Ltd. All rights reserved.

Mesoscale dynamics and boundary-layer structure in topographically forced low-level jets

Two types of mesoscale wind-speed jet and their effects on boundary-layer structure were studied. The first is a coastal jet off the northern California coast, and the second is a katabatic jet over Vatnajökull, Iceland. Coastal regions are highly populated, and studies of coastal meteorology are of general interest for environmental protection, fishing industry, and for air and sea transportation. Not so many people live in direct contact with glaciers but properties of katabatic flows are important for understanding glacier response to climatic changes. Hence, the two jets can potentially influence a vast number of people. Flow response to terrain forcing, transient behavior in time and space, and adherence to simplified theoretical models were examined. The turbulence structure in these stably stratified boundary layers was also investigated. Numerical modeling is the main tool in this thesis; observations are used primarily to ensure a realistic model behavior. Simple shallow-water theory provides a useful framework for analyzing high-velocity flows along mountainous coastlines, but for an unexpected reason. Waves are trapped in the inversion by the curvature of the wind-speed profile, rather than by an infinite stability in the inversion separating two neutral layers, as assumed in the theory. In the absence of blocking terrain, observations of steady-state supercritical flows are not likely, due to the diurnal variation of flow criticality. In many simplified models, non-local processes are neglected. In the flows studied here, we showed that this is not always a valid approximation. Discrepancies between simulated katabatic flow and that predicted by an analytical model are hypothesized to be due to non-local effects, such as surface inhomogeneity and slope geometry, neglected in the theory. On a different scale, a reason for variations in the shape of local similarity scaling functions between studies is suggested to be differences in non-local contributions to the velocity variance budgets.

Modeling the gas and aqueous phase chemistry of the marine boundary layer

Ein eindimensionales numerisches Modell der maritimenGrenzschicht (MBL) wurde erweitert, um chemische Reaktionenin der Gasphase, von Aerosolpartikeln und Wolkentropfen zu beschreiben. Ein Schwerpunkt war dabei die Betrachtung derReaktionszyklen von Halogenen. Soweit Ergebnisse vonMesskampagnen zur Verfuegung standen, wurden diese zurValidierung des Modells benutzt. Die Ergebnisse von frueheren Boxmodellstudien konntenbestaetigt werden. Diese zeigten die saeurekatalysierteAktivierung von Brom aus Seesalzaerosolen, die Bedeutung vonHalogenradikalen fuer die Zerstoerung von O3, diepotentielle Rolle von BrO bei der Oxidation von DMS und dievon HOBr und HOCl in der Oxidation von S(IV). Es wurde gezeigt, dass die Beruecksichtigung derVertikalprofile von meteorologischen und chemischen Groessenvon grosser Bedeutung ist. Dies spiegelt sich darin wider,dass Maxima des Saeuregehaltes von Seesalzaerosolen und vonreaktiven Halogenen am Oberrand der MBL gefunden wurden.Darueber hinaus wurde die Bedeutung von Sulfataerosolen beidem aktiven Recyceln von weniger aktiven zu photolysierbarenBromspezies gezeigt. Wolken haben grosse Auswirkungen auf die Evolution und denTagesgang der Halogene. Dies ist nicht auf Wolkenschichtenbeschraenkt. Der Tagesgang der meisten Halogene ist aufgrundeiner erhoehten Aufnahme der chemischen Substanzen in die Fluessigphase veraendert. Diese Ergebnisse betonen dieWichtigkeit der genauen Dokumentation der meteorologischenBedingungen bei Messkampagnen (besonders Wolkenbedeckungsgrad und Fluessigwassergehalt), um dieErgebnisse richtig interpretieren und mit Modellresultatenvergleichen zu koennen. Dieses eindimensionale Modell wurde zusammen mit einemBoxmodell der MBL verwendet, um die Auswirkungen vonSchiffemissionen auf die MBL abzuschaetzen, wobei dieVerduennung der Abgasfahne parameterisiert wurde. DieAuswirkungen der Emissionen sind am staerksten, wenn sie insauberen Gebieten stattfinden, die Hoehe der MBL gering istund das Einmischen von Hintergrundluft schwach ist.Chemische Reaktionen auf Hintergrundaerosolen spielen nureine geringe Rolle. In Ozeangebieten mit schwachemSchiffsverkehr sind die Auswirkungen auf die Chemie der MBL beschraenkt. In staerker befahrenen Gebieten ueberlappensich die Abgasfahnen mehrerer Schiffe und sorgen fuerdeutliche Auswirkungen. Diese Abschaetzung wurde mitSimulationen verglichen, bei denen die Emissionen alskontinuierliche Quellen behandelt wurden, wie das inglobalen Chemiemodellen der Fall ist. Wenn die Entwicklungder Abgasfahne beruecksichtigt wird, sind die Auswirkungendeutlich geringer da die Lebenszeit der Abgase in der erstenPhase nach Emission deutlich reduziert ist.

Detection and Modeling of Boundary Layer Height

This thesis tackles the problem of the automated detection of the atmospheric boundary layer (BL) height, h, from aerosol lidar/ceilometer observations. A new method, the Bayesian Selective Method (BSM), is presented. It implements a Bayesian statistical inference procedure which combines in an statistically optimal way different sources of information. Firstly atmospheric stratification boundaries are located from discontinuities in the ceilometer back-scattered signal. The BSM then identifies the discontinuity edge that has the highest probability to effectively mark the BL height. Information from the contemporaneus physical boundary layer model simulations and a climatological dataset of BL height evolution are combined in the assimilation framework to assist this choice. The BSM algorithm has been tested for four months of continuous ceilometer measurements collected during the BASE:ALFA project and is shown to realistically diagnose the BL depth evolution in many different weather conditions. Then the BASE:ALFA dataset is used to investigate the boundary layer structure in stable conditions. Functions from the Obukhov similarity theory are used as regression curves to fit observed velocity and temperature profiles in the lower half of the stable boundary layer. Surface fluxes of heat and momentum are best-fitting parameters in this exercise and are compared with what measured by a sonic anemometer. The comparison shows remarkable discrepancies, more evident in cases for which the bulk Richardson number turns out to be quite large. This analysis supports earlier results, that surface turbulent fluxes are not the appropriate scaling parameters for profiles of mean quantities in very stable conditions. One of the practical consequences is that boundary layer height diagnostic formulations which mainly rely on surface fluxes are in disagreement to what obtained by inspecting co-located radiosounding profiles.

Modelling the stable, nocturnal boundary layer: different approaches to turbulence closure problem

Il lavoro è dedicato all'analisi fisica e alla modellizzazione dello strato limite atmosferico in condizioni stabili. L'obiettivo principale è quello di migliorare i modelli di parametrizzazione della turbulenza attualmente utilizzati dai modelli meteorologici a grande scala. Questi modelli di parametrizzazione della turbolenza consistono nell' esprimere gli stress di Reynolds come funzioni dei campi medi (componenti orizzontali della velocità e temperatura potenziale) usando delle chiusure. La maggior parte delle chiusure sono state sviluppate per i casi quasi-neutrali, e la difficoltà è trattare l'effetto della stabilità in modo rigoroso. Studieremo in dettaglio due differenti modelli di chiusura della turbolenza per lo strato limite stabile basati su assunzioni diverse: uno schema TKE-l (Mellor-Yamada,1982), che è usato nel modello di previsione BOLAM (Bologna Limited Area Model), e uno schema sviluppato recentemente da Mauritsen et al. (2007). Le assunzioni delle chiusure dei due schemi sono analizzate con dati sperimentali provenienti dalla torre di Cabauw in Olanda e dal sito CIBA in Spagna. Questi schemi di parametrizzazione della turbolenza sono quindi inseriti all'interno di un modello colonnare dello strato limite atmosferico, per testare le loro predizioni senza influenze esterne. Il confronto tra i differenti schemi è effettuato su un caso ben documentato in letteratura, il "GABLS1". Per confermare la validità delle predizioni, un dataset tridimensionale è creato simulando lo stesso caso GABLS1 con una Large Eddy Simulation. ARPS (Advanced Regional Prediction System) è stato usato per questo scopo. La stratificazione stabile vincola il passo di griglia, poichè la LES deve essere ad una risoluzione abbastanza elevata affinchè le tipiche scale verticali di moto siano correttamente risolte. Il confronto di questo dataset tridimensionale con le predizioni degli schemi turbolenti permettono di proporre un insieme di nuove chiusure atte a migliorare il modello di turbolenza di BOLAM. Il lavoro è stato compiuto all' ISAC-CNR di Bologna e al LEGI di Grenoble.

Iodine speciation in atmospheric aerosols in the marine boundary layer

Iodine chemistry plays an important role in the tropospheric ozone depletion and the new particle formation in the Marine Boundary Layer (MBL). The sources, reaction pathways, and the sinks of iodine are investigated using lab experiments and field observations. The aims of this work are, firstly, to develop analytical methods for iodine measurements of marine aerosol samples especially for iodine speciation in the soluble iodine; secondly, to apply the analytical methods in field collected aerosol samples, and to estimate the characteristics of aerosol iodine in the MBL. Inductively Coupled Plasma – Mass Spectrometry (ICP-MS) was the technique used for iodine measurements. Offline methods using water extraction and Tetra-methyl-ammonium-hydroxide (TMAH) extraction were applied to measure total soluble iodine (TSI) and total insoluble iodine (TII) in the marine aerosol samples. External standard calibration and isotope dilution analysis (IDA) were both conducted for iodine quantification and the limits of detection (LODs) were both 0.1 μg L-1 for TSI and TII measurements. Online couplings of Ion Chromatography (IC)-ICP-MS and Gel electrophoresis (GE)-ICP-MS were both developed for soluble iodine speciation. Anion exchange columns were adopted for IC-ICP-MS systems. Iodide, iodate, and unknown signal(s) were observed in these methods. Iodide and iodate were separated successfully and the LODs were 0.1 and 0.5 μg L-1, respectively. Unknown signals were soluble organic iodine species (SOI) and quantified by the calibration curve of iodide, but not clearly identified and quantified yet. These analytical methods were all applied to the iodine measurements of marine aerosol samples from the worldwide filed campaigns. The TSI and TII concentrations (medians) in PM2.5 were found to be 240.87 pmol m-3 and 105.37 pmol m-3 at Mace Head, west coast of Ireland, as well as 119.10 pmol m-3 and 97.88 pmol m-3 in the cruise campaign over the North Atlantic Ocean, during June – July 2006. Inorganic iodine, namely iodide and iodate, was the minor iodine fraction in both campaigns, accounting for 7.3% (median) and 5.8% (median) in PM2.5 iodine at Mace Head and over the North Atlantic Ocean, respectively. Iodide concentrations were higher than iodate in most of the samples. In the contrast, more than 90% of TSI was SOI and the SOI concentration was correlated significantly with the iodide concentration. The correlation coefficients (R2) were both higher than 0.5 at Mace Head and in the first leg of the cruise. Size fractionated aerosol samples collected by 5 stage Berner impactor cascade sampler showed similar proportions of inorganic and organic iodine. Significant correlations were obtained in the particle size ranges of 0.25 – 0.71 μm and 0.71 – 2.0 μm between SOI and iodide, and better correlations were found in sunny days. TSI and iodide existed mainly in fine particle size range (< 2.0 μm) and iodate resided in coarse range (2.0 – 10 μm). Aerosol iodine was suggested to be related to the primary iodine release in the tidal zone. Natural meteorological conditions such as solar radiation, raining etc were observed to have influence on the aerosol iodine. During the ship campaign over the North Atlantic Ocean (January – February 2007), the TSI concentrations (medians) ranged 35.14 – 60.63 pmol m-3 among the 5 stages. Likewise, SOI was found to be the most abundant iodine fraction in TSI with a median of 98.6%. Significant correlation also presented between SOI and iodide in the size range of 2.0 – 5.9 μm. Higher iodate concentration was again found in the higher particle size range, similar to that at Mace Head. Airmass transport from the biogenic bloom region and the Antarctic ice front sector was observed to play an important role in aerosol iodine enhancement. The TSI concentrations observed along the 30,000 km long cruise round trip from East Asia to Antarctica during November 2005 – March 2006 were much lower than in the other campaigns, with a median of 6.51 pmol m-3. Approximately 70% of the TSI was SOI on average. The abundances of inorganic iodine including iodine and iodide were less than 30% of TSI. The median value of iodide was 1.49 pmol m-3, which was more than four fold higher than that of iodate (median, 0.28 pmol m-3). Spatial variation indicated highest aerosol iodine appearing in the tropical area. Iodine level was considerably lower in coastal Antarctica with the TSI median of 3.22 pmol m-3. However, airmass transport from the ice front sector was correlated with the enhance TSI level, suggesting the unrevealed source of iodine in the polar region. In addition, significant correlation between SOI and iodide was also shown in this campaign. A global distribution in aerosol was shown in the field campaigns in this work. SOI was verified globally ubiquitous due to the presence in the different sampling locations and its high proportion in TSI in the marine aerosols. The correlations between SOI and iodide were obtained not only in different locations but also in different seasons, implying the possible mechanism of iodide production through SOI decomposition. Nevertheless, future studies are needed for improving the current understanding of iodine chemistry in the MBL (e.g. SOI identification and quantification as well as the update modeling involving organic matters).