Modelling of urban sensible heat flux at multiple spatial scales : a demonstration using airborne hyperspectral imagery of Shanghai and a temperature–emissivity separation approach

Urbanization related alterations to the surface energy balance impact urban warming (‘heat islands’), the growth of the boundary layer, and many other biophysical processes. Traditionally, in situ heat flux measures have been used to quantify such processes, but these typically represent only a small local-scale area within the heterogeneous urban environment. For this reason, remote sensing approaches are very attractive for elucidating more spatially representative information. Here we use hyperspectral imagery from a new airborne sensor, the Operative Modular Imaging Spectrometer (OMIS), along with a survey map and meteorological data, to derive the land cover information and surface parameters required to map spatial variations in turbulent sensible heat flux (QH). The results from two spatially-explicit flux retrieval methods which use contrasting approaches and, to a large degree, different input data are compared for a central urban area of Shanghai, China: (1) the Local-scale Urban Meteorological Parameterization Scheme (LUMPS) and (2) an Aerodynamic Resistance Method (ARM). Sensible heat fluxes are determined at the full 6 m spatial resolution of the OMIS sensor, and at lower resolutions via pixel aggregation and spatial averaging. At the 6 m spatial resolution, the sensible heat flux of rooftop dominated pixels exceeds that of roads, water and vegetated areas, with values peaking at ∼ 350 W m− 2, whilst the storage heat flux is greatest for road dominated pixels (peaking at around 420 W m− 2). We investigate the use of both OMIS-derived land surface temperatures made using a Temperature–Emissivity Separation (TES) approach, and land surface temperatures estimated from air temperature measures. Sensible heat flux differences from the two approaches over the entire 2 × 2 km study area are less than 30 W m− 2, suggesting that methods employing either strategy maybe practica1 when operated using low spatial resolution (e.g. 1 km) data. Due to the differing methodologies, direct comparisons between results obtained with the LUMPS and ARM methods are most sensibly made at reduced spatial scales. At 30 m spatial resolution, both approaches produce similar results, with the smallest difference being less than 15 W m− 2 in mean QH averaged over the entire study area. This is encouraging given the differing architecture and data requirements of the LUMPS and ARM methods. Furthermore, in terms of mean study QH, the results obtained by averaging the original 6 m spatial resolution LUMPS-derived QH values to 30 and 90 m spatial resolution are within ∼ 5 W m− 2 of those derived from averaging the original surface parameter maps prior to input into LUMPS, suggesting that that use of much lower spatial resolution spaceborne imagery data, for example from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) is likely to be a practical solution for heat flux determination in urban areas.

An assessment of observed vertical flux divergence in long-term eddy-covariance measurements over two Midwestern forest ecosystems

Vertical divergence of CO2 fluxes is observed over two Midwestern AmeriFlux forest sites. The differences in ensemble averaged hourly CO2 fluxes measured at two heights above canopy are relatively small (0.2–0.5 μmol m−2 s−1), but they are the major contributors to differences (76–256 g C m−2 or 41.8–50.6%) in estimated annual net ecosystem exchange (NEE) in 2001. A friction velocity criterion is used in these estimates but mean flow advection is not accounted for. This study examines the effects of coordinate rotation, averaging time period, sampling frequency and co-spectral correction on CO2 fluxes measured at a single height, and on vertical flux differences measured between two heights. Both the offset in measured vertical velocity and the downflow/upflow caused by supporting tower structures in upwind directions lead to systematic over- or under-estimates of fluxes measured at a single height. An offset of 1 cm s−1 and an upflow/downflow of 1° lead to 1% and 5.6% differences in momentum fluxes and nighttime sensible heat and CO2 fluxes, respectively, but only 0.5% and 2.8% differences in daytime sensible heat and CO2 fluxes. The sign and magnitude of both offset and upflow/downflow angle vary between sonic anemometers at two measurement heights. This introduces a systematic and large bias in vertical flux differences if these effects are not corrected in the coordinate rotation. A 1 h averaging time period is shown to be appropriate for the two sites. In the daytime, the absolute magnitudes of co-spectra decrease with height in the natural frequencies of 0.02–0.1 Hz but increase in the lower frequencies (<0.01 Hz). Thus, air motions in these two frequency ranges counteract each other in determining vertical flux differences, whose magnitude and sign vary with averaging time period. At night, co-spectral densities of CO2 are more positive at the higher levels of both sites in the frequency range of 0.03–0.4 Hz and this vertical increase is also shown at most frequencies lower than 0.03 Hz. Differences in co-spectral corrections at the two heights lead to a positive shift in vertical CO2 flux differences throughout the day at both sites. At night, the vertical CO2 flux differences between two measurement heights are 20–30% and 40–60% of co-spectral corrected CO2 fluxes measured at the lower levels of the two sites, respectively. Vertical differences of CO2 flux are relatively small in the daytime. Vertical differences in estimated mean vertical advection of CO2 between the two measurement heights generally do not improve the closure of the 1D (vertical) CO2 budget in the air layer between the two measurement heights. This may imply the significance of horizontal advection. However, a reliable assessment of mean advection contributions in annual NEE estimate at these two AmeriFlux sites is currently an unsolved problem.

