ALBERT

Description: A sun-tracking method to improve the pointing accuracy of weather radar Atmospheric Measurement Techniques, 5, 547-555, 2012 Author(s): X. Muth, M. Schneebeli, and A. Berne Accurate positioning of data collected by a weather radar is of primary importance for their appropriate georeferencing, which in turn makes it possible to combine those with additional sources of information (topography, land cover maps, meteorological simulations from numerical weather models to list a few). This issue is especially acute for mobile radar systems, for which accurate and stable leveling might be difficult to ensure. The sun is a source of microwave radiation, which can be detected by weather radars and used for accurate positioning of radar data. This paper presents a technique based on the similarity between theodolites and radar systems as well as on the sun echoes to quantify and hence correct the instrumental errors which can affect the pointing accuracy of radar antenna. The proposed method is applied to data collected in the Swiss Alps using a mobile X-band radar system. The obtained instrumental bias values are evaluated by comparing the locations of the ground echoes predicted using these bias estimates with the observed ground echo locations. The very good agreement between the two confirms the accuracy of the proposed method.

Description: Smoke aerosol and its radiative effects during extreme fire event over Central Russia in summer 2010 Atmospheric Measurement Techniques, 5, 557-568, 2012 Author(s): N. Chubarova, Ye. Nezval', I. Sviridenkov, A. Smirnov, and I. Slutsker Different microphysical, optical and radiative properties of aerosol were analyzed during the severe fires in summer 2010 over Central Russia using ground measurements at two AERONET sites in Moscow (Meteorological Observatory of Moscow State University – MSU MO) and Zvenigorod (Moscow Region) and radiative measurements at the MSU MO. Volume aerosol size distribution in smoke conditions had a bimodal character with the significant prevalence of fine mode particles, for which effective radius was shifted to higher values ( r eff-fine = 0.24 μm against approximately 0.15 μm in typical conditions). For smoke aerosol, the imaginary part of refractive index (REFI) in the visible spectral region was lower than that for typical aerosol (REFI λ =675 nm = 0.006 against REFI λ =675 nm = 0.01), while single scattering albedo (SSA) was significantly higher (SSA λ =675 nm = 0.95 against SSA λ =675 nm ~ 0.9). Extremely high aerosol optical thickness at 500 nm (AOT500) was observed on 6–8 August reaching the absolute maximum on 7 August in Moscow (AOT500 = 6.4) and at Zvenigorod (AOT500 = 5.9). A dramatic attenuation of solar irradiance at ground was also recorded. Maximum irradiance loss had reached 64% for global shortwave irradiance, 91% for UV radiation 300–380 nm, and 97% for erythemally-weighted UV irradiance at relatively high solar elevation 47°. Significant spectral dependence in attenuation of solar irradiance in smoky conditions was mainly explained by higher AOT and smaller SSA in UV (0.8–0.9) compared with SSA in the visible region of spectrum. The assessments of radiative forcing effect (RFE) at the TOA indicated a significant cooling of the smoky atmosphere. Instant RFE reached −167 Wm −2 at AOT500 = 6.4, climatological RFE calculated with August 2010 monthly mean AOT was about −65 Wm −2 , compared with −20 Wm −2 for typical aerosol according to the 10 yr period of measurements in Moscow.

Description: Towards a regional ocean forecasting system for the IBI (Iberia-Biscay-Ireland area): developments and improvements within the ECOOP project framework Ocean Science, 8, 143-159, 2012 Author(s): S. Cailleau, J. Chanut, J.-M. Lellouche, B. Levier, C. Maraldi, G. Reffray, and M. G. Sotillo The regional ocean operational system remains a key element in downscaling from large scale (global or basin scale) systems to coastal ones. It enables the transition between systems in which the resolution and the resolved physics are quite different. Indeed, coastal applications need a system to predict local high frequency events (inferior to the day) such as storm surges, while deep sea applications need a system to predict large scale lower frequency ocean features. In the framework of the ECOOP project, a regional system for the Iberia-Biscay-Ireland area has been upgraded from an existing V0 version to a V2. This paper focuses on the improvements from the V1 system, for which the physics are close to a large scale basin system, to the V2 for which the physics are more adapted to shelf and coastal issues. Strong developments such as higher regional physics resolution in the NEMO Ocean General Circulation Model for tides, non linear free surface and adapted vertical mixing schemes among others have been implemented in the V2 version. Thus, regional thermal fronts due to tidal mixing now appear in the latest version solution and are quite well positioned. Moreover, simulation of the stratification in shelf areas is also improved in the V2.

Description: The Nexus Land-Use model version 1.0, an approach articulating biophysical potentials and economic dynamics to model competition for land-use Geoscientific Model Development Discussions, 5, 571-638, 2012 Author(s): F. Souty, T. Brunelle, P. Dumas, B. Dorin, P. Ciais, R. Crassous, C. Müller, and A. Bondeau Interactions between food demand, biomass energy and forest preservation are driving both food prices and land-use changes, regionally and globally. This study presents a new model called Nexus Land-Use version 1.0 which describes these interactions through a generic representation of agricultural intensification mechanisms. The Nexus Land-Use model equations combine biophysics and economics into a single coherent framework to calculate crop yields, food prices, and resulting pasture and cropland areas within 12 regions inter-connected with each other by international trade. The representation of cropland and livestock production systems in each region relies on three components: (i) a biomass production function derived from the crop yield response function to inputs such as industrial fertilisers; (ii) a detailed representation of the livestock production system subdivided into an intensive and an extensive component, and (iii) a spatially explicit distribution of potential (maximal) crop yields prescribed from the Lund-Postdam-Jena global vegetation model for managed Land (LPJmL). The economic principles governing decisions about land-use and intensification are adapted from the Ricardian rent theory, assuming cost minimisation for farmers. The land-use modelling approach described in this paper entails several advantages. Firstly, it makes it possible to explore interactions among different types of biomass demand for food and animal feed, in a consistent approach, including indirect effects on land-use change resulting from international trade. Secondly, yield variations induced by the possible expansion of croplands on less suitable marginal lands are modelled by using regional land area distributions of potential yields, and a calculated boundary between intensive and extensive production. The model equations and parameter values are first described in details. Then, idealised scenarios exploring the impact of forest preservation policies or rising energy price on agricultural intensification are described, and their impacts on pasture and cropland areas are investigated.

Description: Evaluation of Arctic broadband surface radiation measurements Atmospheric Measurement Techniques, 5, 429-438, 2012 Author(s): N. Matsui, C. N. Long, J. Augustine, D. Halliwell, T. Uttal, D. Longenecker, O. Niebergall, J. Wendell, and R. Albee The Arctic is a challenging environment for making in-situ surface radiation measurements. A standard suite of radiation sensors is typically designed to measure incoming and outgoing shortwave (SW) and thermal infrared, or longwave (LW), radiation. Enhancements may include various sensors for measuring irradiance in narrower bandwidths. Many solar radiation/thermal infrared flux sensors utilize protective glass domes and some are mounted on complex mechanical platforms (solar trackers) that keep sensors and shading devices trained on the sun along its diurnal path. High quality measurements require striking a balance between locating stations in a pristine undisturbed setting free of artificial blockage (such as from buildings and towers) and providing accessibility to allow operators to clean and maintain the instruments. Three significant sources of erroneous data in the Arctic include solar tracker malfunctions, rime/frost/snow deposition on the protective glass domes of the radiometers and operational problems due to limited operator access in extreme weather conditions. In this study, comparisons are made between the global and component sum (direct [vertical component] + diffuse) SW measurements. The difference between these two quantities (that theoretically should be zero) is used to illustrate the magnitude and seasonality of arctic radiation flux measurement problems. The problem of rime/frost/snow deposition is investigated in more detail for one case study utilizing both SW and LW measurements. Solutions to these operational problems that utilize measurement redundancy, more sophisticated heating and ventilation strategies and a more systematic program of operational support and subsequent data quality protocols are proposed.

Description: Modelling sub-grid wetland in the ORCHIDEE global land surface model: evaluation against river discharges and remotely sensed data Geoscientific Model Development Discussions, 5, 683-735, 2012 Author(s): B. Ringeval, B. Decharme, S. L. Piao, P. Ciais, F. Papa, N. de Noblet-Ducoudré, C. Prigent, P. Friedlingstein, I. Gouttevin, C. Koven, and A. Ducharne The quality of the global hydrological simulations performed by Land Surface Models (LSMs) strongly depends on processes that occur at unresolved spatial scales. Approaches such as TOPMODEL have been developed, which allow soil moisture redistribution within each grid-cell, based upon sub-grid scale topography. Moreover, the coupling between TOPMODEL and a LSM appears as a potential way to simulate wetland extent dynamic and its sensitivity to climate, a recently identified research problem for biogeochemical modelling, including methane emissions. Global evaluation of the coupling between TOPMODEL and an LSM is difficult, and prior attempts have been indirect, based on the evaluation of the simulated river flow. This study presents a new way to evaluate this coupling, within the ORCHIDEE LSM, using remote sensing data of inundated areas. Because of differences in nature between the satellite derived information – inundation extent – and the variable diagnosed by TOPMODEL/ORCHIDEE – area at maximum soil water content –, the evaluation focuses on the spatial distribution of these two quantities as well as on their temporal variation. Despite some difficulties in exactly matching observed localized inundated events, we obtain a rather good agreement in the distribution of these two quantities at a global scale. Floodplains are not accounted for in the model, and this is a major limitation. The difficulty of reproducing the year-to-year variability of the observed inundated area (for instance, the decreasing trend by the end of 90s) is also underlined. Classical indirect evaluation based on comparison between simulated and observed riverflow is also performed and underlines difficulties to simulate riverflow after coupling with TOPMODEL. The relationship between inundation and river flow at the basin scale in the model is analyzed, using both methods (evaluation against remote sensing data and riverflow). Finally, we discuss the potential of the TOPMODEL/LSM coupling to simulate wetland areas. A major limitation of the coupling for this purpose is linked to its ability to simulate a global wetland coverage consistent with the commonly used datasets. However, it seems to be a good opportunity to account for the wetland areas sensitivity to the climate and thus to simulate its temporal variability.

Description: Technical Note: Detection of gas bubble leakage via correlation of water column multibeam images Ocean Science, 8, 175-181, 2012 Author(s): J. Schneider von Deimling and C. Papenberg Hydroacoustic detection of natural gas release from the seafloor has been conducted in the past by using singlebeam echosounders. In contrast, modern multibeam swath mapping systems allow much wider coverage, higher resolution, and offer 3-D spatial correlation. Up to the present, the extremely high data rate hampers water column backscatter investigations and more sophisticated visualization and processing techniques are needed. Here, we present water column backscatter data acquired with a 50 kHz prototype multibeam system over a period of 75 seconds. Display types are of swath-images as well as of a "re-sorted" singlebeam presentation. Thus, individual and/or groups of gas bubbles rising from the 24 m deep seafloor clearly emerge in the acoustic images, making it possible to estimate rise velocities. A sophisticated processing scheme is introduced to identify those rising gas bubbles in the hydroacoustic data. We apply a cross-correlation technique adapted from particle imaging velocimetry (PIV) to the acoustic backscatter images. Temporal and spatial drift patterns of the bubbles are assessed and are shown to match very well to measured and theoretical rise patterns. The application of this processing to our field data gives clear results with respect to unambiguous bubble detection and remote bubble rise velocimetry. The method can identify and exclude the main source of misinterpretations, i.e. fish-mediated echoes. Although image-based cross-correlation techniques are well known in the field of fluid mechanics for high resolution and non-inversive current flow field analysis, we present the first application of this technique as an acoustic bubble detector.

Description: Black Sea coastal forecasting system Ocean Science, 8, 183-196, 2012 Author(s): A. I. Kubryakov, G. K. Korotaev, V. L. Dorofeev, Y. B. Ratner, A. Palazov, N. Valchev, V. Malciu, R. Matescu, and T. Oguz The Black Sea coastal nowcasting and forecasting system was built within the framework of EU FP6 ECOOP (European COastalshelf sea OPerational observing and forecasting system) project for five regions: the south-western basin along the coasts of Bulgaria and Turkey, the north-western shelf along the Romanian and Ukrainian coasts, coastal zone around of the Crimea peninsula, the north-eastern Russian coastal zone and the coastal zone of Georgia. The system operates in the real-time mode during the ECOOP project and afterwards. The forecasts include temperature, salinity and current velocity fields. Ecosystem model operates in the off-line mode near the Crimea coast.

Description: Sea surface temperature anomalies, seasonal cycle and trend regimes in the Eastern Pacific coast Ocean Science, 8, 81-90, 2012 Author(s): A. Ramos-Rodríguez, D. B. Lluch-Cota, S. E. Lluch-Cota, and A. Trasviña-Castro We used the extended reconstruction of sea surface temperature (ERSST) to analyze the variation of surface temperature and the seasonal cycle along the coast of the eastern Pacific (60° N–60° S, 61 pixels alongshore) from 1950 to 2010 (732 months). First, we analyzed the monthly anomalies and looked for a relationship of such anomalies with total solar irradiance (TSI) and then the Regime Shift Detector (RSD) was applied to detect possible temperature regimes in the series. Afterwards, we calculated a yearly temperature range per pixel (amplitude of seasonal cycle) and through the subtraction of a latitudinal theoretical curve of temperature based on solar irradiance, the residuals of the seasonal cycle were obtained. The results showed an almost complete spatial synchrony and dominance of negative anomalies from 1950 to mid-late 1970's, with a switch to near-zero and positive anomalies that lasted up to late 1990's when a new shift to negative values was detected. Such a shift lasted until the early 2000's when positive anomalies appeared again but there was a change to negative anomalies in the late 2000's. These results were supported by the RSD. The TSI variability shows a clear relationship with that of sea surface temperature anomalies and with the regime changes. This is probably due to a difference in the amount of energy received from the sun. Comparing the "cool regime" versus the "warm regime", the second one received 0.39% more energy (approximately 3 × 10 8 J m −2 ) from the sun. Seasonal cycles show larger ranges at northern latitudes (〉40° N), northern tropical-temperate transition zone (20°–26° N) and in the tropical-equatorial band (0°–30° S). The smallest ranges occur at 0°–16° N and 50°–60° S. The residuals (seasonal minus the theoretical curve) indicated a clear modulation due to advection by ocean currents.

Description: Profiles of CH 4 , HDO, H 2 O, and N 2 O with improved lower tropospheric vertical resolution from Aura TES radiances Atmospheric Measurement Techniques, 5, 397-411, 2012 Author(s): J. Worden, S. Kulawik, C. Frankenberg, V. Payne, K. Bowman, K. Cady-Peirara, K. Wecht, J.-E. Lee, and D. Noone Thermal infrared (IR) radiances measured near 8 microns contain information about the vertical distribution of water vapor (H 2 O), the water isotopologue HDO, and methane (CH 4 ), key gases in the water and carbon cycles. Previous versions (Version 4 or less) of the TES profile retrieval algorithm used a "spectral-window" approach to minimize uncertainty from interfering species at the expense of reduced vertical resolution and sensitivity. In this manuscript we document changes to the vertical resolution and uncertainties of the TES version 5 retrieval algorithm. In this version (Version 5), joint estimates of H 2 O, HDO, CH 4 and nitrous oxide (N 2 O) are made using radiances from almost the entire spectral region between 1100 cm −1 and 1330 cm −1 . The TES retrieval constraints are also modified in order to better use this information. The new H 2 O estimates show improved vertical resolution in the lower troposphere and boundary layer, while the new HDO/H 2 O estimates can now profile the HDO/H 2 O ratio between 925 hPa and 450 hPa in the tropics and during summertime at high latitudes. The new retrievals are now sensitive to methane in the free troposphere between 800 and 150 mb with peak sensitivity near 500 hPa; whereas in previous versions the sensitivity peaked at 200 hPa. However, the upper troposphere methane concentrations are biased high relative to the lower troposphere by approximately 4% on average. This bias is likely related to temperature, calibration, and/or methane spectroscopy errors. This bias can be mitigated by normalizing the CH 4 estimate by the ratio of the N 2 O estimate relative to the N 2 O prior, under the assumption that the same systematic error affects both the N 2 O and CH 4 estimates. We demonstrate that applying this ratio theoretically reduces the CH 4 estimate for non-retrieved parameters that jointly affect both the N 2 O and CH 4 estimates. The relative upper troposphere to lower troposphere bias is approximately 2.8% after this bias correction. Quality flags based upon the vertical variability of the methane and N 2 O estimates can be used to reduce this bias further. While these new CH 4 , HDO/H 2 O, and H 2 O estimates are consistent with previous TES retrievals in the altitude regions where the sensitivities overlap, future comparisons with independent profile measurement will be required to characterize the biases of these new retrievals and determine if the calculated uncertainties using the new constraints are consistent with actual uncertainties.

Description: Two instruments based on differential optical absorption spectroscopy (DOAS) to measure accurate ammonia concentrations in the atmosphere Atmospheric Measurement Techniques, 5, 413-427, 2012 Author(s): H. Volten, J. B. Bergwerff, M. Haaima, D. E. Lolkema, A. J. C. Berkhout, G. R. van der Hoff, C. J. M. Potma, R. J. Wichink Kruit, W. A. J. van Pul, and D. P. J. Swart We present two Differential Optical Absorption Spectroscopy (DOAS) instruments built at RIVM: the RIVM DOAS and the miniDOAS. Both instruments provide virtually interference-free measurements of NH 3 concentrations in the atmosphere, since they measure over an open path, without suffering from inlet problems or interference problems by ammonium aerosols dissociating on tubes or filters. They measure concentrations up to at least 200 μg m −3 , have a fast response, low maintenance demands, and a high up-time. The RIVM DOAS has a high accuracy of typically 0.15 μg m −3 for ammonia for 5-min averages and over a total light path of 100 m. The miniDOAS has been developed for application in measurement networks such as the Dutch National Air Quality Monitoring Network (LML). Compared to the RIVM DOAS it has a similar accuracy, but is significantly reduced in size, costs, and handling complexity. The RIVM DOAS and miniDOAS results showed excellent agreement ( R 2 = 0.996) during a field measurement campaign in Vredepeel, the Netherlands. This measurement site is located in an agricultural area and is characterized by highly variable, but on average high ammonia concentrations in the air. The RIVM-DOAS and miniDOAS results were compared to the results of the AMOR instrument, a continuous-flow wet denuder system, which is currently used in the LML. Averaged over longer time spans of typically a day, the (mini)DOAS and AMOR results agree reasonably well, although an offset of the AMOR values compared to the (mini)DOAS results exists. On short time scales, the (mini)DOAS shows a faster response and does not show the memory effects due to inlet tubing and transport of absorption fluids encountered by the AMOR. Due to its high accuracy, high uptime, low maintenance and its open path, the (mini)DOAS shows a good potential for flux measurements by using two (or more) systems in a gradient set-up and applying the aerodynamic gradient technique.

