The Entrainment Rate for Buoyant Plumes in a Crossflow

We consider large-eddy simulations (LES) of buoyant plumes from a circular source with initial buoyancy flux F ₀ released into a stratified environment with constant buoyancy frequency N and a uniform crossflow with velocity U. We make a systematic comparison of the LES results with the mathematical theory of plumes in a crossflow. We pay particular attention to the limits [(U)\tilde] << 1{\tilde{U}\ll1} and [(U)\tilde] >> 1{\tilde{U}\gg 1}, where [(U)\tilde]=U/(F₀ N)^1/4{\tilde{U}=U/(F_0 N)^{1/4}}, for which analytical results are possible. For [(U)\tilde] >> 1{\tilde{U}\gg 1}, the LES results show good agreement with the well-known two-thirds law for the rise in height of the plume. Sufficiently far above the source, the centreline vertical velocity of the LES plumes is consistent with the analytical z ^−1/3 and z ^−1/2 scalings for respectively [(U)\tilde] << 1{\tilde{U}\ll 1} and [(U)\tilde] >> 1{\tilde{U}\gg 1}. In the general case, where the entrainment is assumed to be the sum of the contributions from the horizontal and vertical velocity components, we find that the discrepancy between the LES data and numerical solutions of the plume equations is largest for [(U)\tilde]=O(1){\tilde{U}=O(1)}. We propose a modified additive entrainment assumption in which the contributions from the horizontal and vertical velocity components are not equally weighted. We test this against observations of the plume generated by the Buncefield fire in the U.K. in December 2005 and find that the results compare favourably. We also show that the oscillations of the plume as it settles down to its final rise height may be attenuated by the radiation of gravity waves. For [(U)\tilde] << 1{\tilde{U}\ll 1} the oscillations decay rapidly due to the transport of energy away from the plume by gravity waves. For ${\tilde{U}>rsim 1}${\tilde{U}>rsim 1} the gravity waves travel in the same direction and at the same speed as the flow. In this case, the oscillations of the plume do not decay greatly by radiation of gravity waves.     

Ocean–atmosphere coupling and the boreal winter MJO

The influence of ocean–atmosphere coupling on the simulation and prediction of the boreal winter Madden–Julian Oscillation (MJO) is examined using the Seoul National University coupled general circulation model (CGCM) and atmospheric—only model (AGCM). The AGCM is forced with daily SSTs interpolated from pentad mean CGCM SSTs. Forecast skill is examined using serial extended simulations spanning 26 different winter seasons with 30-day forecasts commencing every 5 days providing a total of 598 30-day simulations. By comparing both sets of experiments, which share the same atmospheric components, the influence of coupled ocean–atmosphere processes on the simulation and prediction of MJO can be studied. The mean MJO intensity possesses more realistic amplitude in the CGCM than in AGCM. In general, the ocean–atmosphere coupling acts to improve the simulation of the spatio-temporal evolution of the eastward propagating MJO and the phase relationship between convection (OLR) and SST over the equatorial Indian Ocean and the western Pacific. Both the CGCM and observations exhibit a near-quadrature relationship between OLR and SST, with the former lagging by about two pentads. However, the AGCM shows a less realistic phase relationship. As the initial conditions are the same in both models, the additional forcing by SST anomalies in the CGCM extends the prediction skill beyond that of the AGCM. To test the applicability of the CGCM to real-time prediction, we compute the Real-time Multivariate MJO (RMM) index and compared it with the index computed from observations. RMM1 (RMM2) falls away rapidly to 0.5 after 17–18 (15–16) days in the AGCM and 18–19 (16–17) days in the CGCM. The prediction skill is phase dependent in both the CGCM and AGCM.     

Relative roles of climate sensitivity and forcing in defining the ocean circulation response to climate change

The response of the ocean’s meridional overturning circulation (MOC) to increased greenhouse gas forcing is examined using a coupled model of intermediate complexity, including a dynamic 3-D ocean subcomponent. Parameters are the increase in CO₂ forcing (with stabilization after a specified time interval) and the model’s climate sensitivity. In this model, the cessation of deep sinking in the north “Atlantic” (hereinafter, a “collapse”), as indicated by changes in the MOC, behaves like a simple bifurcation. The final surface air temperature (SAT) change, which is closely predicted by the product of the radiative forcing and the climate sensitivity, determines whether a collapse occurs. The initial transient response in SAT is largely a function of the forcing increase, with higher sensitivity runs exhibiting delayed behavior; accordingly, high CO₂-low sensitivity scenarios can be assessed as a recovering or collapsing circulation shortly after stabilization, whereas low CO₂-high sensitivity scenarios require several hundred additional years to make such a determination. We also systemically examine how the rate of forcing, for a given CO₂ stabilization, affects the ocean response. In contrast with previous studies based on results using simpler ocean models, we find that except for a narrow range of marginally stable to marginally unstable scenarios, the forcing rate has little impact on whether the run collapses or recovers. In this narrow range, however, forcing increases on a time scale of slow ocean advective processes results in weaker declines in overturning strength and can permit a run to recover that would otherwise collapse.     

