Search

Biomass burning particulate matter (BBPM) affects regional air quality and global climate, with impacts expected to continue to grow over the coming years. We show that studies of North American fires have a systematic altitude dependence in measured BBPM normalized excess mixing ratio (NEMR; ΔPM/ΔCO), with airborne and high-altitude studies showing a factor of 2 higher NEMR than ground-based measurements. We report direct airborne measurements of BBPM volatility that partially explain the difference in the BBPM NEMR observed across platforms. We find that when heated to 40-45 °C in an airborne thermal denuder, 19% of lofted smoke PM 1 evaporates. Thermal denuder measurements are consistent with evaporation observed when a single smoke plume was sampled across a range of temperatures as the plume descended from 4 to 2 km altitude. We also demonstrate that chemical aging of smoke and differences in PM emission factors can not fully explain the platform-dependent differences. When the measured PM volatility is applied to output from the High Resolution Rapid Refresh Smoke regional model, we predict a lower PM NEMR at the surface compared to the lofted smoke measured by aircraft. These results emphasize the significant role that gas-particle partitioning plays in determining the air quality impacts of wildfire smoke.

Carbonaceous emissions from wildfires are a dynamic mixture of gases and particles that have important impacts on air quality and climate. Emissions that feed atmospheric models are estimated using burned area and fire radiative power (FRP) methods that rely on satellite products. These approaches show wide variability and have large uncertainties, and their accuracy is challenging to evaluate due to limited aircraft and ground measurements. Here, we present a novel method to estimate fire plume-integrated total carbon and speciated emission rates using a unique combination of lidar remote sensing aerosol extinction profiles and in situ measured carbon constituents. We show strong agreement between these aircraft-derived emission rates of total carbon and a detailed burned area-based inventory that distributes carbon emissions in time using Geostationary Operational Environmental Satellite FRP observations (Fuel2Fire inventory, slope = 1.33 ± 0.04, r 2 = 0.93, and RMSE = 0.27). Other more commonly used inventories strongly correlate with aircraft-derived emissions but have wide-ranging over- and under-predictions. A strong correlation is found between carbon monoxide emissions estimated in situ with those derived from the TROPOspheric Monitoring Instrument (TROPOMI) for five wildfires with coincident sampling windows (slope = 0.99 ± 0.18; bias = 28.5%). Smoke emission coefficients (g MJ -1 ) enable direct estimations of primary gas and aerosol emissions from satellite FRP observations, and we derive these values for many compounds emitted by temperate forest fuels, including several previously unreported species.

Accurate fire emissions inventories are crucial to predict the impacts of wildland fires on air quality and atmospheric composition. Two traditional approaches are widely used to calculate fire emissions: a satellite-based top-down approach and a fuels-based bottom-up approach. However, these methods often considerably disagree on the amount of particulate mass emitted from fires. Previously available observational datasets tended to be sparse, and lacked the statistics needed to resolve these methodological discrepancies. Here, we leverage the extensive and comprehensive airborne in situ and remote sensing measurements of smoke plumes from the recent Fire Influence on Regional to Global Environments and Air Quality (FIREX-AQ) campaign to statistically assess the skill of the two traditional approaches. We use detailed campaign observations to calculate and compare emission rates at an exceptionally high-resolution using three separate approaches: top-down, bottom-up, and a novel approach based entirely on integrated airborne in situ measurements. We then compute the daily average of these high-resolution estimates and compare with estimates from lower resolution, global top-down and bottom-up inventories. We uncover strong, linear relationships between all of the high-resolution emission rate estimates in aggregate, however no single approach is capable of capturing the emission characteristics of every fire. Global inventory emission rate estimates exhibited weaker correlations with the high-resolution approaches and displayed evidence of systematic bias. The disparity between the low-resolution global inventories and the high-resolution approaches is likely caused by high levels of uncertainty in essential variables used in bottom-up inventories and imperfect assumptions in top-down inventories., Competing Interests: The authors declare no conflicts of interest relevant to this study., (© 2021. The Authors.)

