Model-based remote sensing algorithms for particulate organic carbon (POC) in the Northeastern Gulf of Mexico

Hydrographic data, including particulate organic carbon (POC) from the Northeastern Gulf of Mexico (NEGOM) study, were combined with remotely-sensed SeaWiFS data to estimate POC concentration using principal component analysis (PCA). The spectral radiance was extracted at each NEGOM station, digitized, and averaged. The mean value and spurious trends were removed from each spectrum. De-trended data included six wavelengths at 58 stations. The correlation between the weighting factors of the first six eigenvectors and POC concentration were applied using multiple linear regression. PCA algorithms based on the first three, four, and five modes accounted for 90, 95, and 98% of total variance and yielded significant correlations with POC with R ² = 0.89, 0.92, and 0.93. These full waveband approaches provided robust estimates of POC in various water types. Three different analyses (root mean square error, mean ratio and standard deviation) showed similar error estimates, and suggest that spectral variations in the modes defined by just the first four characteristic vectors are closely correlated with POC concentration, resulting in only negligible loss of spectral information from additional modes. The use of POC algorithms greatly increases the spatial and temporal resolution for interpreting POC cycling and can be extrapolated throughout and perhaps beyond the area of shipboard sampling.     

The Ashima/MIT Mars GCM and argon in the martian atmosphere

We investigate the ability of modern general circulation models (GCMs) to simulate transport in the martian atmosphere using measurements of argon as a proxy for the transport processes. Argon provides the simplest measure of transport as it is a noble gas with no sinks or sources on seasonal timescales. Variations in argon result solely from ‘freeze distillation’, as the atmosphere condenses at the winter poles, and from atmospheric transport. Comparison of all previously published models when rescaled to a common definition of the argon enhancement factor (EF) suggest that models generally do a poor job in predicting the peak enhancement in southern winter over the winter pole – the time when the capability of the model transport approaches are most severely tested. Despite observed peak EF values of ～6, previously published model predictions peaked at EF values of only 2–3. We introduce a new GCM that provides a better treatment of mass conservation within the dynamical core, includes more sophisticated tracer transport approaches, and utilizes a cube–sphere grid structure thus avoiding the grid-point convergence problem at the pole that exists for most current Mars GCMs. We describe this model – the Ashima Research/Massachusetts Institute of Technology Mars General Circulation Model (Ashima/MIT Mars GCM) and use it to demonstrate the significant sensitivity of peak EF to the choices of transport approach for both tracers and heat. We obtain a peak EF of 4.75 which, while over 50% higher than any prior model, remains well short of the observed value. We show that the polar EF value in winter is primarily determined by the competition between two processes: (1) mean meridional import of lower-latitude air not enriched in argon and (2) the leakage of enriched argon out of the polar column by eddies in the lowest atmospheric levels. We suggest possibilities for improving GCM representation of the CO₂ cycle and the general circulation that may further improve the simulation of the argon cycle. We conclude that current GCMs may be insufficient for detailed simulation of transport-sensitive problems like the water cycle and potentially also the dust cycle.     

N-body simulations of cohesion in dense planetary rings: A study of cohesion parameters

We present results from a large suite of simulations of Saturn’s dense A and B rings using a new model of particle sticking in local simulations (Perrine, R.P., Richardson, D.C., Scheeres, D.J. [2011]. Icarus 212, 719–735). In this model, colliding particles can be incorporated into or help fragment rigid aggregations on the basis of certain user-specified parameters that can represent van der Waals forces or interlocking surface frost layers.Our investigation is motivated by laboratory results that show that interpenetration of surface layers can allow impacting frost-covered ice spheres to stick together. In these experiments, cohesion only occurs below specific impact speeds, which happen to be characteristic of impact speeds in Saturn’s rings. Our goal is to determine if weak bonding is consistent with ring observations, to constrain cohesion parameters in light of existing ring observations, to make predictions about particle populations throughout the rings, and to discover other diagnostics that may constrain bonding parameters.We considered the effects of five parameters on the equilibrium characteristics of our ring simulations: speed-based merge and fragmentation limits, bond strength, ring surface density, and patch orbital distance (i.e., the A or B ring), some with both monodisperse and polydisperse comparison cases. In total, we present data from 95 simulations.We find that weak cohesion is consistent with observations of the A and B rings (e.g., French, R.G., Nicholson, P.D. [2000]. Icarus 145, 502–523), and we present a range of simulation parameters that reproduce the observed size distribution and maximum particle size. It turns out that the parameters that match observations differ between the A and B rings, and we discuss the potential implications of this result. We also comment on other observable consequences of cohesion for the rings, such as optical depth and scale height effects, and discuss whether very large objects (e.g., “propeller” source objects) are grown bottom-up from cohesion of smaller ring particles.     

Numerical simulations of granular dynamics II: Particle dynamics in a shaken granular material

Surfaces of planets and small bodies of our Solar System are often covered by a layer of granular material that can range from a fine regolith to a gravel-like structure of varying depths. Therefore, the dynamics of granular materials are involved in many events occurring during planetary and small-body evolution thus contributing to their geological properties.We demonstrate that the new adaptation of the parallel N-body hard-sphere code pkdgrav has the capability to model accurately the key features of the collective motion of bidisperse granular materials in a dense regime as a result of shaking. As a stringent test of the numerical code we investigate the complex collective ordering and motion of granular material by direct comparison with laboratory experiments. We demonstrate that, as experimentally observed, the scale of the collective motion increases with increasing small-particle additive concentration.We then extend our investigations to assess how self-gravity and external gravity affect collective motion. In our reduced-gravity simulations both the gravitational conditions and the frequency of the vibrations roughly match the conditions on asteroids subjected to seismic shaking, though real regolith is likely to be much more heterogeneous and less ordered than in our idealised simulations. We also show that collective motion can occur in a granular material under a wide range of inter-particle gravity conditions and in the absence of an external gravitational field. These investigations demonstrate the great interest of being able to simulate conditions that are to relevant planetary science yet unreachable by Earth-based laboratory experiments.     

Litterfall,litter and associated chemistry in a dry sclerophyll eucalypt forest and a pine plantation in south-eastern Australia: 1. Litterfall and litter

Litterfall was measured in a dry schlerophyll eucalypt forest and a nearby Pinus radiata plantation of similar tree density and basal area near Canberra in south-eastern Australia. Total annual litterfall for the eucalypts was 329 g m⁻², compared with 180 g m⁻² for the pines, with the bark component being 52 g m⁻² for eucalypts and zero for pines. Barkfall did not occur for the eucalypts during the drought of 1982–1983 but complete bark shedding occurred during the subsequent very wet year when barkfall was 177 g m⁻² for Eucalyptus rossii and 146 g m⁻² for Eucalyptus mannifera (9·3 and 7·6 g m⁻² of basal area, respectively). Barkfall of E. rossii responded to rainfall in the period autumn to early summer, whereas E. mannifera responded to summer rainfall. In the eucalypt forest floor-litter was stratified into a surface layer where the components were substantially intact, and a cohesive layer where the components were fragmented and bound together by fungal hyphae. The amount and residence times of loose and cohesive floor-litter were 1056 g m⁻² and 3·2 years, respectively, for the loose litter layer; and 1164 g m⁻² and 3·5 years for the cohesive layer. The litter biomass represented 17% of the estimated total above-ground biomass of 127 tonnes ha⁻¹. A previous study showed roots to be 25% of total biomass, suggesting a total biomass of 167 tonnes ha⁻¹. © 1998 John Wiley & Sons, Ltd.