Idealized Mesoscale Model Simulations of Open Cellular Convection Over the Sea

The atmospheric conditions during an observed case of open cellular convection over the North Sea were simulated using the Weather Research and Forecasting (WRF) numerical model. Wind, temperature and water vapour mixing ratio profiles from the WRF simulation were used to initialize an idealized version of the model, which excluded the effects of topography, surface inhomogeneities and large-scale weather forcing. Cells with an average diameter of 17.4 km developed. Simulations both with and without a capping inversion were made, and the cell-scale kinetic energy budget was calculated for each case. By considering all sources of explicit diffusion in the model, the budgets were balanced. In comparison with previous work based on observational studies, the use of three-dimensional, gridded model data afforded the possibility of calculating all terms in the budgets, which showed that the important terms in the budgets were buoyancy, pressure balance and inter-scale transfer to subgrid scales. Cells were also composited to calculate the average cell-scale flow and each of the budget terms on two-dimensional cross-sections through the cells, parallel and perpendicular to the mean wind direction.     

Estimating spatial and temporal patterns of urban anthropogenic heat fluxes for UK cities: the case of Manchester

A model is proposed for determining the temporal and spatial patterns of anthropogenic heat fluxes in UK urban areas. It considers buildings, traffic, and metabolic heat flux sources and has been evaluated to a good accuracy against alternative data for the Greater Manchester area in the UK. Results are presented at spatial resolution of 200 ×?200 m although the model itself is scalable depending on data availability. In this paper, results are generated using a set of urban morphology units so that detailed and time-consuming accounting of individual building and road emissions is not required. The model estimates a mean heat emission of 6.12 Wm^-2 across Greater Manchester, with values in the region of 10 Wm^-2 for non central urbanized areas and 23 Wm^-2 in city center areas. Despite this difference, the results are not described by a simple distance decay function, as has been reported for other cities, due to the influence of satellite towns and the influence of the road network. Buildings are the dominant emitter, contributing some 60% of total emissions across the city compared to around 32% for road traffic and 8% for metabolic sources.     

Water age in stormwater management ponds and stormwater management pond-treated catchments

Expansion of impervious surface cover results in “flashy” hydrologic response, elevated flood risk, and degraded water quality in urban watersheds. Stormwater management ponds (SWMPs) are often engineered into stream networks to mitigate these issues. A clearer understanding of how water is stored and released from SWMPs and SWMP-treated catchments is required to better represent these engineered systems in hydrological and water quality models of urban and urbanizing watersheds. Stable water isotopes were used to compare water age in SWMPs and SWMP-treated catchments in an urbanizing watershed. We sampled water biweekly from two SWMPs and five stream sites with varying land cover and stormwater control in their catchments. Two inverse transit time proxies (damping ratio and young water fraction) were computed along with the mean transit time (MTT) by sine–wave fitting for each SWMP and stream site using the δ¹⁸O and δ²H data. Water entering the SWMPs was consistently older (224 and 177 days) than water in or exiting the ponds (ranging from 46 to 91 days and 39 to 67 days, respectively). This finding is likely due to a combination of groundwater infiltration into broken sewer pipes that transport water into the ponds and a bias toward baseflow sampling. At the catchment scale, detention provided by SWMPs was not found to be more significant than the interactive effects of impervious cover, surficial geology, land use proportions, and catchment size in determining MTT. Overall, surficial geology explained the most variation in MTT among the seven sites. This study illustrates the potential for isotope-based approaches of water age to provide information on individual SWMP functioning and the influence of SWMPs on catchment-scale water movement.     

Quantifying chloride retention and release in urban stormwater management ponds using a mass balance approach