Wildfires are a substantial but poorly quantified source of tropospheric ozone (O 3 ). Here, to investigate the highly variable O 3 chemistry in wildfire plumes, we exploit the in situ chemical characterization of western wildfires during the FIREX-AQ flight campaign and show that O 3 production can be predicted as a function of experimentally constrained OH exposure, volatile organic compound (VOC) reactivity, and the fate of peroxy radicals. The O 3 chemistry exhibits rapid transition in chemical regimes. Within a few daylight hours, the O 3 formation substantially slows and is largely limited by the abundance of nitrogen oxides (NO x ). This finding supports previous observations that O 3 formation is enhanced when VOC-rich wildfire smoke mixes into NO x -rich urban plumes, thereby deteriorating urban air quality. Last, we relate O 3 chemistry to the underlying fire characteristics, enabling a more accurate representation of wildfire chemistry in atmospheric models that are used to study air quality and predict climate.

We present a novel method, the Gaussian observational model for edge to center heterogeneity (GOMECH), to quantify the horizontal chemical structure of plumes. GOMECH fits observations of short-lived emissions or products against a long-lived tracer (e.g., CO) to provide relative metrics for the plume width ( w i / w CO ) and center ( b i / w CO ). To validate GOMECH, we investigate OH and NO 3 oxidation processes in smoke plumes sampled during FIREX-AQ (Fire Influence on Regional to Global Environments and Air Quality, a 2019 wildfire smoke study). An analysis of 430 crosswind transects demonstrates that nitrous acid (HONO), a primary source of OH, is narrower than CO ( w HONO / w CO = 0.73-0.84 ± 0.01) and maleic anhydride (an OH oxidation product) is enhanced on plume edges ( w maleicanhydride / w CO = 1.06-1.12 ± 0.01). By contrast, NO 3 production [P(NO 3 )] occurs mainly at the plume center ( w P(NO 3 ) / w CO = 0.91-1.00 ± 0.01). Phenolic emissions, highly reactive to OH and NO 3 , are narrower than CO ( w phenol / w CO = 0.96 ± 0.03, w catechol / w CO = 0.91 ± 0.01, and w methylcatechol / w CO = 0.84 ± 0.01), suggesting that plume edge phenolic losses are the greatest. Yet, nitrophenolic aerosol, their oxidation product, is the greatest at the plume center ( w nitrophenolicaerosol / w CO = 0.95 ± 0.02). In a large plume case study, GOMECH suggests that nitrocatechol aerosol is most associated with P(NO 3 ). Last, we corroborate GOMECH with a large eddy simulation model which suggests most (55%) of nitrocatechol is produced through NO 3 in our case study.

Oceans emit large quantities of dimethyl sulfide (DMS) to the marine atmosphere. The oxidation of DMS leads to the formation and growth of cloud condensation nuclei (CCN) with consequent effects on Earth's radiation balance and climate. The quantitative assessment of the impact of DMS emissions on CCN concentrations necessitates a detailed description of the oxidation of DMS in the presence of existing aerosol particles and clouds. In the unpolluted marine atmosphere, DMS is efficiently oxidized to hydroperoxymethyl thioformate (HPMTF), a stable intermediate in the chemical trajectory toward sulfur dioxide (SO 2 ) and ultimately sulfate aerosol. Using direct airborne flux measurements, we demonstrate that the irreversible loss of HPMTF to clouds in the marine boundary layer determines the HPMTF lifetime ( τ HPMTF < 2 h) and terminates DMS oxidation to SO 2 When accounting for HPMTF cloud loss in a global chemical transport model, we show that SO 2 production from DMS is reduced by 35% globally and near-surface (0 to 3 km) SO 2 concentrations over the ocean are lowered by 24%. This large, previously unconsidered loss process for volatile sulfur accelerates the timescale for the conversion of DMS to sulfate while limiting new particle formation in the marine atmosphere and changing the dynamics of aerosol growth. This loss process potentially reduces the spatial scale over which DMS emissions contribute to aerosol production and growth and weakens the link between DMS emission and marine CCN production with subsequent implications for cloud formation, radiative forcing, and climate., Competing Interests: The authors declare no competing interest.

