Solubility and Diffusivity of Water in Lunar Basalt, and Chemical Zonation in Olivine-Hosted Melt Inclusions

I report the solubility and diffusivity of water in lunar basalt and an iron-free basaltic analogue at 1 atm and 1350 °C. Such parameters are critical for understanding the degassing histories of lunar pyroclastic glasses. Solubility experiments have been conducted over a range of fO₂ conditions from three log units below to five log units above the iron-wüstite buffer (IW) and over a range of pH₂/pH₂O from 0.03 to 24. Quenched experimental glasses were analyzed by Fourier transform infrared spectroscopy (FTIR) and secondary ionization mass spectrometry (SIMS) and were found to contain up to ~420 ppm water. Results demonstrate that, under the conditions of our experiments: (1) hydroxyl is the only H-bearing species detected by FTIR; (2) the solubility of water is proportional to the square root of pH₂O in the furnace atmosphere and is independent of fO₂ and pH₂/pH₂O; (3) the solubility of water is very similar in both melt compositions; (4) the concentration of H₂ in our iron-free experiments is <3 ppm, even at oxygen fugacities as low as IW-2.3 and pH₂/pH₂O as high as 24; and (5) SIMS analyses of water in iron-rich glasses equilibrated under variable fO₂ conditions can be strongly influenced by matrix effects, even when the concentrations of water in the glasses are low. Our results can be used to constrain the entrapment pressure of the lunar melt inclusions of Hauri et al. (2011).

Diffusion experiments were conducted over a range of fO₂ conditions from IW-2.2 to IW+6.7 and over a range of pH₂/pH₂O from nominally zero to ~10. The water concentrations measured in our quenched experimental glasses by SIMS and FTIR vary from a few ppm to ~430 ppm. Water concentration gradients are well described by models in which the diffusivity of water (D^*_water) is assumed to be constant. The relationship between D^*_water and water concentration is well described by a modified speciation model (Ni et al. 2012) in which both molecular water and hydroxyl are allowed to diffuse. The success of this modified speciation model for describing our results suggests that we have resolved the diffusivity of hydroxyl in basaltic melt for the first time. Best-fit values of D^*_water for our experiments on lunar basalt vary within a factor of ~2 over a range of pH₂/pH₂O from 0.007 to 9.7, a range of fO₂ from IW-2.2 to IW+4.9, and a water concentration range from ~80 ppm to ~280 ppm. The relative insensitivity of our best-fit values of D^*_water to variations in pH₂ suggests that H₂ diffusion was not significant during degassing of the lunar glasses of Saal et al. (2008). D^*_water during dehydration and hydration in H₂/CO₂ gas mixtures are approximately the same, which supports an equilibrium boundary condition for these experiments. However, dehydration experiments into CO₂ and CO/CO₂ gas mixtures leave some scope for the importance of kinetics during dehydration into H-free environments. The value of D^*_water chosen by Saal et al. (2008) for modeling the diffusive degassing of the lunar volcanic glasses is within a factor of three of our measured value in our lunar basaltic melt at 1350 °C.

In Chapter 4 of this thesis, I document significant zonation in major, minor, trace, and volatile elements in naturally glassy olivine-hosted melt inclusions from the Siqueiros Fracture Zone and the Galapagos Islands. Components with a higher concentration in the host olivine than in the melt (MgO, FeO, Cr₂O₃, and MnO) are depleted at the edges of the zoned melt inclusions relative to their centers, whereas except for CaO, H₂O, and F, components with a lower concentration in the host olivine than in the melt (Al₂O₃, SiO₂, Na₂O, K₂O, TiO₂, S, and Cl) are enriched near the melt inclusion edges. This zonation is due to formation of an olivine-depleted boundary layer in the adjacent melt in response to cooling and crystallization of olivine on the walls of the melt inclusions concurrent with diffusive propagation of the boundary layer toward the inclusion center.

Concentration profiles of some components in the melt inclusions exhibit multicomponent diffusion effects such as uphill diffusion (CaO, FeO) or slowing of the diffusion of typically rapidly diffusing components (Na₂O, K₂O) by coupling to slow diffusing components such as SiO₂ and Al₂O₃. Concentrations of H2O and F decrease towards the edges of some of the Siqueiros melt inclusions, suggesting either that these components have been lost from the inclusions into the host olivine late in their cooling histories and/or that these components are exhibiting multicomponent diffusion effects.