A critical revision of the estimation of the latent heat flux from two-wavelength scintillometry

Simultaneous scintillometer measurements at multiple wavelengths (pairing visible or infrared with millimetre or radio waves) have the potential to provide estimates of path-averaged surface fluxes of sensible and latent heat. Traditionally, the equations to deduce fluxes from measurements of the refractive index structure parameter at the two wavelengths have been formulated in terms of absolute humidity. Here, it is shown that formulation in terms of specific humidity has several advantages. Specific humidity satisfies the requirement for a conserved variable in similarity theory and inherently accounts for density effects misapportioned through the use of absolute humidity. The validity and interpretation of both formulations are assessed and the analogy with open-path infrared gas analyser density corrections is discussed. Original derivations using absolute humidity to represent the influence of water vapour are shown to misrepresent the latent heat flux. The errors in the flux, which depend on the Bowen ratio (larger for drier conditions), may be of the order of 10%. The sensible heat flux is shown to remain unchanged. It is also verified that use of a single scintillometer at optical wavelengths is essentially unaffected by these new formulations. Where it may not be possible to reprocess two-wavelength results, a density correction to the latent heat flux is proposed for scintillometry, which can be applied retrospectively to reduce the error.

The gravity wave momentum flux in hydrostatic flow with directional shear over elliptical mountains

Semi-analytical expressions for the momentum flux associated with orographic internal gravity waves, and closed analytical expressions for its divergence, are derived for inviscid, stationary, hydrostatic, directionally-sheared flow over mountains with an elliptical horizontal cross-section. These calculations, obtained using linear theory conjugated with a third-order WKB approximation, are valid for relatively slowly-varying, but otherwise generic wind profiles, and given in a form that is straightforward to implement in drag parametrization schemes. When normalized by the surface drag in the absence of shear, a quantity that is calculated routinely in existing drag parametrizations, the momentum flux becomes independent of the detailed shape of the orography. Unlike linear theory in the Ri → ∞ limit, the present calculations account for shear-induced amplification or reduction of the surface drag, and partial absorption of the wave momentum flux at critical levels. Profiles of the normalized momentum fluxes obtained using this model and a linear numerical model without the WKB approximation are evaluated and compared for two idealized wind profiles with directional shear, for different Richardson numbers (Ri). Agreement is found to be excellent for the first wind profile (where one of the wind components varies linearly) down to Ri = 0.5, while not so satisfactory, but still showing a large improvement relative to the Ri → ∞ limit, for the second wind profile (where the wind turns with height at a constant rate keeping a constant magnitude). These results are complementary, in the Ri > O(1) parameter range, to Broad’s generalization of the Eliassen–Palm theorem to 3D flow. They should contribute to improve drag parametrizations used in global weather and climate prediction models.

