The role of local heating in the 2015 Indian Heat Wave

India faced a major heat wave during the summer of 2015. Temperature anomalies peaked in the dry period before the onset of the summer monsoon, suggesting that local land-atmosphere feedbacks involving desiccated soils and vegetation might have played a role in driving the heat extreme. Upon examination of in situ data, reanalysis, satellite observations, and land surface models, we find that the heat wave included two distinct peaks: one in late May, and a second in early June. During the first peak we find that clear skies led to a positive net radiation anomaly at the surface, but there is no significant sensible heat flux anomaly within the core of the heat wave affected region. By the time of the second peak, however, soil moisture had dropped to anomalously low levels in the core heat wave region, net surface radiation was anomalously high, and a significant positive sensible heat flux anomaly developed. This led to a substantial local forcing on air temperature that contributed to the intensity of the event. The analysis indicates that the highly agricultural landscape of North and Central India can reinforce heat extremes under dry conditions.

Seasonal range test run with Global Eta Framework

Abstract. Global Eta Framework (GEF) is a global atmospheric model developed in general curvilinear coordinates and capable of running on arbitrary rectangular quasi-uniform spherical grids, using stepwise (Eta) representation of the terrain. In this study, the model is run on a cubed-sphere grid topology, in a version with uniform Jacobians (UJ), which provides equal-area grid cells, and a smooth transition of coordinate lines across the edges of the cubed-sphere. Within a project at the Brazilian Center for Weather Forecasts and Climate Studies (CPTEC), a nonhydrostatic version of this model is under development and will be applied for seasonal prediction studies. This note describes preliminary tests with the GEF on the UJ cubed-sphere in which model performance is evaluated in seasonal simulations at a horizontal resolution of approximately 25 km, running in the hydrostatic mode. Comparison of these simulations with the ERA-Interim reanalyses shows that the 850 hPa temperature is underestimated, while precipitation pattern is mostly underestimated in tropical continental regions and overestimated in tropical oceanic regions. Nevertheless, the model is still able to well capture the main seasonal climate characteristics. These results will be used as a control run in further tests with the nonhydrostatic version of the model.

Tradeoff between cost and accuracy in large-scale surface water dynamic modeling

Recent efforts have led to the development of the local inertia formulation (INER) for an accurate but still cost-efficient representation of surface water dynamics, compared to the widely used kinematic wave equation (KINE). In this study, both formulations are evaluated over the Amazon basin in terms of computational costs and accuracy in simulating streamflows and water levels through synthetic experiments and comparisons against ground-based observations. Varying time steps are considered as part of the evaluation and INER at 60-second time step is adopted as the reference for synthetic experiments. Five hybrid (HYBR) realizations are performed based on maps representing the spatial distribution of the two formulations that physically represent river reach flow dynamics within the domain. Maps have fractions of KINE varying from 35.6% to 82.8%. KINE runs show clear deterioration along the Amazon river and main tributaries, with maximum RMSE values for streamflow and water level reaching 7827m3.s-1 and 1379cm near the basin's outlet. However, KINE is at least 25% more efficient than INER with low model sensitivity to longer time steps. A significant improvement is achieved with HYBR, resulting in maximum RMSE values of 3.9-292m3.s-1 for streamflows and 1.1-28.5cm for water levels, and cost reduction of 6-16%, depending on the map used. Optimal results using HYBR are obtained when the local inertia formulation is used in about one third of the Amazon basin, reducing computational costs in simulations while preserving accuracy. However, that threshold may vary when applied to different regions, according to their hydrodynamics and geomorphological characteristics.

Effects of life-history requirements on the distribution of a threatened reptile