A model has been developed of the time-dependent evolution of MgO concentration profiles in melt inclusions due to simultaneous depletion of MgO at the inclusion walls due to olivine growth and diffusion of MgO in the melt inclusions in response to this depletion. Observed concentration profiles were fit to this model to constrain their thermal histories. Cooling rates determined by a single-stage linear cooling model are 150–13,000 °C hr^-1 from the liquidus down to ~1000 °C, consistent with previously determined cooling rates for basaltic glasses; compositional trends with melt inclusion size observed in the Siqueiros melt inclusions are described well by this simple single-stage linear cooling model. Despite the overall success of the modeling of MgO concentration profiles using a single-stage cooling history, MgO concentration profiles in some melt inclusions are better fit by a two-stage cooling history with a slower-cooling first stage followed by a faster-cooling second stage; the inferred total duration of cooling from the liquidus down to ~1000 °C is 40 s to just over one hour.

Based on our observations and models, compositions of zoned melt inclusions (even if measured at the centers of the inclusions) will typically have been diffusively fractionated relative to the initially trapped melt; for such inclusions, the initial composition cannot be simply reconstructed based on olivine-addition calculations, so caution should be exercised in application of such reconstructions to correct for post-entrapment crystallization of olivine on inclusion walls. Off-center analyses of a melt inclusion can also give results significantly fractionated relative to simple olivine crystallization.

All melt inclusions from the Siqueiros and Galapagos sample suites exhibit zoning profiles, and this feature may be nearly universal in glassy, olivine-hosted inclusions. If so, zoning profiles in melt inclusions could be widely useful to constrain late-stage syneruptive processes and as natural diffusion experiments.

River Lune processes. A study of change in the River Lune catchment and recommendations for flood defence management

There has been a perception of increasing river channel instability in north west rivers and the River Lune in particular in recent decades. This has been attributed variously to: (a) long-term trends in precipitation-runoff regime; (b) changes in land-use such as moor-draining and sub-soil draining such that the river is more flashy than previously, and (c) a change in the magnitude-frequency relationships of flow such that high discharges are occurring with increased frequency. Resources are available in the form of rainfall and runoff records, archived information on channel planform, land use statistics and local engineering experience which have not been jointly and fully evaluated. Effective interpretation of the nature of channel change through time with respect to this resource may enhance the Environment Agency's ability to manage the river channel efficiently in the future and will aid the development of effective policy. The results of this study will for the first time, provide robust guidance with respect to long-term channel adjustment and the appropriate management options. The research provides suggestions as to how policy might be developed taking account of other pertinent factors.

Implications of using On-Farm Flood Flow Capture to recharge groundwater and mitigate flood risks along the Kings River, CA

Two large hydrologic issues face the Kings Basin, severe and chronic overdraft of about 0.16M ac-ft annually, and flood risks along the Kings River and the downstream San Joaquin River. Since 1983, these floods have caused over $1B in damage in today’s dollars. Capturing flood flows of sufficient volume could help address these two pressing issues which are relevant to many regions of the Central Valley and will only be exacerbated with climate change. However, the Kings River has high variability associated with flow magnitudes which suggests that standard engineering approaches and acquisition of sufficient acreage through purchase and easements to capture and recharge flood waters would not be cost effective. An alternative approach investigated in this study, termed On-Farm Flood Flow Capture, involved leveraging large areas of private farmland to capture flood flows for both direct and in lieu recharge. This study investigated the technical and logistical feasibility of best management practices (BMPs) associated with On-Farm Flood Flow Capture. The investigation was conducted near Helm, CA, about 20 miles west of Fresno, CA. The experimental design identified a coordinated plan to determine infiltration rates for different soil series and different crops; develop a water budget for water applied throughout the program and estimate direct and in lieu recharge; provide a preliminary assessment of potential water quality impacts; assess logistical issues associated with implementation; and provide an economic summary of the program. At check locations, we measured average infiltration rates of 4.2 in/d for all fields and noted that infiltration rates decreased asymptotically over time to about 2 – 2.5 in/d. Rates did not differ significantly between the different crops and soils tested, but were found to be about an order of magnitude higher in one field. At a 2.5 in/d infiltration rate, 100 acres are required to infiltrate 10 CFS of captured flood flows. Water quality of applied flood flows from the Kings River had concentrations of COC (constituents of concern; i.e. nitrate, electrical conductivity or EC, phosphate, ammonium, total dissolved solids or TDS) one order of magnitude or more lower than for pumped groundwater at Terranova Ranch and similarly for a broader survey of regional groundwater. Applied flood flows flushed the root zone and upper vadose zone of nitrate and salts, leading to much lower EC and nitrate concentrations to a depth of 8 feet when compared to fields in which more limited flood flows were applied or for which drip irrigation with groundwater was the sole water source. In demonstrating this technology on the farm, approximately 3,100 ac-ft was diverted, primarily from April through mid-July, with about 70% towards in lieu and 30% towards direct recharge. Substantial flood flow volumes were applied to alfalfa, wine grapes and pistachio fields. A subset of those fields, primarily wine grapes and pistachios, were used primarily to demonstrate direct recharge. For those fields about 50 – 75% of water applied was calculated going to direct recharge. Data from the check studies suggests more flood flows could have been applied and infiltrated, effectively driving up the amount of water towards direct recharge. Costs to capture flood flows for in lieu and direct recharge for this project were low compared to recharge costs for other nearby systems and in comparison to irrigating with groundwater. Moreover, the potentially high flood capture capacity of this project suggests significant flood avoidance costs savings to downstream communities along the Kings and San Joaquin Rivers. Our analyses for Terranova Ranch suggest that allocating 25% or more flood flow water towards in lieu recharge and the rest toward direct recharge will result in an economically sustainable recharge approach paid through savings from reduced groundwater pumping. Two important issues need further consideration. First, these practices are likely to leach legacy salts and nitrates from the unsaturated zone into groundwater. We develop a conceptual model of EC movement through the unsaturated zone and estimated through mass balance calculations that approximately 10 kilograms per square meter of salts will be flushed into the groundwater through displacing 12 cubic meters per square meter of unsaturated zone pore water. This flux would increase groundwater salinity but an equivalent amount of water added subsequently is predicted as needed to return to current groundwater salinity levels. All subsequent flood flow capture and recharge is expected to further decrease groundwater salinity levels. Second, the project identified important farm-scale logistical issues including irrigator training; developing cropping plans to integrate farming and recharge activities; upgrading conveyance; and quantifying results. Regional logistical issues also exist related to conveyance, integration with agricultural management, economics, required acreage and Operation and Maintenance (O&M).