The ACT-America project is a NASA Earth Venture Suborbital-2 mission designed to study the transport and fluxes of greenhouse gases. The open and freely available ACT-America data sets provide airborne in situ measurements of atmospheric carbon dioxide, methane, trace gases, aerosols, clouds, and meteorological properties, airborne remote sensing measurements of aerosol backscatter, atmospheric boundary layer height and columnar content of atmospheric carbon dioxide, tower-based measurements, and modeled atmospheric mole fractions and regional carbon fluxes of greenhouse gases over the Central and Eastern United States. We conducted 121 research flights during five campaigns in four seasons during 2016-2019 over three regions of the US (Mid-Atlantic, Midwest and South) using two NASA research aircraft (B-200 and C-130). We performed three flight patterns (fair weather, frontal crossings, and OCO-2 underflights) and collected more than 1,140 h of airborne measurements via level-leg flights in the atmospheric boundary layer, lower, and upper free troposphere and vertical profiles spanning these altitudes. We also merged various airborne in situ measurements onto a common standard sampling interval, which brings coherence to the data, creates geolocated data products, and makes it much easier for the users to perform holistic analysis of the ACT-America data products. Here, we report on detailed information of data sets collected, the workflow for data sets including storage and processing of the quality controlled and quality assured harmonized observations, and their archival and formatting for users. Finally, we provide some important information on the dissemination of data products including metadata and highlights of applications of ACT-America data sets., (© 2021. The Authors. Earth and Space Science published by Wiley Periodicals LLC on behalf of American Geophysical Union.)

A particle-into-liquid sampler (PILS) was coupled to a flow-through conductivity cell to provide a continuous, nondestructive, online measurement in support of offline ion chromatography analysis. The conductivity measurement provides a rapid assessment of the total ion concentration augmenting slower batch-sample data from offline analysis and is developed primarily to assist airborne measurements, where fast time-response is essential. A conductivity model was developed for measured ions and excellent closure was derived for laboratory-generated aerosols (97% conductivity explained, R 2 > 0.99). The PILS-conductivity measurement was extensively tested throughout the NASA Cloud, Aerosol and Monsoon Processes: Philippines Experiment (CAMP 2 Ex) during nineteen research flights. A diverse range of ambient aerosol was sampled from biomass burning, fresh and aged urban pollution, and marine sources. Ambient aerosol did not exhibit the same degree of closure as the laboratory aerosol, with measured ions only accountable for 43% of the conductivity. The remaining fraction of the conductivity was examined in combination with ion charge balance and found to provide additional supporting information for diagnosing and modeling particle acidity. An urban plume case study was used to demonstrate the utility of the measurement for supplementing compositional data and augmenting the temporal capability of the PILS., Competing Interests: Disclosure statement The authors declare no conflicts of interest, financial interests or benefits.

Nucleation and growth of aerosol particles from atmospheric vapors constitutes a major source of global cloud condensation nuclei (CCN). The fraction of newly formed particles that reaches CCN sizes is highly sensitive to particle growth rates, especially for particle sizes <10 nm, where coagulation losses to larger aerosol particles are greatest. Recent results show that some oxidation products from biogenic volatile organic compounds are major contributors to particle formation and initial growth. However, whether oxidized organics contribute to particle growth over the broad span of tropospheric temperatures remains an open question, and quantitative mass balance for organic growth has yet to be demonstrated at any temperature. Here, in experiments performed under atmospheric conditions in the Cosmics Leaving Outdoor Droplets (CLOUD) chamber at the European Organization for Nuclear Research (CERN), we show that rapid growth of organic particles occurs over the range from [Formula: see text]C to [Formula: see text]C. The lower extent of autoxidation at reduced temperatures is compensated by the decreased volatility of all oxidized molecules. This is confirmed by particle-phase composition measurements, showing enhanced uptake of relatively less oxygenated products at cold temperatures. We can reproduce the measured growth rates using an aerosol growth model based entirely on the experimentally measured gas-phase spectra of oxidized organic molecules obtained from two complementary mass spectrometers. We show that the growth rates are sensitive to particle curvature, explaining widespread atmospheric observations that particle growth rates increase in the single-digit-nanometer size range. Our results demonstrate that organic vapors can contribute to particle growth over a wide range of tropospheric temperatures from molecular cluster sizes onward., Competing Interests: The authors declare no conflict of interest., (Copyright © 2018 the Author(s). Published by PNAS.)

