Indirect air–sea interactions simulated with a coupled turbulence-resolving model

A turbulence-resolving parallelized atmospheric large-eddy simulation model (PALM) has been applied to study turbulent interactions between the humid atmospheric boundary layer (ABL) and the salt water oceanic mixed layer (OML). The most energetic three-dimensional turbulent eddies in the ABL–OML system (convective cells) were explicitly resolved in these simulations. This study considers a case of shear-free convection in the coupled ABL–OML system. The ABL–OML coupling scheme used the turbulent fluxes at the bottom of the ABL as upper boundary conditions for the OML and the sea surface temperature at the top of the OML as lower boundary conditions for the ABL. The analysis of the numerical experiment confirms that the ABL–OML interactions involve both the traditional direct coupling mechanism and much less studied indirect coupling mechanism (Garrett Dyn Atmos Ocean 23:19–34, 1996). The direct coupling refers to a common flux-gradient representation of the air–sea exchange, which is controlled by the temperature difference across the air–water interface. The indirect coupling refers to thermal instability of the Rayleigh–Benard convection, which is controlled by the temperature difference across the entire mixed layer through formation of the large convective eddies or cells. The indirect coupling mechanism in these simulations explained up to 45 % of the ABL–OML co-variability on the turbulent scales. Despite relatively small amplitude of the sea surface temperature fluctuations, persistence of the OML cells organizes the ABL convective cells. Water downdrafts in the OML cells tend to be collocated with air updrafts in the ABL cells. The study concludes that the convective structures in the ABL and the OML are co-organized. The OML convection controls the air–sea turbulent exchange in the quasi-equilibrium convective ABL–OML system.     

The nested superblock approach for regional-scale analytic element models

A new approach is presented for improving the computational efficiency of regional-scale ground water models based on the analytic element method (AEM). The algorithm is an extension of the existing "superblock" algorithm, which combines the effects of multiple analytic elements into Laurent series and Taylor series (superblock expansions). With the new "nested superblock" formulation, Laurent series are nested in a hierarchical (quad-tree) data structure with direct mathematical relationships between parent and child superblock coefficients. Nested superblocks significantly accelerate the evaluation of the complex potential and discharge function in models that contain a large number of analytic elements at multiple scales. This evaluation process, the primary computational cost of AEM models, is required to determine the element coefficients, generate contour plots, and trace pathlines. The performance of the nested superblocks is demonstrated with a simplified model based on the Lake Ontario watershed geometry comprising thousands of hydrogeologic features at multiple geographic scales.     

Vertical structure of recent arctic warming from observed data and reanalysis products

Spatiotemporal patterns of recent (1979–2008) air temperature trends are evaluated using three reanalysis datasets and radiosonde data. Our analysis demonstrates large discrepancies between the reanalysis datasets, possibly due to differences in the data assimilation procedures as well as sparseness and inhomogeneity of high-latitude observations. We test the robustness of arctic tropospheric warming based on the ERA-40 dataset. ERA-40 Arctic atmosphere temperatures tend to be closer to the observed ones in terms of root mean square error compared to other reanalysis products used in the article. However, changes in the ERA-40 data assimilation procedure produce unphysical jumps in atmospheric temperatures, which may be the likely reason for the elevated tropospheric warming trend in 1979–2002. NCEP/NCAR Reanalysis data show that the near-surface upward temperature trend over the same period is greater than the tropospheric trend, which is consistent with direct radiosonde observations and inconsistent with ERA-40 results. A change of sign in the winter temperature trend from negative to positive in the late 1980s is documented in the upper troposphere/lower stratosphere with a maximum over the Canadian Arctic, based on radiosonde data. This change from cooling to warming tendency is associated with weakening of the stratospheric polar vortex and shift of its center toward the Siberian coast and possibly can be explained by the changes in the dynamics of the Arctic Oscillation. This temporal pattern is consistent with multi-decadal variations of key arctic climate parameters like, for example, surface air temperature and oceanic freshwater content. Elucidating the mechanisms behind these changes will be critical to understanding the complex nature of high-latitude variability and its impact on global climate change.     