On-Farm Flood Flow Capture – addressing flood risks and groundwater overdraft in the Kings Basin, with potential applications throughout the Central Valley

Project fact sheet prepared in cooperation with the USDA Natural Resources Conservation Service and the Kings River Conservation District.

The great February 1986 flood [abstract]

EXTRACT (SEE PDF FOR FULL ABSTRACT): There were many similarities between the February 1986 storm and that of December 1964 and also December 1955. The 1964 storm hit hardest a little further north and the North Coast took the brunt of that storm. December 1955 also produced higher north coastal area runoff. December 1955 produced greater peaks in the central part of the state than the 1964 flood and is perhaps more comparable south of the Lake Tahoe-American River area. But the real surprise this time was the volume. Four reservoirs, Folsom, Black Butte, Pardee, and Comanche, were filled completely and became surcharged (storing more water than the designed capacity). The 10 day total rainfall amounted to half the normal annual totals at many precipitation stations. The February 1986 flood is a vivid reminder of the extremes of California climate and the value of the extensive system of flood control works in the state. Before the storm, especially in January, there was much concern about the dryness of the water year. Then with the deluge, California's flood control systems were tested. By and large the system worked preventing untold damage and misery for most dwellers in the flat lands.

Climatic effects on flood frequency: an example from southern Arizona

EXTRACT (SEE PDF FOR FULL ABSTRACT): After 1960, the Santa Cruz River at Tucson, Arizona, an ephemeral stream normally dominated by summer floods, experienced an apparent increased frequency of flooding coincident with an increased percentage of annual floods occurring in fall and winter. This shift reflects large-scale and low-frequency changes in the eastern Pacific Ocean, in part associated with El Niño-Southern Oscillation (ENSO) phenomena. ... Questions are raised about the validity of standard methods of flood-frequency analysis to estimate regulatory and designed floods.

Holocene history of the El Niño phenomenon as recorded in flood sediments of northern coastal Peru

EXTRACT (SEE PDF FOR FULL ABSTRACT): The Holocene history of flooding in northern coastal Peru is believed to be a proxy record for the El Niño phenomenon. A recently completed set of 30 radiocarbon dates on overbank flood deposits and a tsunami deposit from the Casma region (Figure 1 and Table 1) establishes a chronology for the largest events that have occurred during the last 3500 years. ... The data presented here indicate that events much larger than the one in 1982-1983 may occur with a frequency of about once every 1000 years.

Assessment of agrobiodiversity resources in the Borotse flood plain, Zambia

