Assumption Of Horizontal Homogeneity Research Articles

Comparing triple and single Doppler lidar wind measurements with sonic anemometer data based on a new filter strategy for virtual tower measurements

Abstract. In this study, we compare the wind measurements of a virtual tower triple Doppler lidar setup to those of a sonic anemometer located at a height of 90 m above ground on an instrumented tower and with those of two single Doppler lidars to evaluate the effect of the horizontal homogeneity assumption used for single Doppler lidar applications on the measurement accuracy. The triple lidar setup was operated in a 90 m stare and a step–stare mode at six heights between 90 and 500 m above ground, while the single lidars were operated in a continuous scan velocity–azimuth display (VAD) mode where one of them had a zenith angle of 54.7° and the other one of 28.0°. The instruments were set up at the boundary-layer field site of the German Meteorological Service (DWD) in July and August of 2020 during the FESST@MOL (Field Experiment on sub-mesoscale spatiotemporal variability at the Meteorological Observatory Lindenberg) 2020 campaign. Overall, we found good agreement of the lidar methods for the whole study period for different averaging times and scan modes compared to the sonic anemometer. For the step–stare mode wind speed measurements, the comparability between the triple lidar and the sonic anemometer was 0.47 m s−1 at an averaging time of 30 min with a bias value of −0.34 m s−1. For wind speed measured by one single lidar setup for the same period with an averaging time of 30 min, we found a comparability of 0.32 m s−1 at an averaging time of 30 min and a bias value of −0.07 m s−1 as well as values of 0.47 and −0.34 m s−1 for the other one, respectively. We also compared the wind velocity measurements of the single and triple lidars at different heights and found decreasing agreement between them with increasing measurement height up to 495 m above ground for the single lidar systems. We found that the single Doppler lidar with the increased zenith angle produced poorer agreement with the triple Doppler lidar setup than the one with the lower zenith angle, especially at higher altitudes. At a height of 495 m above ground and with an averaging time of 30 min the comparability and bias for the larger zenith angle were 0.71 and −0.50 m s−1, respectively, compared to values of 0.57 and −0.28 m s−1 for the smaller zenith angle. Our results confirm that a single Doppler lidar provides reliable wind speed and direction data over heterogeneous but basically flat terrain in different scan configurations. For the virtual tower scanning strategies, we developed a new filtering approach based on a median absolute deviation (MAD) filter combined with a relatively relaxed filtering criterion for the signal-to-noise ratio output by the instrument.

Comparing triple and single Doppler lidar wind measurements with sonic anemometer data based on a new filter strategy for virtual tower measurements

Abstract. In this study, we compare the wind measurements of a virtual tower triple Doppler lidar setup to those of a sonic anemometer located at a height of 90 m above ground on an instrumented tower and with those of two single Doppler lidars to evaluate the effect of the horizontal homogeneity assumption used for single Doppler lidar applications on the measurement accuracy. The triple lidar setup was operated in a 90 m stare and a step–stare mode at six heights between 90 and 500 m above ground, while the single lidars were operated in a continuous scan velocity–azimuth display (VAD) mode where one of them had a zenith angle of 54.7° and the other one of 28.0°. The instruments were set up at the boundary-layer field site of the German Meteorological Service (DWD) in July and August of 2020 during the FESST@MOL (Field Experiment on sub-mesoscale spatiotemporal variability at the Meteorological Observatory Lindenberg) 2020 campaign. Overall, we found good agreement of the lidar methods for the whole study period for different averaging times and scan modes compared to the sonic anemometer. For the step–stare mode wind speed measurements, the comparability between the triple lidar and the sonic anemometer was 0.47 m s−1 at an averaging time of 30 min with a bias value of −0.34 m s−1. For wind speed measured by one single lidar setup for the same period with an averaging time of 30 min, we found a comparability of 0.32 m s−1 at an averaging time of 30 min and a bias value of −0.07 m s−1 as well as values of 0.47 and −0.34 m s−1 for the other one, respectively. We also compared the wind velocity measurements of the single and triple lidars at different heights and found decreasing agreement between them with increasing measurement height up to 495 m above ground for the single lidar systems. We found that the single Doppler lidar with the increased zenith angle produced poorer agreement with the triple Doppler lidar setup than the one with the lower zenith angle, especially at higher altitudes. At a height of 495 m above ground and with an averaging time of 30 min the comparability and bias for the larger zenith angle were 0.71 and −0.50 m s−1, respectively, compared to values of 0.57 and −0.28 m s−1 for the smaller zenith angle. Our results confirm that a single Doppler lidar provides reliable wind speed and direction data over heterogeneous but basically flat terrain in different scan configurations. For the virtual tower scanning strategies, we developed a new filtering approach based on a median absolute deviation (MAD) filter combined with a relatively relaxed filtering criterion for the signal-to-noise ratio output by the instrument.