Interpretation of Solar Magnetic Field Strength Observations

This study based on longitudinal Zeeman effect magnetograms and spectral line scans investigates the dependence of solar surface magnetic fields on the spectral line used and the way the line is sampled to estimate the magnetic flux emerging above the solar atmosphere and penetrating to the corona from magnetograms of the Mt. Wilson 150-foot tower synoptic program (MWO). We have compared the synoptic program λ5250 Å line of Fe?i to the line of Fe?i at λ5233 Å since this latter line has a broad shape with a profile that is nearly linear over a large portion of its wings. The present study uses five pairs of sampling points on the λ5233 Å line. Line profile observations show that the determination of the field strength from the Stokes V parameter or from line bisectors in the circularly polarized line profiles lead to similar dependencies on the spectral sampling of the lines, with the bisector method being the less sensitive. We recommend adoption of the field determined with the line bisector method as the best estimate of the emergent photospheric flux and further recommend the use of a sampling point as close to the line core as is practical. The combination of the line profile measurements and the cross-correlation of fields measured simultaneously with λ5250 Å and λ5233 Å yields a formula for the scale factor δ ^?1 that multiplies the MWO synoptic magnetic fields. By using ρ as the center-to-limb angle (CLA), a fit to this scale factor is δ ^?1=4.15?2.82sin?²(ρ). Previously δ ^?1=4.5?2.5sin?²(ρ) had been used. The new calibration shows that magnetic fields measured by the MDI system on the SOHO spacecraft are equal to 0.619±0.018 times the true value at a center-to-limb position 30°. Berger and Lites (2003, Solar Phys. 213, 213) found this factor to be 0.64±0.013 based on a comparison using the Advanced Stokes Polarimeter.     

Travel times for snowmelt-dominated headwater catchments: Influences of wetlands and forest harvesting,and linkages to stream water quality

The time it takes water to travel through a catchment, from when it enters as rain and snow to when it leaves as streamflow, may influence stream water quality and catchment sensitivity to environmental change. Most studies that estimate travel times do so for only a few, often rain-dominated, catchments in a region and use relatively short data records (<10 years). A better understanding of how catchment travel times vary across a landscape may help diagnose inter-catchment differences in water quality and response to environmental change. We used comprehensive and long-term observations from the Turkey Lakes Watershed Study in central Ontario to estimate water travel times for 12 snowmelt-dominated headwater catchments, three of which were impacted by forest harvesting. Chloride, a commonly used water tracer, was measured in streams, rain, snowfall and as dry atmospheric deposition over a 31 year period. These data were used with a lumped convolution integral approach to estimate mean water travel times. We explored relationships between travel times and catchment characteristics such as catchment area, slope angle, flowpath length, runoff ratio and wetland coverage, as well as the impact of harvesting. Travel time estimates were then used to compare differences in stream water quality between catchments. Our results show that mean travel times can be variable for small geographic areas and are related to catchment characteristics, in particular flowpath length and wetland cover. In addition, forest harvesting appeared to decrease mean travel times. Estimated mean travel times had complex relationships with water quality patterns. Results suggest that biogeochemical processes, particularly those present in wetlands, may have a greater influence on water quality than catchment travel times.     

Reducing dissolved inorganic nitrogen in surface runoff water from sugarcane production systems

Nitrogen (N) lost from farms, especially as the highly bioavailable dissolved inorganic form, may be damaging Australia's Great Barrier Reef (GBR). As sugarcane is the dominant cropping system in GBR catchments, its N management practises are coming under increasing scrutiny. This study measured dissolved inorganic N lost in surface runoff water and sugarcane productivity over 3 years. The experiment compared the conventional fertiliser N application rate to sugarcane (average 180kg N/ha/year) and a rate based on replacing N exported in the previous crop (average 94kg N/ha/year). Dissolved inorganic N losses in surface water were 72%, 48% and 66% lower in the three monitored years in the reduced N fertiliser treatment. There was no significant difference in sugarcane yield between the two fertiliser N treatments, nor any treatment difference in soil mineral N - both of these results are indicators of the sustainability of the lower fertiliser N applications.     