Concerns about perceived loss of indigenous materials emerged from multiple stakeholders during consultations to plan and design the CGIAR Research Program on Aquatic Agricultural Systems for the Borotse hub in Zambia’s Western Province. To come to grips with and address the concerns, the AAS Borotse hub program of work included an assessment of agrobiodiversity to inform community-level and program initiatives and actions. The agrobiodiversity assessment comprised three components: key informant and expert surveys complemented by review of grey and published literature, focus group discussions in the communities, and individual household surveys. This working paper reports the findings from assessments of agrobiodiversity resources in the Borotse hub by key informants and local experts working in government ministries, departments and agencies, and non-governmental organizations operating in the communities. This working paper covers the following topics: agriculture in the Borotse flood plain; major agricultural land types in the Borotse flood plain; soils and their uses; production systems; crops, including the seed sector and ex-situ resources; indigenous materials collected from the wild, including non-perennial and perennial plants, aquatic plants, and forest biodiversity; fish resources, including both capture fisheries and aquaculture; livestock resources; dietary diversity; and indigenous and local knowledge on management systems.

Reconstructing 20th century flood patterns in Havasu Creek, Arizona, using historical and dendrochronologic data [abstract]

EXTRACT (SEE PDF FOR FULL ABSTRACT): Havasu Creek is the second largest tributary of the Colorado River in Grant Canyon. Perennial streamflow in the creek seldom exceeds 2 cubic meters per second, but it supports an important riparian habitat as well as unique travertine pools and waterfalls that attract over 20,000 tourists annually. Havasu Canyon is also home to over 400 members of the Havasu Tribe. Despite a long history of habitation and recreation in Havasu Canyon, streamflow records for Havasu Creek are extremely limited, making flood prediction difficult.

Community-led resettlement: From a flood affected slum to a new society in Pune, India

This paper describes the resettlement process of a community devastated by annual floods, to newly constructed housing in Pune, India. The relocation from Kamgar Putala slum to a housing society at Hadapsar was organized by a community-led NGO partnership in 2004. The housing development was coordinated by the local NGO Shelter Associates with significant community participation. The housing has been revisited in 2010 to evaluate the sustainability of the resettlement project’s delivery model via stakeholder perception. The process of organizing for resettlement after natural disaster is described along with the implementation and evaluation of the new housing nearly six years after initial occupation. The strong partnership approach overcame a series of political and financial hurdles at various stages of the relocation project. The story of resettling Kamgar Putala is detailed alongside an outline of the current political climate for an alternative slum-upgrading policy in India and Pune. The advantages of an empowered community supported by an influential local NGO demonstrate a commendable team effort which has tackled the threat of floods. The paper highlights the merits of a community-led partnership approach to housing development for achieving sustainable urban development as well as the alleviation of poverty in a developing context.

Community-led resettlement: From a flood-affected slum to a new society in Pune, India

This paper describes the resettlement process of a community devastated by annual floods, to newly constructed housing in Pune, India. The relocation from Kamgar Putala slum to a housing society at Hadapsar was organized by a community-led NGO partnership in 2004. The housing development was coordinated by the local NGO Shelter Associates with significant community participation. The housing has been revisited in 2010 to evaluate the sustainability of the resettlement project's delivery model via stakeholder perception. The process of organizing for resettlement after natural disaster is described along with the implementation and evaluation of the new housing nearly six years after initial occupation. The strong partnership approach overcame a series of political and financial hurdles at various stages of the relocation project. The story of resettling Kamgar Putala is detailed alongside an outline of the current political climate for an alternative slum-upgrading policy in India and Pune. The advantages of an empowered community supported by an influential local NGO demonstrate a commendable team effort which has tackled the threat of floods. The paper highlights the merits of a community-led partnership approach to housing development for achieving sustainable urban development as well as the alleviation of poverty in a developing context. © 2011 Taylor & Francis.

Digital elevation models and GIS for watershed modelling and flood prediction-A case study of Accra Ghana

Geographical Information Systems (GIS) and Digital Elevation Models (DEM) can be used to perform many geospatial and hydrological modelling including drainage and watershed delineation, flood prediction and physical development studies of urban and rural settlements. This paper explores the use of contour data and planimetric features extracted from topographic maps to derive digital elevation models (DEMs) for watershed delineation and flood impact analysis (for emergency preparedness) of part of Accra, Ghana in a GIS environment. In the study two categories of DEMs were developed with 5 m contour and planimetric topographic data; bare earth DEM and built environment DEM. These derived DEMs were used as terrain inputs for performing spatial analysis and obtaining derivative products. The generated DEMs were used to delineate drainage patterns and watershed of the study area using ArcGIS desktop and its ArcHydro extension tool from Environmental Systems Research Institute (ESRI). A vector-based approach was used to derive inundation areas at various flood levels. The DEM of built-up areas was used as inputs for determining properties which will be inundated in a flood event and subsequently generating flood inundation maps. The resulting inundation maps show that about 80% areas which have perennially experienced extensive flooding in the city falls within the predicted flood extent. This approach can therefore provide a simplified means of predicting the extent of inundation during flood events for emergency action especially in less developed economies where sophisticated technologies and expertise are hard to come by. © 2009 Springer Netherlands.