| 1. |
- Niemi, MEK, et al.
(författare)
-
- 2021
-
swepub:Mat__t
|
|
| 2. |
- Romagnoni, A, et al.
(författare)
-
Comparative performances of machine learning methods for classifying Crohn Disease patients using genome-wide genotyping data
- 2019
-
Ingår i: Scientific reports. - : Springer Science and Business Media LLC. - 2045-2322. ; 9:1, s. 10351-
-
Tidskriftsartikel (refereegranskat)abstract
- Crohn Disease (CD) is a complex genetic disorder for which more than 140 genes have been identified using genome wide association studies (GWAS). However, the genetic architecture of the trait remains largely unknown. The recent development of machine learning (ML) approaches incited us to apply them to classify healthy and diseased people according to their genomic information. The Immunochip dataset containing 18,227 CD patients and 34,050 healthy controls enrolled and genotyped by the international Inflammatory Bowel Disease genetic consortium (IIBDGC) has been re-analyzed using a set of ML methods: penalized logistic regression (LR), gradient boosted trees (GBT) and artificial neural networks (NN). The main score used to compare the methods was the Area Under the ROC Curve (AUC) statistics. The impact of quality control (QC), imputing and coding methods on LR results showed that QC methods and imputation of missing genotypes may artificially increase the scores. At the opposite, neither the patient/control ratio nor marker preselection or coding strategies significantly affected the results. LR methods, including Lasso, Ridge and ElasticNet provided similar results with a maximum AUC of 0.80. GBT methods like XGBoost, LightGBM and CatBoost, together with dense NN with one or more hidden layers, provided similar AUC values, suggesting limited epistatic effects in the genetic architecture of the trait. ML methods detected near all the genetic variants previously identified by GWAS among the best predictors plus additional predictors with lower effects. The robustness and complementarity of the different methods are also studied. Compared to LR, non-linear models such as GBT or NN may provide robust complementary approaches to identify and classify genetic markers.
|
|
| 3. |
|
|
| 4. |
|
|
| 5. |
- Birney, Ewan, et al.
(författare)
-
Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project
- 2007
-
Ingår i: Nature. - : Springer Science and Business Media LLC. - 0028-0836 .- 1476-4687. ; 447:7146, s. 799-816
-
Tidskriftsartikel (refereegranskat)abstract
- We report the generation and analysis of functional data from multiple, diverse experiments performed on a targeted 1% of the human genome as part of the pilot phase of the ENCODE Project. These data have been further integrated and augmented by a number of evolutionary and computational analyses. Together, our results advance the collective knowledge about human genome function in several major areas. First, our studies provide convincing evidence that the genome is pervasively transcribed, such that the majority of its bases can be found in primary transcripts, including non-protein-coding transcripts, and those that extensively overlap one another. Second, systematic examination of transcriptional regulation has yielded new understanding about transcription start sites, including their relationship to specific regulatory sequences and features of chromatin accessibility and histone modification. Third, a more sophisticated view of chromatin structure has emerged, including its inter-relationship with DNA replication and transcriptional regulation. Finally, integration of these new sources of information, in particular with respect to mammalian evolution based on inter- and intra-species sequence comparisons, has yielded new mechanistic and evolutionary insights concerning the functional landscape of the human genome. Together, these studies are defining a path for pursuit of a more comprehensive characterization of human genome function.
|
|
| 6. |
- Craddock, Nick, et al.
(författare)
-
Genome-wide association study of CNVs in 16,000 cases of eight common diseases and 3,000 shared controls
- 2010
-
Ingår i: Nature. - : Springer Science and Business Media LLC. - 0028-0836 .- 1476-4687. ; 464:7289, s. 713-720
-
Tidskriftsartikel (refereegranskat)abstract
- Copy number variants (CNVs) account for a major proportion of human genetic polymorphism and have been predicted to have an important role in genetic susceptibility to common disease. To address this we undertook a large, direct genome-wide study of association between CNVs and eight common human diseases. Using a purpose-designed array we typed,19,000 individuals into distinct copy-number classes at 3,432 polymorphic CNVs, including an estimated similar to 50% of all common CNVs larger than 500 base pairs. We identified several biological artefacts that lead to false-positive associations, including systematic CNV differences between DNAs derived from blood and cell lines. Association testing and follow-up replication analyses confirmed three loci where CNVs were associated with disease-IRGM for Crohn's disease, HLA for Crohn's disease, rheumatoid arthritis and type 1 diabetes, and TSPAN8 for type 2 diabetes-although in each case the locus had previously been identified in single nucleotide polymorphism (SNP)-based studies, reflecting our observation that most common CNVs that are well-typed on our array are well tagged by SNPs and so have been indirectly explored through SNP studies. We conclude that common CNVs that can be typed on existing platforms are unlikely to contribute greatly to the genetic basis of common human diseases.
|
|
| 7. |
- Jones, Benedict C, et al.
(författare)
-
To which world regions does the valence-dominance model of social perception apply?
- 2021
-
Ingår i: Nature Human Behaviour. - : Springer Science and Business Media LLC. - 2397-3374. ; 5:1, s. 159-169
-
Tidskriftsartikel (refereegranskat)abstract
- Over the past 10 years, Oosterhof and Todorov's valence-dominance model has emerged as the most prominent account of how people evaluate faces on social dimensions. In this model, two dimensions (valence and dominance) underpin social judgements of faces. Because this model has primarily been developed and tested in Western regions, it is unclear whether these findings apply to other regions. We addressed this question by replicating Oosterhof and Todorov's methodology across 11 world regions, 41 countries and 11,570 participants. When we used Oosterhof and Todorov's original analysis strategy, the valence-dominance model generalized across regions. When we used an alternative methodology to allow for correlated dimensions, we observed much less generalization. Collectively, these results suggest that, while the valence-dominance model generalizes very well across regions when dimensions are forced to be orthogonal, regional differences are revealed when we use different extraction methods and correlate and rotate the dimension reduction solution. PROTOCOL REGISTRATION: The stage 1 protocol for this Registered Report was accepted in principle on 5 November 2018. The protocol, as accepted by the journal, can be found at https://doi.org/10.6084/m9.figshare.7611443.v1 .
|
|
| 8. |
- Gaulton, Kyle J, et al.
(författare)
-
Genetic fine mapping and genomic annotation defines causal mechanisms at type 2 diabetes susceptibility loci.
- 2015
-
Ingår i: Nature Genetics. - : Springer Science and Business Media LLC. - 1546-1718 .- 1061-4036. ; 47:12, s. 1415-1415
-
Tidskriftsartikel (refereegranskat)abstract
- We performed fine mapping of 39 established type 2 diabetes (T2D) loci in 27,206 cases and 57,574 controls of European ancestry. We identified 49 distinct association signals at these loci, including five mapping in or near KCNQ1. 'Credible sets' of the variants most likely to drive each distinct signal mapped predominantly to noncoding sequence, implying that association with T2D is mediated through gene regulation. Credible set variants were enriched for overlap with FOXA2 chromatin immunoprecipitation binding sites in human islet and liver cells, including at MTNR1B, where fine mapping implicated rs10830963 as driving T2D association. We confirmed that the T2D risk allele for this SNP increases FOXA2-bound enhancer activity in islet- and liver-derived cells. We observed allele-specific differences in NEUROD1 binding in islet-derived cells, consistent with evidence that the T2D risk allele increases islet MTNR1B expression. Our study demonstrates how integration of genetic and genomic information can define molecular mechanisms through which variants underlying association signals exert their effects on disease.
|
|
| 9. |
- Regev, A, et al.
(författare)
-
The Human Cell Atlas
- 2017
-
Ingår i: eLife. - : ELIFE SCIENCES PUBLICATIONS LTD. - 2050-084X. ; 6
-
Tidskriftsartikel (refereegranskat)
|
|
| 10. |
- Dunham, I, et al.
(författare)
-
The DNA sequence of human chromosome 22
- 1999
-
Ingår i: Nature. - : Springer Science and Business Media LLC. - 0028-0836 .- 1476-4687. ; 402:6761, s. 489-495
-
Tidskriftsartikel (refereegranskat)
|
|
| 11. |
- Tobin, John J., et al.
(författare)
-
The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars. II. A Statistical Characterization of Class 0 and Class i Protostellar Disks
- 2020
-
Ingår i: Astrophysical Journal. - : American Astronomical Society. - 1538-4357 .- 0004-637X. ; 890:2
-
Tidskriftsartikel (refereegranskat)abstract
- We have conducted a survey of 328 protostars in the Orion molecular clouds with the Atacama Large Millimeter/submillimeter Array at 0.87 mm at a resolution of ∼0.″1 (40 au), including observations with the Very Large Array at 9 mm toward 148 protostars at a resolution of ∼0.″08 (32 au). This is the largest multiwavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and flat-spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and flat-spectrum protostars are 44.9-3.4+5.8, 37.0-3.0+4.9, and 28.5-2.3+3.7 au, respectively, and the mean protostellar dust disk masses are 25.9-4.0+7.7, 14.9-2.2+3.8, 11.6-1.9+3.5 M⊙, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% of 379 detected 0.87 mm continuum sources plus 42 nondetections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than those of Class II disks by a factor of >4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase.
|
|
| 12. |
- Archer, E., et al.
(författare)
-
Biodiversity and ecosystem services on the African continent - What is changing, and what are our options?
- 2021
-
Ingår i: Environmental Development. - : Elsevier BV. - 2211-4645 .- 2211-4653. ; 37
-
Tidskriftsartikel (refereegranskat)abstract
- Throughout the world, biodiversity and nature's contributions to people are under threat, with clear changes evident. Biodiversity and ecosystem services have particular value in Africa- yet they are negatively impacted by a range of drivers, including land use and climate change. In this communication, we show evidence of changing biodiversity and ecosystem services in Africa, as well as the current most significant drivers of change. We then consider five plausible futures for the African continent, each underlain by differing assumptions. In three out of the five futures under consideration, negative impacts on biodiversity and ecosystem services are likely to persist. Those two plausible futures prioritizing environment and sustainability, however, are shown as the most likely paths to achieving long term development objectives without compromising the continent's biodiversity and ecosystem services. Such a finding shows clearly that achievement of such objectives cannot be separated from full recognition of the value of such services.
|
|
| 13. |
- Guerriero, S., et al.
(författare)
-
Systematic approach to sonographic evaluation of the pelvis in women with suspected endometriosis, including terms, definitions and measurements : A consensus opinion from the International Deep Endometriosis Analysis (IDEA) group
- 2016
-
Ingår i: Ultrasound in Obstetrics & Gynecology. - : Wiley. - 0960-7692. ; 48:3, s. 318-332
-
Tidskriftsartikel (refereegranskat)abstract
- The IDEA (International Deep Endometriosis Analysis group) statement is a consensus opinion on terms, definitions and measurements that may be used to describe the sonographic features of the different phenotypes of endometriosis. Currently, it is difficult to compare results between published studies because authors use different terms when describing the same structures and anatomical locations. We hope that the terms and definitions suggested herein will be adopted in centers around the world. This would result in consistent use of nomenclature when describing the ultrasound location and extent of endometriosis. We believe that the standardization of terminology will allow meaningful comparisons between future studies in women with an ultrasound diagnosis of endometriosis and should facilitate multicenter research.
|
|
| 14. |
- Moncalvo, J. M., et al.
(författare)
-
The cantharelloid clade: dealing with incongruent gene trees and phylogenetic reconstruction methods
- 2006
-
Ingår i: Mycologia. - : Informa UK Limited. - 0027-5514 .- 1557-2536. ; 98:6, s. 937-948
-
Tidskriftsartikel (refereegranskat)abstract
- We reassessed the circumscription of the cantharelloid clade and identified monophyletic groups by using nLSU, nSSU, mtSSU and RPB2 sequence data. Results agreed with earlier studies that placed the genera Cantharellus, Craterellus, Hydnum, Clavulina, Membranomyces, Multiclavula, Sistotrema, Botryobasidium and the family Ceratobasidiaceae in that clade. Phylogenetic analyses support monophyly of all genera except Sistotrema, which was highly polyphyletic. Strongly supported monophyletic groups were: (i) Cantharellus-Craterellus, Hydnum, and the Sistotrema confluens group; (ii) Clavulina-Membranomyces and the S. brinkmannii-oblongisporum group, with Multiclavula being possibly sister of that clade; (iii) the Sistotrema eximum-octosporum group; (iv) Sistotrema adnatum and S. coronilla. Positions of Sistotrema raduloides and S. athelioides were unresolved, as were basal relationships. Botryobasidium was well supported as the sister taxon of all the above taxa, while Ceratobasidiaceae was the most basal lineage. The relationship between Tulasnella and members of the cantharelloid clade will require further scrutiny, although there is cumulative evidence that they are probably sister groups. The rates of molecular evolution of both the large and small nuclear ribosomal RNA genes (nuc-rDNA) are much higher in Cantharellus, Craterellus and Tulasnella than in the other cantharelloid taxa, and analyses of nuc-rDNA sequences strongly placed Tulasnella close to Cantharellus-Craterellus. In contrast analyses with RPB2 and mtSSU sequences placed Tulasnella at the base of the cantharelloid clade. Our attempt to reconstruct a "supertree" from tree topologies resulting from separate analyses that avoided phylogenetic reconstruction problems associated with missing data and/or unalignable sequences proved unsuccessful.
|
|
| 15. |
- Moshontz, Hannah, et al.
(författare)
-
The Psychological Science Accelerator: Advancing Psychology Through a Distributed Collaborative Network
- 2018
-
Ingår i: Advances in Methods and Practices in Psychological Science. - : SAGE Publications. - 2515-2459 .- 2515-2467. ; 1:4, s. 501-515
-
Tidskriftsartikel (refereegranskat)abstract
- Concerns about the veracity of psychological research have been growing. Many findings in psychological science are based on studies with insufficient statistical power and nonrepresentative samples, or may otherwise be limited to specific, ungeneralizable settings or populations. Crowdsourced research, a type of large-scale collaboration in which one or more research projects are conducted across multiple lab sites, offers a pragmatic solution to these and other current methodological challenges. The Psychological Science Accelerator (PSA) is a distributed network of laboratories designed to enable and support crowdsourced research projects. These projects can focus on novel research questions or replicate prior research in large, diverse samples. The PSA’s mission is to accelerate the accumulation of reliable and generalizable evidence in psychological science. Here, we describe the background, structure, principles, procedures, benefits, and challenges of the PSA. In contrast to other crowdsourced research networks, the PSA is ongoing (as opposed to time limited), efficient (in that structures and principles are reused for different projects), decentralized, diverse (in both subjects and researchers), and inclusive (of proposals, contributions, and other relevant input from anyone inside or outside the network). The PSA and other approaches to crowdsourced psychological science will advance understanding of mental processes and behaviors by enabling rigorous research and systematic examination of its generalizability.
|
|
| 16. |
- Tobin, J. J., et al.
(författare)
-
A triple protostar system formed via fragmentation of a gravitationally unstable disk
- 2016
-
Ingår i: Nature. - : Springer Science and Business Media LLC. - 0028-0836 .- 1476-4687. ; 538:7626, s. 483-
-
Tidskriftsartikel (refereegranskat)abstract
- Binary and multiple star systems are a frequent outcome of the star formation process(1,2) and as a result almost half of all stars with masses similar to that of the Sun have at least one companion star(3). Theoretical studies indicate that there are two main pathways that can operate concurrently to form binary/multiple star systems: large-scale fragmentation of turbulent gas cores and filaments(4,5) or smaller-scale fragmentation of a massive protostellar disk due to gravitational instability(6,7). Observational evidence for turbulent fragmentation on scales of more than 1,000 astronomical units has recently emerged(8,9). Previous evidence for disk fragmentation was limited to inferences based on the separations of more-evolved pre-main sequence and protostellar multiple systems(10-13). The triple protostar system L1448 IRS3B is an ideal system with which to search for evidence of disk fragmentation as it is in an early phase of the star formation process, it is likely to be less than 150,000 years old(14) and all of the protostars in the system are separated by less than 200 astronomical units. Here we report observations of dust and molecular gas emission that reveal a disk with a spiral structure surrounding the three protostars. Two protostars near the centre of the disk are separated by 61 astronomical units and a tertiary protostar is coincident with a spiral arm in the outer disk at a separation of 183 astronomical units(13). The inferred mass of the central pair of protostellar objects is approximately one solar mass, while the disk surrounding the three protostars has a total mass of around 0.30 solar masses. The tertiary protostar itself has a minimum mass of about 0.085 solar masses. We demonstrate that the disk around L1448 IRS3B appears susceptible to disk fragmentation at radii between 150 and 320 astronomical units, overlapping with the location of the tertiary protostar. This is consistent with models for a protostellar disk that has recently undergone gravitational instability, spawning one or two companion stars.
|
|
| 17. |
- Tobin, John J., et al.
(författare)
-
The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars. I. Identifying and Characterizing the Protostellar Content of the OMC-2 FIR4 and OMC-2 FIR3 Regions
- 2019
-
Ingår i: Astrophysical Journal. - : American Astronomical Society. - 1538-4357 .- 0004-637X. ; 886:1
-
Tidskriftsartikel (refereegranskat)abstract
- We present Atacama Large Millimeter/submillimeter Array (0.87 mm) and Very Large Array (9 mm) observations toward OMC-2 FIR4 and OMC-2 FIR3 within the Orion integral-shaped filament, thought to be two of the nearest regions of intermediate-mass star formation. We characterize the continuum sources within these regions on ?40 au (01) scales and associated molecular line emission at a factor of ?30 better resolution than previous observations at similar wavelengths. We identify six compact continuum sources within OMC-2 FIR4, four in OMC-2 FIR3, and one additional source just outside OMC-2 FIR4. This continuum emission is tracing the inner envelope and/or disk emission on less than 100 au scales. HOPS-108 is the only protostar in OMC-2 FIR4 that exhibits emission from high-excitation transitions of complex organic molecules (e.g., methanol and other lines) coincident with the continuum emission. HOPS-370 in OMC-2 FIR3, with L;?;360 L, also exhibits emission from high-excitation methanol and other lines. The methanol emission toward these two protostars is indicative of temperatures high enough to thermally evaporate it from icy dust grains; overall, these protostars have characteristics similar to hot corinos. We do not identify a clear outflow from HOPS-108 in (CO)-C-12, but we find evidence of interaction between the outflow/jet from HOPS-370 and the OMC-2 FIR4 region. A multitude of observational constraints indicate that HOPS-108 is likely a low- to intermediate-mass protostar in its main mass accretion phase and is the most luminous protostar in OMC-2 FIR4. The high-resolution data presented here are essential for disentangling the embedded protostars from their surrounding dusty environments and characterizing them.
|
|
| 18. |
- Svoboda, B. E., et al.
(författare)
-
THE BOLOCAM GALACTIC PLANE SURVEY. XIV. PHYSICAL PROPERTIES of MASSIVE STARLESS and STAR-FORMING CLUMPS
- 2016
-
Ingår i: Astrophysical Journal. - 0004-637X .- 1538-4357. ; 822:2
-
Tidskriftsartikel (refereegranskat)abstract
- We sort 4683 molecular clouds between 10° < ℓ < 65° from the Bolocam Galactic Plane Survey based on observational diagnostics of star formation activity: compact 70 μm sources, mid-IR color-selected YSOs, H2O and CH3OH masers, and UCH ii regions. We also present a combined NH3-derived gas kinetic temperature and H2O maser catalog for 1788 clumps from our own GBT 100 m observations and from the literature. We identify a subsample of 2223 (47.5%) starless clump candidates (SCCs), the largest and most robust sample identified from a blind survey to date. Distributions of flux density, flux concentration, solid angle, kinetic temperature, column density, radius, and mass show strong (>1 dex) progressions when sorted by star formation indicator. The median SCC is marginally subvirial (α ∼ 0.7) with >75% of clumps with known distance being gravitationally bound (α < 2). These samples show a statistically significant increase in the median clump mass of ΔM ∼ 170-370 M o from the starless candidates to clumps associated with protostars. This trend could be due to (i) mass growth of the clumps at M o Myr-1 for an average freefall 0.8 Myr timescale, (ii) a systematic factor of two increase in dust opacity from starless to protostellar phases, and/or (iii) a variation in the ratio of starless to protostellar clump lifetime that scales as ∼M -0.4. By comparing to the observed number of CH3OH maser containing clumps, we estimate the phase lifetime of massive (M > 103 M o) starless clumps to be 0.37 ± 0.08 Myr (M/103 M o)-1; the majority (M < 450 M o) have phase lifetimes longer than their average freefall time. © 2016. The American Astronomical Society. All rights reserved.
|
|
| 19. |
|
|
| 20. |
|
|
| 21. |
- Rada-Iglesias, Alvaro, et al.
(författare)
-
Butyrate mediates decrease of histone acetylation centered on transcription start sites and down-regulation of associated genes
- 2007
-
Ingår i: Genome Research. - : Cold Spring Harbor Laboratory. - 1088-9051 .- 1549-5469. ; 17:6, s. 708-719
-
Tidskriftsartikel (refereegranskat)abstract
- Butyrate is a histone deacetylase inhibitor (HDACi) with anti-neoplastic properties, which theoretically reactivates epigenetically silenced genes by increasing global histone acetylation. However, recent studies indicate that a similar number or even more genes are down-regulated than up-regulated by this drug. We treated hepatocarcinoma HepG2 cells with butyrate and characterized the levels of acetylation at DNA-bound histones H3 and H4 by ChIP-chip along the ENCODE regions. In contrast to the global increases of histone acetylation, many genomic regions close to transcription start sites were deacetylated after butyrate exposure. In order to validate these findings, we found that both butyrate and trichostatin A treatment resulted in histone deacetylation at selected regions, while nucleosome loss or changes in histone H3 lysine 4 trimethylation (H3K4me3) did not occur in such locations. Furthermore, similar histone deacetylation events were observed when colon adenocarcinoma HT-29 cells were treated with butyrate. In addition, genes with deacetylated promoters were down-regulated by butyrate, and this was mediated at the transcriptional level by affecting RNA polymerase II (POLR2A) initiation/elongation. Finally, the global increase in acetylated histones was preferentially localized to the nuclear periphery, indicating that it might not be associated to euchromatin. Our results are significant for the evaluation of HDACi as anti-tumourogenic drugs, suggesting that previous models of action might need to be revised, and provides an explanation for the frequently observed repression of many genes during HDACi treatment.
|
|
| 22. |
|
|
| 23. |
- Almquist, Martin, et al.
(författare)
-
Elastic wave propagation in anisotropic solids using energy-stable finite differences with weakly enforced boundary and interface conditions
- 2021
-
Ingår i: Journal of Computational Physics. - : Elsevier. - 0021-9991 .- 1090-2716. ; 424
-
Tidskriftsartikel (refereegranskat)abstract
- Summation-by-parts (SBP) finite difference methods have several desirable properties for second-order wave equations. They combine the computational efficiency of narrow-stencil finite difference operators with provable stability on curvilinear multiblock grids. While several techniques for boundary and interface conditions exist, weak imposition via simultaneous approximation terms (SATs) is perhaps the most flexible one. Although SBP methods have been applied to elastic wave equations many times, an SBP-SAT method for general anisotropic elastic wave equations has not yet been presented in the literature. We fill this gap by deriving energy-stable self-adjoint SBP-SAT methods for general anisotropic materials on curvilinear multiblock grids. The methods are based on fully compatible SBP operators. Although this paper focuses on classical SBP finite difference operators, the presented boundary and interface treatments are general and apply to a range of methods that satisfy an SBP property. We demonstrate the stability and accuracy properties of a particular set of fully compatible SBP-SAT schemes using the method of manufactured solutions. We also demonstrate the utility of the new method in elastodynamic cloaking and seismic imaging in mountainous regions.
|
|
| 24. |
- Almquist, Martin, et al.
(författare)
-
Non-stiff boundary and interface penalties for narrow-stencil finite difference approximations of the Laplacian on curvilinear multiblock grids
- 2020
-
Ingår i: Journal of Computational Physics. - : Elsevier BV. - 0021-9991 .- 1090-2716. ; 408
-
Tidskriftsartikel (refereegranskat)abstract
- The Laplacian appears in several partial differential equations used to model wave propagation. Summation-by-parts simultaneous approximation term (SBP-SAT) finite difference methods are often used for such equations, as they combine computational efficiency with provable stability on curvilinear multiblock grids. However, the existing SBP-SAT discretization of the Laplacian quickly becomes prohibitively stiff as grid skewness increases. The stiffness stems from the SATs that impose inter-block couplings and Dirichlet boundary conditions. We resolve this issue by deriving stable SATs whose stiffness is almost insensitive to grid skewness. The new discretization thus allows for large time steps in explicit time integrators, even on very skewed grids. It also applies to the variable-coefficient generalization of the Laplacian. We demonstrate the efficacy and versatility of the new SATs by applying them to acoustic wave propagation problems inspired by marine seismic exploration and infrasound monitoring of volcanoes.
|
|
| 25. |
- Bader, Milad, et al.
(författare)
-
Modeling and inversion in acoustic-elastic coupled media using energy-stable summation-by-parts operators
- 2023
-
Ingår i: Geophysics. - : SOC EXPLORATION GEOPHYSICISTS - SEG. - 0016-8033 .- 1942-2156. ; 88:3
-
Tidskriftsartikel (refereegranskat)abstract
- Seismic data acquired at the seafloor are valuable in charac-terizing the subsurface and monitoring producing hydrocarbon fields. To fully use such data, a stable, accurate, and efficient numerical scheme is needed that accounts for acoustic and elas-tic wave propagation, their interaction at the seafloor interface, and for sources and receivers placed on either side of that inter-face. Existing methods either make incorrect assumptions or have high implementation and computational costs. We have developed a high-order finite-difference summation-by-parts framework for the acoustic-elastic wave equations in second -or-der form (in terms of displacements in the solid and an extension of velocity potential in the fluid). Our modified discretization of the elastic operator overcomes the dispersion errors known to plague displacement-based schemes in the high VP/VS limit. We weakly impose boundary and interface conditions using simulta-neous approximation terms, leading to an energy-stable numeri-cal scheme that rigorously handles point injection and extraction. The fully discrete system is self-adjoint after a time reversal and a sign flip and is furthermore a high-order accurate discretization of the continuous problem. The self-adjointness ensures that for-ward and adjoint wavefields are computed with similar accuracy and simplifies gradient computation for inversion purposes. We find that our numerical scheme achieves accuracy comparable to a more computationally expensive spectral-element method and demonstrate its application to full-waveform inversion using the Marmousi2 model.
|
|