Survival and reproduction are the two primary life-history traits essential for species' persistence; however, the environmental conditions that support each of these traits may not be the same. Despite this, reproductive requirements are seldom considered when estimating species' potential distributions. We sought to examine potentially limiting environmental factors influencing the distribution of an oviparous reptile of conservation concern with respect to the species' survival and reproduction and to assess the implications of the species' predicted climatic constraints on current conservation practices. We used ecological niche modeling to predict the probability of environmental suitability for the alligator snapping turtle (Macrochelys temminckii). We built an annual climate model to examine survival and a nesting climate model to examine reproduction. We combined incubation temperature requirements, products of modeled soil temperature data, and our estimated distributions to determine whether embryonic development constrained the northern distribution of the species. Low annual precipitation constrained the western distribution of alligator snapping turtles, whereas the northern distribution was constrained by thermal requirements during embryonic development. Only a portion of the geographic range predicted to have a high probability of suitability for alligator snapping turtle survival was estimated to be capable of supporting successful embryonic development. Historic occurrence records suggest adult alligator snapping turtles can survive in regions with colder climes than those associated with consistent and successful production of offspring. Estimated egg-incubation requirements indicated that current reintroductions at the northern edge of the species' range are within reproductively viable environmental conditions. Our results highlight the importance of considering survival and reproduction when estimating species' ecological niches, implicating conservation plans, and benefits of incorporating physiological data when evaluating species' distributions.

Utility of Satellite Remote Sensing for Land-Atmosphere Coupling and Drought Metrics

Feedbacks between the land and the atmosphere can play an important role in the water cycle and a number of studies have quantified Land-Atmosphere (L-A) interactions and feedbacks through observations and prediction models. Due to the complex nature of L-A interactions, the observed variables are not always available at the needed temporal and spatial scales. This work derives the Coupling Drought Index (CDI) solely from satellite data and evaluates the input variables and the resultant CDI against in-situ data and reanalysis products. NASA's AQUA satellite and retrievals of soil moisture and lower tropospheric temperature and humidity properties are used as input. Overall, the AQUA-based CDI and its inputs perform well at a point, spatially, and in time (trends) compared to in-situ and reanalysis products. In addition, this work represents the first time that in-situ observations were utilized for the coupling classification and CDI. The combination of in-situ and satellite remote sensing CDI is unique and provides an observational tool for evaluating models at local and large scales. Overall, results indicate that there is sufficient information in the signal from simultaneous measurements of the land and atmosphere from satellite remote sensing to provide useful information for applications of drought monitoring and coupling metrics.

The role of low-level terrain-induced jets in rainfall variability in Tigris-Euphrates Headwaters

Rainfall variability in the Tigris-Euphrates Headwaters is a result of interaction between topography and meteorological features at a range of spatial scales. Here, we have implemented the Weather Research and Forecasting (WRF) model, driven by NCEP/DOE R2, to better understand these interactions. Simulations were performed over a domain covering most of the Middle-East. The extended simulation period (1983-2013) enables us to study seasonality, interannual variability, spatial variability and extreme events of rainfall. Results showed that the annual cycle of precipitation produced by WRF agrees much more closely with observations than does R2. This was particularly evident during the transition months of April and October, which were further examined to study the underlying physical mechanisms. In both months, WRF improves representation of interannual variability relative to R2, with a substantially larger benefit in April. This improvement results primarily from WRF's ability to resolve two low-level terrain-induced flows in the region that are either absent or weak in NCEP/DOE: one parallel to western edge of the Zagros Mountains, and one along the East Turkish Highlands. The first shows a complete reversal in its direction during wet and dry days: when flowing southeasterly it transports moisture from the Persian Gulf to the region, and when flowing northwesterly it blocks moisture and transports it away from the region. The second is more directly related to synoptic-scale systems and carries moist, warm air from the Mediterranean and Red Seas toward the region. The combined contribution of these flows explains about 50% of interannual variability in both WRF and observations for April and October precipitation.

A land data assimilation system for sub-Saharan Africa food and water security applications