Chloride (Cl⁻) in urban waterways largely originates from runoff containing deicing salts. Cl⁻ is retained in watersheds after deicing ends, resulting in deleterious effects on aquatic biota. Stormwater management ponds (SWMPs), designed to mitigate ‘flashy’ urban runoff response, are known to impact pollutant transport. However, there is little information on what role SWMPs play in the timing and magnitude of Cl⁻ transport over different timescales. This study quantifies the mass of Cl⁻ retained in two SWMPs over varying timescales. Both ponds are in an urbanizing watershed in south-central Ontario; one drains a commercial area, the other, a residential area. High frequency measurements of water level and specific conductivity, from which flow and Cl⁻ concentration were derived, were taken with sensors at pond inlets and outlets. For one SWMP, data were also collected upstream and downstream of the confluence of the pond outflow and the receiving creek to quantify the in-stream response to Cl⁻-laden pond outflows. The findings suggest that SWMPs likely play a role in watershed-scale Cl⁻ retention; one SWMP consistently retained Cl⁻ while the other had variable retention and release of Cl⁻. In the receiving creek, Cl⁻ concentrations downstream of the pond exceeded the acute toxicity threshold for aquatic organisms twice as often as concentrations upstream of the pond, and Cl⁻ pulses corresponded to Cl⁻ release events from the pond. The results of this study suggest that SWMPs concentrate spatially distributed salt inputs and modify the timing and magnitude of their release to receiving streams. Stream reaches that receive water inputs from SWMPs may be more vulnerable to Cl⁻ toxicity than reaches that do not receive flow via SWMPs. The results of this study will help parameterize the role of SWMPs in watershed-scale Cl⁻ transport models and geospatial models of salt vulnerable areas.     

Determination of bulk density for small meteorite fragments via visible light 3‐D laser imaging

Abstract– Bulk density is an important intrinsic property of meteorites, but the necessary bulk volume measurement is difficult to do in a truly nondestructive way. Archimedean methods involving the displacement of a 40–100 μm beads “fluid” are commonly applied, but can encounter systematic errors. Herein, we report a visible light laser imaging technique for the nondestructive measurement of meteorite surface features, allowing for the subsequent assembly of 3‐D volumetric models; the method is particularly applicable to small meteorite fragments and to fragile specimens. We have acquired laser image data for 24 fragments from 18 ordinary chondrites, carbonaceous chondrites, and achondrites, with masses ranging from 265.0 to 1.2 g. Laser imaging bulk density is consistent between sister fragments of meteorites down to sizes of about 0.5 cm³, an order of magnitude smaller than can be reliably measured with Archimedean beads techniques. Uncertainty is less than 2% for fragments >4 cm³, and typically between 2 and 4% for small fragments <4 cm³. For 10 fragments, 3‐D laser imaging volumes are on average 1.3% smaller than those obtained with Archimedean beads. In a wider comparison using 21 meteorite fragments, 3‐D laser imaging bulk densities are on average 2.14 ± 2.36% greater than the corresponding Archimedean method literature values for these meteorites. Difficulties in the procedure of 3‐D image alignment may lead to a slight overestimation of meteorite bulk density, and so laser imaging‐based bulk densities are maximum estimates that can be viewed as being complementary to the minimum bulk density estimates obtained using Archimedean beads methods.     

Atmospheric modeling of Mars methane surface releases

We utilize the MarsWRF general circulation model (GCM) to address the behavior of gas plumes in the Martian atmosphere, with the specific goal of characterizing the source of the recently identified methane detection in the Martian atmosphere. These observations have been interpreted as the release of methane from localized surface sources with spatial and temporal variabilities. Due to the limited temporal coverage of ground-based observations, we use a GCM to simulate the development of passive atmospheric plumes over relevant timescales. The observations can be reproduced best if the release occurred just before the time of observation—no more than 1–2 sols earlier—and if this release were nearly instantaneous rather than a slow, steady emission. Furthermore, it requires a source region spanning a broad latitudinal range rather than a point emission. While the accuracy of our conclusions about this specific methane release scenario is limited by the uncertainties inherent in GCM simulations of the Martian atmosphere, our findings regarding generalized plume behavior are robust, and illustrate the potential power of numerical modeling for constraining plume source conditions.     

Stratospheric superrotation in the TitanWRF model

TitanWRF general circulation model simulations performed without sub-grid-scale horizontal diffusion of momentum produce roughly the observed amount of superrotation in Titan’s stratosphere. We compare these results to Cassini-Huygens measurements of Titan’s winds and temperatures, and predict temperature and winds at future seasons. We use angular momentum and transformed Eulerian mean diagnostics to show that equatorial superrotation is generated during episodic angular momentum ‘transfer events’ during model spin-up, and maintained by similar (yet shorter) events once the model has reached steady state. We then use wave and barotropic instability analysis to suggest that these transfer events are produced by barotropic waves, generated at low latitudes then propagating poleward through a critical layer, thus accelerating low latitudes while decelerating the mid-to-high latitude jet in the late fall through early spring hemisphere. Finally, we identify the dominant waves responsible for the transfers of angular momentum close to northern winter solstice during spin-up and at steady state. Problems with our simulations include peak latitudinal temperature gradients and zonal winds occurring ∼60 km lower than observed by Cassini CIRS, and no reduction in zonal wind speed around 80 km, as was observed by Huygens. While the latter may have been due to transient effects (e.g. gravity waves), the former suggests that our low (∼420 km) model top is adversely affecting the circulation near the jet peak, and/or that we require active haze transport in order to correctly model heating rates and thus the circulation. Future work will include running the model with a higher top, and including advection of a haze particle size distribution.     

