High-resolution climate simulation of the last glacial maximum

The climate of the last glacial maximum (LGM) is simulated with a high-resolution atmospheric general circulation model, the NCAR CCM3 at spectral truncation of T170, corresponding to a grid cell size of roughly 75 km. The purpose of the study is to assess whether there are significant benefits from the higher resolution simulation compared to the lower resolution simulation associated with the role of topography. The LGM simulations were forced with modified CLIMAP sea ice distribution and sea surface temperatures (SST) reduced by 1°C, ice sheet topography, reduced CO₂, and 21,000 BP orbital parameters. The high-resolution model captures modern climate reasonably well, in particular the distribution of heavy precipitation in the tropical Pacific. For the ice age case, surface temperature simulated by the high-resolution model agrees better with those of proxy estimates than does the low-resolution model. Despite the fact that tropical SSTs were only 2.1°C less than the control run, there are many lowland tropical land areas 4–6°C colder than present. Comparison of T170 model results with the best constrained proxy temperature estimates (noble gas concentrations in groundwater) now yield no significant differences between model and observations. There are also significant upland temperature changes in the best resolved tropical mountain belt (the Andes). We provisionally attribute this result in part as resulting from decreased lateral mixing between ocean and land in a model with more model grid cells. A longstanding model-data discrepancy therefore appears to be resolved without invoking any unusual model physics. The response of the Asian summer monsoon can also be more clearly linked to local geography in the high-resolution model than in the low-resolution model; this distinction should enable more confident validation of climate proxy data with the high-resolution model. Elsewhere, an inferred salinity increase in the subtropical North Atlantic may have significant implications for ocean circulation changes during the LGM. A large part of the Amazon and Congo Basins are simulated to be substantially drier in the ice age—consistent with many (but not all) paleo data. These results suggest that there are considerable benefits derived from high-resolution model regarding regional climate responses, and that observationalists can now compare their results with models that resolve geography at a resolution comparable to that which the proxy data represent.     

Advances in understanding the “Perfect Monsoon-influenced Typhoon”: Summary from International Conference on Typhoon Morakot (2009)

Typhoon Morakot (2009) produced 2855 mm of rain and was the deadliest typhoon to impact Taiwan with 619 deaths and 76 missing persons, including a landslide that wiped out an entire village. While Morakot did not exceed the heaviest 24-h rain record, the combination of heavy rain and long duration that led to the record accumulation is attributed to the southwest summer monsoon influence on the typhoon. Thus, a special combination of factors was involved in the Morakot disaster: (i) Strong southwesterly monsoon winds; (ii) Convergence between the typhoon circulation and monsoon flow to form an east-west oriented convective band over the Taiwan Strait that was quasi-stationary and long-lasting; (iii) A typhoon in a specific location relative to the Central Mountain Range and moving slowly; and (iv) Steep topography that provided rapid lifting of the moist air stream. The contributions of each of these four factors in leading to the Morakot disaster are reviewed primarily based on new research presented at the International Conference on Typhoon Morakot (2009). Historical data sets, new Doppler radar observations, and numerical modeling have advanced the understanding of the special conditions of monsoon-influenced typhoons such as Morakot. This research is also leading to modifications of existing and development of new forecasting tools. Gaps in scientific understanding, limits to the predictability, and requirements for advanced forecast guidance tools are described that are challenges to improved warnings of these extreme precipitation and flooding events in monsoon-influenced typhoons.     

Evaluation of the vegetated urban canopy model (VUCM) and its impacts on urban boundary layer simulation

The vegetated urban canopy model (VUCM) is implemented in a meteorological model, the Regional Atmospheric Modeling System (RAMS), for urban atmospheric modeling. The VUCM includes various urban physical processes such as in-canyon radiative transfer, turbulent energy exchanges, substrate heat conduction, and in-canyon momentum drag. The coupled model RAMS/VUCM is evaluated and then used to examine its impacts on the dynamic and thermodynamic structure of the urban boundary layer (UBL) in the Seoul metropolitan area. The spatial pattern of the nocturnal urban heat island (UHI) in Seoul is quite well simulated by the RAMS/VUCM. A statistical evaluation of 2-m air temperature reveals a significant improvement in model performance, especially in the nighttime. The RAMS/VUCM simulates the diurnal variations of surface energy balance fluxes realistically. This contributes to a reasonable UBL formation. A weakly unstable UBL is formed in the nighttime with UBL heights of about 100–200 m. When urban surfaces are represented in the RAMS using a land surface model of the Land Ecosystem-Atmosphere Feedback (LEAF), the RAMS/LEAF produces strong cold biases and thus fails to simulate UHI formation. This is due to the poor representation or absence of important urban physical processes in the RAMS/LEAF. This study implies that urban physical processes should be included in numerical models in order to reasonably simulate meteorology and air quality in urban areas and that the VUCM is one of the promising urban canopy models.     