Interannual variability of Mediterranean evaporation and its relation to regional climate

Gridded monthly evaporation data for 1958–2006 from the Woods Hole Oceanographic Institution data set are used to investigate interannual variability of Mediterranean evaporation during cold and hot seasons and its relation to regional atmospheric dynamics, sea surface temperature and atmospheric elements of the hydrological cycle. The first EOF mode of Mediterranean evaporation, explaining more than 50% of its total variance, is characterized by the monopole pattern both in winter and summer. However, despite structural similarity, the EOF-1 of Mediterranean evaporation is affected by different climate signals in cold and hot seasons. During winter the EOF-1 is associated with the East Atlantic teleconnection pattern. In summer, there is indication of tropical influence on the EOF-1 of Mediterranean evaporation (presumably from Asian monsoon). Both in winter and summer, principal components of EOF-1 demonstrate clear interdecadal signals (with a stronger signature in summer) associated with large sea surface temperature anomalies. The results of a sensitivity analysis suggest that in winter both the meridional wind and the vertical gradient of saturation specific humidity (GSSH) near the sea surface contribute to the interdecadal evaporation signal. In summer, however, it is likely that the signal is more related to GSSH. Our analysis did not reveal significant links between the Mediterranean evaporation and the North Atlantic Oscillation in any season. The EOF-2 of evaporation accounts for 20% (11%) of its total variance in winter (in summer). Both in winter and summer the EOF-2 is characterized by a zonal dipole with opposite variations of evaporation in western and eastern parts of the Mediterranean Sea. This mode is associated presumably with smaller scale (i.e., local) effects of atmospheric dynamics. Seasonality of the leading modes of the Mediterranean evaporation is also clearly seen in the character of their links to atmospheric elements of the regional hydrological cycle. In particular, significant links to precipitation in some regions have been found in winter, but not in summer.     

秘鲁西北部早白垩世Oyotún组火山岩层序、年代学及地球化学特征——对古Farallon板块俯冲作用细节的揭示

Oyotún组火山岩为中安第斯西北缘中生代出露范围最广的一套火山岩,对研究中安第斯早白垩世壳幔相互作用具有重要意义。研究区西部Oyotún组火山岩出露面积小,厚度薄,岩性单一,而东部出露范围大,厚度厚,岩性复杂。LA-ICP-MS锆石U-Pb定年表明,Oyotún组火山岩时代集中在139.7~131.0Ma之间。主量元素SiO_2含量介于54.44%~72.33%之间,中等富MgO(0.56%~1.55%),偏铝-强过铝,A/CNK=0.98~3.51,富集Cs、Rb、Th、U、Pb和轻稀土元素,无明显Eu异常,亏损Nb、Ta、Ti、P等,具明显的火山弧特征。Oyotún组火山岩εHf(t)值在3.37~8.42之间,且底部εHf(t)值显著低于顶部,指示其可能形成于活动陆缘弧后伸展背景下,由亏损地幔衍生的新生玄武质下地壳,经板片流体作用而产生的低程度部分熔融过程,并受到太古宙古老陆壳物质的混染,其混染程度随时代变新而逐渐减小。     

Soil compaction under different management practices in a Croatian vineyard

The use of machinery in vineyards is increasing soil compaction and erosion. However, there is a lack of information about the impacts of different management practices on soil conditions of vineyards. Therefore, the aim of this study was to assess soil compaction in Croatian vineyards under four different management systems: no-tillage (NT) system, conventional tillage (CT), yearly inversed grass covered (INV-GC) and tillage managed (INV-T) treatments. Four key top-soil (0–20 cm) parameters were assessed in the different land uses: bulk density (BD), penetration resistance (PR), soil water content (SWC) and carbon dioxide (CO₂) fluxes. Tractor traffic increased the BD and PR in all treatments, with exception of CT treatment, as consequence of tillage. SWC showed higher values in INV-GC treatment during the dry period; meanwhile, it was similar during the wet season in every management type. Lower CO₂ fluxes were found in INV-GC and NT than in the CT and INV-T treatments.     

Osmium isotope and highly siderophile element systematics of the lunar crust

