ASU Scholarship Showcase

We develop a general framework to analyze the controllability of multiplex networks using multiple-relation networks and multiple-layer networks with interlayer couplings as two classes of prototypical systems. In the former, networks associated with different physical variables share the same set of nodes and in the latter, diffusion processes take place. We find that, for a multiple-relation network, a layer exists that dominantly determines the controllability of the whole network and, for a multiple-layer network, a small fraction of the interconnections can enhance the controllability remarkably. Our theory is generally applicable to other types of multiplex networks as well, leading to significant insights into the control of complex network systems with diverse structures and interacting patterns.

We investigate near-field radiative heat transfer between Indium Tin Oxide (ITO) nanowire arrays which behave as type 1 and 2 hyperbolic metamaterials. Using spatial dispersion dependent effective medium theory to model the dielectric function of the nanowires, the impact of filling fraction on the heat transfer is analyzed. Depending on the filling fraction, it is possible to achieve both types of hyperbolic modes. At 150 nm vacuum gap, the heat transfer between the nanowires with 0.5 filling fraction can be 11 times higher than that between two bulk ITOs. For vacuum gaps less than 150 nm the heat transfer increases as the filling fraction decreases. Results obtained from this study will facilitate applications of ITO nanowires as hyperbolic metamaterials for energy systems.

A film-coupled metamaterial structure is numerically investigated for enhancing the light absorption in an ultrathin photovoltaic layer of crystalline gallium arsenide (GaAs). The top subwavelength concave grating and the bottom metallic film could not only effectively trap light with the help of wave interference and magnetic resonance effects excited above the bandgap, but also practically serve as electrical contacts for photon-generated charge collection. The energy absorbed by the active layer is greatly enhanced with the help of the film-coupled metamaterial structure, resulting in significant improvement on the short-circuit current density by three times over a free-standing GaAs layer at the same thickness. The performance of the proposed light trapping structure is demonstrated to be little affected by the grating ridge width considering the geometric tolerance during fabrication. The optical absorption at oblique incidences also shows direction-insensitive behavior, which is highly desired for efficiently converting off-normal sunlight to electricity. The results would facilitate the development of next-generation ultrathin solar cells with lower cost and higher efficiency.

In this work, we report the design of a wavelength-tunable infrared metamaterial by tailoring magnetic resonance condition with the phase transition of vanadium dioxide (VO₂). Numerical simulation based on the finite-difference time-domain method shows a broad absorption peak at the wavelength of 10.9 μm when VO₂ is a metal, but it shifts to 15.1 μm when VO₂ changes to dielectric phase below its phase transition temperature of 68 °C. The large tunability of 38.5% in the resonance wavelength stems from the different excitation conditions of magnetic resonance mediated by plasmon in metallic VO₂ but optical phonons in dielectric VO₂. The physical mechanism is elucidated with the aid of electromagnetic field distribution at the resonance wavelengths. A hybrid magnetic resonance mode due to the plasmon-phonon coupling is also discussed. The results here would be beneficial for active control of thermal radiation in novel electronic, optical, and thermal devices.

In this work, we numerically demonstrate an infrared (IR) frequency-tunable selective thermal emitter made of graphene-covered silicon carbide (SiC) gratings. Rigorous coupled-wave analysis shows temporally-coherent emission peaks associated with magnetic polariton (MP), whose resonance frequency can be dynamically tuned within the phonon absorption band of SiC by varying graphene chemical potential. An analytical inductor–capacitor circuit model is introduced to quantitatively predict the resonance frequency and further elucidate the mechanism for the tunable emission peak. The effects of grating geometric parameters, such as grating height, groove width and grating period, on the selective emission peak are explored. The direction-independent behavior of MP and associated coherent emission are also demonstrated. Moreover, by depositing four layers of graphene sheets onto the SiC gratings, a large tunability of 8.5% in peak frequency can be obtained to yield the coherent emission covering a broad frequency range from 820 to 890 cm^-1. The novel tunable metamaterial could pave the way to a new class of tunable thermal sources in the IR region.

