The Correlation Structure of Matheron's Classical Variogram Estimator Under Elliptically Contoured Distributions

The classical variogram estimator proposed by Matheron can be written as a quadratic form of the observations. When data have an elliptically contoured distribution with constant mean, the correlation between the classical variogram estimator at two different lags is a function of the spatial design matrix, the covariance matrix, and the kurtosis. Several specific cases are studied closely. A subclass of elliptically contoured distributions with a particular family of covariance matrices is shown to possess exactly the same correlation structure for the classical variogram estimator as the multivariate independent Gaussian distribution. The consequences on variogram fitting by generalized least squares are discussed.     

Park Forest (L5) and the asteroidal source of shocked L chondrites

The Park Forest (L5) meteorite fell in a suburb of Chicago, Illinois (USA) on March 26, 2003. It is one of the currently 25 meteorites for which photographic documentation of the fireball enabled the reconstruction of the meteoroid orbit. The combination of orbits with pre‐atmospheric sizes, cosmic‐ray exposure (CRE), and radiogenic gas retention ages (“cosmic histories”) is significant because they can be used to constrain the meteoroid's “birth region,” and test models of meteoroid delivery. Using He, Ne, Ar, ¹⁰Be, and ²⁶Al, as well as a dynamical model, we show that the Park Forest meteoroid had a pre‐atmospheric size close to 180 g cm^?2, 0–40% porosity, and a pre‐atmospheric mass range of ~2–6 tons. It has a CRE age of 14 ± 2 Ma, and (U, Th)‐He and K‐Ar ages of 430 ± 90 and 490 ± 70 Ma, respectively. Of the meteorites with photographic orbits, Park Forest is the second (after Novato) that was shocked during the L chondrite parent body (LCPB) break‐up event approximately 470 Ma ago. The suggested association of this event with the formation of the Gefion family of asteroids has recently been challenged and we suggest the Ino family as a potential alternative source for the shocked L chondrites. The location of the LCPB break‐up event close to the 5:2 resonance also allows us to put some constraints on the possible orbital migration paths of the Park Forest meteoroid.     

Scientific inference and experiment in Ecosystem Based Fishery Management,with application to Steller sea lions in the Bering Sea and Western Gulf of Alaska

Learning about ecosystem processes and patterns is an essential component of Ecosystem Based Fishery Management and the sustainable use of natural resources.     

Remote sensing-based estimation of rice yields using various models:A critical review

Reliable estimation of region-wide rice yield is vital for food security and agricultural management. Field-scale models have increased our understanding of ric...     

Star-like activity from a very young 'isolated planet'

The discovery of isolated bodies of planetary mass has challenged the paradigm that planets form only as companions to stars. To determine whether 'isolated planets', brown dwarfs and stars can have a common origin, we have made deep submillimetre observations of part of the ρ Oph B star formation region. Spectroscopy of the 9-Jupiter-mass core Oph B-11 has revealed carbon monoxide line wings such as those of a protostar. Moreover, the estimated mass of outflowing gas lies on the force versus core-mass relation for protostars and protobrown dwarfs. This is evidence for a common process that can form any object between planetary and stellar masses in a molecular cloud. In a submillimetre continuum map, six compact cores in ρ Oph B were found to have masses presently below the deuterium-burning limit, extending the core mass function down to  0.01 M_⊙  with the approximate form  d N /d M ∝ M ^−3/2  . If these lowest-mass cores are not transient and can collapse under gravity, then isolated planets should be very common in ρ Oph in the future, as is the case in the Orion star formation region. In fact, the isolated planetary objects that may form from these cores would outnumber the massive planets that have been found as companions to stars.     

Dynamics of photosynthetic photon flux density (PPFD) and estimates in coastal northern California

Plants require solar radiation for photosynthesis and their growth is directly related to the amount received, assuming that other environmental parameters are not limiting. Therefore, precise estimation of photosynthetically active radiation (PAR) is necessary to enhance overall accuracies of plant growth models. This study aimed to explore the PAR radiant flux in the San Francisco Bay Area of northern California. During the growing season (March through August) for 2?years 2007?C2008, the on-site magnitudes of photosynthetic photon flux densities (PPFD) were investigated and then processed at both the hourly and daily time scales. Combined with global solar radiation (R _S) and simulated extraterrestrial solar radiation, five PAR-related values were developed, i.e., flux density-based PAR (PPFD), energy-based PAR (PARE), from-flux-to-energy conversion efficiency (fFEC), and the fraction of PAR energy in the global solar radiation (fE), and a new developed indicator??lost PARE percentages (LPR)??when solar radiation penetrates from the extraterrestrial system to the ground. These PAR-related values indicated significant diurnal variation, high values occurring at midday, with the low values occurring in the morning and afternoon hours. During the entire experimental season, the overall mean hourly value of fFEC was found to be 2.17???mol?J^?1, while the respective fE value was 0.49. The monthly averages of hourly fFEC and fE at the solar noon time ranged from 2.15 in March to 2.39???mol?J^?1 in August and from 0.47 in March to 0.52 in July, respectively. However, the monthly average daily values were relatively constant, and they exhibited a weak seasonal variation, ranging from 2.02?mol?MJ^?1 and 0.45 (March) to 2.19?mol?MJ^?1 and 0.48 (June). The mean daily values of fFEC and fE at the solar noon were 2.16?mol?MJ^?1 and 0.47 across the entire growing season, respectively. Both PPFD and the ever first reported LPR showed strong diurnal patterns. However, they had opposite trends. PPFD was high around noon, resulting in low values of LPR during the same time period. Both were found to be highly correlated with global solar radiation R _S, solar elevation angle h, and the clearness index K _t. Using the best subset selection of variables, two parametric models were developed for estimating PPFD and LPR, which can easily be applied in radiometric sites, by recording only global solar radiation measurements. These two models were found to be involved with the most commonly measured global solar radiation (R _S) and two large-scale geometric parameters, i.e., extraterrestrial solar radiation and solar elevation. The models were therefore insensitive to local weather conditions such as temperature. In particular, with two test data sets collected in USA and Greece, it was verified that the models could be extended across different geographical areas, where they performed well. Therefore, these two hourly based models can be used to provide precise PAR-related values, such as those required for developing precise vegetation growth models.     