A simplified atmospheric boundary layer model for an improved representation of air–sea interactions in eddying oceanic models: implementation and first evaluation in NEMO (4.0)

Abstract. A simplified model of the atmospheric boundary layer (ABL) of intermediate complexity between a bulk parameterization and a three-dimensional atmospheric model is developed and integrated to the Nucleus for European Modelling of the Ocean (NEMO) general circulation model. An objective in the derivation of such a simplified model, called ABL1d, is to reach an apt representation in ocean-only numerical simulations of some of the key processes associated with air–sea interactions at the characteristic scales of the oceanic mesoscale. In this paper we describe the formulation of the ABL1d model and the strategy to constrain this model with large-scale atmospheric data available from reanalysis or real-time forecasts. A particular emphasis is on the appropriate choice and calibration of a turbulent closure scheme for the atmospheric boundary layer. This is a key ingredient to properly represent the air–sea interaction processes of interest. We also provide a detailed description of the NEMO-ABL1d coupling infrastructure and its computational efficiency. The resulting simplified model is then tested for several boundary-layer regimes relevant to either ocean–atmosphere or sea-ice–atmosphere coupling. The coupled system is also tested with a realistic 0.25∘ resolution global configuration. The numerical results are evaluated using standard metrics from the literature to quantify the wind–sea-surface-temperature (a.k.a. thermal feedback effect), wind–current (a.k.a. current feedback effect), and ABL–sea-ice couplings. With respect to these metrics, our results show very good agreement with observations and fully coupled ocean–atmosphere models for a computational overhead of about 9 % in terms of elapsed time compared to standard uncoupled simulations. This moderate overhead, largely due to I/O operations, leaves room for further improvement to relax the assumption of horizontal homogeneity behind ABL1d and thus to further improve the realism of the coupling while keeping the flexibility of ocean-only modeling.

A Theoretical Analysis of Mixing Length for Atmospheric Models From Micro to Large Scales

A new mixing length adapted to the constraints of the hectometric-scale gray zone of turbulence for neutral and convective boundary layers is proposed. It combines a mixing length for mesoscale simulations, where the turbulence is fully subgrid and a mixing length for Large-Eddy Simulations, where the coarsest turbulent eddies are explicitly resolved. The mixing length is built for isotropic turbulence schemes, as well as schemes using the horizontal homogeneity assumption. This mixing length is tested over three boundary layer cases: a free convective case, a neutral case and a cold air outbreak case. The later combines turbulence from thermal and dynamical origins as well as presence of clouds. With this new mixing length, the turbulence scheme produces the right proportion between subgrid and resolved turbulent exchanges in Large Eddy Simulations, in the gray zone and at the mesoscale. This opens the way of using a single mixing length whatever the grid mesh of the atmospheric model, the evolution stage or the depth of the boundary layer.

One-dimensional models of radiation transfer in heterogeneous canopies: a review, re-evaluation, and improved model