Influence of leaf water potential on diurnal changes in CO2 and water vapour fluxes

Mass and energy fluxes between the atmosphere and vegetation are driven by meteorological variables, and controlled by plant water status, which may change more markedly diurnally than soil water. We tested the hypothesis that integration of dynamic changes in leaf water potential may improve the simulation of CO₂ and water fluxes over a wheat canopy. Simulation of leaf water potential was integrated into a comprehensive model (the ChinaAgrosys) of heat, water and CO₂ fluxes and crop growth. Photosynthesis from individual leaves was integrated to the canopy by taking into consideration the attenuation of radiation when penetrating the canopy. Transpiration was calculated with the Shuttleworth-Wallace model in which canopy resistance was taken as a link between energy balance and physiological regulation. A revised version of the Ball-Woodrow-Berry stomatal model was applied to produce a new canopy resistance model, which was validated against measured CO₂ and water vapour fluxes over winter wheat fields in Yucheng (36°57′ N, 116°36′ E, 28 m above sea level) in the North China Plain during 1997, 2001 and 2004. Leaf water potential played an important role in causing stomatal conductance to fall at midday, which caused diurnal changes in photosynthesis and transpiration. Changes in soil water potential were less important. Inclusion of the dynamics of leaf water potential can improve the precision of the simulation of CO₂ and water vapour fluxes, especially in the afternoon under water stress conditions.     

Observing and Simulating Wind-Turbine Wakes During the Evening Transition

Wind-turbine-wake evolution during the evening transition introduces variability to wind-farm power production at a time of day typically characterized by high electricity demand. During the evening transition, the atmosphere evolves from an unstable to a stable regime, and vertical stratification of the wind profile develops as the residual planetary boundary layer decouples from the surface layer. The evolution of wind-turbine wakes during the evening transition is examined from two perspectives: wake observations from single turbines, and simulations of multiple turbine wakes using the mesoscale Weather Research and Forecasting (WRF) model. Throughout the evening transition, the wake’s wind-speed deficit and turbulence enhancement are confined within the rotor layer when the atmospheric stability changes from unstable to stable. The height variations of maximum upwind-downwind differences of wind speed and turbulence intensity gradually decrease during the evening transition. After verifying the WRF-model-simulated upwind wind speed, wind direction and turbulent kinetic energy profiles with observations, the wind-farm-scale wake evolution during the evening transition is investigated using the WRF-model wind-farm parametrization scheme. As the evening progresses, due to the presence of the wind farm, the modelled hub-height wind-speed deficit monotonically increases, the relative turbulence enhancement at hub height grows by 50%, and the downwind surface sensible heat flux increases, reducing surface cooling. Overall, the intensifying wakes from upwind turbines respond to the evolving atmospheric boundary layer during the evening transition, and undermine the power production of downwind turbines in the evening.     

Effects of climate variability on the distribution and fishing conditions of yellowfin tuna (Thunnus albacares) in the western Indian Ocean