Seasonal agricultural drought monitoring systems, which rely on satellite remote sensing and land surface models (LSMs), are important for disaster risk reduction and famine early warning. These systems require the best available weather inputs, as well as a long-term historical record to contextualize current observations. This article introduces the Famine Early Warning Systems Network (FEWS NET) Land Data Assimilation System (FLDAS), a custom instance of the NASA Land Information System (LIS) framework. The FLDAS is routinely used to produce multi-model and multi-forcing estimates of hydro-climate states and fluxes over semi-arid, food insecure regions of Africa. These modeled data and derived products, like soil moisture percentiles and water availability, were designed and are currently used to complement FEWS NET's operational remotely sensed rainfall, evapotranspiration, and vegetation observations. The 30+ years of monthly outputs from the FLDAS simulations are publicly available from the NASA Goddard Earth Science Data and Information Services Center (GES DISC) and recommended for use in hydroclimate studies, early warning applications, and by agro-meteorological scientists in Eastern, Southern, and Western Africa.

Comparing potential recharge estimates from three Land Surface Models across the Western US

Groundwater is a major source of water in the western US. However, there are limited recharge estimates available in this region due to the complexity of recharge processes and the challenge of direct observations. Land surface Models (LSMs) could be a valuable tool for estimating current recharge and projecting changes due to future climate change. In this study, simulations of three LSMs (Noah, Mosaic and VIC) obtained from the North American Land Data Assimilation System (NLDAS-2) are used to estimate potential recharge in the western US. Modeled recharge was compared with published recharge estimates for several aquifers in the region. Annual recharge to precipitation ratios across the study basins varied from 0.01-15% for Mosaic, 3.2-42% for Noah, and 6.7-31.8% for VIC simulations. Mosaic consistently underestimates recharge across all basins. Noah captures recharge reasonably well in wetter basins, but overestimates it in drier basins. VIC slightly overestimates recharge in drier basins and slightly underestimates it for wetter basins. While the average annual recharge values vary among the models, the models were consistent in identifying high and low recharge areas in the region. Models agree in seasonality of recharge occurring dominantly during the spring across the region. Overall, our results highlight that LSMs have the potential to capture the spatial and temporal patterns as well as seasonality of recharge at large scales. Therefore, LSMs (specifically VIC and Noah) can be used as a tool for estimating future recharge rates in data limited regions.

Impacts of uncertainties in European gridded precipitation observations on regional climate analysis

Gridded precipitation data sets are frequently used to evaluate climate models or to remove model output biases. Although precipitation data are error prone due to the high spatio-temporal variability of precipitation and due to considerable measurement errors, relatively few attempts have been made to account for observational uncertainty in model evaluation or in bias correction studies. In this study, we compare three types of European daily data sets featuring two Pan-European data sets and a set that combines eight very high-resolution station-based regional data sets. Furthermore, we investigate seven widely used, larger scale global data sets. Our results demonstrate that the differences between these data sets have the same magnitude as precipitation errors found in regional climate models. Therefore, including observational uncertainties is essential for climate studies, climate model evaluation, and statistical post-processing. Following our results, we suggest the following guidelines for regional precipitation assessments. (1) Include multiple observational data sets from different sources (e.g. station, satellite, reanalysis based) to estimate observational uncertainties. (2) Use data sets with high station densities to minimize the effect of precipitation undersampling (may induce about 60% error in data sparse regions). The information content of a gridded data set is mainly related to its underlying station density and not to its grid spacing. (3) Consider undercatch errors of up to 80% in high latitudes and mountainous regions. (4) Analyses of small-scale features and extremes are especially uncertain in gridded data sets. For higher confidence, use climate-mean and larger scale statistics. In conclusion, neglecting observational uncertainties potentially misguides climate model development and can severely affect the results of climate change impact assessments.

Does temperature nudging overwhelm aerosol radiative effects in regional integrated climate models?