Abstract. Despite recent advances in the development of detailed plant radiative transfer models, large-scale canopy models generally still rely on simplified one-dimensional (1-D) radiation models based on assumptions of horizontal homogeneity, including dynamic ecosystem models, crop models, and global circulation models. In an attempt to incorporate the effects of vegetation heterogeneity or “clumping” within these simple models, an empirical clumping factor, commonly denoted by the symbol Ω, is often used to effectively reduce the overall leaf area density and/or index value that is fed into the model. While the simplicity of this approach makes it attractive, Ω cannot in general be readily estimated for a particular canopy architecture and instead requires radiation interception data in order to invert for Ω. Numerous simplified geometric models have been previously proposed, but their inherent assumptions are difficult to evaluate due to the challenge of validating heterogeneous canopy models based on field data because of the high uncertainty in radiative flux measurements and geometric inputs. This work provides a critical review of the origin and theory of models for radiation interception in heterogeneous canopies and an objective comparison of their performance. Rather than evaluating their performance using field data, where uncertainty in the measured model inputs and outputs can be comparable to the uncertainty in the model itself, the models were evaluated by comparing against simulated data generated by a three-dimensional leaf-resolving model in which the exact inputs are known. A new model is proposed that generalizes existing theory and is shown to perform very well across a wide range of canopy types and ground cover fractions.

Detecting wind turbine wakes with nacelle lidars

Because the horizontal homogeneity assumption is violated in wakes flows, lidars face difficulties when reconstructing wind fields. Further, small-scale turbulence which is prevalent in wake flows causes Doppler spectrum widths to be broader than in the free stream. In this study the Doppler peak variance is used as a detection parameter for wakes. A one month long measurement campaign, where a continuous-wave lidar on a turbine has been exposed to multiple wake situations, is used to test the detection capabilities. The results show that it is possible to identify situation where a downstream turbine is in wake by comparing the peak widths. The used lidar is inexpensive and brings instalments on every turbine within economical reach. Thus, the information gathered by the lidars can be used for improved control at wind farm level.

Errors induced by different approximations in handling horizontal atmospheric inhomogeneities in MIPAS/ENVISAT retrievals

Abstract. MIPAS (Michelson Interferometer for Passive Atmospheric Sounding) is a mid-infrared limb emission sounder that operated on board the polar satellite ENVISAT from 2002 to 2012. The retrieval algorithm used by the European Space Agency to process MIPAS measurements exploits the assumption that the atmosphere is horizontally homogeneous. However, previous studies highlighted how this assumption causes significant errors on the retrieved profiles of some MIPAS target species.In this paper we quantify the errors induced by this assumption and evaluate the performances of three different algorithms that can be used to mitigate the problem. We generate synthetic observations with a high spatial resolution atmospheric model and carry out the retrievals with four alternative methods. The first assumes horizontal homogeneity (1-D retrieval), the second includes a model of the horizontal gradient of atmospheric temperature (1-D plus temperature gradient retrieval), the third accounts for an horizontal gradient of temperature and composition (1-D plus temperature and composition gradient retrieval), while the fourth is the full two-dimensional (2-D) inversion approach.Our results highlight that the 1-D retrieval implies errors that are significant for averages of profiles. Furthermore, for some targets (e.g. T, CH4 and N2O below 10 hPa) the error induced by the 1-D approximation also becomes visible in the individual retrieved profiles. The inclusion of any kind of horizontal variability model improves all the targets with respect to the horizontal homogeneity assumption. For temperature, HNO3 and CFC-11, the inclusion of an horizontal temperature gradient leads to a significant reduction of the error. For other targets, such as H2O, O3, N2O, CH4, the improvements due to the inclusion of an horizontal temperature gradient are minor. In these cases, the inclusion of a gradient in the target volume mixing ratio leads to significant improvements. Among all the methods tested in this work, the 2-D approach, as expected, implies the smallest errors for almost all the target parameters. This residual error of the 2-D approach is the smoothing caused by the retrieval grid, which is coarser than that of the atmospheric model.

Accuracy of Rotational and Divergent Kinetic Energy Spectra Diagnosed from Flight-Track Winds

