The full-spectrum correlated-k method for longwave atmospheric radiative transfer using an effective Planck function

The correlated k-distribution (CKD) method is widely used in the radiative transfer schemes of atmospheric models and involves dividing the spectrum into a number of bands and then reordering the gaseous absorption coefficients within each one. The fluxes and heating rates for each band may then be computed by discretizing the reordered spectrum into of order 10 quadrature points per major gas and performing a monochromatic radiation calculation for each point. In this presentation it is shown that for clear-sky longwave calculations, sufficient accuracy for most applications can be achieved without the need for bands: reordering may be performed on the entire longwave spectrum. The resulting full-spectrum correlated k (FSCK) method requires significantly fewer monochromatic calculations than standard CKD to achieve a given accuracy. The concept is first demonstrated by comparing with line-by-line calculations for an atmosphere containing only water vapor, in which it is shown that the accuracy of heating-rate calculations improves approximately in proportion to the square of the number of quadrature points. For more than around 20 points, the root-mean-squared error flattens out at around 0.015 K/day due to the imperfect rank correlation of absorption spectra at different pressures in the profile. The spectral overlap of m different gases is treated by considering an m-dimensional hypercube where each axis corresponds to the reordered spectrum of one of the gases. This hypercube is then divided up into a number of volumes, each approximated by a single quadrature point, such that the total number of quadrature points is slightly fewer than the sum of the number that would be required to treat each of the gases separately. The gaseous absorptions for each quadrature point are optimized such that they minimize a cost function expressing the deviation of the heating rates and fluxes calculated by the FSCK method from line-by-line calculations for a number of training profiles. This approach is validated for atmospheres containing water vapor, carbon dioxide, and ozone, in which it is found that in the troposphere and most of the stratosphere, heating-rate errors of less than 0.2 K/day can be achieved using a total of 23 quadrature points, decreasing to less than 0.1 K/day for 32 quadrature points. It would be relatively straightforward to extend the method to include other gases

Equitability revisited: why the “equitable threat score” is not equitable

In the forecasting of binary events, verification measures that are “equitable” were defined by Gandin and Murphy to satisfy two requirements: 1) they award all random forecasting systems, including those that always issue the same forecast, the same expected score (typically zero), and 2) they are expressible as the linear weighted sum of the elements of the contingency table, where the weights are independent of the entries in the table, apart from the base rate. The authors demonstrate that the widely used “equitable threat score” (ETS), as well as numerous others, satisfies neither of these requirements and only satisfies the first requirement in the limit of an infinite sample size. Such measures are referred to as “asymptotically equitable.” In the case of ETS, the expected score of a random forecasting system is always positive and only falls below 0.01 when the number of samples is greater than around 30. Two other asymptotically equitable measures are the odds ratio skill score and the symmetric extreme dependency score, which are more strongly inequitable than ETS, particularly for rare events; for example, when the base rate is 2% and the sample size is 1000, random but unbiased forecasting systems yield an expected score of around −0.5, reducing in magnitude to −0.01 or smaller only for sample sizes exceeding 25 000. This presents a problem since these nonlinear measures have other desirable properties, in particular being reliable indicators of skill for rare events (provided that the sample size is large enough). A potential way to reconcile these properties with equitability is to recognize that Gandin and Murphy’s two requirements are independent, and the second can be safely discarded without losing the key advantages of equitability that are embodied in the first. This enables inequitable and asymptotically equitable measures to be scaled to make them equitable, while retaining their nonlinearity and other properties such as being reliable indicators of skill for rare events. It also opens up the possibility of designing new equitable verification measures

New near-IR observations of mesospheric CO2 and H2O clouds on Mars

Carbon dioxide clouds, which are speculated by models on solar and extra-solar planets, have been recently observed near the equator of Mars. The most comprehensive identification of Martian CO2 ice clouds has been obtained by the near-IR imaging spectrometer OMEGA. CRISM, a similar instrument with a higher spatial resolution, cannot detect these clouds with the same method due to its shorter wavelength range. Here we present a new method to detect CO2 clouds using near-IR data based on the comparison of H2O and CO2 ice spectral properties. The spatial and seasonal distributions of 54 CRISM observations containing CO2 clouds are reported, in addition to 17 new OMEGA observations. CRISM CO2 clouds are characterized by grain size in the 0.5-2\mum range and optical depths lower than 0.3. The distributions of CO2 clouds inferred from OMEGA and CRISM are consistent with each other and match at first order the distribution of high altitude (>60km) clouds derived from previous studies. At second order, discrepancies are observed. We report the identification of H2O clouds extending up to 80 km altitude, which could explain part of these discrepancies: both CO2 and H2O clouds can exist at high, mesospheric altitudes. CRISM observations of afternoon CO2 clouds display morphologies resembling terrestrial cirrus, which generalizes a previous result to the whole equatorial clouds season. Finally, we show that morning OMEGA observations have been previously misinterpreted as evidence for cumuliform, and hence potentially convective, CO2 clouds.Comment: Vincendon, M., C. Pilorget, B. Gondet, S. Murchie, and J.-P. Bibring (2011), New near-IR observations of mesospheric CO2 and H2O clouds on Mars, J. Geophys. Res., 116, E00J0