Reconstruction of geomagnetic activity and near-Earth interplanetary conditions over the past 167 yr - Part 4: Near-Earth solar wind speed, IMF, and open solar flux

In the concluding paper of this tetralogy, we here use the different geomagnetic activity indices to reconstruct the near-Earth interplanetary magnetic field (IMF) and solar wind flow speed, as well as the open solar flux (OSF) from 1845 to the present day. The differences in how the various indices vary with near-Earth interplanetary parameters, which are here exploited to separate the effects of the IMF and solar wind speed, are shown to be statistically significant at the 93% level or above. Reconstructions are made using four combinations of different indices, compiled using different data and different algorithms, and the results are almost identical for all parameters. The correction to the aa index required is discussed by comparison with the Ap index from a more extensive network of mid-latitude stations. Data from the Helsinki magnetometer station is used to extend the aa index back to 1845 and the results confirmed by comparison with the nearby St Petersburg observatory. The optimum variations, using all available long-term geomagnetic indices, of the near-Earth IMF and solar wind speed, and of the open solar flux, are presented; all with ±2sigma� uncertainties computed using the Monte Carlo technique outlined in the earlier papers. The open solar flux variation derived is shown to be very similar indeed to that obtained using the method of Lockwood et al. (1999).

Factors contributing to fecal incontinence in older people and outcome of routine management in home, hospital and nursing home settings

Objective: Fecal loading, cognitive impairment, loose stools, functional disability, comorbidity and anorectal incontinence are recognized as factors contributing to loss of fecal continence in older adults. The objective of this project was to assess the relative distribution of these factors in a variety of settings along with the outcome of usual management. Methods: One hundred and twenty adults aged 65 years and over with fecal incontinence recruited by convenience sampling from four different settings were studied. They were either living at home or in a nursing home or receiving care on an acute or rehabilitation elderly care ward. A structured questionnaire was used to elicit which factors associated with fecal incontinence were present from subjects who had given written informed consent or for whom assent for inclusion in the study had been obtained. Results: Fecal loading (Homes 6 [20%]; Acute care wards 17 [57%]; Rehabilitation wards 19 [63%]; Nursing homes 21 [70%]) and functional disability (Homes 5 [17%]; Acute care wards 25 [83%]; Rehabilitation wards 25 [83%]; Nursing homes 20 [67%]) were significantly more prevalent in the hospital and nursing home settings than in those living at home (P < 0.01). Loose stools were more prevalent in the hospital setting than in the other settings (Homes 11 [37%]; Acute care wards 20 [67%]; Rehabilitation wards 17 [57%]; Nursing homes 6 [20%]) (P < 0.01). Cognitive impairment was significantly more common in the nursing home than in the other settings (Nursing homes 26 [87%], Homes 5 [17%], Acute care wards 13 [43%], Rehabilitation wards 14 [47%]) (P < 0.01). Loose stools were the most prevalent factor present at baseline in 13 of the 19 (68%) subjects whose fecal incontinence had resolved at 3 months. Conclusion: The distribution of the factors contributing to fecal incontinence in older people living at home differs from those cared for in nursing home and hospital wards settings. These differences need to be borne in mind when assessing people in different settings. Management appears to result in a cure for those who are not significantly disabled with loose stools as a cause for their fecal incontinence, but this would need to be confirmed by further research.

Use of pelleted sericea lespedeza (Lespedeza cuneata) for natural control of coccidia and gastrointestinal nematodes in weaned goats