Abstract Under assumptions of horizontal homogeneity and isotropy, one may derive relations between rotational or divergent kinetic energy spectra and velocities along one-dimensional tracks, such as might be measured by aircraft. Two recent studies, differing in details of their implementation, have applied these relations to the Measurement of Ozone and Water Vapor by Airbus In-Service Aircraft (MOZAIC) dataset and reached different conclusions with regard to the mesoscale ratio of divergent to rotational kinetic energy. In this study the accuracy of the method is assessed using global atmospheric simulations performed with the Model for Prediction Across Scales, where the exact decomposition of the horizontal winds into divergent and rotational components may be easily computed. For data from the global simulations, the two approaches yield similar and very accurate results. Errors are largest for the divergent component on synoptic scales, which is shown to be related to a very dominant rotational mode. The errors are, in particular, sufficiently small so that the mesoscale ratio of divergent to rotational kinetic energy can be derived correctly. The proposed technique thus provides a strong observational check of model results with existing large commercial aircraft datasets. The results do, however, show a significant dependence on the height and latitude ranges considered, and the disparate conclusions drawn from previous applications to MOZAIC data may result from the use of different subsets of the data.

Limitations of One-Dimensional Mesoscale PBL Parameterizations in Reproducing Mountain-Wave Flows

Abstract Mesoscale models are considered to be the state of the art in modeling mountain-wave flows. Herein, the authors investigate the role and accuracy of planetary boundary layer (PBL) parameterizations in handling the interaction between large-scale mountain waves and the atmospheric boundary layer. To that end, recent large-eddy simulation (LES) results of mountain waves over a symmetric two-dimensional bell-shaped hill are used and compared to four commonly used PBL schemes. It is found that one-dimensional PBL parameterizations produce reasonable agreement with the LES results in terms of vertical wavelength, amplitude of velocity, and turbulent kinetic energy distribution in the downhill shooting-flow region. However, the assumption of horizontal homogeneity in PBL parameterizations does not hold in the context of these complex flow configurations. This inappropriate modeling assumption results in a vertical wavelength shift, producing errors of approximately 10 m s−1 at downstream locations because of the presence of a coherent trapped lee wave that does not mix with the atmospheric boundary layer. In contrast, horizontally integrated momentum flux derived from these PBL schemes displays a realistic pattern. Therefore, results from mesoscale models using ensembles of one-dimensional PBL schemes can still potentially be used to parameterize drag effects in general circulation models. Nonetheless, three-dimensional PBL schemes must be developed in order for mesoscale models to accurately represent complex terrain and other types of flows where one-dimensional PBL assumptions are violated.

Quantifying error of lidar and sodar Doppler beam swinging measurements of wind turbine wakes using computational fluid dynamics

Abstract. Wind-profiling lidars are now regularly used in boundary-layer meteorology and in applications such as wind energy and air quality. Lidar wind profilers exploit the Doppler shift of laser light backscattered from particulates carried by the wind to measure a line-of-sight (LOS) velocity. The Doppler beam swinging (DBS) technique, used by many commercial systems, considers measurements of this LOS velocity in multiple radial directions in order to estimate horizontal and vertical winds. The method relies on the assumption of homogeneous flow across the region sampled by the beams. Using such a system in inhomogeneous flow, such as wind turbine wakes or complex terrain, will result in errors. To quantify the errors expected from such violation of the assumption of horizontal homogeneity, we simulate inhomogeneous flow in the atmospheric boundary layer, notably stably stratified flow past a wind turbine, with a mean wind speed of 6.5 m s−1 at the turbine hub-height of 80 m. This slightly stable case results in 15° of wind direction change across the turbine rotor disk. The resulting flow field is sampled in the same fashion that a lidar samples the atmosphere with the DBS approach, including the lidar range weighting function, enabling quantification of the error in the DBS observations. The observations from the instruments located upwind have small errors, which are ameliorated with time averaging. However, the downwind observations, particularly within the first two rotor diameters downwind from the wind turbine, suffer from errors due to the heterogeneity of the wind turbine wake. Errors in the stream-wise component of the flow approach 30% of the hub-height inflow wind speed close to the rotor disk. Errors in the cross-stream and vertical velocity components are also significant: cross-stream component errors are on the order of 15% of the hub-height inflow wind speed (1.0 m s−1) and errors in the vertical velocity measurement exceed the actual vertical velocity. By three rotor diameters downwind, DBS-based assessments of wake wind speed deficits based on the stream-wise velocity can be relied on even within the near wake within 1.0 m s−1 (or 15% of the hub-height inflow wind speed), and the cross-stream velocity error is reduced to 8% while vertical velocity estimates are compromised. Measurements of inhomogeneous flow such as wind turbine wakes are susceptible to these errors, and interpretations of field observations should account for this uncertainty.