Nudging (data assimilation) is used in many regional integrated meteorology-air quality models to reduce biases in simulated climatology. However, in such modeling systems, temperature changes due to nudging could compete with temperature changes induced by radiatively active and hygroscopic short-lived tracers leading to two interesting dilemmas: when nudging is continuously applied, what are the relative sizes of these two radiative forces at regional and local scales? How do these two forces present in the free atmosphere differ from those present at the surface? This work studies these two issues by converting temperature changes due to nudging into pseudo radiative effects (PRE) at the surface (PRE_sfc), in troposphere (PRE_atm), and at the top of atmosphere (PRE_toa), and comparing PRE with the reported aerosol radiative effects (ARE). Results show that the domain-averaged PRE_sfc is smaller than ARE_sfc estimated in previous studies and this work, but could be significantly larger than ARE_sfc at local scales. PRE_atm is also much smaller than ARE_atm. These results indicate that appropriate nudging methodology could be applied to the integrated models to study aerosol radiative effects at continental/regional scales, but it should be treated with caution for local scale applications.

Assimilation of SMOS Retrievals in the Land Information System

The Soil Moisture and Ocean Salinity (SMOS) satellite provides retrievals of soil moisture in the upper 5 cm with a 30-50 km resolution and a mission accuracy requirement of 0.04 cm3 cm-3. These observations can be used to improve land surface model soil moisture states through data assimilation. In this paper, SMOS soil moisture retrievals are assimilated into the Noah land surface model via an Ensemble Kalman Filter within the NASA Land Information System. Bias correction is implemented using Cumulative Distribution Function (CDF) matching, with points aggregated by either land cover or soil type to reduce sampling error in generating the CDFs. An experiment was run for the warm season of 2011 to test SMOS data assimilation and to compare assimilation methods. Verification of soil moisture analyses in the 0-10 cm upper layer and root zone (0-1 m) was conducted using in situ measurements from several observing networks in the central and southeastern United States. This experiment showed that SMOS data assimilation significantly increased the anomaly correlation of Noah soil moisture with station measurements from 0.45 to 0.57 in the 0-10 cm layer. Time series at specific stations demonstrate the ability of SMOS DA to increase the dynamic range of soil moisture in a manner consistent with station measurements. Among the bias correction methods, the correction based on soil type performed best at bias reduction but also reduced correlations. The vegetation-based correction did not produce any significant differences compared to using a simple uniform correction curve.

Comparisons of Modeled and Observed Reflectivities and Fall Speeds for Snowfall of Varied Riming Degrees during Winter Storms on Long Island, New York

Derived radar reflectivities and fall speeds for four Weather Research and Forecasting (WRF) Model bulk microphysical parameterizations (BMPs) run at 1.33-km grid spacing are compared with ground-based, vertically pointing Ku-band radar, scanning S-band radar, and in situ measurements at Stony Brook, New York. Simulations were partitioned into periods of observed riming degree as determined manually using a stereo microscope and camera during nine winter storms. Simulations were examined to determine whether the selected BMPs captured the effects of varying riming intensities, provided a reasonable match to the vertical structure of radar reflectivity or fall speed, and whether they produced reasonable surface fall speed distributions. Schemes assuming nonspherical mass-diameter relationships yielded reflectivity distributions closer to observed values. All four schemes examined in this study provided a better match to the observed, vertical structure of reflectivity during moderate riming than light riming periods. The comparison of observed and simulated snowfall speeds had mixed results. One BMP produced episodes of excessive cloud water at times, resulting in fall speeds that were too large. However, most schemes had frequent periods of little or no cloud water during moderate riming periods and thus underpredicted the snowfall speeds at lower levels. Short, 1-4-h periods with relatively steady snow conditions were used to compare BMP and observed size and fall speed distributions. These limited data suggest the examined BMPs underpredict fall speeds of cold-type snow habits and underrepresent aggregates larger than 4-mm diameter.

Physiologically grounded metrics of model skill: a case study estimating heat stress in intertidal populations