Infection with Eimeria spp. (coccidia) can be devastating in goats, particularly for young, recently-weaned kids, resulting in diarrhea, dehydration, and even death. Feeding dried sericea lespedeza [SL; Lespedeza cuneata (Dum.-Cours.) G. Don.] to young goats has been reported to reduce the effects of internal parasites, including gastrointestinal nematodes (GIN) but there have been no reports of the effects of feeding this forage on Eimeria spp. in goats. Two confinement feeding experiments were completed on recently-weaned intact bucks (24 Kiko-cross, Exp. 1; 20 Spanish, Exp. 2) to determine effects of SL pellets on an established infection of GIN and coccidia. The bucks were assigned to 1 of 2 (Exp. 1) or 3 (Exp. 2) treatment groups based upon the number of Eimeria spp. oocysts per gram (OPG) of feces. In Exp. 1, the kids were fed 1 of 2 pelleted rations ad libitum; 90% SL leaf meal + 10% of a liquid molasses/lignin binder mix and a commercial pellet with 12% crude protein (CP) and 24% acid detergent fiber (n = 12/treatment group, 2 animals/pen). For Exp. 2, treatment groups were fed 1) 90% SL leaf meal pellets from leaves stored 3 years (n = 7), 2) 90% SL pellets from leaf meal stored less than 6 months, (n = 7), and the commercial pellets (n = 6) ad libitum. For both trials, fecal and blood samples were taken from individual animals every 7 days for 28 days to determine OPG and GIN eggs per gram (EPG) and packed cell volume (PCV), respectively. In Exp. 2, feces were scored for consistency (1 = solid pellets, 5 = slurry) as an indicator of coccidiosis. In Exp. 1, EPG (P < 0.001) and OPG (P < 0.01) were reduced by 78.7 and 96.9%, respectively, 7 days after initiation of feeding in goats on the SL pellet diet compared with animals fed the control pellets. The OPG and EPG remained lower in treatment than control animals until the end of the trial. In Exp. 2, goats fed new and old SL leaf meal pellets had 66.2 and 79.2% lower (P < 0.05) EPG and 92.2 and 91.2% lower (P < 0.05) OPG, respectively, than control animals within 7 days, and these differences were maintained or increased throughout the trial. After 4 weeks of pellet feeding in Exp. 2, fecal scores were lower (P < 0.01) in both SL-fed groups compared with control animals, indicating fewer signs of coccidiosis. There was no effect of diet on PCV values throughout either experiment. Dried, pelleted SL has excellent potential as a natural anti-coccidial feed for weaned goats.

Centennial variations in sunspot number, open solar flux and streamer belt width: 1. correction of the sunspot number record since 1874

We analyse the widely-used international/ Zürich sunspot number record, R, with a view to quantifying a suspected calibration discontinuity around 1945 (which has been termed the “Waldmeier discontinuity” [Svalgaard, 2011]). We compare R against the composite sunspot group data from the Royal Greenwich Observatory (RGO) network and the Solar Optical Observing Network (SOON), using both the number of sunspot groups, N{sub}G{\sub}, and the total area of the sunspots, A{sub}G{\sub}. In addition, we compare R with the recently developed interdiurnal variability geomagnetic indices IDV and IDV(1d). In all four cases, linearity of the relationship with R is not assumed and care is taken to ensure that the relationship of each with R is the same before and after the putative calibration change. It is shown the probability that a correction is not needed is of order 10{sup}−8{\sup} and that R is indeed too low before 1945. The optimum correction to R for values before 1945 is found to be 11.6%, 11.7%, 10.3% and 7.9% using A{sub}G{\sub}, N{sub)G{\sub}, IDV, and IDV(1d), respectively. The optimum value obtained by combining the sunspot group data is 11.6% with an uncertainty range 8.1-14.8% at the 2σ level. The geomagnetic indices provide an independent yet less stringent test but do give values that fall within the 2σ uncertainty band with optimum values are slightly lower than from the sunspot group data. The probability of the correction needed being as large as 20%, as advocated by Svalgaard [2011], is shown to be 1.6 × 10{sup}−5{\sup}.

Centennial variations in sunspot number, open solar flux and streamer belt width: 2. Comparison with the geomagnetic data