Position error in profiles retrieved from MIPAS observations with a 1-D algorithm

Abstract. The information load (IL) analysis, first introduced for the two-dimensional approach (Carlotti and Magnani, 2009), is applied to the inversion of MIPAS (Michelson Interferometer for Passive Atmospheric Sounding) observations operated with a 1-dimensional (1-D) retrieval algorithm. The IL distribution of MIPAS spectra is shown to be often asymmetrical with respect to the tangent points of the observations and permits us to define the preferential latitude where the profiles retrieved with a 1-D algorithm should be geo-located. Therefore, defining the geo-location of the retrieved profile by means of the tangent points leads to a "position error". We assess the amplitude of the position error for some of the MIPAS main products and we show that the IL analysis can also be used as a tool for the selection of spectral intervals that, when analyzed, minimize the position error of the retrieved profile. When the temperature (T) profiles are used for the retrieval of volume mixing ratio (VMR) of atmospheric constituents, the T-position error (of the order of 1.5 degrees of latitude) induces a VMR error that is directly connected with the horizontal T gradients. Temperature profiles can be externally-provided or determined in a previous step of the retrieval process. In the first case, the IL analysis shows that a meaningful fraction (often exceeding 50%) of the VMR error deriving from the 1-D approximation is to be attributed to the mismatch between the position assigned to the external T profile and the positions where T is required by the analyzed observations. In the second case the retrieved T values suffer by an error of 1.5–2 K due to neglecting the horizontal variability of T; however the error induced on VMRs is of minor concern because of the generally small mismatch between the IL distribution of the observations analyzed to retrieve T and those analyzed to retrieve the VMR target. An estimate of the contribution of the T-position error to the error budget is provided for MIPAS main products. This study shows that the information load analysis can be successfully exploited in a 1-D context that makes the assumption of horizontal homogeneity of the analyzed portion of atmosphere. The analysis that we propose can be extended to the 1-D inversion of other limb-sounding experiments.

Estimate of the Arctic Convective Boundary Layer Height from Lidar Observations: A Case Study

A new automated small size lidar system (microlidar or MULID) has been developed and employed to perform aerosol measurements since March 2010 at Ny Ålesund (78.9°N,11.9°E), Svalbard. The lidar observations have been used to estimate the PBL height by using the gradient method based on abrupt changes in the vertical aerosol profile and monitor its temporal evolution. The scope of the present study is to compare several approaches to estimate the PBL height, by using lidar observations, meteorological measurements by radio soundings, and a zero-order one-dimensional model based on a parameterization of the turbulent kinetic energy budget within the mixing layer, under the assumptions of horizontal homogeneity, and neglecting radiation and latent heat effects. A case study is presented here for a convective PBL, observed in June 2010 in order to verify whether the Gradient Method can be applied to lidar measurements in the Arctic region to obtain the PBL height. The results obtained are in good agreement with the PBL height estimated by the analysis of thermodynamic measurements obtained from radio sounding and with the model.

Use of a High-Resolution Sodar to Study Surface-layer Turbulence at Night

