Simulations and mechanisms of subtropical low-cloud response to climate change

This thesis focuses on improving the simulation skills and the theoretical understanding of the subtropical low cloud response to climate change.

First, an energetically consistent forcing framework is designed and implemented for the large eddy simulation (LES) of the low-cloud response to climate change. The three representative current-day subtropical low cloud regimes of cumulus (Cu), cumulus-over-stratocumulus, and stratocumulus (Sc) are all well simulated with this framework, and results are comparable to the conventional fixed-SST approach. However, the cumulus response to climate warming subject to energetic constraints differs significantly from the conventional approach with fixed SST. Under the energetic constraint, the subtropics warm less than the tropics, since longwave (LW) cooling is more efficient with the drier subtropical free troposphere. The surface latent heat flux (LHF) also increases only weakly subject to the surface energetic constraint. Both factors contribute to an increased estimated inversion strength (EIS), and decreased inversion height. The decreased Cu-depth contributes to a decrease of liquid water path (LWP) and weak positive cloud feedback. The conventional fixed-SST approach instead simulates a strong increase in LHF and deepening of the Cu layer, leading to a weakly negative cloud feedback. This illustrates the importance of energetic constraints to the simulation and understanding of the sign and magnitude of low-cloud feedback.

Second, an extended eddy-diffusivity mass-flux (EDMF) closure for the unified representation of sub-grid scale (SGS) turbulence and convection processes in general circulation models (GCM) is presented. The inclusion of prognostic terms and the elimination of the infinitesimal updraft fraction assumption makes it more flexible for implementation in models across different scales. This framework can be consistently extended to formulate multiple updrafts and downdrafts, as well as variances and covariances. It has been verified with LES in different boundary layer regimes in the current climate, and further development and implementation of this closure may help to improve our simulation skills and understanding of low-cloud feedback through GCMs.

Estimating relative abundance from catch and effort data, using neural networks

We develop and test a method to estimate relative abundance from catch and effort data using neural networks. Most stock assessment models use time series of relative abundance as their major source of information on abundance levels. These time series of relative abundance are frequently derived from catch-per-unit-of-effort (CPUE) data, using general linearized models (GLMs). GLMs are used to attempt to remove variation in CPUE that is not related to the abundance of the population. However, GLMs are restricted in the types of relationships between the CPUE and the explanatory variables. An alternative approach is to use structural models based on scientific understanding to develop complex non-linear relationships between CPUE and the explanatory variables. Unfortunately, the scientific understanding required to develop these models may not be available. In contrast to structural models, neural networks uses the data to estimate the structure of the non-linear relationship between CPUE and the explanatory variables. Therefore neural networks may provide a better alternative when the structure of the relationship is uncertain. We use simulated data based on a habitat based-method to test the neural network approach and to compare it to the GLM approach. Cross validation and simulation tests show that the neural network performed better than nominal effort and the GLM approach. However, the improvement over GLMs is not substantial. We applied the neural network model to CPUE data for bigeye tuna (Thunnus obesus) in the Pacific Ocean.

The annual course of precipitation over much of the United States: observed versus GCM simulation

General Circulation Models (GCMs) may be useful in estimating the ecological impacts of global climatic change. We analyzed seasonal weather patterns over the conterminous United States and determined that regional patterns of rainfall seasonality appear to control the distributions of the Nation's major biomes. These regional patterns were compared to the output from three GCMs for validation. The models appear to simulate the appropriate seasonal climates in the northern tier of states. However, the spatial extent of these regions is distorted. None of the models accurately portrayed rainfall seasonalities in the southern tier of states, where biomes are primarily influenced by the Bermuda High.

GCM-data intercomparison: the good news and the bad