We investigate the relationship between interdiurnal variation geomagnetic activity indices, IDV and IDV(1d), corrected sunspot number, R{sub}C{\sub}, and the group sunspot number R{sub}G{\sub}. R{sub}C{\sub} uses corrections for both the “Waldmeier discontinuity”, as derived in Paper 1 [Lockwood et al., 2014c], and the “Wolf discontinuity” revealed by Leussu et al. [2013]. We show that the simple correlation of the geomagnetic indices with R{sub}C{\sub}{sup}n{\sup} or R{sub}G{\sub}{sup}n{\sup} masks a considerable solar cycle variation. Using IDV(1d) or IDV to predict or evaluate the sunspot numbers, the errors are almost halved by allowing for the fact that the relationship varies over the solar cycle. The results indicate that differences between R{sub}C{\sub} and R{sub}G{\sub} have a variety of causes and are highly unlikely to be attributable to errors in either R{sub}G{\sub} alone, as has recently been assumed. Because it is not known if R{sub}C{\sub} or R{sub}G{\sub} is a better predictor of open flux emergence before 1874, a simple sunspot number composite is suggested which, like R{sub}G{\sub}, enables modelling of the open solar flux for 1610 onwards in Paper 3, but maintains the characteristics of R{sub}C{\sub}.

Centennial variations in sunspot number, open solar flux and streamer belt width: 3. Modelling

From the variation of near-Earth interplanetary conditions, reconstructed for the mid-19th century to the present day using historic geomagnetic activity observations, Lockwood and Owens [2014] have suggested that Earth remains within a broadened streamer belt during solar cycles when the Open Solar Flux (OSF) is low. From this they propose that the Earth was immersed in almost constant slow solar wind during the Maunder minimum (c. 1650-1710). In this paper, we extend continuity modelling of the OSF to predict the streamer belt width using both group sunspot numbers and corrected international sunspot numbers to quantify the emergence rate of new OSF. The results support the idea that the solar wind at Earth was persistently slow during the Maunder minimum because the streamer belt was broad.

Attribution of the spatial pattern of CO2-forced sea level change to ocean surface flux changes

Climate models taking part in the coupled model intercomparison project phase 5 (CMIP5) all predict a global mean sea level rise for the 21st century. Yet the sea level change is not spatially uniform and differs among models. Here we evaluate the role of air–sea fluxes of heat, water and momentum (windstress) to find the spatial pattern associated to each of them as well as the spread they can account for. Using one AOGCM to which we apply the surface flux changes from other AOGCMs, we show that the heat flux and windstress changes dominate both the pattern and the spread, but taking the freshwater flux into account as well yields a sea level change pattern in better agreement with the CMIP5 ensemble mean. Differences among the CMIP5 control ocean temperature fields have a smaller impact on the sea level change pattern.

Effect of fall-grazed sericea lespedeza (Lespedeza cuneata) on gastrointestinal nematode infections of growing goats