Measurements in the atmospheric surface layer are generally made with point sensors located in the first few tens of metres. In most cases, however, these measurements are not representative of the whole surface layer. Standard Doppler sodars allow a continuous display of the turbulent thermal structure and wind profiles in the boundary layer up to 1000 m, with a few points, if any, in the surface layer. To overcome these limitations a new sodar configuration is proposed that allows for a higher resolution in the surface layer. Because of its capabilities (echo recording starting at 2 m, echo intensity vertical resolution of approximately 2 m, temporal resolution of 1 s) this sodar is called the surface-layer mini-sodar (SLM-sodar). Features and capabilities of the SLM-sodar are described and compared with the sodar. The comparison of the thermal vertical structure given by the SLM-sodar and the sodar provides evidence that, in most cases, the surface layer presents a level of complexity comparable to that of the entire boundary layer. Considering its high vertical resolution, the SLM-sodar is a promising system for the study of the nocturnal surface layer. The nocturnal SLM-sodar measurements have shown that, depending on wind speed, the structure of the surface layer may change substantially within a short time period. At night, when the wind speed is greater than 3 m s−1, mechanical mixing destroys the wavy structure present in the nocturnal layer. Sonic anemometer measurements have shown that, in such cases, also the sensible heat flux varies with height, reaching a peak in correspondence with the wind speed peak. Under these conditions the assumption of horizontal homogeneity of the surface layer and the choice of the averaging time need to be carefully treated.

Iodine-filter-based mobile Doppler lidar to make continuous and full-azimuth-scanned wind measurements: data acquisition and analysis system, data retrieval methods, and error analysis

An incoherent Doppler wind lidar based on iodine edge filters has been developed at the Ocean University of China for remote measurements of atmospheric wind fields. The lidar is compact enough to fit in a minivan for mobile deployment. With its sophisticated and user-friendly data acquisition and analysis system (DAAS), this lidar has made a variety of line-of-sight (LOS) wind measurements in different operational modes. Through carefully developed data retrieval procedures, various wind products are provided by the lidar, including wind profile, LOS wind velocities in plan position indicator (PPI) and range height indicator (RHI) modes, and sea surface wind. Data are processed and displayed in real time, and continuous wind measurements have been demonstrated for as many as 16 days. Full-azimuth-scanned wind measurements in PPI mode and full-elevation-scanned wind measurements in RHI mode have been achieved with this lidar. The detection range of LOS wind velocity PPI and RHI reaches 8-10 km at night and 6-8 km during daytime with range resolution of 10 m and temporal resolution of 3 min. In this paper, we introduce the DAAS architecture and describe the data retrieval methods for various operation modes. We present the measurement procedures and results of LOS wind velocities in PPI and RHI scans along with wind profiles obtained by Doppler beam swing. The sea surface wind measured for the sailing competition during the 2008 Beijing Olympics is also presented. The precision and accuracy of wind measurements are estimated through analysis of the random errors associated with photon noise and the systematic errors introduced by the assumptions made in data retrieval. The three assumptions of horizontal homogeneity of atmosphere, close-to-zero vertical wind, and uniform sensitivity are made in order to experimentally determine the zero wind ratio and the measurement sensitivity, which are important factors in LOS wind retrieval. Deviations may occur under certain meteorological conditions, leading to bias in these situations. Based on the error analyses and measurement results, we point out the application ranges of this Doppler lidar and propose several paths for future improvement.

An Attempt to Close the Daytime Surface Energy Balance Using Spatially-Averaged Flux Measurements