Comparison of precipitation in the regional climate model BALTIMOS to radar observations

Observational data and simulations of the regional climate system Baltic integrated model system (BALTIMOS) were used to study precipitation in the Baltic Sea and its drainage basin with a special focus on the diurnal cycle. The study includes a general evaluation of BALTIMOS precipitation, showing that BALTIMOS has too many light rain events causing an overestimation of the total annual precipitation amount. The diurnal cycle as well as its spatial distribution was analysed. BALTIMOS captures the broad characteristics: a significant diurnal variability with an afternoon peak above land and weak variability with a nocturnal peak above sea. An algorithm to distinguish between frontal and convective precipitation was applied to examine the diurnal cycle more thoroughly. The local solar time of maximum rain in summer is about 1 to 2 h earlier in BALTIMOS than in radar observations of precipitation.     

An Experimental Comparison of Ordinary and Universal Kriging and Inverse Distance Weighting

A factorial, computational experiment was conducted to compare the spatial interpolation accuracy of ordinary and universal kriging and two types of inverse squared-distance weighting. The experiment considered, in addition to these four interpolation methods, the effects of four data and sampling characteristics: surface type, sampling pattern, noise level, and strength of small-scale spatial correlation. Interpolation accuracy was measured by the natural logarithm of the mean squared interpolation error. Main effects of all five factors, all two-factor interactions, and several three-factor interactions were highly statistically significant. Among numerous findings, the most striking was that the two kriging methods were substantially superior to the inverse distance weighting methods over all levels of surface type, sampling pattern, noise, and correlation.     

Precipitation and temperature space–time variability and extremes in the Mediterranean region: evaluation of dynamical and statistical downscaling methods

This study evaluates how statistical and dynamical downscaling models as well as combined approach perform in retrieving the space–time variability of near-surface temperature and rainfall, as well as their extremes, over the whole Mediterranean region. The dynamical downscaling model used in this study is the Weather Research and Forecasting (WRF) model with varying land-surface models and resolutions (20 and 50 km) and the statistical tool is the Cumulative Distribution Function-transform (CDF-t). To achieve a spatially resolved downscaling over the Mediterranean basin, the European Climate Assessment and Dataset (ECA&D) gridded dataset is used for calibration and evaluation of the downscaling models. In the frame of HyMeX and MED-CORDEX international programs, the downscaling is performed on ERA-I reanalysis over the 1989–2008 period. The results show that despite local calibration, CDF-t produces more accurate spatial variability of near-surface temperature and rainfall with respect to ECA&D than WRF which solves the three-dimensional equation of conservation. This first suggests that at 20–50 km resolutions, these three-dimensional processes only weakly contribute to the local value of temperature and precipitation with respect to local one-dimensional processes. Calibration of CDF-t at each individual grid point is thus sufficient to reproduce accurately the spatial pattern. A second explanation is the use of gridded data such as ECA&D which smoothes in part the horizontal variability after data interpolation and damps the added value of dynamical downscaling. This explains partly the absence of added-value of the 2-stage downscaling approach which combines statistical and dynamical downscaling models. The temporal variability of statistically downscaled temperature and rainfall is finally strongly driven by the temporal variability of its forcing (here ERA-Interim or WRF simulations). CDF-t is thus efficient as a bias correction tool but does not show any added-value regarding the time variability of the downscaled field. Finally, the quality of the reference observation dataset is a key issue. Comparison of CDF-t calibrated with ECA&D dataset and WRF simulations to local measurements from weather stations not assimilated in ECA&D, shows that the temporal variability of the downscaled data with respect to the local observations is closer to the local measurements than to ECA&D data. This highlights the strong added-value of dynamical downscaling which improves the temporal variability of the atmospheric dynamics with regard to the driving model. This article highlights the benefits and inconveniences emerging from the use of both downscaling techniques for climate research. Our goal is to contribute to the discussion on the use of downscaling tools to assess the impact of climate change on regional scales.     

Surface modes in the infrared spectrum of hydrous minerals: the OH stretching modes of bayerite

The theoretical infrared (IR) and Raman spectra of bayerite (β-Al(OH)₃) are computed in the density functional theory framework, using the linear response theory. The results are consistent with the occurrence of six non-equivalent OH groups in a bayerite structure with space group P2₁/n. Similar to gibbsite, the transmission powder IR spectrum of bayerite in the region of the OH stretching bands is found to depend on the shape of particles. In particular, the broadening of the strong band observed at about 3,460 cm⁻¹ in the spectrum of Al hydroxides is related to the electrostatic charges occurring at the surface of the polarized dielectric particles. The experimental correlation observed between the shape of this band and morphological parameters has therefore a physical, instead of chemical, origin.