Scheduling satellite-based SAR acquisition for sequential assimilation of water level observations into flood modelling

Satellite-based Synthetic Aperture Radar (SAR) has proved useful for obtaining information on flood extent, which, when intersected with a Digital Elevation Model (DEM) of the floodplain, provides water level observations that can be assimilated into a hydrodynamic model to decrease forecast uncertainty. With an increasing number of operational satellites with SAR capability, information on the relationship between satellite first visit and revisit times and forecast performance is required to optimise the operational scheduling of satellite imagery. By using an Ensemble Transform Kalman Filter (ETKF) and a synthetic analysis with the 2D hydrodynamic model LISFLOOD-FP based on a real flooding case affecting an urban area (summer 2007,Tewkesbury, Southwest UK), we evaluate the sensitivity of the forecast performance to visit parameters. We emulate a generic hydrologic-hydrodynamic modelling cascade by imposing a bias and spatiotemporal correlations to the inflow error ensemble into the hydrodynamic domain. First, in agreement with previous research, estimation and correction for this bias leads to a clear improvement in keeping the forecast on track. Second, imagery obtained early in the flood is shown to have a large influence on forecast statistics. Revisit interval is most influential for early observations. The results are promising for the future of remote sensing-based water level observations for real-time flood forecasting in complex scenarios.

The tropopause inversion layer in models and analyses

Recent high-resolution radiosonde climatologies have revealed a tropopause inversion layer (TIL) in the extratropics: temperature strongly increases just above a sharp local cold point tropopause. Here, it is asked to what extent a TIL exists in current general circulation models (GCMs) and meteorological analyses. Only a weak hint of a TIL exists in NCEP/NCAR reanalysis data. In contrast, the Canadian Middle Atmosphere Model (CMAM), a comprehensive GCM, exhibits a TIL of realistic strength. However, in data assimilation mode CMAM exhibits a much weaker TIL, especially in the Southern Hemisphere where only coarse satellite data are available. The discrepancy between the analyses and the GCM is thus hypothesized to be mainly due to data assimilation acting to smooth the observed strong curvature in temperature around the tropopause. This is confirmed in the reanalysis where the stratification around the tropopause exhibits a strong discontinuity at the start of the satellite era.

Programmable data acquisition system for research measurements from meteorological radiosondes

A programmable data acquisition system to allow novel use of meteorological radiosondes for atmospheric science measurements is described. In its basic form it supports four analogue inputs at 16 bit resolution, and up to two further inputs at lower resolution configurable instead for digital instruments. It also provides multiple instrument power supplies (+8V, +16V, +5V and -8V) from the 9V radiosonde battery. During a balloon flight encountering air temperatures from +17°C to -66°C, the worst case voltage drift in the 5V unipolar digitisation circuitry was 20mV. The system liberates a new range of low cost atmospheric research measurements, by utilising radiosondes routinely launched internationally for weather forecasting purposes. No additional receiving equipment is required. Comparisons between the specially instrumented and standard meteorological radiosondes show negligible effect of the additional instrumentation on the standard meteorological data.

The role of baroclinic waves in the initiation of tropical cyclones across the southern Indian Ocean

Cases where tropical storms are initiated simultaneously along one latitude are investigated. It is argued that such structure arises as part of a baroclinic wave. A case from February 2008 is examined using European Centre for Medium-Range Forecasts (ECMWF) analyses; the birth of three tropical cyclones in the low-level cyclonic regions to the east of upper-level troughs suggests that the wave was instrumental for initiation. Archived satellite imagery and storm warnings reveal that baroclinic waves over the southern Indian Ocean accompany tropical cyclogenesis twice a season on average, mainly in late summer, when breaking Rossby waves on the subtropical westerly jet are closest to the Intertropical Convergence Zone (ITCZ). Copyright © 2012 Royal Meteorological Society

The European storm Kyrill in January 2007: synoptic evolution, meteorological impacts and some considerations with respect to climate change