Models of ecological responses to climate change fundamentally assume that predictor variables, which are often measured at large scales, are to some degree diagnostic of the smaller-scale biological processes that ultimately drive patterns of abundance and distribution. Given that organisms respond physiologically to stressors, such as temperature, in highly non-linear ways, small modelling errors in predictor variables can potentially result in failures to predict mortality or severe stress, especially if an organism exists near its physiological limits. As a result, a central challenge facing ecologists, particularly those attempting to forecast future responses to environmental change, is how to develop metrics of forecast model skill (the ability of a model to predict defined events) that are biologically meaningful and reflective of underlying processes. We quantified the skill of four simple models of body temperature (a primary determinant of physiological stress) of an intertidal mussel, Mytilus californianus, using common metrics of model performance, such as root mean square error, as well as forecast verification skill scores developed by the meteorological community. We used a physiologically grounded framework to assess each model's ability to predict optimal, sub-optimal, sub-lethal and lethal physiological responses. Models diverged in their ability to predict different levels of physiological stress when evaluated using skill scores, even though common metrics, such as root mean square error, indicated similar accuracy overall. Results from this study emphasize the importance of grounding assessments of model skill in the context of an organism's physiology and, especially, of considering the implications of false-positive and false-negative errors when forecasting the ecological effects of environmental change.

Scale Dependence of Land-Atmosphere Interactions in Wet and Dry Regions as Simulated with NU-WRF over the Southwestern and South-Central United States

Large-scale forcing and land-atmosphere interactions on precipitation are investigated with NASA-Unified WRF (NU-WRF) simulations during fast transitions of ENSO phases from spring to early summer of 2010 and 2011. The model is found to capture major precipitation episodes in the 3-month simulations without resorting to nudging. However, the mean intensity of the simulated precipitation is underestimated by 46% and 57% compared with the observations in dry and wet regions in the southwestern and south-central United States, respectively. Sensitivity studies show that large-scale atmospheric forcing plays a major role in producing regional precipitation. A methodology to account for moisture contributions to individual precipitation events, as well as total precipitation, is presented under the same moisture budget framework. The analysis shows that the relative contributions of local evaporation and large-scale moisture convergence depend on the dry/wet regions and are a function of temporal and spatial scales. While the ratio of local and large-scale moisture contributions vary with domain size and weather system, evaporation provides a major moisture source in the dry region and during light rain events, which leads to greater sensitivity to soil moisture in the dry region and during light rain events. The feedback of land surface processes to large-scale forcing is well simulated, as indicated by changes in atmospheric circulation and moisture convergence. Overall, the results reveal an asymmetrical response of precipitation events to soil moisture, with higher sensitivity under dry than wet conditions. Drier soil moisture tends to suppress further existing below-normal precipitation conditions via a positive soil moisture-land surface flux feedback that could worsen drought conditions in the southwestern United States.

Verification of the isotopic composition of precipitation simulated by a regional isotope circulation model over Japan

The isotopic composition (δ(18)O and δ(2)H) of precipitation simulated by a regional isotope circulation model with a horizontal resolution of 10, 30 and 50 km was compared with observations at 56 sites over Japan in 2013. All simulations produced reasonable spatio-temporal variations in δ(18)O in precipitation over Japan, except in January. In January, simulated δ(18)O values in precipitation were higher than observed values on the Pacific side of Japan, especially during an explosively developing extratropical cyclone event. This caused a parameterisation of precipitation formulation about the large fraction of precipitated water to liquid detrained water in the lower troposphere. As a result, most water vapour that transported from the Sea of Japan precipitated on the Sea of Japan side. The isotopic composition of precipitation was a useful verification tool for the parameterisation of precipitation formulation as well as large-scale moisture transport processes in the regional isotope circulation model.

Air pollutant emissions from Chinese households: A major and underappreciated ambient pollution source