Description: Evaluation of a near-global eddy-resolving ocean model Geoscientific Model Development Discussions, 5, 4305-4354, 2012 Author(s): P. R. Oke, D. A. Griffin, A. Schiller, R. J. Matear, R. Fiedler, J. Mansbridge, A. Lenton, M. Cahill, M. A. Chamberlain, and K. Ridgway Analysis of the variability in an 18-yr run of a near-global, eddy-resolving ocean general circulation model coupled with biogeochemistry is presented. Comparisons between modelled and observed mean sea level (MSL), mixed-layer depth (MLD), sea-level anomaly (SLA), sea-surface temperature (SST), and Chlorophyll a indicate that the model variability is realistic. We find some systematic errors in the modelled MLD, with the model generally deeper than observations, that results in errors in the Chlorophyll a , owing to the strong biophysical coupling. We evaluate several other metrics in the model, including the zonally-averaged seasonal cycle of SST, meridional overturning, volume transports through key Straits and passages, zonal averaged temperature and salinity, and El Nino-related SST indices. We find that the modelled seasonal cycle in SST is 0.5–1.5 °C weaker than observed; volume transports of the Antarctic Circumpolar Current, the East Australian Current, and Indonesian Throughflow are in good agreement with observational estimates; and the correlation between the modelled and observed NINO SST indices exceed 0.91. Most aspects of the model circulation are realistic. We conclude that the model output is suitable for broader analysis to better understand ocean dynamics and ocean variability.

Description: Technical Note: The single particle soot photometer fails to reliably detect PALAS soot nanoparticles Atmospheric Measurement Techniques, 5, 3099-3107, 2012 Author(s): M. Gysel, M. Laborde, A. A. Mensah, J. C. Corbin, A. Keller, J. Kim, A. Petzold, and B. Sierau The single particle soot photometer (SP2) uses laser-induced incandescence (LII) for the measurement of atmospheric black carbon (BC) particles. The BC mass concentration is obtained by combining quantitative detection of BC mass in single particles with a counting efficiency of 100% above its lower detection limit. It is commonly accepted that a particle must contain at least several tenths of a femtogram BC in order to be detected by the SP2. Here we show the result that most BC particles from a PALAS spark discharge soot generator remain undetected by the SP2, even if their BC mass, as independently determined with an aerosol particle mass analyser (APM), is clearly above the typical lower detection limit of the SP2. Comparison of counting efficiency and effective density data of PALAS soot with flame generated soot (combustion aerosol standard burner, CAST), fullerene soot and carbon black particles (Cabot Regal 400R) reveals that particle morphology can affect the SP2's lower detection limit. PALAS soot particles are fractal-like agglomerates of very small primary particles with a low fractal dimension, resulting in a very low effective density. Such loosely packed particles behave like "the sum of individual primary particles" in the SP2's laser. Accordingly, most PALAS soot particles remain undetected as the SP2's laser intensity is insufficient to heat the primary particles to their vaporisation temperature because of their small size ( D pp ≈ 5–10 nm). Previous knowledge from pulsed laser-induced incandescence indicated that particle morphology might have an effect on the SP2's lower detection limit, however, an increase of the lower detection limit by a factor of ∼5–10, as reported here for PALAS soot, was not expected. In conclusion, the SP2's lower detection limit at a certain laser power depends primarily on the total BC mass per particle for compact particles with sufficiently high effective density. By contrast, the BC mass per primary particle mainly determines whether fractal-like particles with low fractal dimension and very small primary particles are detectable, while their total BC mass has only a minor influence. This effect shifts the lower detection limit to much higher BC mass, or makes them completely undetectable. Consequently, care has to be taken when using the SP2 in applications dealing with loosely packed particles that have very small primary particles as building blocks.

Description: Evaluating a lightning parameterization based on cloud-top height for mesoscale numerical model simulations Geoscientific Model Development Discussions, 5, 3493-3531, 2012 Author(s): J. Wong, M. C. Barth, and D. Noone The Price and Rind lightning parameterization based on cloud-top height is a commonly used method for predicting flash rate in global chemistry models. As mesoscale simulations begin to implement flash rate predictions at resolutions that partially resolve convection, it is necessary to validate and understand the behavior of this method within such regime. In this study, we tested the flash rate parameterization, intra-cloud/cloud-to-ground (IC:CG) partitioning parameterization, and the associated resolution dependency "calibration factor" by Price and Rind using the Weather Research and Forecasting (WRF) model running at 36 km, 12 km, and 4 km grid spacings within the continental United States. Our results show that while the integrated flash count is consistent with observation when model biases in convection are taken into account, an erroneous frequency distribution is simulated. When the spectral characteristics of lightning flash rate is a concern, we recommend the use of prescribed IC:CG values. In addition, using cloud-top from convective parameterization, the "calibration factor" is also shown to be insufficient in reconciling the resolution dependency at the tested grid spacing used in this study. We recommend scaling by areal ratio relative to a base-case grid spacing determined by convective core density.

Description: Using model reduction to predict the soil-surface C 18 OO flux: an example of representing complex biogeochemical dynamics in a computationally efficient manner Geoscientific Model Development Discussions, 5, 3469-3491, 2012 Author(s): W. J. Riley Earth System Models (ESMs) must calculate large-scale interactions between the land and atmosphere while accurately characterizing fine-scale spatial heterogeneity in water, carbon, and nutrient dynamics. We present here a high-dimensional model representation (HDMR) approach that allows detailed process representation of a coupled carbon and water tracer (the δ 18 O value of the soil-surface CO 2 flux (δ F s )) in a computationally tractable manner. δ F s depends on the δ 18 O value of soil water, soil moisture, soil temperature, and soil CO 2 production (all of which are depth-dependent), and the δ 18 O value of above-surface CO 2 . We tested the HDMR approach over a growing season in a C 4 -dominated pasture using two vertical soil discretizations. The difference between the HDMR approach and the full model solution in the three-month integrated isoflux was less than 0.2% (0.5 mol m −2 ‰), and the approach is up to 100 times faster than the full numerical solution. This type of model reduction approach allows representation of complex coupled biogeochemical processes in regional and global climate models and can be extended to characterize subgrid-scale spatial heterogeneity.

Description: ECOCLIMAP-II/Europe: a twofold database of ecosystems and surface parameters at 1-km resolution based on satellite information for use in land surface, meteorological and climate models Geoscientific Model Development Discussions, 5, 3573-3620, 2012 Author(s): S. Faroux, A. T. Kaptué Tchuenté, J.-L. Roujean, V. Masson, E. Martin, and P. Le Moigne The overall objective of the present study is to introduce the new ECOCLIMAP-II database for Europe, which is an upgrade for this region of the former initiative, ECOCLIMAP-I, already implemented at global scale. The ECOCLIMAP programme is a dual database at 1-km resolution that includes an ecosystem classification and a coherent set of land surface parameters that are primarily mandatory in meteorological modelling (notably leaf area index and albedo). Hence, the aim of this innovative physiography is to enhance the quality of initialisation and impose some surface attributes within the scope of weather forecasting and climate related studies. The strategy for implementing ECOCLIMAP-II is to depart from prevalent land cover products such as CLC2000 (Corine Land Cover) and GLC2000 (Global Land Cover) by splitting existing classes into new classes that possess a better regional character by virtue of the climatic environment (latitude, proximity to the sea, topography). The leaf area index (LAI) from MODIS and NDVI from SPOT/Vegetation yield the two proxy variables that were considered here in order to perform a multi-year trimmed analysis between 1999 and 2005 using the K-means method. Further, meteorological applications require each land cover type to appear as a partition of fractions of 4 main surface types or tiles (nature, water bodies, sea, urban areas) and, inside the nature tile, fractions of 12 Plant Functional Types (PFTs) representing generic vegetation types – principally broadleaf forest, needleleaf forest, C3 and C4 crops, grassland and bare land – as incorporated by the SVAT model ISBA developed at Météo France. This landscape division also forms the cornerstone of a validation exercise. The new ECOCLIMAP-II can be verified with auxiliary land cover products at very fine and coarse resolutions by means of versatile land occupation nomenclatures.

Description: Microstructure observations during the spring 2011 STRATIPHYT-II cruise in the northeast Atlantic Ocean Science, 8, 945-957, 2012 Author(s): E. Jurado, H. A. Dijkstra, and H. J. van der Woerd Small-scale temperature and conductivity variations have been measured in the upper 100 m of the northeast Atlantic during the STRATIPHYT-II cruise (Las Palmas–Reykjavik, 6 April–3 May 2011). The measurements were done at midday and comprised 2 to 15 vertical profiles at each station. The derived turbulent quantities show a transition between weakly-stratified (mixed layer depth, MLD, 100), which was centered at about 48° N. The temperature eddy diffusivities, K T , range from 10 −5 to 10 0 m 2 s −1 in the weakly-stratified stations, and range from 3 × 10 −4 to 2 × 10 0 m 2 s −1 in the well-mixed stations. The turbulent kinetic energy dissipation rates, ε, range from 3 × 10 −8 to 2 × 10 −6 m 2 s −3 south of the transition zone, and from 10 −7 to 10 −5 m 2 s −3 north of the transition zone. The station-averaged K T values throughout the mixed layer increase exponentially with the wind speed. The station-averaged ε values throughout the mixed layer scale with the wind stress similarity variable with a scaling factor of about 1.8 in the wind-dominated stations (ε ≈ 1.8 u ☆ 3 /(−κ z )). The values of K T and ε are on average 10 times higher compared to the values measured at the same stations in July 2009. The results presented here constitute a unique data set giving large spatial coverage of upper ocean spring turbulence quantities.

Description: Intercomparison of temperature trends in IPCC CMIP5 simulations with observations, reanalyses and CMIP3 models Geoscientific Model Development Discussions, 5, 3621-3645, 2012 Author(s): J. Xu and A. M. Powell Jr. On the basis of the fifth Coupled Model Intercomparison Project (CMIP5) and the climate model simulations covering 1979 through 2005, the temperature trends and their uncertainties have been examined to note the similarities or differences compared to the radiosonde observations, reanalyses and the third Coupled Model Intercomparison Project (CMIP3) simulations. The results show noticeable discrepancies for the estimated temperature trends in the four data groups (Radiosonde, Reanalysis, CMIP3 and CMIP5) although similarities can be observed. Compared to the CMIP3 model simulations, the simulation in some of CMIP5 models were improved. The CMIP5 models displayed a negative temperature trend in the stratosphere closer to the strong negative trend seen in the observations. However, the positive tropospheric trend in the tropics is overestimated by the CMIP5 models relative to CMIP3 models. While some of the models produce temperature trend patterns more highly correlated with the observed patterns in CMIP5, the other models (such as CCSM4 and IPSL_CM5A-LR) exhibit the reverse tendency. The CMIP5 temperature trend uncertainty was significantly reduced in most areas, especially in the Arctic and Antarctic stratosphere, compared to the CMIP3 simulations. Similar to the CMIP3, the CMIP5 simulations overestimated the tropospheric warming in the tropics and Southern Hemisphere and underestimated the stratospheric cooling. The crossover point where tropospheric warming changes into stratospheric cooling occurred near 100 hPa in the tropics, which is higher than in the radiosonde and reanalysis data. The result is likely related to the overestimation of convective activity over the tropical areas in both the CMIP3 and CMIP5 models. Generally, for the temperature trend estimates associated with the numerical models including the reanalyses and global climate models, the uncertainty in the stratosphere is much larger than that in the troposphere, and the uncertainty in the Antarctic is the largest. In addition, note that the reanalyses show the largest uncertainty in the lower tropical stratosphere, and the CMIP3 simulations show the largest uncertainty in both the south and north polar regions.