The synoptic evolution and some meteorological impacts of the European winter storm Kyrill that swept across Western, Central, and Eastern Europe between 17 and 19 January 2007 are investigated. The intensity and large storm damage associated with Kyrill is explained based on synoptic and mesoscale environmental storm features, as well as on comparisons to previous storms. Kyrill appeared on weather maps over the US state of Arkansas about four days before it hit Europe. It underwent an explosive intensification over the Western North Atlantic Ocean while crossing a very intense zonal polar jet stream. A superposition of several favourable meteorological conditions west of the British Isles caused a further deepening of the storm when it started to affect Western Europe. Evidence is provided that a favourable alignment of three polar jet streaks and a dry air intrusion over the occlusion and cold fronts were causal factors in maintaining Kyrill's low pressure very far into Eastern Europe. Kyrill, like many other strong European winter storms, was embedded in a pre-existing, anomalously wide, north-south mean sea-level pressure (MSLP) gradient field. In addition to the range of gusts that might be expected from the synoptic-scale pressure field, mesoscale features associated with convective overturning at the cold front are suggested as the likely causes for the extremely damaging peak gusts observed at many lowland stations during the passage of Kyrill's cold front. Compared to other storms, Kyrill was by far not the most intense system in terms of core pressure and circulation anomaly. However, the system moved into a pre-existing strong MSLP gradient located over Central Europe which extended into Eastern Europe. This fact is considered determinant for the anomalously large area affected by Kyrill. Additionally, considerations of windiness in climate change simulations using two state-of-the-art regional climate models driven by ECHAM5 indicate that not only Central, but also Eastern Central Europe may be affected by higher surface wind speeds at the end of the 21st century. These changes are partially associated with the increased pressure gradient over Europe which is identified in the ECHAM5 simulations. Thus, with respect to the area affected, as well as to the synoptic and mesoscale storm features, it is proposed that Kyrill may serve as an interesting study case to assess future storm impacts.

On the 'downward control' of extratropical diabatic circulations by eddy-induced mean zonal forces

The situation considered is that of a zonally symmetric model of the middle atmosphere subject to a given quasi-steady zonal force F̄, conceived to be the result of irreversible angular momentum transfer due to the upward propagation and breaking of Rossby and gravity waves together with any other dissipative eddy effects that may be relevant. The model's diabatic heating is assumed to have the qualitative character of a relaxation toward some radiatively determined temperature field. To the extent that the force F̄ may be regarded as given, and the extratropical angular momentum distribution is realistic, the extratropical diabatic mass flow across a given isentropic surface may be regarded as controlled exclusively by the F̄ distribution above that surface (implying control by the eddy dissipation above that surface and not, for instance, by the frequency of tropopause folding below). This “downward control” principle expresses a critical part of the dynamical chain of cause and effect governing the average rate at which photochemical products like ozone become available for folding into, or otherwise descending into, the extratropical troposphere. The dynamical facts expressed by the principle are also relevant, for instance, to understanding the seasonal-mean rate of upwelling of water vapor to the summer mesopause, and the interhemispheric differences in stratospheric tracer transport. The robustness of the principle is examined when F̄ is time-dependent. For a global-scale, zonally symmetric diabatic circulation with a Brewer-Dobson-like horizontal structure given by the second zonally symmetric Hough mode, with Rossby height HR = 13 km in an isothermal atmosphere with density scale height H = 7 km, the vertical partitioning of the unsteady part of the mass circulation caused by fluctuations in F̄ confined to a shallow layer LF̄ is always at least 84% downward. It is 90% downward when the force fluctuates sinusoidally on twice the radiative relaxation timescale and 95% if five times slower. The time-dependent adjustment when F̄ is changed suddenly is elucidated, extending the work of Dickinson (1968), when the atmosphere is unbounded above and below. Above the forcing, the adjustment is characterized by decay of the meridional mass circulation cell at a rate proportional to the radiative relaxation rate τr−1 divided by {1 + (4H2/HR2)}. This decay is related to the boundedness of the angular momentum that can be taken up by the finite mass of air above LF̄ without causing an ever-increasing departure from thermal wind balance. Below the forcing, the meridional mass circulation cell penetrates downward at a speed τr−1 HR2/H. For the second Hough mode, the time for downward penetration through one density scale height is about 6 days if the radiative relaxation time is 20 days, the latter being representative of the lower stratosphere. At any given altitude, a steady state is approached. The effect of a rigid lower boundary on the time-dependent adjustment is also considered. If a frictional planetary boundary layer is present then a steady state is ultimately approached everywhere, with the mass circulation extending downward from LF̄ and closing via the boundary layer. Satellite observations of temperature and ozone are used in conjunction with a radiative transfer scheme to estimate the altitudes from which the lower stratospheric diabatic vertical velocity is controlled by the effective F̄ in the real atmosphere. The data appear to indicate that about 80% of the effective control is usually exerted from below 40 km but with significant exceptions up to 70 km (in the high latitude southern hemispheric winter). The implications for numerical modelling of chemical transport are noted.