As part of the 12th Five-Year Plan, the Chinese government has developed air pollution prevention and control plans for key regions with a focus on the power, transport, and industrial sectors. Here, we investigate the contribution of residential emissions to regional air pollution in highly polluted eastern China during the heating season, and find that dramatic improvements in air quality would also result from reduction in residential emissions. We use the Weather Research and Forecasting model coupled with Chemistry to evaluate potential residential emission controls in Beijing and in the Beijing, Tianjin, and Hebei (BTH) region. In January and February 2010, relative to the base case, eliminating residential emissions in Beijing reduced daily average surface PM2.5 (particulate mater with aerodynamic diameter equal or smaller than 2.5 micrometer) concentrations by 14 ± 7 μg⋅m(-3) (22 ± 6% of a baseline concentration of 67 ± 41 μg⋅m(-3); mean ± SD). Eliminating residential emissions in the BTH region reduced concentrations by 28 ± 19 μg⋅m(-3) (40 ± 9% of 67 ± 41 μg⋅m(-3)), 44 ± 27 μg⋅m(-3) (43 ± 10% of 99 ± 54 μg⋅m(-3)), and 25 ± 14 μg⋅m(-3) (35 ± 8% of 70 ± 35 μg⋅m(-3)) in Beijing, Tianjin, and Hebei provinces, respectively. Annually, elimination of residential sources in the BTH region reduced emissions of primary PM2.5 by 32%, compared with 5%, 6%, and 58% achieved by eliminating emissions from the transportation, power, and industry sectors, respectively. We also find air quality in Beijing would benefit substantially from reductions in residential emissions from regional controls in Tianjin and Hebei, indicating the value of policies at the regional level.

Evaluation of NU-WRF rainfall forecasts for IFloodS

The Iowa Flood Studies (IFloodS) campaign was conducted in eastern Iowa as a pre-GPM-launch campaign from 1 May to 15 June 2013. During the campaign period, real time forecasts are conducted utilizing NASA-Unified Weather Research and Forecasting (NU-WRF) model to support the everyday weather briefing. In this study, two sets of the NU-WRF rainfall forecasts are evaluated with Stage IV and Multi-Radar Multi-Sensor (MRMS) Quantitative Precipitation Estimation (QPE), with the objective to understand the impact of Land Surface initialization on the predicted precipitation. NU-WRF is also compared with North American Mesoscale Forecast System (NAM) 12 km forecast. In general, NU-WRF did a good job at capturing individual precipitation events. NU-WRF is also able to replicate a better rainfall spatial distribution compare with NAM. Further sensitivity tests show that the high-resolution makes a positive impact on rainfall forecast. The two sets of NU-WRF simulations produce very close rainfall characteristics. The Land surface initialization do not show significant impact on short-term rainfall forecast, and it is largely due to the soil conditions during the field campaign period.

Confronting weather and climate models with observational data from soil moisture networks over the United States

Four land surface models in uncoupled and coupled configurations are compared to observations of daily soil moisture from 19 networks in the conterminous United States to determine the viability of such comparisons and explore the characteristics of model and observational data. First, observations are analyzed for error characteristics and representation of spatial and temporal variability. Some networks have multiple stations within an area comparable to model grid boxes; for those we find that aggregation of stations before calculation of statistics has little effect on estimates of variance, but soil moisture memory is sensitive to aggregation. Statistics for some networks stand out as unlike those of their neighbors, likely due to differences in instrumentation, calibration and maintenance. Buried sensors appear to have less random error than near-field remote sensing techniques, and heat dissipation sensors show less temporal variability than other types. Model soil moistures are evaluated using three metrics: standard deviation in time, temporal correlation (memory) and spatial correlation (length scale). Models do relatively well in capturing large-scale variability of metrics across climate regimes, but poorly reproduce observed patterns at scales of hundreds of kilometers and smaller. Uncoupled land models do no better than coupled model configurations, nor do reanalyses outperform free-running models. Spatial decorrelation scales are found to be difficult to diagnose. Using data for model validation, calibration or data assimilation from multiple soil moisture networks with different types of sensors and measurement techniques requires great caution. Data from models and observations should be put on the same spatial and temporal scales before comparison.

Calibration to improve forward model simulation of microwave emissivity at GPM frequencies over the U.S. Southern Great Plains