Description: Effect of air composition (N 2 , O 2 , Ar, and H 2 O) on CO 2 and CH 4 measurement by wavelength-scanned cavity ring-down spectroscopy: calibration and measurement strategy Atmospheric Measurement Techniques, 5, 2689-2701, 2012 Author(s): H. Nara, H. Tanimoto, Y. Tohjima, H. Mukai, Y. Nojiri, K. Katsumata, and C. W. Rella We examined potential interferences from water vapor and atmospheric background gases (N 2 , O 2 , and Ar), and biases by isotopologues of target species, on accurate measurement of atmospheric CO 2 and CH 4 by means of wavelength-scanned cavity ring-down spectroscopy (WS-CRDS). Changes of the background gas mole fractions in the sample air substantially impacted the CO 2 and CH 4 measurements: variation of CO 2 and CH 4 due to relative increase of each background gas increased as Ar 〈 O 2 〈 N 2 , suggesting similar relation for the pressure-broadening effects (PBEs) among the background gas. The pressure-broadening coefficients due to variations in O 2 and Ar for CO 2 and CH 4 are empirically determined from these experimental results. Calculated PBEs using the pressure-broadening coefficients are linearly correlated with the differences between the mole fractions of O 2 and Ar and their ambient abundances. Although the PBEs calculation showed that impact of natural variation of O 2 is negligible on the CO 2 and CH 4 measurements, significant bias was inferred for the measurement of synthetic standard gases. For gas standards balanced with purified air, the PBEs were estimated to be marginal (up to 0.05 ppm for CO 2 and 0.01 ppb for CH 4 ) although the PBEs were substantial (up to 0.87 ppm for CO 2 and 1.4 ppb for CH 4 ) for standards balanced with synthetic air. For isotopic biases on CO 2 measurements, we compared experimental results and theoretical calculations, which showed excellent agreement within their uncertainty. We derived instrument-specific water correction functions empirically for three WS-CRDS instruments (Picarro EnviroSense 3000i, G-1301, and G-2301), and evaluated the transferability of the water correction function from G-1301 among these instruments. Although the transferability was not proven, no significant difference was found in the water vapor correction function for the investigated WS-CRDS instruments as well as the instruments reported in the past studies within the typical analytical precision at sufficiently low water concentrations (

Description: Atmospheric effect on the ground-based measurements of broadband surface albedo Atmospheric Measurement Techniques, 5, 2675-2688, 2012 Author(s): T. Manninen, A. Riihelä, and G. de Leeuw Ground-based pyranometer measurements of the (clear-sky) broadband surface albedo are affected by the atmospheric conditions (mainly by aerosol particles, water vapour and ozone). A new semi-empirical method for estimating the magnitude of the effect of atmospheric conditions on surface albedo measurements in clear-sky conditions is presented. Global and reflected radiation and/or aerosol optical depth (AOD) at two wavelengths are needed to apply the method. Depending on the aerosol optical depth and the solar zenith angle values, the effect can be as large as 20%. For the cases we tested using data from the Cabauw atmospheric test site in the Netherlands, the atmosphere caused typically up to 5% overestimation of surface albedo with respect to corresponding black-sky surface albedo values.

Description: Improved Micro Rain Radar snow measurements using Doppler spectra post-processing Atmospheric Measurement Techniques, 5, 2661-2673, 2012 Author(s): M. Maahn and P. Kollias The Micro Rain Radar 2 (MRR) is a compact Frequency Modulated Continuous Wave (FMCW) system that operates at 24 GHz. The MRR is a low-cost, portable radar system that requires minimum supervision in the field. As such, the MRR is a frequently used radar system for conducting precipitation research. Current MRR drawbacks are the lack of a sophisticated post-processing algorithm to improve its sensitivity (currently at +3 dBz), spurious artefacts concerning radar receiver noise and the lack of high quality Doppler radar moments. Here we propose an improved processing method which is especially suited for snow observations and provides reliable values of effective reflectivity, Doppler velocity and spectral width. The proposed method is freely available on the web and features a noise removal based on recognition of the most significant peak. A dynamic dealiasing routine allows observations even if the Nyquist velocity range is exceeded. Collocated observations over 115 days of a MRR and a pulsed 35.2 GHz MIRA35 cloud radar show a very high agreement for the proposed method for snow, if reflectivities are larger than −5 dBz. The overall sensitivity is increased to −14 and −8 dBz, depending on range. The proposed method exploits the full potential of MRR's hardware and substantially enhances the use of Micro Rain Radar for studies of solid precipitation.

Description: Comparison of satellite microwave backscattering (ASCAT) and visible/near-infrared reflectances (PARASOL) for the estimation of aeolian aerodynamic roughness length in arid and semi-arid regions Atmospheric Measurement Techniques, 5, 2703-2712, 2012 Author(s): C. Prigent, C. Jiménez, and J. Catherinot Previous studies examined the possibility to estimate the aeolian aerodynamic roughness length from satellites, either from visible/near-infrared observations or from microwave backscattering measurements. Here we compare the potential of the two approaches and propose to merge the two sources of information to benefit from their complementary aspects, i.e. the high spatial resolution of the visible/near-infrared (6 km for PARASOL that is part of the A-Train) and the independence from atmospheric contamination of the active microwaves (ASCAT on board MetOp with a lower spatial resolution of 25 km). A global map of the aeolian aerodynamic roughness length at 6 km resolution is derived, for arid and semi-arid regions. It shows very good consistency with the existing information on the properties of these surfaces. The dataset is available to the community, for use in atmospheric dust transport models.

Description: Anisotropy of small-scale stratospheric irregularities retrieved from scintillations of a double star α-Cru observed by GOMOS/ENVISAT Atmospheric Measurement Techniques, 5, 2713-2722, 2012 Author(s): V. Kan, V. F. Sofieva, and F. Dalaudier In this paper, we discuss estimating anisotropy of air density irregularities (ratio of characteristic horizontal and vertical scales) from satellite observations of bi-chromatic scintillations of a double star whose components are not resolved by the detector. The analysis is based on fitting experimental auto- and cross-spectra of scintillations by those computed using the 3-D spectral model of atmospheric irregularities consisting of anisotropic and isotropic components. Application of the developed method to the scintillation measurements of the double star α-Cru by GOMOS (Global Ozone Monitoring by Occultation of Stars) fast photometers results in estimates of anisotropy coefficient of ~15–20 at altitudes 30–38 km, as well as other parameters of atmospheric irregularities. The obtained estimates of the anisotropy coefficient correspond to small-scale irregularities, close to the buoyancy scale.

Description: Development of a new data-processing method for SKYNET sky radiometer observations Atmospheric Measurement Techniques, 5, 2723-2737, 2012 Author(s): M. Hashimoto, T. Nakajima, O. Dubovik, M. Campanelli, H. Che, P. Khatri, T. Takamura, and G. Pandithurai In order to reduce uncertainty in the estimation of Direct Aerosol Radiative Forcing (DARF), it is important to improve the estimation of the single scattering albedo (SSA). In this study, we propose a new data processing method to improve SSA retrievals for the SKYNET sky radiometer network, which is one of the growing number of networks of sun-sky photometers, such as NASA AERONET and others. There are several reports that SSA values from SKYNET have a bias compared to those from AERONET, which is regarded to be the most accurate due to its rigorous calibration routines and data quality and cloud screening algorithms. We investigated possible causes of errors in SSA that might explain the known biases through sensitivity experiments using a numerical model, and also using real data at the SKYNET sites at Pune (18.616° N/73.800° E) in India and Beijing (39.586° N/116.229° E) in China. Sensitivity experiments showed that an uncertainty of the order of ±0.03 in the SSA value can be caused by a possible error in the ground surface albedo or solid view angle assumed for each observation site. Another candidate for possible error in the SSA was found in cirrus contamination generated by imperfect cloud screening in the SKYNET data processing. Therefore, we developed a new data quality control method to get rid of low quality or cloud contamination data, and we applied this method to the real observation data at the Pune site in SKYNET. After applying this method to the observation data, we were able to screen out a large amount of cirrus-contaminated data and to reduce the deviation in the SSA value from that of AERONET. We then estimated DARF using data screened by our new method. The result showed that the method significantly reduced the difference of 5 W m −2 that existed between the SKYNET and AERONET values of DARF before screening. The present study also suggests the necessity of preparing suitable a priori information on the distribution of coarse particles ranging in radius between 10 μm and 30 μm for the analysis of heavily dust-laden atmospheric cases.

Description: Calibration and intercomparison of acetic acid measurements using proton-transfer-reaction mass spectrometry (PTR-MS) Atmospheric Measurement Techniques, 5, 2739-2750, 2012 Author(s): K. B. Haase, W. C. Keene, A. A. P. Pszenny, H. R. Mayne, R. W. Talbot, and B. C. Sive Acetic acid is one of the most abundant organic acids in the ambient atmosphere, with maximum mixing ratios reaching into the tens of parts per billion by volume (ppbv) range. The identities and associated magnitudes of the major sources and sinks for acetic acid are poorly characterized, due in part to the limitations of available measurement techniques. This paper demonstrates that, when properly calibrated, proton-transfer-reaction mass spectrometry (PTR-MS) can be a valuable technique for fast response, accurate quantification of acetic acid in ambient air. Three different PTR-MS configurations were calibrated at low ppbv mixing ratios using permeation tubes, which yielded calibration factors between 7.0 and 10.9 normalized counts per second per ppbv (ncps ppbv −1 ) at a drift tube field strength of 132 Townsend (Td). Detection limits ranged from 0.06 to 0.32 ppbv with dwell times of 5 s. These calibration factors showed negligible humidity dependence. Acetic acid was measured with PTR-MS on Appledore B Island, ME, during the International Consortium for Atmospheric Research on Transport and Transformation (ICARTT) campaign and validated based on acetic acid measured in parallel using tandem mist chambers coupled with ion chromatography (MC/IC). Mixing ratios ranged from a minimum of 0.075 ± 0.004 ppbv to 3.555 ± 0.171 ppbv, with a median mixing ratio of 0.530 ± 0.025 ppbv. An orthogonal least squares linear regression of paired data yielded a slope of 1.14 ± 0.06 (2σ), an intercept of 0.049 ± 0.020 (2σ) ppbv, and an R 2 of 0.78.

Description: Airborne observations of formic acid using a chemical ionization mass spectrometer Atmospheric Measurement Techniques, 5, 3029-3039, 2012 Author(s): M. Le Breton, M. R. McGillen, J. B. A. Muller, A. Bacak, D. E. Shallcross, P. Xiao, L. G. Huey, D. Tanner, H. Coe, and C. J. Percival The first airborne measurements of formic acid mixing ratios over the United Kingdom were measured on the FAAM BAe-146 research aircraft on 16 March 2010 with a chemical ionization mass spectrometer using I − reagent ions. The I − ionization scheme was able to measure formic acid mixing ratios at 1 Hz in the boundary layer. In-flight standard addition calibrations from a formic acid source were used to determine the instrument sensitivity of 35 ± 6 ion counts pptv −1 s −1 and a limit of detection of 25 pptv. Routine measurements were made through a scrubbed inlet to determine the instrumental background. Three plumes of formic acid were observed over the UK, originating from London, Humberside and Tyneside. The London plume had the highest formic acid mixing ratio throughout the flight, peaking at 358 pptv. No significant correlations of formic acid with NO x and ozone were found, but a positive correlation was observed between CO and HCOOH within the two plumes where coincident data were recorded. A trajectory model was employed to determine the sources of the plumes and compare modelled mixing ratios with measured values. The model underestimated formic acid concentrations by up to a factor of 2. This is explained by missing sources in the model, which were considered to be both primary emissions of formic acid of mainly anthropogenic origin and a lack of precursor emissions, such as isoprene, from biogenic sources, whose oxidation in situ would lead to formic acid formation.

Description: Applying spaceborne reflectivity measurements for calculation of the solar ultraviolet radiation at ground level Atmospheric Measurement Techniques, 5, 3041-3054, 2012 Author(s): P. N. den Outer, A. van Dijk, H. Slaper, A. V. Lindfors, H. De Backer, A. F. Bais, U. Feister, T. Koskela, and W. Josefsson Long-term analysis of cloud effects on ultraviolet (UV) radiation on the ground using spaceborne observations requires the use of instruments that have operated consecutively. The longest data record can be built from the reflectivity measurements produced by the instruments Total Ozone Mapping Spectrometers (TOMS) flown on Nimbus 7 from 1979 to 1992, TOMS on Earth Probe from 1996 to 2005, and the Ozone Monitoring Instrument (OMI) flown on EOS Aura since 2004. The reflectivity data produced by TOMS on Earth Probe is only included until 2002. A comparison is made with cloud effects inferred from ground-based pyranometer measurements at over 83 World Radiation Data Centre stations. Modelled UV irradiances utilizing the standard reflectivity are compared with measurements of UV irradiances at eight European low-elevation stations. The reflectivity data of the two TOMS instruments shows a consistent agreement, and the required corrections are of low percentage, i.e. 2–3%. In contrast, the reflectivity product of OMI requires correction of 7–10%, and a solar angle dependency therein is more pronounced. These corrections were inferred from a comparison with pyranometer data, and tested using the UV measurements. The average reduction of UV radiation due to clouds for all sites together indicates a small trend: a diminishing cloudiness, in line with ground-based UV observations. Uncorrected implementation of the reflectivity data would have indicated the opposite. An optimal area was established for reflectivity data for the calculation of daily sums of UV radiation. It measures approximately 1.25° in latitudinal direction for square-shaped areas overhead the ground-based UV stations. Such an area can be traversed within 5 to 7 h at the average wind speeds found for the West European continent.

Description: Present state of global wetland extent and wetland methane modelling: methodology of a model intercomparison project (WETCHIMP) Geoscientific Model Development Discussions, 5, 4071-4136, 2012 Author(s): R. Wania, J. R. Melton, E. L. Hodson, B. Poulter, B. Ringeval, R. Spahni, T. Bohn, C. A. Avis, G. Chen, A. V. Eliseev, P. O. Hopcroft, W. J. Riley, Z. M. Subin, H. Tian, V. Brovkin, P. M. van Bodegom, T. Kleinen, Z. C. Yu, J. S. Singarayer, S. Zürcher, D. P. Lettenmaier, D. J. Beerling, S. N. Denisov, C. Prigent, F. Papa, and J. O. Kaplan The Wetland and Wetland CH 4 Intercomparison of Models Project (WETCHIMP) was created to evaluate our present ability to simulate large-scale wetland characteristics and corresponding methane (CH 4 ) emissions. A multi-model comparison is essential to evaluate the key uncertainties in the mechanisms and parameters leading to methane emissions. Ten modelling groups joined WETCHIMP to run eight global and two regional models with a common experimental protocol using the same climate and atmospheric carbon dioxide (CO 2 ) forcing datasets. We reported the main conclusions from the intercomparison effort in a companion paper (Melton et al., 2012). Here we provide technical details for the six experiments, which included an equilibrium, a transient, and an optimized run plus three sensitivity experiments (temperature, precipitation, and atmospheric CO 2 concentration). The diversity of approaches used by the models is summarized through a series of conceptual figures, and is used to evaluate the wide range of wetland extents and CH 4 fluxes predicted by the models in the equilibrium run. We discuss relationships among the various approaches and patterns in consistencies of these model predictions. Within this group of models, there are three broad classes of methods used to estimate wetland extent: prescribed based on wetland distribution maps, prognostic relationships between hydrological states based on satellite observations, and explicit hydrological mass balances. A larger variety of approaches was used to estimate the net CH 4 fluxes from wetland systems. Even though modelling of wetland extents and CH 4 emissions has progressed significantly over recent decades, large uncertainties still exist when estimating CH 4 emissions: there is little consensus on model structure or complexity due to knowledge gaps, different aims of the models, and the range of temporal and spatial resolutions of the models.

Description: Improving the representation of secondary organic aerosol (SOA) in the MOZART-4 global chemical transport model Geoscientific Model Development Discussions, 5, 4187-4232, 2012 Author(s): A. Mahmud and K. C. Barsanti The secondary organic aerosol (SOA) module in the Model for Ozone and Related chemical Tracers, version 4 (MOZART-4) has been updated by replacing existing two-product (2p) parameters with those obtained from two-product volatility basis set (2p-VBS) fits, and by treating SOA formation from the following volatile organic compounds (VOCs): isoprene, propene and lumped alkenes. Strong seasonal and spatial variations in global SOA distributions were demonstrated, with significant differences in the predicted concentrations between the base-case and updated model versions. The base-case MOZART-4 predicted annual average SOA of 0.36 ± 0.50 μg m −3 in South America, 0.31 ± 0.38 μg m −3 in Indonesia, 0.09 ± 0.05 μg m −3 in the USA, and 0.12 ± 0.07 μg m −3 in Europe. Concentrations from the updated versions of the model showed a~marked increase in annual average SOA. Using the updated set of parameters alone (MZ4-v1) increased annual average SOA by ~8%, ~16%, ~56%, and ~108% from the base-case in South America, Indonesia, USA, and Europe, respectively. Treatment of additional parent VOCs (MZ4-v2) resulted in an even more dramatic increase of ~178–406% in annual average SOA for these regions over the base-case. The increases in predicted SOA concentrations further resulted in increases in corresponding SOA contributions to annual average total aerosol optical depth (AOD) by

Description: Modeling agriculture in the Community Land Model Geoscientific Model Development Discussions, 5, 4137-4185, 2012 Author(s): B. Drewniak, J. Song, J. Prell, V. R. Kotamarthi, and R. Jacob The potential impact of climate change on agriculture is uncertain. In addition, agriculture could influence above- and below-ground carbon storage. Development of models that represent agriculture is necessary to address these impacts. We have developed an approach to integrate agriculture representations for three crop types – maize, soybean, and spring wheat – into the coupled carbon-nitrogen version of the Community Land Model (CLM), to help address these questions. Here we present the new model, CLM-Crop, validated against observations from two AmeriFlux sites in the United States, planted with maize and soybean. Seasonal carbon fluxes compared well with field measurements. CLM-Crop yields were comparable with observations in some regions, although the generality of the crop model and its lack of technology and irrigation made direct comparison difficult. CLM-Crop was compared against the standard CLM3.5, which simulates crops as grass. The comparison showed improvement in gross primary productivity in regions where crops are the dominant vegetation cover. Crop yields and productivity were negatively correlated with temperature and positively correlated with precipitation. In case studies with the new crop model looking at impacts of residue management and planting date on crop yield, we found that increased residue returned to the litter pool increased crop yield, while reduced residue returns resulted in yield decreases. Using climate controls to signal planting date caused different responses in different crops. Maize and soybean had opposite reactions: when low temperature threshold resulted in early planting, maize responded with a loss of yield, but soybean yields increased. Our improvements in CLM demonstrate a new capability in the model – simulating agriculture in a realistic way, complete with fertilizer and residue management practices. Results are encouraging, with improved representation of human influences on the land surface and the potentially resulting climate impacts.

Description: A test of numerical instability and stiffness in the parametrizations of the ARPÉGE and ALADIN models Geoscientific Model Development Discussions, 5, 4233-4268, 2012 Author(s): M. Tudor Meteorological numerical weather prediction (NWP) models solve a system of partial differential equations in time and space. Semi-lagrangian advection scheme in the model dynamics allows for long time-steps. These longer time-steps can result in instabilities occurring in the model physics. A system of differential equations in which some solution components decay more rapidly than others is stiff. In this case it is stability rather than accuracy that restricts the time-step. The vertical diffusion parametrization can cause fast non-meteorological oscillations around the slowly evolving true solution (fibrillations). These are treated with an anti-fibrillation scheme. But small oscillations remain in an operational weather forecasts using ARPÉGE and ALADIN models. It is needed to test of the complete model formulation, as implemented in the operational forecast. In this paper, a simple test is designed. The test reveals if the formulation of particular physical parametrization is a stiff problem or potentially numerically unstable in combination with any other part of the model. When the test is applied to a stable scheme, the solution remains stable. But, applying the test to a potentially unstable scheme yields a solution with fibrillations of substantial amplitude. The parametrizations of a NWP model ARPÉGE were tested one by one to see which one may be the source of unstable model behaviour. The test has identified the stratiform precipitation scheme (a diagnostic Kessler type scheme) as a stiff problem, particularly the term that describes the evaporation of snow.

Description: Impact of the sea surface temperature forcing on hindcasts of Madden-Julian Oscillation events using the ECMWF model Ocean Science, 8, 1071-1084, 2012 Author(s): E. de Boisséson, M. A. Balmaseda, F. Vitart, and K. Mogensen This paper explores the sensitivity of hindcasts of the Madden Julian Oscillation (MJO) to the use of different sea surface temperture (SST) products as lower boundary conditions in the European Centre for Medium-range Weather Forecasts (ECMWF) atmospheric model. Three sets of monthly hindcast experiments are conducted, starting from initial conditions from the ERA interim reanalysis. First, as a reference, the atmosphere is forced by the SST used to produce ERA interim. In the second and third experiments, the SST is switched to the OSTIA (Operational Sea Surface Temperature and Sea-Ice Analysis) and the AVHRR-only (Advanced Very High Resolution Radiometer) reanalyses, respectively. Tests on the temporal resolution of the SST show that monthly fields are not optimal, while weekly and daily resolutions provide similar MJO scores. When using either OSTIA or AVHRR, the propagation of the MJO is degraded and the resulting scores are lower than in the reference experiment. Further experiments show that this loss of skill cannot be attributed to either the difference in mean state or temporal variability between the SST products. Additional diagnostics show that the phase relationship between either OSTIA or AVHRR SST and the MJO convection is distorted with respect to satellite observations and the ERA interim reanalysis. This distortion is expected to impact the MJO hindcasts, leading to a relative loss of forecast skill. A realistic representation of ocean–atmosphere interactions is thus needed for MJO hindcasts, but not all SST products – though accurate for other purposes – fulfill this requirement.

Description: Imbalance of energy and momentum source terms of the sea wave transfer equation for fully developed seas Ocean Science, 8, 1085-1098, 2012 Author(s): G. V. Caudal In the concept of full development, the sea wave spectrum is regarded as a nearly stationary solution of the wave transfer equation, where source and sink terms should be in balance with respect to both energy and momentum. Using a two-dimensional empirical sea wave spectral model at full development, this paper performs an assessment of the compatibility of the energy and momentum budgets of sea waves over the whole spectral range. Among the various combinations of model functions for wave breaking and wind source terms tested, not one is found to fulfill simultaneously the energy and momentum balance of the transfer equation. Based on experimental and theoretical grounds, wave breaking is known to contribute to frequency downshift of a narrow-banded wave spectrum when the modulational instability is combined with wave breaking. On those grounds, it is assumed that, in addition to dissipation, wave breaking produces a spectral energy flux directed toward low wavenumbers. I show that it is then possible to remove the energy and momentum budget inconsistency, and correspondingly the required strength of this spectral flux is estimated. Introducing such a downward spectral flux permits fulfilling both energy and momentum balance conditions. Meanwhile, the consistency between the transfer equation and empirical spectra, estimated by means of a cost function K, is either improved or slightly reduced, depending upon the wave breaking and wind source terms chosen. Other tests are performed in which it is further assumed that wave breaking would also be associated with azimuthal diffusion of the spectral energy. This would correspondingly reduce the required downward spectral flux by a factor of up to 5, although it would not be able to remove it entirely.

Description: MESMO 2: a mechanistic marine silica cycle and coupling to a simple terrestrial scheme Geoscientific Model Development Discussions, 5, 2999-3033, 2012 Author(s): K. Matsumoto, K. S. Tokos, A. Huston, and H. Joy-Warren Here we describe the second version of Minnesota Earth System Model for Ocean biogeochemistry (MESMO 2), an earth system model of intermediate complexity, which consists of a dynamical ocean, dynamic-thermodynamic sea ice, and energy moisture balanced atmosphere. The new version has more realistic land ice masks and is driven by seasonal winds. A major aim in version 2 is representing the marine silica cycle mechanistically in order to investigate climate-carbon feedbacks involving diatoms, a critically important class of phytoplankton in terms of carbon export production. This is achieved in part by including iron, on which phytoplankton uptake of silicic acid depends. Also, MESMO 2 is coupled to an existing terrestrial model, which allows for the exchange of carbon, water, and energy between land and the atmosphere. The coupled model, called MESMO 2E, is appropriate for more complete earth system simulations. The new version was calibrated with the goal of preserving reasonable interior ocean ventilation and various biological production rates in the ocean and land, while simulating key features of the marine silica cycle.

Description: High frequency fluctuations in the heat content of an ocean general circulation model Ocean Science, 8, 813-825, 2012 Author(s): A. M. Huerta-Casas and D. J. Webb The transport and storage of heat by the ocean is of crucial importance because of its effect on ocean dynamics and its impact on the atmosphere, climate and climate change. Unfortunately, limits to the amount of data that can be collected and stored mean that many experimental and modelling studies of the heat budget have to make use of mean datasets where the effects of short term fluctuations are lost. In this paper we investigate the magnitude of the resulting errors by making use of data from OCCAM, a high resolution global ocean model. The model carries out a proper heat balance every time step so any imbalances that are found in the analysis must result from the use of mean fields. The study concentrates on two areas of the ocean affecting the El Nino. The first is the region of tropical instability waves north of the Equator. The second is in the upwelling region along the Equator. It is shown that in both cases, processes with a period of less than five days can have a significant impact on the heat budget. Thus, analyses using data averaged over five days or more are likely to have significant errors. It is also shown that if a series of instantaneous values is available, reasonable estimates can be made of the size of the errors. In model studies, such values are available in the form of the datasets used to restart the model. In experimental studies they may be in the form of individual unaveraged observations.

Description: A 20-year reanalysis experiment in the Baltic Sea using three-dimensional variational (3DVAR) method Ocean Science, 8, 827-844, 2012 Author(s): W. Fu, J. She, and M. Dobrynin A 20-year retrospective reanalysis of the ocean state in the Baltic Sea is constructed by assimilating available historical temperature and salinity profiles into an operational numerical model with three-dimensional variational (3DVAR) method. To determine the accuracy of the reanalysis, the authors present a series of comparisons to independent observations on a monthly mean basis. In the reanalysis, temperature (T) and salinity (S) fit better with independent measurements than the free run at different depths. Overall, the mean biases of temperature and salinity for the 20 year period are reduced by 0.32 °C and 0.34 psu, respectively. Similarly, the mean root mean square error (RMSE) is decreased by 0.35 °C for temperature and 0.3 psu for salinity compared to the free run. The modeled sea surface temperature, which is mainly controlled by the weather forcing, shows the least improvements due to sparse in situ observations. Deep layers, on the other hand, witness significant and stable model error improvements. In particular, the salinity related to saline water intrusions into the Baltic Proper is largely improved in the reanalysis. The major inflow events such as in 1993 and 2003 are captured more accurately as the model salinity in the bottom layer is increased by 2–3 psu. Compared to independent sea level at 14 tide gauge stations, the correlation between model and observation is increased by 2%–5%, while the RMSE is generally reduced by 10 cm. It is found that the reduction of RMSE comes mainly from the reduction of mean bias. In addition, the changes in density induced by the assimilation of T/S contribute little to the barotropic transport in the shallow Danish Transition zone. The mixed layer depth exhibits strong seasonal variations in the Baltic Sea. The basin-averaged value is about 10 m in summer and 30 m in winter. By comparison, the assimilation induces a change of 20 m to the mixed layer depth in deep waters and wintertime, whereas small changes of about 2 m occur in summer and shallow waters. It is related to the strong heating in summer and the dominant role of the surface forcing in shallow water, which largely offset the effect of the assimilation.

Description: Corrigendum to "Greenhouse gas measurements over a 144 km open path in the Canary Islands" published in Atmos. Meas. Tech., 5, 2309–2319, 2012 Atmospheric Measurement Techniques, 5, 2349-2349, 2012 Author(s): J. S. A. Brooke, P. F. Bernath, G. Kirchengast, C. B. Thomas, J.-G. Wang, K. A. Tereszchuk, G. González Abad, R. J. Hargreaves, C. A. Beale, J. J. Harrison, S. Schweitzer, V. Proschek, P. A. Martin, V. L. Kasyutich, C. Gerbig, O. Kolle, and A. Loescher No abstract available.

Description: NO 2 observations over the western Pacific and Indian Ocean by MAX-DOAS on Kaiyo , a Japanese research vessel Atmospheric Measurement Techniques, 5, 2351-2360, 2012 Author(s): H. Takashima, H. Irie, Y. Kanaya, and F. Syamsudin Nitrogen dioxide (NO 2 ) profile retrievals were performed by ship-borne Multi-Axis Differential Optical Absorption Spectroscopy (MAX-DOAS) using a compact/low-power spectrometer on the Japanese research vessel Kaiyo during two ocean cruises around Japan and Japan–Bali (Indonesia)–Indian Ocean. DOAS analysis using a 425–450 nm fitting window revealed a clear land–ocean contrast in NO 2 differential slant column densities (DSCDs) but poor fitting results and negative values, especially at low elevation angles at low latitudes (

Description: A new method to diagnose the contribution of anthropogenic activities to temperature: temperature tagging Geoscientific Model Development Discussions, 5, 3183-3215, 2012 Author(s): V. Grewe This study presents a new methodology, called temperature tagging. It keeps track of the contributions of individual processes to temperature within a climate model simulation. As a first step and as a test bed a simple box climate model is regarded. The model consists of an atmosphere, which absorbs and emits radiation and of a surface, which reflects, absorbs and emits radiation. The tagging methodology is used to investigate the impact of the atmosphere on surface temperature. Four processes are investigated in more detail and their contribution to the surface temperature quantified: (i) shortwave influx and shortwave atmospheric absorption ("sw"), (ii) longwave atmospheric absorption due to non-CO 2 greenhouse gases ("nC"), (iii) due to a base case CO 2 concentration ("bC"), and (iv) due to an enhanced CO 2 concentration ("eC"). The differential equation for the temperature in the box climate model is decomposed into four equations for the tagged temperatures. This method is applied to investigate the contribution of longwave absorption to the surface temperature (greenhouse effect), which is calculated to be 68 K. This estimate contrasts an alternative calculation of the greenhouse effect of slightly more than 30 K based on the difference of the surface temperature with and without an atmosphere. The difference of the two estimates is due to a shortwave cooling effect and a reduced contribution of the shortwave to the total downward flux: The shortwave absorption of the atmosphere results in a reduced net shortwave flux at the surface of 192 W m −2 , leading to a cooling of the surface by 14 K. Introducing an atmosphere results in a downward longwave flux at the surface due to atmospheric absorption of 189 W m −2 , which roughly equals the net shortwave flux of 192 W m −2 . This longwave flux is a result of both, the radiation due to atmospheric temperatures and its longwave absorption. Hence the longwave absorption roughly accounts for 91 W m −2 out of a total of 381 W m −2 (roughly 25%) and therefore accounts for a temperature of 68 K. In a second experiment, the CO 2 concentration is doubled, which leads to an increase in surface temperature of 1.2 K, resulting from an temperature increase due to CO 2 of 1.9 K, due to non-CO 2 greenhouse gases of 0.6 K and a cooling of 1.3 K due to a reduced importance of the solar heating for the surface and atmospheric temperatures. These two experiments show the feasibility of temperature tagging and its potential as a diagnostic for climate simulations.

Description: Evaluation of Release-05 GRACE time-variable gravity coefficients over the ocean Ocean Science, 8, 859-868, 2012 Author(s): D. P. Chambers and J. A. Bonin The latest release of GRACE (Gravity Recovery and Climate Experiment) gravity field coefficients (Release-05, or RL05) are evaluated for ocean applications. Data have been processed using the current methodology for Release-04 (RL04) coefficients, and have been compared to output from two different ocean models. Results indicate that RL05 data from the three Science Data Centers – the Center for Space Research (CSR), GeoForschungsZentrum (GFZ), and Jet Propulsion Laboratory (JPL) – are more consistent among themselves than the previous RL04 data. Moreover, the variance of residuals with the output of an ocean model is 50–60% lower for RL05 data than for RL04 data. A more optimized destriping algorithm is also tested, which improves the results slightly. By comparing the GRACE maps with two different ocean models, we can better estimate the uncertainty in the RL05 maps. We find the standard error to be about 1 cm (equivalent water thickness) in the low- and mid-latitudes, and between 1.5 and 2 cm in the polar and subpolar oceans, which is comparable to estimated uncertainty for the output from the ocean models.

Description: Greenhouse gas measurements over a 144 km open path in the Canary Islands Atmospheric Measurement Techniques, 5, 2309-2319, 2012 Author(s): J. S. A. Brooke, P. F. Bernath, G. Kirchengast, C. B. Thomas, J.-G. Wang, K. A. Tereszchuk, G. González Abad, R. J. Hargreaves, C. A. Beale, J. J. Harrison, S. Schweitzer, V. Proschek, P. A. Martin, V. L. Kasyutich, C. Gerbig, O. Kolle, and A. Loescher A new technique for the satellite remote sensing of greenhouse gases in the atmosphere via the absorption of short-wave infrared laser signals transmitted between counter-rotating satellites in low Earth orbit has recently been proposed; this would enable the acquisition of a long-term, stable, global set of altitude-resolved concentration measurements. We present the first ground-based experimental demonstration of this new infrared-laser occultation method, in which the atmospheric absorption of CO 2 near 2.1 μm was measured over a ~144 km path length between two peaks in the Canary Islands (at an altitude of ~2.4 km), using relatively low power diode lasers (~4 to 10 mW). The retrieved CO 2 volume mixing ratio of 400 ppm (±15 ppm) is consistent within experimental uncertainty with simultaneously recorded in situ validation measurements. We conclude that the new method has a sound basis for monitoring CO 2 in the free atmosphere; other greenhouse gases such as methane, nitrous oxide and water vapour can be monitored in the same way.

Description: An empirical model of optical and radiative characteristics of the tropospheric aerosol over West Siberia in summer Atmospheric Measurement Techniques, 5, 1513-1527, 2012 Author(s): M. V. Panchenko, T. B. Zhuravleva, S. A. Terpugova, V. V. Polkin, and V. S. Kozlov An empirical model of the vertical profiles of aerosol optical characteristics is described. This model was developed based on data acquired from multi-year airborne sensing of optical and microphysical characteristics of the tropospheric aerosol over West Siberia. The main initial characteristics for the creation of the model were measurement data of the vertical profiles of the aerosol angular scattering coefficients in the visible wavelength range, particle size distribution functions and mass concentrations of black carbon (BC). The proposed model allows us to retrieve the aerosol optical and radiative characteristics in the visible and near-IR wavelength range, using the season, air mass type and time of day as input parameters. The columnar single scattering albedo and asymmetry factor of the aerosol scattering phase function, calculated using the average vertical profiles, are in good agreement with data from the AERONET station located in Tomsk. For solar radiative flux calculations, this empirical model has been tested for typical summer conditions. The available experimental database obtained for the regional features of West Siberia and the model developed on this basis are shown to be sufficient for performing these calculations.

Description: Calibration of an all-sky camera for obtaining sky radiance at three wavelengths Atmospheric Measurement Techniques, 5, 2013-2024, 2012 Author(s): R. Román, M. Antón, A. Cazorla, A. de Miguel, F. J. Olmo, J. Bilbao, and L. Alados-Arboledas This paper proposes a method to obtain spectral sky radiances, at three wavelengths (464, 534 and 626 nm), from hemispherical sky images. Images are registered with the All-Sky Imager installed at the Andalusian Center for Environmental Research (CEAMA) in Granada (Spain). The methodology followed in this work for the absolute calibration in radiance of this instrument is based on the comparison of its output measurements with modelled sky radiances derived from the LibRadtran/UVSPEC radiative transfer code under cloud-free conditions. Previously, in order to check the goodness of the simulated radiances, these are compared with experimental values recorded by a CIMEL sunphotometer. In general, modelled radiances are in agreement with experimental data, showing mean differences lower than 20% except for the pixels located next to the Sun position that show larger errors. The relationship between the output signal of the All-Sky Imager and the modelled sky radiances provides a calibration matrix for each image. The variability of the matrix coefficients is analyzed, showing no significant changes along a period of 5 months. Therefore, a unique calibration matrix per channel is obtained for all selected images (a total of 705 images per channel). Camera radiances are compared with CIMEL radiances, finding mean absolute differences between 2% and 15% except for pixels near to the Sun and high scattering angles. We apply these calibration matrices to three images in order to study the sky radiance distributions for three different sky conditions: cloudless, overcast and partially cloudy. Horizon brightening under cloudless conditions has been observed together with the enhancement effect of individual clouds on sky radiance.

Description: Aircraft measurements of carbon dioxide and methane for the calibration of ground-based high-resolution Fourier Transform Spectrometers and a comparison to GOSAT data measured over Tsukuba and Moshiri Atmospheric Measurement Techniques, 5, 2003-2012, 2012 Author(s): T. Tanaka, Y. Miyamoto, I. Morino, T. Machida, T. Nagahama, Y. Sawa, H. Matsueda, D. Wunch, S. Kawakami, and O. Uchino Aircraft measurements of carbon dioxide and methane over Tsukuba (36.05° N, 140.12° E) (February 2010) and Moshiri (44.36° N, 142.26° E) (August 2009) were made to calibrate ground-based high-resolution Fourier Transform Spectrometers (g-b FTSs) and to compare with the Greenhouse gases Observing SATellite (GOSAT). The aircraft measurements over Tsukuba in February 2010 were successful in synchronizing with both the g-b FTS and GOSAT for the first time. Airborne in situ and flask-sampling instruments were mounted on the aircraft, and measurements were carried out between altitudes of 0.5 and 7 km to obtain vertical profiles of carbon dioxide (CO 2 ), methane (CH 4 ), and other gaseous species. By comparing the g-b FTS measurements with the airborne measurements, the column-averaged dry air mole fractions of CO 2 (X CO 2 ) and CH 4 (X CH 4 ) retrieved from the g-b FTS measurements at Tsukuba were biased low by 0.33 ± 0.11% for X CO 2 and 0.69 ± 0.29% for X CH 4 . The g-b FTS values at Moshiri were biased low by 1.24% for X CO 2 and 2.11% for X CH 4 . The GOSAT data show biases that are 3.1% ± 1.7% lower for X CO 2 and 2.5% ± 0.8% lower for X CH 4 than the aircraft measurements obtained over Tsukuba.

Description: Quantitative bias estimates for tropospheric NO 2 columns retrieved from SCIAMACHY, OMI, and GOME-2 using a common standard for East Asia Atmospheric Measurement Techniques, 5, 2403-2411, 2012 Author(s): H. Irie, K. F. Boersma, Y. Kanaya, H. Takashima, X. Pan, and Z. F. Wang For the intercomparison of tropospheric nitrogen dioxide (NO 2 ) vertical column density (VCD) data from three different satellite sensors (SCIAMACHY, OMI, and GOME-2), we use a common standard to quantitatively evaluate the biases for the respective data sets. As the standard, a regression analysis using a single set of collocated ground-based Multi-Axis Differential Optical Absorption Spectroscopy (MAX-DOAS) observations at several sites in Japan and China from 2006–2011 is adopted. Examinations of various spatial coincidence criteria indicates that the slope of the regression line can be influenced by the spatial distribution of NO 2 over the area considered. While the slope varies systematically with the distance between the MAX-DOAS and satellite observation points around Tokyo in Japan, such a systematic dependence is not clearly seen and correlation coefficients are generally higher in comparisons at sites in China. On the basis of these results, we focus mainly on comparisons over China and estimate the biases in SCIAMACHY, OMI, and GOME-2 data (TM4NO2A and DOMINO version 2 products) against the MAX-DOAS observations to be −5 ± 14%, −10 ± 14%, and +1 ± 14%, respectively, which are all small and insignificant. We suggest that these small biases now allow for analyses combining these satellite data for air quality studies, which are more systematic and quantitative than previously possible.

Description: The Mediterranean Ocean Colour Observing System – system development and product validation Ocean Science, 8, 869-883, 2012 Author(s): G. Volpe, S. Colella, V. Forneris, C. Tronconi, and R. Santoleri This paper presents the Mediterranean Ocean Colour Observing System in the framework of the growing demand of near real-time data emerging within the operational oceanography international context. The main issues related to the satellite operational oceanography are tied to the following: (1) the near real-time ability to track data flow uncertainty sources; (2) in case of failure, to provide backup solutions to end-users; and (3) to scientifically assess the product quality. We describe the major scientific and technological steps made to develop, maintain and improve the operational system and its products. A method for assessing the near real-time product quality is developed and its limitation discussed. Main results are concerned with the degradation, starting from mid-2010, of the MODIS Aqua channel at 443 nm with its successive recovery thanks to the new calibration scheme implemented in the recently released SeaDAS version 6.4. The product validation analysis highlights that SeaWiFS chlorophyll product over the Mediterranean Sea is the best performing in comparison with those of MODIS and MERIS. Despite their general good agreement with in situ observations, MODIS- and MERIS-derived chlorophyll present a slight and systematic underestimation of the in situ counter part. The most relevant implications induced by these results are discussed from an operational point of view.

Description: Implementation of the Fast-JX Photolysis scheme into the UKCA component of the MetUM chemistry climate model Geoscientific Model Development Discussions, 5, 3217-3260, 2012 Author(s): P. J. Telford, N. L. Abraham, A. T. Archibald, P. Braesicke, M. Dalvi, O. Morgenstern, F. M. O'Connor, N. A. D. Richards, and J. A. Pyle Atmospheric chemistry is driven by photolytic reactions, making their modelling a crucial component of atmospheric models. We describe the implementation and validation of Fast-JX, a state of the art model of interactive photolysis, into the MetUM chemistry climate model. This allows for interactive photolysis frequencies to be calculated in the troposphere and augments the calculation of the frequencies in the stratosphere by accounting for clouds and aerosols in addition to ozone. In order to demonstrate the effectiveness of this new photolysis scheme we employ new methods of validating the model, including techniques for sampling the model to compare to flight track and satellite data.

Description: COSTRICE – three model online coupling using OASIS: problems and solutions Geoscientific Model Development Discussions, 5, 3261-3310, 2012 Author(s): H. T. M. Ho, B. Rockel, H. Kapitza, B. Geyer, and E. Meyer The coupled system COSTRICE is developed for the first time in order to reproduce the interactions and feedbacks between atmosphere, ocean and sea-ice in a two-way online coupled model system containing three component models for regional climate simulations over Baltic Sea and North Sea regions. The regional climate model CCLM 1 is coupled to the regional ocean model TRIMNP 1 and the sea ice model CICE 1 via the coupler OASIS3. In this study, CCLM is setup with a horizontal grid mesh size of 50 km and 32 vertical atmosphere layers and driven by the 6-h ERA-interim reanalysis data as initial and boundary conditions. TRIMNP is setup with a horizontal grid mesh size of 12.8 km and 50 vertical ocean levels. CICE calculates ice in 5 categories and runs with the same horizontal resolution as TRIMNP but only over the Baltic Sea and the Kattegat Bay of the North Sea. In a two-way online coupling process, CCLM is linked to TRIMNP through sea surface temperature (SST) as lower boundary condition every 3 h and TRIMNP is driven by 1-h atmospheric state variables and fluxes of CCLM. The data exchange processes between TRIMNP and CICE as well as from CCLM to CICE take place with an interval of 3 h. The coupled model is applied in a study for climate simulations over Baltic Sea and North Sea regions in 1997. The coupled system is set up to run in parallel on the super computing system IBM-power 6 at the German Climate Computing Center (DKRZ). 1 See Table A1.

Description: A generalized tagging method Geoscientific Model Development Discussions, 5, 3311-3324, 2012 Author(s): V. Grewe The understanding of causes of changes in climate-chemistry simulations is an important, but often challenging task. In atmospheric chemistry, one approach is to tag species according to their origin (e.g. emission categories) and to inherit these tags to other species during subsequent reactions. This concept was recently employed to calculate the contribution of atmospheric processes to temperature. Here a new concept for tagging any state variable is presented. This generalized tagging method results from a sensitivity analysis of the forcing terms of the right hand side of the governing differential equations. In a couple of examples, the consistency with previous approaches is shown. Since the method is based on a ratio describing relative sensitivities, singularities occur where the method is not applicable. For some applications, like in atmospheric chemistry, these singularities can easily be removed. However, one theoretical example is given, where this method is not applicable at all.

Description: Combined wind measurements by two different lidar instruments in the Arctic middle atmosphere Atmospheric Measurement Techniques, 5, 2433-2445, 2012 Author(s): J. Hildebrand, G. Baumgarten, J. Fiedler, U.-P. Hoppe, B. Kaifler, F.-J. Lübken, and B. P. Williams During a joint campaign in January 2009, the Rayleigh/Mie/Raman (RMR) lidar and the sodium lidar at the ALOMAR Observatory (69° N, 16° E) in Northern Norway were operated simultaneously for more than 40 h, collecting data for wind measurements in the middle atmosphere from 30 up to 110 km altitude. As both lidars share the same receiving telescopes, the upper altitude range of the RMR lidar and the lower altitude range of the sodium lidar overlap in the altitude region of ≈80–85 km. For this overlap region we are thus able to present the first simultaneous wind measurements derived from two different lidar instruments. The comparison of winds derived by RMR and sodium lidar is excellent for long integration times of 10 h as well as shorter ones of 1 h. Combination of data from both lidars allows identifying wavy structures between 30 and 110 km altitude, whose amplitudes increase with height. We have also performed vertical wind measurements and measurements of the same horizontal wind component using two independent lasers and telescopes of the RMR lidar and show how to use this data to calibrate and validate the wind retrieval. For the latter configuration we found a good agreement of the results but also identified inhomogeneities in the horizontal wind at about 55 km altitude of up to 20 ms −1 for an integration time of nearly 4 h. Such small-scale inhomogeneities in the horizontal wind field are an essential challenge when comparing data from different instruments.

Description: Technical Note: Improving computational efficiency in large linear inverse problems: an example from carbon dioxide flux estimation Geoscientific Model Development Discussions, 5, 3325-3342, 2012 Author(s): V. Yadav and A. M. Michalak Addressing a variety of questions within Earth science disciplines entails the inference of the spatio-temporal distribution of parameters of interest based on observations of related quantities. Such estimation problems often represent inverse problems that are formulated as linear optimization problems. Computational limitations arise when the number of observations and/or the size of the discretized state space become large, especially if the inverse problem is formulated in a probabilistic framework and therefore aims to assess the uncertainty associated with the estimates. This work proposes two approaches to lower the computational costs and memory requirements for large linear space-time inverse problems, taking the Bayesian approach for estimating carbon dioxide (CO 2 ) emissions and uptake (a.k.a. fluxes) as a prototypical example. The first algorithm can be used to efficiently multiply two matrices, as long as one can be expressed as a Kronecker product of two smaller matrices, a condition that is typical when multiplying a sensitivity matrix by a covariance matrix in the solution of inverse problems. The second algorithm can be used to compute a posteriori uncertainties directly at aggregated spatio-temporal scales, which are the scales of most interest in many inverse problems. Both algorithms have significantly lower memory requirements and computational complexity relative to direct computation of the same quantities (O( n 2.5 ) vs. O( n 3 )). For an examined benchmark problem, the two algorithms yielded a three and six order of magnitude increase in computational efficiency, respectively, relative to direct computation of the same quantities. Sample computer code is provided for assessing the computational and memory efficiency of the proposed algorithms for matrices of different dimensions.

Description: Quality assessment concept of the World Data Center for Climate and its application to CMIP5 data Geoscientific Model Development Discussions, 5, 781-802, 2012 Author(s): M. Stockhause, H. Höck, F. Toussaint, and M. Lautenschlager The preservation of data in a high state of quality and suitable for interdisciplinary use is one of the most pressing and challenging current issues in long-term archiving. For high volume data such as climate model data, the data and data replica are no longer stored centrally but distributed over several local data repositories, e.g. the data of the Climate Model Intercomparison Project No. 5 (CMIP5). The most important part of the data is to be published as DOI according to the World Data Center for Climate's (WDCC) application of the DataCite regulations. The integrated part of WDCC's data publication process, the data quality assessment, was adapted to the requirements of a federated data infrastructure. A concept of a distributed and federated quality assessment procedure was developed, in which the work load and responsibility for quality control is shared between the three primary CMIP5 data centers: Program for Climate Model Diagnosis and Intercomparison (PCMDI), British Atmospheric Data Centre (BADC), and WDCC. This distributed quality control concept, its pilot implementation for CMIP5, and first experiences are presented.

Description: Global distributions of C 2 H 6 , C 2 H 2 , HCN, and PAN retrieved from MIPAS reduced spectral resolution measurements Atmospheric Measurement Techniques, 5, 723-734, 2012 Author(s): A. Wiegele, N. Glatthor, M. Höpfner, U. Grabowski, S. Kellmann, A. Linden, G. Stiller, and T. von Clarmann Vertical profiles of mixing ratios of C 2 H 6 , C 2 H 2 , HCN, and PAN were retrieved from MIPAS reduced spectral resolution nominal mode limb emission measurements. The retrieval strategy follows that of the analysis of MIPAS high resolution measurements, with occasional adjustments to cope with the reduced spectral resolution under which MIPAS is operated since 2005. MIPAS measurements from January 2005 to January 2010 have been analyzed with special emphasis on October 2007. Largest mixing ratios are found in the troposphere, and reach 1.2 ppbv for C 2 H 6 , 1 ppbv for HCN, 600 pptv for PAN, and 450 pptv for C 2 H 2 . The estimated precisions in case of significantly enhanced mixing ratios (including measurement noise and propagation of uncertain parameters randomly varying in the time domain) and altitude resolution are typically 10%, 3–4.5 km for C 2 H 6 , 15%, 4–6 km for HCN, 6%, 2.5–3.5 km for PAN, and 7%, 2.5–4 km for C 2 H 2 .

Description: DOAS measurements of NO 2 from an ultralight aircraft during the Earth Challenge expedition Atmospheric Measurement Techniques, 5, 2057-2068, 2012 Author(s): A. Merlaud, M. Van Roozendael, J. van Gent, C. Fayt, J. Maes, X. Toledo-Fuentes, O. Ronveaux, and M. De Mazière We report on airborne Differential Optical Absorption Spectroscopy (DOAS) measurements of NO 2 tropospheric columns above South Asia, the Arabic peninsula, North Africa, and Italy in November and December 2009. The DOAS instrument was installed on an ultralight aircraft involved in the Earth Challenge project, an expedition of seven pilots flying on four ultralight aircraft between Australia and Belgium. The instrument recorded spectra in limb geometry with a large field of view, a set-up which provides a high sensitivity to the boundary layer NO 2 while minimizing the uncertainties related to the attitude variations. We compare our measurements with OMI (Ozone Monitoring Instrument) and GOME-2 (Global Ozone Monitoring Experiment 2) tropospheric NO 2 products when the latter are available. Above Rajasthan and the Po Valley, two areas where the NO 2 field is homogeneous, data sets agree very well. Our measurements in these areas are 0.1 ± 0.1 to 3 ± 1 × 10 15 molec cm −2 and 2.6 ± 0.8 × 10 16 molec cm −2 , respectively. Flying downwind of Riyadh, our NO 2 measurements show the structure of the megacity's exhaust plume with a higher spatial resolution than OMI. Moreover, our measurements are larger (up to 40%) than those seen by satellites. We also derived tropospheric columns when no satellite data were available if it was possible to get information on the visibility from satellite measurements of aerosol optical thickness. This experiment also provides a confirmation for the recent finding of a soil signature above desert.

Description: SPITFIRE-2: an improved fire module for Dynamic Global Vegetation Models Geoscientific Model Development Discussions, 5, 2347-2443, 2012 Author(s): M. Pfeiffer and J. O. Kaplan Fire is the primary disturbance factor in many terrestrial ecosystems. Wildfire alters vegetation structure and composition, affects carbon storage and biogeochemical cycling, and results in the release of climatically relevant trace gases, including CO 2 , CO, CH 4 , NO x , and aerosols. Assessing the impacts of global wildfire on centennial to multi-millennial timescales requires the linkage of process-based fire modeling with vegetation modeling using Dynamic Global Vegetation Models (DGVMs). Here we present a new fire module, SPITFIRE-2, and an update to the LPJ-DGVM that includes major improvements to the way in which fire occurrence, behavior, and the effect of fire on vegetation is simulated. The new fire module includes explicit calculation of natural ignitions, the representation of multi-day burning and coalescence of fires and the calculation of rates of spread in different vegetation types, as well as a simple scheme to model crown fires. We describe a new representation of anthropogenic biomass burning under preindustrial conditions that distinguishes the way in which the relationship between humans and fire are different between hunter-gatherers, obligate pastoralists, and farmers. Where and when available, we evaluate our model simulations against remote-sensing based estimates of burned area. While wildfire in much of the modern world is largely influenced by anthropogenic suppression and ignitions, in those parts of the world where natural fire is still the dominant process, e.g. in remote areas of the boreal forest, our results demonstrate a significant improvement in simulated burned area over previous models. With its unique properties of being able to simulate preindustrial fire, the new module we present here is particularly well suited for the investigation of climate-human-fire relationships on multi-millennial timescales.

Description: The Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP): overview and description of models, simulations and climate diagnostics Geoscientific Model Development Discussions, 5, 2445-2502, 2012 Author(s): J.-F. Lamarque, D. T. Shindell, B. Josse, P. J. Young, I. Cionni, V. Eyring, D. Bergmann, P. Cameron-Smith, W. J. Collins, R. Doherty, S. Dalsoren, G. Faluvegi, G. Folberth, S. J. Ghan, L. W. Horowitz, Y. H. Lee, I. A. MacKenzie, T. Nagashima, V. Naik, D. Plummer, M. Righi, S. Rumbold, M. Schulz, R. B. Skeie, D. S. Stevenson, S. Strode, K. Sudo, S. Szopa, A. Voulgarakis, and G. Zeng The Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP) consists of a series of timeslice experiments targeting the long-term changes in atmospheric composition between 1850 and 2100, with the goal of documenting radiative forcing and the associated composition changes. Here we introduce the various simulations performed under ACCMIP and the associated model output. The ACCMIP models have a wide range of horizontal and vertical resolutions, vertical extent, chemistry schemes and interaction with radiation and clouds. While anthropogenic and biomass burning emissions were specified for all time slices in the ACCMIP protocol, it is found that the natural emissions lead to a significant range in emissions, mostly for ozone precursors. The analysis of selected present-day climate diagnostics (precipitation, temperature, specific humidity and zonal wind) reveals biases consistent with state-of-the-art climate models. The model-to-model comparison of changes in temperature, specific humidity and zonal wind between 1850 and 2000 and between 2000 and 2100 indicates mostly consistent results, but with outliers different enough to possibly affect their representation of climate impact on chemistry.