The San Joaquin Valley (SJV) of California experiences high concentrations of particulate matter NH 4 NO 3 during episodes of meteorological stagnation in winter. A rich data set of observations related to NH 4 NO 3 formation was acquired during multiple periods of elevated NH 4 NO 3 during the Deriving Information on Surface Conditions from Column and Vertically Resolved Observations Relevant to Air Quality (DISCOVER-AQ) field campaign in SJV in January and February 2013. Here NH 4 NO 3 is simulated during the SJV DISCOVER-AQ study period with the Community Multiscale Air Quality (CMAQ) model, diagnostic model evaluation is performed using the DISCOVER-AQ data set, and integrated reaction rate analysis is used to quantify HNO 3 production rates. Simulated NO 3 - generally agrees well with routine monitoring of 24-hr average NO 3 - , but comparisons with hourly average NO 3 - measurements in Fresno revealed differences at higher time resolution. Predictions of gas-particle partitioning of total nitrate (HNO 3 + NO 3 - ) and NHx (NH 3 + NH 4 + ) generally agree well with measurements in Fresno, although partitioning of total nitrate to HNO 3 is sometimes overestimated at low relative humidity in afternoon. Gas-particle partitioning results indicate that NH 4 NO 3 formation is limited by HNO 3 availability in both the model and ambient. NH 3 mixing ratios are underestimated, particularly in areas with large agricultural activity, and additional work on the spatial allocation of NH 3 emissions is warranted. During a period of elevated NH 4 NO 3 , the model predicted that the OH + NO 2 pathway contributed 46% to total HNO 3 production in SJV and the N 2 O 5 heterogeneous hydrolysis pathway contributed 54%. The relative importance of the OH + NO 2 pathway for HNO 3 production is predicted to increase as NOx emissions decrease.

The evolution of atmospheric organic carbon as it undergoes oxidation has a controlling influence on concentrations of key atmospheric species, including particulate matter, ozone and oxidants. However, full characterization of organic carbon over hours to days of atmospheric processing has been stymied by its extreme chemical complexity. Here we study the multigenerational oxidation of α-pinene in the laboratory, characterizing products with several state-of-the-art analytical techniques. Although quantification of some early generation products remains elusive, full carbon closure is achieved (within measurement uncertainty) by the end of the experiments. These results provide new insights into the effects of oxidation on organic carbon properties (volatility, oxidation state and reactivity) and the atmospheric lifecycle of organic carbon. Following an initial period characterized by functionalization reactions and particle growth, fragmentation reactions dominate, forming smaller species. After approximately one day of atmospheric aging, most carbon is sequestered in two long-lived reservoirs-volatile oxidized gases and low-volatility particulate matter.

This study discusses an analysis of combined airborne and ground observations of particulate nitrate (NO 3 - (p) ) concentrations made during the wintertime DISCOVER-AQ study at one of the most polluted cities in the United States, Fresno, CA in the San Joaquin Valley (SJV) and focuses on development of understanding of the various processes that impact surface nitrate concentrations during pollution events. The results provide an explicit case-study illustration of how nighttime chemistry can influence daytime surface-level NO 3 - (p) concentrations, complementing previous studies in the SJV. The observations exemplify the critical role that nocturnal chemical production of NO 3 - (p) aloft in the residual layer (RL) can play in determining daytime surface-level NO 3 - (p) concentrations. Further, they indicate that nocturnal production of NO 3 - (p) in the RL, along with daytime photochemical production, can contribute substantially to the build-up and sustaining of severe pollution episodes. The exceptionally shallow nocturnal boundary layer heights characteristic of wintertime pollution events in the SJV intensifies the importance of nocturnal production aloft in the residual layer to daytime surface concentrations. The observations also demonstrate that dynamics within the RL can influence the early-morning vertical distribution of NO 3 - (p) , despite low wintertime wind speeds. This overnight reshaping of the vertical distribution above the city plays an important role in determining the net impact of nocturnal chemical production on local and regional surface-level NO 3 - (p) concentrations. Entrainment of clean free tropospheric air into the boundary layer in the afternoon is identified as an important process that reduces surface-level NO 3 - (p) and limits build-up during pollution episodes. The influence of dry deposition of HNO 3 gas to the surface on daytime particulate nitrate concentrations is important but limited by an excess of ammonia in the region, which leads to only a small fraction of nitrate existing in the gas-phase even during the warmer daytime. However, in late afternoon, when diminishing solar heating leads to a rapid fall in the mixed boundary layer height, the impact of surface deposition is temporarily enhanced and can lead to a substantial decline in surface-level particulate nitrate concentrations; this enhanced deposition is quickly arrested by a decrease in surface temperature, which drops the gas-phase fraction to near zero. The overall importance of enhanced late afternoon gas-phase loss to the multiday build-up of pollution events is limited by the very shallow nocturnal boundary layer. The case study here demonstrates that mixing down of NO 3 - (p) from the RL can contribute a majority of the surface-level NO 3 - (p) in the morning (here, ~80%), and a strong influence can persist into the afternoon even when photochemical production is maximum. The particular day-to-day contribution of aloft nocturnal NO 3 - (p) production to surface concentrations will depend on prevailing chemical and meteorological conditions. Although specific to the SJV, the observations and conceptual framework further developed here provide general insights into the evolution of pollution episodes in wintertime environments., Competing Interests: Author Contributions The authors declare that they have no conflict of interest.