Better estimation of land surface microwave emissivity promises to improve over-land precipitation retrievals in the GPM era. Forward models of land microwave emissivity are available but have suffered from poor parameter specification and limited testing. Here, forward models are calibrated and the accompanying change in predictive power is evaluated. With inputs (e.g., soil moisture) from the Noah land surface model and applying MODIS LAI data, two microwave emissivity models are tested, the Community Radiative Transfer Model (CRTM) and Community Microwave Emission Model (CMEM). The calibration is conducted with the NASA Land Information System (LIS) parameter estimation subsystem using AMSR-E based emissivity retrievals for the calibration dataset. The extent of agreement between the modeled and retrieved estimates is evaluated using the AMSR-E retrievals for a separate 7-year validation period. Results indicate that calibration can significantly improve the agreement, simulating emissivity with an across-channel average root-mean-square-difference (RMSD) of about 0.013, or about 20% lower than if relying on daily estimates based on climatology. The results also indicate that calibration of the microwave emissivity model alone, as was done in prior studies, results in as much as 12% higher across-channel average RMSD, as compared to joint calibration of the land surface and microwave emissivity models. It remains as future work to assess the extent to which the improvements in emissivity estimation translate into improvements in precipitation retrieval accuracy.

The role of vegetation-microclimate feedback in promoting shrub encroachment in the northern Chihuahuan desert

Many arid and semi-arid landscapes around the world are affected by a shift from grassland to shrubland vegetation, presumably induced by climate warming, increasing atmospheric CO2 concentrations, and/or changing land use. This major change in vegetation cover is likely sustained by positive feedbacks with the physical environment. Recent research has focused on a feedback with microclimate, whereby cold intolerant shrubs increase the minimum nocturnal temperatures in their surroundings. Despite the rich literature on the impact of land cover change on local climate conditions, changes in microclimate resulting from shrub expansion into desert grasslands have remained poorly investigated. It is unclear to what extent such a feedback can affect the maximum extent of shrub expansion and the configuration of a stable encroachment front. Here, we focus on the case of the northern Chihuahuan desert, where creosotebush (Larrea tridentata) has been replacing grasslands over the past 100-150 years. We use a process-based coupled atmosphere-vegetation model to investigate the role of this feedback in sustaining shrub encroachment in the region. Simulations indicate that the feedback allows juvenile shrubs to establish in the grassland during average years and, once established, reduce their vulnerability to freeze-induced mortality by creating a warmer microclimate. Such a feedback is crucial in extreme cold winters as it may reduce shrub mortality. We identify the existence of a critical zone in the surroundings of the encroachment front, in which vegetation dynamics are bistable: in this zone, vegetation can be stable both as grassland and as shrubland. The existence of these alternative stable states explains why in most cases the shift from grass to shrub cover is found to be abrupt and often difficult to revert.

Directional analysis of the storm surge from Hurricane Sandy 2012, with applications to Charleston, New Orleans, and the Philippines

Hurricane Sandy in late October 2012 drove before it a storm surge that rose to 4.28 meters above mean lower low water at The Battery in lower Manhattan, and flooded the Hugh L. Carey automobile tunnel between Brooklyn and The Battery. This study examines the surge event in New York Harbor using the Weather Research and Forecasting (WRF) atmospheric model and the Coupled-Ocean-Atmosphere-Wave- Sediment Transport / Regional Ocean Modeling System (COAWST/ROMS). We present a new technique using directional analysis to calculate and display maps of a coastline's potential for storm surge; these maps are constructed from wind fields blowing from eight fixed compass directions. This analysis approximates the surge observed during Hurricane Sandy. The directional analysis is then applied to surge events at Charleston, South Carolina, New Orleans, Louisiana, and Tacloban City, the Philippines. Emergency managers could use these directional maps to prepare their cities for an approaching storm, on planning horizons from days to years.

Preliminary evaluation of a regional atmospheric chemical data assimilation system for environmental surveillance