The sensitivity of Earth’s wetlands to observed shifts in global precipitation and temperature patterns and their ability to produce large quantities of methane gas are key global change questions. We present a microwave satellite-based approach for mapping fractional surface water (FW) globally at 25-km resolution. The approach employs a land cover-supported, atmospherically-corrected dynamic mixture model applied to 20+ years (1992–2013) of combined, daily, passive/active microwave remote sensing data. The resulting product, known as Surface Water Microwave Product Series (SWAMPS), shows strong microwave sensitivity to sub-grid scale open water and inundated wetlands comprising open plant canopies. SWAMPS’ FW compares favorably (R² = 91%–94%) with higher-resolution, global-scale maps of open water from MODIS and SRTM-MOD44W. Correspondence of SWAMPS with open water and wetland products from satellite SAR in Alaska and the Amazon deteriorates when exposed wetlands or inundated forests captured by the SAR products were added to the open water fraction reflecting SWAMPS’ inability to detect water underneath the soil surface or beneath closed forest canopies. Except for a brief period of drying during the first 4 years of observation, the inundation extent for the global domain excluding the coast was largely stable. Regionally, inundation in North America is advancing while inundation is on the retreat in Tropical Africa and North Eurasia. SWAMPS provides a consistent and long-term global record of daily FW dynamics, with documented accuracies suitable for hydrologic assessment and global change-related investigations.

A warming climate is altering land-atmosphere exchanges of carbon, with a potential for increased vegetation productivity as well as the mobilization of permafrost soil carbon stores. Here we investigate land-atmosphere carbon dioxide (CO₂) cycling through analysis of net ecosystem productivity (NEP) and its component fluxes of gross primary productivity (GPP) and ecosystem respiration (ER) and soil carbon residence time, simulated by a set of land surface models (LSMs) over a region spanning the drainage basin of Northern Eurasia. The retrospective simulations cover the period 1960–2009 at 0.5° resolution, which is a scale common among many global carbon and climate model simulations. Model performance benchmarks were drawn from comparisons against both observed CO₂ fluxes derived from site-based eddy covariance measurements as well as regional-scale GPP estimates based on satellite remote-sensing data.

The site-based comparisons depict a tendency for overestimates in GPP and ER for several of the models, particularly at the two sites to the south. For several models the spatial pattern in GPP explains less than half the variance in the MODIS MOD17 GPP product. Across the models NEP increases by as little as 0.01 to as much as 0.79 g C m^-2 yr^-2, equivalent to 3 to 340 % of the respective model means, over the analysis period. For the multimodel average the increase is 135 % of the mean from the first to last 10 years of record (1960–1969 vs. 2000–2009), with a weakening CO₂ sink over the latter decades. Vegetation net primary productivity increased by 8 to 30 % from the first to last 10 years, contributing to soil carbon storage gains. The range in regional mean NEP among the group is twice the multimodel mean, indicative of the uncertainty in CO₂ sink strength.

The models simulate that inputs to the soil carbon pool exceeded losses, resulting in a net soil carbon gain amid a decrease in residence time. Our analysis points to improvements in model elements controlling vegetation productivity and soil respiration as being needed for reducing uncertainty in land-atmosphere CO₂ exchange. These advances will require collection of new field data on vegetation and soil dynamics, the development of benchmarking data sets from measurements and remote-sensing observations, and investments in future model development and intercomparison studies.

Dynamical processes occurring on the edges in complex networks are relevant to a variety of real-world situations. Despite recent advances, a framework for edge controllability is still required for complex networks of arbitrary structure and interaction strength. Generalizing a previously introduced class of processes for edge dynamics, the switchboard dynamics, and exploit- ing the exact controllability theory, we develop a universal framework in which the controllability of any node is exclusively determined by its local weighted structure. This framework enables us to identify a unique set of critical nodes for control, to derive analytic formulas and articulate efficient algorithms to determine the exact upper and lower controllability bounds, and to evaluate strongly structural controllability of any given network. Applying our framework to a large number of model and real-world networks, we find that the interaction strength plays a more significant role in edge controllability than the network structure does, due to a vast range between the bounds determined mainly by the interaction strength. Moreover, transcriptional regulatory networks and electronic circuits are much more strongly structurally controllable (SSC) than other types of real-world networks, directed networks are more SSC than undirected networks, and sparse networks are typically more SSC than dense networks.