EXTRACT (SEE PDF FOR FULL ABSTRACT): General circulation models (GCMs) are probably the most sophisticated theoretical tools we have to simulate possible climatic effects of increasing carbon dioxide and other greenhouse gases. ... As will be illustrated here using a variety of examples, although the models do simulate "reality" very well on the "grand" scale (e.g., global, hemispheric, zonal), substantial differences are more apparent as the scale is reduced to areas particularly relevant to regional planners. It is particularly important that workers more clearly recognize the potential dangers in relying too heavily on simple summary statistics such as averages estimated over large regional scales. Many shortcomings are apparent in the model simulations of the present climate, indicating that further model improvements are needed to achieve reliable regional and seasonal projections of the future climatic conditions.

Design and process selection for joints in flexible materials

The application of high performance textiles has grown significantly in the last 10 to 15 years. Various research groups throughout the United Kingdom, such as the Department of Trade and Industry, have identified technical textiles as a field for future development. There is little design guidance for joining of flexible materials or general property models that can be applied to theses materials. This lack is due to the large diversity of properties, structures and resulting behaviours of the materials that are classified as "Flexible Materials". This dissertation explores the issues that are involved in characterising the materials at the fibre, bulk and textile levels. Different units of measurement are used for each stage of the manufacturing process of flexible materials and this disparity creates problems when trying to make general comparisons (e.g. comparing textiles to polymer films). Thus, a possible solution to this is to create selection charts that allow designers to compare the strength of materials for a given mass per unit area. A design tool was created using the Cambridge Engineering Selector (CES) software to enable the selection of joining processes for material. The tool is effective in selecting a reduced number of viable joining processes. Through case studies it was shown that designers are required to examine the selected processes (identified by the software) in greater detail - in particular the economics and geometry of the joint - in order to identify the optimum joining process.

On the use of combined compact scheme for numerical solution of the two-layer oceanic shallow water equations

Many types of oceanic physical phenomena have a wide range in both space and time. In general, simplified models, such as shallow water model, are used to describe these oceanic motions. The shallow water equations are widely applied in various oceanic and atmospheric extents. By using the two-layer shallow water equations, the stratification effects can be considered too. In this research, the sixth-order combined compact method is investigated and numerically implemented as a high-order method to solve the two-layer shallow water equations. The second-order centered, fourth-order compact and sixth-order super compact finite difference methods are also used to spatial differencing of the equations. The first part of the present work is devoted to accuracy assessment of the sixth-order super compact finite difference method (SCFDM) and the sixth-order combined compact finite difference method (CCFDM) for spatial differencing of the linearized two-layer shallow water equations on the Arakawa's A-E and Randall's Z numerical grids. Two general discrete dispersion relations on different numerical grids, for inertia-gravity and Rossby waves, are derived. These general relations can be used for evaluation of the performance of any desired numerical scheme. For both inertia-gravity and Rossby waves, minimum error generally occurs on Z grid using either the sixth-order SCFDM or CCFDM methods. For the Randall's Z grid, the sixth-order CCFDM exhibits a substantial improvement , for the frequency of the barotropic and baroclinic modes of the linear inertia-gravity waves of the two layer shallow water model, over the sixth-order SCFDM. For the Rossby waves, the sixth-order SCFDM shows improvement, for the barotropic and baroclinic modes, over the sixth-order CCFDM method except on Arakawa's C grid. In the second part of the present work, the sixth-order CCFDM method is used to solve the one-layer and two-layer shallow water equations in their nonlinear form. In one-layer model with periodic boundaries, the performance of the methods for mass conservation is compared. The results show high accuracy of the sixth-order CCFDM method to simulate a complex flow field. Furthermore, to evaluate the performance of the method in a non-periodic domain the sixth-order CCFDM is applied to spatial differencing of vorticity-divergence-mass representation of one-layer shallow water equations to solve a wind-driven current problem with no-slip boundary conditions. The results show good agreement with published works. Finally, the performance of different schemes for spatial differencing of two-layer shallow water equations on Z grid with periodic boundaries is investigated. Results illustrate the high accuracy of combined compact method.

The adoption of sustainable remediation behaviour in the US and UK: a cross country comparison and determinant analysis.

