Testing variograms for positive-definiteness

A method based on Bochner's theorem is described for demonstrating the positive-definiteness of variogram models and for generating classes of valid variogram functions.Work performed while on leave at the Centre de Géostatistique et de Morphologie Mathématique.     

Measuring Urban Sustainability: the potential and pitfalls of city rankings

Measuring the sustainability of a city, and ranking the city against other cities, is increasingly seen as a way to achieve urban sustainability. A number of studies have highlighted how Sydney is unsustainable, with two recent ranking exercises appearing to confirm these findings by placing Sydney below many other Australian cities for urban sustainability. These two exercises are situated within international moves towards ranking cities. The article explores the assumptions and criteria used by Mercer in its Eco-City Index and the Australian Conservation Foundation in its Sustainable Cities Index. The conclusion discusses how Sydney could become more sustainable, and identifies why normative urban sustainability rankings should be recognised for their potential while simultaneously being treated with caution.     

Observable consequences of planet formation models in systems with close-in terrestrial planets

To date, two planetary systems have been discovered with close-in, terrestrial-mass planets     . Many more such discoveries are anticipated in the coming years with radial velocity and transit searches. Here we investigate the different mechanisms that could form 'hot Earths' and their observable predictions. Models include: (1) in situ accretion; (2) formation at larger orbital distance followed by inward 'type 1' migration; (3) formation from material being 'shepherded' inward by a migrating gas giant planet; (4) formation from material being shepherded by moving secular resonances during dispersal of the protoplanetary disc; (5) tidal circularization of eccentric terrestrial planets with close-in perihelion distances and (6) photoevaporative mass-loss of a close-in giant planet. Models 1–4 have been validated in previous work. We show that tidal circularization can form hot Earths, but only for relatively massive planets     with very close-in perihelion distances (≲0.025 au), and even then the net inward movement in orbital distance is at most only 0.1–0.15 au. For planets of less than     , photoevaporation can remove the planet's envelope and leave behind the solid core on a Gyr time-scale, but only for planets inside 0.025–0.05 au. Using two quantities that are observable by current and upcoming missions, we show that these models each produce unique signatures, and can be observationally distinguished. These observables are the planetary system architecture (detectable with radial velocities, transits and transit timing) and the bulk composition of transiting close-in terrestrial planets (measured by transits via the planet's radius).     

Nonlinear Force-Free Modeling of Coronal Magnetic Fields. II. Modeling a Filament Arcade and Simulated Chromospheric and Photospheric Vector Fields

We compare a variety of nonlinear force-free field (NLFFF) extrapolation algorithms, including optimization, magneto-frictional, and Grad – Rubin-like codes, applied to a solar-like reference model. The model used to test the algorithms includes realistic photospheric Lorentz forces and a complex field including a weakly twisted, right helical flux bundle. The codes were applied to both forced “photospheric” and more force-free “chromospheric” vector magnetic field boundary data derived from the model. When applied to the chromospheric boundary data, the codes are able to recover the presence of the flux bundle and the field’s free energy, though some details of the field connectivity are lost. When the codes are applied to the forced photospheric boundary data, the reference model field is not well recovered, indicating that the combination of Lorentz forces and small spatial scale structure at the photosphere severely impact the extrapolation of the field. Preprocessing of the forced photospheric boundary does improve the extrapolations considerably for the layers above the chromosphere, but the extrapolations are sensitive to the details of the numerical codes and neither the field connectivity nor the free magnetic energy in the full volume are well recovered. The magnetic virial theorem gives a rapid measure of the total magnetic energy without extrapolation though, like the NLFFF codes, it is sensitive to the Lorentz forces in the coronal volume. Both the magnetic virial theorem and the Wiegelmann extrapolation, when applied to the preprocessed photospheric boundary, give a magnetic energy which is nearly equivalent to the value derived from the chromospheric boundary, but both underestimate the free energy above the photosphere by at least a factor of two. We discuss the interpretation of the preprocessed field in this context. When applying the NLFFF codes to solar data, the problems associated with Lorentz forces present in the low solar atmosphere must be recognized: the various codes will not necessarily converge to the correct, or even the same, solution. On 07/07/2007, the NLFFF team was saddened by the news that Tom Metcalf had died as the result of an accident. We remain grateful for having had the opportunity to benefit from his unwavering dedication to the problems encountered in attempting to understand the Sun’s magnetic field; Tom had completed this paper several months before his death, leading the team through the many steps described above.     

Resolving the Azimuthal Ambiguity in Vector Magnetogram Data with the Divergence-Free Condition: Theoretical Examination

We demonstrate that the azimuthal ambiguity that is present in solar vector magnetogram data can be resolved with line-of-sight and horizontal heliographic derivative information by using the divergence-free property of magnetic fields without additional assumptions. We discuss the specific derivative information that is sufficient to resolve the ambiguity away from disk centre, with particular emphasis on the line-of-sight derivative of the various components of the magnetic field. Conversely, we also show cases where ambiguity resolution fails because sufficient line-of-sight derivative information is not available. For example, knowledge of only the line-of-sight derivative of the line-of-sight component of the field is not sufficient to resolve the ambiguity away from disk centre.     