Description: Evaluating nighttime CALIOP 0.532 μm aerosol optical depth and extinction coefficient retrievals Atmospheric Measurement Techniques, 5, 2143-2160, 2012 Author(s): J. R. Campbell, J. L. Tackett, J. S. Reid, J. Zhang, C. A. Curtis, E. J. Hyer, W. R. Sessions, D. L. Westphal, J. M. Prospero, E. J. Welton, A. H. Omar, M. A. Vaughan, and D. M. Winker NASA Cloud Aerosol Lidar with Orthogonal Polarization (CALIOP) Version 3.01 5-km nighttime 0.532 μm aerosol optical depth (AOD) datasets from 2007 are screened, averaged and evaluated at 1° × 1° resolution versus corresponding/co-incident 0.550 μm AOD derived using the US Navy Aerosol Analysis and Prediction System (NAAPS), featuring two-dimensional variational assimilation of quality-assured NASA Moderate Resolution Imaging Spectroradiometer (MODIS) and Multi-angle Imaging Spectroradiometer (MISR) AOD. In the absence of sunlight, since passive radiometric AOD retrievals rely overwhelmingly on scattered radiances, the model represents one of the few practical global estimates available from which to attempt such a validation. Daytime comparisons, though, provide useful context. Regional-mean CALIOP vertical profiles of night/day 0.532 μm extinction coefficient are compared with 0.523/0.532 μm ground-based lidar measurements to investigate representativeness and diurnal variability. In this analysis, mean nighttime CALIOP AOD are mostly lower than daytime (0.121 vs. 0.126 for all aggregated data points, and 0.099 vs. 0.102 when averaged globally per normalised 1° × 1° bin), though the relationship is reversed over land and coastal regions when the data are averaged per normalised bin (0.134/0.108 vs. 0140/0.112, respectively). Offsets assessed within single bins alone approach ±20%. CALIOP AOD, both day and night, are higher than NAAPS over land (0.137 vs. 0.124) and equal over water (0.082 vs. 0.083) when averaged globally per normalised bin. However, for all data points inclusive, NAAPS exceeds CALIOP over land, coast and ocean, both day and night. Again, differences assessed within single bins approach 50% in extreme cases. Correlation between CALIOP and NAAPS AOD is comparable during both day and night. Higher correlation is found nearest the equator, both as a function of sample size and relative signal magnitudes inherent at these latitudes. Root mean square deviation between CALIOP and NAAPS varies between 0.1 and 0.3 globally during both day/night. Averaging of CALIOP along-track AOD data points within a single NAAPS grid bin improves correlation and RMSD, though day/night and land/ocean biases persist and are believed systematic. Vertical profiles of extinction coefficient derived in the Caribbean compare well with ground-based lidar observations, though potentially anomalous selection of a priori lidar ratios for CALIOP retrievals is likely inducing some discrepancies. Mean effective aerosol layer top heights are stable between day and night, indicating consistent layer-identification diurnally, which is noteworthy considering the potential limiting effects of ambient solar noise during day.