Synthesis of a spinifex‐textured basalt as an analog to Gusev crater basalts,Mars

Abstract– Analyses by the Mars Exploration Rover (MER), Spirit, of Martian basalts from Gusev crater show that they are chemically very different from terrestrial basalts, being characterized in particular by high Mg‐ and Fe‐contents. To provide suitable analog basalts for the International Space Analogue Rockstore (ISAR), a collection of analog rocks and minerals for preparing in situ space missions, especially, the upcoming Mars mission MSL‐2011 and the future international Mars‐2018 mission, it is necessary to synthesize Martian basalts. The aim of this study was therefore to synthesize Martian basalt analogs to the Gusev crater basalts, based on the geochemical data from the MER rover Spirit. We present the results of two experiments, one producing a quench‐cooled basalt (<1 h) and one producing a more slowly cooled basalt (1 day). Pyroxene and olivine textures produced in the more slowly cooled basalt were surprisingly similar to spinifex textures in komatiites, a volcanic rock type very common on the early Earth. These kinds of ultramafic rocks and their associated alteration products may have important astrobiological implications when associated with aqueous environments. Such rocks could provide habitats for chemolithotrophic microorganisms, while the glass and phyllosilicate derivatives can fix organic compounds.     

Evaporation and recondensation of sodium in Semarkona Type II chondrules

We have investigated the Na distributions in Semarkona Type II chondrules by electron microprobe, analyzing olivine and melt inclusions in it, mesostasis and bulk chondrule, to see whether they indicate interactions with an ambient gas during chondrule formation. Sodium concentrations of bulk chondrule liquids, melt inclusions and mesostases can be explained to a first approximation by fractional crystallization of olivine ± pyroxene. The most primitive olivine cores in each chondrule are mostly between Fa₈ and Fa₁₃, with 0.0022–0.0069 ± 0.0013 wt.% Na₂O. Type IIA chondrule olivines have consistently higher Na contents than olivines in Type IIAB chondrules. We used the dependence of olivine–liquid Na partitioning on FeO in olivine as a measure of equilibration. Extreme olivine rim compositions are ～Fa₃₅ and 0.03 wt.% Na₂O and are close to being in equilibrium with the mesostasis glass. Olivine cores compared with the bulk chondrule compositions, particularly in IIA chondrules, show very high apparent D_Na, indicating disequilibrium and suggesting that chondrule initial melts were more Na-rich than present chondrule bulk compositions. The apparent D_Na values correlate with the Na concentrations of the olivine, but not with concentrations in the bulk melt. We use equilibrium D_Na to find the Na content of the true parent liquid and estimate that Type IIA chondrules lost more than half their Na and recondensation was incomplete, whereas Type IIAB chondrules recovered most of theirs in their mesostases.Glass inclusions in olivine have lower Na than expected from fractionation of bulk composition liquids, and mesostases have higher Na than expected in calculated daughter liquids formed by fractional crystallization alone. These observations also require open system behavior of chondrules, specifically evaporation of Na before formation of melt inclusions followed by recondensation of Na in mesostases. Within this record of evaporation followed by recondensation, there is no indication of a stage with zero Na in the chondrules, which is predicted by models for shock wave cooling at canonical nebular pressures, suggesting high P_T.The high Na concentrations in olivine and mesostases indicate very high P_Na while chondrules were molten. This may be explained by local, very high particle densities where Type II chondrules formed. The high P_T, P_Na and number densities of chondrules implied suggest formation in debris clouds after protoplanetary collisions as an alternative to formation after passage of shock waves through large particle-rich clumps in the disk. Encounters of partially molten chondrules should have been frequent in these dense swarms. However, in many ordinary chondrites like Semarkona, “cluster chondrites”, compound chondrules are not abundant but instead chondrules aggregated into clusters. Chondrule melting, cooling and clustering in dense swarms contributed to rapid accretion, possibly after collision, by fallback on the grandparent body and by reaccretion as a new body downrange.