We report enhancements of glyoxal and methylglyoxal relative to carbon monoxide and formaldehyde in agricultural biomass burning plumes intercepted by the NOAA WP-3D aircraft during the 2013 Southeast Nexus and 2015 Shale Oil and Natural Gas Nexus campaigns. Glyoxal and methylglyoxal were measured using broadband cavity enhanced spectroscopy, which for glyoxal provides a highly selective and sensitive measurement. While enhancement ratios of other species such as methane and formaldehyde were consistent with previous measurements, glyoxal enhancements relative to carbon monoxide averaged 0.0016 ± 0.0009, a factor of 4 lower than values used in global models. Glyoxal enhancements relative to formaldehyde were 30 times lower than previously reported, averaging 0.038 ± 0.02. Several glyoxal loss processes such as photolysis, reactions with hydroxyl radicals, and aerosol uptake were found to be insufficient to explain the lower measured values of glyoxal relative to other biomass burning trace gases, indicating that glyoxal emissions from agricultural biomass burning may be significantly overestimated. Methylglyoxal enhancements were three to six times higher than reported in other recent studies, but spectral interferences from other substituted dicarbyonyls introduce an estimated correction factor of 2 and at least a 25% uncertainty, such that accurate measurements of the enhancements are difficult.

Organic compounds in the atmosphere vary widely in their molecular composition and chemical properties, so no single instrument can reasonably measure the entire range of ambient compounds. Over the past decade, a new generation of in situ, field-deployable mass spectrometers has dramatically improved our ability to detect, identify, and quantify these organic compounds, but no systematic approach has been developed to assess the extent to which currently available tools capture the entire space of chemical identity and properties that is expected in the atmosphere. Reduced-parameter frameworks that have been developed to describe atmospheric mixtures are exploited here to characterize the range of chemical properties accessed by a suite of instruments. Multiple chemical spaces (e.g. oxidation state of carbon vs. volatility, and oxygen number vs. carbon number) were populated with ions measured by several mass spectrometers, with gas- and particle-phase α-pinene oxidation products serving as the test mixture of organic compounds. Few gaps are observed in the coverage of the parameter spaces by the instruments employed in this work, though the full extent to which comprehensive measurement was achieved is difficult to assess due to uncertainty in the composition of the mixture. Overlaps between individual ions and regions in parameter space were identified, both between gas- and particle-phase measurements, and within each phase. These overlaps were conservatively found to account for little (<10%) of the measured mass. However, challenges in identifying overlaps and in accurately converting molecular formulas into chemical properties (such as volatility or reactivity) highlight a continued need to incorporate structural information into atmospheric measurements.

Natural emissions of ozone-and-aerosol-precursor gases such as isoprene and monoterpenes are high in the southeast of the US. In addition, anthropogenic emissions are significant in the Southeast US and summertime photochemistry is rapid. The NOAA-led SENEX (Southeast Nexus) aircraft campaign was one of the major components of the Southeast Atmosphere Study (SAS) and was focused on studying the interactions between biogenic and anthropogenic emissions to form secondary pollutants. During SENEX, the NOAA WP-3D aircraft conducted 20 research flights between 27 May and 10 July 2013 based out of Smyrna, TN. Here we describe the experimental approach, the science goals and early results of the NOAA SENEX campaign. The aircraft, its capabilities and standard measurements are described. The instrument payload is summarized including detection limits, accuracy, precision and time resolutions for all gas-and-aerosol phase instruments. The inter-comparisons of compounds measured with multiple instruments on the NOAA WP-3D are presented and were all within the stated uncertainties, except two of the three NO 2 measurements. The SENEX flights included day- and nighttime flights in the Southeast as well as flights over areas with intense shale gas extraction (Marcellus, Fayetteville and Haynesville shale). We present one example flight on 16 June 2013, which was a daytime flight over the Atlanta region, where several crosswind transects of plumes from the city and nearby point sources, such as power plants, paper mills and landfills, were flown. The area around Atlanta has large biogenic isoprene emissions, which provided an excellent case for studying the interactions between biogenic and anthropogenic emissions. In this example flight, chemistry in and outside the Atlanta plumes was observed for several hours after emission. The analysis of this flight showcases the strategies implemented to answer some of the main SENEX science questions.