Description: The scientific basis for a satellite mission to retrieve CCN concentrations and their impacts on convective clouds Atmospheric Measurement Techniques, 5, 2039-2055, 2012 Author(s): D. Rosenfeld, E. Williams, M. O. Andreae, E. Freud, U. Pöschl, and N. O. Rennó The cloud-mediated aerosol radiative forcing is widely recognized as the main source of uncertainty in our knowledge of the anthropogenic forcing on climate. The current challenges for improving our understanding are (1) global measurements of cloud condensation nuclei (CCN) in the cloudy boundary layer from space, and (2) disentangling the effects of aerosols from the thermodynamic and meteorological effects on the clouds. Here, we present a new conceptual framework to help us overcome these two challenges, using relatively simple passive satellite measurements in the visible and infared (IR). The idea is to use the clouds themselves as natural CCN chambers by retrieving simultaneously the number of activated aerosols at cloud base, N a , and the cloud base updraft speed. The N a is obtained by analyzing the distribution of cloud drop effective radius in convective elements as a function of distance above cloud base. The cloud base updraft velocities are estimated by double stereoscopic viewing and tracking of the evolution of cloud surface features just above cloud base. In order to resolve the vertical dimension of the clouds, the field of view will be 100 m for the microphysical retrievals, and 50 m for the stereoscopic measurements. The viewing geometry will be eastward and 30 degrees off nadir, with the Sun in the back at 30 degrees off zenith westward, requiring a Sun-synchronous orbit at 14 LST. Measuring simultaneously the thermodynamic environment, the vertical motions of the clouds, their microstructure and the CCN concentration will allow separating the dynamics from the CCN effects. This concept is being applied in the proposed satellite mission named Clouds, Hazards and Aerosols Survey for Earth Researchers (CHASER).

Description: Evaluation of continuous water vapor δD and δ 18 O measurements by off-axis integrated cavity output spectroscopy Atmospheric Measurement Techniques, 5, 2069-2080, 2012 Author(s): N. Kurita, B. D. Newman, L. J. Araguas-Araguas, and P. Aggarwal Recent commercially available laser spectroscopy systems enabled us to continuously and reliably measure the δD and δ 18 O of atmospheric water vapor. The use of this new technology is becoming popular because of its advantages over the conventional approach based on cold trap collection. These advantages include much higher temporal resolution/continuous monitoring and the ability to make direct measurements of both isotopes in the field. Here, we evaluate the accuracy and precision of the laser based water vapor isotope instrument through a comparison of measurements with those found using the conventional cold trap method. A commercially available water vapor isotope analyzer (WVIA) with the vaporization system of a liquid water standard (Water Vapor Isotope Standard Source, WVISS) from Los Gatos Research (LGR) Inc. was used for this study. We found that the WVIA instrument can provide accurate results if (1) correction is applied for time-dependent isotope drift, (2) normalization to the VSMOW/SLAP scale is implemented, and (3) the water vapor concentration dependence of the isotopic ratio is also corrected. In addition, since the isotopic value of water vapor generated by the WVISS is also dependent on the concentration of water vapor, this effect must be considered to determine the true water vapor concentration effect on the resulting isotope measurement. To test our calibration procedure, continuous water vapor isotope measurements using both a laser instrument and a cold trap system were carried out at the IAEA Isotope Hydrology Laboratory in Vienna from August to December 2011. The calibrated isotopic values measured using the WVIA agree well with those obtained via the cold trap method. The standard deviation of the isotopic difference between both methods is about 1.4‰ for δD and 0.28‰ for δ 18 O. This precision allowed us to obtain reliable values for d -excess. The day-to-day variation of d -excess measured by WVIA also agrees well with that found using the cold trap method. These results demonstrate that a coupled system, using commercially available WVIA and WVISS instruments can provide continuous and accurate isotope data, with results achieved similar to those obtained using the conventional method, but with drastically improved temporal resolution.

Description: Setup of the PMIP3 paleoclimate experiments conducted using an Earth System Model, MIROC-ESM Geoscientific Model Development Discussions, 5, 2527-2569, 2012 Author(s): T. Sueyoshi, R. Ohgaito, A. Yamamoto, M. O. Chikamoto, T. Hajima, H. Okajima, M. Yoshimori, M. Abe, R. O'ishi, F. Saito, S. Watanabe, M. Kawamiya, and A. Abe-Ouchi The importance of climate model evaluation using paleoclimate simulations for better future climate projections has been recognized by the Intergovernmental Panel on Climate Change. In recent years, Earth System Models (ESMs) were developed to investigate carbon-cycle climate feedback, as well as to project the future climate. Paleoclimate events, especially those associated with the variations in atmospheric CO 2 level or land vegetation, provide suitable benchmarks to evaluate ESMs. Here we present implementations of the paleoclimate experiments proposed by the Coupled Model Intercomparison Project phase 5/Paleoclimate Modelling Intercomparison Project phase 3 (CMIP5/PMIP3) using an Earth System Model, MIROC-ESM. In this paper, experimental settings and procedures of the mid-Holocene, the Last Glacial Maximum, and the Last Millennium experiments are explained. The first two experiments are time slice experiments and the last one is a transient experiment. The complexity of the model requires various steps to correctly configure the experiments. Several basic outputs are also shown.

Description: A disjunct eddy accumulation system for the measurement of BVOC fluxes: instrument characterizations and field deployment Atmospheric Measurement Techniques, 5, 2115-2132, 2012 Author(s): G. D. Edwards, D. K. Martins, T. Starn, K. Pratt, and P. B. Shepson Biological volatile organic compounds (BVOCs), such as isoprene and monoterpenes, are emitted in large amounts from forests. Quantification of the flux of BVOCs is critical in the evaluation of the impact of these compounds on the concentrations of atmospheric oxidants and on the production of secondary organic aerosol. A disjunct eddy accumulation (DEA) sampler system was constructed for the measurement of speciated BVOC fluxes. Unlike traditional eddy covariance (EC), the relatively new technique of disjunct sampling differs by taking short, discrete samples that allow for slower sampling frequencies. Disjunct sample airflow is directed into cartridges containing sorbent materials at sampling rates proportional to the magnitude of the vertical wind. Compounds accumulated on the cartridges are then quantified by thermal desorption and gas chromatography. Herein, we describe our initial tests to evaluate the disjunct sampler including the application of vertical wind measurements to create optimized sampling thresholds. Measurements of BVOC fluxes obtained from DEA during its deployment above a mixed hardwood forest at the University of Michigan Biological Station (Pellston, MI) during the 2009 CABINEX field campaign are reported. Daytime (09:00 a.m. to 05:00 p.m. LT) isoprene fluxes, when averaged over the footprint of the tower, were 1.31 mg m −2 h −1 which are comparable to previous flux measurements at this location. Speciated monoterpene fluxes are some of the first to be reported from this site. Daytime averages were 26.7 μg m −2 h −1 for α-pinene and 10.6 μg m −2 h −1 for β-pinene. These measured concentrations and fluxes were compared to the output of an atmospheric chemistry model, and were found to be consistent with our knowledge of the variables that control BVOCs fluxes at this site.