High prevalence of anthelmintic-resistant gastrointestinal nematodes (GIN) in goats has increased pressure to find effective, alternative non-synthetic control methods, one of which is adding forage of the high condensed tannin (CT) legume sericea lespedeza (SL; Lespedeza cuneata) to the animal's diet. Previous work has demonstrated good efficacy of dried SL (hay, pellets) against small ruminant GIN, but information is lacking on consumption of fresh SL, particularly during the late summer–autumn period in the southern USA when perennial warm-season grass pastures are often low in quality. A study was designed to determine the effects of autumn (September–November) consumption of fresh SL forage, grass pasture (predominantly bermudagrass, BG; Cynodon dactylon), or a combination of SL + BG forage by young goats [intact male Spanish kids, 9 months old (20.7 ± 1.1 kg), n = 10/treatment group] on their GIN infection status. Three forage paddocks (0.40 ha) were set up at the Fort Valley State University Agricultural Research Station (Fort Valley, GA) for an 8-week trial. The goats in each paddock were supplemented with a commercial feed pellet at 0.45 kg/head/d for the first 4 weeks of the trial, and 0.27 kg/head/d for the final 4 weeks. Forage samples taken at the start of the trial were analyzed for crude protein (CP), neutral detergent fiber (NDF), and acid detergent fiber (ADF) content, and a separate set of SL samples was analyzed for CT in leaves, stems, and whole plant using the benzyl mercaptan thiolysis method. Animal weights were taken at the start and end of the trial, and fecal and blood samples were collected weekly for determination of fecal egg counts (FEC) and packed cell volume (PCV), respectively. Adult GIN was recovered from the abomasum and small intestines of all goats at the end of the experiment for counting and speciation. The CP levels were highest for SL forage, intermediate for SL + BG, and lowest for BG forage samples, while NDF and ADF values were the opposite, with highest levels in BG and lowest in SL forage samples. Sericea lespedeza leaves had more CT than stems (16.0 g vs. 3.3 g/100 g dry weight), a slightly higher percentage of PDs (98% vs. 94%, respectively) and polymers of larger mean degrees of polymerization (42 vs. 18, respectively). There were no differences in average daily gain or blood PCV between the treatment groups, but SL goats had lower FEC (P < 0.05) than the BG or SL + BG forage goats throughout most of the trial. The SL + BG goats had lower FEC than the BG forage animals by the end of the trial (week 8, P < 0.05). The SL goats had lower numbers (P < 0.05) of male Haemonchus contortus and tended to have fewer female (P < 0.10) and total (P < 0.07) H. contortus compared with the BG goats. The predominant GIN in all the goats was Trichostrongylus colubriformis (73% of total GIN). As a low-input forage with activity against pathogenic GIN (H. contortus), SL has a potential to reduce producers’ dependence upon synthetic anthelmintics and also to fill the autumn ‘window’ in good-quality fresh forages for goat grazing in the southern USA.

Centennial changes in the solar wind speed and in the open solar flux

We use combinations of geomagnetic indices, based on both variation range and hourly means, to derive the solar wind flow speed, the interplanetary magnetic field strength at 1 AU and the total open solar flux between 1895 and the present. We analyze the effects of the regression procedure and geomagnetic indices used by adopting four analysis methods. These give a mean interplanetary magnetic field strength increase of 45.1 ± 4.5% between 1903 and 1956, associated with a 14.4 ± 0.7% rise in the solar wind speed. We use averaging timescales of 1 and 2 days to allow for the difference between the magnetic fluxes threading the coronal source surface and the heliocentric sphere at 1 AU. The largest uncertainties originate from the choice of regression procedure: the average of all eight estimates of the rise in open solar flux is 73.0 ± 5.0%, but the best procedure, giving the narrowest and most symmetric distribution of fit residuals, yields 87.3 ± 3.9%.

An optimal inverse method using Doppler lidar measurements to estimate the surface sensible heat flux

Inverse methods are widely used in various fields of atmospheric science. However, such methods are not commonly used within the boundary-layer community, where robust observations of surface fluxes are a particular concern. We present a new technique for deriving surface sensible heat fluxes from boundary-layer turbulence observations using an inverse method. Doppler lidar observations of vertical velocity variance are combined with two well-known mixed-layer scaling forward models for a convective boundary layer (CBL). The inverse method is validated using large-eddy simulations of a CBL with increasing wind speed. The majority of the estimated heat fluxes agree within error with the proscribed heat flux, across all wind speeds tested. The method is then applied to Doppler lidar data from the Chilbolton Observatory, UK. Heat fluxes are compared with those from a mast-mounted sonic anemometer. Errors in estimated heat fluxes are on average 18 %, an improvement on previous techniques. However, a significant negative bias is observed (on average −63%) that is more pronounced in the morning. Results are improved for the fully-developed CBL later in the day, which suggests that the bias is largely related to the choice of forward model, which is kept deliberately simple for this study. Overall, the inverse method provided reasonable flux estimates for the simple case of a CBL. Results shown here demonstrate that this method has promise in utilizing ground-based remote sensing to derive surface fluxes. Extension of the method is relatively straight-forward, and could include more complex forward models, or other measurements.