Jet characterization in the upper troposphere/lower stratosphere (UTLS): applications to climatology and transport studies

A method of classifying the upper tropospheric/lower stratospheric (UTLS) jets has been developed that allows satellite and aircraft trace gas data and meteorological fields to be efficiently mapped in a jet coordinate view. A detailed characterization of multiple tropopauses accompanies the jet characterization. Jet climatologies show the well-known high altitude subtropical and lower altitude polar jets in the upper troposphere, as well as a pattern of concentric polar and subtropical jets in the Southern Hemisphere, and shifts of the primary jet to high latitudes associated with blocking ridges in Northern Hemisphere winter. The jet-coordinate view segregates air masses differently than the commonly-used equivalent latitude (EqL) coordinate throughout the lowermost stratosphere and in the upper troposphere. Mapping O3 data from the Aura Microwave Limb Sounder (MLS) satellite and the Winter Storms aircraft datasets in jet coordinates thus emphasizes different aspects of the circulation compared to an EqL-coordinate framework: the jet coordinate reorders the data geometrically, thus highlighting the strong PV, tropopause height and trace gas gradients across the subtropical jet, whereas EqL is a dynamical coordinate that may blur these spatial relationships but provides information on irreversible transport. The jet coordinate view identifies the concentration of stratospheric ozone well below the tropopause in the region poleward of and below the jet core, as well as other transport features associated with the upper tropospheric jets. Using the jet information in EqL coordinates allows us to study trace gas distributions in regions of weak versus strong jets, and demonstrates weaker transport barriers in regions with less jet influence. MLS and Atmospheric Chemistry Experiment-Fourier Transform Spectrometer trace gas fields for spring 2008 in jet coordinates show very strong, closely correlated, PV, tropopause height and trace gas gradients across the jet, and evidence of intrusions of stratospheric air below the tropopause below and poleward of the subtropical jet; these features are consistent between instruments and among multiple trace gases. Our characterization of the jets is facilitating studies that will improve our understanding of upper tropospheric trace gas evolution.

Assessment of intercalibration methods for satellite microwave humidity sounders

Three methods for intercalibrating humidity sounding channels are compared to assess their merits and demerits. The methods use the following: (1) natural targets (Antarctica and tropical oceans), (2) zonal average brightness temperatures, and (3) simultaneous nadir overpasses (SNOs). Advanced Microwave Sounding Unit-B instruments onboard the polar-orbiting NOAA 15 and NOAA 16 satellites are used as examples. Antarctica is shown to be useful for identifying some of the instrument problems but less promising for intercalibrating humidity sounders due to the large diurnal variations there. Owing to smaller diurnal cycles over tropical oceans, these are found to be a good target for estimating intersatellite biases. Estimated biases are more resistant to diurnal differences when data from ascending and descending passes are combined. Biases estimated from zonal-averaged brightness temperatures show large seasonal and latitude dependence which could have resulted from diurnal cycle aliasing and scene-radiance dependence of the biases. This method may not be the best for channels with significant surface contributions. We have also tested the impact of clouds on the estimated biases and found that it is not significant, at least for tropical ocean estimates. Biases estimated from SNOs are the least influenced by diurnal cycle aliasing and cloud impacts. However, SNOs cover only relatively small part of the dynamic range of observed brightness temperatures.