Gas-phase low volatility organic compounds (LVOC), produced from oxidation of isoprene 4-hydroxy-3-hydroperoxide (4,3-ISOPOOH) under low-NO conditions, were observed during the FIXCIT chamber study. Decreases in LVOC directly correspond to appearance and growth in secondary organic aerosol (SOA) of consistent elemental composition, indicating that LVOC condense (at OA below 1 μg m(-3)). This represents the first simultaneous measurement of condensing low volatility species from isoprene oxidation in both the gas and particle phases. The SOA formation in this study is separate from previously described isoprene epoxydiol (IEPOX) uptake. Assigning all condensing LVOC signals to 4,3-ISOPOOH oxidation in the chamber study implies a wall-loss corrected non-IEPOX SOA mass yield of ∼4%. By contrast to monoterpene oxidation, in which extremely low volatility VOC (ELVOC) constitute the organic aerosol, in the isoprene system LVOC with saturation concentrations from 10(-2) to 10 μg m(-3) are the main constituents. These LVOC may be important for the growth of nanoparticles in environments with low OA concentrations. LVOC observed in the chamber were also observed in the atmosphere during SOAS-2013 in the Southeastern United States, with the expected diurnal cycle. This previously uncharacterized aerosol formation pathway could account for ∼5.0 Tg yr(-1) of SOA production, or 3.3% of global SOA.

During the Deepwater Horizon (DWH) oil spill, a wide range of gas and aerosol species were measured from an aircraft around, downwind, and away from the DWH site. Additional hydrocarbon measurements were made from ships in the vicinity. Aerosol particles of respirable sizes were on occasions a significant air quality issue for populated areas along the Gulf Coast. Yields of organic aerosol particles and emission factors for other atmospheric pollutants were derived for the sources from the spill, recovery, and cleanup efforts. Evaporation and subsequent secondary chemistry produced organic particulate matter with a mass yield of 8 ± 4% of the oil mixture reaching the water surface. Approximately 4% by mass of oil burned on the surface was emitted as soot particles. These yields can be used to estimate the effects on air quality for similar events as well as for this spill at other times without these data. Whereas emission of soot from burning surface oil was large during the episodic burns, the mass flux of secondary organic aerosol to the atmosphere was substantially larger overall. We use a regional air quality model to show that some observed enhancements in organic aerosol concentration along the Gulf Coast were likely due to the DWH spill. In the presence of evaporating hydrocarbons from the oil, NO(x) emissions from the recovery and cleanup operations produced ozone.

An ultraviolet (UV) optical absorption system has been developed for absolute calibrations of nitric acid (HNO3) and ammonia (NH3) permeation tube emission rates. Using this technique, dilute mixtures containing NH3 or HNO3, both of which interact strongly with many surfaces, are accurately measured at levels below a part per million by volume. This compact and portable instrument operates continuously and autonomously to rapidly (<1 h) quantify the emission of trace gases from permeation devices that are commonly used to calibrate air-monitoring instruments. The output from several HNO3 and NH3 permeation tubes, with emission rates that ranged between 13 and 150 ng/min, was examined as a function of temperature, pressure, and carrier gas flow. Absorptions of 0.015% can be detected which allows a precision (3sigma) of +/-1 ng/min for the HNO3 and NH3 permeation tubes studied here. The accuracy of the measurements, which relies on published UV absorption cross sections, is estimated to be +/-10%. Measurements of permeation tube emission rates using ion chromatography analysis are made to further assess measurement accuracy. The output from the HNO3 and NH3 permeation tubes examined here was stable over the study period, which ranged between 3 months and 1 year for each permeation tube.