The sustainable remediation concept, aimed at maximizing the net environmental, social, and economic benefits in contaminated site remediation, is being increasingly recognized by industry, governments, and academia. However, there is limited understanding of actual sustainable behaviour being adopted and the determinants of such sustainable behaviour. The present study identified 27 sustainable practices in remediation. An online questionnaire survey was used to rank and compare them in the US (n=112) and the UK (n=54). The study also rated ten promoting factors, nine barriers, and 17 types of stakeholders' influences. Subsequently, factor analysis and general linear models were used to determine the effects of internal characteristics (i.e. country, organizational characteristics, professional role, personal experience and belief) and external forces (i.e. promoting factors, barriers, and stakeholder influences). It was found that US and UK practitioners adopted many sustainable practices to similar extents. Both US and UK practitioners perceived the most effectively adopted sustainable practices to be reducing the risk to site workers, protecting groundwater and surface water, and reducing the risk to the local community. Comparing the two countries, we found that the US adopted innovative in-situ remediation more effectively; while the UK adopted reuse, recycling, and minimizing material usage more effectively. As for the overall determinants of sustainable remediation, the country of origin was found not to be a significant determinant. Instead, organizational policy was found to be the most important internal characteristic. It had a significant positive effect on reducing distant environmental impact, sustainable resource usage, and reducing remediation cost and time (p<0.01). Customer competitive pressure was found to be the most extensively significant external force. In comparison, perceived stakeholder influence, especially that of primary stakeholders (site owner, regulator, and primary consultant), did not appear to have as extensive a correlation with the adoption of sustainability as one would expect.

未来气候变化对黄土高原黑河流域水资源的影响

气候变化对黄土高原的水资源有重要影响,对其影响进行评估可以为区域发展提供重要的决策依据。基于分布式水文模型SWAT和4种全球环流模式的各3种排放情景,评估了2010～2039年黄土高塬沟壑区黑河流域水资源对气候变化的潜在响应。结果表明,黑河流域2010～2039年的年均降水变化-2.3%～7.8%,年均最高和最低温度分别升高0.7～2.2℃和1.2～2.8℃,年均径流量变化-19.8%～37.0%,1.2m剖面年均土壤水分含量变化-5.5%～17.2%,年均蒸散量普遍增长0.1%～5.9%;水文气象变量变化趋势复杂,但T检验表明年降水、径流、土壤水分和蒸散增长的概率较大。对于季节变化,降水可能在12～7月份和9月份增长,8月份和10～11月份减少;径流在4～7月份和9～10月份增加,11～3月份和8月份减少;土壤水分在各月都增长;蒸散11～6月份普遍增长,7～10月份减少的可能性较大。未来气候将发生显著变化并对水资源有重要影响,需采取必要的措施来减缓其不利影响。

A global shallow water model using high order multi-moment constrained finite volume method and icosahedral grid

A novel accurate numerical model for shallow water equations on sphere have been developed by implementing the high order multi-moment constrained finite volume (MCV) method on the icosahedral geodesic grid. High order reconstructions are conducted cell-wisely by making use of the point values as the unknowns distributed within each triangular cell element. The time evolution equations to update the unknowns are derived from a set of constrained conditions for two types of moments, i.e. the point values on the cell boundary edges and the cell-integrated average. The numerical conservation is rigorously guaranteed. in the present model, all unknowns or computational variables are point values and no numerical quadrature is involved, which particularly benefits the computational accuracy and efficiency in handling the spherical geometry, such as coordinate transformation and curved surface. Numerical formulations of third and fourth order accuracy are presented in detail. The proposed numerical model has been validated by widely used benchmark tests and competitive results are obtained. The present numerical framework provides a promising and practical base for further development of atmospheric and oceanic general circulation models. (C) 2009 Elsevier Inc. All rights reserved.

土壤侵蚀高维指数非线性模型及其参数辨识