Description: Validation of six years of SCIAMACHY carbon monoxide observations using MOZAIC CO profile measurements Atmospheric Measurement Techniques, 5, 2133-2142, 2012 Author(s): A. T. J. de Laat, R. Dijkstra, H. Schrijver, P. Nédélec, and I. Aben This paper presents a validation study of SCanning Imaging Absorption SpectroMeter for Atmospheric CartograpHY (SCIAMACHY) carbon monoxide (CO) total column measurements from the Iterative Maximum Likelihood Method (IMLM) algorithm using vertically integrated profile aircraft measurements obtained within the MOZAIC project for the six year time period of 2003–2008. Overall we find a good agreement between SCIAMACHY and airborne measurements for both mean values – also on a year-to-year basis – as well as seasonal variations. Several locations show large biases that are attributed to local effects like orography and proximity of large emission sources. Differences were detected for individual years: 2003, 2004 and 2006 have larger biases than 2005, 2007 and 2008, which appear to be related to SCIAMACHY instrumental issues but require more research. Results from this study are consistent with, and complementary to, findings from a previous validation study using ground-based measurements (de Laat et al., 2010b). According to this study, the SCIAMACHY data, if individual measurements are of sufficient quality – good signal-to-noise, can be used to determine the spatial distribution and seasonal cycles of CO total columns over clean areas. Biases found over areas with strong emissions (Africa, China) could be explained by low sensitivity of the instrument in the boundary layer and users are recommended to avoid using the SCIAMACHY data while trying to quantify CO burden and/or retrieve CO emissions in such areas.

Description: Corrigendum to "Trend analysis of aerosol optical thickness and Ångström exponent derived from the global AERONET spectral observations" published in Atmos. Meas. Tech., 5, 1271–1299, 2012 Atmospheric Measurement Techniques, 5, 2113-2113, 2012 Author(s): J. Yoon, W. von Hoyningen-Huene, A. A. Kokhanovsky, M. Vountas, and J. P. Burrows No abstract available.

Description: Using sonic anemometer temperature to measure sensible heat flux in strong winds Atmospheric Measurement Techniques, 5, 2095-2111, 2012 Author(s): S. P. Burns, T. W. Horst, L. Jacobsen, P. D. Blanken, and R. K. Monson Sonic anemometers simultaneously measure the turbulent fluctuations of vertical wind ( w ') and sonic temperature ( T s '), and are commonly used to measure sensible heat flux ( H ). Our study examines 30-min heat fluxes measured with a Campbell Scientific CSAT3 sonic anemometer above a subalpine forest. We compared H calculated with T s to H calculated with a co-located thermocouple and found that, for horizontal wind speed ( U ) less than 8 m s −1 , the agreement was around ±30 W m −2 . However, for U ≈ 8 m s −1 , the CSAT H had a generally positive deviation from H calculated with the thermocouple, reaching a maximum difference of ≈250 W m −2 at U ≈ 18 m s −1 . With version 4 of the CSAT firmware, we found significant underestimation of the speed of sound and thus T s in high winds (due to a delayed detection of the sonic pulse), which resulted in the large CSAT heat flux errors. Although this T s error is qualitatively similar to the well-known fundamental correction for the crosswind component, it is quantitatively different and directly related to the firmware estimation of the pulse arrival time. For a CSAT running version 3 of the firmware, there does not appear to be a significant underestimation of T s ; however, a T s error similar to that of version 4 may occur if the CSAT is sufficiently out of calibration. An empirical correction to the CSAT heat flux that is consistent with our conceptual understanding of the T s error is presented. Within a broader context, the surface energy balance is used to evaluate the heat flux measurements, and the usefulness of side-by-side instrument comparisons is discussed.

Description: Aerosol-climate interactions in the Norwegian Earth System Model – NorESM Geoscientific Model Development Discussions, 5, 2599-2685, 2012 Author(s): A. Kirkevåg, T. Iversen, Ø. Seland, C. Hoose, J. E. Kristjánsson, H. Struthers, A. M. L. Ekman, S. Ghan, J. Griesfeller, E. D. Nilsson, and M. Schulz The objective of this study is to document and evaluate recent changes and updates to the module for aerosols and aerosol-cloud-radiation interactions in the atmospheric module CAM4-Oslo of the Norwegian Earth System Model (NorESM). Particular attention is paid to the role of natural organics, sea salt, and mineral dust in determining the gross aerosol properties as well as the anthropogenic contribution to these properties and the associated direct and indirect radiative forcing. The aerosol module is extended from earlier versions that have been published, and includes life-cycling of sea-salt, mineral dust, particulate sulphate, black carbon, and primary and secondary organics. The impacts of most of the numerous changes since previous versions are thoroughly explored by sensitivity experiments. The most important changes are: modified prognostic sea salt emissions; updated treatment of precipitation scavenging and gravitational settling; inclusion of biogenic primary organics and methane sulphonic acid (MSA) from oceans; almost doubled production of land-based biogenic secondary organic aerosols (SOA); and increased ratio of organic matter to organic carbon (OM / OC) for biomass burning aerosols from 1.4 to 2.6. Compared with in-situ measurements and remotely sensed data, the new treatments of sea salt and dust aerosols give smaller biases in near surface mass concentrations and aerosol optical depth than in the earlier model version. The model biases for mass concentrations are approximately unchanged for sulphate and BC. The enhanced levels of modeled OM yield improved overall statistics, even though OM is still underestimated in Europe and over-estimated in North America. The global direct radiative forcing (DRF) at the top of the atmosphere has changed from a small positive value to −0.08 W m −2 in CAM4-Oslo. The sensitivity tests suggest that this change can be attributed to the new treatment of biomass burning aerosols and gravitational settling. Although it has not been a goal in this study, the new DRF estimate is closer both to the median model estimate from the AeroCom inter-comparison and the best estimate in IPCC AR4. Estimated DRF at the ground surface has increased by ca. 60%, to −1.89 W m −2 . We show that this can be explained by new emission data and omitted mixing of constituents between updrafts and downdrafts in convective clouds. The increased abundance of natural OM and the introduction of a cloud droplet spectral dispersion formulation are the most important contributions to a considerably decreased estimate of the indirect radiative forcing (IndRF). The IndRF is also found to be sensitive to assumptions about the coating of insoluble aerosols by sulphate and OM. The IndRF of −1.2 W m −2 , which is closer to the IPCC AR4 estimates than the previous estimate of −1.9 W m −2 , has thus been obtained without imposing unrealistic artificial lower bounds on cloud droplet number concentrations.

Description: Inclusion of Ash and SO 2 emissions from volcanic eruptions in WRF-CHEM: development and some applications Geoscientific Model Development Discussions, 5, 2571-2597, 2012 Author(s): M. Stuefer, S. R. Freitas, G. Grell, P. Webley, S. Peckham, and S. A. McKeen We describe a new functionality within the Weather Research and Forecasting model with coupled Chemistry (WRF-Chem) that allows simulating emission, transport, dispersion, transformation and sedimentation of pollutants released during volcanic activities. Emissions from both an explosive eruption case and relatively calm degassing situation are considered using the most recent volcanic emission databases. A preprocessor tool provides emission fields and additional information needed to establish the initial three-dimensional cloud umbrella/vertical distribution within the transport model grid, as well as the timing and duration of an eruption. From this source condition, the transport, dispersion and sedimentation of the ash-cloud can be realistically simulated by WRF-Chem using its own dynamics, physical parameterization as well as data assimilation. Examples of model validation include a comparison of tephra fall deposits from the 1989 eruption of Mount Redoubt (Alaska), and the dispersion of ash from the 2010 Eyjafjallajökull eruption in Iceland. Both model applications show good coincidence between WRF-Chem and observations.

Description: Remote sensing of near-infrared chlorophyll fluorescence from space in scattering atmospheres: implications for its retrieval and interferences with atmospheric CO 2 retrievals Atmospheric Measurement Techniques, 5, 2081-2094, 2012 Author(s): C. Frankenberg, C. O'Dell, L. Guanter, and J. McDuffie With the advent of dedicated greenhouse gas space-borne spectrometers sporting high resolution spectra in the O 2 A-band spectral region (755–774 nm), the retrieval of chlorophyll fluorescence has become feasible on a global scale. If unaccounted for, however, fluorescence can indirectly perturb the greenhouse gas retrievals as it perturbs the oxygen absorption features. As atmospheric CO 2 measurements are used to invert net fluxes at the land–atmosphere interface, a bias caused by fluorescence can be crucial as it will spatially correlate with the fluxes to be inverted. Avoiding a bias and retrieving fluorescence accurately will provide additional constraints on both the net and gross fluxes in the global carbon cycle. We show that chlorophyll fluorescence, if neglected, systematically interferes with full-physics multi-band X CO 2 retrievals using the O 2 A-band. Systematic biases in X CO 2 can amount to +1 ppm if fluorescence constitutes 1% to the continuum level radiance. We show that this bias can be largely eliminated by simultaneously fitting fluorescence in a full-physics based retrieval. If fluorescence is the primary target, a dedicated but very simple retrieval based purely on Fraunhofer lines is shown to be more accurate and very robust even in the presence of large scattering optical depths. We find that about 80% of the surface fluorescence is retained at the top-of-atmosphere, even for cloud optical thicknesses around 2–5. We further show that small instrument modifications to future O 2 A-band spectrometer spectral ranges can result in largely reduced random errors in chlorophyll fluorescence, paving the way towards a more dedicated instrument exploiting solar absorption features only.

Description: Net primary productivity, upwelling and coastal currents in the Gulf of Ulloa, Baja California, México Ocean Science, 8, 703-711, 2012 Author(s): E. González-Rodríguez, A. Trasviña-Castro, G. Gaxiola-Castro, L. Zamudio, and R. Cervantes-Duarte The Gulf of Ulloa, a highly productive area off the western coast of the Baja California Peninsula, is examined for five successive years (2003–2007) by using satellite data and seasonal net primary productivity (NPP) estimates obtained from a vertical generalised production model. The results identify that northwestern winds blow parallel to the coast throughout the year. However, highest NPP occurs from March to June. During this period, an equatorward coastal current transports water from neighbouring upwelling areas to the northern Gulf of Ulloa and in combination with local upwelling, which injects nutrients into the euphotic zone, produce the observed increase in NPP. The opposite situation occurs in late summer when a warm poleward current of tropical characteristics arrives and inhibits the productivity in the whole region and generates the yearly lowest NPP levels. Our findings reveal the importance of lateral advection in the modulation of the primary productivity in this subtropical upwelling region.

Description: Assessment of a physical-biogeochemical coupled model system for operational service in the Baltic Sea Ocean Science, 8, 683-701, 2012 Author(s): Z. Wan, J. She, M. Maar, L. Jonasson, and J. Baasch-Larsen Thanks to the abundant observation data, we are able to deploy the traditional point-to-point comparison and statistical measures in combination with a comprehensive model validation scheme to assess the skills of the biogeochemical model ERGOM in providing an operational service for the Baltic Sea. The model assessment concludes that the operational products can resolve the main observed seasonal features for phytoplankton biomass, dissolved inorganic nitrogen, dissolved inorganic phosphorus and dissolved oxygen in euphotic layers as well as their vertical profiles. This assessment reflects that the model errors of the operational system at the current stage are mainly caused by insufficient light penetration, excessive organic particle export downward, insufficient regional adaptation and some from improper initialization. This study highlights the importance of applying multiple schemes in order to assess model skills rigidly and identify main causes for major model errors.

Description: A web-based software tool to estimate unregulated daily streamflow at ungauged rivers Geoscientific Model Development Discussions, 5, 2503-2526, 2012 Author(s): S. A. Archfield, P. A. Steeves, J. D. Guthrie, and K. G. Ries III Streamflow information is critical for solving any number of hydrologic problems. Often times, streamflow information is needed at locations which are ungauged and, therefore, have no observations on which to base water management decisions. Furthermore, there has been increasing need for daily streamflow time series to manage rivers for both human and ecological functions. To facilitate negotiation between human and ecological demands for water, this paper presents the first publically-available, map-based, regional software tool to interactively estimate daily streamflow time series at any user-selected ungauged river location. The map interface allows users to locate and click on a river location, which then returns estimates of daily streamflow for the location selected. For the demonstration region in the northeast United States, daily streamflow was shown to be reliably estimated by the software tool, with efficiency values computed from observed and estimated streamflows ranging from 0.69 to 0.92. The software tool provides a general framework that can be applied to other regions for which daily streamflow estimates are needed.

Description: An operational model for the West Iberian coast: products and services Ocean Science, 8, 713-732, 2012 Author(s): M. Mateus, G. Riflet, P. Chambel, L. Fernandes, R. Fernandes, M. Juliano, F. Campuzano, H. de Pablo, and R. Neves This paper presents the structure and application of a regional scale operational modelling tool for the West Iberian coast, and discusses its potential for products and services for both scientific and coastal management activities. The forecasting suite includes nested hydrodynamic models forced with up-to-date meteorological forecast data and large-scale model results. The present status of the system and its recent upgrades are reviewed, offering a general description of the main components of the system: the forcing data, the circulation model, the model outputs and the validation methodology of model results. Seasonal differences in temperature, salinity and current velocity fields are illustrated and show satisfactory reproduction of the top and deep layer thermodynamics. The system provides boundary forcing for a number of local-scale model applications via downscaling of the solution and enables potential products and services from which civil society will benefit.

Description: On the shelf resonances of the Gulf of Carpentaria and the Arafura Sea Ocean Science, 8, 733-750, 2012 Author(s): D. J. Webb A numerical model is used to investigate the resonances of the Gulf of Carpentaria and the Arafura Sea, and the additional insights that come from extending the analysis into the complex angular velocity plane. When the model is forced at the shelf edge with physically realistic real values of the angular velocity, the response functions at points within the region show maxima and other behaviour which imply that resonances are involved but provide little additional information. The study is then extended to complex angular velocities, and the results then show a clear pattern of gravity wave and Rossby wave like resonances. The properties of the resonances are investigated and used to reinterpret the response at real values of angular velocity. It is found that in some regions the response is dominated by modes trapped between the shelf edge and the coast or between opposing coastlines. In other regions the resonances show cooperative behaviour, possibly indicating the importance of other physical processes.

Description: The analysis of large-scale turbulence characteristics in the Indonesian seas derived from a regional model based on the Princeton Ocean Model Ocean Science, 8, 615-631, 2012 Author(s): K. O'Driscoll and V. Kamenkovich Turbulence characteristics in the Indonesian seas on the horizontal scale of order of 100 km were calculated with a regional model of the Indonesian seas circulation in the area based on the Princeton Ocean Model (POM). As is well known, the POM incorporates the Mellor–Yamada turbulence closure scheme. The calculated characteristics are: twice the turbulence kinetic energy per unit mass, q 2 ; the turbulence master scale, ℓ; mixing coefficients of momentum, K M ; and temperature and salinity, K H ; etc. The analyzed turbulence has been generated essentially by the shear of large-scale ocean currents and by the large-scale wind turbulence. We focused on the analysis of turbulence around important topographic features, such as the Lifamatola Sill, the North Sangihe Ridge, the Dewakang Sill, and the North and South Halmahera Sea Sills. In general, the structure of turbulence characteristics in these regions turned out to be similar. For this reason, we have carried out a detailed analysis of the Lifamatola Sill region because dynamically this region is very important and some estimates of mixing coefficients in this area are available. Briefly, the main results are as follows. The distribution of q 2 is quite adequately reproduced by the model. To the north of the Lifamatola Sill (in the Maluku Sea) and to the south of the Sill (in the Seram Sea), large values of q 2 occur in the deep layer extending several hundred meters above the bottom. The observed increase of q 2 near the very bottom is probably due to the increase of velocity shear and the corresponding shear production of q 2 very close to the bottom. The turbulence master scale, ℓ, was found to be constant in the main depth of the ocean, while ℓ rapidly decreases close to the bottom, as one would expect. However, in deep profiles away from the sill, the effect of topography results in the ℓ structure being unreasonably complicated as one moves towards the bottom. Values of 15 to 20 × 10 −4 m 2 s −1 were obtained for K M and K H in deep water in the vicinity of the Lifamatola Sill. These estimates agree well with basin-scale averaged values of 13.3 × 10 −4 m 2 s −1 found diagnostically for K H in the deep Banda and Seram Seas (Gordon et al., 2003) and a value of 9.0 × 10 −4 m 2 s −1 found diagnostically for K H for the deep Banda Sea system (van Aken et al., 1988). The somewhat higher simulated values can be explained by the presence of steep topography around the sill.

Description: Airborne intercomparison of HO x measurements using laser-induced fluorescence and chemical ionization mass spectrometry during ARCTAS Atmospheric Measurement Techniques, 5, 2025-2037, 2012 Author(s): X. Ren, J. Mao, W. H. Brune, C. A. Cantrell, R. L. Mauldin III, R. S. Hornbrook, E. Kosciuch, J. R. Olson, J. H. Crawford, G. Chen, and H. B. Singh The hydroxyl (OH) and hydroperoxyl (HO 2 ) radicals, collectively called HO x , play central roles in tropospheric chemistry. Accurate measurements of OH and HO 2 are critical to examine our understanding of atmospheric chemistry. Intercomparisons of different techniques for detecting OH and HO 2 are vital to evaluate their measurement capabilities. Three instruments that measured OH and/or HO 2 radicals were deployed on the NASA DC-8 aircraft throughout Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) in the spring and summer of 2008. One instrument was the Penn State Airborne Tropospheric Hydrogen Oxides Sensor (ATHOS) for OH and HO 2 measurements based on Laser-Induced Fluorescence (LIF) spectroscopy. A second instrument was the NCAR Selected-Ion Chemical Ionization Mass Spectrometer (SI-CIMS) for OH measurement. A third instrument was the NCAR Peroxy Radical Chemical Ionization Mass Spectrometer (PeRCIMS) for HO 2 measurement. Formal intercomparison of LIF and CIMS was conducted for the first time on a same aircraft platform. The three instruments were calibrated by quantitative photolysis of water vapor by ultraviolet (UV) light at 184.9 nm with three different calibration systems. The absolute accuracies were ±32% (2σ) for the LIF instrument, ±65% (2σ) for the SI-CIMS instrument, and ±50% (2σ) for the PeRCIMS instrument. In general, good agreement was obtained between the CIMS and LIF measurements of both OH and HO 2 measurements. Linear regression of the entire data set yields [OH] CIMS = 0.89 × [OH] LIF + 2.8 × 10 4 cm −3 with a correlation coefficient r 2 = 0.72 for OH, and [HO 2 ] CIMS = 0.86 × [HO 2 ] LIF + 3.9 parts per trillion by volume (pptv, equivalent to pmol mol −1 ) with a correlation coefficient r 2 = 0.72 for HO 2 . In general, the difference between CIMS and LIF instruments for OH and HO 2 measurements can be explained by their combined measurement uncertainties. Comparison with box model results shows some similarities for both the CIMS and LIF measurements. First, the observed-to-modeled HO 2 ratio increases greatly for higher NO mixing ratios, indicating that the model may not properly account for HO x sources that correlate with NO. Second, the observed-to-modeled OH ratio increases with increasing isoprene mixing ratios, suggesting either incomplete understanding of isoprene chemistry in the model or interferences in the measurements in environments where biogenic emissions dominate ambient volatile organic compounds.