Single-tower eddy-covariance measurements represent the complete surface flux of a scalar only under idealized conditions. Therefore, we often find an underestimation of energy fluxes expressed as a lack of energy balance closure at many sites. In this study, a multi-tower approach to measure atmospheric energy fluxes based on spatial averaging is evaluated and possible mechanisms causing a lack of energy balance closure are analysed, focussing on daytime data only. It is shown that the multi-tower technique is also unable to measure the entire flux for our site, likely because the assumption of horizontal homogeneity is violated. Heterogeneity-induced and buoyancy-driven quasi-stationary circulations are probably the dominant processes causing the underestimation of energy fluxes. A dependence of the energy balance residual on stability is found, with residuals close to zero for stable stratification, a maximum under unstable to near-neutral conditions and still relatively large residuals for stronger instability. Assuming the processes transporting energy and CO2 are similar, the implications on long-term CO2 flux measurements are analysed. Accordingly, the resulting selective systematic error of cumulative net ecosystem exchange estimates for agricultural regions such as ours can be of the order of more than 100%, since mainly the fluxes during periods of net CO2 uptake are underestimated while periods of net CO2 release are much less affected by this bias. Further investigations about this issue are highly warranted.

Flux Footprints Over an Undulating Surface

The flux footprint probability distribution (FPD) functions for near-surface receptors over an idealised undulating surface are evaluated using a backward Lagrangian stochastic model. The wind and turbulence fields employed to drive the stochastic model are derived from large-eddy simulations, in which the horizontal wind aligns with the surface-elevation-varying direction. The flux FPD for a receptor is affected by flow divergence or convergence, and varies with the receptor’s location. The widest crosswind-integrated FPD (CIFPD) curve with the smallest peak value appears when the receptor is located in the crest area, while the narrowest CIFPD curve with the largest peak value appears when the receptor is located in the windward area. Experiments are designed to highlight the impact of the horizontal homogeneity assumption on the estimation of the FPD. When the receptor is located in the area with surface-wind convergence, the peak value of the CIFPD is larger than its counterpart under assumed horizontally homogeneous flow conditions, with the peak position being closer to the receptor. The case is reversed when the receptor is located in the area with surface-wind divergence. Similar results are obtained when the CIFPD derived from an analytical footprint model (developed under the assumption of horizontally homogeneous flow conditions) is compared with that from the stochastic model over the undulating surface. The analytical model fails to simulate the CIFPD in the local downwind area under weak wind conditions due to the longitudinal wind fluctuation not being considered.

Assimilation of MIPAS limb radiances in the ECMWF system. II: Experiments with a 2‐dimensional observation operator and comparison to retrieval assimilation

AbstractThis study applies and compares three different ways of assimilating MIPAS data in a global numerical weather prediction system. The three methods are: direct assimilation of emitted infrared limb radiances with a 1‐dimensional radiative transfer model that assumes local horizontal homogeneity; direct assimilation of radiances with a 2‐dimensional radiative transfer model that takes into account horizontal gradients in the atmosphere; and assimilation of retrieved profiles of temperature, humidity and ozone, derived under the assumption of horizontal homogeneity. The resulting analyses are compared against each other and against independent retrievals.The use of a 2‐dimensional radiative transfer model in the radiance assimilation leads to relatively small differences in the analyses compared to using a 1‐dimensional observation operator in our experiments. Nevertheless, the results show that the 2‐dimensional operator correctly takes into account the effect of tangent point drift, is capable of extracting a limited amount of horizontal structure from a single MIPAS scan, and leads to smaller first‐guess departures for lower tangent altitudes and more strongly absorbing channels. Consequently, humidity and ozone increments from a 2‐dimensional operator are smaller in the lower stratosphere and upper troposphere in areas where considerable horizontal gradients prevail, and forecasts of humidity and ozone are improved in these regions.In these first trials with assimilation of limb radiances, both the radiance and the retrieval assimilation appear capable of incorporating useful information from MIPAS in the analyses. Both methods introduce broadly similar changes to the mean analyses, with little indication of which method should be favoured. Nevertheless, results from the retrieval assimilation compare better to independent ozone data in the Tropics and over Antarctica. For the radiance assimilation, there are many areas with scope for improvement; these are discussed at the end of the paper. Copyright © 2007 Royal Meteorological Society

A fast radiative‐transfer model for the assimilation of MIPAS limb radiances: Accounting for horizontal gradients