The GALFA-HI Compact Cloud Catalog

We present a catalog of 1964 isolated, compact neutral hydrogen clouds from the Galactic Arecibo L-Band Feed Array Survey Data Release One (GALFA-HI DR1). The clouds were identified by a custom machine-vision algorithm utilizing Difference of Gaussian kernels to search for clouds smaller than 20'. The clouds have velocities typically between |VLSR| = 20-400 km/s, linewidths of 2.5-35 km/s, and column densities ranging from 1 - 35 x 10^18 cm^-2. The distances to the clouds in this catalog may cover several orders of magnitude, so the masses may range from less than a Solar mass for clouds within the Galactic disc, to greater than 10^4 Solar Masses for HVCs at the tip of the Magellanic Stream. To search for trends, we separate the catalog into five populations based on position, velocity, and linewidth: high velocity clouds (HVCs); galaxy candidates; cold low velocity clouds (LVCs); warm, low positive-velocity clouds in the third Galactic Quadrant; and the remaining warm LVCs. The observed HVCs are found to be associated with previously-identified HVC complexes. We do not observe a large population of isolated clouds at high velocities as some models predict. We see evidence for distinct histories at low velocities in detecting populations of clouds corotating with the Galactic disc and a set of clouds that is not corotating.Comment: 34 Pages, 9 Figures, published in ApJ (2012, ApJ, 758, 44), this version has the corrected fluxes and corresponding flux histogram and masse

A Large-scale CO Imaging of the Galactic Center. II. Dynamical Properties of Molecular Clouds

The data from the Nobeyama Radio Observatory 45 m telescope Galactic Center CO survey have been analyzed to generate a compilation of molecular clouds with intense CO emission in this region. Clouds are identified in an automated manner throughout the main part of the survey data for all CO emission peaks exceeding 10 K ($T_R^*$). Correlations between the size, velocity dispersion, virial mass, and the CO luminosity, for the molecular clouds in the Galactic center were shown. We diagnosed gravitational stabilities of identified clouds assuming that the disk clouds are nearly at the onset of gravitational instability. Most of the clouds and cloud complexes in the Galactic center are gravitationally stable, while some clouds with intense CO emission are gravitationally unstable.Comment: 4 pages, 4 figures, to appear in the Proceedings of the 32nd COSPAR Scientific Assembl

A Herschel [C II] Galactic plane survey II: CO-dark H2 in clouds

ABRIDGED: Context: HI and CO large scale surveys of the Milky Way trace the diffuse atomic clouds and the dense shielded regions of molecular hydrogen clouds. However, until recently, we have not had spectrally resolved C+ surveys to characterize the photon dominated interstellar medium, including, the H2 gas without C, the CO-dark H2, in a large sample of clouds. Aims: To use a sparse Galactic plane survey of the 1.9 THz [C II] spectral line from the Herschel Open Time Key Programme, Galactic Observations of Terahertz C+ (GOT C+), to characterize the H2 gas without CO in a statistically significant sample of clouds. Methods: We identify individual clouds in the inner Galaxy by fitting [CII] and CO isotopologue spectra along each line of sight. We combine these with HI spectra, along with excitation models and cloud models of C+, to determine the column densities and fractional mass of CO-dark H2 clouds. Results: We identify 1804 narrow velocity [CII] interstellar cloud components in different categories. About 840 are diffuse molecular clouds with no CO, 510 are transition clouds containing [CII] and 12CO, but no 13CO, and the remainder are dense molecular clouds containing 13CO emission. The CO-dark H2 clouds are concentrated between Galactic radii 3.5 to 7.5 kpc and the column density of the CO-dark H2 layer varies significantly from cloud-to-cloud with an average 9X10^(20) cm-2. These clouds contain a significant fraction of CO-dark H2 mass, varying from ~75% for diffuse molecular clouds to ~20% for dense molecular clouds. Conclusions: We find a significant fraction of the warm molecular ISM gas is invisible in HI and CO, but is detected in [CII]. The fraction of CO-dark H2 is greatest in the diffuse clouds and decreases with increasing total column density, and is lowest in the massive clouds.Comment: 21 pages, 19 figures, accepted for publication in A&A (2014