Soil temperature (T_s) change is a key indicator of the dynamics of permafrost. On seasonal and interannual timescales, the variability of T_s determines the active-layer depth, which regulates hydrological soil properties and biogeochemical processes. On the multi-decadal scale, increasing T_s not only drives permafrost thaw/retreat but can also trigger and accelerate the decomposition of soil organic carbon. The magnitude of permafrost carbon feedbacks is thus closely linked to the rate of change of soil thermal regimes. In this study, we used nine process-based ecosystem models with permafrost processes, all forced by different observation-based climate forcing during the period 1960–2000, to characterize the warming rate of T_s in permafrost regions. There is a large spread of T_s trends at 20 cm depth across the models, with trend values ranging from 0.010 ± 0.003 to 0.031 ± 0.005 °C yr^-1. Most models show smaller increase in T_s with increasing depth. Air temperature (Tsub>a) and longwave downward radiation (LWDR) are the main drivers of T_s trends, but their relative contributions differ amongst the models. Different trends of LWDR used in the forcing of models can explain 61 % of their differences in T_s trends, while trends of T_a only explain 5 % of the differences in T_s trends. Uncertain climate forcing contributes a larger uncertainty in T_s trends (0.021 ± 0.008 °C yr^-1, mean ± standard deviation) than the uncertainty of model structure (0.012 ± 0.001 °C yr^-1), diagnosed from the range of response between different models, normalized to the same forcing. In addition, the loss rate of near-surface permafrost area, defined as total area where the maximum seasonal active-layer thickness (ALT) is less than 3 m loss rate, is found to be significantly correlated with the magnitude of the trends of T_s at 1 m depth across the models (R = −0.85, P = 0.003), but not with the initial total near-surface permafrost area (R = −0.30, P = 0.438). The sensitivity of the total boreal near-surface permafrost area to T_s at 1 m is estimated to be of −2.80 ± 0.67 million km²°C^-1. Finally, by using two long-term LWDR data sets and relationships between trends of LWDR and T_s across models, we infer an observation-constrained total boreal near-surface permafrost area decrease comprising between 39 ± 14  ×  10³ and 75 ± 14  ×  10³km²yr^-1 from 1960 to 2000. This corresponds to 9–18 % degradation of the current permafrost area.

A realistic simulation of snow cover and its thermal properties are important for accurate modelling of permafrost. We analyse simulated relationships between air and near-surface (20  cm) soil temperatures in the Northern Hemisphere permafrost region during winter, with a particular focus on snow insulation effects in nine land surface models, and compare them with observations from 268 Russian stations. There are large cross-model differences in the simulated differences between near-surface soil and air temperatures (ΔT; 3 to 14 °C), in the sensitivity of soil-to-air temperature (0.13 to 0.96 °C °^C-1), and in the relationship between ΔT and snow depth. The observed relationship between ΔT and snow depth can be used as a metric to evaluate the effects of each model's representation of snow insulation, hence guide improvements to the model's conceptual structure and process parameterisations. Models with better performance apply multilayer snow schemes and consider complex snow processes. Some models show poor performance in representing snow insulation due to underestimation of snow depth and/or overestimation of snow conductivity. Generally, models identified as most acceptable with respect to snow insulation simulate reasonable areas of near-surface permafrost (13.19 to 15.77 million  km²). However, there is not a simple relationship between the sophistication of the snow insulation in the acceptable models and the simulated area of Northern Hemisphere near-surface permafrost, because several other factors, such as soil depth used in the models, the treatment of soil organic matter content, hydrology and vegetation cover, also affect the simulated permafrost distribution.