Description: TOPAZ4: an ocean-sea ice data assimilation system for the North Atlantic and Arctic Ocean Science, 8, 633-656, 2012 Author(s): P. Sakov, F. Counillon, L. Bertino, K. A. Lisæter, P. R. Oke, and A. Korablev We present a detailed description of TOPAZ4, the latest version of TOPAZ – a coupled ocean-sea ice data assimilation system for the North Atlantic Ocean and Arctic. It is the only operational, large-scale ocean data assimilation system that uses the ensemble Kalman filter. This means that TOPAZ features a time-evolving, state-dependent estimate of the state error covariance. Based on results from the pilot MyOcean reanalysis for 2003–2008, we demonstrate that TOPAZ4 produces a realistic estimate of the ocean circulation in the North Atlantic and the sea-ice variability in the Arctic. We find that the ensemble spread for temperature and sea-level remains fairly constant throughout the reanalysis demonstrating that the data assimilation system is robust to ensemble collapse. Moreover, the ensemble spread for ice concentration is well correlated with the actual errors. This indicates that the ensemble statistics provide reliable state-dependent error estimates – a feature that is unique to ensemble-based data assimilation systems. We demonstrate that the quality of the reanalysis changes when different sea surface temperature products are assimilated, or when in-situ profiles below the ice in the Arctic Ocean are assimilated. We find that data assimilation improves the match to independent observations compared to a free model. Improvements are particularly noticeable for ice thickness, salinity in the Arctic, and temperature in the Fram Strait, but not for transport estimates or underwater temperature. At the same time, the pilot reanalysis has revealed several flaws in the system that have degraded its performance. Finally, we show that a simple bias estimation scheme can effectively detect the seasonal or constant bias in temperature and sea-level.

Description: Using multi-model averaging to improve the reliability of catchment scale nitrogen predictions Geoscientific Model Development Discussions, 5, 2289-2310, 2012 Author(s): J.-F. Exbrayat, N. R. Viney, H.-G. Frede, and L. Breuer Hydro-biogeochemical models are used to foresee the impact of mitigation measures on water quality. Usually, scenario-based studies rely on single model applications. This is done in spite of the widely acknowledged advantage of ensemble approaches to cope with structural model uncertainty issues. As an attempt to demonstrate the reliability of such multi-model efforts in the hydro-biogeochemical context, this methodological contribution proposes an adaptation of the Reliability Ensemble Averaging (REA) philosophy to nitrogen losses predictions. A total of 4 models are used to predict the total nitrogen (TN) losses from the well-monitored Ellen Brook catchment in Western Australia. Simulations include re-predictions of current conditions and a set of straightforward management changes targeting fertilization scenarios. Results show that, in spite of good calibration metrics, one of the models provides a very different response to management changes. This behaviour leads the simple average of the ensemble members to also predict reductions in TN export that are not in agreement with the other models. However, considering the convergence of model predictions in the more sophisticated REA approach assigns more weight to previously less well calibrated models that are more in agreement with each other. This method also avoids having to disqualify any of the ensemble members, which is always sensible.

Description: Consistency between Fourier transform and small-volume few-wave decomposition for spectral and spatial variability of gravity waves above a typhoon Atmospheric Measurement Techniques, 5, 1637-1651, 2012 Author(s): C. I. Lehmann, Y.-H. Kim, P. Preusse, H.-Y. Chun, M. Ern, and S.-Y. Kim Convective gravity wave (GW) sources are spatially localized and emit at the same time waves with a wide spectrum of phase speeds. Any wave analysis therefore compromises between spectral and spatial resolution. Future satellite borne limb imagers will for a first time provide real 3-D volumes of observations. These volumes will be however limited which will impose further constraints on the analysis technique. In this study a three dimensional few-wave approach fitting sinusoidal waves to limited 3-D volumes is introduced. The method is applied to simulated GWs above typhoon Ewiniar and GW momentum flux is estimated from temperature fluctuations. Phase speed spectra as well as average profiles of positive, negative and net momentum fluxes are compared to momentum flux estimated by Fourier transform as well as spatial averaging of wind fluctuations. The results agree within 10–20%. The few-wave method can also reveal the spatial orientation of the GWs with respect to the source. The relevance of the results for different types of measurements as well as its applicability to model data is discussed.

Description: Satellite retrieval of the liquid water fraction in tropical clouds between −20 and −38 °C Atmospheric Measurement Techniques, 5, 1683-1698, 2012 Author(s): D. L. Mitchell and R. P. d'Entremont This study describes a satellite remote sensing method for directly retrieving the liquid water fraction in mixed phase clouds, and appears unique in this respect. The method uses MODIS split-window channels for retrieving the liquid fraction from cold clouds where the liquid water fraction is less than 50% of the total condensate. This makes use of the observation that clouds only containing ice exhibit effective 12-to-11 μm absorption optical thickness ratios (β eff ) that are quasi-constant with retrieved cloud temperature T . This observation was made possible by using two CO 2 channels to retrieve T and then using the 12 and 11 μm channels to retrieve emissivities and β eff . Thus for T 〈 −40 °C, β eff is constant, but for T 〉 −40 °C, β eff slowly increases due to the presence of liquid water, revealing mean liquid fractions of ~ 10% around −22 °C from tropical clouds identified as cirrus by the cloud mask. However, the uncertainties for these retrievals are large, and extensive in situ measurements are needed to refine and validate these retrievals. Such liquid levels are shown to reduce the cloud effective diameter D e such that cloud optical thickness will increase by more than 50% for a given water path, relative to D e corresponding to pure ice clouds. Such retrieval information is needed for validation of the cloud microphysics in climate models. Since low levels of liquid water can dominate cloud optical properties, tropical clouds between −25 and −20 °C may be susceptible to the first aerosol indirect effect.

Description: Atlantic transport variability at 25° N in six hydrographic sections Ocean Science, 8, 497-523, 2012 Author(s): C. P. Atkinson, H. L. Bryden, S. A. Cunningham, and B. A. King In January and February 2010, a sixth transatlantic hydrographic section was completed across 25° N, extending the hydrographic record at this latitude to over half a century. In combination with continuous transport measurements made since 2004 at 26.5° N by the Rapid-WATCH project, we reassess transport variability in the 25° N hydrographic record. Past studies of transport variability at this latitude have assumed transport estimates from each hydrographic section to represent annual average conditions. In this study the uncertainty in this assumption is assessed through use of Rapid-WATCH observations to quantify sub-seasonal and seasonal transport variability. Whilst in the upper-ocean no significant interannual or decadal transport variability are identified in the hydrographic record, in the deep ocean transport variability in both depth and potential temperature classes suggests some interannual or decadal variability may have occurred. This is particularly striking in the lower North Atlantic Deep Water where southward transports prior to 1998 were greater than recent transports by several Sverdrups. Whilst a cooling and freshening of Denmark Straits Overflow Water has occurred which is coincident with these transport changes, these water mass changes appear to be density compensated. Transport changes are the result of changing velocity shear in the vicinity of the Deep Western Boundary Current.

Description: Propagation and dissipation of internal tides in the Oslofjord Ocean Science, 8, 525-543, 2012 Author(s): A. Staalstrøm, E. Aas, and B. Liljebladh Observations of velocity, pressure, temperature and salinity in the inner Oslofjord have been analysed to provide new information about the relationships between internal tides generated by tidal currents across the Drøbak Sill and dissipation and diffusivity in the fjord. The most energetic vertical displacement of density surfaces inside the sill is associated with the first internal mode that has maximum amplitude around sill depth. The amplitude of the vertical displacement around sill depth correlates with the amplitude of the surface elevation, and, at a distance of 1 km inside the sill, the ratio between the amplitudes is 38, decreasing to 11 at a distance of 10 km. The greatest vertical displacements inside the sill, however, are found at 40 m depth. These latter internal waves are not associated with a first-mode internal tide, but are rather associated with higher internal modes controlled by stratification. The energy flux of the internal wave propagating from the Drøbak Sill into the inner fjord on the east side of the Håøya Island is estimated to vary in the range 155–430 kW. This is the same order of magnitude as the estimated barotropic energy loss over the Drøbak Sill (250 kW), but only 4–10% of the total barotropic flux. Approximately 40–70% of the internal energy flux is lost within a distance of 10 km from the sill. The mean diffusivity below 90 m depth in this area (~20 cm 2 s −1 ) is more than four times higher than in the rest of the fjord (~5 cm 2 s −1 or less).

Description: Carbon monoxide measurements onboard the CARIBIC passenger aircraft using UV resonance fluorescence Atmospheric Measurement Techniques, 5, 1753-1760, 2012 Author(s): D. Scharffe, F. Slemr, C. A. M. Brenninkmeijer, and A. Zahn Goal of the project CARIBIC (Civil Aircraft for the Regular Investigation of the atmosphere Based on an Instrumented Container) is to carry out regular and detailed observations of atmospheric composition (particles and gases) at cruising altitudes of passenger aircraft, i.e. at 9–12 km. Continuous, fast measurement of CO is indispensable for the chemical characterization of encountered air masses, for the detection of plumes of polluted air and for studying troposphere-stratosphere transport. CO is measured by a commercial resonance fluorescence UV instrument modified for the use onboard passenger aircraft. Modifications were necessary to optimize the instrument reliability allowing unattended operation for several days. The instrument has a precision of 1–2 ppbv at an integration time of 1 s. The response time to reach 63.2% signal strength is 2 s. We describe the modifications of the instrument, the experiences made during its operation since December 2004, the quality control of CO measurements onboard CARIBIC, and suggest a regular service routine that guarantees long-term high-quality data.

Description: Development of high resolution land surface parameters for the Community Land Model Geoscientific Model Development Discussions, 5, 1435-1481, 2012 Author(s): Y. Ke, L. R. Leung, M. Huang, A. M. Coleman, H. Li, and M. S. Wigmosta There is a growing need for high-resolution land surface parameters as land surface models are being applied at increasingly higher spatial resolution offline as well as in regional and global models. The default land surface parameters for the most recent version of the Community Land Model (i.e. CLM 4.0) are at 0.5° or coarser resolutions, released with the model from the National Center for Atmospheric Research (NCAR). Plant Functional Types (PFTs), vegetation properties such as Leaf Area Index (LAI), Stem Area Index (SAI), and non-vegetated land covers were developed using remotely-sensed datasets retrieved in late 1990's and the beginning of this century. In this study, we developed new land surface parameters for CLM 4.0, specifically PFTs, LAI, SAI and non-vegetated land cover composition, at 0.05° resolution globally based on the most recent MODIS land cover and improved MODIS LAI products. Compared to the current CLM 4.0 parameters, the new parameters produced a decreased coverage by bare soil and trees, but an increased coverage by shrub, grass, and cropland. The new parameters result in a decrease in global seasonal LAI, with the biggest decrease in boreal forests; however, the new parameters also show a large increase in LAI in tropical forest. Differences between the new and the current parameters are mainly caused by changes in the sources of remotely sensed data and the representation of land cover in the source data. The new high-resolution land surface parameters have been used in a coupled land-atmosphere model (WRF-CLM) applied to the western US to demonstrate their use in high-resolution modeling. Future work will include global offline CLMsimulations to examine the impacts of source data resolution and subsequent land parameter changes on simulated land surface processes.

Description: The Model of Emissions of Gases and Aerosols from Nature version 2.1 (MEGAN2.1): an extended and updated framework for modeling biogenic emissions Geoscientific Model Development Discussions, 5, 1503-1560, 2012 Author(s): A. B. Guenther, X. Jiang, C. L. Heald, T. Sakulyanontvittaya, T. Duhl, L. K. Emmons, and X. Wang The Model of Emissions of Gases and Aerosols from Nature version 2.1 (MEGAN2.1) is a modeling framework for estimating fluxes of 147 biogenic compounds between terrestrial ecosystems and the atmosphere using simple mechanistic algorithms to account for the major known processes controlling biogenic emissions. It is available as an offline code and has also been coupled into land surface models and atmospheric chemistry models. MEGAN2.1 is an update from the previous versions including MEGAN2.0 for isoprene emissions and MEGAN2.04, which estimates emissions of 138 compounds. Isoprene comprises about half of the estimated total global biogenic volatile organic compound (BVOC) emission of 1 Pg (1000 Tg or 10 15 g). Another 10 compounds including methanol, ethanol, acetaldehyde, acetone, α-pinene, β-pinene, t −β-ocimene, limonene, ethene, and propene together contribute another 30% of the estimated emission. An additional 20 compounds (mostly terpenoids) are associated with another 17% of the total emission with the remaining 3% distributed among 125 compounds. Emissions of 41 monoterpenes and 32 sesquiterpenes together comprise about 15% and 3%, respectively, of the total global BVOC emission. Tropical trees cover about 18% of the global land surface and are estimated to be responsible for 60% of terpenoid emissions and 48% of other VOC emissions. Other trees cover about the same area but are estimated to contribute only about 10% of total emissions. The magnitude of the emissions estimated with MEGAN2.1 are within the range of estimates reported using other approaches and much of the differences between reported values can be attributed to landcover and meteorological driving variables. The offline version of MEGAN2.1 source code and driving variables is available from http://acd.ucar.edu/~guenther/MEGAN/MEGAN.htm and the version integrated into the Community Land Model version 4 (CLM4) can be downloaded from http://www.cesm.ucar.edu/ .

Description: Evaluation of turbulent dissipation rate retrievals from Doppler Cloud Radar Atmospheric Measurement Techniques, 5, 1375-1385, 2012 Author(s): M. D. Shupe, I. M. Brooks, and G. Canut Turbulent dissipation rate retrievals from cloud radar Doppler velocity measurements are evaluated using independent, in situ observations in Arctic stratocumulus clouds. In situ validation data sets of dissipation rate are derived using sonic anemometer measurements from a tethered balloon and high frequency pressure variation observations from a research aircraft, both flown in proximity to stationary, ground-based radars. Modest biases are found among the data sets in particularly low- or high-turbulence regimes, but in general the radar-retrieved values correspond well with the in situ measurements. Root mean square differences are typically a factor of 4–6 relative to any given magnitude of dissipation rate. These differences are no larger than those found when comparing dissipation rates computed from tethered-balloon and meteorological tower-mounted sonic anemometer measurements made at spatial distances of a few hundred meters. Temporal lag analyses suggest that approximately half of the observed differences are due to spatial sampling considerations, such that the anticipated radar-based retrieval uncertainty is on the order of a factor of 2–3. Moreover, radar retrievals are clearly able to capture the vertical dissipation rate structure observed by the in situ sensors, while offering substantially more information on the time variability of turbulence profiles. Together these evaluations indicate that radar-based retrievals can, at a minimum, be used to determine the vertical structure of turbulence in Arctic stratocumulus clouds.

Description: CH 4 , CO, and H 2 O spectroscopy for the Sentinel-5 Precursor mission: an assessment with the Total Carbon Column Observing Network measurements Atmospheric Measurement Techniques, 5, 1387-1398, 2012 Author(s): A. Galli, A. Butz, R. A. Scheepmaker, O. Hasekamp, J. Landgraf, P. Tol, D. Wunch, N. M. Deutscher, G. C. Toon, P. O. Wennberg, D. W. T. Griffith, and I. Aben The TROPOspheric Monitoring Instrument (TROPOMI) will be part of ESA's Sentinel-5 Precursor (S5P) satellite platform scheduled for launch in 2015. TROPOMI will monitor methane and carbon monoxide concentrations in the Earth's atmosphere by measuring spectra of back-scattered sunlight in the short-wave infrared (SWIR). S5P will be the first satellite mission to rely uniquely on the spectral window at 4190–4340 cm −1 (2.3 μm) to retrieve CH 4 and CO. In this study, we investigated if the absorption features of the three relevant molecules CH 4 , CO, and H 2 O are adequately known. To this end, we retrieved total columns of CH 4 , CO, and H 2 O from absorption spectra measured by two ground-based Fourier transform spectrometers that are part of the Total Carbon Column Observing Network (TCCON). The retrieval results from the 4190–4340 cm −1 range at the TROPOMI resolution (0.45 cm −1 ) were then compared to the CH 4 results obtained from the 6000 cm −1 region, and the CO results obtained from the 4190–4340 cm −1 region at the higher TCCON resolution (0.02 cm −1 ). For TROPOMI-like settings, we were able to reproduce the CH 4 columns to an accuracy of 0.3% apart from a constant bias of 1%. The CO retrieval accuracy was, through interference, systematically influenced by the shortcomings of the CH 4 and H 2 O spectroscopy. In contrast to CH 4 , the CO column error also varied significantly with atmospheric H 2 O content. Unaddressed, this would introduce seasonal and latitudinal biases to the CO columns retrieved from TROPOMI measurements. We therefore recommend further effort from the spectroscopic community to be directed at the H 2 O and CH 4 spectroscopy in the 4190–4340 cm −1 region.

Description: New Aura Microwave Limb Sounder observations of BrO and implications for Br y Atmospheric Measurement Techniques, 5, 1741-1751, 2012 Author(s): L. Millán, N. Livesey, W. Read, L. Froidevaux, D. Kinnison, R. Harwood, I. A. MacKenzie, and M. P. Chipperfield This paper introduces a new inversion algorithm for retrievals of stratospheric BrO from the Aura Microwave Limb Sounder. This version is based on the algorithm described by Livesey et al. (2006a) but uses a more realistic atmospheric state to constrain the retrieval. A description of the methodology and an error analysis are presented. Single daily profile precision uncertainty, when taking the ascending-descending (day-night) difference, was found to be up to 40 pptv while systematic error biases were estimated to be less than about 3 pptv. Monthly mean comparisons show broad agreement with other measurements as well as with state-of-the-art numerical models. We infer a 2005 yearly total inorganic Br y using the measured MLS BrO to be 20.3 ± 4.5 pptv, which implies a contribution from very short lived substances to the stratospheric bromine budget of 5 ± 4.5 pptv.

Description: A simulation study of the ensemble-based data assimilation of satellite-borne lidar aerosol observations Geoscientific Model Development Discussions, 5, 1877-1947, 2012 Author(s): T. T. Sekiyama, T. Y. Tanaka, and T. Miyoshi A four-dimensional ensemble-based data assimilation system was assessed by observing system simulation experiments (OSSEs), in which the CALIPSO satellite was emulated via simulated satellite-borne lidar aerosol observations. Its performance over athree-month period was validated according to the Method for Object-based Diagnostic Evaluation (MODE), using aerosol optical thickness (AOT) distributions in East Asia as the objects of analysis. Consequently, this data assimilation system demonstrated the ability to produce better analyses of sulfate and dust aerosols in comparison to a free-running simulation model. For example, the mean centroid distance (from the truth) over a three-month collection period of aerosol plumes was improved from 2.15 grids (≈ 600 km) to 1.45 grids (≈ 400 km) for sulfate aerosols and from 2.59 grids (≈ 750 km) to 1.14 grids (≈ 330 km) for dust aerosols; the mean area ratio (to the truth) over a three-month collection period of aerosol plumes was improved from 0.49 to 0.76 for sulfate aerosols and from 0.51 to 0.72 for dust aerosols. The satellite-borne lidar data assimilation successfully improved the aerosol plume analysis and the dust emission estimation in the OSSEs. These results present great possibilities for the beneficial use of lidar data, whose distribution is vertically/temporally dense but horizontally sparse, when coupled with a four-dimensional data assimilation system. In addition, sensitivity tests were conducted, and their results indicated that the degree of freedom to control the aerosol variables was probably limited in the data assimilation because the meteorological field in the system was constrained to weather reanalysis using Newtonian relaxation. Further improvements to the aerosol analysis can be performed through the simultaneous assimilation of aerosol observations with meteorological observations. The OSSE results strongly suggest that the use of real CALIPSO data will have a beneficial effect on obtaining more accurate sulfate and dust aerosol analyses. Furthermore, the use of the same OSSE technique will allow us to perform a prior assessment of the next-generation lidar satellite EarthCARE, which will be launched in 2015.