Improving a convection-permitting model simulation of a cold air outbreak

A convection-permitting local-area model was used to simulate a cold air outbreak crossing from the Norwegian Sea into the Atlantic Ocean near Scotland. A control model run based on an operational configuration of the Met Office UKV high-resolution (1.5 km grid spacing) NWP model was compared to satellite, aircraft and radar data. While the control model captured the large-scale features of the synoptic situation, it was not able to reproduce the shallow (<1.5 km) stratiform layer to the north of the open cellular convection. Liquid water paths were found to be too low in both the stratiform and convective cloud regions. Sensitivity analyses including a modified boundary-layer diagnosis to generate a more well-mixed boundary layer and inhibition of ice formation to lower temperatures improved cloud morphology and comparisons with observational data. Copyright © 2013 Royal Meteorological Society and British Crown Copyright, the Met Office

The global distribution of phytoplankton size spectrum and size classes from their light-absorption spectra derived from satellite data

The absorption spectra of phytoplankton in the visible domain hold implicit information on the phytoplankton community structure. Here we use this information to retrieve quantitative information on phytoplankton size structure by developing a novel method to compute the exponent of an assumed power-law for their particle-size spectrum. This quantity, in combination with total chlorophyll-a concentration, can be used to estimate the fractional concentration of chlorophyll in any arbitrarily-defined size class of phytoplankton. We further define and derive expressions for two distinct measures of cell size of mixed populations, namely, the average spherical diameter of a bio-optically equivalent homogeneous population of cells of equal size, and the average equivalent spherical diameter of a population of cells that follow a power-law particle-size distribution. The method relies on measurements of two quantities of a phytoplankton sample: the concentration of chlorophyll-a, which is an operational index of phytoplankton biomass, and the total absorption coefficient of phytoplankton in the red peak of visible spectrum at 676 nm. A sensitivity analysis confirms that the relative errors in the estimates of the exponent of particle size spectra are reasonably low. The exponents of phytoplankton size spectra, estimated for a large set of in situ data from a variety of oceanic environments (~ 2400 samples), are within a reasonable range; and the estimated fractions of chlorophyll in pico-, nano- and micro-phytoplankton are generally consistent with those obtained by an independent, indirect method based on diagnostic pigments determined using high-performance liquid chromatography. The estimates of cell size for in situ samples dominated by different phytoplankton types (diatoms, prymnesiophytes, Prochlorococcus, other cyanobacteria and green algae) yield nominal sizes consistent with the taxonomic classification. To estimate the same quantities from satellite-derived ocean-colour data, we combine our method with algorithms for obtaining inherent optical properties from remote sensing. The spatial distribution of the size-spectrum exponent and the chlorophyll fractions of pico-, nano- and micro-phytoplankton estimated from satellite remote sensing are in agreement with the current understanding of the biogeography of phytoplankton functional types in the global oceans. This study contributes to our understanding of the distribution and time evolution of phytoplankton size structure in the global oceans.

Modelling the time-evolution of phytoplankton size spectra from satellite remote sensing

A dynamic size-structured model is developed for phytoplankton and nutrients in the oceanic mixed layer and applied to extract phytoplankton biomass at discrete size fractions from remotely sensed, ocean-colour data. General relationships between cell size and biophysical processes (such as sinking, grazing, and primary production) of phytoplankton were included in the model through a bottom–up approach. Time-dependent, mixed-layer depth was used as a forcing variable, and a sequential data-assimilation scheme was implemented to derive model trajectories. From a given time-series, the method produces estimates of size-structured biomass at every observation, so estimates seasonal succession of individual phytoplankton size, derived here from remote sensing for the first time. From these estimates, normalized phytoplankton biomass size spectra over a period of 9 years were calculated for one location in the North Atlantic. Further analysis demonstrated that strong relationships exist between the seasonal trends of the estimated size spectra and the mixed-layer depth, nutrient biomass, and total chlorophyll. The results contain useful information on the time-dependent biomass flux in the pelagic ecosystem.