XUV Photometer System (XPS): Improved Solar Irradiance Algorithm Using CHIANTI Spectral Models

Solar soft X-ray (XUV) radiation is highly variable on all time scales and strongly affects Earth’s ionosphere and upper atmosphere; consequently, the solar XUV irradiance is important for atmospheric studies and for space weather applications. Although there have been several recent measurements of the solar XUV irradiance, detailed understanding of the solar XUV irradiance, especially its variability during flares, has been hampered by the broad bands measured in the XUV range. In particular, the simple conversion of the XUV photometer signal into irradiance, in which a static solar spectrum is assumed, overestimates the flare variations by more than a factor of two as compared to the atmospheric response to the flares. To address this deficiency in the simple conversion, an improved algorithm using CHIANTI spectral models has been developed to process the XUV Photometer System (XPS) measurements with its broadband photometers. Model spectra representative of quiet Sun, active region, and flares are combined to match the signals from the XPS and produce spectra from 0.1 to 40 nm in 0.1-nm intervals for the XPS Level 4 data product. The two XPS instruments are aboard NASA’s Solar Radiation and Climate Experiment (SORCE) and Thermosphere, Ionosphere, Mesosphere, Energetics, and Dynamics (TIMED) satellites. In addition, the XPS responsivities have been updated for the latest XPS data processing version. The new XPS results are consistent with daily variations from the previous simple conversion technique used for XPS and are also consistent with spectral measurements made at wavelengths longer than 27 nm. Most importantly, the XPS flare variations are reduced by factors of 2 – 4 at wavelengths shorter than 14 nm and are more consistent, for the first time, with atmospheric response to solar flares. Along with the details of the new XPS algorithm, several comparisons to dayglow and photoelectron measurements and model results are also presented to help verify the accuracy of the new XUV irradiance spectra.     

Spectroscopy, morphometry, and photoclinometry of Titan's dunefields from Cassini/VIMS

Fine-resolution (500 m/pixel) Cassini Visual and Infrared Mapping Spectrometer (VIMS) T20 observations of Titan resolve that moon's sand dunes. The spectral variability in some dune regions shows that there are sand-free interdune areas, wherein VIMS spectra reveal the exposed dune substrate. The interdunes from T20 are, variously, materials that correspond to the equatorial bright, 5-μm-bright, and dark blue spectral units. Our observations show that an enigmatic “dark red” spectral unit seen in T5 in fact represents a macroscopic mixture with 5-μm-bright material and dunes as its spectral endmembers. Looking more broadly, similar mixtures of varying amounts of dune and interdune units of varying composition can explain the spectral and albedo variability within the dark brown dune global spectral unit that is associated with dunes. The presence of interdunes indicates that Titan's dunefields are both mature and recently active. The spectrum of the dune endmember reveals the sand to be composed of less water ice than the rest of Titan; various organics are consistent with the dunes' measured reflectivity. We measure a mean dune spacing of 2.1 km, and find that the dunes are oriented on the average in an east-west direction, but angling up to 10° from parallel to the equator in specific cases. Where no interdunes are present, we determine the height of one set of dunes photoclinometrically to be between 30 and 70 m. These results pave the way for future exploration and interpretation of Titan's sand dunes.     

Short communication: public interest in rip currents relative to other natural hazards: evidence from Google search data

Natural Hazards - Informing and warning beach users about the danger posed by the rip current hazard remain a difficult challenge. Recent evidence suggests that warning signs and flags alone have...     

Spectral analysis of continuous redox data reveals geochemical dynamics near the stream–aquifer interface

Changes in streamflow and water table elevation influence oxidation–reduction (redox) conditions near river–aquifer interfaces, with potentially important consequences for solute fluxes and biogeochemical reaction rates. Although continuous measurements of groundwater chemistry can be arduous, in situ sensors reveal chemistry dynamics across a wide range of timescales. We monitored redox potential in an aquifer adjacent to a tidal river and used spectral and wavelet analyses to link redox responses to hydrologic perturbations within the bed and banks. Storms perturb redox potential within both the bed and banks over timescales of days to weeks. Tides drive semidiurnal oscillations in redox potential within the streambed that are absent in the banks. Wavelet analysis shows that tidal redox oscillations in the bed are greatest during late summer (wavelet magnitude of 5.62 mV) when river stage fluctuations are on the order of 70 cm and microbial activity is relatively high. Tidal redox oscillations diminish during the winter (wavelet magnitude of 2.73 mV) when river stage fluctuations are smaller (on the order of 50 cm) and microbial activity is presumably low. Although traditional geochemical observations are often limited to summer baseflow conditions, in situ redox sensing provides continuous, high‐resolution chemical characterization of the subsurface, revealing transport and reaction processes across spatial and temporal scales in aquifers.