AbstractWe investigate the influence of horizontal structure in the atmosphere on the simulation of emitted infrared limb radiances from the Envisat‐MIPAS instrument, and examine a pragmatic extension of the regression‐based RTMIPAS fast radiative‐transfer model to deal with horizontal atmospheric structure. The extension of RTMIPAS takes into account the horizontal structure at the ray‐tracing step which is used to determine scene temperatures, path lengths, and other conditions along the ray. But the transmittance parametrization used in the extension is based on regression models previously derived from a set of horizontally homogeneous atmospheres. This reduces the computational effort required for the derivation of the regression models. The effects of horizontal gradients in the atmosphere on MIPAS radiances are examined by comparing radiance spectra simulated with RTMIPAS to line‐by‐line equivalents and to real MIPAS observations.It is shown that neglecting horizontal structure in the atmosphere can introduce errors in the radiance simulation that exceed by far the instrument noise, particularly for lower tangent heights and for strongly absorbing spectral regions. The extension of the fast radiative‐transfer model can substantially reduce this error, and for tangent pressures less than 350 hPa (tangent heights greater than about 8 km) almost all error due to horizontal gradients is eliminated. In agreement with this, fast‐model radiances simulated from cross‐sections of the European Centre for Medium‐Range Weather Forecasts model compare significantly better to observed MIPAS limb radiances than radiances simulated under the assumption of horizontal homogeneity. This further confirms a considerable reduction in forward‐model error when horizontal atmospheric structure is taken into account. Copyright © 2006 Royal Meteorological Society

Comparison of high-latitude line-of-sight ozone column density with derived ozone fields and the effects of horizontal inhomogeneity

Abstract. Extensive ozone measurements were made during the second SAGE III Ozone Loss and Validation Experiment (SOLVE II). We compare high-latitude line-of-sight (LOS) slant column ozone measurements from the NASA DC-8 to ozone simulated by forward integration of measurement-derived ozone fields constructed both with and without the assumption of horizontal homogeneity. The average bias and rms error of the simulations assuming homogeneity are relatively small (−6 and 10%, respectively) in comparison to the LOS measurements. The comparison improves significantly (−2% bias; 8% rms error) using forward integrations of three-dimensional proxy ozone fields reconstructed from potential vorticity-O3 correlations. The comparisons provide additional verification of the proxy fields and quantify the influence of large-scale ozone inhomogeneity. The spatial inhomogeneity of the atmosphere is a source of error in the retrieval of trace gas vertical profiles and column abundance from LOS measurements, as well as a complicating factor in intercomparisons that include LOS measurements at large solar zenith angles.

Perturbed profiles of oxygen nightglow emissions as observed by WINDII on UARS

Perturbed profiles having multiple peaks in the O( 1S) green line, O 2(b 1∑ g +) Atmospheric (0–0) band, and OH (8–3) Meinel band P 1(3) line nightglow emissions in the mesopause region are studied using observations of the Wind Imaging Interferometer (WINDII) onboard the Upper Atmosphere Research Satellite (UARS). These profiles were frequently observed by WINDII over the study period of 1991–1997, comprising up to 40% of the observations in bins having certain latitudes and local times. The preferred locations and times for their occurrences coincide with steep emission rate gradients associated with the diurnal tide. The assumption of horizontal homogeneity in the WINDII data processing implies that the perturbed profiles may be caused by horizontal as well as vertical variations. A 1-D model is developed to simulate the horizontal variation of volume emission rate embedded in a Gaussian distributed emission profile. The model results show that the secondary peak in the vertical emission profile could be the reflection of a wave-like horizontal structure with wavelengths between 100 and 400 km. The features of the horizontal structures correspond to localized regions of gravity waves with extents up to 2000 km and lifetimes of one or more days. It is concluded that the WINDII perturbed emission profiles could be produced by horizontal variations induced by gravity waves with wavelengths between 100 and 400 km, and that these occur in preferred regions of steep emission rate gradients associated with the diurnal tide.