Generalized Bayesian cloud detection for satellite imagery. part 1: technique and validation for night-time imagery over land and sea

Numerical Weather Prediction (NWP) fields are used to assist the detection of cloud in satellite imagery. Simulated observations based on NWP are used within a framework based on Bayes' theorem to calculate a physically-based probability of each pixel with an imaged scene being clear or cloudy. Different thresholds can be set on the probabilities to create application-specific cloud-masks. Here, this is done over both land and ocean using night-time (infrared) imagery. We use a validation dataset of difficult cloud detection targets for the Spinning Enhanced Visible and Infrared Imager (SEVIRI) achieving true skill scores of 87% and 48% for ocean and land, respectively using the Bayesian technique, compared to 74% and 39%, respectively for the threshold-based techniques associated with the validation dataset.

Generalised Bayesian cloud detection for satellite imagery. part 2: technique and validation for day-time imagery

Numerical Weather Prediction (NWP) fields are used to assist the detection of cloud in satellite imagery. Simulated observations based on NWP are used within a framework based on Bayes' theorem to calculate a physically-based probability of each pixel with an imaged scene being clear or cloudy. Different thresholds can be set on the probabilities to create application-specific cloud masks. Here, the technique is shown to be suitable for daytime applications over land and sea, using visible and near-infrared imagery, in addition to thermal infrared. We use a validation dataset of difficult cloud detection targets for the Spinning Enhanced Visible and Infrared Imager (SEVIRI) achieving true skill scores of 89% and 73% for ocean and land, respectively using the Bayesian technique, compared to 90% and 70%, respectively for the threshold-based techniques associated with the validation dataset.

Cloud optical and microphysical properties derived from ground-based and satellite sensors over a site in the Yangtze Delta region

Comprehensive surface-based retrievals of cloud optical and microphysical properties were made at Taihu, a highly polluted site in the central Yangtze Delta region, during a research campaign from May 2008 to December 2009. Cloud optical depth (COD), effective radius (Re), and liquid water path (LWP) were retrieved from measurements made with a suite of ground-based and spaceborne instruments, including an Analytical Spectral Devices spectroradiometer, a multi␣lter rotating shadowband radiometer, a multichannel microwave radiometer profiler, and the Moderate Resolution Imaging Spectroradiometer (MODIS) on Terra and Aqua satellites. Retrievals from zenith radiance measurements capture better the temporal variation of cloud properties than do retrievals from hemispherical fluxes. Annual mean LWP, COD, and Re are 115.8 ± 90.8 g/m2, 28.5 ± 19.2, and 6.9 ± 4.2 microns. Over 90% of LWP values are less than 250 g/m2. Most of the COD values (>90%) fall between 5 and 60, and ~80% of Re values are less than 10 microns. Maximum (minimum) values of LWP and Re occur in summer (winter); COD is highest in winter and spring. Raining and nonraining clouds have signi␣cant differences in LWP, COD, and Re. Rainfall frequency is best correlated with LWP, followed by COD and Re. Cloud properties retrieved from multiple ground-based instruments are also compared with those from satellite retrievals. On average, relative to surface retrievals, mean differences of satellite retrievals in cloud LWP, COD, and Re were -33.6 g/m2 (-26.4%), -5.8 (-31.4%), and 2.9 ␣m (29.3%) for 11 MODIS-Terra overpasses and -43.3 g/m2 (-22.3%), -3.0 (-10.0%), and -1.3 ␣m (-12.0%) for 8 MODIS-Aqua overpasses, respectively. These discrepancies indicate that MODIS cloud products still suffer from large uncertainties in this region.