Airborne Laser-Induced Fluorescence (LIF) Light Detection and Ranging (LiDAR) for the Quantification of Dissolved Organic Matter Concentration in Natural Waters

Analysis of dissolved organic matter (DOM) concentration and composition is essential to quantifying biological and chemical oxygen demand and atmosphere?Cocean heat flux exchange in natural waters. However, manual water sampling is costly and time consuming over large areas. The purpose of this research was to analyze the applicability of airborne laser-induced fluorescence light detection and ranging (LiDAR) for the detection of DOM in estuarine ecosystems impacted by agriculture. A fluorescence LiDAR system (Airborne Marine) (FLS-AM) was used to assess the DOM concentration of the Annapolis River and Basin, Nova Scotia, Canada, as well as three rivers and their estuaries in Prince Edward Island, Canada. Two FLS-AM flight missions were conducted in the summers of 2008 and 2009 and the resulting datasets were compared with spectral fluorescence signature (SFS DOM) and dissolved organic carbon (DOC) analysis of in situ water samples. Significant positive correlations were found at five of seven sites between the FLS-AM DOM and SFS DOM relationship which indicates that the FLS-AM sensor is a good surrogate for traditional sample collection of DOM data in estuaries in this region. Positive correlations were also found at all sites between FLS-AM DOM values and DOC. FLS-AM DOM patterns show that DOM values are significantly higher in rivers and estuaries that drain watersheds which are heavily impacted by agricultural practices. The results of this study show that the FLS-AM can be used efficiently as a general indicator for how estuaries are affected by runoff from agricultural watersheds in real time and thus reduce the requirement for traditional water sample collection and laboratory analysis methods.     

Be-daughter minerals in fluid and melt inclusions: implications for the enrichment of Be in granite–pegmatite systems

The study of re-homogenized melt inclusions in the same growth planes of quartz of pegmatites genetically linked to the Variscan granite of the Ehrenfriedersdorf complex, Erzgebirge, Germany, by ion microprobe analyses has determined high concentrations of Be, up to 10,000 ppm, in one type of melt inclusion, as well as moderate concentrations in the 100 ppm range in a second type of melt inclusion. Generally, the high Be concentrations are associated with the H₂O- and other volatile-rich type-B melt inclusions, and the lower Be concentration levels are connected to H₂O-poor type-A melt inclusions. Both inclusion types, representing conjugate melt pairs, are formed by a liquid–liquid immiscibility separation process. This extremely strong and very systematic scattering in Be provides insights into the origin of Be concentration and transport mechanisms in pegmatite-forming melts. In this contribution, we present more than 250 new analytical data and show with ion microprobe and fs-LA-ICPMS studies on quenched glasses, as well as with confocal Raman spectroscopy of daughter minerals in unheated melt inclusions, that the concentrations of Be may achieve such extreme levels during melt–melt immiscibility of H₂O-, B-, F-, P-, ± Li-enriched pegmatite-forming magmas. Starting from host granite with about 10 ppm Be, melt inclusions with 10,000 ppm Be correspond to enrichment by a factor of over 1,000. This strong enrichment of Be is the result of processes of fractional crystallization and further enrichment in melt patches of pegmatite bodies due to melt–melt immiscibility at fluid saturation. We also draw additional conclusions regarding the speciation of Be in pegmatite-forming melt systems from investigation of the Be-bearing daughter mineral phases in the most H₂O-rich melt inclusions. In the case of evolved volatile and H₂O-rich pegmatite systems, B, P, and carbonates are important for the enrichment and formation of stable Be complexes.     

Tsunami simulation with Green–Naghdi theory

Green–Naghdi (GN) theory is a fully nonlinear wave theory which has been used with success to simulate nonlinear water waves. In previous applications of GN theory to water wave problems the ocean bottom was assumed to be time invariant. In this work no such restriction is made and GN theory is used to simulate tsunami caused by bottom fluctuation. As first test cases we simulate two-dimensional nonlinear surface waves generated by positive bottom movements. The results in the generation region for three different seabed movements compare well against earlier experimental data. The results in the downstream region for impulsive seabed movements show some discrepancies in wave phase and amplitude compared with earlier experimental values. It is suspected that the viscous effects may have played a role. The GN theory is then used to study three-dimensional near-field tsunami amplitudes caused by submarine landslides and slumps spreading in two orthogonal directions. The GN results agree with previous linear solution very well when the ratio of the velocities is v₁/v₂=1.0. But GN theory give more believable results for the case of v_T/v=0.1 and v₁/v₂=0.1.