Variations in the abundance and distribution of pelagic tuna populations have been associated with large-scale climate indices such as the Southern Oscillation Index in the Pacific Ocean and the North Atlantic Oscillation in the Atlantic Ocean. Similarly to the Pacific and Atlantic, variability in the distribution and catch rates of tuna species have also been observed in association with the Indian Ocean Dipole (IOD), a basin-scale pattern of sea surface and subsurface temperatures that affect climate in the Indian Ocean. The environmental processes associated with the IOD that drive variability in tuna populations, however, are largely unexplored. To better understand these processes, we investigated longline catch rates of yellowfin tuna and their distributions in the western Indian Ocean in relation to IOD events, sea surface water temperatures (SST) and estimates of net primary productivity (NPP). Catch per unit effort (CPUE) was observed to be negatively correlated to the IOD with a periodicity centred around 4 years. During positive IOD events, SSTs were relatively higher, NPP was lower, CPUE decreased and catch distributions were restricted to the northern and western margins of the western Indian Ocean. During negative IOD events, lower SSTs and higher NPP were associated with increasing CPUE, particularly in the Arabian Sea and seas surrounding Madagascar, and catches expanded into central regions of the western Indian Ocean. These findings provide preliminary insights into some of the key environmental features driving the distribution of yellowfin tuna in the western Indian Ocean and associated variability in fisheries catches.     

The influence of coupled physical swelling and chemical reactions on deformable geomaterials

Coupled thermo‐hydro‐mechanical‐chemical modelling has attracted attention in past decades due to many contemporary geotechnical engineering applications (e.g., waste disposal, carbon capture and storage). However, molecular‐scale interactions within geomaterials (e.g., swelling and dissolution/precipitation) have a significant influence on the mechanical behaviour, yet are rarely incorporated into existing Thermal‐Hydro‐Mechanical‐Chemical (THMC) frameworks. This paper presents a new coupled hydro‐mechanical‐chemical constitutive model to bridge molecular‐scale interactions with macro‐physical deformation by combining the swelling and dissolution/precipitation through an extension of the new mixture‐coupling theory. Entropy analysis of the geomaterial system provides dissipation energy, and Helmholtz free energy gives the relationship between solids and fluids. Numerical simulation is used to compare with the selected recognized models, which demonstrates that the swelling and dissolution/precipitation processes may have a significant influence on the mechanical deformation of the geomaterials.     

Water source dynamics of high Arctic river basins

Arctic river basins are amongst the most vulnerable to climate change. However, there is currently limited knowledge of the hydrological processes that govern flow dynamics in Arctic river basins. We address this research gap using natural hydrochemical and isotopic tracers to identify water sources that contributed to runoff in river basins spanning a gradient of glacierization (0–61%) in Svalbard during summer 2010 and 2011. Spatially distinct hydrological processes operating over diurnal, weekly and seasonal timescales were characterized by river hydrochemistry and isotopic composition. Two conceptual water sources (‘meltwater’ and ‘groundwater’) were identified and used as a basis for end‐member mixing analyses to assess seasonal and year‐to‐year variability in water source dynamics. In glacier‐fed rivers, meltwater dominated flows at all sites (typically >80%) with the highest contributions observed at the beginning of each study period in early July when snow cover was most extensive. Rivers in non‐glacierized basins were sourced initially from snowmelt but became increasingly dependent on groundwater inputs (up to 100% of total flow volume) by late summer. These hydrological changes were attributed to the depletion of snowpacks and enhanced soil water storage capacity as the active layer expanded throughout each melt season. These findings provide insight into the processes that underpin water source dynamics in Arctic river systems and potential future changes in Arctic hydrology that might be expected under a changing climate. Copyright © 2013 John Wiley & Sons, Ltd.     

Quantitative estimation of submarine groundwater discharge using airborne thermal infrared data acquired at two different tidal heights

Submarine groundwater discharge (SGD) plays an important role in coastal biogeochemical processes and hydrological cycles, particularly off volcanic islands in oligotrophic oceans. However, the spatial and temporal variations of SGD are still poorly understood owing to difficulty in taking rapid SGD measurements over a large scale. In this study, we used four airborne thermal infrared surveys (twice each during high and low tides) to quantify the spatiotemporal variations of SGD over the entire coast of Jeju Island, Korea. On the basis of an analytical model, we found a linear positive correlation between the thermal anomaly and squares of the groundwater discharge velocity and a negative exponential correlation between the anomaly and water depth (including tide height and bathymetry). We then derived a new equation for quantitatively estimating the SGD flow rates from thermal anomalies acquired at two different tide heights. The proposed method was validated with the measured SGD flow rates using a current meter at Gongcheonpo Beach. We believe that the method can be effectively applied for rapid estimation of SGD over coastal areas, where fresh groundwater discharge is significant, using airborne thermal infrared surveys.