模型验证和数据库组建基础上,用WinEPIC模型定量模拟研究了黄土高原半湿润区长武、半干旱区固原和半干旱偏旱区海原20～30年内苜蓿草地水分生产潜力、10m土层土壤有效含水量和土壤湿度剖面分布特征的动态变化.结果表明:长武、固原和海原苜蓿草地水分生产潜力模拟值随降水量变化而呈现波动性降低趋势,其平均值分别为8.81、3.83和2.48t.hm-2;长武、固原和海原苜蓿草地10m土层逐月土壤有效含水量模拟值均呈现明显的波动性降低趋势,模拟初期,4～8年生苜蓿草地土壤干燥化趋势十分强烈,此后,随降水量变化长期在较低水平上波动;随着苜蓿生长年限的延长,苜蓿草地土壤干层逐年加深、加厚,长武、固原和海原土壤干层分布深度达到10m所需时间依次为6、6和4年,此后苜蓿草地降水渗深以下土层长期维持较为稳定的干燥化状态;苜蓿草地水分持续利用的合理年限为半湿润区8～10年,半干旱区6～8年,半干旱偏旱区4～6年.

Intraseasonal variability of Indian Ocean sea surface temperature during boreal winter: Madden-Julian Oscillation versus submonthly forcing and processes

[ 1] Intraseasonal variability of Indian Ocean sea surface temperature (SST) during boreal winter is investigated by analyzing available data and a suite of solutions to an ocean general circulation model for 1998 - 2004. This period covers the QuikSCAT and Tropical Rainfall Measuring Mission (TRMM) observations. Impacts of the 30 - 90 day and 10 - 30 day atmospheric intraseasonal oscillations (ISOs) are examined separately, with the former dominated by the Madden-Julian Oscillation (MJO) and the latter dominated by convectively coupled Rossby and Kelvin waves. The maximum variation of intraseasonal SST occurs at 10 degrees S - 2 degrees S in the wintertime Intertropical Convergence Zone (ITCZ), where the mixed layer is thin and intraseasonal wind speed reaches its maximum. The observed maximum warming ( cooling) averaged over ( 60 degrees E - 85 degrees E, 10 degrees S - 3 degrees S) is 1.13 degrees C ( - 0.97 degrees C) for the period of interest, with a standard deviation of 0.39 degrees C in winter. This SST change is forced predominantly by the MJO. While the MJO causes a basin-wide cooling ( warming) in the ITCZ region, submonthly ISOs cause a more complex SST structure that propagates southwestward in the western-central basin and southeastward in the eastern ocean. On both the MJO and submonthly timescales, winds are the deterministic factor for the SST variability. Short-wave radiation generally plays a secondary role, and effects of precipitation are negligible. The dominant role of winds results roughly equally from wind speed and stress forcing. Wind speed affects SST by altering turbulent heat fluxes and entrainment cooling. Wind stress affects SST via several local and remote oceanic processes.

太阳辐射对热层和电离层变化性的影响