We report the progress of an ongoing effort by the Air Resources Laboratory, NOAA to build a prototype regional Chemical Analysis System (ARLCAS). The ARLCAS focuses on providing long-term analysis of the three dimensional (3D) air-pollutant concentration fields over the continental U.S. It leverages expertise from the NASA Earth Science Division-sponsored Air Quality Applied Science Team (AQAST) for the state-of-science knowledge in atmospheric and data assimilation sciences. The ARLCAS complies with national operational center requirement protocols and aims to have the modeling system to be maintained by a national center. Meteorology and chemistry observations consist of land-, air- and space-based observed and quality-assured data. We develop modularized testing to investigate the efficacies of the various components of the ARLCAS. The sensitivity testing of data assimilation schemes showed that with the increment of additional observational data sets, the accuracy of the analysis chemical fields also increased incrementally in varying margins. The benefit is especially noted for additional data sets based on a different platform and/or a different retrieval algorithm. We also described a plan to apply the analysis chemical fields in environmental surveillance at the Centers for Disease Control and Prevention.

Streamflow impacts of biofuel policy-driven landscape change

Likely changes in precipitation (P) and potential evapotranspiration (PET) resulting from policy-driven expansion of bioenergy crops in the United States are shown to create significant changes in streamflow volumes and increase water stress in the High Plains. Regional climate simulations for current and biofuel cropping system scenarios are evaluated using the same atmospheric forcing data over the period 1979-2004 using the Weather Research Forecast (WRF) model coupled to the NOAH land surface model. PET is projected to increase under the biofuel crop production scenario. The magnitude of the mean annual increase in PET is larger than the inter-annual variability of change in PET, indicating that PET increase is a forced response to the biofuel cropping system land use. Across the conterminous U.S., the change in mean streamflow volume under the biofuel scenario is estimated to range from negative 56% to positive 20% relative to a business-as-usual baseline scenario. In Kansas and Oklahoma, annual streamflow volume is reduced by an average of 20%, and this reduction in streamflow volume is due primarily to increased PET. Predicted increase in mean annual P under the biofuel crop production scenario is lower than its inter-annual variability, indicating that additional simulations would be necessary to determine conclusively whether predicted change in P is a response to biofuel crop production. Although estimated changes in streamflow volume include the influence of P change, sensitivity results show that PET change is the significantly dominant factor causing streamflow change. Higher PET and lower streamflow due to biofuel feedstock production are likely to increase water stress in the High Plains. When pursuing sustainable biofuels policy, decision-makers should consider the impacts of feedstock production on water scarcity.

The role of temperature and humidity on seasonal influenza in tropical areas: Guatemala, El Salvador and Panama, 2008-2013

Background

The role of meteorological factors on influenza transmission in the tropics is less defined than in the temperate regions. We assessed the association between influenza activity and temperature, specific humidity and rainfall in 6 study areas that included 11 departments or provinces within 3 tropical Central American countries: Guatemala, El Salvador and Panama.

Method/Findings

Logistic regression was used to model the weekly proportion of laboratory-confirmed influenza positive samples during 2008 to 2013 (excluding pandemic year 2009). Meteorological data was obtained from the Tropical Rainfall Measuring Mission satellite and the Global Land Data Assimilation System. We found that specific humidity was positively associated with influenza activity in El Salvador (Odds Ratio (OR) and 95% Confidence Interval of 1.18 (1.07–1.31) and 1.32 (1.08–1.63)) and Panama (OR = 1.44 (1.08–1.93) and 1.97 (1.34–2.93)), but negatively associated with influenza activity in Guatemala (OR = 0.72 (0.6–0.86) and 0.79 (0.69–0.91)). Temperature was negatively associated with influenza in El Salvador's west-central departments (OR = 0.80 (0.7–0.91)) whilst rainfall was positively associated with influenza in Guatemala's central departments (OR = 1.05 (1.01–1.09)) and Panama province (OR = 1.10 (1.05–1.14)). In 4 out of the 6 locations, specific humidity had the highest contribution to the model as compared to temperature and rainfall. The model performed best in estimating 2013 influenza activity in Panama and west-central El Salvador departments (correlation coefficients: 0.5–0.9).

Conclusions/Significance

The findings highlighted the association between influenza activity and specific humidity in these 3 tropical countries. Positive association with humidity was found in El Salvador and Panama. Negative association was found in the more subtropical Guatemala, similar to temperate regions. Of all the study locations, Guatemala had annual mean temperature and specific humidity that were lower than the others.