School of Earth Sciences - Research Publications

Permanent URI for this collection

Search Results

Now showing 1 - 10 of 13
  • Item
    No Preview Available
    Simultaneous nested modeling from the synoptic scale to the LES scale for wind energy applications
    Liu, Y ; Warner, T ; Liu, Y ; Vincent, C ; Wu, W ; Mahoney, B ; Swerdlin, S ; Parks, K ; Boehnert, J (ELSEVIER, 2011-04)
  • Item
    Thumbnail Image
    Non-stationarity in daily and sub-daily intense rainfall – Part 2: regional assessment for sites in south-east Australia
    Jakob, D. ; Karoly, D. J. ; Seed, A. (Copernicus Publications on behalf of the European Geoscience Union, 2011)
    Using data for a common period (1976–2005) for a set of 31 sites located in south-east Australia, variations in frequency and magnitude of intense rainfall events across durations from 6 min to 72 h were assessed. This study was driven by a need to clarify how variations in climate might affect intense rainfall and the potential for flooding. Sub-daily durations are of particular interest for urban applications. Worldwide, few such observation-based studies exist, which is mainly due to limitations in data. Analysis of seasonality in frequency and magnitude of events revealed considerable variation across the set of sites, implying different dominating rainfall-producing mechanisms and/or interactions with local topography. Both these factors are relevant when assessing the potential effects of climate variations on intense rainfall events. The set of sites was therefore split into groups ("north cluster" and "south cluster") according to the characteristics of intense rainfall events. There is a strong polarisation in the nature of changes found for the north cluster and south cluster. While sites in the north cluster typically exhibit decrease in frequency of events, particularly in autumn and at durations of 1 h and longer; sites in the south cluster experience an increase in frequency of events, particularly for summer and sub-hourly durations. Non-stationarity found in historical records has the potential to significantly affect design rainfall estimates. An assessment of quantile estimates derived using a standard regionalisation technique and periods representative of record lengths available for practical applications show that such estimates may not be representative of long-term conditions, so alternative approaches need to be considered, particularly where short records are concerned. Additional rainfall information, in particular radar data, could be used for an in-depth spatial analysis of intense rainfall events.
  • Item
    Thumbnail Image
    Non-stationarity in daily and sub-daily intense rainfall – Part 1: Sydney, Australia
    Jakob, D. ; Karoly, D. J. ; Seed, A. (Copernicus Publications on behalf of the European Geoscience Union, 2011)
    This study was driven by a need to clarify how variations in climate might affect intense rainfall and the potential for flooding. Sub-daily durations are of particular interest for urban applications. Worldwide, few such observation-based studies exist, which is mainly due to limitations in data. While there are still large discrepancies between precipitation data sets from observations and models, both show that there is a tendency for moist regions to become wetter and for dry regions to become drier. However, changes in extreme conditions may show the opposite sign to those in average conditions. Where changes in observed intense precipitation have been studied, this has typically been for daily durations or longer. The purpose of this two-part study is to examine daily and sub-daily rainfall extremes for evidence of non-stationarity. Here the problem was addressed by supplementing one long record (Part 1) by a set of shorter records for a 30-yr concurrent period (Part 2). Variations in frequency and magnitude of rainfall extremes across durations from 6 min to 72 h were assessed using data from sites in the south-east of Australia. For the analyses presented in this paper, a peaks-over-threshold approach was chosen since it allows investigating changes in frequency as well as magnitude. Non-parametric approaches were used to assess changes in frequency, magnitude, and quantile estimates as well as the statistical significance of changes for one station (Sydney Observatory Hill) for the period 1921 to 2005. Deviations from the long-term average vary with season, duration, and threshold. The effects of climate variations are most readily detected for the highest thresholds. Deviations from the long-term average tend to be larger for frequencies than for magnitudes, and changes in frequency and magnitude may have opposite signs. Investigations presented in this paper show that variations in frequency and magnitude of events at daily durations are a poor indicator of changes at sub-daily durations. Studies like the one presented here should be undertaken for other regions to allow the identification of regions with significant increase/decrease in intense rainfall, whether there are common features with regards to duration and season exhibiting most significant changes (which in turn could lead to establishing a theoretical framework), and assist in validation of projections of rainfall extremes.
  • Item
    Thumbnail Image
    The Greenland Ice Sheet Response to Transient Climate Change
    Ren, D ; Fu, R ; Leslie, LM ; Chen, J ; Wilson, CR ; Karoly, DJ (AMER METEOROLOGICAL SOC, 2011-07)
    Abstract This study applies a multiphase, multiple-rheology, scalable, and extensible geofluid model to the Greenland Ice Sheet (GrIS). The model is driven by monthly atmospheric forcing from global climate model simulations. Novel features of the model, referred to as the scalable and extensible geofluid modeling system (SEGMENT-Ice), include using the full Navier–Stokes equations to account for nonlocal dynamic balance and its influence on ice flow, and a granular sliding layer between the bottom ice layer and the lithosphere layer to provide a mechanism for possible large-scale surges in a warmer future climate (granular basal layer is for certain specific regions, though). Monthly climate of SEGMENT-Ice allows an investigation of detailed features such as seasonal melt area extent (SME) over Greenland. The model reproduced reasonably well the annual maximum SME and total ice mass lost rate when compared observations from the Special Sensing Microwave Imager (SSM/I) and Gravity Recovery and Climate Experiment (GRACE) over the past few decades. The SEGMENT-Ice simulations are driven by projections from two relatively high-resolution climate models, the NCAR Community Climate System Model, version 3 (CCSM3) and the Model for Interdisciplinary Research on Climate 3.2, high-resolution version [MIROC3.2(hires)], under a realistic twenty-first-century greenhouse gas emission scenario. They suggest that the surface flow would be enhanced over the entire GrIS owing to a reduction of ice viscosity as the temperature increases, despite the small change in the ice surface topography over the interior of Greenland. With increased surface flow speed, strain heating induces more rapid heating in the ice at levels deeper than due to diffusion alone. Basal sliding, especially for granular sediments, provides an efficient mechanism for fast-glacier acceleration and enhanced mass loss. This mechanism, absent from other models, provides a rapid dynamic response to climate change. Net mass loss estimates from the new model should reach ~220 km3 yr−1 by 2100, significantly higher than estimates by the Intergovernmental Panel on Climate Change (IPCC) Assessment Report 4 (AR4) of ~50–100 km3 yr−1. By 2100, the perennial frozen surface area decreases up to ~60%, to ~7 × 105 km2, indicating a massive expansion of the ablation zone. Ice mass change patterns, particularly along the periphery, are very similar between the two climate models.
  • Item
    Thumbnail Image
    A regional modeling study of climate change impacts on warm-season precipitation in the Central United States
    Bukovsky, Melissa S. ; Karoly, David J. (American Meteorological Society, 2011)
    In this study, the Weather Research and Forecasting (WRF) model is employed as a nested regional climate model to dynamically downscale output from the National Center for Atmospheric Research’s (NCAR’s) Community Climate SystemModel (CCSM) version 3 and the National Centers for Environmental Prediction (NCEP)–NCARglobal reanalysis (NNRP). The latter is used for verification of late-twentieth-century climate simulations from the WRF.
  • Item
    Thumbnail Image
    A three-dimensional synthesis inversion of the molecular hydrogen cycle: sources and sinks budget and implications for the soil uptake
    Bousquet, P. ; Yver, C. ; Pison, I. ; Li, Y. S. ; Fortems, A. ; Hauglustaine, D. ; Szopa, S. ; Rayner, P. J. ; Novelli, P. ; Langenfelds, R. ; STEELE, P. ; Ramonet, M. ; Schmidt, M. ; Foster, P. ; Morfopoulos, C. ; Ciais, P. (American Geophysical Union, 2011)
    Our understanding of the global budget of atmospheric hydrogen (H2) contains large uncertainties. An atmospheric Bayesian inversion of H2 sources and sinks is presented for the period 1991-2004, based on a two networks of flask measurement stations. The types of fluxes and the spatial scales potentially resolvable by the inversion are first estimated from an analysis of the correlations of errors between the different processes and regions emitting or absorbing H2. Then, the estimated budget of H2 and its uncertainties is presented and discussed, for five groups of fluxes and three groups of large regions, in terms of mean fluxes, seasonal and interannual variations, and long-term trends. One main focus of the study is the improvement of the estimate of H2 soil uptake, which is the largest sink of H2. Various sensitivity tests are performed defining an ensemble of more than 20 inversions. We show that inferring a robust estimate of the H2 soil uptake requires to prescribe the prior magnitude of some other sources and sinks with a small uncertainty. Doing so an estimate of the H2 soil uptake of -62 ± 3 Tg y−1 is inferred for the period 1991-2004 (the uncertainty is the residual error after inversion). The inferred soil H2 sink presents a negative long-term trend that is qualitatively consistent with a bottom-up process-based model.
  • Item
    Thumbnail Image
    A severe thunderstorm climatology for Australia and associated thunderstorm environments
    Allen, John Terrence ; Karoly, David J. ; Mills, Graham A. (Australian Meteorological and Oceanographic Society, 2011)
    Severe thunderstorms can present a significant threat to life and property in Australia. A unique and broad database of severe thunderstorm reports has been constructed for the Australian region for 2003–2010 from observer reports of hailstones, winds in excess of 90 km h–1 and, less frequently, tornadoes. Based on this database, a climatology of atmospheric environments associated with the occurrence of severe thunderstorms in Australia was developed using pseudo-proximity soundings from the MesoLAPS numerical weather prediction model simulations. Observed soundings have been used to verify derived soundings from MesoLAPS simulations, with a reasonable performance over much of the continent. Proximity rawinsonde soundings from the MesoLAPS simulations were identified for each of the severe thunderstorm reports to develop the climatology of environments. This climatology was then used to derive discriminants between environments with an increased likelihood of severe thunderstorm occurrence and other thunderstorm environments. This appears to be the best way to produce a long-term climatology of severe thunderstorm environment occurrence in a sparsely populated continent without considering the complex problem of initiation.
  • Item
    Thumbnail Image
    Atmospheric CO2 inversion validation using vertical profile measurements: analysis of four independent inversion models
    Peylin, P. ; Maksyutov, S. ; Marshall, J. ; Rödenbeck, C. ; Langenfelds, R.L. ; Steele, L.P. ; Francey, R.J. ; Tans P. ; Sweeney C. ; Pickett-Heaps, C. A. ; Rayner, P. J. ; Law, R. M. ; Ciais, P. ; Patra, P. K. ; Bousquet, P. (American Geophysical Union, 2011)
    We present the results of a validation of atmospheric inversions of CO2 fluxes using four transport models. Each inversion uses data primarily from surface stations, combined with an atmospheric transport model, to estimate surface fluxes. The validation (or model evaluation) consists of running these optimized fluxes through the forward model and comparing the simulated concentrations with airborne concentration measurements. We focus on profiles from Cape Grim, Tasmania, and Carr, Colorado, while using other profile sites to test the generality of the comparison. Fits to the profiles are generally worse than to the surface data from the inversions and worse than the expected model-data mismatch. Thus inversion estimates are generally not consistent with the profile measurements. The TM3 model does better by some measures than the other three models. Models perform better over Tasmania than Colorado, and other profile sites bear out a general improvement from north to south and from continental to marine locations. There are also errors in the interannual variability of the fit, consistent in time and common across models. This suggests real variations in sources visible to the profile but not the surface measurements.
  • Item
    Thumbnail Image
    Optimal representation of source-sink fluxes for mesoscale carbon dioxide inversion with synthetic data
    Wu, Lin ; Bocquet, Marc ; Lauvaux, Thomas ; Chevallier, Frédéric ; RAYNER, PETER ; Davis, Kenneth (American Geophysical Union, 2011)
    The inversion of CO2 surface fluxes from atmospheric concentration measurements involves discretizing the flux domain in time and space. The resolution choice is usually guided by technical considerations despite its impact on the solution to the inversion problem. In our previous studies, a Bayesian formalism has recently been introduced to describe the discretization of the parameter space over a large dictionary of adaptive multiscale grids. In this paper, we exploit this new framework to construct optimal space-time representations of carbon fluxes for mesoscale inversions. Inversions are performed using synthetic continuous hourly CO2 concentration data in the context of the Ring 2 experiment in support of the North American Carbon Program Mid Continent Intensive (MCI). Compared with the regular grid at finest scale, optimal representations can have similar inversion performance with far fewer grid cells. These optimal representations are obtained by maximizing the number of degrees of freedom for the signal (DFS) that measures the information gain from observations to resolve the unknown fluxes. Consequently information from observations can be better propagated within the domain through these optimal representations. For the Ring 2 network of eight towers, in most cases, the DFS value is relatively small compared to the number of observations d (DFS/d < 20%). In this multiscale setting, scale-dependent aggregation errors are identified and explicitly formulated for more reliable inversions. It is recommended that the aggregation errors should be taken into account, especially when the correlations in the errors of a priori fluxes are physically unrealistic. The optimal multiscale grids allow to adaptively mitigate the aggregation errors.
  • Item
    Thumbnail Image
    A European summertime CO2 biogenic flux inversion at mesoscale from continuous in situ mixing ratio measurements
    Broquet, Grégoire ; Chevallier, Frédéric ; RAYNER, PETER ; Aulagnier, Céline ; Pison, Isabelle ; Ramonet, Michel ; Schmidt, Martina ; Vermeulen, Alex T. ; Ciais, Philippe (American Geophysical Union, 2011)
    A regional variational inverse modeling system for the estimation of European biogenic CO2 fluxes is presented. This system is based on a 50 km horizontal resolution configuration of a mesoscale atmospheric transport model and on the adjoint of its tracer transport code. It exploits hourly CO2 in situ data from 15 CarboEurope-Integrated Project stations. Particular attention in the inversion setup is paid to characterizing the transport model error and to selecting the observations to be assimilated as a function of this error. Comparisons between simulations and data of CO2 and 222Rn concentrations indicate that the model errors should have a standard deviation which is less than 7 ppm when simulating the hourly variability of CO2 at low altitude during the afternoon and evening or at high altitude at night. Synthetic data are used to estimate the uncertainty reduction for the fluxes using this inverse modeling system. The improvement brought by the inversion to the prior estimate of the fluxes for both the mean diurnal cycle and the monthly to synoptic variability in the fluxes and associated mixing ratios are checked against independent atmospheric data and eddy covariance flux measurements. Inverse modeling is conducted for summers 2002 - 2007 which should reduce the uncertainty in the biogenic fluxes by ∼60% during this period. The trend in the mean flux corrections between June and September is to increase the uptake of CO2 by ∼12 gCm−2. Corrections at higher resolution are also diagnosed that reveal some limitations of the underlying prior model of the terrestrial biosphere.