Solar ultraviolet (UV) radiation at wavelengths less than 400 nm is an important source of energy for aeronomic processes throughout the solar system. Solar UV photons are absorbed in planetary atmospheres, as well as throughout the heliosphere, via photodissociation of molecules, photoionization of molecules and atoms, and photoexcitation toexcitation including resonance scattering. In this paper, the solar irradiances data measured by TIMED SEE, as well as the solar proxies such as F10.7 and Mg II, thermosphere neutral density of CHAMP measurements and topside ionospheric plasmas densities from DMSP, are used to analyze solar irradiance effects on the variabilities of the thermosphere and the ionosphere. First, thermosphere densities near 410 km altitude are analyzed for solar irradiance variability effects during the period 2002-2004. Correlations between the densities and the solar irradiances for different spectral lines and wavelength ranges reveal significantly different characteristics. The density correlates remarkably well with all the selected solar irradiances except the lower chromospheric O I (130.4 nm) emission. Among the chosen solar proxies, the Mg II core-to-wing ratio index, EUV (30-120 nm) and F10.7 show the highest correlations with the density for short-term (< ~27 days) variations. For both long- (> ~27 days) and short-term variations, linear correlation coefficients exhibit a decreasing trend from low latitudes towards high latitudes. The density variability can be effectively modeled (capturing 71% of the variance) using multiple solar irradiance indices, including F10.7, SEUV (the EUV 30-120 nm index), and SFUV (the FUV 120-193 nm index), in which a lag time of 1 day was used for both F10.7 and SEUV, and 5 days for SFUV. In our regression formulation SEUV has the largest contribution to the density variation (40%), with the F10.7 having the next largest contribution (32%) and SFUV accounting for the rest (28%). Furthermore, a pronounced period of about 27.2 days (mean period of the Sun's rotation) is present in both density and solar irradiance data of 2003 and 2004, and a pronounced period of about 54.4 days (doubled period of the solar rotation) is also revealed in 2004. However, soft X-ray and FUV irradiances did not present a pronounced 54.4 day period in 2004, in spite of their high correlation with the densities. The Ap index also shows 54-day periodicities in 2004, and magnetic activity, together with solar irradiance, affects the 54-day variation in density significantly. In addition, NRLMSISE00, DTM-2000 and JB2006 model predictions are compared with density measurements from CHAMP to assess their accuracy, and the results show that these models underestimate the response of the thermosphere to variations induced by solar rotation. Next, the equatorial topside ionospheric plasmas densities Ni are analyzed for solar irradiance variability effects during the period 2002-2005. Linear correlations between Ni and the solar irradiances for different wavelength ranges reveal significantly different characteristics. XUV (0-35 nm) and EUV (115-130 nm) show higher correlation with Ni for the long-term variations, whereas EUV (35-115 nm) show higher correlation for the short-term variations. Moreover, partial correlation analysis shows that the long-term variations of Ni are affected by both XUV (0-35 nm) and EUV (35-115 nm), whereas XUV (0-35 nm) play a more important role; the short-term variations of Ni are mostly affected by EUV (35-115 nm). Furthermore, a pronounced period of about 27 days is present in both Ni and solar irradiance data of 2003 and 2004, and a pronounced period of about 54 days is also revealed in 2004. Finally, prompted by previous studies that have suggested solar EUV radiation as a means of driving the semiannual variation, we investigate the intra-annual variation in thermosphere neutral density near 400 km during 2002-2005. The intra-annual variation, commonly referred to as the ‘semiannual variation’, is characterized by significant latitude structure, hemispheric asymmetries, and inter-annual variability. The magnitude of the maximum yearly difference, from the yearly minimum to the yearly maximum, varies by as much as 60% from year to year, and the phases of the minima and maxima also change by 20-40 days from year to year. Each annual harmonic of the intra-annual variation, namely, annual, semiannual, ter-annual and quatra-annual, exhibits a decreasing trend from 2002 through 2005 that is correlated with the decline in solar activity. In addition, some variations in these harmonics are correlated with geomagnetic activity, as represented by the daily mean value of Kp. Recent empirical models of the thermosphere are found to be deficient in capturing most of the latitude dependencies discovered in our data. In addition, the solar flux and geomagnetic activity proxies that we have employed do not capture some latitude and inter-annual variations detected in our data. It is possible that these variations are partly due to other effects, such as seasonal-latitudinal variations in turbopause altitude (and hence O/N2 composition) and ionosphere coupling processes that remain to be discovered in the context of influencing the intra-annual variations depicted here. Our results provide a new dataset to challenge and validate thermosphere-ionosphere general circulation models that seek to delineate the thermosphere intra-annual variation and to understand the various competing mechanisms that may contribute to its existence and variability. We furthermore suggest that the term “intra-annual” variation be adopted to describe the variability in thermosphere and ionosphere parameters that is well-captured through a superposition of annual, semiannual, ter-annual, and quatra-annual harmonic terms, and that “semiannual’ be used strictly in reference to a pure 6-monthly sinusoidal variation. Moreover, we propose the term “intra-seasonal” to refer to those shorter-term variations that arise as residuals from the above Fourier representation.

Victims of infanticide and conspecific bite wounding in a female-dominant primate: a long-term study.

The aggression animals receive from conspecifics varies between individuals across their lifetime. As poignantly evidenced by infanticide, for example, aggression can have dramatic fitness consequences. Nevertheless, we understand little about the sources of variation in received aggression, particularly in females. Using a female-dominant species renowned for aggressivity in both sexes, we tested for potential social, demographic, and genetic patterns in the frequency with which animals were wounded by conspecifics. Our study included 243 captive, ring-tailed lemurs (Lemur catta), followed from infancy to adulthood over a 35-year time span. We extracted injury, social, and life-history information from colony records and calculated neutral heterozygosity for a subset of animals, as an estimate of genetic diversity. Focusing on victims rather than aggressors, we used General Linear Models to explain bite-wound patterns at different life stages. In infancy, maternal age best predicted wounds received, as infants born to young mothers were the most frequent infanticide victims. In adulthood, sex best predicted wounds received, as males were three times more likely than females to be seriously injured. No relation emerged between wounds received and the other variables studied. Beyond the generally expected costs of adult male intrasexual aggression, we suggest possible additive costs associated with female-dominant societies - those suffered by young mothers engaged in aggressive disputes and those suffered by adult males aggressively targeted by both sexes. We propose that infanticide in lemurs may be a costly by-product of aggressively mediated, female social dominance. Accordingly, the benefits of female behavioral 'masculinization' accrued to females through priority of access to resources, may be partially offset by early costs in reproductive success. Understanding the factors that influence lifetime patterns of conspecific wounding is critical to evaluating the fitness costs associated with social living; however, these costs may vary substantially between societies.

A shared care model pilot for palliative home care in a rural area:Impact on symptoms, distress, and place of death

Context: Shared care models integrating family physician services with interdisciplinary palliative care specialist teams are critical to improve access to quality palliative home care and address multiple domains of end-of-life issues and needs. Objectives: To examine the impact of a shared care pilot program on the primary outcomes of symptom severity and emotional distress (patient and family separately) over time and, secondarily, the concordance between patient preferences and place of death. Methods: An inception cohort of patients (n = 95) with advanced, progressive disease, expected to die within six months, were recruited from three rural family physician group practices (21 physicians) and followed prospectively until death or pilot end. Serial measurement of symptoms, emotional distress (patient and family), and preferences for place of death was performed, with analysis of changes in distress outcomes assessed using t-tests and general linear models. Results: Symptoms trended toward improvement, with a significant reduction in anxiety from baseline to 14 days noted. Symptom and emotional distress were maintained below high severity (7-10), and a high rate of home death compared with population norms was observed. Conclusion: Future controlled studies are needed to examine outcomes for shared care models with comparison groups. Shared care models build on family physician capacity and as such are promising in the development of palliative home care programs to improve access to quality palliative home care and foster health system integration. © 2011 U.S. Cancer Pain Relief Committee. Published by Elsevier Inc. All rights reserved.

A comparative study of Type II-P and II-L supernova rise times as exemplified by the case of LSQ13cuw

We report on our findings based on the analysis of observations of the Type II-L supernova LSQ13cuw within the framework of currently accepted physical predictions of core-collapse supernova explosions. LSQ13cuw was discovered within a day of explosion, hitherto unprecedented for Type II-L supernovae. This motivated a comparative study of Type II-P and II-L supernovae with relatively well-constrained explosion epochs and rise times to maximum (optical) light. From our sample of twenty such events, we find evidence of a positive correlation between the duration of the rise and the peak brightness. On average, SNe II-L tend to have brighter peak magnitudes and longer rise times than SNe II-P. However, this difference is clearest only at the extreme ends of the rise time versus peak brightness relation. Using two different analytical models, we performed a parameter study to investigate the physical parameters that control the rise time behaviour. In general, the models qualitatively reproduce aspects of the observed trends. We find that the brightness of the optical peak increases for larger progenitor radii and explosion energies, and decreases for larger masses. The dependence of the rise time on mass and explosion energy is smaller than the dependence on the progenitor radius. We find no evidence that the progenitors of SNe II-L have significantly smaller radii than those of SNe II-P.