Description: Multi-wavelength Raman lidar, sun photometric and aircraft measurements in combination with inversion models for the estimation of the aerosol optical and physico-chemical properties over Athens, Greece Atmospheric Measurement Techniques, 5, 1793-1808, 2012 Author(s): R. E. Mamouri, A. Papayannis, V. Amiridis, D. Müller, P. Kokkalis, S. Rapsomanikis, E. T. Karageorgos, G. Tsaknakis, A. Nenes, S. Kazadzis, and E. Remoundaki A novel procedure has been developed to retrieve, simultaneously, the optical, microphysical and chemical properties of tropospheric aerosols with a multi-wavelength Raman lidar system in the troposphere over an urban site (Athens, Greece: 37.9° N, 23.6° E, 200 m a.s.l.) using data obtained during the European Space Agency (ESA) THERMOPOLIS project, which took place between 15–31 July 2009 over the Greater Athens Area (GAA). We selected to apply our procedure for a case study of intense aerosol layers that occurred on 20–21 July 2009. The National Technical University of Athens (NTUA) EOLE 6-wavelength Raman lidar system has been used to provide the vertical profiles of the optical properties of aerosols (extinction and backscatter coefficients, lidar ratio) and the water vapor mixing ratio. An inversion algorithm was used to derive the mean aerosol microphysical properties (mean effective radius ( r eff ), single-scattering albedo ω ) and mean complex refractive index ( m )) at selected heights in the 2–3 km height region. We found that r eff was 0.14–0.4 (±0.14) μm, ω was 0.63–0.88 (±0.08) (at 532 nm) and m ranged from 1.44 (±0.10) + 0.01 (±0.01) i to 1.55 (±0.12) + 0.06 (±0.02) i , in good agreement (only for the r eff values) with in situ aircraft measurements. The water vapor and temperature profiles were incorporated into the ISORROPIA II model to propose a possible in situ aerosol composition consistent with the retrieved m and ω values. The retrieved aerosol chemical composition in the 2–3 km height region gave a variable range of sulfate (0–60%) and organic carbon (OC) content (0–50%), although the OC content increased (up to 50%) and the sulfate content dropped (up to 30%) around 3 km height; the retrieved low ω value (0.63), indicates the presence of absorbing biomass burning smoke mixed with urban haze. Finally, the retrieved aerosol microphysical properties were compared with column-integrated sun photometer CIMEL data.

Description: The RAMNI airborne lidar for cloud and aerosol research Atmospheric Measurement Techniques, 5, 1779-1792, 2012 Author(s): F. Cairo, G. Di Donfrancesco, L. Di Liberto, and M. Viterbini We describe an airborne lidar for the characterization of atmospheric aerosol. The system has been set up in response to the need to monitor extended regions where the air traffic may be posed at risk by the presence of potentially harmful volcanic ash, and to study the characteristics of volcanic emissions both near the source region and when transported over large distances. The lidar provides backscatter and linear depolarization profiles at 532 nm, from which aerosol and cloud properties can be derived. The paper presents the characteristics and capabilities of the lidar system and gives examples of its airborne deployment. Observations from three flights, aimed at assessing the system capabilities in unperturbed atmospheric conditions, and at characterizing the emissions near a volcanic ash source (Mt. Etna) and transported far away from the source, are presented and discussed.

Description: Sensitivity analysis and calibration of a soil carbon model (SoilGen2) in two contrasting loess forest soils Geoscientific Model Development Discussions, 5, 1817-1849, 2012 Author(s): Y. Y. Yu, P. A. Finke, H. B. Wu, and Z. T. Guo To accurately estimate past terrestrial carbon pools is the key to understand the global carbon cycle and its relationship with the climate system. SoilGen2 is a useful tool to obtain aspects of soil properties (including carbon content) by simulating soil formation processes; thus it offers an opportunity for past soil carbon pool reconstruction. In order to apply it to various environmental conditions, parameters related to carbon cycle process in SoilGen2 are calibrated based on 6 soil pedons from two typical loess deposition regions (Belgium and China). Sensitivity analysis using Morris' method shows that decomposition rate of humus ( k HUM ), fraction of incoming plant material as leaf litter (fr ecto ) and decomposition rate of resistant plant material ( k RPM ) are 3 most sensitive parameters that would cause the greatest uncertainty in simulated change of soil organic carbon in both regions. According to the principle of minimizing the difference between simulated and measured organic carbon by comparing quality indices, the suited values of k HUM , fr ecto and k RPM in the model are deduced step by step. The difference of calibrated parameters between Belgium and China may be attributed to their different vegetation types and climate conditions. This calibrated model is improved for better simulation of carbon change in the whole pedon and has potential for future modeling of carbon cycle in paleosols.

Description: Activation of the operational ecohydrodynamic model (3-D CEMBS) – the hydrodynamic part Geoscientific Model Development Discussions, 5, 1851-1875, 2012 Author(s): L. Dzierzbicka-Głowacka, J. Jakacki, M. Janecki, and A. Nowicki The paper presents a description of the hydrodynamic part of the coupled ice-ocean model that also includes ecosystem predictive model for evaluation of the condition of the marine environment and the Baltic ecosystem, as well as a preliminary empirical verification of the operational hydrodynamic model based on the POP code in order to determine the consistence between the results obtained from the model and experimental results for the sea surface temperature. The current Baltic Sea model is based on the Community Earth System Model (CESM from NCAR – National Center for Atmospheric Research). CESM was adopted for the Baltic Sea as a coupled sea-ice model. It consists of the Community Ice Code (CICE model, version 4.0) and the Parallel Ocean Program (POP, version 2.1). The models are coupled through the coupler (CPL7), which is based on the Model Coupling Toolkit (MCT) routines. The current horizontal resolution is about 2 km (1/48 degrees). The ocean model has 21 vertical levels. The driver time step is 1440 s and it is also coupling the time step. The ocean model time step is about 480 s (8 min). Currently, the model is forced by fields from the European Center for Medium Weather Forecast. In the operational mode, 48-h atmospheric forecasts are used, which are supplied by the UM model of the Interdisciplinary Centre for Mathematical and Computational Modelling of the Warsaw University. The model of the marine ecosystem is the right tool for monitoring the state and bioproductivity of the marine ecosystem and forecasting the physical and ecological changes in the studied basin.

Description: Operational SAR-based sea ice drift monitoring over the Baltic Sea Ocean Science, 8, 473-483, 2012 Author(s): J. Karvonen An algorithm for computing ice drift from pairs of synthetic aperture radar (SAR) images covering a common area has been developed at FMI. The algorithm has been developed based on the C-band SAR data over the Baltic Sea. It is based on phase correlation in two scales (coarse and fine) with some additional constraints. The algorithm has been running operationally in the Baltic Sea from the beginning of 2011, using Radarsat-1 ScanSAR wide mode and Envisat ASAR wide swath mode data. The resulting ice drift fields are publicly available as part of the MyOcean EC project. The SAR-based ice drift vectors have been compared to the drift vectors from drifter buoys in the Baltic Sea during the first operational season, and also these validation results are shown in this paper. Also some navigationally useful sea ice quantities, which can be derived from ice drift vector fields, are presented.

Description: Characterization of atmospheric aerosol in the US Southeast from ground- and space-based measurements over the past decade Atmospheric Measurement Techniques, 5, 1667-1682, 2012 Author(s): E. J. Alston, I. N. Sokolik, and O. V. Kalashnikova This study examines how aerosols measured from the ground and space over the US Southeast change temporally over a regional scale during the past decade. PM 2.5 (particulate matter with aerodynamic diameter 〉2.5 micrometers) data consist of two datasets that represent the measurements that are used for regulatory purposes by the US EPA (Environmental Protection Agency) and continuous measurements used for quickly disseminating air quality information. AOD (aerosol optical depth) data come from three NASA sensors: the MODIS sensors onboard Terra and Aqua satellites and the MISR sensor onboard the Terra satellite. We analyze all available data over the state of Georgia from 2000–2009 of both types of aerosol data. The analysis reveals that during the summer the large metropolitan area of Atlanta has average PM 2.5 concentrations that are 50% more than the remainder of the state. Strong seasonality is detected in both the AOD and PM 2.5 datasets, as evidenced by a threefold increase of AOD from mean winter values to mean summer values, and the increase in PM 2.5 concentrations is almost twofold over the same period. Additionally, there is agreement between MODIS and MISR onboard the Terra satellite during the spring and summer, having correlation coefficients of 0.64 and 0.71, respectively. Monthly anomalies were used to determine the presence of a trend in all considered aerosol datasets. We found negative linear trends for both the monthly AOD anomalies from MODIS onboard Terra and the PM 2.5 datasets, which are statistically significant. Decreasing trends were also found for MISR onboard Terra and MODIS onboard Aqua, but those trends were not statistically significant. The observed decrease in AOD and PM 2.5 concentrations may be indicative of the brightening over the study region during the past decade.

Description: Eddy-covariance flux measurements with a weight-shift microlight aircraft Atmospheric Measurement Techniques, 5, 1699-1717, 2012 Author(s): S. Metzger, W. Junkermann, M. Mauder, F. Beyrich, K. Butterbach-Bahl, H. P. Schmid, and T. Foken The objective of this study is to assess the feasibility and quality of eddy-covariance flux measurements from a weight-shift microlight aircraft (WSMA). Firstly, we investigate the precision of the wind measurement (σ u,v ≤ 0.09 m s −1 , σ w = 0.04 m s −1 ), the lynchpin of flux calculations from aircraft. From here, the smallest resolvable changes in friction velocity (0.02 m s −1 ), and sensible- (5 W m −2 ) and latent (3 W m −2 ) heat flux are estimated. Secondly, a seven-day flight campaign was performed near Lindenberg (Germany). Here we compare measurements of wind, temperature, humidity and respective fluxes between a tall tower and the WSMA. The maximum likelihood functional relationship (MLFR) between tower and WSMA measurements considers the random error in the data, and shows very good agreement of the scalar averages. The MLFRs for standard deviations (SDs, 2–34%) and fluxes (17–21%) indicate higher estimates of the airborne measurements compared to the tower. Considering the 99.5% confidence intervals, the observed differences are not significant, with exception of the temperature SD. The comparison with a large-aperture scintillometer reveals lower sensible heat flux estimates at both tower (−40 to −25%) and WSMA (−25–0%). We relate the observed differences to (i) inconsistencies in the temperature and wind measurement at the tower and (ii) the measurement platforms' differing abilities to capture contributions from non-propagating eddies. These findings encourage the use of WSMA as a low cost and highly versatile flux measurement platform.

Description: Tagged ozone mechanism for MOZART-4, CAM-chem, and other chemical transport models Geoscientific Model Development Discussions, 5, 1949-1985, 2012 Author(s): L. K. Emmons, P. G. Hess, J.-F. Lamarque, and G. G. Pfister A procedure for tagging ozone produced from NO sources through updates to an existing chemical mechanism is described, and results from its implementation in the Model for Ozone and Related chemical Tracers (MOZART-4), a global chemical transport model, are presented. Artificial tracers are added to the mechanism, thus not affecting the standard chemistry. The results are linear in the troposphere, i.e., the sum of ozone from individual tagged sources equals the ozone from all sources to within 3% in zonal mean monthly averages. The stratospheric ozone contribution to the troposphere determined from the difference between total ozone and ozone from all tagged sources is significantly less than estimates using a traditional stratospheric ozone tracer (8 vs 20 ppbv at the surface). The commonly used technique of perturbing NO emissions by 20% in a region to determine its ozone contribution is compared to the tagging technique, showing that the tagged ozone is 2–4 times the ozone contribution that was deduced from perturbing emissions.

Description: Global and regional evaluation of over-land spectral aerosol optical depth retrievals from SeaWiFS Atmospheric Measurement Techniques, 5, 1761-1778, 2012 Author(s): A. M. Sayer, N. C. Hsu, C. Bettenhausen, M.-J. Jeong, B. N. Holben, and J. Zhang This study evaluates a new spectral aerosol optical depth (AOD) dataset derived from Sea-viewing Wide Field-of-view Sensor (SeaWiFS) measurements over land. First, the data are validated against Aerosol Robotic Network (AERONET) direct-sun AOD measurements and found to compare well on a global basis. If only data with the highest quality flag are used, the correlation is 0.86 and 72% of matchups fall within an expected absolute uncertainty of 0.05 + 20% (for the wavelength of 550 nm). The quality is similar at other wavelengths and stable over the 13-yr (1997–2010) mission length. Performance tends to be better over vegetated, low-lying terrain with typical AOD of 0.3 or less, such as found over much of North America and Eurasia. Performance tends to be poorer for low-AOD conditions near backscattering geometries, where SeaWiFS overestimates AOD, or optically-thick cases of absorbing aerosol, where SeaWiFS tends to underestimate AOD. Second, the SeaWiFS data are compared with midvisible AOD derived from the Moderate Resolution Imaging Spectrometer (MODIS) and Multiangle Imaging Spectroradiometer (MISR). All instruments show similar spatial and seasonal distributions of AOD, although there are regional and seasonal offsets between them. At locations where AERONET data are available, these offsets are largely consistent with the known validation characteristics of each dataset. With the results of this study in mind, the SeaWiFS over-land AOD record is suitable for quantitative scientific use.

Description: Coupling technologies for Earth System Modelling Geoscientific Model Development Discussions, 5, 1987-2006, 2012 Author(s): S. Valcke, V. Balaji, A. Craig, C. DeLuca, R. Dunlap, R. W. Ford, R. Jacob, J. Larson, R. O'Kuinghttons, G. D. Riley, and M. Vertenstein This paper presents a review of the software currently used in climate modelling in general and in CMIP5 in particular to couple the numerical codes representing the different components of the Earth system. The coupling technologies presented show common features, such as the ability to communicate and regrid data, but also offer different functions and implementations. Design characteristics of the different approaches are discussed as well as future challenges arising from the increasing complexity of scientific problems and computing platforms.

Description: Downscale cascades in tracer transport test cases: an intercomparison of the dynamical cores in the Community Atmosphere Model CAM5 Geoscientific Model Development Discussions, 5, 1781-1816, 2012 Author(s): J. Kent, C. Jablonowski, J. P. Whitehead, and R. B. Rood The accurate modelling of cascades to unresolved scales is an important part of the tracer transport component of dynamical cores of weather and climate models. This paper aims to investigate the ability of the advection schemes in the National Center for Atmospheric Research's Community Atmosphere Model version 5 (CAM5) to model this cascade. In order to quantify the effects of the different advection schemes in CAM5, four two-dimensional tracer transport test cases are presented. Three of the tests stretch the tracer below the scale of coarse resolution grids to ensure the downscale cascade of tracer variance. These results are compared with a high resolution reference solution, which is simulated on a resolution fine enough to resolve the tracer during the test. The fourth test has two separate flow cells, and is designed so that any tracer in the Western Hemisphere should not pass into the Eastern Hemisphere. This is to test whether the diffusion in transport schemes, often in the form of explicit hyper-diffusion terms or implicit through monotonic limiters, contains unphysical mixing. An intercomparison of three of the dynamical cores of the National Center for Atmospheric Research's Community Atmosphere Model version 5 is performed. The results show that the finite-volume (CAM-FV) and spectral element (CAM-SE) dynamical cores model the downscale cascade of tracer variance better than the semi-Lagrangian transport scheme of the Eulerian spectral transform core (CAM-EUL). Each scheme tested produces unphysical mass in the Eastern Hemisphere of the separate cells test.

Description: Critical surface albedo and its implications to aerosol remote sensing Atmospheric Measurement Techniques, 5, 1653-1665, 2012 Author(s): F. C. Seidel and C. Popp We analyse the critical surface albedo (CSA) and its implications to aerosol remote sensing. CSA is defined as the surface albedo where the reflectance at top-of-atmosphere (TOA) does not depend on aerosol optical depth (AOD). AOD retrievals are therefore inaccurate at the CSA. The CSA is obtained by derivatives of the TOA reflectance with respect to AOD using a radiative transfer code. We present the CSA and the effect of surface albedo uncertainties on AOD retrieval and atmospheric correction as a function of aerosol single-scattering albedo, illumination and observation geometry, wavelength and AOD. In general, increasing aerosol absorption and increasing scattering angles lead to lower CSA. In contrast to the strict definition of the CSA, we show that the CSA can also slightly depend on AOD and therefore rather represent a small range of surface albedo values. This was often neglected in previous studies. The following implications to aerosol remote sensing applications were found: (i) surface albedo uncertainties result in large AOD retrieval errors, particularly close to the CSA; (ii) AOD retrievals of weakly or non-absorbing aerosols require dark surfaces, while strongly absorbing aerosols can be retrieved more accurately over bright surfaces; (iii) the CSA may help to estimate aerosol absorption; and (iv) the presented sensitivity of the reflectance at TOA to AOD provides error estimations to optimise AOD retrieval algorithms.

Description: Impact of SLA assimilation in the Sicily Channel Regional Model: model skills and mesoscale features Ocean Science, 8, 485-496, 2012 Author(s): A. Olita, S. Dobricic, A. Ribotti, L. Fazioli, A. Cucco, C. Dufau, and R. Sorgente The impact of the assimilation of MyOcean sea level anomalies along-track data on the analyses of the Sicily Channel Regional Model was studied. The numerical model has a resolution of 1/32° degrees and is capable to reproduce mesoscale and sub-mesoscale features. The impact of the SLA assimilation is studied by comparing a simulation (SIM, which does not assimilate data) with an analysis (AN) assimilating SLA along-track multi-mission data produced in the framework of MyOcean project. The quality of the analysis was evaluated by computing RMSE of the misfits between analysis background and observations (sea level) before assimilation. A qualitative evaluation of the ability of the analyses to reproduce mesoscale structures is accomplished by comparing model results with ocean colour and SST satellite data, able to detect such features on the ocean surface. CTD profiles allowed to evaluate the impact of the SLA assimilation along the water column. We found a significant improvement for AN solution in terms of SLA RMSE with respect to SIM (the averaged RMSE of AN SLA misfits over 2 years is about 0.5 cm smaller than SIM). Comparison with CTD data shows a questionable improvement produced by the assimilation process in terms of vertical features: AN is better in temperature while for salinity it gets worse than SIM at the surface. This suggests that a better a-priori description of the vertical error covariances would be desirable. The qualitative comparison of simulation and analyses with synoptic satellite independent data proves that SLA assimilation allows to correctly reproduce some dynamical features (above all the circulation in the Ionian portion of the domain) and mesoscale structures otherwise misplaced or neglected by SIM. Such mesoscale changes also infer that the eddy momentum fluxes (i.e. Reynolds stresses) show major changes in the Ionian area. Changes in Reynolds stresses reflect a different pumping of eastward momentum from the eddy to the mean flow, in turn influencing transports through the channel.