Coupled ¹⁸⁷Os/¹⁸⁸Os and highly siderophile element (HSE: Os, Ir, Ru, Pt, Pd, and Re) abundance data are reported for pristine lunar crustal rocks 60025, 62255, 65315 (ferroan anorthosites, FAN) and 76535, 78235, 77215 and a norite clast in 15455 (magnesian-suite rocks, MGS). Osmium isotopes permit more refined discrimination than previously possible of samples that have been contaminated by meteoritic additions and the new results show that some rocks, previously identified as pristine, contain meteorite-derived HSE. Low HSE abundances in FAN and MGS rocks are consistent with derivation from a strongly HSE-depleted lunar mantle. At the time of formation, the lunar floatation crust, represented by FAN, had 1.4 ± 0.3 pg g^? 1 Os, 1.5 ± 0.6 pg g^? 1 Ir, 6.8 ± 2.7 pg g^? 1 Ru, 16 ± 15 pg g^? 1 Pt, 33 ± 30 pg g^? 1 Pd and 0.29 ± 0.10 pg g^? 1 Re (～ 0.00002 × CI) and Re/Os ratios that were modestly elevated (¹⁸⁷Re/¹⁸⁸Os = 0.6 to 1.7) relative to CI chondrites. MGS samples are, on average, characterised by more elevated HSE abundances (～ 0.00007 × CI) compared with FAN. This either reflects contrasting mantle-source HSE characteristics of FAN and MGS rocks, or different mantle–crust HSE fractionation behaviour during production of these lithologies. Previous studies of lunar impact-melt rocks have identified possible elevated Ru and Pd in lunar crustal target rocks. The new results provide no supporting evidence for such enrichments.If maximum estimates for HSE in the lunar mantle are compared with FAN and MGS averages, crust–mantle concentration ratios (D-values) must be ≤ 0.3. Such D-values are broadly similar to those estimated for partitioning between the terrestrial crust and upper mantle, with the notable exception of Re. Given the presumably completely different mode of origin for the primary lunar floatation crust and tertiary terrestrial continental crust, the potential similarities in crust–mantle HSE partitioning for the Earth and Moon are somewhat surprising. Low HSE abundances in the lunar crust, coupled with estimates of HSE concentrations in the lunar mantle implies there may be a ‘missing component’ of late-accreted materials (as much as 95%) to the Moon if the Earth/Moon mass-flux estimates are correct and terrestrial mantle HSE abundances were established by late accretion.     

Moveout approximation for horizontal transversely isotropic and vertical transversely isotropic layered medium. Part II: effective model‡

We use residual moveouts measured along continuous full azimuth reflection angle gathers, in order to obtain effective horizontal transversely isotropic model parameters. The angle gathers are generated through a special angle domain imaging system, for a wide range of reflection angles and full range of phase velocity azimuths. The estimation of the effective model parameters is performed in two stages. First, the background horizontal transversely isotropic (HTI)/vertical transversely isotropic (VTI) layered model is used, along with the values of reflection angles, for converting the measured residual moveouts (or traveltime errors) into azimuthally dependent normal moveout (NMO) velocities. Then we apply a digital Fourier transform to convert the NMO velocities into azimuthal wavenumber domain, in order to obtain the effective HTI model parameters: vertical time, vertical compression velocity, Thomsen parameter delta and the azimuth of the medium axis of symmetry. The method also provides a reliability criterion of the HTI assumption. The criterion shows whether the medium possesses the HTI type of symmetry, or whether the azimuthal dependence of the residual traveltime indicates to a more complex azimuthal anisotropy. The effective model used in this approach is defined for a 1D structure with a set of HTI, VTI and isotropic layers (with at least one HTI layer). We describe and analyse the reduction of a multi‐layer structure into an equivalent effective HTI model. The equivalent model yields the same NMO velocity and the same offset azimuth on the Earth's surface as the original layered structure, for any azimuth of the phase velocity. The effective model approximates the kinematics of an HTI/VTI layered structure using only a few parameters. Under the hyperbolic approximation, the proposed effective model is exact.